ALBERT

Description: The deflection and the control of the effects of the complex urban seismic wavefield on the built environment is a major challenge in earthquake engineering. The interactions between the soil and the structures and between the structures strongly modify the lateral variability of ground motion seen in connection to earthquake damage. Here we investigate the idea that flexural and compressional resonances of tall turbines in a wind farm strongly influence the propagation of the seismic wavefield. A large-scale geophysical experiment demonstrates that surface waves are strongly damped in several distinct frequency bands when interacting at the resonances of a set of wind turbines. The ground-anchored arrangement of these turbines produces unusual amplitude and phase patterns in the observed seismic wavefield, in the intensity ratio between stations inside and outside the wind farm and in surface wave polarization while there is no metamaterial-like complete extinction of the wavefield. This demonstration is done by setting up a dense grid of 400 geophones and another set of radial broadband stations outside the wind farm to study the properties of the seismic wavefield propagating through the wind farm. Additional geophysical equipment (e.g., an optical fiber, rotational and barometric sensors) was used to provide essential explanatory and complementary measurements. A numerical model of the turbine also confirms the mechanical resonances that are responsible for the strong coupling between the wind turbines and the seismic wavefield observed in certain frequency ranges of engineering interest.

Description: Teleseismic back-projection imaging has emerged as a powerful tool for understanding the rupture propagation of large earthquakes. However, its application often suffers from artifacts related to the receiver array geometry. We developed a teleseismic back-projection technique that can accommodate data from multiple arrays. Combined processing of P and pP waveforms may further improve the resolution. The method is suitable for defining arrays ad-hoc to achieve a good azimuthal distribution for most earthquakes. We present a catalog of short-period rupture histories (0.5–2.0 Hz) for all earthquakes from 2010 to 2022 with MW ≥ 7.5 and depth less than 200 km (56 events). The method provides automatic estimates of rupture length, directivity, speed, and aspect ratio, a proxy for rupture complexity. We obtained short-period rupture length scaling relations that are in good agreement with previously published relations based on estimates of total slip. Rupture speeds were consistently in the sub-Rayleigh regime for thrust and normal earthquakes, whereas a tenth of strike-slip events propagated at supershear speeds. Many rupture histories exhibited complex behaviors, for example, rupture on conjugate faults, bilateral propagation, and dynamic triggering by a P wave. For megathrust earthquakes, ruptures encircling asperities were frequently observed, with downdip, updip, and balanced patterns. Although there is a preference for short-period emissions to emanate from central and downdip parts of the megathrust, emissions updip of the main asperity are more frequent than suggested by earlier results.

Description: The proximity of fast‐slipping crustal faults to urban areas may result in pulse‐like ground motions from rupture directivity, which can contribute to increased levels of damage even for engineered structures. Systematic modeling of directivity within probabilistic seismic hazard analysis (PSHA) remains challenging to implement at the regional scale, despite the availability of directivity models in the literature. In the process of developing the 2022 National Seismic Hazard Model for New Zealand (2022 NSHM), we explored the feasibility and impact of modeling directivity for PSHA at a national scale using the previous generation 2010 NSHM. The results of this analysis allowed us to quantify the impact of directivity on the resulting seismic hazard maps for New Zealand and gain insights into the factors that contribute to the expected increases (and decreases) in ground‐motion level. For the 2022 NSHM, the earthquake rupture forecast (ERF) seismogenic source models introduced enormous challenges for directivity modeling due to the abundance of large multisegment or multifault ruptures with complex geometries. To overcome these challenges, we applied a machine learning‐based strategy to “overfit” an artificial neural network to capture the distributions of directivity amplification and its variability for each unique rupture in the earthquake rupture forecast. This produces a compact representation of the spatial fields of amplification that are computationally efficient to generate within a complete PSHA calculation for the 2022 NSHM. This flexible and reproducible framework facilitates the implementation of directivity in PSHA at a regional scale for complex ERF source models and opens up the possibility of more complex characterization of epistemic uncertainties for near‐source ground motion in practice.

Description: The morphology of coastal sequences provides fundamental observations to unravel past sea level (SL) variations. For that purpose, converting morphometric observations into a SL datum requires understanding their morphogenesis. The long-lasting sequence of coral reef terraces (CRTs) at Cape Laundi (Sumba Island, Indonesia) could serve as a benchmark. Yet, it epitomizes a pitfall that challenges the ultimate goal: the overall chronology of its development remains poorly constrained. The polycyclic nature of the terraces, involving marine erosion and reoccupation of old coral colonies by more recent ones hinders any clear assignment of Marine Isotope Stages (MIS) to specific terraces, in particular the reference datum corresponding to the last Interglacial maximum (i.e., MIS 5e). Thus, to overcome these obstacles, we numerically model the genesis of the sequence, testing a range of eustatic SL (ESL) reconstructions and uplift rates, as well as exploring the parameter space to address reef growth, erosion and sedimentation. A total of 625 model runs allowed us to improve the morpho-chronological constraints of the coastal sequence and, more particularly, to explain the morphogenesis of the several CRTs associated with MIS 5e. Our results suggest that the lowermost main terrace was first constructed during the marine transgression of MIS 5e and was later reshaped during the marine regression of MIS 5e, as well as during the MIS 5c and MIS 5a highstands. Finally, we discuss the general morphology of the sequence and the implications it may have on SL reconstructions. At Cape Laundi, as elsewhere, we emphasize the necessity of addressing the development of CRT sequences with a dynamic approach, that is, considering that a CRT is a landform built continuously throughout the history of SL oscillations, and not simply during a singular SL maximum.

Description: Hydraulic fracturing has been widely used to enhance reservoir permeability during the extraction of shale gas. As one of the external input parameters, injection rate has a significant impact on formation breakdown pressure and the complexity of hydraulic fractures. To gain deeper insights into the effect of injection rate on breakdown pressure and fracture morphology, we conducted five hydraulic fracturing experiments on Changning shale in the laboratory. We used five different injection rates between 3 and 30 mL/min to fracture cylindrical core samples with 50 mm in diameter and 100 mm in length. We monitored acoustic emissions and surface displacements during the tests, and analyzed the fracture pattern post mortem by using a fluorescent tracer. We find a semi-logarithmic relationship between the breakdown pressure and the injection rates. Second, we find that it is the injection rate that dictates sample deformation and crack formation during breakdown rather than the fluid volume injected during the whole process. The analysis of amplitudes and frequency of acoustic signals indicates that hydraulic fracturing of Changning shale is overall dominated by tensile fractures (〉 60%). However, at low injection rates, shear events are facilitated before rock breakdown. On the other hand, high injection rates result in reducing fracture tortuosity and surface roughness due to limited fluid infiltration in the relatively short injection window. We close this study with a conceptual model to explain the difference between fluid infiltration (low injection rates) and the loading rate effect (high injection rate) in low-permeability shale rocks. The findings obtained in this study can help to adjust injection rates in the field to economically and safely produce gas from shale.

Description: Secondary ion mass spectrometry was used to test the d18O and d34S nanogram-scale homogeneity of a suite of candidate sulfate minerals, ultimately selecting three barite, two anhydrite, and two gypsum samples from the Royal Ontario Museum that have repeatabilities for their SIMS measurements of better than 0.39‰ and 0.37‰ (1s) for oxygen and sulfur isotope ratios, respectively. Metrological splits of each of the seven materials were sent to multiple gas source isotope ratio mass spectrometry laboratories in order to establish their absolute 18O/16O and 34S/32S ratios. The inter-laboratory results of GS-IRMS analyses yielded reasonably narrow ranges in d18OVSMOW, whereas larger variations in d34SVCDT values were found between the results from the gas source laboratories. All samples have good reproducibility within laboratories of GS-IRMS 103d18O values of between 0.24‰ and 0.44‰ (1s). The reproducibility within laboratories of GS-IRMS 103d34S values range from 0.07‰ to 0.99‰ (1s). Here we also discuss some of the current analytical limitations affecting these isotope-mineral systems. A total of 256 metrological splits have been prepared from each of these seven materials; these aliquots will be made available to the global geochemical community.

Description: The goals of the Paris Agreement (PA) on collectively managing climate change can only be reached if all parties to the United Nations Framework Convention on Climate Change (UNFCCC) commit to actions supporting their Nationally Determined Contributions (NDCs). Developing-economy nations play a crucial role in reaching the PA targets, particularly in the Agriculture, Forest, and Other Land Uses (AFOLU) sector. However, developing country Parties also face several constraints in tracking and communicating progress towards their climate policy targets and implementation of their NDCs. The operationalization of Biennial Transparency Report (BTR) and Enhanced Transparency Framework (ETF) under the PA will bring stricter reporting timeframes and advanced transparency for all parties. With these requirements rapidly coming into force, addressing reporting gaps is now a pressing priority. The present study analyzes the NDCs, and Biennial Update Reports (BURs) submitted by developing country Parties under the UNFCCC. In an illustrative exercise, our in-depth analysis concentrates on reporting on the AFOLU sector and identifies issues impeding a comprehensive and comparable Global Stock Take (GST): (i) issues of consistency in reporting timeframes (ii) issues in transparency of reporting on mitigation sectors and on relevant progress indicators (iii) incomparability of methodological approaches proposed and used, and (iv) the implications of limited national capacity for transparent reporting. The UNFCCC and developed country Parties now have the opportunity of providing specialized support for developing country Parties. This could include tailored guidance to address gaps in both greenhouse-gas (GHG) emissions accounting, and reporting challenges, to ensure consistent, comprehensive, and transparent reporting to reinforce capacities moving forward following the next GST.

Description: The increasing demand for fertilizers and their rising prices has led to the search for new nutrient sources, especially in rural areas where family farming predominates. In this study, we assessed the potential of reusing sediment deposited in surface reservoirs as a soil conditioner in a semiarid region, focusing on two features: the characterization of sediment physicochemical properties at the regional scale and the effect of the substrate containing sediment on the growth and physiology of maize.

Description: The crustal motions throughout Germany have not yet been fully understood because the research scope of previous studies often focuses only on some active grabens. Thus, we investigate it in detail to identify the neotectonic motion characteristics and specific deformation-ongoing regions. High accuracy for monitoring and data analyses is required because the expected crustal deformation in Germany is small. For this reason, we use high-precision GNSS time series processing techniques and interdisciplinary data to reflect actual motions and determine the causes of deformation. Also, an advanced technique of discontinuity correction is introduced to unify the fragments of the GNSS coordinate time series for better velocity field reliability. Our findings show that the crustal motions in Germany tend to increase at a maximum speed of +1.0 mm/year. Meanwhile, local subsidence of around 0.8 mm/year is concentrated in the river basins (e.g., the Rhine, Ems, Elbe, Northern Oder, and Danube) and extensive mining regions. The Earth’s crust here also behaves with noticeable compressions. The intra-plate motion in Germany is 0.8 mm/year. A special region with an extension rate of +4.3 nstrain/year is observed along the North–South trending Regensburg-Leipzig-Rostock shear zone. Machine Learning clusters the 3D plate velocity field in Germany into three distinct regions with increasing speeds: Northwest, East, and Southwest. Significant surface deformations are detected mainly in the Upper Rhine graben, Eifel volcanic field, and Thuringian-Vogtland slate mountains. The harmonic motions of the Earth’s crust in Germany have an amplitude of 4.7 mm, in which the surface loads contribute half to this type of motion. The findings will contribute to the overall picture of neotectonics here.

Description: To determine the relationships between the functional trait composition of forest communities and environmental gradients across scales and biomes and the role of species relative abundances in these relationships.

Description: Microbial blooms colonize the Greenland Ice Sheet bare ice surface during the ablation season and significantly reduce its albedo. On the ice surface, microbes are exposed to high levels of irradiance, freeze–thaw cycles, and low nutrient concentrations. It is well known that microorganisms secrete metabolites to maintain homeostasis, communicate with other microorganisms, and defend themselves. Yet, the exometabolome of supraglacial microbial blooms, dominated by the pigmented glacier ice algae Ancylonema alaskanum and Ancylonema nordenskiöldii, remains thus far unstudied. Here, we use a high-resolution mass spectrometry-based untargeted metabolomics workflow to identify metabolites in the exometabolome of microbial blooms on the surface of the southern tip of the Greenland Ice Sheet. Samples were collected every 6 h across two diurnal cycles at 5 replicate sampling sites with high similarity in community composition, in terms of orders and phyla present. Time of sampling explained 46% (permutational multivariate analysis of variance [PERMANOVA], pseudo-F = 3.7771, p = 0.001) and 27% (PERMANOVA, pseudo-F = 1.8705, p = 0.001) of variance in the exometabolome across the two diurnal cycles. Annotated metabolites included riboflavin, lumichrome, tryptophan, and azelaic acid, all of which have demonstrated roles in microbe–microbe interactions in other ecosystems and should be tested for potential roles in the development of microbial blooms on bare ice surfaces.

Description: The analysis of Coulomb stress changes has become an important tool for seismic hazard evaluation because such stress changes may trigger or delay subsequent earthquakes. Processes that can cause significant Coulomb stress changes include coseismic slip and transient postseismic processes such as poroelastic effects and viscoelastic relaxation. However, the combined influence of poroelastic effects and viscoelastic relaxation on co- and postseismic Coulomb stress changes has not been systematically studied so far. Here, we use three-dimensional finite-element models with arrays of normal and thrust faults to investigate how pore fluid pressure changes and viscoelastic relaxation overlap during the postseismic phase. In different experiments, we vary the permeability of the upper crust and the viscosity of the lower crust or lithospheric mantle while keeping the other parameters constant. In addition, we perform experiments in which we combine a high (low) permeability of the upper crust with a low (high) viscosity of the lower crust. Our results show that the coseismic (i.e., static) Coulomb stress changes are altered by the signal from poroelastic effects and viscoelastic relaxation during the first month after the earthquake. For sufficiently low viscosities, the Coulomb stress change patterns show a combined signal from poroelastic and viscoelastic effects already during the first postseismic year. For sufficiently low permeabilities, Coulomb stress changes induced by poroelastic effects overlap with the signals from viscoelastic relaxation and interseismic stress accumulation for decades. Our results imply that poroelastic and viscoelastic effects have a strong impact on postseismic Coulomb stress changes and should therefore be considered together when analyzing Coulomb stress transfer between faults.

