ALBERT

Description: Grain boundary networks of quartz, plagioclase and olivine crystal aggregates in metamorphic rocks have been investigated from the nanometer to the millimeter scale by polarized-light microscopy, SEM, and TEM. The studied materials show different grain sizes and experienced different retrograde P-T histories. The aggregates of quartz and plagioclase are traversed by networks of ∼90% continuously open boundaries with μm-sized cavities along the boundaries or at triple junctions. The boundaries are up to ∼500 nm wide open with typically parallel opposing grain faces. Olivine boundaries are filled with serpentine that does not replace olivine but fills the initially open space homogeneously and mostly with random orientation. For quartz there is no correlation between the crystallographic orientation of grain boundaries and their widths. Amongst all samples analyzed, a weak positive correlation exists between grain size and width of open grain boundaries. The application of measured volume changes and elasticity data from the literature to the cooling-decompression paths of the analyzed materials suggests that fracturing with subsequent widening of the grain boundaries starts at temperatures recognizably below the transition from crystal-plastic to brittle behavior of quartz, plagioclase and olivine but not only under surface conditions. The high amount of open boundaries causes an extensive permeability.

Description: The Cretaceous provides us with an excellent case history of ocean-climate-biota system perturbations. Such perturbations occurred several times during the Cretaceous, such as oceanic anoxic events and the end-Cretaceous mass extinction, which have been the subject of an abundant literature. Other perturbations, such as the mid-Maastrichtian Event (MME) remain poorly understood. The MME was associated with global sea-level rise, changes in climate and deep-water circulation that were accompanied by biotic extinctions including ‘true inoceramids’ and the demise of the Caribbean-Tethyan rudist reef ecosystems. So far, the context and causes behind the MME remain poorly studied. We conducted high-resolution integrated biotic, petrological and geochemical studies in order to fill this knowledge gap. We studied, in particular, carbonate Nd and Os isotopes, whole-rock Hg, C and N content, C and N isotopes in organic matter, S isotopes in carbonate-associated sulfate, along with C and O isotopes in foraminifera from the European Chalk Sea: the Polanówka UW-1 core from Poland and the Stevns-1 core from Denmark. Our data showed that sea-level rise of ∼50–100 m lasted around ∼2 Ma and co-occurred with anomalously high mercury concentration in seawater. Along with previously published data, our results strongly suggest that the MME was driven by intense volcanic–tectonic activity, likely related to the production of vast oceanic plateaus (LIP, Large Igneous Province). The collapse of reef ecosystems could have been the consequence of LIP-related environmental stress factors, including climate warming, presumably caused by emission of greenhouse gases, modification of the oceanic circulation, oceanic acidification and/or toxic metal input. The disappearance of the foraminifer Stensioeina lineage on the European shelf was likely caused by the collapse of primary production triggered by sea-level rise and limited amount of nutrient input. Nd isotopes and foraminiferal assemblages attest for changes in sea-water circulation in the European Shelf and the increasing contribution of North Atlantic water masses

Description: The structural response to compression of the synthetic high-pressure hydroxide perovskite MgSi(OH)6, the so-called “3.65 Å phase,” has been determined to 8.4 GPa at room temperature using single-crystal XRD in the diamond-anvil cell. Two very similar structures have been determined in space groups P21 and P21/n, for which differences in oxygen donor-acceptor distances indicate that the non-centrosymmetric structure is likely the correct one. This structure has six nonequivalent H sites, of which two are fully occupied and four are half-occupied. Half-occupied sites are associated with a well-defined crankshaft of hydrogen-bonded donor-acceptor oxygens extending parallel to c. Half occupancy of these sites arises from the averaging of two orientations of the crankshaft H atoms (|| ±c) in equal proportions. The P21 and P21/n structures are compared. It is shown that the former is likely the correct space group, which is also consistent with recent spectroscopic studies that recognize six nonequivalent O-H. The structure of MgSi(OH)6 at pressures up to 8.4 GPa was refined in both space groups to see how divergent the two models are. There is a very close correspondence between the responses of the two structures implying that, at least to 8.4 GPa, non-centrosymmetry does not affect compressional behavior. The very different compressional behavior of MgO6 and SiO6 octahedra observed in this study suggests that structural phase transformations or discontinuities likely occur in MgSi(OH)6 above 9 GPa.

Description: The Combination Service for Time-variable Gravity fields (COST-G) operationally provides combinations of monthly Earth gravity field models derived from observations of the microwave ranging instrument of the GRACE Follow-on (GRACE-FO) satellite mission, applying the quality control and combination methodology originally developed by the Horizon 2020 project European Gravity Service for Improved Emergency Management for the data of the GRACE satellites. In the frame of the follow-up Horizon 2020 project Global Gravity-based Groundwater Product (G3P), the GRACE-FO combination is used to derive global grids of groundwater storage anomalies. To meet the user requirements and achieve optimal signal-to-noise ratio, the combination has been further developed and extended to incorporate: • new time-series based on the alternative accelerometer transplant product generated in the frame of the project by the Institute of Geodesy at the Graz University of Technology, which specifically improves the estimation of the C30 coefficient and also reduces the noise at medium to short wavelengths, and • the new time-series AIUB–GRACE-FO–RL02 of monthly GRACE-FO gravity fields, which is derived at the Astronomical Institute of the University of Bern by applying empirical noise modelling techniques. The COST-G quality control confirms the consistency of the contributing GRACE-FO time-series concerning the signal amplitude of seasonal hydrology in large river basins and the secular mass change in polar regions, but it also indicates rather diverse noise characteristics. The difference in the noise levels is taken into account in the combination process by relative weights derived by variance component estimation on the solution level. The weights are expected to be inverse proportional to the noise levels of the individual gravity field solutions. However, this expectation is violated when applying the weighting scheme as developed for the GRACE combination. The reason is found in the high-order coefficients of the gravity field, which are poorly determined from the low–low range-rate observations due to the observation geometry and suffer from aliasing due to the malfunctioning accelerometer onboard one of the GRACE-FO satellites. Hence, for the final G3P-combination a revised weighting scheme is applied where the gravity field coefficients beyond order 60 are excluded from the determination of the weights. The quality of the combined gravity fields is assessed by comparison of the noise content and the signal-to-noise ratio with the individual time-series. Independent validation is provided by the COST-G validation centre at the GFZ German Research Centre for Geosciences, where orbit fits of the low-flying Gravity and steady-state Ocean Circulation Explorer satellite are performed that confirm the high quality of the combined GRACE-FO gravity fields. By the end of the G3P project, the new combination scheme is implemented by COST-G as the new COST-G–GRACE-FO–RL02 and continued to be used for the operational GRACE-FO combination.

Description: With the ongoing deployment of Global Navigation Satellite Systems (GNSS) ground stations and the modernization of satellite signal systems, the utilization of various augmentation technologies enables the realization of Precise Point Positioning (PPP) in real-time. Augmentation technology, which introduces precise atmospheric and signal-related delays, has become an essential component of high-precision real-time services and is attracting growing interest in scientific research, disaster monitoring, autopilot, etc. Previous studies have dedicated significant efforts to enhance the generation and dissemination of augmentation information on the service side and improve real-time positioning algorithms on the user side. The real-time atmosphere augmentation information with sufficient accuracy and proper constraint, and reliable Ambiguity Resolution (AR) for this purpose is the main focus of current GNSS research. However, these efforts have primarily been concentrated on small or medium-sized regions with the capability for transmitting massive data volumes. Alternatively, they have focused on larger areas, but with slow convergence due to the imprecise nature of atmosphere information. To address the challenge posed by the trade-offs among service area size, correction volume, and the precision of represented correction, a new augmentation strategy is proposed. This approach integrates the advantages of atmospheric delay fitting models, unmodeled residuals, and uncertainty information to achieve rapid and high-precision positioning, all while reducing data transmission volume for larger areas. It also allows users to implement different positioning modes depending on their communication capacity. Additionally, all deviations among different types of receivers and satellite signals are calibrated in this study for reliable AR can be achieved on all reference stations. The main contribution of this thesis is summarized as follows. With the real-time precise orbit, clock, and Uncalibrated Phase Delay (UPD) products, precise atmospheric delay corrections relying on reliable AR can be derived for large-areas augmentation services. To address the challenge of achieving reliable AR across different receiver types and various satellite signals, this thesis proposes a comprehensive method for calibrating receiver-type-related satellite-specific deviations and analyzes the impact of satellite signal bias corrections in data processing. The primary objective is to enhance the reliability of AR, enabling the utilization of all available signals and receiver types in large-area services. Subsequently, new tropospheric and ionospheric delay fitting models applied for large-area are carried out according to the properties of their propagation paths. In addition, the corresponding atmospheric delay uncertainty for large areas is introduced based on the fitting residuals. Finally, a hierarchical mode is developed for augmentation services, leveraging the advantages of the fitting model and uncertainty grid to reduce data volume and incorporating regional fitting residuals using the interpolation model and ionospheric delay error function, depending on the network capability. Based on hierarchical augmentation, positioning in large areas can not only achieve rapid/instantaneous high-precision convergence but also overcome the conflict among correction volume, represented precision, and coverage size. In order to derive precise atmospheric delay and accelerate positioning, implementing reliable and robust AR across all types of receivers and satellite signals is essential. It also demonstrates and discusses the advantages of calibrating satellite-signal and receiver-type-related satellite-specific deviations in AR solutions. The deviations related to receivers in terms of UPD products are assessed and calibrated, confirming that a 0.03 cycle consistency in wide lane UPD can be achieved. The effectiveness of the proposed approach is demonstrated using GPS satellite signals, which can improve the AR rate by at least 10% and produce more reliable results. In addition, the impact of different signal settings and corrections on orbit, clock, and UPD generation, as well as positioning and pseudo-range signal systematic and stochastic residuals, is analyzed. These processing strategies provide flexible observation selections, allowing the utilization of all available satellite signals and receiver types, thereby enabling reliable AR and a higher fixing rate. As a result, an AR fixing rate exceeding 95% is achievable across all stations in large-area services. For precise atmospheric delay modeling over large areas, new models are proposed, including a tropospheric Zenith Wet Delay (ZWD) model and a satellite-wise ionospheric slant delay fitting model. The tropospheric delay model takes the exponential function of water vapor vertical changes into consideration, addressing model anomalies in areas with large altitude differences. The new ionospheric delay fitting model introduces the trigonometric functions to describe differences in slant path delays between the optimal reference propagation path and others, achieving superior modeling performance in large areas. The precision of the fitting model, utilizing a 200 km station-spacing network, demonstrates tropospheric ZWD and ionospheric slant delays of 1.3 cm and 8.9 cm, respectively, with smaller standard deviations. These new fitting models overcome the challenge of handling massive information for providing station-wise corrections and avoid an increase in the number of coefficients. In addition to the function model, the stochastic model, i.e., uncertainty information, is essential for describing the quality of corrections. The atmospheric delay uncertainty for the large-area fitting model is generated based on the fitting residuals and represented in forms of grid-point. Additionally, regional ionosphere unmodeled residual uncertainty is represented by the form of liner function, which is established by the relationship between distance and interpolation precision through inter-satellite cross-verification among all reference stations. The differences between uncertainty value and real delays are 2.5 cm and 0.5 cm for grid and function forms, respectively. For real-time applications in large areas, the fitting model and grid-based atmosphere uncertainty serve as the essential information, satisfying the requirement of rapid positioning. By further incorporating unmodeled residuals and ionosphere error function, a hierarchical augmentation model is provided. Based on the fitting model established for large areas, unmodeled residuals are further introduced as optional compensation for specific areas, depending on the magnitude of fitting residuals. This approach results in a 97% reduction in tropospheric delay and a 65% reduction in ionospheric delay transmission volume. Furthermore, leveraging the regional high capability of communication, 85.3% of all solutions can achieve instantaneous convergence at the first epoch with the aid of corresponding regional compensation. This thesis proposes a large areas augmentation service to overcome the conflict among correction data volume, represented precision, and coverage size. It demonstrates the benefits of an augmentation mode that integrates regional information into large-area services. Under these conditions, a more reliable and rapid AR solution can be easily achieved based on precise atmospheric delay correction and uncertainty in large areas with fewer data volume requirements. This is beneficial for actual real-time services and applications.

Description: Mit der laufenden Bereitstellung von Bodenstationen für globale Navigationssatellitensysteme (GNSS) und der Modernisierung von Satellitensignal-Systemen ermöglicht die Nutzung verschiedener Augmentationstechnologien die Realisierung der Präzisen Punkt-Positionierung (PPP) in Echtzeit. Augmentationstechnologie, die präzise atmosphärische und signalbezogene Verzögerungen einführt, ist zu einem wesentlichen Bestandteil hochpräziser Echtzeitdienste geworden und findet wachsendes Interesse in wissenschaftlicher Forschung, Katastrophenüberwachung, Autopiloten usw. Frühere Studien haben erhebliche Anstrengungen darauf verwendet, die Erzeugung und Verbreitung von Augmentationsinformationen auf der Dienstseite zu verbessern und Echtzeit-Positionierungsalgorithmen auf der Benutzerseite zu optimieren. Die Echtzeit-Atmosphärenaugmentationsinformationen mit ausreichender Genauigkeit und angemessener Einschränkung sowie zuverlässige Ambiguitätsauflösung (AR) für diesen Zweck stehen im Mittelpunkt der aktuellen GNSS-Forschung. Diese Bemühungen konzentrierten sich jedoch hauptsächlich auf kleine oder mittelgroße Regionen mit der Fähigkeit zur Übertragung großer Datenmengen. Alternativ richteten sie sich auf größere Gebiete, jedoch mit langsamer Konvergenz aufgrund der ungenauen Natur der Atmosphäreninformation. Um der Herausforderung durch die Abwägung zwischen Größe des Dienstleistungsgebiets, Korrekturvolumen und Präzision der dargestellten Korrektur zu begegnen, wird eine neue Augmentationsstrategie vorgeschlagen. Dieser Ansatz integriert die Vorteile atmosphärischer Verzögerungsanpassungsmodelle, nicht modellierter Reste und Unsicherheitsinformationen, um eine schnelle und hochpräzise Positionierung zu erreichen, und das bei gleichzeitiger Reduzierung der Datenübertragungsvolumina für größere Gebiete. Es ermöglicht den Benutzern auch, verschiedene Positionierungsmodi je nach ihrer Kommunikationskapazität zu implementieren. Zusätzlich werden in dieser Studie alle Abweichungen zwischen verschiedenen Typen von Empfängern und Satellitensignalen kalibriert, um eine zuverlässige AR an allen Referenzstationen zu erreichen. Die Hauptbeiträge dieser Arbeit werden wie folgt zusammengefasst. Mit den Echtzeit-Präzbitbahnen, Uhren und Uncalibrated Phase Delay (UPD)-Produkten können präzise atmosphärische Verzögerungskorrekturen für großflächige Augmentationsdienste abgeleitet werden, die auf zuverlässiger AR basieren. Um die Herausforderung zu bewältigen, eine zuverlässige AR über verschiedene Empfängertypen und verschiedene Satellitensignale hinweg zu erreichen, schlägt diese Arbeit eine umfassende Methode zur Kalibrierung von empfängertypbezogenen satellspezifischen Abweichungen vor und analysiert die Auswirkungen von Korrekturen für Satellitensignalverzerrungen in der Datenverarbeitung. Das Hauptziel besteht darin, die Zuverlässigkeit der AR zu verbessern und die Nutzung aller verfügbaren Signale und Empfängertypen in großflächigen Diensten zu ermöglichen. Anschließend werden neue troposphärische und ionosphärische Verzögerungsanpassungsmodelle für großflächige Anwendungen gemäß den Eigenschaften ihrer Ausbreitungspfade durchgeführt. Darüber hinaus wird die entsprechende atmosphärische Verzögerungsunsicherheit für große Gebiete auf der Grundlage der Anpassungsreste eingeführt. Schließlich wird ein hierarchischer Modus für Augmentationsdienste entwickelt, der die Vorteile des Anpassungsmodells und des Unsicherheitsgitters nutzt, um das Datenvolumen zu reduzieren und regionale Anpassungsreste unter Verwendung des Interpolationsmodells und der ionosphärischen Verzögerungsfehlerfunktion, abhängig von der Netzwerkfähigkeit, zu integrieren. Basierend auf der hierarchischen Augmentation kann die Positionierung in großen Gebieten nicht nur eine schnelle/instantane hochpräzise Konvergenz erreichen, sondern auch den Konflikt zwischen Korrekturvolumen, dargestellter Präzision und Abdeckungsgröße überwinden. Um präzise atmosphärische Verzögerungen abzuleiten und die Positionierung zu beschleunigen, ist es entscheidend, eine zuverlässige und robuste AR über alle Arten von Empfängern und Satellitensignalen zu implementieren. Es zeigt auch die Vorteile der Kalibrierung von satellitensignal- und empfängertypbezogenen satellspezifischen Abweichungen in AR-Lösungen auf. Die Abweichungen im Zusammenhang mit Empfängern in Bezug auf UPD-Produkte werden bewertet und kalibriert, wobei bestätigt wird, dass eine Konsistenz von 0,03 Zyklen bei Wide-Lane-UPD erreicht werden kann. Die Wirksamkeit des vorgeschlagenen Ansatzes wird unter Verwendung von GPS-Satellitensignalen demonstriert, die die AR-Rate um mindestens 10% verbessern und zu zuverlässigeren Ergebnissen führen können. Darüber hinaus wird der Einfluss unterschiedlicher Signalparameter und Korrekturen auf die Erzeugung von Orbit, Uhr und UPD sowie auf die Positionierung und systematische und stochastische Reste der Pseudo-Range-Signale analysiert. Diese Verarbeitungsstrategien bieten flexible Auswahlmöglichkeiten bei der Beobachtung und ermöglichen die Nutzung aller verfügbaren Satellitensignale und Empfängertypen, wodurch eine zuverlässige AR und eine höhere Fixierungsrate ermöglicht wird. Als Ergebnis ist eine AR-Fixierungsrate von über 95% bei allen Stationen in großflächigen Diensten erreichbar. Für eine präzise Modellierung atmosphärischer Verzögerungen über großen Gebieten werden neue Modelle vorgeschlagen, darunter ein troposphärisches Zenith Wet Delay (ZWD)-Modell und ein satellitenweises ionosphärisches Schrägverzögerungsanpassungsmodell. Das troposphärische Verzögerungsmodell berücksichtigt die exponentielle Funktion der vertikalen Änderungen des Wasserdampfs und behebt Modellanomalien in Gebieten mit großen Höhendifferenzen. Das neue ionosphärische Verzögerungsanpassungsmodell verwendet trigonometrische Funktionen, um Unterschiede in den Schrägpfadverzögerungen zwischen dem optimalen Referenzausbreitungspfad und anderen zu beschreiben und erreicht so eine überlegene Modellierungsleistung in großen Gebieten. Die Präzision des Anpassungsmodells, unter Verwendung eines 200 km-Stationen-Netzwerks, zeigt troposphärische ZWD- und ionosphärische Schrägverzögerungen von jeweils 1,3 cm und 8,9 cm mit kleineren Standardabweichungen. Diese neuen Anpassungsmodelle überwinden die Herausforderung, massive Informationen für die Bereitstellung stationsspezifischer Korrekturen zu verarbeiten, und vermeiden eine Zunahme der Anzahl der Koeffizienten. Neben dem Funktionsmodell ist das stochastische Modell, d. h. Unsicherheitsinformationen, entscheidend für die Beschreibung der Qualität der Korrekturen. Die Unsicherheit der atmosphärischen Verzögerung für das großflächige Anpassungsmodell wird auf der Grundlage der Anpassungsreste generiert und in Form von Gitterpunkten dargestellt. Zusätzlich wird die regionale ionosphärische nicht modellierte Restunsicherheit durch die Form einer linearen Funktion repräsentiert, die durch die Beziehung zwischen Entfernung und Interpolationsgenauigkeit durch inter-satellitenkreuz-Verifikation zwischen allen Referenzstationen etabliert wird. Die Unterschiede zwischen Unsicherheitswert und realen Verzögerungen betragen 2,5 cm bzw. 0,5 cm für Gitter- und Funktionsformen. Für Echtzeitanwendungen in großen Gebieten dienen das Anpassungsmodell und die gitterbasierte Atmosphärenunsicherheit als wesentliche Informationen, die die Anforderungen an schnelle Positionierung erfüllen. Durch die weitere Integration von nicht modellierten Resten und Ionosphärenfehlerfunktion wird ein hierarchisches Augmentationsmodell bereitgestellt. Basierend auf dem für große Gebiete etablierten Anpassungsmodell werden nicht modellierte Reste zusätzlich als optionale Kompensation für spezifische Bereiche eingeführt, abhängig von der Größenordnung der Anpassungsreste. Dieser Ansatz führt zu einer Reduktion von 97% der troposphärischen Verzögerung und einer Reduktion von 65% des ionosphärischen Verzögerungsvolumens. Darüber hinaus können unter Nutzung der regionalen hohen Kommunikationsfähigkeit 85,3% aller Lösungen mit Hilfe entsprechender regionaler Kompensation eine sofortige Konvergenz beim ersten Epochenzeitpunkt erreichen. Diese Dissertation schlägt einen großflächigen Augmentationsdienst vor, um den Konflikt zwischen Korrekturvolumen, dargestellter Präzision und Abdeckungsgröße zu überwinden. Sie zeigt die Vorteile eines Augmentationsmodus, der regionale Informationen in großflächige Dienste integriert. Unter diesen Bedingungen kann eine zuverlässigere und schnellere AR-Lösung basierend auf präziser atmosphärischer Verzögerungskorrektur und Unsicherheit in großen Gebieten mit geringeren Anforderungen an das Datenvolumen leicht erreicht werden. Dies ist vorteilhaft für tatsächliche Echtzeitdienste und Anwendungen.