Description: The impact of faults on the contemporary stress field in the upper crust has been discussed in various studies. Data and models clearly show that there is an effect, but so far, a systematic study quantifying the impact as a function of distance from the fault is lacking. In the absence of data, here we use a series of generic 3-D models to investigate which component of the stress tensor is affected at which distance from the fault. Our study concentrates on the far field, lo- cated hundreds of metres from the fault zone. The models assess various techniques to represent faults, different mate- rial properties, different boundary conditions, variable orien- tation, and the fault’s size. The study findings indicate that most of the factors tested do not have an influence on ei- ther the stress tensor orientation or principal stress magni- tudes in the far field beyond 1000 m from the fault. Only in the case of oblique faults with a low static friction coeffi- cient of μ = 0.1 can noteworthy stress perturbations be seen up to 2000 m from the fault. However, the changes that we detected are generally small and of the order of lateral stress variability due to rock property variability. Furthermore, only in the first hundreds of metres to the fault are variations large enough to be theoretically detected by borehole-based stress data when considering their inherent uncertainties. This find- ing agrees with robust stress magnitude measurements and stress orientation data. Thus, in areas where high-quality and high-resolution data show gradual and continuous stress ten- sor rotations of 〉 20◦ observed over lateral spatial scales of 10 km or more, we infer that these rotations cannot be at- tributed to faults. We hypothesize that most stress orienta- tion changes attributed to faults may originate from different sources such as density and strength contrasts.

Description: The Bakreswar geothermal province represents a medium enthalpy geothermal system with its Bakreswar and Tantloie hot springs. It lies within the Chotanagpur Granite Gneissic Complex in the eastern part of the Indian Peninsula. The province has a high heat flow and a high geothermal gradient of 90°C/km. Magnetotelluric data from 95 sites in a frequency range of 10 kHz–10 Hz were acquired over the Bakreswar geothermal province to obtain an electrical conductivity model and map the geothermal reservoir with its fluid pathways and related geological structures. Subsurface conductivity models obtained from three-dimensional inversions of the Magnetotelluric data exhibit several prominent anomalies, which are supplemented by gravity results. The conductivity model maps three features which act as a conduit (a) a northwest–southeast trending feature, (b) an east–west trending feature to the south of the northwest–southeast trending feature (which lies 1 km north of the Oil and Natural Gas Corporation fault marked by previous studies) and (c) shallow conducting features close to Bakreswar hot spring. The northwest–southeast trending feature coincides with the boundary of the high-density intrusive block. This northwest–southeast trending feature provides the pathway for the meteoric water to reach a maximum depth of 2.7 km, where it gets heated by interacting with deep-seated structures and then it rises towards the surface. The radiogenic process occurring within the granites of Chotanagpur Granite Gneissic Complex provides the heat responsible for heating the meteoric water. The northwest–southeast and east–west trending features are responsible for the transport of meteoric water to deeper depths and then towards the shallow regions of the Earth. The near surface features close to the Bakreswar hot spring are responsible for carrying the water further towards the hot spring. The resistivity of these structures plotted as a function of salinity and temperatures for saline crustal fluids suggests the involvement of meteoric water. Further, applying Archie's law to this resistivity suggests that the conduit path has a porosity greater than 10%. This study successfully maps the anomalous structures which might foster the migration of geothermal fluid in Bakreswar geothermal province.

Description: The Samail Ophiolite in the Oman Mountains formed at a Cretaceous subduction zone that was part of a wider Neo-Tethys plate-boundary system. The original configuration and evolution of this plate-boundary system is hidden in a structurally and metamorphically complex nappe stack below the Samail Ophiolite. Previous work provided evidence for high-temperature metamorphism high in the nappe pile (in the metamorphic sole of the Samail Ophiolite), and high-pressure metamorphism in the deepest part of the nappe pile (Saih Hatat window), possibly reflecting a downward younging, progressive accretion history at the Samail subduction zone. However, there is evidence that the two subduction-related metamorphic events are disparate, but temporally overlapping during the mid-Cretaceous. We present the first geochronologic dataset across the entire high-pressure nappe stack below the Samail Ophiolite, and the shear zones between the high-pressure nappes. Our 22 new Rbsingle bondSr multimineral isochron ages from the Saih Hatat window, along with independent new field mapping and kinematic reconstructions, constrain the timing and geometry of tectonometamorphic events. Our work indicates the existence of a high-pressure metamorphic event in the nappes below the ophiolite that was synchronous with the high-temperature conditions in the metamorphic sole. We argue that the thermal conditions of these synchronous metamorphic events can only be explained through the existence of two Cretaceous subduction zones/segments that underwent distinctly different thermal histories during subduction infancy. We infer that these two subduction zones initially formed at two perpendicular subduction segments at the Arabian margin and subsequently rotated relative to each other and, as a consequence, their records became juxtaposed: (1) The high-temperature metamorphic sole and the Samail Ophiolite both formed above the structurally higher, outboard, ‘hot’ and rotating Samail subduction zone and, (2) the high-pressure nappes developed within the structurally lower, inboard, ‘cold’ Ruwi subduction zone. We conclude that the formation and evolution of both subduction zones were likely controlled by the density structure of the mafic-rock-rich Arabian rifted margin and outermost Arabian Platform, and the subsequent arrival of the buoyant, largely mafic-rock-free, full-thickness Arabian lithosphere, which eventually halted subduction at the southern margin of Neo-Tethys. Previous article in issue

Description: This data set contains the results from a 2023 GFZ Innovative Research Expedition project to explore for natural hydrogen gas (H2) occurrences in the NW Pyrenean foreland, near the town of Biarritz in France. The data represent in-situ measurements of soil and spring water gas, as well as in-situ spring water property measurements, complemented with laboratory analysis results of gas contents and noble gas isotopic compositions of gas and spring water samples collected during the expedition. This GFZ Innovative Research Expedition was inspired by previous exploration efforts in the region by Lefeuvre et al. (2021, 2022). These authors detected elevated concentrations of natural H2 gas in the soil and interpreted this natural H2 to be derived from serpentinizing mantle rocks below the Pyrenees. The main aims of this expedition were the following: (1) in-situ measuring soil gas contents and taking soil gas samples for laboratory analysis at a site near the town of Peyrehorade in the NW of the general study area of Lefeuvre et al. (2021), thus improving the soil gas data coverage along the NW end of the North Pyrenean Frontal Thrust (NPFT); (2) taking gas samples from degassing springs (or water samples from non-degassing springs to be degassed in the lab) in the general Lefeuvre et al. (2021) study area for additional laboratory analysis of gas contents and noble gas isotopic compositions, which may be indicative of (deep) gas origins; and (3) performing a detailed soil gas analysis by means of a portable mass spectrometer at Sauveterre-de-Béarn, a site along the NPFT where Lefeuvre et al. (2022) measured elevated concentrations of natural H2 in the soil. Furthermore, we also measured the properties of the visited springs (temperature, pH, conductivity) while on site, and performed additional in-situ soil gas measurements from manual drillholes. Details on the measurement and sampling methods, on the laboratory analyses, as well as the results of these measurements and analyses are provided in the data description file The expedition involved six field days in July 2023, during which a total of 26 sites were visited. These sites were selected for their vicinity near a major geological contact or fault zone that could have facilitated upward circulation of gas or (thermal) water from the (deep) subsurface (i.e., potentially from the mantle).

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: The gravity anomalies reflect density perturbations at different depths, which control the physical state and dynamics of the lithosphere and sub-lithospheric mantle. However, the gravity effect of the crust masks the mantle signals. In this study, we develop two frameworks (correction with density contrasts and actual densities) to calculate the gravity anomalies generated by the layered crust. We apply the proposed approaches to evaluate the global mantle gravity disturbances based on the new crustal models. Consistent patterns and an increasing linear trend of the mantle gravity disturbances with lithospheric thickness and Vs velocities at 150 km depth are obtained. Our results indicate denser lithospheric roots in most cratons and lighter materials in the oceanic mantle. Furthermore, our gravity map corresponds well to regional geological features, providing new insights into mantle structure and dynamics. Specifically, (1) reduced anomalies associated with the Superior and Rae cratons indicate more depleted roots compared with other cratons of North America. (2) Negative anomalies along the Cordillera (western North America) suggest mass deficits owing to the buoyant hot mantle. (3) Positive anomalies in the Baltic, East European, and Siberian cratons support thick, dense lithosphere with significant density heterogeneities, which could result from thermo-chemical modifications of the cratonic roots. (4) Pronounced positive anomalies correspond to stable blocks, e.g., Arabian Platform, Indian Craton, and Tarim basin, indicating a thick, dense lithosphere. (5) Low anomalies in the active tectonic units and back-arc basins suggest local mantle upwellings. (6) The cold subducting/detached plates may result in the high anomalies observed in the Zagros and Tibet.

Description: The 1400 km long North Anatolian Fault Zone in Türkiye runs through numerous densely populated regions, including the city of Düzce that was recently hit by an Mw 6.1 earthquake on 23 November 2022. This was the first moderate event in the region after the devastating Mw 7.2 earthquake in 1999, which cost the lives of over 700 people. Despite its moderate size, the earthquake caused unexpected severe damage to a significant number of buildings, as reported by local institutions (Disaster and Emergency Management Presidency, AFAD). It is well established that ground motions in the near field can lead to increased damage due to near-field domain effects, such as groundmotion pulses and directivity effects (i.e., when the site is aligned with rupture propagation). We examine potential near-field effects using the strong ground motion database of AFAD-Turkish Accelerometric Database and Analysis Systems. To achieve this, we first analyze the behavior of the ground-motion intensities in terms of their spatial distribution and observe higher peak ground velocity than expected by ground-motion models in spatially constrained azimuthal ranges. Furthermore, we find that the majority of the near-fault recordings contain velocity pulses that are primary concentrated on the fault-parallel component. This outcome questions the widely accepted understanding from the previous studies, which mainly suggested that impulsive ground motions that are associated with directivity effects primarily occur on the fault-normal component of large-magnitude events.

Description: We report uplift and shortening rates from a late Neogene–Pleistocene deformation stage of the frontal fold-thrust belt and adjacent wedge-top in the Principal Cordillera of the southern Central Andes (33-39° SL). A structural model is presented based on integration of surface field data and subsurface 2D seismic sections. Shortening, uplift, and sedimentation rates were calculated from different steps of kinematic modelling. Our structural interpretations and modelling are integrated with new detrital zircon U-Pb geochronology to define a previously overlooked Pleistocene period of orogenic shortening and syntectonic sedimentation in the Malargüe basin. This task was possible due to the dating of three samples yielding between ∼12 and 1 Ma obtained from a 900 m deep well located in the foreland. From stratigraphic correlations, our data records an active Plio-Pleistocene wedge-top depozone coeval with retreat of the volcanism, and the emplacement of retroarc basalts. Structural modelling, together with detrital zircon U-Pb provenance data register shortening producing a foredeep to wedge-top Plio-Pleistocene transition, adjusting and completing the knowledge of the frontal fold-thrust belt and foreland basin in the southern Central Andes.

Description: Arsenic and silica are known inhibitors of the crystallization of iron minerals from poorly ordered precursor phases. However, little is known about the effects of co-existing As and Si on the crystallization and long-term stability of mixed-valence Fe minerals such as green rust (GR). GR usually forms in anoxic, Fe2+-rich, near-neutral pH environments, where they influence the speciation and mobility of trace elements, nutrients and contaminants. In this work, the Fe2+-induced transformation of As- and/or Si-bearing ferrihydrite (FHY) was monitored at pH 8 ([As]initial = 100 μM, Si/As = 10) over 720 h. Our results showed that in the presence of As(III) + Si or As(V) + Si, GR sulfate (GRSO4) formation from FHY was up to four times slower compared to single species system containing only As(III), As(V) or Si. Co-existing As(III) + Si and As(V) + Si also inhibited GRSO4 transformation to magnetite, contrary to systems with only Si or As(V). Overall, our findings demonstrate the synergistic inhibitory effect of co-existing Si on the crystallization and solid-phase stability of As-bearing GRSO4, establishing an inhibitory effect ladder: As(III) + Si 〉 As(V) + Si 〉 As(III) 〉 Si 〉 As(V). This further highlights the importance of GR in potentially controlling the fate and mobility of As in ferruginous, Si-rich groundwater and sediments such as those in South and Southeast Asia.

Description: Lava flows form important fluid reservoirs and have been extensively exploited for water aquifers, geothermal energy, hydrocarbon production and, more recently, for carbon storage. Effusive subaerial mafic to intermediate lava flows account for vast rock volumes globally, and form reservoirs with properties dictated by well-known lava flow facies ranging from pāhoehoe through several transitional forms to ‘a’ā lava. These variations in flow type lead to critical differences in the pore structure, distribution, connectivity, strength and fracturing of individual lava flows, which, alongside lava flow package architectures, determine primary reservoir potential. Lava flow margins with vesicular, fracture and often autobreccia-hosted pore structures can have porosities commonly exceeding 40% and matrix permeabilities over 10−11 m2 (〉10 D) separated by much lower porosity and permeability flow interiors. Secondary post-emplacement physicochemical changes related to fracturing, meteoric, diagenetic and hydrothermal alteration can significantly modify reservoir potential through a complex interplay of mineral transformation, pore-clogging secondary minerals and dissolution, which must be carefully characterized and assessed during exploration and appraisal. Within this contribution, a review of selected global lava flow-hosted reservoir occurrences is presented, followed by a discussion of the factors that influence lava flow reservoir potential.