Description: Inland water bodies play a vital role at all scales in the terrestrial water balance and Earth’s climate variability. Thus, an inventory of inland waters is crucially important for hydrologic and ecological studies and management. Therefore, the main aim of this study was to develop a deep learning-based method for inventorying and mapping inland water bodies using the RGB band of high-resolution satellite imagery automatically and accurately. The Sentinel-2 Harmonized dataset, together with ZABAGED-validated ground truth, was used as the main dataset for the model training step. Three different deep learning algorithms based on U-Net architecture were employed to segment inland waters, including a simple U-Net, Residual Attention U-Net, and VGG16-U-Net. All three algorithms were trained using a combination of Sentinel-2 visible bands (Red [B04; 665nm], Green [B03; 560nm], and Blue [B02; 490 nm]) at a 10-meter spatial resolution. The Residual Attention U-Net achieved the highest computational cost due to the increased number of trainable parameters. The VGG16-U-Net had the shortest run time and the lowest number of trainable parameters, attributed to its architecture compared to the simple and Residual Attention U-Net architectures, respectively. As a result, the VGG16-U-Net provided the best segmentation results with a mean-IoU score of 0.9850, a slight improvement compared to other proposed U-Net-based architectures. Although the accuracy of the model based on VGG16-U-Net does not make a difference from Residual Attention U-Net, the computation costs for training VGG16-U-Net were dramatically lower than Residual Attention U-Net.

Description: Northeastern Eurasia is one of the least explored regions in the world. Very little geophysical data is available for this inaccessible area. Even the exact location of the plate boundary between Eurasia and North America remains a subject of ongoing debate. The effective elastic thickness (EET) of the lithosphere is a proxy for lithospheric strength and can provide insight into the thermal regime and tectonic processes. We have computed a high-resolution map of the EET for northeastern Eurasia using the fan wavelet coherence technique applied to the Bouguer gravity anomalies and topography/bathymetry data, appropriately adjusted to account for the influence of density variations within sediments. The results obtained provide insights into different tectonic regimes within this predominantly understudied region. In particular, we identify the boundary between the Eurasian and North American plates in Siberia as a rheologically weak diffusive zone extending from the Verkhoyansk and Sette-Daban Ranges to the eastern boundary of the Chersky Range. Unlike the Sette-Daban and Verkhoyansk Ranges, which were formed by plate collision and have an EET of 30–50 km, other mountainous regions have much lower EET values, usually less than 15 km. These areas have recently experienced tectonic activity that has weakened the lithosphere.

Description: This study introduces a new approach for in situ Rb–Sr dating that utilizes rapid line scans instead of static spot ablation, enabling the creation of two-dimensional 87Rb/86Sr and 87Sr/86Sr isotope ratio and Rb–Sr age maps. The data acquisition is conducted utilizing an ICP-MS/MS instrument with N2O as the reaction gas, coupled to a 193 nm excimer laser via a low-aerosol-dispersion interface. This configuration allows for high repetition rates (〉100 Hz) and sensitivities, enabling data acquisition at a high scanning speed and small laser beam size (3–4 μm). Notably, this approach requires just about 1/30 of the sample volume typically utilized in conventional spot ablation mode, while achieving similar levels of precision and accuracy. Line scan ablation is tested and compared to spot ablation on age-homogeneous crystalline muscovite and biotite, for which reference Rb–Sr age data is acquired through ID-TIMS. Results show that a key requirement for accurate Rb–Sr ages based on line scan analyses is matrix correction using chemically matched crystalline mica. By presenting Rb–Sr age maps of three naturally deformed mica samples, we highlight the potential of Rb–Sr mapping for extracting age data from rocks that exhibit complex metamorphic-metasomatic histories and microscale dynamic recrystallization. Additionally, we show that quantitative elemental information (Al, Fe, Si, Li) can be collected alongside Rb–Sr isotope data. This advancement offers a distinctly more insightful assessment of isotope mobility in natural systems, the timing of element enrichment processes and enables, in high-Rb/Sr rock systems, precise and accurate isotopic dating of intricate geological processes at small scales.

Description: High-resolution flood maps are needed for more effective flood risk assessment and management. Producing these directly with hydrodynamic models is slow and computationally prohibitive at large scales. Here we demonstrate a new algorithm for post-processing low-resolution inundation layers by using high-resolution terrain models to disaggregate or downscale. The new algorithm is roughly 8 times faster than state-of-the-art algorithms and shows a slight improvement in accuracy when evaluated against observations of a recent flood using standard performance metrics. Qualitatively, the algorithm generates more physically coherent flood maps in some hydraulically challenging regions compared to the state of the art. The algorithm developed here is open source and can be applied in conjunction with a low-resolution hydrodynamic model and a high-resolution DEM to rapidly produce high-resolution inundation maps. For example, in our case study with a river reach of 20 km, the proposed algorithm generated a 4 m resolution inundation map from 32 m hydrodynamic model outputs in 33 s compared to a 4 m hydrodynamic model runtime of 34 min. This 60-fold improvement in runtime is associated with a 25 % increase in RMSE when compared against the 4 m hydrodynamic model results and observations of a recent flood. Substituting downscaling into flood risk model chains for high-resolution modelling has the potential to drastically improve the efficiency of inundation map production and increase the lead time of impact-based forecasts, helping more at-risk communities prepare for and mitigate flood damages.

Description: To design user-centred and scientifically high-quality outreach products to inform about earthquake-related hazards and the associated risk, a close collaboration between the model developers and communication experts is needed. In this contribution, we present the communication strategy developed to support the public release of the first openly available European Seismic Risk Model and the updated European Seismic Hazard Model. The backbone of the strategy was the communication concept in which the overall vision, communication principles, target audiences (including personas), key messages, and products were defined. To fulfil the end-users' needs, we conducted two user testing surveys: one for the interactive risk map viewer and one for the risk poster with a special emphasis on the European earthquake risk map. To further ensure that the outreach products are not only understandable and attractive for different target groups but also adequate from a scientific point of view, a two-fold feedback mechanism involving experts in the field was implemented. Through a close collaboration with a network of communication specialists from other institutions supporting the release, additional feedback and exchange of knowledge was enabled. Our insights, gained as part of the release process, can support others in developing user-centred products reviewed by experts in the field to inform about hazard and risk models.

Description: Distributed Acoustic Sensing (DAS) is becoming a powerful tool for earthquake monitoring, providing continuous strain-rate records of seismic events along fiber optic cables. However, the use of standard seismological techniques for earthquake source characterization requires the conversion of data in ground motion quantities. In this study we provide a new formulation for far-field strain radiation emitted by a seismic rupture, which allows to directly analyze DAS data in their native physical quantity. This formulation naturally accounts for the complex directional sensitivity of the fiber to body waves and to the shallow layering beneath the cable. In this domain, we show that the spectral amplitude of the strain integral is related to the Fourier transform of the source time function, and its modeling allows to determine the source parameters. We demonstrate the validity of the technique on two case-studies, where source parameters are consistent with estimates from standard seismic instruments in magnitude range 2.0–4.3. When analyzing events from a 1-month DAS survey in Chile, moment-corner frequency distribution shows scale invariant stress drop estimates, with an average of Δσ = (0.8 ± 0.6) MPa. Analysis of DAS data acquired in the Southern Apennines shows a dominance of the local attenuation that masks the effective corner frequency of the events. After estimating the local attenuation coefficient, we were able to retrieve the corner frequencies for the largest magnitude events in the catalog. Overall, this approach shows the capability of DAS technology to depict the characteristic scales of seismic sources and the released moment.

Description: The bacterial strains Brochothrix thermosphacta DH-B18 and Rathayibacter sp. DH-RSZ4 were isolated from raw sausage and escalope samples and grown in a CO2-rich modified atmosphere. Here, we present both circular genomes obtained by nanopore sequencing.

Description: Digital infrastructures have become indispensable in the field of modern research and science. These technological frameworks play a crucial role for the entire research cycle, supporting literature searches, aiding in data collection and analysis, facilitating the creation and publication of scholarly works, and ensuring the thorough documentation and long-term storage of research findings. Additionally, these infrastructures serve as a vital means for networking and communication among peers, creating the essential foundation of an open and transparent science and research ecosystem. Helmholtz employees were invited to join the Helmholtz Open Science Forum "Towards Open Digital Research Ecosystems - Interconnection Infrastructures" on February 14, 2024, where options for the seamless integration of these digital infrastructures have been discussed. Speakers presented insights into diverse efforts to the provision of open infrastructure structures and how their interconnection offers new possibilities for seamless and integrated workflows within the increasingly digitized research. Further, it was examined how such an integrated ecosystem can support open science practices and vice versa

Description: On January 22, 2024 the Helmholtz Open Science Office hosted the Second Helmholtz Open Science Forum on the topic of Open Science and Transfer. The online event addressed various aspects and issues around the interplay of Open Science, Technology Transfer, Knowledge Transfer and Citizen Science at the Helmholtz Association. Together with the participants important overlaps were identified against the backdrop of the digital transformation. Open science as a standard for scientific work creates foundations for successful transfer - and both topics can complement each other very well. The Helmholtz-internal event gave insights into current projects and initiatives relating to transfer to society, business and industry. Moreover, Helmholtz initiatives for Citizen Science and the successful practical implementation of Open Hardware were presented. The event also offered opportunities for networking and the exchange of ideas. This report documents the Second Helmholtz Open Science Forum on Open Science and Transfer.

Description: This study conducts mineralogical and chemical investigations on the oldest achondrite, Erg Chech 002 (∼4565 million yr old). This meteorite exhibits a disequilibrium igneous texture characterized by high-Mg-number (atomic Mg/(Mg + Fe2+)) orthopyroxene xenocrysts (Mg number = 60–80) embedded in an andesitic groundmass. Our research reveals that these xenocrysts were early formed crystals, loosely accumulated or scattered in the short-period magma ocean on the parent body. Subsequently, these crystals underwent agitation due to the influx of external materials. The assimilation of these materials enriched the 16O component of the magma ocean and induced a relatively reduced state. Furthermore, this process significantly cooled the magma ocean and inhibited the evaporation of alkali elements, leading to elevated concentrations of Na and K within the meteorite. Our findings suggest that the introduced materials are probably sourced from the reservoirs of CR clan meteorites, indicating extensive transport and mixing of materials within the early solar system.

Description: Continental crust at temperatures 〉 400 °C and depths 〉 10–20 km normally deforms in a ductile manner, but can become brittle and permeable in response to changes in temperature or stress state induced by fluid injection. In this study, we quantify the theoretical power generation potential of an enhanced geothermal system (EGS) at 15–17 km depth using a numerical model considering the dynamic response of the rock to injection-induced pressurization and cooling. Our simulations suggest that an EGS circulating 80 kg s−1 of water through initially 425 ℃ hot rock can produce thermal energy at a rate of ~ 120 MWth (~ 20 MWe) for up to two decades. As the fluid temperature decreases (less than 400 ℃), the corresponding thermal energy output decreases to around 40 MWth after a century of fluid circulation. However, exploiting these resources requires that temporal embrittlement of nominally ductile rock achieves bulk permeability values of ~ 10–15–10–14 m2 in a volume of rock with dimensions ~ 0.1 km3, as lower permeabilities result in unreasonably high injection pressures and higher permeabilities accelerate thermal drawdown. After cooling of the reservoir, the model assumes that the rock behaves in a brittle manner, which may lead to decreased fluid pressures due to a lowering of thresholds for failure in a critically stressed crust. However, such an evolution may also increase the risk for short-circuiting of fluid pathways, as in regular EGS systems. Although our theoretical investigation sheds light on the roles of geologic and operational parameters, realizing the potential of the ductile crust as an energy source requires cost-effective deep drilling technology as well as further research describing rock behavior at elevated temperatures and pressures.

Description: Primary granitic melt inclusions are trapped in garnets of eclogites in the garnet peridotite body of Pfaffenberg, Granulitgebirge (Bohemian Massif, Germany). These polycrystalline inclusions, based on their nature and composition, can be called nanogranitoids and contain mainly phlogopite/biotite, kumdykolite, quartz/rare cristobalite, a phase with the main Raman peak at 412 cm-1, a phase with the main Raman peak at 430 cm-1, osumilite and plagioclase. The melt is hydrous, peraluminous and granitic and significantly enriched in large ion lithophile elements (LILE), Th, U, Li, B and Pb. The melt major element composition resembles that of melts produced by the partial melting of metasediments, as also supported by its trace element signature characterized by elements (LILE, Pb, Li and B) typical of the continental crust. These microstructural and geochemical features suggest that the investigated melt originated in the subducted continental crust and interacted with the mantle to produce the Pfaffenberg eclogite. Moreover, in situ analyses and calculations based on partition coefficients between apatite and melt show that the melt was also enriched in Cl and F, pointing toward the presence of a brine during melting. The melt preserved in inclusions can thus be regarded as an example of a metasomatizing agent present at depth and responsible for the interaction between the crust and the mantle. Chemical similarities between this melt and other metasomatizing melts measured in other eclogites from the Granulitgebirge and Erzgebirge, in addition to the overall similar enrichment in trace elements observed in other metasomatized mantle rocks from central Europe, suggest an extended crustal contamination of the mantle beneath the Bohemian Massif during the Variscan orogeny.