Description: For an ω2-source model, moment-based estimates of the stress drop are obtained by combining corner frequency and seismic moment source parameters. Therefore, the moment-based estimates of the stress drop are informative about the amount of energy radiated at high frequencies by dynamic rupture processes. This study aims to systematically estimate such stress drop from the harmonized dataset at the European scale and to characterize the distributions of the stress drop for application in future stochastic simulations. We analyze the seismological records associated with shallow crustal seismic events that occurred in Western Europe between January 1990 and May 2020. We processed 220,000 high-quality records and isolated the contributions of the source, site, and path contributions using the Generalized Inversion Technique. The source parameters, including the corner frequency, moment magnitude, and stress drop, of 6135 seismic events are calculated. The events processed are mainly tectonic events (e.g., earthquakes of the central Italy 2009–2016 sequence), although non-tectonic events associated with the Groningen gas field and mining activities in Western Europe are also included in the analysis. The impact of different attenuation models and reference site choices are evaluated. Most of the obtained source spectra follow the standard ω2-model except for a few events where the data sampling considered does not allow an effective spectral decomposition. The resulting stress drop shows a positive correlation with moment magnitude between 3 and 4, and a self-similarity for magnitudes greater than 4 with a mean stress drop of 13.8 MPa.

Description: During flood events near real-time, synthetic aperture radar (SAR) satellite imagery has proven to be an efficient management tool for disaster management authorities. However, one of the challenges is accurate classification and segmentation of flooded water. A common method of SAR-based flood mapping is binary segmentation by thresholding, but this method is limited due to the effects of backscatter, geographical area, and surface characterstics. Recent advancements in deep learning algorithms for image segmentation have demonstrated excellent potential for improving flood detection. In this paper, we present a deep learning approach with a nested UNet architecture based on a backbone of EfficientNet-B7 by leveraging a publicly available Sentinel‑1 dataset provided jointly by NASA and the IEEE GRSS Committee. The performance of the nested UNet model was compared with several other UNet-based convolutional neural network architectures. The models were trained on flood events from Nebraska and North Alabama in the USA, Bangladesh, and Florence, Italy. Finally, the generalization capacity of the trained nested UNet model was compared to the other architectures by testing on Sentinel‑1 data from flood events of varied geographical regions such as Spain, India, and Vietnam. The impact of using different polarization band combinations of input data on the segmentation capabilities of the nested UNet and other models is also evaluated using Shapley scores. The results of these experiments show that the UNet model architectures perform comparably to the UNet++ with EfficientNet-B7 backbone for both the NASA dataset as well as the other test cases. Therefore, it can be inferred that these models can be trained on certain flood events provided in the dataset and used for flood detection in other geographical areas, thus proving the transferability of these models. However, the effect of polarization still varies across different test cases from around the world in terms of performance; the model trained with the combinations of individual bands, VV and VH, and polarization ratios gives the best results.

Description: To trace Critical Zone processes and to quantify Li fluxes from one Critical Zone compartment into another, we investigated the Li concentration and isotopic composition (δ7Li) of time-series water samples (including subsurface flow, groundwater and creek water), vegetation, bedrock (including separated minerals from bedrock), and regolith (including exchangeable fraction and clay-sized fraction of regolith) in a temperate forested headwater catchment in the Black Forest, Conventwald, Germany. Our estimation of the Li budget shows that atmospheric deposition and biological processes played minor roles in the Li cycle relative to chemical weathering. All water samples (δ 7Li value of 6.5 to 20.4 ‰) were enriched in 7Li compared to bedrock (-1.3 ‰) and regolith (∼-1.7 ‰), but δ7Li differed between water pathways: δ7Li variations in subsurface flow, creek water and groundwater were controlled by conservative mixing, exchangeable pool buffering and Li incorporation/adsorption, respectively. Fractionated heavy Li isotopes in water samples resulted from the formation of secondary solids which preferentially incorporated 6Li, with the separated clay-sized fraction of the regolith exhibiting more negative δ7Li values (-5.4 to −3.5 ‰) than the bulk regolith (∼-1.7 ‰). However, Li in secondary solids only accounted for 8 ± 6 % of the total Li hosted in bulk regolith, and consequently δ7Li in soil did not differ significantly from δ 7Li in bedrock. This is unexpected considering water is continuously removing 7Li in preference over 6Li from regolith. Mass balance calculations applied at the catchment scale point to an irreconcilable imbalance with our data. On one hand, the regolith’s δ7Li values are not negative enough to balance the 7Li export by river water, and on the other hand Li in the riverine dissolved load only accounts for ∼ 30 % of the Li solubilized from regolith. Therefore, we suggest that there might be a “hidden export pathway” for Li at our site, possibly subsurface removal of fine particles enriched in 6Li. In light of increasingly frequent observations of such isotopic imbalances in the Critical Zone this phenomenon deserves increased attention.

Description: Estuarine processes are key in modulating the riverine input of particle-reactive trace elements to the ocean. An important, but still under-utilized member of these elements is beryllium-9 (9Be) that together with cosmogenic 10Be has been suggested to serve as a quantitative tracer of present and past continental weathering flux. This study investigates different pathways of terrigenous 9Be through coastal areas into the ocean, based on dissolved 9Be concentrations in surface and bottom waters together with corresponding particulate 9Be concentrations along the salinity gradient in the Changjiang Estuary. Dissolved 9Be in the Changjiang Estuary shows a non-conservative behavior: At low to mid-salinity where water is well-mixed, 9Be is removed from both surface and bottom waters at low salinity and then released back into the water column at mid-salinity. At high salinity where water is stratified, dissolved 9Be is removed from surface waters, but is released back into bottom waters. In combination with hydrochemical (e.g., dissolved oxygen) and particulate 9Be data obtained from different extracted phases, we attribute the removal of dissolved 9Be at low salinity to salt-induced colloidal flocculation, whereas in surface waters at high salinity, we ascribe the removal to biological scavenging facilitated by phytoplankton blooms. The release of 9Be into mid- and high-salinity bottom waters is likely dominated by benthic processes, including porewater diffusion and/or submarine groundwater discharge. The contribution from desorption of 9Be from suspended particulate matter is negligible throughout the entire estuary. We propose that the release of 9Be through benthic processes potentially presents the most important contributor to the marine 9Be budget, where this benthic flux of 9Be is likely enhanced by hypoxic conditions in coastal bottom waters.

Description: The recent rapid improvement of machine learning techniques had a large impact on the way seismological data can be processed. During the last years several machine learning algorithms determining seismic onset times have been published facilitating the automatic picking of large data sets. Here we apply the deep neural network PhaseNet to a network of over 900 permanent and temporal broad band stations that were deployed as part of the AlpArray research initiative in the Greater Alpine Region (GAR) during 2016-2020. We selected 384 well distributed earthquakes with ML ≥ 2.5 for our study and developed a purely data-driven pre-inversion pick selection method to consistently remove outliers from the automatic pick catalog. This allows us to include observations throughout the crustal triplication zone resulting in 39,599 P and 13,188 S observations. Using the established VELEST and the recently developed McMC codes we invert for the 1D P- and S-wave velocity structure including station correction terms while simultaneously relocating the events. As a result we present two separate models differing in the maximum included observation distance and therefore their suggested usage. The model AlpsLocPS is based on arrivals from ≤ 130 km and therefore should be used to consistently (re)-locate seismicity based on P & S observations. The model GAR1D_PS includes the entire observable distance range of up to 1000 km and for the first time provides consistent P- & S-phase synthetic travel times for the entire Alpine orogen. Comparing our relocated seismicity with hypocentral parameters from other studies in the area we quantify the absolute horizontal and vertical accuracy of event locations as ≈ 2.0 km and ≈ 6.0 km, respectively.

Description: Significance Continental glaciers and ice sheets are excellent indicators of ongoing and past climate changes. The Patagonian ice sheet (PIS) was the largest extrapolar ice sheet in the Southern Hemisphere. Many studies have investigated the advances of the PIS on its eastern side, but there are only a few PIS records on the Pacific side. We show that three active intervals occurred during the last ~140 ka, with an extended PIS that contributed to the release of large amounts of freshwater and sediment into the Pacific. Active intervals during the last glacial period occurred from ~70 to 60 ka and from ~40 to 18 ka, with four and five phases of increased ice discharge, respectively, most likely driven by precipitation changes. Abstract Terrestrial glacial records from the Patagonian Andes and New Zealand Alps document quasi-synchronous Southern Hemisphere–wide glacier advances during the late Quaternary. However, these records are inherently incomplete. Here, we provide a continuous marine record of western–central Patagonian ice sheet (PIS) extent over a complete glacial–interglacial cycle back into the penultimate glacial (~140 ka). Sediment core MR16-09 PC03, located at 46°S and ~150 km offshore Chile, received high terrestrial sediment and meltwater input when the central PIS extended westward. We use biomarkers, foraminiferal oxygen isotopes, and major elemental data to reconstruct terrestrial sediment and freshwater input related to PIS variations. Our sediment record documents three intervals of general PIS marginal fluctuations, during Marine Isotope Stage (MIS) 6 (140 to 135 ka), MIS 4 (~70 to 60 ka), and late MIS 3 to MIS 2 (~40 to 18 ka). These higher terrigenous input intervals occurred during sea-level low stands, when the western PIS covered most of the Chilean fjords, which today retain glaciofluvial sediments. During these intervals, high-amplitude phases of enhanced sediment supply occur at millennial timescales, reflecting increased ice discharge most likely due to a growing PIS. We assign the late MIS 3 to MIS 2 phases and, by inference, older advances to Antarctic cold stages. We conclude that the increased sediment/meltwater release during Southern Hemisphere millennial-scale cold phases was likely related to higher precipitation caused by enhanced westerly winds at the northwestern margin of the PIS. Our records complement terrestrial archives and provide evidence for PIS climate sensitivity.

Description: Postglacial flooding of the Persian Gulf (PG) was important in shaping human history and driving landscape changes in the region. However, there is a paucity of data regarding the postglacial transgression. The position of the PG at the edge of major synoptic systems of the Indian Ocean Summer Monsoon (IOSM) and Mid-latitude Westerlies (MLW) makes the environment particularly sensitive to Holocene climate shifts. To investigate the timing of the flooding and to detect the impacts of significant climate shifts on the regional environment during the Holocene, a multiproxy study was conducted on three short sediment cores from two deep sites in the PG. Sedimentological, palynological and geochemical analyses were performed on the cores. The results show that inundation of the western part of the PG that started from ca. 11.5 ka bp continued with successive prominent phases of transgression centered on 10.4 and 9.2 ka cal bp, and definitive marine conditions were established around 8.8 ka cal bp. The IOSM was the dominant system in the region until about 9 to ~6.3 ka cal bp. After that time, the intensity of the IOSM declined, as MLW dominated the region after ~6.3 ka cal bp. These climatic shifts induced significant changes in regional vegetation and hydrology, and possibly triggered socio-cultural transformations.

Description: The Caucasus region is key for understanding early human dispersal and evolution in Eurasia, and characterizing the environmental contrast between Last Glacial Maximum and Holocene is crucial for investigating human adaptation strategies to large climatic shifts. However, a paucity of high-resolution paleoclimate records leave this context largely unknown for early human populations in the Caucasus region. Based on our model-proxy comparison of high- and low-resolution records of 24 stalagmites from three caves, we find spatially distinct changes in vegetation and seasonality of precipitation, especially under glacial conditions. Supported by modern oxygen-isotope data and climate modeling, we identify a supraregional cold-season temperature control for oxygen isotopes in Black Sea speleothems, which previously had been interpreted as a local moisture-source signal. Carbon-isotope and trace-element data further suggest disproportionate changes in vegetation cover and soil dynamics at high altitudes, which would have resulted in a reduction but not a disappearance of human refugia during the Last Glacial Maximum, relative to the current interglacial. Our findings imply that abrupt climatic pressures from harsh conditions were overcome by adaptive strategies in the past.

Description: Continuous monitoring of volcanic gas emissions is crucial for understanding volcanic activity and potential eruptions. However, emissions of volcanic gases underwater are infrequently studied or quantified. This study explores the potential of Distributed Acoustic Sensing (DAS) technology to monitor underwater volcanic degassing. DAS converts fiber-optic cables into high-resolution vibration recording arrays, providing measurements at unprecedented spatio-temporal resolution. We conducted an experiment at Laacher See volcano in Germany, immersing a fiber-optic cable in the lake and interrogating it with a DAS system. We detected and analyzed numerous acoustic signals that we associated with bubble emissions in different lake areas. Three types of text-book bubbles exhibiting characteristic waveforms are all found from our detections, indicating different nucleation processes and bubble sizes. Using clustering algorithms, we classified bubble events into four distinct clusters based on their temporal and spectral characteristics. The temporal distribution of the events provided insights into the evolution of gas seepage patterns. This technology has the potential to revolutionize underwater degassing monitoring and provide valuable information for studying volcanic processes and estimating gas emissions. Furthermore, DAS can be applied to other applications, such as monitoring underwater carbon capture and storage operations or methane leaks associated with climate change.

Description: Understanding the temporal and spatial environmental response to past climate change during the Last Glacial-Interglacial Transition (LGIT, 16-8 ka) across Europe relies on precise chronologies for palaeoenvironmental records. Tephra layers (volcanic ash) are a powerful chronological tool to synchronise disparate records across the continent. Yet, some regions remain overlooked in terms of cryptotephra investigations. Building on earlier work at the same sites, we present the first complete LGIT high-resolution cryptotephra investigation of two lake records in the Carpathian Mountains in Romania, Lake Brazi and Lake Lia. Numerous volcanic glass shards have been recognised as originating from various volcanic regions, including: Iceland (Katla, Askja, and Torfajokull), Italy (Campi Flegrei, Ischia, Lipari, and Pantelleria), and central Anatolia (Acigol and Ericyes). In total, four distinct tephra horizons have now been identified in these records: 1) an LGIT Lipari tephra (11,515–12,885 cal BP, 95.4% range); 2) Askja-S (11,070–10,720 cal BP, 95.4% range); 3) an Early Holocene Lipari tephra,(12,590–10,845 cal BP, 95.4% range) and; 4) an Early Holocene Ischia tephra (11,120–10,740 cal BP, 95.4% range). The use of trace element analysis on selected cryptotephra layers provided additional important information in identifying volcanic source and facilitating correlations. These tephra layers, along with numerous other discrete cryptotephra layers, offer promise as significant future isochrons for comprehending the spatial and temporal fluctuations in past climate change throughout Europe and the Mediterranean area. This research has emphasized the significance of the Carpathian region in expanding the European and Mediterranean tephra lattice and establishing it as a keystone area within the framework.

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: The eastern Alpine crust has been shaped by the continental collision of the European and Adriatic plates beginning at 35 Ma and was affected by a major reorganization after 20 Ma. To better understand how the eastern Alpine surface structures link with deep seated processes, we analyze the depth-dependent seismic anisotropy based on Rayleigh wave propagation. Ambient noise recordings are evaluated to extract Rayleigh wave phase dispersion measurements. These are inverted in a two step approach for the azimuthally anisotropic shear velocity structure. Both steps are performed with a reversible jump Markov chain Monte Carlo (rj-McMC) approach that estimates data errors and propagates the modeled uncertainties from the phase velocity maps into the depth inversion. A two layer structure of azimuthal anisotropy is imaged in the Alpine crust, with an orogen-parallel upper crust and approximately orogen-perpendicular layer in the lower crust and the uppermost mantle. In the upper layer, the anisotropy tends to follow major fault lines and may thus be an apparent, structurally driven anisotropy. The main foliation and fold axis orientations might contribute to the anisotropy. In the lower crust, the N-S orientation of the fast axis is mostly confined to regions north of the Periadriatic Fault and may be related to European subduction. Outside the orogen, no clearly layered structure is identified. The anisotropy pattern in the northern Alpine foreland is found to be similar compared to SKS studies which is an indication of very homogeneous fast axis directions throughout the crust and the upper mantle.

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: Although many collisional orogens form after subduction of oceanic lithosphere between two continents, some orogens result from strain localization within a continent via inversion of structures inherited from continental rifting. Intracontinental rift-inversion orogens exhibit a range of structural styles, but the underlying causes of such variability have not been extensively explored. We use numerical models of intracontinental rift inversion to investigate the impact of parameters including rift structure, rift duration, post-rift cooling, and convergence velocity on orogen structure. Our models reproduce the natural variability of rift-inversion orogens and can be categorized using three endmember styles: asymmetric underthrusting (AU), distributed thickening (DT), and localized polarity flip (PF). Inversion of narrow rifts tends to produce orogens with more localized deformation (styles AU and PF) than those resulting from wide rifts. However, multiple combinations of the parameters we investigated can produce the same structural style. Thus, our models indicate no unique relationship between orogenic structure and the conditions prior to and during inversion. Because the style of rift-inversion orogenesis is highly contingent upon the rift history prior to inversion, knowing the geologic history that preceded rift inversion is essential for translating orogenic structure into the processes that produced that structure.

Description: The distinction between the shear behavior of tensile- and shear-induced fractures is critical to understanding the deformation and failure of geologic discontinuities at different scales. To investigate these differences, a series of direct shear tests were performed on sandstone specimens with a continuous fracture created by either splitting or shearing. The acoustic emission (AE) technique was used to examine variations in grain-size cracking behavior between specimens with tensile- and shear-induced fractures. An increase in normal stress for both fracture types correlates with increased microcrack density and energy release. However, there are notable differences: during the shear process, tensile-induced fractures produce AE sequences similar to the seismic patterns observed along natural tectonic faults, with foreshocks, mainshocks, and aftershocks. In contrast, the AE sequence for shear-induced fractures during the shear process lacks prominent mainshocks and deviates progressively from the power-law function with time as normal stress increases. In addition, the AE b-value for tension-induced fractures initially shows a gradual decrease as the mainshock approaches and then slowly increases during the aftershock period. In contrast, the b-value remains nearly constant for shear-induced fractures due to the low roughness and heterogeneity of the fracture surface. These differences highlight the strong correlation between AE responses and fault heterogeneity, paving the way for fault characterization and risk assessment in subsurface energy extraction.

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: Although the general impacts of zooplankton grazing on phytoplankton communities are clear, we know comparatively less about how specific grazing strategies interact with environmental conditions to shape the size structure of phytoplankton communities. Here, we present a new data-driven, size-based model that describes changes in the size composition of lake phytoplankton under various environmental constraints. The model includes an ecological trade-off emerging from observed allometric relationships between (1) phytoplankton cell size and phytoplankton growth and (2) phytoplankton cell size and zooplankton grazing. In our model, phytoplankton growth is nutrient-dependent and zooplankton grazing varies according to specific grazing strategies, namely, specialists (targeting a narrow range of the size-feeding spectrum) vs. generalists (targeting a wide range of the size-feeding spectrum). Our results indicate that grazing strategies shape the size composition of the phytoplankton community in nutrient-rich conditions, whereas inorganic nutrient concentrations govern phytoplankton biomass. Under oligotrophic regimes, the phytoplankton community is dominated by small cell sizes and the grazers have little to no impact. Under eutrophic regimes, dominating specialist grazers push phytoplankton towards small cells, whereas dominating generalist grazers push phytoplankton towards large cells. Our work highlights that trait-based modeling, based on realistic eco-physiological trade-offs, represents a valuable tool for disentangling the interactive roles played by nutrient regimes and grazing strategies in determining the size compositions of lake phytoplankton. Ultimately, our study offers a quantitative basis for understanding how communities of lake phytoplankton may reorganize in the future in response to changes in nutrient levels and zooplankton grazing strategies.

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: Plankton community modeling is a critical tool for understanding the processes that shape marine ecosystems and their impacts on global biogeochemical cycles. These models can be of variable ecological, physiological, and physical complexity. Many published models are either not publicly available or implemented in static and inflexible code, thus hampering adoption, collaboration, and reproducibility of results. Here we present Phydra, an open-source library for plankton community modeling, and Xarray-simlab-ODE (XSO), a modular framework for efficient, flexible, and reproducible model development based on ordinary differential equations. Both tools are written in Python. Phydra provides pre-built models and model components that can be modified and assembled to develop plankton community models of various levels of ecological complexity. The components can be created, adapted, and modified using standard variable types provided by the XSO framework. XSO is embedded in the Python scientific ecosystem and is integrated with tools for data analysis and visualization. To demonstrate the range of applicability and how Phydra and XSO can be used to develop and execute models, we present three applications: (1) a highly simplified nutrient–phytoplankton (NP) model in a chemostat setting, (2) a nutrient–phytoplankton–zooplankton–detritus (NPZD) model in a zero-dimensional pelagic ocean setting, and (3) a size-structured plankton community model that resolves 50 phytoplankton and 50 zooplankton size classes with functional traits determined by allometric relationships. The applications presented here are available as interactive Jupyter notebooks and can be used by the scientific community to build, modify, and run plankton community models based on differential equations for a diverse range of scientific pursuits.

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: The data set was collected to identify hydrological processes and their evolution over it time. It consists of several individual files in tabstop delimeted text format. The data set contains the data obtained from deuterium and brilliant blue tracer experiments at two chronosequence studies in the glacier forefield of the Stone Glacier and the Griessfirn in the central Alps, Switzerland. Each chronosequence consisted of four moraines of different ages (from 30 to 13500 years). At each forefield sprinkling experiments with deuterium and dye tracer experiments with blue dye (Brilliant Blue) were conducted on three plots per moraine. The moraines at the forefield of the Stone Glacier developed from siliceous parent material and at the forefield of the Griessfirn from calcareous parent material. Data from the siliceous forefield are marked with (S) and data from the calcareous forefield are marked with (C). The data set consist of soil moisture time series and soil water isotope profiles of the sprinkling experiments with deuterium, as well as trinary images of stained vertical subsurface flow paths from the dye tracer experiment. The individual plots per moraine are distinguished via their position relative to one another on the moraine (left, middle, and right, looking upslope). The plots used for the sprinkling experiments were located in close vicinity to the plots used for the dye tracer experiments. For the sprinkling experiments with deuterium each plot (4m x 6m) per age class was equipped with 6 soil moisture sensors. Three of these sensors were installed as a sensor profile at one side of the plot about one meter downslope from the upper plot boundary. The sensors were installed at 10, 30, and 50 cm soil depth. On the other side of the plot, two sensors were placed in 10 cm depth, one opposite to the sensor profile and the second sensor one meter upslope from the lower plot boundary. The sixth sensor was placed at 10 cm depth in the center of the plot. The plots were irrigated on three consecutive days with three different irrigation intensities and deuterium concentrations. Per forefield, the soil moisture data are listed in one file per age class. The file contains for each plot, the time stamp and the soil moisture values of the 6 sensors.

Description: Silicate weathering and organic carbon (OC) burial in soil regulate atmospheric CO2, but their influence on each other remains unclear. Generally, OC oxidation can generate acids that drive silicate weathering, yet clay minerals that form during weathering can protect OC and limit oxidation. This poses a conundrum where clay formation and OC preservation either compete or cooperate. Debate remains about their relative contributions because quantitative tools to simultaneously probe these processes are lacking while those that exist are often not measured in concert. Here we demonstrate that Li isotope ratios of sediment, commonly used to trace clay formation, can help constrain OC cycling. Measurements of river suspended sediment from two watersheds of varying physiography and analysis of published data from Hawaii soil profiles show negative correlations between solid-phase values and OC content, indicating the association of clay mineral formation with OC accumulation. Yet, the localities differ in their ranges of values and OC contents, which we interpret with a model of soil formation. We find that temporal trends of Li isotopes and OC are most sensitive to mineral dissolution/clay formation rates, where higher rates yield greater OC stocks and lower values. Whereas OC-enhanced dissolution primarily dictates turnover times of OC and silicate minerals, clay protection distinctly modifies soil formation pathways and is likely required to explain the range of observations. These findings underscore clay mineral formation, driven primarily by bedrock chemistry and secondarily by climate, as a principal modulator of weathering fluxes and OC accumulation in soil.

Description: This data repository for the Southern Caribbean and NW South America contains a 3D thermal model computed down to 75 km depth, the modelled hypocentral temperatures and geothermal gradients at the locations of crustal earthquakes, and the crustal seismogenic depths calculated from earthquake statistics, as well as the associated modelled temperatures. We used the uppermost 75 km of the gravity-constrained structural and density model of Gómez-García et al. (2020, 2021) to derive the 3D thermal configuration of the study area (5°-15° N, 63°-82° W). A steady-state approach was followed, in which upper and lower boundary conditions were set to run the thermal calculations using the software GOLEM (Cacace & Jacquey, 2017; Jacquey & Cacace, 2017). A catalogue of earthquakes occurred within the study area and surroundings was compiled from public sources. In the database archived here, we provide data of the best located crustal earthquakes within the boundaries of this area, from January 1980 to June 2021. Earthquakes below the magnitude of completeness, or with poorly determined depths, were disregarded. Earthquakes were deemed crustal if their hypocentres were located between the topo-bathymetry from the GEBCO relief (Weatherall et al., 2015) and the Moho depth from the GEMMA model (Reguzzoni & Sampietro, 2015). We computed the crustal seismogenic depth as the 90th and 95th percentiles (D90 and D95), respectively, of the crustal hypocentral depths. These percentiles were mapped on a latitude-longitude grid, using for each grid node at least the 20 closest earthquakes as sample. The hypocentral temperatures, the geothermal gradient at the earthquake locations, and the temperatures at the D90 and D95 surfaces were calculated from the lithospheric-scale thermal model. For more details about the modelling approach and interpretation of the results, we kindly ask the reader to refer to the main publication: Gomez-Garcia et al. (2024).

Description: Introduction: Urban pluvial flooding is a growing concern worldwide as consequence of rising urban population and climate change induced increases in heavy rainfall. Easy-to-implement and fast simulation tools are needed to cope with this challenge. Methods: This study describes the development of the parsimonious, GPU-accelerated hydraulic model RIM2D for urban pluvial flood simulations. This is achieved by considering the built-up urban area as flow obstacles, and by introducing capacity-based approaches to consider urban drainage by infiltration on pervious surfaces and sewer drainage from roofs and sealed surfaces. The model performance was analyzed by simulating 8 heavy rainfall events in a test area in the city of Dresden, Germany. For these events detailed discharge measurements of sewer discharge are available, providing a unique dataset for evaluating the sewer drainage simulation, which is of high importance for realistic pluvial inundation simulations in urban areas. Results and discussion: We show that the model simulates the temporal dynamics of the sewer discharge and the sewer volume within acceptable ranges. Moreover, the erratic variation of the simulated to measured sewer discharge suggests that the deviations from the measurements are caused by the precipitation input rather than the model simplifications. We conclude that RIM2D is a valid tool for urban inundation simulation. Its short simulation runtimes allow probabilistic flood risk assessments and operational flood forecasts.

Description: Repeating earthquakes, or repeaters, affecting overlapping rupture patches with a similar focal mechanism, have important implications to track fault slip rates, aseismic deformation, slow earthquakes and earthquake nucleation processes. They are often detected based on highly similar waveforms. Here, we discuss earthquakes with highly anti-correlated waveforms, denoting a reversed seismogenic process at the same or a neighbouring location, which we refer to as true and quasi anti-repeaters. We first report a range such observations in different environments, including volcano seismicity, intermediate depth seismicity and injection-induced microseismicity. Then, we review conceptual models proposed to explain them. True and quasi anti-repeaters can be robustly identified via a three-component single station or distributed network data. They are key indicators for stress perturbation transients or local stress heterogeneities. Since most of these observations were explained as the response to fluid migration processes, they may help to identify and track fluid movements in the subsurface.