Description: This dataset contains 92 estimates from individual studies for groundwater recharge rates on the Arabian Peninsula. Following information is sorted for each study: Location information: Country, Latitude*, Longitude* Estimated groundwater recharge rate: Representative value, Lower/Upper estimate range (all in mm/yr) Estimating methods** Scale of study: Aquifer scale, Study period, Study years Credibility***: Confidence, Confidence criteria (*) Location information was set as the middle point of the study area, in case that spatial coordinates are not given by the authors. (**) If more than 1 methods were used for the estimation, additional methods were written in "Method_2" (***) Confidence of estimates was evaluated by the same criteria used in another meta-study for the African continent (MacDonald et al. 2021; https://doi.org/10.1088/1748-9326/abd661) This dataset has been used to train the neural network model targeting global-scale estimation of groundwater recharge rate together with datasets used in other meta-studies. More detailed information is provided in the paper "Can eXplainable AI offer a new perspective for groundwater recharge estimation? – Global-scale modeling using neural network“.

Description: The fifth meeting of this series took place 5-6 February 2024 at the Helmholtz-Centre for Environmental Research – UFZ. A key topic of this meeting was the formation of Open Source Program Offices. The term is used to describe a concerted effort for an overarching structure in research organizations where research software engineers, computer departments, research data management units, technology/knowledge transfer units, libraries and legal departments cooperatively aim to build an environment conductive to excellent research, excellent research software, open science and technology/knowledge transfer.

Description: Flood-prone people and decision-makers are often unwilling to discuss and prepare for exceptional events, as such events are hard to perceive and out of experience for most people. Once an exceptional flood occurs, affected people and decision-makers are able to learn from this event. However, this learning is often focussed narrowly on the specific disaster experienced, thus missing an opportunity to explore and prepare for even more severe, or different, events. We propose spatial counterfactual floods as a means to motivate society to discuss exceptional events and suitable risk management strategies. We generate a set of extreme floods across Germany by shifting observed rainfall events in space and then propagating these shifted fields through a flood model. We argue that the storm tracks that caused past floods could have developed several tens of km away from the actual tracks. The set of spatial counterfactual floods generated contains events which are more than twice as severe as the most disastrous flood since 1950 in Germany. Moreover, regions that have been spared from havoc in the past should not feel safe, as they could have been badly hit as well. We propose spatial counterfactuals as a suitable approach to overcome society's unwillingness to think about and prepare for exceptional floods expected to occur more frequently in a warmer world.

Description: Excess heat (i.e., Urban Heat Island; UHI) and other urban conditions affect tree physiology with outcomes from enhanced growth to mortality. Resilient urban forests in the face of climate change require species-specific understanding of growth responses. However, previous studies assessing growth dynamics were primarily based on remote sensing of communities rather than individuals, or relied on labor-intensive methods that can limit the spatial coverage necessary to account for highly variable urban growing conditions. Here, we analyze growth dynamics of common urban street tree species over time and across space for Berlin (Germany) combining dendroecological (temporal) and inventory assessments (spatial). First, we show annual increments increased across the 20th century for early (i.e., young) growth. Second, we use an approach relying on open inventory data to identify growth potential in relation to excess heat while accounting for age, potential management effects, and the urban fabric (i.e., planting area; building density, height; available soil nutrients) with generalized additive models for the ten most abundant species. Our analyses showed that younger trees may benefit from increased temperatures, while older individuals feature lower growth at greater UHI magnitudes. Furthermore, planting area as well as building density modulate growth responses to temperature. Lastly, we discuss management implications in the context of climate change mitigation, considering that younger trees are predominantly located at UHI “hot spots” and will undergo the observed age-dependent shift in temperature-growth sensitivity. By relying on increasingly available open data, our approach here is or will be transferable to other urban regions.

Description: We report systematic first-principles results of structural properties and compression behavior based on density functional theory (DFT) and an exchange-correlation functional for solids, of Al-bearing garnets of general compositions in the pyrope-almandine-grossular solid solution system. The combination of DFT and a simple solid solution model is able to produce a compositional dependence of the compression curve consistent with trends observed in experimental studies. Using end-member properties extrapolated from our computations and perturbing an extant thermodynamic model we observe only marginal effects on the bulk sound velocity of pyrolite and MORB along relevant geothermal paths. However, this could hide important effects on the elemental partitioning between garnet and other major phases which should be further investigated both experimentally and computationally. We also present simulations of the effect of combined Fe and Ca substitutions for Mg on the elastic tensor of Al-bearing garnets, our simplified modeling shows only partial agreement with the trends observed in experiments. Therefore, further computational investigations, especially of the effect of Fe-Mg substitution on the tensor, are needed.

Description: Rift-Rift-Rift triple junctions are regions where three plates interact, generating complex networks of variably oriented faults. While the geometry of the fault networks is easily constrained from their surface expression, what remains unclear is how the kinematics of faults and their interactions vary spatially, and how these relate to the unusual crustal motions that result from three plates diverging from each other. The Afar depression lies at the triple junction between the African, Arabian, and Somalian plates (in the Horn of Africa), where the unique combination of observational data from structural mapping, seismicity, and Global Navigation Satellite System (GNSS) allows us to understand the link between fault kinematics and plate motions. We complement these observations with an analog model to gain insights into how the patterns and directions of faults relate to overall plate motions. A key finding in both the model and nature is that some adjacent normal faults form at high angles and generate T-shaped structures. These purely normal faults are synchronously active, which means that the extension direction varies ∼90° locally. These kinematic contrasts in our model and in nature occur despite the relatively smooth pattern of overall surface motions. The results indicate that normal faults interacting at high angles to form the T-shaped structures can evolve synchronously within a stress field that varies gently in magnitude but dramatically in orientation over a few kilometers.

Description: Rift propagation is a 3D thermo-mechanical process that often precedes continental breakup. Pre-existing microcontinental blocks and the associated lithospheric strength heterogeneities influence the style of rift propagation. Interestingly, some rifts propagate into pre-existing blocks and eventually cut through them (e.g., the Zhongsha Block and the Reed Bank), while others bypass these microcontinental blocks forming distinct overlapping rift branches (e.g., the East African Rift System). In this study, we use 3D numerical models to investigate the interaction between microcontinental blocks and rift propagation under different far-field extension rates. In doing so, we assess the impact of mantle lithospheric thicknesses and lower crustal rheology on the style of rift propagation. Our models reproduce the two types of rift propagation, characterized by propagating rifts that either split or bypass the pre-existing microcontinental blocks. We find that lithospheric thickness exerts dominant control, while lower crustal rheology of microcontinental blocks and the extension rate have less effect on rift propagation. Our model results can explain rift propagation patterns, block rotation and strong lithospheric thinning in the South China Sea, the East African Rift System, and the Woodlark Basin.

Description: Opportunistic sensors are increasingly used for rainfall measurement. However, their raw data are collected by a variety of systems that are often not primarily intended for rainfall monitoring, resulting in a plethora of different data formats and a lack of common standards. This hinders the sharing of opportunistic sensing (OS) data, their automated processing, and, at the end, their practical usage and integration into standard observation systems. This paper summarises the experiences of the more than 100 members of the OpenSense Cost Action involved in the OS of rainfall. We review the current practice of collecting and storing precipitation OS data and corresponding metadata, and propose new common guidelines describing the requirements on data and metadata collection, harmonising naming conventions, and defining human-readable and machine readable file formats for data and metadata storage. We focus on three sensors identified by the OpenSense community as prominent representatives of the OS of precipitation: Commercial microwave links (CML): fixed point-to-point radio links mainly used as backhauling connections in telecommunication networks Satellite microwave links (SML): radio links between geostationary Earth orbit (GEO) satellites and ground user terminals. Personal weather stations (PWS): non-professional meteorological sensors owned by citizens. The conventions presented in this paper are primarily designed for storing, handling, and sharing historical time series and do not consider specific requirements for using OS data in real time for operational purposes. The conventions are already now accepted by the ever growing OpenSense community and represent an important step towards automated processing of OS raw data and community development of joint OS software packages.

Description: Pyrite is a ubiquitous sulfide mineral found in diverse geological settings and holds great significance in the formation of Au deposits as well as the safe utilization of groundwater due to its remarkable ability to incorporate substantial amounts of As. However, despite its importance, there remains a dearth of fundamental data on the partitioning of As between pyrite and fluid, which is key for accurately modeling the As distribution in these environments. Here, we present new insights into the partitioning behavior of As between pyrite and fluid at conditions that mimic natural fluid systems. Pyrite was synthesized by replacement of natural siderite in hydrothermal experiments at 200 °C and pH 5 applying a wide range of fluid As concentrations, spanning from 0.001 to 100 µg/g. The As distribution and concentration in synthetic pyrite was analyzed by quantitative LA-ICP-MS mapping providing a high spatial resolution and sensitivity at 2–3 µm image pixel size at a detection limit of ∼1 µg/g at the single pixel scale. Pyrite-fluid partitioning coefficients (DAs(py/fluid)) between synthetic pyrite and experimental fluid agree with previously published data for high fluid As concentrations of 1 µg/g to 100 µg/g (DAs 〈 2000). However, at low As concentrations in the experimental fluid (〈1 µg/g), a steep increase in the DAs(py/fluid) values of up to ∼30,000 was detected, demonstrating even stronger As partitioning into pyrite. This is confirmed by the analyses of natural pyrite that precipitated from As-poor fluids (0.3–0.4 ng/g) within a deep anoxic aquifer in SE Sweden. The discovery holds significant implications for the mobility and scavenging of As, which in turn is important for understanding the formation and fingerprinting of mineral deposits as well as for the secure utilization of groundwater resources.

Description: Important features of Sn mineralization are the heterogeneous geographic distribution and frequent regional separation from W mineralization in spite some similarities of Sn and W behavior during magmatic processes. Major Sn and W mineralization is often spatially associated with peraluminous granites, which are derived from partial melting of metasediments. Several concepts have been suggested to explain those features, such as a weathering-related Sn-enriched source, Sn redistribution between melts and restite during protolith melting, and extensive fractional crystallization. We demonstrate the importance of protolith composition for the formation of Sn (and W) granites by using a comprehensive bulk-rock composition dataset from Precambrian metasediments of the South China Sn-W province and employing a thermodynamic modeling approach. We used four compositional proxies for phase equilibria calculations, which are the metasediments of the Mengdong, Sibao, Pingbian, and Shuangqiaoshan Groups. It is well documented that those Precambrian metasediments are important protoliths of Sn granites in South China. We present quantitative evaluation of the control of protolith composition in the generation of Sn-enriched granitic melts using South China as example, but our conclusions may also be applicable to worldwide Sn–enriched granites. Our results indicate that the protolith major-element geochemistry controls the anatectic reactions and melt productivity at specific melting conditions, and consequently the partitioning behavior of Sn. Further, pre-enrichment of Sn is crucial to the fertility of granitic melt and may be a prerequisite, particularly for the formation of giant Sn deposits. We propose that the heterogeneous distribution of favorable source rocks is one of the important factors that control the spatial distribution of major Sn (and W) districts in South China and other regions worldwide.

Description: The evolution of topography in forearc regions results from the complex interplay of crustal and mantle processes. The Southern Apennines represent a well‐studied forearc region that experienced several tectonic phases, initially marked by compressional deformation followed by extension and large‐scale uplift. We present a new structural, geomorphic and fluvial analysis of the Pollino Massif and surrounding intermontane basins (Mercure, Campotenese and Castrovillari) to unravel their evolution since the Pliocene. We constrain multiple tectonic transport directions, evolution of the drainage, and magnitude and timing of long‐ term incision following base level falls. Two sets of knickpoints suggest two phases of base level lowering and allow to estimate ∼500 m of long‐term uplift (late Pleistocene), as observed in the Sila Massif. On a smaller spatial scale, the evolution and formation of topographic relief, sedimentation, and opening of intermontane basins is strongly controlled by the recent increase in rock uplift rate and fault activity. At the regional scale, an along‐strike, long‐wavelength uplift pattern from north to south can be explained by progressive lateral slab tearing and inflow of asthenospheric mantle beneath Pollino and Sila, which in turn may have promoted extensional tectonics. The lower uplift of Le Serre Massif may be explained as result of weak plate coupling due to narrowing of the Calabrian slab. The onset of uplift in the Pollino Massif, ranging from 400 to 800 ka, is consistent with that one proposed in the southern Calabrian forearc, suggesting a possible synchronism of uplift, and lateral tearing of the Calabrian slab.

Description: We address the possible link between the age of subducting oceanic lithosphere and growth of elevated cordilleras versus extension-dominated arc regions. Singularity exists in South America: the lowest elevated Andean segments are found in Patagonia where the active Chile Ridge enters the trench. Subduction of active ridge triggers thermal doming, crustal extension and attenuation of former cordilleras. At the Antarctica–South America connection, three active ridge subductions induced the disruption of a former continuous cordillera during the opening of Drake Passage. Active ridge subduction induces lithosphere thermal erosion and related crustal extension in the upper plate. Evolution of regions worldwide experiencing ridge subduction confirms this hypothesis.

Description: Der vorliegende Kolloquiumsband der ”Blauen Bände“ umfasst 41 Beiträge. Dem diesjährigen Kolloquiumsband wurde ein neues Format hinzugefügt. Der sogenannte 1-Pager bietet die Möglichkeit eines einseitigen Beitrages, bestehend aus einer Zusammenfassung, einer Abbildung und einer Referenz zu weiterführender Information. In der Vergangenheit haben sich immer wieder Diskussionen zur Publikation von wissenschaftlichen Arbeiten in den ”Blauen Bänden“ ergeben. Mit diesem neuen niederschwelligen 1-Pager Format wollen wir alle EMTF Kolloquiums Teilnehmer*innen motivieren, ihren Beitrag in diesem traditionellen Kolloquiumsband sichtbar zu machen. Das ist mit einer vergleichsweise großen Anzahl von Beitr¨agen gut gelungen. In diesem Band sind Extended Abstracts (4), Vortragsfolien (8), Postern (19) sowie der neu eingeführte 1-Pager (10) zusammengefasst. Als weitere Neuerung wurde dem Band erstmalig das Tagungsprogramm als Anhang beigefügt. Wir danken allen Teilnehmern, die mit ihren Beiträgen dazu beitragen, die Breite des Kolloquiums auch in den ”Blauen Bänden“ zu dokumentieren.

Description: The Zagros Fold and Thrust Belt (ZFTB) is an outstanding orogen running from eastern Turkey to the Makran area. It is formed as a consequence of the convergence between the Arabian and the Eurasian plates that occurred in the Neogene. This still active and long-lasting process generated a topographic configuration dominated by a series of parallel folding structures which, at places, isolate internal basins. The topographic configuration has, in turn, profoundly influenced the river network evolution, which follows a trellis pattern with the main valleys developed in the synclines and rivers that occasionally cut into anticlines. The peculiar climate, characterised by arid and semi-arid conditions, makes most of the rivers ephemeral, alimented only by short rainfall events. For this reason, the sediments are transported over short distances and deposited in huge alluvial fans. Although the Zagros is one of the most studied belts in the world, its tectonic evolution is far from being fully understood. Debated, for example, are the beginning of collision, the primary deformation mechanism, the evolution of the drainage system, the formation process of the alluvial fans, and the interrelations between landscape, tectonics, and climate. This paper, focusing on the geodynamic, geological, stratigraphic, and topographic configuration of the Zagros belt, is intended to be a compendium of the most up-to-date knowledge on the Zagros and aims to provide the cognitive basis for future research that can find answers to outstanding questions.

Description: Slow slip events (SSEs) have been observed in spatial and temporal proximity to megathrust earthquakes in various subduction zones, including the 2014 Mw 7.3 Guerrero, Mexico earthquake which was preceded by a Mw 7.6 SSE. However, the underlying physics connecting SSEs to earthquakes remains elusive. Here, we link 3D slow‐slip cycle models with dynamic rupture simulations across the geometrically complex flat‐slab Cocos plate boundary. Our physics‐based models reproduce key regional geodetic and teleseismic fault slip observations on timescales from decades to seconds. We find that accelerating SSE fronts transiently increase shear stress at the down‐dip end of the seismogenic zone, modulated by the complex geometry beneath the Guerrero segment. The shear stresses cast by the migrating fronts of the 2014 Mw 7.6 SSE are significantly larger than those during the three previous episodic SSEs that occurred along the same portion of the megathrust. We show that the SSE transient stresses are large enough to nucleate earthquake dynamic rupture and affect rupture dynamics. However, additional frictional asperities in the seismogenic part of the megathrust are required to explain the observed complexities in the coseismic energy release and static surface displacements of the Guerrero earthquake. We conclude that it is crucial to jointly analyze the long‐ and shortterm interactions and complexities of SSEs and megathrust earthquakes across several (a)seismic cycles accounting for megathrust geometry. Our study has important implications for identifying earthquake precursors and understanding the link between transient and sudden megathrust faulting processes.

Description: Fault‐damage zones comprise multiscale fracture networks that may slip dynamically and interact with the main fault during earthquake rupture. Using 3D dynamic rupture simulations and scale‐dependent fracture energy, we examine dynamic interactions of more than 800 intersecting multiscale fractures surrounding a listric fault, emulating a major listric fault and its damage zone. We investigate 10 distinct orientations of maximum horizontal stress, probing the conditions necessary for sustained slip within the fracture network or activating the main fault. Additionally, we assess the feasibility of nucleating dynamic rupture earthquake cascades from a distant fracture and investigate the sensitivity of fracture network cascading rupture to the effective normal stress level. We model either pure cascades or main fault rupture with limited offfault slip. We find that cascading ruptures within the fracture network are dynamically feasible under certain conditions, including: (a) the fracture energy scales with fracture and fault size, (b) favorable relative pre‐stress of fractures within the ambient stress field, and (c) close proximity of fractures. We find that cascading rupture within the fracture network discourages rupture on the main fault. Our simulations suggest that fractures with favorable relative pre‐stress, embedded within a fault damage zone, may lead to cascading earthquake rupture that shadows main fault slip. We find that such off‐fault events may reach moment magnitudes up to Mw ≈ 5.5, comparable to magnitudes that can be otherwise hosted by the main fault. Our findings offer insights into physical processes governing cascading earthquake dynamic rupture within multiscale fracture networks.

Description: Wetlands in Arctic drained lake basins (DLBs) have a high potential for carbon storage in vegetation and peat as well as for elevated greenhouse gas emissions. However, the evolution of vegetation and organic matter is rarely studied in DLBs, making these abundant wetlands especially uncertain elements of the permafrost carbon budget. We surveyed multiple DLB generations in northern Alaska with the goal to assess vegetation, microtopography, and organic matter in surface sediment and pond water in DLBs and to provide the first high-resolution land cover classification for a DLB system focussing on moisture-related vegetation classes for the Teshekpuk Lake region. We associated sediment properties and methane concentrations along a post-drainage succession gradient with remote sensing-derived land cover classes. Our study distinguished five eco-hydrological classes using statistical clustering of vegetation data, which corresponded to the land cover classes. We identified surface wetness and time since drainage as predictors of vegetation composition. Microtopographic complexity increased after drainage. Organic carbon and nitrogen contents in sediment, and dissolved organic carbon (DOC) and dissolved nitrogen (DN) in ponds were high throughout, indicating high organic matter availability and decomposition. We confirmed wetness as a predictor of sediment methane concentrations. Our findings suggest moderate to high methane concentrations independent of drainage age, with particularly high concentrations beneath submerged patches (up to 200 μmol l−1) and in pond water (up to 22 μmol l−1). In our DLB system, wet and shallow submerged patches with high methane concentrations occupied 54% of the area, and ponds with high DOC, DN and methane occupied another 11%. In conclusion, we demonstrate that DLB wetlands are highly productive regarding organic matter decomposition and methane production. Machine learning-aided land cover classification using high-resolution multispectral satellite imagery provides a useful tool for future upscaling of sediment properties and methane emission potentials from Arctic DLBs.

Description: We present a seismic catalog (Bindi et al., 2024, https://doi.org/10.5880/GFZ.2.6.2023.010) including energy magnitude Me estimated from P waves recorded at teleseismic distances in the range 20° 1 98° and for depths shorter than 80 km. The catalog is built starting from the event catalog disseminated by GEOFON (GEOFOrschungsNetz), considering 6349 earthquakes with moment magnitude Mw 5 occurring between 2011 and 2023. Magnitudes are computed using 1 031 396 freely available waveforms archived in EIDA (European Integrated Data Archive) and IRIS (Incorporated Research Institutions for Seismology) repositories, retrieved through the standard International Federation of Digital Seismograph Networks (FDSN) web services (https://www.fdsn.org/webservices/, last access: March 2024). A reduced, high-quality catalog for events with Mw 5〉_8 and from which stations and events with only few recordings were removed forms the basis of a detailed analysis of the residuals of individual station measurements, which are decomposed into station- and event-specific terms and a term accounting for remaining variability. The derived Me values are compared to Mw computed by GEOFON and with the Me values calculated by IRIS. Software and tools developed for downloading and processing waveforms for bulk analysis and an add-on for SeisComP for real-time assessment of Me in a monitoring context are also provided alongside the catalog. The SeisComP add-on has been part of the GEOFON routine processing since December 2021 to compute and disseminate Me for major events via the existing services.

Description: The West Siberian Seaway connected the Tethys to the Arctic Ocean in the Paleogene and played an important role for Eurasian-Arctic biogeography, ocean circulation, and climate. However, the paleogeography and geological mechanisms enabling the seaway are not well constrained, which complicates linking the seaway evolution to paleoenvironmental changes. Here, we investigate the paleogeography of the Peri-Tethys realms for the Cenozoic time (66–0 Ma), including the West Siberian Seaway, and quantify the influence of mantle convection and corresponding dynamic topography. We start by generating continuous digital elevation models for Eurasia, Arabia, and Northern Africa, by digitizing regional paleogeographic maps and additional geological information and incorporate them in a global paleogeography model with nominal million-year resolution. Then we compute time-dependent dynamic topography for the same time interval and find a clear correlation between changes in dynamic topography and the paleogeographic evolution of Central Eurasia and the West Siberian Seaway. Our results suggest that mantle convection played a greater role in Eurasian paleogeography than previously recognized. Mantle flow may have influenced oceanic connections between the Arctic and global ocean providing a link between deep mantle convection, surface evolution, and environmental changes. Our reconstructions also indicate that the Arctic Ocean may have been isolated from the global ocean in the Eocene, even if the West Siberian Seaway was open, as the Peri-Tethys – Tethys connection was limited, and the Greenland-Scotland Ridge was a landbridge.

Description: To reach their net-zero targets, countries will have to compensate hard-to-abate CO2 emissions through carbon dioxide removal (CDR). Yet, current assessments rarely include socio-cultural or institutional aspects or fail to contextualize CDR options for implementation. Here we present a context-specific feasibility assessment of CDR options for the example of Germany. We assess 14 CDR options, including three chemical carbon capture options, six options for bioenergy combined with carbon capture and storage (BECCS), and five options that aim to increase ecosystem carbon uptake. The assessment addresses technological, economic, environmental, institutional, social-cultural and systemic considerations using a traffic-light system to evaluate implementation opportunities and hurdles. We find that in Germany CDR options like cover crops or seagrass restoration currently face comparably low implementation hurdles in terms of technological, economic, or environmental feasibility and low institutional or social opposition but show comparably small CO2 removal potentials. In contrast, some BECCS options that show high CDR potentials face significant techno-economic, societal and institutional hurdles when it comes to the geological storage of CO2. While a combination of CDR options is likely required to meet the net-zero target in Germany, the current climate protection law includes a limited set of options. Our analysis aims to provide comprehensive information on CDR hurdles and possibilities for Germany for use in further research on CDR options, climate, and energy scenario development, as well as an effective decision support basis for various actors.

Description: In this study, we investigated four different ferroan dolomite samples using in situ Raman spectroscopy and powder X-ray diffraction (XRD) at high pressures up to 48 GPa and at room temperature. Our results show that the transition from Dolomite-I (Dol-I) to Dolomite-II (Dol-II) occurs above 13-16 GPa, and the transition pressure depends on the composition of the solid solution. Compression above 32-35 GPa results in the appearance of the Dolomite-IIIc (Dol-IIIc) or Dolomite-IIIb (Dol-IIIb). In the high-pressure XRD study, we found that the XRD patterns of the Ca0.97(Mg0.77,Fe0.23Mn0.03)(CO3)2 (xFe = 0.23, Ank23) can be indexed as Dol-IIIc at 44 GPa, while the rhombohedral Dol-IIIb structure matches better with the XRD patterns of the xFe = 0.40 (Ank40) and 0.64 (Ank64) solid solutions. Additionally, in the Raman spectra of the Fe-richest sample (Ca0.99(Mg0.33Fe0.64Mn0.05) (CO3)2 (Ank64), we observed an abrupt frequency downshift of the CO3-stretching vibrations between 40 – 42 GPa, which may reflect a pressure-induced Fe2+ spin transition. We further investigated two samples with xFe= 0.19 (Ank19) and 0.23 (Ank23) at high pressure and high temperatures, up to at least 2600 K. The experiments revealed that the unquenchable Dol-IIIc structure could be a stable high-pressure/high-temperature polymorph in ferroan dolomite up to at least 2600 K.

Description: Global climate warming is accelerating permafrost degradation. The large amounts of soil organic matter in permafrost-affected soils are prone to increased microbial decomposition in a warming climate. Along with permafrost degradation, changes to the soil microbiome play a crucial role in enhancing our understanding and in predicting the feedback of permafrost carbon. In this article, we review the current state of knowledge of carbon-cycling microbial ecology in permafrost regions. Microbiomes in degrading permafrost exhibit variations across spatial and temporal scales. Among the short-term, rapid degradation scenarios, thermokarst lakes have distinct biogeochemical conditions promoting emission of greenhouse gases. Additionally, extreme climatic events can trigger drastic changes in microbial consortia and activity. Notably, environmental conditions appear to exert a dominant influence on microbial assembly in permafrost ecosystems. Furthermore, as the global climate is closely connected to various permafrost regions, it will be crucial to extend our understanding beyond local scales, for example by conducting comparative and integrative studies between Arctic permafrost and alpine permafrost on the Qinghai–Tibet Plateau at global and continental scales. These comparative studies will enhance our understanding of microbial functioning in degrading permafrost ecosystems and help inform effective strategies for managing and mitigating the impacts of climate change on permafrost regions.

Description: Geomechanics has a marked impact on the safe and sustainable use of the subsurface. This Special Publication contains contributions detailing the latest efforts in present-day in-situ stress characterization, prediction and modelling from the borehole to plate-tectonic scale. A particular emphasis is on the uncertainties that are often associated with geomechanics.

Description: Desert environments constitute one of the largest and yet most fragile ecosystems on Earth. Under the absence of regular precipitation, microorganisms are the main ecological component mediating nutrient fluxes by using soil components, like minerals and salts, and atmospheric gases as a source for energy and water. While most of the previous studies on microbial ecology of desert environments have focused on surface environments, little is known about microbial life in deeper sediment layers. Our study is extending the limited knowledge about microbial communities within the deeper subsurface of the hyperarid core of the Atacama Desert. By employing intracellular DNA extraction and subsequent 16S rRNA sequencing of samples collected from a soil pit in the Yungay region of the Atacama Desert, we unveiled a potentially viable microbial subsurface community residing at depths down to 4.20 m. In the upper 80 cm of the playa sediments, microbial communities were dominated by Firmicutes taxa showing a depth-related decrease in biomass correlating with increasing amounts of soluble salts. High salt concentrations are possibly causing microbial colonization to cease in the lower part of the playa sediments between 80 and 200 cm depth. In the underlying alluvial fan deposits, microbial communities reemerge, possibly due to gypsum providing an alternative water source. The discovery of this deeper subsurface community is reshaping our understanding of desert soils, emphasizing the need to consider subsurface environments in future explorations of arid ecosystems.

Description: The current crustal stress field is of key importance to understand geodynamic processes and to assess stability aspects during subsurface usage. To provide a 3-D continuous description of the stress state, linear elastic forward geomechanical-numerical models are used. These models solve the equilibrium of forces between gravitational volume forces and surfaces forces im- posed mainly by plate tectonics. The latter are responsible for the horizontal stress anisotropy and impose the inverse problem to estimate horizontal displacement boundary conditions that provide a fit best to horizontal stress magnitude data within the model volume. Ho wever , horizontal stress magnitude data have high uncertainties and they are sparse, clustered and not necessaril y representati ve for a larger rock v olume. Even w hen Bay esian statistics are incor - porated and additional stress information such as borehole failure observations or formation integrity test are used to further constrain the solution space, this approach may result in a low accuracy of the model results, that is the result is not correct. Here, we present an alternative approach that removes the dependence of the solution space based on stress magnitude data to avoid potential low accuracy . Initially , a solution space that contains all stress states that are physically reasonable is defined. Stress magnitude data and the additional stress information are then used in a Bayesian framework to e v aluate which solutions are more likely than others. We first show and validate our approach with a generic truth model and then apply it to a case study of the Molasse foreland basin of the Alps in Southern Germany. The results show that the model’s ability to predict a reliable stress state is increasing while the number of likely solutions may also increase, and that outlier of stress magnitude data can be identified. This alternative approach results in a substantial increase in computational speed as we perform most of the calculations anal yticall y.

Description: Cyanobacteria are major contributors to algal blooms in inland waters, threatening ecosystem function and water uses, especially when toxin-producing strains dominate. Here, we examine 140 hyperspectral (HS) images of five representatives of the widespread, potentially toxin-producing and bloom-forming genera Microcystis, Planktothrix, Aphanizomenon, Chrysosporum and Dolichospermum, to determine the potential of utilizing visible and near-infrared (VIS/NIR) reflectance for their discrimination. Cultures were grown under various light and nutrient conditions to induce a wide range of pigment and spectral variability, mimicking variations potentially found in natural environments. Importantly, we assumed a simplified scenario where all spectral variability was derived from cyanobacteria. Throughout the cyanobacterial life cycle, multiple HS images were acquired along with extractions of chlorophyll a and phycocyanin. Images were calibrated and average spectra from the region of interest were extracted using k-means algorithm. The spectral data were pre-processed with seven methods for subsequent integration into Random Forest models, whose performances were evaluated with different metrics on the training, validation and testing sets. Successful classification rates close to 90 % were achieved using either the first or second derivative along with spectral smoothing, identifying important wavelengths in both the VIS and NIR. Microcystis and Chrysosporum were the genera achieving the highest accuracy (〉95 %), followed by Planktothrix (79 %), and finally Dolichospermum and Aphanizomenon (〉50 %). The potential of HS imagery to discriminate among toxic cyanobacteria is discussed in the context of advanced monitoring, aiming to enhance remote sensing capabilities and risk predictions for water bodies affected by cyanobacterial harmful algal blooms.

Description: This archive disseminated through the GFZ-Data Service includes both results and information as-sociated to Bindi et al. (2023). In particular, the archive includes a seismic catalogue reporting ener-gy magnitude Me estimated form vertical P-waves recorded at teleseismic distances in the range 20°≤ D ≤ 98°, following Di Giacomo et al (2008, 2010). The catalogue is built considering 6349 earth-quakes included in the GEOFON (Quinteros et al, 2021) catalogue with moment magnitude Mw larger than 5 and occurring after 2011. Tools used to compute the energy magnitude are free available. In particular, we used stream2segment (Zaccarelli, 2018) to download data from IRIS (https://ds.iris.edu/ds) and EIDA (Strollo et al., 2021) repositories, and me-compute [Zaccarelli, 2023) to process waveforms and compute Me. The methodology applied to me-compute is also implemented as add-on for SeicomP (GFZ and Gempa, 2020) in order to allow the real time computation of Me (https://github.com/SeisComP/scmert).

Description: The EU funded project CRM-geothermal aims to establish an overview of the potential for critical raw materials (CRM) in geothermal fluids across the EU and third countries (Ref). Within this framework, the geothermal sites of Tuzla, Seferihisar and Dikili in eastern Turkey have been visited in March 2023. To estimate the potential of CRM at these sites, a comprehensive sampling program was performed. Rock samples (drill gravel) of the production borehole and scaling from gas-water separators were obtained. Furthermore, sampling of geothermal fluids (gas and brine) and precipitates (salt) along the production line was performed. Here, the results of the geochemical analyses of solid sample materials (drill gravel, scales and salt) are presented. All analyses were performed in the ElMiE-Lab (Elements and Minerals of the Earth Laboratory) at German Research Centre for Geosciences Potsdam, Germany (https://labinfrastructure.geo-x.net/laboratories/8). For their major and minor element compositions, bulk samples of drill gravel and scales were analyzed with XRF and ICP-MS, respectively. Salt precipitates were analyzed for dry loss and mineral composition using XRD

Description: It is well known that triaxial deformation is a common feature of continental tectonics, and is accommodated by complex polymodal fault networks. Field investigations confirm that multiple phases involving time-dependent three-dimensional strain conditions (e.g. constriction, plane, and flattening strain) affect the spatial and temporal interaction of polymodal fault systems. However, a key question remains: How do changing strain conditions affect the reactivation of fault systems that formed during a previous deformation phase? Here, we conduct scaled analogue models with time-dependent boundary conditions to investigate how fault networks evolve under changing boundary conditions, including reactivation and formation of new faults. We have developed a setup in which a basal rubber sheet is stretched in one direction, so that longitudinal extension and layer thinning are accompanied by lateral shortening, hence producing triaxial deformation (Liu et al. in revision). According to previous brittle-viscous experiments with this set-up, an increase in longitudinal extension velocity results in a higher coupling between the rubber base and brittle layer, generating increasing transmission of lateral shortening from the base into the brittle layer. We thus induce constriction-to-plane strain conditions in the brittle layer as a function of longitudinal extension velocity by varying the magnitude of lateral contraction. In a new set of experiments, by varying extension velocity either stepwise or continuously, we realize time-dependent kinematic boundary conditions including deformation phases and secular changes, respectively. Digital image correlation (DIC) and photogrammetry (structure from motion, SFM) are employed to track the 3D kinematic surface and topography evolution, respectively. Preliminary observations show both the formation of new faults and the reactivation of early phase faults through a change from plane to constriction strain. Conversely, a change from constriction to plane strain conditions results in the abandonment of the early phase fault network as it becomes overprinted by fault systems of the subsequent phase. Moreover, early-phase fault systems influence the propagation and linkage of fault populations in subsequent phases. Our analogue models highlight the impact of strain conditions on the overall plan-view geometry of fault populations, providing alternative explanations for complex fault patterns and interactions (e.g. the Jeanne d’Arc basin, the North Træna Basin, and the Beagle Platform).

Description: In autumn 2022, an expedition to Tanzania was undertaken within the framework of the research project “CRM-geothermal” and Scintific Priority Program (SPP) 2238 “Dynamics of Ore Metal Enrichment”. Within „CRM-geothermal“ we are looking for an environmentally friendly co-production of critical raw materials together with the provision of geothermal energy. In the EARS, high levels of rare earth elements (REE), Sr, Ba and Mg are expected in waters and solids in areas with alkaline volcanic rocks, while other critical elements, including helium, have been sought in other localities. In particular, the eastern branch is the most juvenile sector and has increased geothermal potential related to hot fluids migrating along permeable faults. Tanzania was crossed from north to south, along the eastern arm of the EARS, to collect gas, water, rock and sediment samples associated with natural hot springs, lakes and vents. On site, physical and chemical parameters were measured in-situ and documented together with the geology, infrastructure and the domestic use of the hot site. In the south, existing drill sites and geothermal development areas were visited and gas and water samples were taken from boreholes and rocks sampled from drill cores. The survey covered 13 sites, from Lake Natron in the north to Lake Malawi in the south (see map).

Description: Present day system Earth research utilizes the tool ‘Scientific Drilling’ to access samples and to monitor deep Earth processes that cannot be tackled by other scientific means. Unlike most laboratory experiments or computer modelling, drilling projects are massive field endeavours requiring intense collaboration of researchers with engineers and service providers. In the framework of the International Continental Scientific Drilling Program, ICDP, more than seventy drilling projects have been conducted, from multiyear big research programs to short, smallscale deployments such as lake drilling projects. ICDP has supported these projects not only through grants covering field-related costs, but also through a variety of scientific-technical services and support, as well as active help in data management, outreach and publication. These services are described in this booklet. Due to its instructional character, we call it the ICDP Primer.

Description: In this paper, we anticipate geospatial population distributions to quantify the future number of people living in earthquake-prone and tsunami-prone areas of Lima and Callao, Peru. We capitalize upon existing gridded population time series data sets, which are provided on an open-source basis globally, and implement machine learning models tailored for time series analysis, i.e., based on long short-term memory (LSTM) networks, for prediction of future time steps. Specifically, we harvest WorldPop population data and teach LSTM and convolutional LSTM models equipped with both unidirectional and bidirectional learning mechanisms, which are derived from different feature sets, i.e., driving factors. To gain insights regarding the competitive performance of LSTM-based models in this application context, we also implement multilinear regression and random forest models for comparison. The results clearly underline the value of the LSTM-based models for forecasting gridded population data; the most accurate prediction obtained with an LSTM equipped with a bidirectional learning scheme features a root-mean-squared error of 3.63 people per 100 × 100 m grid cell while maintaining an excellent model fit (R2= 0.995). We deploy this model for anticipation of population along a 3-year interval until the year 2035. Especially in areas of high peak ground acceleration of 207–210 cm s−2, the population is anticipated to experience growth of almost 30 % over the forecasted time span, which simultaneously corresponds to 70 % of the predicted additional inhabitants of Lima. The population in the tsunami inundation area is anticipated to grow by 61 % until 2035, which is substantially more than the average growth of 35 % for the city. Uncovering those relations can help urban planners and policymakers to develop effective risk mitigation strategies.

Description: Natural gas hydrates form at elevated pressure and low temperatures in the presence of sufficient quantities of gas and water and have therefore been discovered on all continental margins and in permafrost regions. In the marine hydrate-bearing sediments, gas hydrates, depending on their content, can transform a loose sediment into a consolidated rock with a strongly increased strength. In permafrost regions the hydrate stability zone can extent deep into the ice-bearing permafrost and, therefore, both, ice and hydrate can consolidate the sediment. However, the strength of methane hydrate is much higher than that of ice, which behaves much more ductile. Consequently, the resulting strength of a sediment, containing both components, strongly depends on the ice to hydrate ratio. Conversely, the decomposition of natural gas hydrates in marine or permafrost sediments leads to a reduction in the mechanical strength of the host sediment. In addition, the release of gas can create overpressure in the pore spaces, reducing the effective stress and leading to instabilities in the sediment structure. Since both continental margins and permafrost regions are used by humans for various activities that largely depend on the mechanical stability of the sediments, knowledge of the main factors and processes that determine the stability of weakly consolidated sediments is crucial. Both the thawing of ice and the decomposition of gas hydrates in permafrost soils lead to a change in the geo-mechanical properties of the host sediment. The residual and peak shear strengths of ice- and hydrate-bearing sediments were investigated using a ring shear cell developed at the GFZ. Based on literature data and our results, we discuss the dependence of the geo-mechanical properties of sediments on ice and hydrate saturation and the possible consequences if their proportion diminishes.

Description: Natural gas hydrates are crystalline compounds that are formed from hydrogen-bonded water molecules and gas molecules. They mainly contain climate-active CH4, but also other light hydrocarbons, CO2 or H2S They exhibit a high sensitivity to variations in temperature and pressure, mainly driven by environmental changes. The oceanic or atmospheric warming resulting from climate change may trigger the decompositions of hydrates, potentially releasing significant amounts of CH4. To assess the potential risks associated with CH4 release from destabilized hydrate deposits, a precise understanding of the dissociation behaviour of gas hydrates becomes crucial. In this study, a systematic investigation on the dissociation process of sI CH4 hydrates, sII CH4+C3H8 hydrates, and sII multi-component CH4+C2H6+C3H8+CO2 mixed hydrates was reported. We employed a combination of experimental and molecular dynamics (MD) simulations to provide a more nuanced understanding of the hydrate dissociation behaviours, which primarily shed light on the molecular aspects. The dissociation was induced through thermal stimulation to mimic climate warming. Both in situ and ex situ Raman spectroscopic measurements were performed continuously to characterize the hydrate phase. Throughout the dissociation process, hydrate composition, surface morphology, and the large-to-small cavity ratios were determined. MD simulations were carried out under similar conditions, providing advanced insights and perspectives that couldn't be readily extracted from experimental observations alone. Both experimental and simulation outcomes indicate that intrinsic kinetics likely govern the early stage of hydrate dissociation. A significant development in the dissociation process is the hindrance caused by the formation of a quasi-liquid or amorphous phase at the surface of the hydrate particles after the breakup of the outer layer of hydrate cavities. The unstable (partial) hydrate cavities that form within this quasi-liquid phase are oversaturated with gas molecules. Consequently, gas hydrates undergo a cycle of decomposition-reformation-continuing decomposition until the crystal eventually disappears. With decomposition dominating the process, both experimental and numerical simulation results demonstrate that the breakup of large cavities (51262) is faster than that of small ones (512) in sI hydrates. Conversely, a faster breakdown of small 512 cavities in sII hydrates is observed. Additionally, during the dissociation process of sII CH4-C3H8 hydrate, the cavities occupied by CH4 preferentially collapse compared to those filled with C3H8. Similarly, over the dissociation of multi-component hydrate, cavities filled with CH4 exhibit a preferential collapse compared to those filled with C3H8, C2H6, and CO2. These findings show the complexity and differences in the dissociation behavior of natural gas hydrates depending on their composition and structure and can therefore make an important contribution to an accurate assessment of CH4 release from destabilized hydrate deposits in response to climate change.

Description: This dataset is the result of an experimental series that was carried out in September/October 2022 at GFZ German Research Centre for Geosciences, Potsdam, Germany to observe biosorption of lead under extreme conditions. Synthetic solutions, simulating the geothermal fluids from the Heemskerk geothermal power plant were were prepared in 30 ml glass vials (Rotalibo screw neck ND24 EPA). To prepare the stock solutions, sodium chloride (NaCl, 99.8 %, Cellpure, Merck, DE) was added at 265 g/L and Pb(II), in form of lead nitrate (Pb(NO3 )2 , Merck, DE), at 1 g/L to ultrapure water. To assess the impact of acetic acid on lead biosorption, two treatments were done: one without acetic acid and one where acetic acid (100 %, Merck, DE) was added at 60 mg/L. Finally, dead biomass of the fungus Penicillium citrinum was added in the samples at a concentration of 4 g/L (Wahab et al., 2017). The samples were incubated in an autoclave at a pressure of 8 bars on a rotative shaker. The temperature was set at 25 °C, 60 °C or 98 °C with three contact times (1, 2 and 3 h). All treatments were performed in triplicates. For each treatment, two controls without biomass were done. Control samples without the addition of NaCl were done in duplicate, at 25 °C and for 2 h. After incubation, samples were filtered through a 0.22 µm nitrocellulose filter (Sartorius Stedim Biotech, FR) to separate the biomass from the liquid. The biomass on the filters was dried for 24 h at 45 °C before being scraped from the filter and kept in a Falcon tube at room temperature.

Description: Pit lakes in the ‘anthropogenic lake district’ in the Muskau Arch (western Poland; central Europe) are strongly affected by acid mine drainage (AMD). The studied acidic pit lake, ŁK-61 (pH 〈3), is also exposed to floods due to its location in the flood hazard area, which may significantly influence the geochemical behavior of elements. The elemental compositions of water and lake sediment samples were measured with ICP–OES and ICP–MS. The sediment profile was also examined for 137Cs and 210Po activity concentrations using gamma and alpha spectrometry, respectively. Grain size distribution, mineralogical composition, diatoms, and organic matter content in the collected core were also determined. The key factors responsible for the distribution of selected heavy metals (e.g., Cu, Ni, Pb, Zn) and radioisotopes (137Cs and 210Po) in the bottom sediments of Lake ŁK-61 are their coprecipitation/precipitation with Fe and Al secondary minerals and their sorption onto authigenic and allogenic phases. These processes are likely driven by the lake tributary, which is an important source of dissolved elements. The data also showed that the physiochemical parameters of Lake ŁK-61 water changed during an episodic depositional event, i.e., the flood of the Nysa Łużycka River in the summer of 2010. The flood caused an increase in the water pH, as interpreted from the subfossil diatom studies. The down-core profiles of the studied heavy metal and radionuclide (HMRs) contents were probably affected by this depositional event, which prevented a detailed age determination of the collected lake sediments with 137Cs and 210Pb dating methods. Geochemical modeling indicates that the flood-related shift in the physicochemical parameters of the lake water could have caused the scavenging of dissolved elements by the precipitation of fresh secondary minerals. Moreover, particles contaminated with HMRs have also possibly been delivered by the river, along with the nutrients (e.g., phosphorus and nitrogen).