Description: The Hipercorig Hallstatt History (H3) project aims to unravel the entire Late Glacial to Holocene sedimentary succession recording past climate, environment, natural hazard impacts, human–environment interactions, and prehistoric mining history. We successfully cored 51 m of the sedimentary succession of Lake Hallstatt, revealing a high-resolution Late Pleistocene to Holocene sediment record, overcoming the previous coring limit of 15.63 m (dated to ∼ 2.3 ka cal BP). The novel drilling platform Hipercorig allows the recovery of undisturbed long cores and the acquisition of borehole logging data of deep lakes. The sedimentary record is spliced to a composite core profile, and for the first time borehole logging data are linked to a core–log seismic correlation of an intra-mountainous lake of the Eastern Alps. The recovered sequence consists of two major lithostratigraphic units: (i) Unit 1 (Holocene, 0–41.7 m below lake floor) with 10 (up to 5.1 m thick) instantaneous deposits and (ii) Unit 2 (Late Pleistocene, 〉 41.7 m below lake floor). The Late Pleistocene sediments comprise the Younger Dryas and the deepest recovered sediments likely date back to the Allerød interstadial. Within the Holocene, six different periods are observed in the core and borehole logging data, showing distinct physical property fluctuations and an overall increase in sedimentation rate upcore. Lake Hallstatt provides a unique prehistoric archive, being located within the UNESCO World Heritage area Hallstatt–Dachstein/Salzkammergut, a region with a rich history of human salt mining dating back to 3400 cal BP (Middle/Late Bronze Age) and one of the oldest documented cultural landscapes worldwide.

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: We discuss data of three laboratory stick-slip experiments on Westerly Granite samples performed at elevated confining pressure and constant displacement rate on rough fracture surfaces. The experiments produced complex slip patterns including fast and slow ruptures with large and small fault slips, as well as failure events on the fault surface producing acoustic emission bursts without externally-detectable stress drop. Preparatory processes leading to large slips were tracked with an ensemble of ten seismo-mechanical and statistical parameters characterizing local and global damage and stress evolution, localization and clustering processes, as well as event interactions. We decompose complex spatio-temporal trends in the lab-quake characteristics and identify persistent effects of evolving fault roughness and damage at different length scales, and local stress evolution approaching large events. The observed trends highlight labquake localization processes on different spatial and temporal scales. The preparatory process of large slip events includes smaller events marked by confined bursts of acoustic emission activity that collectively prepare the fault surface for a system-wide failure by conditioning the large-scale stress field. Our results are consistent overall with an evolving process of intermittent criticality leading to large failure events, and may contribute to improved forecasting of large natural earthquakes.

Description: The empirical Båth’s law states that the average magnitude difference (⁠⁠) between a mainshock and its strongest aftershock is ∼1.2, independent of the size of the mainshock. Although this observation can generally be explained by a scaling of aftershock productivity with mainshock magnitude in combination with a Gutenberg–Richter frequency–magnitude distribution, estimates of may be preferable because they are directly related to the most interesting information, namely the magnitudes of the main events, without relying on assumptions. However, a major challenge in calculating this value is the bias introduced by missing data points when the strongest aftershock is below the observed cut‐off magnitude. Ignoring missing values leads to a systematic error because the data points removed are those with particularly large magnitude differences ⁠. The error can be minimized by restricting the statistics to mainshocks that are at least 2 magnitude units above the cutoff, but then the sample size is strongly reduced. This work provides an innovative approach for modeling by adapting methods for time‐to‐event data, which often suffer from incomplete observations (censoring). In doing so, we adequately account for unobserved values and estimate a fully parametric distribution of the magnitude differences for mainshocks in a global earthquake catalog. Our results suggest that magnitude differences are best modeled by the Gompertz distribution and that larger are expected at increasing depths and higher heat flows.

Description: Using Global Navigation Satellite System-Reflectometry (GNSS-R) for soil moisture (SM) retrieval has recently gained importance due to its high temporal-spatial resolution. However, the current methods, i.e., constructing a single machine learning (ML)-based model, have large model uncertainty resulting from ML networks and input schemes. Moreover, traditional Normalized Difference Vegetation Index (NDVI) cannot capture the rapid vegetation changes well. In this paper, a new SM retrieval method of constructing a hybrid model based on Bayesian model averaging (BMA) is employed to reduce the model uncertainty. Meanwhile, novel Sun-induced fluorescence (SIF) data is used as ancillary data to represent the rapid change of vegetation. We validate the proposed method at point and regional scales using in-situ data and the Global Land Data Assimilation System (GLDAS) product. The results demonstrate that our method has high accuracy and low uncertainty in SM retrieval. At the point scale, as accuracy indices, the average R () of BMA increases from 0.90 to 0.93 and the average root-mean-square-error () decreases from 0.034 to 0.029 ; as indices of uncertainty, the standard deviations of R and RMSE ( and ) decrease by 32 % and 9 % compared to the single ML-based model. For the regional scale, the increases from 0.79 to 0.81, the decreases from 0.024 to 0.023 , and the decreases by 19 %. Moreover, we take the point-scale experiment as an example for comparison to compare the performance of SIF with that of NDVI. The of BMA trained by SIF is 0.03 higher than that trained by NDVI and the decreases by 0.002 ; and decrease by 25 % and 6 %. Based on these results, the proposed method can reduce the uncertainty and the advantage of SIF has potential for improving the SM retrieval.

Description: The Lacq area in southwest France has been associated with continuous moderate induced seismic activity since 1969. However, the mechanisms driving this induced seismicity are not fully understood: reservoir depletion has been proposed as the main factor, and more recently wastewater injection has been suggested to play a more important role (Grasso et al., 2021). The interpretation of these mechanisms relies heavily on the quality of earthquake locations, which we prove to be weak due to a lack of local instrumentation for several years. In order to provide the most complete and reliable induced event catalog for the studies of the Lacq induced seismicity mechanisms & seismic hazard, we made an exhaustive compilation, analysis and improvement of all available catalogs. We also provided new earthquake detections & relocations in a 3D velocity model from past and present temporary deployments never used for studying the Lacq area. Important remaining location uncertainties lead us to also carefully sort the events according to their location confidence, defining 3 classes of events (unconstrained location, location constrained within 2-3 km and 1-2 km respectively). This new harmonized catalog and the identification of well-constrained events, covering 50 years of induced seismicity, allow us to propose that wastewater injection is almost certainly the main mechanism driving the seismicity, with (i) most of the constrained events located within the reservoir boundaries and (ii) the released seismic energy variations following variations in injection operations at different scales. In particular, we have also highlighted a change in the injection-seismicity relationship around 2010–2013. From 2013, despite lower injection volumes, seismicity remained persistent and some clusters of earthquakes were detected predominantly in spring, summer, and early autumn, except in winter periods. From 2016, we observed a strong temporal relationship between days with higher rate/volume injections (approximately above 400m3/day) and both clustered events and higher magnitude earthquakes (greater than 2.4).

Description: This dataset reports measurements from a laboratory incubation of soils sourced from a boreal peatland and surrounding habitats (Siikaneva Bog, Finland). In August 2021, soil cores were collected from three habitat zones: a well-drained upland forest, an intermediate margin ecotone, and a Sphagnum moss bog. The cores from each habitat were taken from surface to approximately 50cm below surface using an Eijelkamp peat corer and subdivided by soil horizon. The samples were then incubated anaerobically for 140 days in three temperature treatment groups (0, 4, 20°C). Subsamples of the incubations headspace (250 µL) were measured on a gas chromatograph (7890A, Agilent Technologies, USA) with flame ionization detection (FID) for CO2 and CH4 concentrations. The rate of respiration from the samples were calculated per gram carbon and per gram soil as described in the method of Robertson., et al. (1999) and reported here, along with other relevant parameters.

Description: The 2022 revision of Aotearoa New Zealand National Seismic Hazard Model (NZ NSHM2022) has involved significant revision of all datasets and model components. In this article, we present a subset of many results from the model as well as an overview of the governance, scientific, and review processes followed by the NZ NSHM team. The calculated hazard from the NZ NSHM 2022 has increased for most of New Zealand when compared with the previous models. The NZ NSHM 2022 models and results are available online.

Description: Underwater Distributed Acoustic Sensing (DAS) utilizes optical fiber as a continuous sensor array. It enables high-resolution data collection over long distances and holds promise to enhance tsunami early warning capabilities. This research focuses on detecting infragravity and tsunami waves associated with earthquakes and understanding their origin and dispersion characteristics through frequency-wavenumber domain transformations and beamforming techniques. We propose a velocity correction method based on adjusting the apparent channel spacing according to water depth to overcome the challenge of detecting long-wavelength and long-period tsunami signals. Experimental results demonstrate the successful retrieval of infragravity and tsunami waves using a subsea optical fiber in offshore Oregon. These findings underscore the potential of DAS technology to complement existing infragravity waves detection systems, enhance preparedness, and improve response efforts in coastal communities. Further research and development in this field are crucial to fully utilize the capabilities of DAS for enhanced tsunami monitoring and warning systems.

Description: A seismicity rate model (SRM) has been developed as part of the 2022 Aotearoa New Zealand National Seismic Hazard Model revision. The SRM consists of many component models, each of which falls into one of two classes: (1) inversion fault model (IFM); or (2) distributed seismicity model (DSM). Here we provide an overview of the SRM and a brief description of each of the component models. The upper plate IFM forecasts the occurrence rate for hundreds of thousands of potential ruptures derived from the New Zealand Community Fault Model version 1.0 and utilizing either geologic- or geodeticbased fault-slip rates. These ruptures are typically less than a couple of hundred kilometers long, but can exceed 1500 km and extend along most of the length of the country (albeit with very low probabilities of exceedance [PoE]). We have also applied the IFM method to the two subduction zones of New Zealand and forecast earthquake magnitudes of up to ∼Mw 9.4, again with very low PoE. The DSM combines a hybrid model developed using multiple datasets with a non-Poisson uniform rate zone model for lower seismicity regions of New Zealand. Forecasts for 100 yr are derived that account for overdispersion of the rate variability when compared with Poisson. Finally, the epistemic uncertainty has been modeled via the range of models and parameters implemented in an SRM logic tree. Results are presented, which indicate the sensitivity of hazard results to the logic tree branches and that were used to reduce the overall complexity of the logic tree.

Description: We present a kinematic model developed from geodetic observations, topography analysis and analogue tectonic modelling results, which reveals a striking similarity between the rotational tectonic settings of the Gakkel Ridge-Chersky Range system in the Arctic, and the Central Indian Tectonic Zone within the Indian subcontinent. A crucial aspect of large-scale extensional rift systems is the gradual variation of extension along the rift axis, due to plate rotation about a Euler pole, which may lead to contraction on the opposite side of the Euler pole to form a rotational tectonic system. Our geodetic and topographic analysis, combined with the reanalysis of analogue tectonic modelling results demonstrates such rotational tectonic plate motion in both the Arctic and Indian case. However, the plate boundary between the North American and Eurasian Plates as represented by the Arctic Gakkel Ridge-Chersky Range system is strongly localized, whereas the Central Indian Tectonic Zone that separates the North and South India Plates involves diffuse deformation instead. Furthermore, in both the Arctic and Central Indian we find that the relative Euler rotation pole is located near an indenter-like feature, which possibly controls the present-day rotational tectonics and contrasting topography on opposite sides of the Euler pole.

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: To test whether a globally inferred sediment thickness value from geomorphological studies can be used as a proxy to predict earthquake site amplification, we derive site-amplification models from the relation between empirical amplification for sites in Europe and Türkiye and the geomorphological sediment thickness. The new site-amplification predictions are then compared to predictions from site-amplification models derived using the traditional site proxies, VS30 inferred from slope, slope itself, and geological era and slope combined. The ability of each proxy to capture the site amplification is evaluated based on the reduction in site-to-site variability caused by each proxy. The results show that the highest reduction is caused by geological era and slope combined, while the geomorphological sediment thickness shows a slightly larger or equal reduction in site-to-site variability as inferred VS30 and slope. We therefore argue that including geology and geomorphology in site-amplification modelling on regional scale can give an important added value and that globally or regionally inferred models for soil and sediment thickness from fields beyond engineering seismology can have a great potential in regional seismic hazard and risk assessments. Furthermore, the differences between the site-amplification maps derived from different proxies capture the epistemic uncertainty of site-amplification modelling. While the different proxies predict similar features on a large scale, local differences can be large. This shows that using only one proxy when predicting site amplification does not capture the full epistemic uncertainty, which is demonstrated by looking into detail on the site-amplification maps predicted for eastern Türkiye and Syria, where the devastating Kahramanmaraş earthquake sequence occurred in February 2023.

Description: Assessing alternative agricultural water management strategies requires long-term field trials or vast data collection for model calibration and simulation. This work aims to assess whether an uncalibrated agro-hydrological model using global input datasets for climate, soil and crop information can serve as a decision support tool for crop water management under data scarcity. This study employs the Cool Farm Tool Water (CFTW) at eight eddy covariance sites of the FLUXNET2015 dataset. CFTW is tested using global (CFTWglobal) and local (CFTWlocal) input datasets under current and alternative management scenarios. Results show that the use of global datasets for estimating daily evapotranspiration had little effect on the median Root Mean Square Error ( ) (CFTWglobal: 1.70 mm, CFTWlocal: 1.79 mm), while, however, the median model is much greater (CFTWglobal: −18.6%, CFTWlocal: −4.3%). Furthermore, the periods of water stress were little affected by the use of local or global data (median accuracy: 0.84), whereas the use of global data inputs led to a significant overestimation of irrigation water requirements (median difference: 110 mm). The model performance improves predominantly through the use of more representative local precipitation data, followed by local reference evapotranspiration and soil for some European growing seasons. We identify model outputs that can support decision-making when relying on global data, such as periods of water stress and the daily dynamics of water use. However, our findings also emphasize the difficulty of overcoming data scarcity in decision-making in agricultural water management. Furthermore, we provide recommendations for enhancing model performance and thus may increase the accessibility of reliable decision support tools in the future.

Description: The death toll andmonetary damages from landslides continue to rise despite advancements in predictive modeling. These models’ performances are limited as landslide databases used in developing them often miss crucial information, e.g., underlying movement types. This study introduces a method of discerning landslide movements, such as slides, flows, and falls, by analyzing landslides’ 3D shapes. By examining landslide topological properties, we discover distinct patterns in their morphology, indicating different movements including complex ones with multiple coupled movements. We achieve 80- 94% accuracy by applying topological properties in identifying landslide movements across diverse geographical and climatic regions, including Italy, the US Pacific Northwest, Denmark, Turkey, and Wenchuan in China. Furthermore, we demonstrate a real-world application on undocumented datasets from Wenchuan. Our work introduces a paradigm for studying landslide shapes to understand their underlying movements through the lens of landslide topology, which could aid landslide predictive models and risk evaluations.

Description: We present transient simulations of the last glacial inception using the Earth system model CLIMBER-X with dynamic vegetation, interactive ice sheets, and visco-elastic solid Earth responses. The simulations are initialized at the middle of the Eemian interglacial (125 kiloyears before present, ka) and run until 100 ka, driven by prescribed changes in Earth's orbital parameters and greenhouse gas concentrations from ice core data. CLIMBER-X simulates a rapid increase in Northern Hemisphere ice sheet area through MIS5d, with ice sheets expanding over northern North America and Scandinavia, in broad agreement with proxy reconstructions. While most of the increase in ice sheet area occurs over a relatively short period between 119 and 117 ka, the larger part of the increase in ice volume occurs afterwards with an almost constant ice sheet extent. We show that the vegetation feedback plays a fundamental role in controlling the ice sheet expansion during the last glacial inception. In particular, with prescribed present-day vegetation the model simulates a global sea level drop of only ∼ 20 m, compared with the ∼ 35 m decrease in sea level with dynamic vegetation response. The ice sheet and carbon cycle feedbacks play only a minor role during the ice sheet expansion phase prior to ∼ 115 ka but are important in limiting the deglaciation during the following phase characterized by increasing summer insolation. The model results are sensitive to climate model biases and to the parameterization of snow albedo, while they show only a weak dependence on changes in the ice sheet model resolution and the acceleration factor used to speed up the climate component. Overall, our simulations confirm and refine previous results showing that climate–vegetation–cryosphere feedbacks play a fundamental role in the transition from interglacial to glacial states characterizing Quaternary glacial cycles.

Description: The Atmosphere and Ocean non-tidal De-aliasing Level-1B (AOD1B) product is widely used in precise orbit determination and satellite gravimetry to correct for transient effects of atmosphere–ocean mass variability that would otherwise alias into monthly mean global gravity fields. The most recent release is based on the global ERA5 reanalysis and ECMWF operational data together with simulations from the general ocean circulation model MPIOM consistently forced with fields from the corresponding atmospheric dataset. As background models are inevitably imperfect, residual errors will consequently propagate into the resulting geodetic products. Accounting for uncertainties of the background model data in a statistical sense, however, has been shown before to be a useful approach to mitigate the impact of residual errors leading to temporal aliasing artefacts. In light of the changes made in the new release RL07 of AOD1B, previous uncertainty assessments are deemed too pessimistic and thus need to be revisited. We here present an analysis of the residual errors in AOD1B RL07 based on ensemble statistics derived from different atmospheric reanalyses, including ERA5, MERRA2 and JRA55. For the oceans, we investigate the impact of both the forced and intrinsic variability through differences in MPIOM simulation experiments. The atmospheric and oceanic information is then combined to produce a new time-series of true errors, called AOe07, which is applicable in combination with AOD1B RL07. AOe07 is further complemented by a new spatial error variance–covariance matrix. Results from gravity field recovery simulation experiments for the planned Mass-Change and Geosciences International Constellation (MAGIC) based on GFZ’s EPOS software demonstrate improvements that can be expected from rigorously implementing the newly available stochastic information from AOD1B RL07 into the gravity field estimation process.

Description: Bones and teeth are often the only fossil remains of vertebrates that are preserved over geological time in sedimentary rocks. They render valuable archives for geochemical proxies which are commonly used for paleo-reconstructions. However, the fossilization mechanisms of bone are not yet well understood. Crucial processes are the transformation of bioapatite and the replacement of collagen by thermodynamically more stable apatite phases, such as fluorapatite. In the present study, aqueous alteration experiments on cortical bone samples were performed under simulated early diagenetic conditions in order to investigate whether and how an external PO4 source affects post mortem bone mineralization. Additionally, abiotic oxygen isotope exchange mechanisms between bioapatite and aqueous solutions were assessed by using either 18O-enriched water or phosphate as a tracer. The presence of an external sedimentary 18O-labeled PO4 source led to a rapid formation of new fluorapatite crystallites at the sample’s margin that was highly enriched in 18O. Meanwhile, in the interior of the samples carbonate-poor HAp formed through a dissolution-precipitation process without incorporating any of the 18O tracer. These two processes appear to act independently from each other. In samples exposed to 18O-labeled aqueous solutions lacking a PO4 source, no newly grown apatite crystallites were found, however in the interior of these samples, nano-crystalline carbonate-poor or -free hydroxylapatite precipitated. A comparatively low but uniform 18O-enrichment was measured from the sample’s margin towards its interior, which is assumed to have resulted from the adsorption of H218O onto crystallite surface sites and collagen. Overall, our results suggest that a fast incorporation of 18O-doped PO4 from the sediment source accelerated bone mineralization and consequently changed the oxygen isotope composition of the PO4 group in the bone mineral phase more rapidly and to a greater extent than in a diagenetic setting lacking additional external PO4. Dissolved phosphate from the taphonomic setting thus seems to be an important factor fostering bone fossilization and preservation as well as oxygen isotope alteration.

Description: Phosphatized fish fossils occur in various locations worldwide. Although these fossils have been intensively studied over the past decades they remain a matter of ongoing research. The mechanism of the permineralization reaction itself remains still debated in the community. The mineralization in apatite of a whole fish requires a substantial amount of phosphate which is scarce in seawater, so the origin of the excess is unknown. Previous research has shown that alkaline phosphatase, a ubiquitous enzyme, can increase the phosphate content in vitro in a medium to the degree of saturation concerning apatite. We applied this principle to an experimental setup where fish scales were exposed to commercial bovine alkaline phosphatase. We analyzed the samples with SEM and TEM and found that apatite crystals had formed on the remaining soft tissue. A comparison of these newly formed apatite crystals with fish fossils from the Solnhofen and Santana fossil deposits showed striking similarities. Both are made up of almost identically sized and shaped nano-apatites. This suggests a common formation process: the spontaneous precipitation from an oversaturated solution. The excess activity of alkaline phosphatase could explain that effect. Therefore, our findings could provide insight into the formation of well-preserved fossils.

Description: Strain engineering beyond substrate limitation of colossal magnetoresistant thin (La0.6Pr0.4)0.7Ca0.3MnO3 (LPCMO) films on LaAlO3-buffered SrTiO3 (LAO/STO) substrates has been demonstrated using metalorganic aerosol deposition technique. By growing partially relaxed 7–27 nm thick heteroepitaxial LAO buffer layers on STO a perfect lattice matching to the LPCMO has been achieved. As a result, strain-free heteroepitaxial 10–20 nm thick LPCMO/LAO/STO films with bulk-like ferromagnetic metallic ground state were obtained. Without buffer the coherently strained thin LPCMO/STO and LPCMO/LAO films were insulating and weakly magnetic. The reason for the optimized magnetotransport in strain-free LPCMO films was found to be a large octahedral Mn–O–Mn bond angle φOOR ~ 166–168° as compared to the significantly smaller one of φOOR ~ 152–156° determined for the tensile (LPCMO/STO) and compressively (LPCMO/LAO) strained films.

Description: The movement of the Earth's surface mass, including the atmosphere and oceans, as well as hydrology and glacier melting, causes the redistribution of surface loads, deformation of the solid Earth, and fluctuations in the gravity field. Global Navigation Satellite Systems (GNSS) provide useful information about the movement of the Earth's surface mass. The impact of surface loading deformation over 145 GNSS sites in Africa was investigated using vertical height time series analysis. The study investigates and quantifies the impact of surface loading on the GNSS coordinates utilizing GNSS Precise Point Positioning (PPP) approach. The German Research Center for Geosciences (GFZ) EPOS.P8 software was used to process and analyze eleven years of GPS data from all the stations, as well as dedicated hydrological and atmospheric loading correction models given by the Earth System Modeling group at Deutsches GeoForschungsZentrum (ESMGFZ). The results of the hydrological loading corrections arising from the surface-deformation were analysed to determine the extent of station improvements. The results revealed about 40% of the stations showed improvement with an average Root Mean Square Error (RMSE) residual of 7.3 mm before the application of the hydrological loading corrections and 7.1 mm Root Mean Square Error (RMSE) after the application of the hydrological loading corrections. Similarly, the atmospheric loading corrections gave an improvement of about 57%. Furthermore, the amplitude values decreased from 4.1–8.1 mm to 3.5–6.2 mm after atmospheric loading corrections. This finding presupposes that applying loading corrections to the derived time series reduces amplitude in some African regions.

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: Determining the timing of landslide occurrence is crucial for establishing an accurate, comprehensive and systematic landslide inventory while assessing the potential for reducing landslide risk. Unfortunately, many existing landslide inventories lack temporal information such as the precise time of landslide events. Optical and Synthetic Aperture Radar (SAR) sensors are the most commonly used remote sensing technologies for landslide detection. Unlike optical sensors, SAR sensors are not affected by cloudy conditions and provide valuable imagery regardless of sunlight availability. Therefore, SAR-derived parameters, i.e., SAR amplitude, interferometric coherence, and polarimetric features (alpha and entropy), offer a higher temporal resolution for detecting landslide occurrence times compared to optical data. Despite the advantages, there is currently no universally accepted automatic method for determining the time of landslide events using SAR data. This is due to the lack of anomaly labels and the high time-series volatility in detecting landslide occurrence times. Despite advances in deep-learning methods for anomaly detection in time-series, only a few of them can address these challenges in our case. In this paper, we propose an unsupervised multivariate transformed-based deep-learning model to automatically and efficiently estimate landslide occurrence times using multivariate SAR-derived parameters time-series analysis. The designed gated relative position can increase robustness and temporal context information, by learning global temporal trends in the time-series. Subsequently, the time-series of the anomaly score derived from the proposed Transformer model is analyzed using an adaptive thresholding strategy to dynamically and automatically mark anomalies related to the landslide occurrence. Our research focuses on collapsed landslides characterized by dramatic changes in ground surface topography, with a particular attention for the need of a prior knowledge about landslide boundaries. We assess the performance of the proposed methodology for several collapsed landslides including the July 21, 2020 Shaziba and 23 July, 2019 Shuicheng landslides in China, March 19, 2019 Takht landslide in Iran, June 15, 2018 Jalgyz-Jangak and May 25, 2018 Kugart landslides in Kyrgyzstan, July 7, 2018 Hitardalur landslide in Iceland, and January 25, 2019 Brumadinho landslide in Brazil. In comparison to commonly used neural networks like the LSTM algorithm, our proposed framework leads to a more accurate estimate for the time of landslide failure using time-series of SAR-derived parameters. Furthermore, our results suggest the great potential of SAR data to narrow the time period detected from optical data when used in conjunction with them.

Description: Drainage divides separate Earth’s surface into individual river basins. Divide migration impacts the evolution of landforms, regional climate, ecosystems and biodiversity. In this Review, we assess the processes and dynamics of divide migration and offer insights into the impact on climate and biodiversity. Drainage divides are not static: they can move through the processes of gradual migration that is continuous in unsteady landscapes, or sudden through infrequent river capture events. Divides tend to move in the direction of slower erosion, faster uplift or with horizontal tectonic advection, with rates typically ranging between 0.001 and 10 mm year−1, and a global average of 0.6 mm year−1. Evidence of river capture, such as a sharp change in flow direction with an upstream waterfall, can constrain divide migration history. Topographic metrics, such as cross-divide steepness, can predict the migration of drainage divides towards directions with a lower topographic steepness. Divide migration influences the spatial distribution of regional precipitation, temperature and topographic connectivity between species, thereby affecting biodiversity. For example, freshwater fish can migrate into a new drainage basin through river capture, potentially increasing the species richness. Future research should couple advanced landscape evolution models and observations from field and remote sensing to better investigate divide migration dynamics.

Description: Assessing the impact of global change on hydrological systems is one of the greatest hydrological challenges of our time. Changes in land cover, land use, and climate have an impact on water quantity, quality, and temporal availability. There is a widespread consensus that, given the far-reaching effects of global change, hydrological systems can no longer be viewed as static in their structure; instead, they must be regarded as entire ecosystems, wherein hydrological processes interact and coevolve with biological, geomorphological, and pedological processes. To accurately predict the hydrological response under the impact of global change, it is essential to understand this complex coevolution. The knowledge of how hydrological processes, in particular the formation of subsurface (preferential) flow paths, evolve within this coevolution and how they feed back to the other processes is still very limited due to a lack of observational data. At the hillslope scale, this intertwined system of interactions is known as the hillslope… Die Einschätzung der Auswirkungen des globalen Wandels auf die lokale Hydrologie stellt zweifellos eine der bedeutendsten hydrologischen Herausforderungen unserer Zeit dar. Die zuverlässige Vorhersage der zukünftigen Verfügbarkeit, Menge und Qualität des Wassers in Landschaften gewinnt dabei an höchster Bedeutung. Es herrscht weitgehender Konsens darüber, dass hydrologische Systeme aufgrund des globalen Wandels nicht mehr als statische Gebilde betrachtet werden können. Vielmehr sind sie als ganzheitliche Ökosysteme zu verstehen, in denen hydrologische Prozesse mit biologischen, geomorphologischen und pedologischen Faktoren interagieren und sich gemeinsam entwickeln. Allerdings ist das Wissen über die Entwicklung und die Rückkopplung hydrologischer Prozesse, insbesondere im Hinblick auf die Entstehung unterirdischer (präferenzieller) Fließwege, aufgrund mangelnder Beobachtungsdaten noch stark begrenzt. Das Hauptziel dieser Studie liegt daher darin, anhand von Untersuchungen auf der Hangskala die Wechselwirkung zwischen Struktur und…

Description: Die Einschätzung der Auswirkungen des globalen Wandels auf die lokale Hydrologie stellt zweifellos eine der bedeutendsten hydrologischen Herausforderungen unserer Zeit dar. Die zuverlässige Vorhersage der zukünftigen Verfügbarkeit, Menge und Qualität des Wassers in Landschaften gewinnt dabei an höchster Bedeutung. Es herrscht weitgehender Konsens darüber, dass hydrologische Systeme aufgrund des globalen Wandels nicht mehr als statische Gebilde betrachtet werden können. Vielmehr sind sie als ganzheitliche Ökosysteme zu verstehen, in denen hydrologische Prozesse mit biologischen, geomorphologischen und pedologischen Faktoren interagieren und sich gemeinsam entwickeln. Allerdings ist das Wissen über die Entwicklung und die Rückkopplung hydrologischer Prozesse, insbesondere im Hinblick auf die Entstehung unterirdischer (präferenzieller) Fließwege, aufgrund mangelnder Beobachtungsdaten noch stark begrenzt. Das Hauptziel dieser Studie liegt daher darin, anhand von Untersuchungen auf der Hangskala die Wechselwirkung zwischen Struktur und hydrologischem Verhalten zu erforschen. Dadurch soll ein besseres Verständnis für den Rückkopplungszyklus gewonnen werden. Die vier in dieser Arbeit präsentierten Studien befassen sich mit der Entwicklung hydrologischer, biologischer und physikalischer Eigenschaften von Hängen sowie der Evolution hydrologischer Fließwege über einen Zeitraum von zehntausend Jahren in zwei unterschiedlichen geologischen Geologien. Durch die Analyse von Moränen-Chronosequenzen in zwei Gletschervorfeldern, die jeweils aus silikatreichem bzw. kalkreichem Geschiebemergel entstanden sind, wird die komplexe Koevolution in Abhängigkeit von der geologischen Beschaffenheit des Ausgangsmaterials beleuchtet. Die Ergebnisse zeigen, dass die verschiedenen Eigenschaften von silikatreichem und kalkhaltigem Ausgangsmaterial unter den gegebenen Klimabedingungen zu Unterschieden in Bodeneigenschaften, Durchlässigkeit und Wasserspeicherung führen. Dies wiederum begünstigt unterschiedliche Vegetationstypen und resultiert in vielfältigen Ökosystemen mit variierenden hydrologischen Dynamiken. Interessanterweise treibt das silikatreiche Ausgangsmaterial die Koevolution stärker an. Hierbei spielt der pH-Wert des Bodens eine Schlüsselrolle, da er Auswirkungen auf Vegetation, Bodenbildung und folglich auch auf die Hydrologie hat. Die saure Verwitterung des silikatreichen Materials begünstigt die Anreicherung organischer Substanz, was zu einer Erhöhung der Wasserspeicherkapazität des Bodens und zur Podsolisierung führt. Innerhalb weniger Jahrhunderte entwickelt sich das unterirdische Wassertransportsystem von einer vertikalen, heterogenen Matrixströmung zu fingerartigen Fließwegen und nach mehreren Jahrtausenden zu einem System aus Makroporen, wobei die Wasserspeicherung und die unterirdische, laterale Strömung überwiegen. Im Gegensatz dazu weist der kalkhaltige Standort eine hohe pH-Pufferkapazität auf. In dem basischen bis neutralen Milieu sammelt sich vergleichsweise wenig abgestorbene, organische Substanz an, was zu einer geringeren Wasserspeicherkapazität führt und die Ansiedlung von Grasvegetation begünstigt. Die Fließwege ändern lediglich ihre Form, nicht jedoch ihre Richtung. Es wurde lediglich ein Übergang von vertikaler, heterogener Matrixströmung zu vertikalen, fingerartigen Fließwegen beobachtet. Der schnelle vertikale Wassertransport mit geringer Speicherung begünstigt die Erhaltung der Grasvegetation. Die durchgeführten Studien verdeutlichen, dass die Entwicklung der Fließwege hauptsächlich von den Eigenschaften des Ausgangsgesteins und seinen Verwitterungsprodukten sowie von deren komplexen Interaktionen mit den ursprünglichen Fließwegen und der Entwicklung der Vegetation geprägt ist. Im Gegensatz dazu spielt die Zeit keine entscheidende Rolle bei der Beschreibung der Evolution der hydrologischen Reaktion. Diese Studie trägt wesentlich dazu bei, die bestehende Lücke in Beobachtungen für die Erforschung der Koevolution von hydrologischen, biologischen, geomorphologischen und pedologischen Prozessen zu schließen. Dies ist von großer Bedeutung, um Vorhersagen hydrologischer Prozesse in sich wandelnden Landschaften zu verbessern. Sie zeigt außerdem die Relevanz interdisziplinärer Studien auf, um den zukünftigen Herausforderungen in der Hydrologie im Zuge des globalen Wandels erfolgreich zu begegnen.

Description: Reusing deep wells whether from oil- and gas or geothermal energy production, after years of inactivity can be economically beneficial for the geothermal industry. This study examines the processes that occur in an abandoned well and their impact on reservoir performance. The Groß Sch¨onebeck geothermal research well Gt GrSk 4/05 was reopened after a seven-year production stop to characterize its hydraulic and hydrochemical properties. It targets a reservoir of Rotliegend sandstone and underlying volcanic rocks at a depth of approximately 4000 m. Prior to the hydraulic test, gas samples were collected at the wellhead outlet. These samples showed a significant increase in methane content from approximately 10–15 vol.-% as measured between 2009 and 2013 to 65 vol.-% as measured currently). The objective of this study was to identify the processes responsible for the increase in methane and to assess their impact on well productivity. Well productivity was evaluated through a slug-withdrawal test, followed by the collection of two fluid samples at depths of 1500 m and 4000 m. The chemical composition of the gas and liquid samples was analyzed along with the microbial community. Isotopes of sulfur (34S) and oxygen (18O) were measured in dissolved sulfate and 13C isotopes were measured in the gas phase of hydrocarbons. The particles separated from the suspension of the collected fluid samples were analyzed by scanning electron microscopy. Compared to previous analyses of fluid samples from this well, the methane content has significantly increased and hydrogen sulfide has been detected in the gas phase. Additionally, sulfides of copper and iron were found in the black-colored suspension from the deep-fluid sample taken at 4000 m. The two samples also differ significantly in their organic carbon content, with the sample from 1500 m containing a high amount (135 mg C/L) and at the sample from 4000 m containing a relatively low amount (5 mg C/L). Remarkably, despite the extreme conditions (150 ◦C; 265 g/L TDS at 4000 m), both samples contained several bacterial and fungal genera). Isotope analysis of hydrocarbons in the gas phase indicated mixed methane production from both, thermogenic and biogenic processes. The strong negative 34S of sulfate in the liquid samples rather suggested a thermogenic process of sulfate reduction. It is assumed that several processes occur along the wellbore, including microbial steel corrosion of the well casing (which provides hydrogen, methane, organic carbon, and iron), thermal sulfate reduction, and inflow of formation gases. However, despite the observation of complex geochemical and microbial reactions in the well, the results from the slug withdrawal test indicated no significant change in well productivity (0.6 m3/h*MPa) since the last test in 2014.

Description: Flood losses have steadily increased in the past and are expected to grow even further owing to climate and socioeconomic change. The reduction of flood vulnerability, for example, through adaptation, plays a key role in the mitigation of future flood risk. However, lacking knowledge about vulnerability dynamics, which arise from the interaction between floods and the ensuing response by society, limits the scope of current risk projections. We present a socio-hydrological method for flood risk assessment that simulates the interaction between society and flooding continuously, including changes in vulnerability through collective (structural) and private (non structural) measures. Our probabilistic approach quantifies uncertainties and exploits empirical data to chart risk dynamics including how society copes with flooding. In a case study for the commercial sector in Dresden, Germany, we show that increased adaptation is necessary to counteract the expected four-fold growth in flood risk due to transient hydroclimatic and socioeconomic boundary conditions. We further use our holistic approach to identify solutions for effective long-term adaptation, demonstrating that integrated adaptation strategies (i.e., combined structural and non structural measures) can reduce the average risk by up to 60% at the study site. Ultimately, our case study highlights the benefit of the model for robust flood risk assessment as it can capture unintended, adverse feedbacks of adaptation measures such as the levee effect. Consequently, our socio-hydrological method contributes to a more systemic and reliable flood risk assessment that can inform adaptation planning by exploring the possible system evolutions comprehensively including unlikely futures.

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: U–Pb geochronology of shocked monazite can be used to date hypervelocity impact events. Impact-induced recrystallisation and formation of mechanical twins in monazite have been shown to result in radiogenic Pb loss and thus constrain impact ages. However, little is known about the effect of porosity on the U–Pb system in shocked monazite. Here we investigate monazite in two impact melt rocks from the Hiawatha impact structure, Greenland by means of nano- and micrometre-scale techniques. Microstructural characterisation by scanning electron and transmission electron microscopy imaging and electron backscatter diffraction reveals shock recrystallisation, microtwins and the development of widespread micrometre- to nanometre-scale porosity. For the first time in shocked monazite, nanophases identified as cubic Pb, Pb3O4, and cerussite (PbCO3) were observed. We also find evidence for interaction with impact melt and fluids, with the formation of micrometre-scale melt-bearing channels, and the precipitation of the Pb-rich nanophases by dissolution–precipitation reactions involving pre-existing Pb-rich high-density clusters. To shed light on the response of monazite to shock metamorphism, high-spatial-resolution U–Pb dating by secondary ion mass spectrometry was completed. Recrystallised grains show the most advanced Pb loss, and together with porous grains yield concordia intercept ages within uncertainty of the previously established zircon U–Pb impact age attributed to the Hiawatha impact structure. Although porous grains alone yielded a less precise age, they are demonstrably useful in constraining impact ages. Observed relatively old apparent ages can be explained by significant retention of radiogenic lead in the form of widespread Pb nanophases. Lastly, we demonstrate that porous monazite is a valuable microtexture to search for when attempting to date poorly constrained impact structures, especially when shocked zircon or recrystallised monazite grains are not present.

Description: The response of evapotranspiration to anthropogenic warming is of critical importance for the water and carbon cycle. Contradictory conclusions about evapotranspiration changes are caused primarily by their brevity in time and sparsity in space, as well as the strong influence of internal variability. Here, we present the first gridded reconstruction of the summer (June, July, and August) vapor pressure deficit (VPD) for the past 4 centuries at the European level. This gridded reconstruction is based on 26 European tree ring oxygen isotope records and is obtained using a random forest approach. According to validation scores obtained with the Nash–Sutcliffe model efficiency, our reconstruction is robust over large parts of Europe since 1600, in particular for the westernmost and northernmost regions, where most tree ring records are located. Based on our reconstruction, we show that from the mid-1700s a trend towards higher summer VPD occurred in central Europe and the Mediterranean region that is related to a simultaneous increase in temperature and decrease in precipitation. This increasing summer VPD trend continues throughout the observational period and in recent times. Moreover, our summer VPD reconstruction helps to visualize the local and regional impacts of the current climate change, as well as to minimize statistical uncertainties of historical VPD variability. This paper provides also new insights into the relationship between summer VPD and large-scale atmospheric circulation, and we show that summer VPD has two preferred modes of variability, namely a NW–SE dipole-like mode and a N–S dipole-like mode. Furthermore, the interdisciplinary use of the data should be emphasized, as summer VPD is a crucial parameter for many climatological feedback processes in the Earth's surface system. The reconstructed summer VPD gridded data over the last 400 years are available at the following link: https://doi.org/10.5281/zenodo.5958836 (Balting et al., 2022).

Description: This study presents the mathematical framework and experiments to achieve an orbit and clock combination using Multi-GNSS products. Over the past years, the International GNSS Service (IGS) has been putting efforts into extending its service by setting up and running the Multi-GNSS Experiment and Pilot Project (MGEX). Several MGEX Analysis Centers (ACs) contribute by providing solutions containing not only GPS and GLONASS but also Galileo, BeiDou, and QZSS. As the current IGS combination software can only handle one constellation at a time, obtaining a consistent MGEX orbit and clock product requires substantial modifications. Therefore is presented a least-squares framework for a Multi-GNSS orbit and clock combination, where the weights used to combine the ACs’ products are determined by least-squares variance component estimation. Two different weighting strategies are proposed and discussed: AC specific weights or AC and constellation specific weights. An automated Z-score test is implemented, yielding a common set of core satellites that are used to determine the weights. Both strategies are tested using MGEX orbit and clock solutions for a period of two and a half years, where the two data processing yielded similar results. The agreement with the ACs’ orbits is around one centimeter level for GPS and up to a few centimeters for the other constellations. The 3D-RMS is generally slightly better with the AC and constellation weighting. As a validation of the strategies, a comparison of our combination approach with the official IGS combination using three years of GPS and GLONASS orbits shows an agreement of better than 5 mm and 12 mm for GPS and GLONASS, respectively. An external validation using Satellite Laser Ranging (SLR) shows that the mean residuals of our combined products are around −3 mm for Galileo, 6 mm for GLONASS, −8 mm for BeiDou, and −31 mm for QZSS. For the clock products, the agreement with the IGS final product is around 32 ps and the Precise Point Positioning (PPP) analysis showed that the combined product agrees with the residuals of the ACs with values of less than 10 mm. The MGEX clock combination showed that most of the ACs agreement is below 80 ps for GPS and Galileo, and for the other constellations, the values are more sparse. In addition, it is noted that SHA and JAX presented big discrepancies compared with the other ACs. The validation is performed again in comparison IGS official combination where the agreement with the presented solutions is around 32 ps. A PPP processing is also performed using the combined solutions, showing the suitability of the products. The ACs products from the IGS third reprocessing campaign are also used as input in the combination and discussed. Finally, an extension of the orbit combination is done for the combination of Low Earth Orbit (LEO) satellites. The results show that the proposed combination can deliver reliable products following the standards of the IGS.

Description: Diese Studie stellt den mathematischen Rahmen und die Experimente vor, um eine Kombination von Umlaufbahn und Uhr mit Hilfe von Multi-GNSS-Produkten zu erreichen. In den letzten Jahren hat der International GNSS Service (IGS) Anstrengungen unternommen, um seinen Dienst durch die Einrichtung und Durchführung des Multi-GNSS-Experiment- and Pilot Project (MGEX) zu erweitern. Mehrere MGEX Analysis Centers (ACs) tragen dazu bei, indem sie Lösungen anbieten, die nicht nur GPS und GLONASS, sondern auch Galileo, BeiDou und QZSS enthalten. Da die derzeitige IGS-Kombinationssoftware jeweils nur eine Konstellation verarbeiten kann, sind für den Erhalt eines konsistenten MGEX-Bahn und Uhrenprodukts erhebliche Änderungen erforderlich. Daher wird ein Least-Squares-Rahmen für eine Multi-GNSS-Bahn- und Uhrenkombination vorgestellt, bei dem die zur Kombination der AC-Produkte verwendeten Gewichte durch eine Least-Squares-Varianzkomponentenschätzung bestimmt werden. Es werden zwei verschiedene Gewichtungsstrategien vorgeschlagen und diskutiert: AC-spezifische Gewichte oder AC- und konstellationsspezifische Gewichte. Es wird ein automatischer Z-Score-Test durchgeführt, der einen gemeinsamen Satz von Kernsatelliten ergibt, die zur Bestimmung der Gewichte verwendet werden. Beide Strategien werden anhand von MGEX-Bahn- und Uhrenlösungen über einen Zeitraum von zweieinhalb Jahren getestet, wobei die beiden Datenverarbeitungen ähnliche Ergebnisse erbrachten. Die Übereinstimmung mit den Bahnen der ACs liegt bei GPS bei etwa einem Zentimeter und bei den anderen Konstellationen bei bis zu einigen Zentimetern. Der 3D-RMS ist im Allgemeinen etwas besser mit der AC- und Konstellationsgewichtung. Zur Validierung der Strategien zeigt ein Vergleich unseres Kombinationsansatzes mit der offiziellen IGS-Kombination unter Verwendung von GPS- und GLONASS-Umlaufbahnen aus drei Jahren eine Übereinstimmung von besser als 5 mm und 12 mm für GPS bzw. GLONASS. Eine externe Validierung mittels Satellite Laser Ranging (SLR) zeigt, dass die mittleren Residuen unserer kombinierten Produkte bei etwa −3 mm für Galileo, 6 mm für GLONASS, −8 mm für BeiDou und −31 mm für QZSS liegen. Bei den Uhrenprodukten liegt die Übereinstimmung mit dem IGS-Endprodukt bei etwa 32 ps, und die Analyse der präzisen Precise Point Positioning (PPP) ergab, dass das kombinierte Produkt mit den Residuen der ACs mit Werten von weniger als 10 mm übereinstimmt. Die MGEX-Uhrenkombination zeigte, dass die meisten ACs unter 80 ps für GPS und Galileo übereinstimmen, und für die anderen Konstellationen sind die Werte spärlicher. Darüber hinaus ist festzustellen, dass SHA und JAX im Vergleich zu den anderen ACs große Diskrepanzen aufweisen. Die Validierung erfolgt wiederum im Vergleich mit der offiziellen IGS-Kombination, bei der die Übereinstimmung mit den vorgestellten Lösungen bei etwa 32 ps liegt. Eine PPP-Verarbeitung wird ebenfalls unter Verwendung der kombinierten Lösungen durchgeführt, was die Eignung der Produkte zeigt. Die ACs-Produkte aus der dritten IGS-Reprocessing-Kampagne werden ebenfalls als Input für die Kombination verwendet und diskutiert. Schließlich wird eine Erweiterung der Orbit-Kombination für die Kombination von Low Earth Orbit (LEO)-Satelliten durchgeführt. Die Ergebnisse zeigen, dass die vorgeschlagene Kombination zuverlässige Produkte gemäß den Standards des IGS liefern kann.

Description: Watershed management requires an understanding of key hydrochemical processes. The Pra Basin is one of the five major river basins in Ghana with a population of over 4.2 million people. Currently, water resources management faces challenges due to surface water pollution caused by the unregulated release of untreated household and industrial waste into aquatic ecosystems and illegal mining activities. This has increased the need for groundwater as the most reliable water supply. Our understanding of groundwater recharge mechanisms and chemical evolution in the basin has been inadequate, making effective management difficult. Therefore, the main objective of this work is to gain insight into the processes that determine the hydrogeochemical evolution of groundwater quality in the Pra Basin. The combined use of stable isotope, hydrochemistry, and water level data provides the basis for conceptualizing the chemical evolution of groundwater in the Pra Basin. For this purpose, the origin and evaporation rates of water infiltrating into the unsaturated zone were evaluated. In addition, Chloride Mass Balance (CMB) and Water Table Fluctuations (WTF) were considered to quantify groundwater recharge for the basin. Indices such as water quality index (WQI), sodium adsorption ratio (SAR), Wilcox diagram, and salinity (USSL) were used in this study to determine the quality of the resource for use as drinking water and for irrigation purposes. Due to the heterogeneity of the hydrochemical data, the statistical techniques of hierarchical cluster and factor analysis were applied to subdivide the data according to their spatial correlation. A conceptual hydrogeochemical model was developed and subsequently validated by applying combinatorial inverse and reaction pathway-based geochemical models to determine plausible mineral assemblages that control the chemical composition of the groundwater. The interactions between water and rock determine the groundwater quality in the Pra Basin. The results underline that the groundwater is of good quality and can be used for drinking water and irrigation purposes. It was demonstrated that there is a large groundwater potential to meet the entire Pra Basin’s current and future water demands. The main recharge area was identified as the northern zone, while the southern zone is the discharge area. The predominant influence of weathering of silicate minerals plays a key role in the chemical evolution of the groundwater. The work presented here provides fundamental insights into the hydrochemistry of the Pra Basin and provides data important to water managers for informed decision-making in planning and allocating water resources for various purposes. A novel inverse modelling approach was used in this study to identify different mineral compositions that determine the chemical evolution of groundwater in the Pra Basin. This modelling technique has the potential to simulate the composition of groundwater at the basin scale with large hydrochemical heterogeneity, using average water composition to represent established spatial groupings of water chemistry.

Description: Die Bewirtschaftung von Wassereinzugsgebieten erfordert ein Verständnis der wichtigsten hydrochemischen Prozesse. Das Pra-Becken ist eines der fünf großen Flusseinzugsgebiete Ghanas mit einer Bevölkerung von über 4,2 Millionen Menschen. Die Bewirtschaftung der Wasserressourcen wird derzeit durch die Verschmutzung der Oberflächengewässer erschwert, die durch die unkontrollierte Einleitung von unbehandelten Haushalts- und Industrieabfällen in die aquatischen Ökosysteme und durch illegale Bergbauaktivitäten entsteht. Dies hat den Bedarf an Grundwasser als zuverlässigste Wasserversorgung erhöht. Unser Verständnis der Mechanismen der Grundwasserneubildung und der chemischen Entwicklung im Einzugsgebiet ist bislang unzureichend, was eine wirksame Bewirtschaftung erschwert. Daher ist das Hauptziel dieser Arbeit Einblicke in die Prozesse zu bekommen, welche die hydrogeochemische Entwicklung der Grundwasserqualität im Pra-Becken bestimmen. Die kombinierte Verwendung von Daten stabiler Isotope, der Hydrochemie und von Wasserständen bildet die Grundlage für die Konzeption der chemischen Entwicklung des Grundwassers im Pra-Becken. Dafür wurden die Herkunft und die Verdunstungsraten des in die ungesättigte Zone infiltrierenden Wassers bewertet. Darüber hinaus wurden die Chlorid-Massenbilanz und die Wasserspiegelschwankungen betrachtet, um die Grundwasserneubildung für das Einzugsgebiet zu quantifizieren. Indizes wie der Wasserqualitätsindex (WQI), das Natriumadsorptionsverhältnis (SAR), das Wilcox-Diagramm und der Salzgehalt (USSL) wurden in dieser Studie verwendet, um die Qualität der Ressource für die Verwendung als Trinkwasser und zu Bewässerungszwecken zu bestimmen. Aufgrund der Heterogenität der hydrochemischen Daten wurden die statistischen Verfahren der hierarchischen Cluster- und Faktorenanalyse angewandt, um die Daten entsprechend ihrer räumlichen Korrelation zu unterteilen. Ein konzeptionelles hydrogeochemisches Modell wurde entwickelt und anschließend durch Anwendung kombinatorischer inverser und reaktionspfadbasierter geochemischer Modelle validiert, um plausible mineralische Assemblagen zu bestimmen, welche die chemische Zusammensetzung des Grundwassers kontrollieren. Die Wechselwirkungen zwischen Wasser und Gestein bestimmen die Grundwasserqualität im Pra-Becken. Die Ergebnisse unterstreichen, dass das Grundwasser eine gute Qualität aufweist und als Trinkwasser und für Bewässerungszwecke genutzt werden kann. Es wurde nachgewiesen, dass ein großes Grundwasserpotenzial vorhanden ist, um den derzeitigen und künftigen Wasserbedarf des gesamten Pra-Beckens zu decken. Als Hauptneubildungsgebiet wurde die nördliche Zone im Gebiet identifiziert, während die südliche Zone das Abflussgebiet ist. Der vorherrschende Einfluss der Verwitterung von Silikatmineralen spielt bei der chemischen Entwicklung des Grundwassers eine zentrale Rolle. Die hier vorgestellte Arbeit gibt grundlegende Einblicke in die Hydrochemie des Pra-Beckens und liefert für das Wassermanagement wichtige Daten für eine fundierte Entscheidungsfindung bei der Planung und Zuweisung von Wasserressourcen für verschiedene Zwecke. Ein neuartiger Ansatz zur inversen Modellierungwurde in dieser Studie eingesetzt, um unterschiedliche Mineralzusammensetzungen zu ermitteln, welche die chemische Entwicklung des Grundwassers im Pra-Becken bestimmen. Diese Modellierungstechnik hat das Potenzial, die Zusammensetzung eines Grundwassers auf der Skala eines Beckens mit großer hydrochemischer Heterogenität zu simulieren, wobei die durchschnittliche Wasserzusammensetzung zur Darstellung der etablierten räumlichen Gruppierungen der Wasserchemie verwendet wird.

Description: Watershed management requires an understanding of key hydrochemical processes. The Pra Basin is one of the five major river basins in Ghana with a population of over 4.2 million people. Currently, water resources management faces challenges due to surface water pollution caused by the unregulated release of untreated household and industrial waste into aquatic ecosystems and illegal mining activities. This has increased the need for groundwater as the most reliable water supply. Our understanding of groundwater recharge mechanisms and chemical evolution in the basin has been inadequate, making effective management difficult. Therefore, the main objective of this work is to gain insight into the processes that determine the hydrogeochemical evolution of groundwater quality in the Pra Basin. The combined use of stable isotope, hydrochemistry, and water level data provides the basis for conceptualizing the chemical evolution of groundwater in the Pra Basin. For this purpose, the origin and evaporation rates of water infiltrating into the unsaturated zone were evaluated. In addition, Chloride Mass Balance (CMB) and Water Table Fluctuations (WTF) were considered to quantify groundwater recharge for the basin. Indices such as water quality index (WQI), sodium adsorption ratio (SAR), Wilcox diagram, and salinity (USSL) were used in this study to determine the quality of the resource for use as drinking water and for irrigation purposes. Due to the heterogeneity of the hydrochemical data, the statistical techniques of hierarchical cluster and factor analysis were applied to subdivide the data according to their spatial correlation. A conceptual hydrogeochemical model was developed and subsequently validated by applying combinatorial inverse and reaction pathway-based geochemical models to determine plausible mineral assemblages that control the chemical composition of the groundwater. The interactions between water and rock determine the groundwater quality in the Pra Basin. The results underline that the groundwater is of good quality and can be used for drinking water and irrigation purposes. It was demonstrated that there is a large groundwater potential to meet the entire Pra Basin’s current and future water demands. The main recharge area was identified as the northern zone, while the southern zone is the discharge area. The predominant influence of weathering of silicate minerals plays a key role in the chemical evolution of the groundwater. The work presented here provides fundamental insights into the hydrochemistry of the Pra Basin and provides data important to water managers for informed decision-making in planning and allocating water resources for various purposes. A novel inverse modelling approach was used in this study to identify different mineral compositions that determine the chemical evolution of groundwater in the Pra Basin. This modelling technique has the potential to simulate the composition of groundwater at the basin scale with large hydrochemical heterogeneity, using average water composition to represent established spatial groupings of water chemistry.