ALBERT

Description: The Paleoproterozoic tectono-magmatic history of the Indian plate has been modelled mostly based on investigations of the Indian parts with limited studies considering the Pakistani domains of the mountain range. In this study we investigate basement rocks and younger intrusions of the Indian crust from within the core of the Indus Syntaxis and the Nanga Parbat Syntaxis in northern Pakistan, with the goal of establishing a chronology for the entire tectono-magmatic progression of northwestern Indian plate margin, including the source rocks that were ultimately involved in the formation of younger magmatic pulses. We present in-situ zircon U-Pb geochronology and O-isotopes complemented by whole rock geochemistry of granitoids and paragneisses that yielded evidence for two distinct, large-scale Paleoproterozoic magmatic events that took place in ∼ 1860 Ma and ∼ 2200 Ma in the Indian crust. Later, leucogranites (6.4 ± 0.1 Ma) intruded into the Nanga Parbat basement that were the result of melting of the basement induced by rapid uplift. The δ18OSMOW compositions of zircons from our basement samples range between 7.5 to 9.1 ‰, indicating the derivation of parental magma from the crustal source while the younger leucogranite is somewhat lighter with δ18O ranging between 7.06 to 8.23 ‰. Geochemical data show extensional tectonic settings for the basement rocks that have an A-type affinity, again pointing to a crustal precursor. We propose little to no δ18O exchange took place between the young anatectic melt and basement rocks during crustal evolution. Moreover, geochemical signatures record the crystallization of basement rocks dominated the northwestern margin of Indian plate in a post-orogenic setting in an overall extensional tectonic realm. Based on these observations, we infer that the northwestern Indian margin has experienced a large-scale magmatism during Paleoproterozoic associated with the amalgamation of the Columbia supercontinent.

Description: Tropospheric delay and multipath effect are two key errors that are difficult to be accurately corrected in Global Navigation Satellite System (GNSS) Precise Point Positioning Ambiguity Resolution (PPP-AR). The tropospheric residuals due to imperfect modeling can be significant in harsh environments, like in regions with complex terrain and during extreme weather. As the tropospheric delay and multipath effect are coupled in the unmodeled errors, the tropospheric residuals could be misunderstood as multipath. Therefore, accurate correction of the tropospheric delay is crucial for estimating the multipath. However, the coupling effect was not properly considered in previous studies. We propose a refined joint troposphere-multipath hemispherical map (TMM), by constructing a refined troposphere hemispherical map (THM) and an improved multipath hemispherical map (C-TMHM). We use ray-tracing and meteorological data to construct THM correction tables, while tropospheric delays in PPP-AR are corrected by retrieving the corresponding satellite tropospheric delay estimates to obtain "cleaner" multipath C-TMHM model values. Results show that the tropospheric THM model reduces GNSS residuals from about 10–2 mm at low-elevation (7°–30°) compared to the Vienna Mapping Functions 3 (VMF3). Because that the topographic complexity and the rapid variations in atmospheric water vapor are not adequately considered by simply using the elevation-dependent mapping function and horizontal gradients. In particular, part of the tropospheric residuals in the low-elevation are likely to be misinterpreted as multipath. Compared with multi-GNSS PPP-AR positioning performance using traditional model (VMF3 and uncorrected multipath), the proposed TMM (THM and C-TMHM) model improves the positioning accuracy by 32.12% and 36.18% under the cases of complex terrain and extreme weather, respectively, while shortens the convergence time by 33.04% and 30.7%.

Description: Above-Ground Biomass (AGB) is an important parameter in the conservation of mangrove ecosystem owing to their ecological and economic benefits. LiDAR technologies in forest studies have become popular, due to its highly accurate 3D spatial data acquisition. In this study, we propose an end-to-end framework for estimating AGB of mangroves from Terrestrial Laser Scanner (TLS) point clouds. The framework includes pre-processing of data, segmenting the wood and foliage at tree level using Weighted Random Forest (WRF) classifier and constructing Quantitative Structure Model (QSM) of the wooden components to estimate its biomass. The flow was extended to AGB estimation of 33 x 33 m plot by integrating tree level framework. The study also finds a unique solution to estimate the contribution of pneumatophores in the AGB. Segmentation of wood/foliage of tree point cloud using WRF yielded better results with an increment of 15.27 % in Balanced accuracy, 0.2 of Cohen’s Kappa coefficient, and 7.45 % in F1score than RF classifier. AGB estimation of mangroves using our approach using TLS data is 47.54 T/ha which has a mean bias of 0.0044 T/ha and RMS variation of 0.026 T/ ha when compared with the allometric methods. A Breadth-first graph-search segmentation approach was used to count the pneumatophores, aerial roots seen in few mangrove species (R2 = 0.94 with manual counting) and estimate its contribution to AGB of mangroves which is first of its kind using TLS point cloud. This outcome could also aid future studies in modeling the underlying root network and estimating the below-ground biomass.

Description: Geothermal energy is one of the most viable sources of renewable heat. However, the potential risk of induced seismicity associated with geothermal operations may slow down the growth of the geothermal sector. Previous research has led to significant progress in understanding fluidinjection- induced seismicity in geothermal reservoirs. However, an in-depth assessment of thermal effects on the seismic risk was generally considered to be of secondary importance. This study aims to investigate the relative influence of temperature and key geological and operational parameters on the slip tendency of pre-existing faults. This is done through coupled thermo-hydro-mechanical simulations of the injection and production processes in synthetic geothermal reservoir models of the most utilized and potentially exploitable Dutch geothermal reservoir formations: Slochteren sandstone, Delft sandstone and Dinantian limestone. In our study, changes in the slip tendency of a fault can largely be attributed to thermo-elastic effects, which confirms the findings of recent studies linking thermal stresses to induced seismicity. While the direct pore pressure effect on slip tendency tends to dominate over the early phase of the operations, once pore pressure equilibrium is established in a doublet system, it is the additional stress change associated with the growing cold-water front around the injection well that has the greatest influence. Therefore, the most significant increase in the slip tendency was observed when this low-temperature front reached the fault zone. The distance between an injection well and a pre-existing fault thus plays a pivotal role in determining the mechanical stability of a fault. A careful selection of a suitable target formation together with an appropriate planning of the operational parameters is also crucial to mitigate the risk of induced seismicity. Besides the well-known relevance of the in situ stress field and local fault geometry, rock-mechanical properties and operation conditions exert a major influence on induced stress changes and therefore on the fault (re)activation potential during geothermal operations.

Description: The CNSC, the Canadian regulator for the nuclear industry, participated in DECOVALEX-2023 Task G that focuses on the thermo (T) - hydro (H)- mechanical (M) behaviour of rock joints. Joints are omnipresent in rock masses and are planes of weakness in the host rock. When deep geological repositories (DGRs) for radioactive waste are being considered in areas where rock joints are present, the joints could be preferential pathways for radionuclide migration. Therefore, their THM behaviour must be better understood to assess the safety of the DGR. Under different possible internal and external perturbations, a joint can move by shear and dilation. If the joint crosses the emplacement area of a waste container, the heat generated from the waste can itself induce shearing of the joint. Excessive shear movement can in turn lead to failure of the container, resulting in earlier release of radionuclides. Furthermore, dilation that might accompany shear, results in an increase in the joint aperture creating a faster flow path for radionuclide transport. Mathematical models are important tools that need to be developed and employed, in order to assess joint shear and dilation under different loading conditions, such as the heat generated from the emplaced waste. The authors have developed such a mathematical model based on a macroscopic formulation within the framework of elasto-plasticity. It is verified against analytical solutions and validated against shear under constant normal load tests and thermal shearing tests of joints in granite.

Description: Premise The rise of angiosperm-dominated tropical rainforests has been proposed to have occurred shortly after the Cretaceous–Paleogene transition. Paleocene fossil wood assemblages are rare yet provide important data for understanding these forests and whether their wood anatomical features can be used to document the changes that occurred during this transition. Methods We used standard techniques to section 11 fossil wood specimens of Paleocene-age, described the anatomy using standard terminology, and investigated their affinities to present-day taxa. Results We report here the first middle Paleocene fossil wood specimens from Myanmar, which at the time was near the equator and anchored to India. Some fossils share affinities with Arecaceae, Sapindales (Anacardiaceae, Meliaceae) and Moraceae and possibly Fabaceae or Lauraceae. One specimen is described as a new species and genus: Compitoxylon paleocenicum gen. et sp. nov. Conclusions This assemblage reveals the long-lasting presence of these aforementioned groups in South Asia and suggests the early presence of multiple taxa of Laurasian affinity in Myanmar and India. The wood anatomical features of the dicotyledonous specimens reveal that both “modern” and “primitive” features (in a Baileyan scheme) are present with proportions similar to features in specimens from Paleocene Indian localities. Their anatomical diversity corroborates that tropical flora display “modern” features early in the history of angiosperms and that their high diversity remained steady afterward.

Description: Scientific drilling expeditions offer a unique op- portunity to characterize microbial communities in the sub- surface that have long been isolated from the surface. With subsurface microbial biomass being low in general, biologi- cal contamination from the drilling fluid, sample processing, or molecular work is a major concern. To address this, char- acterization of the contaminant populations in the drilling fluid and negative extraction controls are essential for assess- ing and evaluating such sequencing data. Here, rock cores down to 2250 m depth, groundwater-bearing fractures, and the drilling fluid were sampled for DNA to characterize the microbial communities using a broad genomic approach. However, even after removing potential contaminant popu- lations present in the drilling fluid, notorious contaminants were abundant and mainly affiliated with the bacterial order Burkholderiales. These contaminant microorganisms likely originated from the reagents used for isolating DNA despite stringent quality standards during the molecular work. The detection of strictly anaerobic sulfate reducers such as Candi- datus Desulforudis audaxviator suggested the presence of au- tochthonous deep biosphere taxa in the sequenced libraries, yet these clades represented only a minor fraction of the se- quence counts (〈 0.1 %), hindering further ecological inter- pretations. The described methods and findings emphasize the importance of sequencing extraction controls and can support experimental design for future microbiological stud- ies in conjunction with continental drilling operations.

Description: Most of in-situ stress data in the Australian continent comes from wellbore stress analysis in deep hydrocarbon reservoirs, and earthquake focal mechanism solutions near the Australian plate boundaries, where geophysical tools facilitate understanding of the present-day stress patterns. This resulted in a paucity of stress information in many other regions such as the northern Bowen Basin, which is an active mining province, but with low seismicity rates and limited deep petroleum exploration. The mining industry runs several hundred kilometres of image logs annually to characterise geotechnical attributes. These logs provide an image from the borehole wall, which facilitates analysis of stress-related borehole deformations for in-situ stress characterisation. This paper examines the orientation of horizontal in-situ stress using different types of image logs in mine boreholes across the northern Bowen Basin. Analyses of 128 km of image logs in 680 vertical boreholes resulted in the interpretation of 9046 pairs of stress-related indicators including 735 drilling induced fractures and 8311 borehole breakouts. Our comprehensive database comprises 890 quality-ranked data records for the orientation of maximum horizontal stress (SHmax), which makes the Bowen Basin as a basin with the highest data density in the world in terms of quality-ranked stress information according to the World Stress Map. Statistical analysis of SHmax orientation reveals that the mean SHmax orientation in northern Bowen Basin is N018◦ ± 16◦. The results show that this orientation is consistent over long distances, which is in contrast with several eastern Australian basins. This uniform stress pattern agrees well with plate-scale geomechanical model predictions, which further highlights the impact of plate boundary forces in the contemporary stress pattern of this region. Detailed image log investigation did not show any systematic rotation of stress; however, some small-scale stress perturbations were observed in the vicinity of sharp stiffness contrasts and geological structures.

Description: The microbiota is attributed to be important for initial soil formation under extreme climate conditions, but experimental evidence for its relevance is scarce. To fill this gap, we investigated the impact of in situ microbial communities and their interrelationship with biocrust and plants compared to abiotic controls on soil formation in initial arid and semiarid soils. Additionally, we assessed the response of bacterial communities to climate change. Topsoil and subsoil samples from arid and semiarid sites in the Chilean Coastal Cordillera were incubated for 16 weeks under diurnal temperature and moisture variations to simulate humid climate conditions as part of a climate change scenario. Our findings indicate that microorganism-plant interaction intensified aggregate formation and stabilized soil structure, facilitating initial soil formation. Interestingly, microorganisms alone or in conjunction with biocrust showed no discernible patterns compared to abiotic controls, potentially due to watermasking effects. Arid soils displayed reduced bacterial diversity and developed a new community structure dominated by Proteobacteria, Actinobacteriota, and Planctomycetota, while semiarid soils maintained a consistently dominant community of Acidobacteriota and Proteobacteria. This highlighted a sensitive and specialized bacterial community in arid soils, while semiarid soils exhibited a more complex and stable community. We conclude that microorganism-plant interaction has measurable impacts on initial soil formation in arid and semiarid regions on short time scales under climate change. Additionally, we propose that soil and climate legacies are decisive for the present soil microbial community structure and interactions, future soil development, and microbial responses.

Description: Anthropogenic climate change drives extreme weather events, leading to significant consequences for both society and the environment. This includes damage to road infrastructure, causing disruptions in transportation, obstructing access to emergency services, and hindering humanitarian organizations after natural disasters. In this study, we develop a novel method for analyzing the impacts of natural hazards on transportation networks rooted in the gravity model of travel, offering a fresh perspective to assess the repercussions of natural hazards on transportation network stability. Applying this approach to the Ahr valley flood of 2021, we discovered that the destruction of bridges and roads caused major bottlenecks, affecting areas considerably distant from the flood’s epicenter. Furthermore, the flood-induced damage to the infrastructure also increased the response time of emergency vehicles, severely impeding the accessibility of emergency services. Our findings highlight the need for targeted road repair and reinforcement, with a focus on maintaining traffic flow for emergency responses. This research provides a new perspective that can aid in prioritizing transportation network resilience measures to reduce the economic and social costs of future extreme weather events.

Description: Analyzing seismic data in a timely manner is essential for potential eruption forecasting and early warning in volcanology. Here, we demonstrate that unsupervised machine learning methods can automatically uncover hidden details from the continuous seismic signals recorded during Iceland’s 2021 Geldingadalir eruption. By pinpointing the eruption’s primary phases, including periods of unrest, ongoing lava extrusion, and varying lava fountaining intensities, we can effectively chart its temporal progress. We detect a volcanic tremor sequence three days before the eruption, which may signify impending eruptive activities. Moreover, the discerned seismicity patterns and their temporal changes offer insights into the shift from vigorous outflows to lava fountaining. Based on the extracted patterns of seismicity and their temporal variations we propose an explanation for this transition. We hypothesize that the emergence of episodic tremors in the seismic data in early May could be related to an increase in the discharge rate in late April.

Description: In modelling atmospheric loading effects for terrestrial gravimetry, state-of-the-art approaches take advantage of numerical weather models to account for the global 3-D distribution of air masses. Deformation effects are often computed assuming the Inverse Barometer (IB) hypothesis to be generally valid over the oceans. By a revision of the IB assumption and its consequences we show that although the seafloor is not deformed by atmospheric pressure changes, there exists a fraction of ocean mass that current modelling schemes are usually not accounting for. This causes an overestimation of the atmospheric attraction effect over oceans, even when the dynamic response of the ocean to atmospheric pressure and wind is accounted through dynamic ocean models. This signal can reach a root mean square variability of a few nm s−2, depending on the location of the station. We therefore test atmospheric and non-tidal ocean loading effects at five superconducting gravimeter (SG) stations, showing that a better representation of the residual gravity variations is found when Newtonian attraction effects due to the IB response of the ocean are correctly considered. A sliding window variance analysis shows that the main reduction takes place for periods between 5 and 10 d, even for stations far away from the oceans. Since periods of non-tidal ocean mass variability closely resemble atmospheric signals recorded by SGs, we recommend to directly incorporate both an ocean component together with the IB into services that provide weather-related corrections for terrestrial gravimetry.

Description: Existing research indicates that to create geothermal reservoirs using CO2 injection, additional stimulation methods are necessary. N, N-bis(carboxymethyl)-L-glutamic acid (GLDA) injection has been predicted to increase the permeability of CO2 injection-induced cloud-fracture networks (CFNs) and could serve as an additional stimulation method. Nevertheless, the influence of differential stress, flow geometry, and scale on the characteristics of permeability enhancement by GLDA injection is yet to be clarified. Accordingly, this study experimentally elucidated the permeability enhancement characteristics of injecting a chelating agent in fractured granite under differential stress conditions as an additional method for creating geothermal reservoirs using CO2 injection. GLDA injection experiments were conducted on fractured-granite samples under conventional- and true-triaxial stress states under varying differential stress and pH conditions. Regardless of the differential stress and pH conditions, rock deformation and acoustic emission (AE) were negligible during the chelating agent flow-through experiments on the fractured samples, whereas similar permeability enhancement factors were achieved within the same duration. Thus, stress did not affect the permeability enhancement by chelating agent injections. The permeability enhancement factors were inferred to be high near the injection borehole because of the high viscosity of the solution. Therefore, reservoir stimulation should be conducted using low-concentration chelating agent solutions at constant injection pressures. The study provides insights into the stimulation strategies for creating geothermal reservoirs using CO2 injection.

Description: The vapor pressure deficit reflects the difference between how much moisture the atmosphere could and actually does hold, a factor that fundamentally affects evapotranspiration, ecosystem functioning, and vegetation carbon uptake. Its spatial variability and long-term trends under natural versus human-influenced climate are poorly known despite being essential for predicting future effects on natural ecosystems and human societies such as crop yield, wildfires, and health. Here we combine regionally distinct reconstructions of pre-industrial summer vapor pressure deficit variability from Europe’s largest oxygen-isotope network of tree-ring cellulose with observational records and Earth system model simulations with and without human forcing included. We demonstrate that an intensification of atmospheric drying during the recent decades across different European target regions is unprecedented in a pre-industrial context and that it is attributed to human influence with more than 98% probability. The magnitude of this trend is largest in Western and Central Europe, the Alps and Pyrenees region, and the smallest in southern Fennoscandia. In view of the extreme drought and compound events of the recent years, further atmospheric drying poses an enhanced risk to vegetation, specifically in the densely populated areas of the European temperate lowlands.

Description: The accurate estimation of flood probability is crucial for designing water storage and flood retention structures. However, the assumption of identical distribution in flood samples is unrealistic, given the influence of various flood mechanisms. To address this challenge, we proposed a novel framework based on flood clustering and data pooling that encompasses the key steps such as 1) flood event separation based on a peak-detection flood separation algorithm, 2) grouping flood events using the k-prototypes algorithm, 3) application of the UNprecedented Simulated Extreme ENsemble (UNSEEN) approach to pool reforecast ensemble datasets, and 4) statistical mixing approach to derive common quantiles from all the flood groups. We applied the framework to the Dresden gauge in the Elbe River for a detailed case study. Various tests have been performed to assess the applicability of the UNSEEN approach and the reforecast dataset consistently shows the potential for data pooling. The proposed methodology outperformed the classical approach in terms of goodness-of-fit. The relative difference between the classical and the proposed approach ((classical-proposed)/proposed) for the 100-year return level is 0.16, with a reduction in root mean square error (RMSE) value from 163 to 98 m3/s. Further, replication of the approach to the gauges in North Germany exhibited a relative difference ranging from −0.3 to +0.15 and produced better estimates in terms of RMSE compared with the traditional model. In summary, the proposed framework offers a better estimation of flood probability by addressing the inherent sample inhomogeneity along with the inclusion of unprecedented flood samples.

Description: African forest are increasingly in decline as a result of land-use conversion due to human activities. However, a consistent and detailed characterization and mapping of land-use change that results in forest loss is not available at the spatial-temporal resolution and thematic levels suitable for decisionmaking at the local and regional scales; so far they have only been provided on coarser scales and restricted to humid forests. Here we present the first high-resolution (5 m) and continental-scale mapping of land use following deforestation in Africa, which covers an estimated 13.85% of the global forest area, including humid and dry forests. We use reference data for 15 different land-use types from 30 countries and implement an active learning framework to train a deep learning model for predicting land-use following deforestation with an F1-score of 84 ± 0.7 for the whole of Africa. Our results show that the causes of forest loss vary by region. In general, small-scale cropland is the dominant driver of forest loss in Africa, with hotspots in Madagascar and DRC. In addition, commodity crops such as cacao, oil palm, and rubber are the dominant drivers of forest loss in the humid forests of western and central Africa, forming an “arc of commodity crops” in that region. At the same time, the hotspots for cashew are found to increasingly dominate in the dry forests of both western and southeastern Africa, while larger hotspots for large-scale croplands were found in Nigeria and Zambia. The increased expansion of cacao, cashew, oil palm, rubber, and large-scale croplands observed in humid and dry forests of western and south-eastern Africa suggests they are vulnerable to future land-use changes by commodity crops, thus creating challenges for achieving the zero deforestation supply chains, support REDD+ initiatives, and towards sustainable development goals.

Description: Detecting phase arrivals and pinpointing the arrival times of seismic phases in seismograms is crucial for many seismological analysis workflows. For land station data, machine learning methods have already found widespread adoption. However, deep learning approaches are not yet commonly applied to ocean bottom data due to a lack of appropriate training data and models. Here, we compiled an extensive and labeled ocean bottom seismometer (OBS) data set from 15 deployments in different tectonic settings, comprising ∼90,000 P and ∼63,000 S manual picks from 13,190 events and 355 stations. We propose PickBlue, an adaptation of the two popular deep learning networks EQTransformer and PhaseNet. PickBlue joint processes three seismometer recordings in conjunction with a hydrophone component and is trained with the waveforms in the new database. The performance is enhanced by employing transfer learning, where initial weights are derived from models trained with land earthquake data. PickBlue significantly outperforms neural networks trained with land stations and models trained without hydrophone data. The model achieves a mean absolute deviation of 0.05 s for P-waves and 0.12 s for S-waves, and we apply the picker on the Hikurangi Ocean Bottom Tremor and Slow Slip OBS deployment offshore New Zealand. We integrate our data set and trained models into SeisBench to enable an easy and direct application in future deployments.

Description: Surface roughness ubiquitously prevails in natural faults across various length scales. Despite extensive studies highlighting the important role of fault geometry in the dynamics of tectonic earthquakes, whether and how fault roughness affects fluid-induced seismicity remains elusive. Here, we investigate the effects of fault geometry and stress heterogeneity on fluid-induced fault slip and associated seismicity characteristics using laboratory experiments and numerical modeling. We perform fluid injection experiments on quartz-rich sandstone samples containing either a smooth or a rough fault. We find that geometrical roughness slows down injection-induced fault slip and reduces macroscopic slip velocities and fault slip-weakening rates. Stress heterogeneity and roughness control hypocenter distribution, frequency–magnitude characteristics, and source mechanisms of injection-induced acoustic emissions (AEs) (analogous to natural seismicity). In contrast to smooth faults where injection-induced AEs are uniformly distributed, slip on rough faults produces spatially localized AEs with pronounced non-double-couple source mechanisms. We demonstrate that these clustered AEs occur around highly stressed asperities where induced local slip rates are higher, accompanied by lower Gutenberg–Richter b-values. Our findings suggest that real-time monitoring of induced microseismicity during fluid injection may allow identifying progressive localization of seismic activity and improve forecasting of runaway events.

Description: Arctic warming increases the degradation of permafrost soils but little is known about floodplain soils in the permafrost region. This study quantifies soil organic carbon (SOC) and soil nitrogen stocks, and the potential CH4 and CO2 production from seven cores in the active floodplains in the Lena River Delta, Russia. The soils were sandy but highly heterogeneous, containing deep, organic rich deposits with 〉60% SOC stored below 30 cm. The mean SOC stocks in the top 1 m were 12.9 ± 6.0 kg C m−2. Grain size analysis and radiocarbon ages indicated highly dynamic environments with sediment re-working. Potential CH4 and CO2 production from active floodplains was assessed using a 1-year incubation at 20°C under aerobic and anaerobic conditions. Cumulative aerobic CO2 production mineralized a mean 4.6 ± 2.8% of initial SOC. The mean cumulative aerobic:anaerobic C production ratio was 2.3 ± 0.9. Anaerobic CH4 production comprised 50 ± 9% of anaerobic C mineralization; rates were comparable or exceeded those for permafrost region organic soils. Potential C production from the incubations was correlated with total organic carbon and varied strongly over space (among cores) and depth (active layer vs. permafrost). This study provides valuable information on the carbon cycle dynamics from active floodplains in the Lena River Delta and highlights the key spatial variability, both among sites and with depth, and the need to include these dynamic permafrost environments in future estimates of the permafrost carbon-climate feedback.

Description: This contribution is presenting a multidisciplinary investigation of heterogeneities in a clay rock formation, based on seismic tomography, logging and core analysis, as a reconnaissance study for a diffusion experiment. Diffusion experiments in clay rock formations provide crucial experimental data on diffusive transport of radionuclides (RN) in extremely low hydraulic conductivity media. Previous diffusion experiments, conducted, for example, in the Mont Terri underground rock laboratory within the relatively homogeneous shaly facies of Opalinus Clay, and modelling studies of these experiments have demonstrated that the clay rock could sufficiently well be described as a homogeneous anisotropic medium. For other lithofacies, characterized by larger heterogeneity, such simplification may be unsuitable, and the description of heterogeneity over a range of scales will be important. The sandy facies of the Opalinus Clay exhibits a significantly more pronounced heterogeneity compared to the shaly facies, and a combined characterization and RN diffusion study has been initiated to investigate various approaches of heterogeneity characterization and subsequent diffusion in a heterogeneous environment. As an initial step, two inclined exploratory boreholes have been drilled to access the margins of the experiment location. These boreholes have been used to acquire a cross-hole tomographic seismic data set. Optical, natural gamma and backscattering logging were applied and rock cores were analysed. The integrated results of these investigations allowed the identification of an anomalous brighter layer within the investigated area of the sandy facies of approximately 1 m thickness and with its upper bound at roughly 10 m depth within the inclined exploratory wells. Mineralogical analyses revealed only slight variations throughout the rock cores and indicated that the anomalous layer exhibited a slightly higher quartz content, and locally significantly higher calcite contents, accompanied by a lower content of clay minerals. The anomalous layer was characterized by reduced natural gamma emissions, due to the lower clay content, and increased neutron backscattering likely indicating an increased porosity. Seismic P-wave velocities, derived from anisotropic tomography, exhibited a maximal gradient near the top of this layer. The transition from the overlaying darker rock matrix into this layer has been identified as an appropriate location for the setup of a tracer diffusion experiment in a heterogeneous environment.

Description: Accurate assessment of the rate and state friction parameters of rocks is essential for producing realistic earthquake rupture scenarios and, in turn, for seismic hazard analysis. Those parameters can be directly measured on samples, or indirectly based on inversion of coseismic or postseismic slip evolution. However, both direct and indirect approaches require assumptions that might bias the results. Aiming to reduce the potential sources of bias, we take advantage of a downscaled analog model reproducing megathrust earthquakes. We couple the simulated annealing algorithm with quasi-dynamic numerical models to retrieve rate and state parameters reproducing the recurrence time, rupture duration and slip of the analog model, in the ensemble. Then, we focus on how the asperity size and the neighboring segments' properties control the seismic cycle characteristics and the corresponding variability of rate and state parameters. We identify a tradeoff between (a–b) of the asperity and (a–b) of neighboring creeping segments, with multiple parameter combinations that allow mimicking the analog model behavior. Tuning of rate and state parameters is required to fit laboratory experiments with different asperity lengths. Poorly constrained frictional properties of neighboring segments are responsible for uncertainties of (a–b) of the asperity in the order of per mille. Roughly one order of magnitude larger uncertainties derive from asperity size. Those results provide a glimpse of the variability that rate and state friction estimates might have when used as a constraint to model fault slip behavior in nature.

Description: The European Geosciences Union (EGU) brings together geoscientists from all over Europe and the rest of the world, covering all disciplines of Earth, planetary and space sciences. The Division on Energy, Resources and the Environment (ERE), as part of the EGU, follows an interdisciplinary approach to serve society and provide solutions to challenges of our time and in the future. One task for humankind, for example, is to provide adequate and reliable supplies of affordable energy and other resources, obtained in environmentally sustainable ways, which will be essential for economic prosperity, environmental quality and political stability around the world. This volume of Advances in Geosciences spans the range of topics of the division and continues a series of ten ERE special issues over the course of the last ten years. We incorporate emerging topics into the division ERE along the line and we advocate that every idea and opportunity should be studied and tested.

Description: The European Geosciences Union (EGU) brings together geoscientists from all over Europe and the rest of the world, covering all disciplines of Earth, planetary and space sciences. The Division on Energy, Resources and the Environment (ERE), as part of the EGU, follows an interdisciplinary approach to serve society and provide solutions to challenges of our time and in the future. One task for humankind, for example, is to provide adequate and reliable supplies of affordable energy and other resources, obtained in environmentally sustainable ways, which will be essential for economic prosperity, environmental quality and political stability around the world. This volume of Advances in Geosciences spans the range of topics of the division and continues a series of ten ERE special issues over the course of the last ten years. We incorporate emerging topics into the division ERE along the line and we advocate that every idea and opportunity should be studied and tested.

Description: Several studies investigated changes in microbial community composition in thawing permafrost landscapes, but microbial assemblages in the transient ecosystems of the Arctic coastline remain poorly understood. Thermokarst lakes, abrupt permafrost thaw features, are widespread along the pan-Arctic coast and transform into thermokarst lagoons upon coastal erosion and sea-level rise. This study looks at the effect of marine water inundation (imposing a sulfate-rich, saline environment on top of former thermokarst lake sediments) on microbial community composition and the processes potentially driving microbial community assembly. In the uppermost lagoon sediment influenced from marine water inflow, the microbial structures were significantly different from those deeper in the lagoon sediment and from those of the lakes. In addition, they became more similar along depth compared with lake communities. At the same time, the diversity of core microbial consortia community decreased compared with the lake sediments. This work provides initial observational evidence that Arctic thermokarst lake to lagoon transitions do not only substantially alter microbial communities but also that this transition has a larger effect than permafrost thaw and lake formation history.

Description: Many observed time series of precipitation and streamflow show heavy-tail behaviour. For heavy-tailed distributions, the occurrence of extreme events has a higher probability than for distributions with an exponentially receding tail. If we neglect heavy-tail behaviour we might underestimate the magnitude of rarely observed, high-impact events. Robust estimation of upper-tail behaviour is often hindered by the limited length of observational records. Using long time series and a better understanding of the relevant process controls can help with achieving more robust tail estimations. Here, a simulation-based approach is used to analyse the effect of precipitation and runoff generation characteristics on the upper tail of flood peak distributions. Long, synthetic precipitation time series with different tail behaviour are produced by a stochastic weather generator. These are used to force a conceptual rainfall–runoff model. In addition, catchment characteristics linked to a threshold process in the runoff generation are varied between model runs. We characterize the upper-tail behaviour of the simulated precipitation and discharge time series with the shape parameter of the generalized extreme value (GEV) distribution. Our analysis shows that runoff generation can strongly modulate the tail behaviour of flood peak distributions. In particular, threshold processes in the runoff generation lead to heavier tails. Beyond a certain return period, the influence of catchment processes decreases and the tail of the rainfall distribution asymptotically governs the tail of the flood peak distribution. Beyond which return period this is the case depends on the catchment storage in relation to the mean annual rainfall amount.

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: 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: 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: 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: Public earthquake early warning systems (PEEWSs) have the potential to save lives by warning people of incoming seismic waves up to tens of seconds in advance. Given the scale and geographical extent of their impact, this potential is greatest for destructive earthquakes, such as the M7.8 Pazarcik (Türkiye) event of 6 February 2023, which killed almost 60,000 people. However, warning people of imminent strong shaking is particularly difficult for large-magnitude earthquakes because the warning must be given before the earthquake has reached its final size. Here, we show that the Earthquake Network (EQN), the first operational smartphone-based PEEWS and apparently the only one operating during this earthquake, issued a cross-border alert within 12 s of the beginning of the rupture. A comparison with accelerometer and macroseismic data reveals that, owing to the EQN alerting strategy, Turkish and Syrian EQN users exposed to intensity IX and above benefitted from a warning time of up to 58 s before the onset of strong ground shaking. If the alert had been extended to the entire population, approximately 2.7 million Turkish and Syrian people exposed to a lifethreatening earthquake would have received a warning ranging from 30 to 66 s in advance.

Description: This study presents the findings of a splitting analysis conducted on core-refracted teleseismic shear waves (SKS, SKKS and PKS, called together as XKS) and local shear waves, obtained from a dense seismological network spanning the Kamchatka Peninsula. The objective of the study is to examine the pattern of mantle flow beneath the study area through the investigation of seismic anisotropy. The peninsula is situated at the northeastern end of the Kuril–Kamchatka subduction zone, where the Kuril trench intersects with the western boundary of the Aleutian trench. The data set utilized in this study comprises waveform data from a dense network of seismic stations (99 broad-band and short-period stations for the local shear wave splitting analysis and 69 broad-band stations for the SKS splitting analysis). The seismograms were downloaded from publicly available data repositories including the IRIS Data Management Center and the GFZ Data Services (GEOFON program). The dense station coverage allows us to investigate the lateral variations in anisotropy, providing insights into the flow patterns within the mantle. The processing of the combined data sets of local shear wave and teleseismic XKS waves allowed us to partially decipher the source of anisotropy in the mantle. Small delay (splitting) times (∼0.35 s) observed from the local-S data suggest that anisotropy in the mantle wedge is relatively weak with lateral variations. Larger splitting times (∼1.1 s) observed for the XKS waves relative to local S suggest that the main part of splitting on the XKS waves occurs in the subslab mantle. On the other hand, the rotational pattern of seismic anisotropy observed by both the local S and XKS waves suggests the presence of a toroidal flow at the NE edge of the subducting slab, which affects both the mantle wedge and subslab mantle. For the regions away from the edge of the slab, the mantle flow seems to be governed mainly by the drag of the lithospheric plate over the underlying asthenosphere.

Description: The joint European Space Agency and Chinese Academy of Sciences Solar wind Magnetosphere Ionosphere Link Explorer (SMILE) mission will explore global dynamics of the magnetosphere under varying solar wind and interplanetary magnetic field conditions, and simultaneously monitor the auroral response of the Northern Hemisphere ionosphere. Combining these large-scale responses with medium and fine-scale measurements at a variety of cadences by additional ground-based and space-based instruments will enable a much greater scientific impact beyond the original goals of the SMILE mission. Here, we describe current community efforts to prepare for SMILE, and the benefits and context various experiments that have explicitly expressed support for SMILE can offer. A dedicated group of international scientists representing many different experiment types and geographical locations, the Ground-based and Additional Science Working Group, is facilitating these efforts. Preparations include constructing an online SMILE Data Fusion Facility, the discussion of particular or special modes for experiments such as coherent and incoherent scatter radar, and the consideration of particular observing strategies and spacecraft conjunctions. We anticipate growing interest and community engagement with the SMILE mission, and we welcome novel ideas and insights from the solar-terrestrial community.

Description: This dataset comprises event peak flows, representing extreme floods at 516 stations in Germany. The data generation process involves several key steps. Initially, observed rainfall events associated with 10 historical flood disasters from 1950 to 2021 are undergone spatial shifts. These shifts involve three distances (20, 50, and 100 km) and eight directions (North, Northeast, East, Southeast, South, Southwest, West, Northwest), resulting in 24 counterfactual precipitation events. Including the factual (no shift) event, a total of 25 distinct shifting events are considered. Subsequently, these shifted fields are used as atmospheric forcing for a mesoscale hydrological model (mHM) set up and calibrated for the entire Germany. The model produces daily stream flows across its domain, from which the event peak flows are derived. This dataset is expected to provide a valuable resource for analyzing and modeling the dynamics extreme flood events in Germany.

Description: We present SeisMIC, a fast, versatile, and adaptable open-source software to estimate seismic velocity changes from ambient seismic noise. SeisMIC includes a broad set of tools and functions to facilitate end-to-end processing of ambient noise data, from data retrieval and raw data analysis via spectrogram computation, over waveform coherence analysis, to post-processing of the final velocity change estimates. A particular highlight of the software is its ability to invert velocity change time series onto a spatial grid, making it possible to create maps of velocity changes. To tackle the challenge of processing large continuous datasets, SeisMIC can exploit multithreading at high efficiency with an about five-time improvement in compute time compared to MSNoise, probably the most widespread ambient noise software. In this manuscript, we provide a short tutorial and tips for users on how to employ SeisMIC most effectively. Extensive and up-to-date documentation is available online. Its broad functionality combined with easy adaptability and high efficiency make SeisMIC a well-suited tool for studies across all scales.

Description: The dataset is the basis for describing a 60-year-long evolution of groundwater dynamics and thermal field in the North German Basin beneath the Federal State of Brandenburg (NE Germany), covering the period between 1953 and 2014 with monthly increments. It was produced by one-way coupling of a near-surface distributed hydrologic model to a 3D basin-scale thermohydraulic groundwater model with the goal of investigating feedbacks between climate-driven forcing (in terms of time- and space-varying recharge and temperature), basin-scale geology, and topographic gradients. Modeled pressure and temperature distributions are validated against published groundwater level and temperature time series from observation wells. Our results indicate the spatio-temporal extent of the groundwater system subjected to nonlinear interactions between local geological variability and climate conditions. The dataset comprises of input files and scripts required to run the groundwater model in GOLEM and output files from the transient thermo-hydraulic simulations in EXODUS format. The input and output data is organized as separate archived folders (*.gz format).

Description: Assessing the potential and extent of earthquake-induced liquefaction is paramount for seismic hazard assessment, for the large ground deformations it causes can result in severe damage to infrastructure and pose a threat to human lives, as evidenced by many contemporary and historical case studies in various tectonic settings. In that regard, numerical modeling of case studies, using state-of-the-art soil constitutive models and numerical frameworks, has proven to be a tailored methodology for liquefaction assessment. Indeed, these simulations allow for the dynamic response of liquefiable soils in terms of effective stresses, large strains, and ground displacements to be captured in a consistent manner with experimental and in-situ observations. Additionally, the impact of soil properties spatial variability in liquefaction response can be assessed, because the system response to waves propagating are naturally incorporated within the model. Considering that, we highlight that the effect of shear-wave velocity Vs spatial variability has not been thoroughly assessed. In a case study in Metropolitan Concepción, Chile, our research addresses the influence of Vs spatial variability on the dynamic response to liquefaction. At the study site, the 2010 Maule Mw 8.8 megathrust Earthquake triggered liquefaction-induced damage in the form of ground cracking, soil ejecta, and building settlements. Using simulated 2D Vs profiles generated from real 1D profiles retrieved with ambient noise methods, along with a PressureDependentMultiYield03 sand constitutive model, we studied the effect of Vs spatial variability on pore pressure generation, vertical settlements, and shear and volumetric strains by performing effective stress site response analyses. Our findings indicate that increased Vs variability reduces the median settlements and strains for soil units that exhibit liquefaction-like responses. On the other hand, no significant changes in the dynamic response are observed in soil units that exhibit non-liquefaction behavior, implying that the triggering of liquefaction is not influenced by spatial variability in Vs. We infer that when liquefaction-like behavior is triggered, an increase of the damping at the shallowest part of the soil domain might be the explanation for the decrease in the amplitude of the strains and settlements as the degree of Vs variability increases.

Description: Drainage basins delineate Earth's land surface into individual water collection units. Basin shape and river sinuosity determine water and sediment dynamics, affecting landscape evolution and connectivity between ecosystems and freshwater species. However, a high-resolution global dataset for the boundaries and geometry of basins is still missing. Using a 90 m resolution digital elevation model, we measured the areas, lengths, widths, aspect ratios, slopes, and elevations of basins over 50 km2 globally. Additionally, we calculated the lengths and sinuosities of the longest river channels within these 0.67 million basins. We built a new global dataset, Basin90m, to present the basins and rivers, as well as their morphological metrics. To highlight the use cases of Basin90m, we explored the correlations among morphological metrics, such as Hack's law. By comparing with HydroSHEDS, HydroATLAS, and Google Earth images, we demonstrated the high accuracy of Basin90m. Basin90m, available in shapefile format, can be used on various GIS platforms, including QGIS, ArcGIS, and GeoPandas. Basin90m has substantial application prospects in geomorphology, hydrology, and ecology. Basin90m is available at https://doi.org/10.5880/GFZ.4.6.2023.004 (He et al., 2023).

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: The distribution of earthquakes in stable intracontinental tectonic settings is typically far more diffuse than along plate boundaries and the causative mechanisms underlying some recognizable clustering are not understood. Here we show that seismicity in intraplate western and central Europe is largely limited to regions that exhibit a low-density layer in the uppermost lithospheric mantle and preferentially clustered above lateral gradients in upper mantle effective viscosity. The basis for these new insights into the thermal and density configuration of the upper mantle is provided by a shear-wave tomographic model. We propose that the spatial correlations between mantle low-density bodies and crustal seismicity reflect the importance of buoyancy forces within the mantle lithosphere. In addition, under the interaction of forces due to mantle gravitational instabilities, plate tectonics and postglacial rebound, the variably hot and strong mantle lithosphere responds by localized deformation which imposes differential loading on the overlying crust.

Description: Oil-source and oil-oil correlations are challenging when the primary chemical composition of petroleum was altered or destroyed by biodegradation. Biomarkers that are bound to macromolecular matrices rather than being in a free form are largely protected from alteration and may therefore offer a solution to the challenges. Using Microscale Sealed Vessel Catalytic Hydrogenation (MSSV-Hy), the bound biomarkers of two North Sea and two Barents Sea oil samples, with varying levels of biodegradation, were released from their respective asphaltenes. Despite biodegradation having changed the relative concentrations of steranes (C27–C30 steranes) and hopanes (C29–C30 hopanes) in the free biomarker fraction to varying degrees, the asphaltene-bound biomarkers were protected and show pristine ratios unaffected by biodegradation, e.g., the relative content of C27–C29 steranes, the C30/C29 sterane ratio, and the C29/C30 hopane ratio. Empirical calibration then allows these bound biomarkers to be converted into regular biomarker ratios. The MSSV-Hy technique has been shown to be an efficient and effective tool for releasing asphaltene-bound biomarkers from biodegraded oil and revealing genetic information when free biomarkers are altered by biodegradation.

Description: A natural gem quality, inclusion-free, U-poor xenotime has been metasomatically altered at 900 °C and 500–1000 MPa in a series of alkali-bearing fluids including 2 N KOH, 2 N NaOH, Na2Si2O5 + H2O, and NaF + H2O to which UO2 and SiO2 have been added. All experiments were buffered to CO2-CO-graphite with the exception of the NaF + H2O experiment at 900 °C and 500 MPa, which was buffered to magnetite-hematite. With the exception of the Na2Si2O5 + H2O experiment, the xenotime reacted with the fluid tried in all of the remaining experiments to varying degrees. The xenotime showed the greatest degree of reactivity in the two NaF + H2O experiments tried. This reactivity takes the form of a coupled dissolution-reprecipitation reaction, which involves enriching reacted areas of the xenotime in U + Si as well as remobilizing the HREE in these areas inherent to the original xenotime. One result was the formation of alternating bands of Y-rich and HREE- + U-rich bands in the NaF + H2O experiment when the system was oxidized, which are presumed to have formed due a combination of simple chemical zoning coupled with diffusion and coupled dissolution-reprecipitation in the altered areas of the xenotime. A total lack of fluid and mineral inclusions on the nano-scale, as well as no disturbances in the crystal lattice in the reacted areas of the xenotime under HRTEM imaging, indicated that total recrystallization of the reacted areas took place during the course of the NaF + H2O experiment. The results of these experiments has important implications for the metasomatic resetting of the xenotime geochromometer in the presence of alkali-bearing fluids, which would potentially allow metasomatised xenotime to time and record metasomatic events affecting the rock.

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 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: Investigating past interglacial climatic and environmental changes can enhance our understanding of the natural rates and ranges of climate variability under interglacial boundary conditions. However, comparing past interglacial palaeoclimate records from different regions and archives is often complicated by differing and uncertain chronologies. For instance, the duration of the Last Interglacial in Europe is still controversial as southern European palaeoclimate records suggest a duration of ~16 500–18 000 years, while a length of only ~11 000 years in northern-central Europe was previously inferred from the analysis of partly annually laminated (varved) palaeolake sediments recovered at Bispingen, northern Germany. To resolve this discrepancy, we here present sediment microfacies, geochemistry and pollen data from a new sediment core from the Bispingen palaeolake sediment succession, covering the entire Last Interglacial (Eemian) and the earliest part of the Last Glacial (Weichselian). In particular, we provide evidence that the duration of the Last Interglacial at Bispingen must have been hitherto underestimated due to the investigation of an incomplete sediment core. Using microscopic varve counting and sedimentation rate estimates for non-varved sections on the new sediment core, we show that the Eemian in northern-central Europe probably lasted at least ~15 000 years, about 4000 years longer than previously thought. This new duration estimate is in much better agreement with results from southern European palaeoclimate records, clarifying the enigma of a steep trans-European vegetation gradient for several millennia at the end of the Last Interglacial.

Description: Robust chronologies and time equivalent tephra markers are essential to better understand spatial palaeoenvironmental response to past abrupt climatic changes. Identification of well-dated and widely dispersed volcanic ash by tephra and cryptotephra (microscopic volcanic ash) provides time synchronous tie-points and strongly reduces chronological uncertainties. Here, we present the major, minor and trace element analyses of cryptotephra shards in the Dead Sea Deep Drilling sedimentary record (DSDDP 5017-1A) matching the Campanian Ignimbrite (CI). This geochemical identification expands the known dispersal range of the CI to the southeastern Mediterranean, over 2300 km from the volcanic source. Due to the CI eruption occurring near-synchronous with North Atlantic ice surge of Heinrich Event 4 (HE4), this tephra provides insights into regional responses to large-scale climatic change in the Mediterranean. In the Dead Sea, the CI layer is associated with wetter climatic conditions. This contrasts with the contemporaneous occurrence of the CI deposition and dry conditions in the central and eastern Mediterranean suggesting a possible climate time-transgressive expansion of HE4. Our finding underscores the temporal and spatial complexity of regional climate responses and emphasises the importance of tephra as a time marker for studying large-scale climatic changes verses regional variations.

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

Description: Beryllium isotopes have emerged as a quantitative tracer of continental weathering, but accurate and precise determination of the cosmogenic 10Be and stable 9Be in seawater is challenging, because seawater contains high concentrations of matrix elements but extremely low concentrations of 9Be and 10Be. In this study, we develop a new, time-efficient procedure for the simultaneous preconcentration of 9Be and 10Be from (coastal) seawater based on the iron co-precipitation method. The concentrations of 9Be, 10Be, and the resulting 10Be/9Be ratio for Changjiang Estuary water derived from the new procedure agree well with those obtained from the conventional procedure requiring separate preconcentration for 9Be and 10Be determinations. By avoiding the separate preconcentration, our newly developed procedure contributes toward more time-efficient handling of samples, less sample cross-contamination, and a more reliable 10Be/9Be ratio. Prior to this, we validated the iron co-precipitation method using artificial seawater and natural water samples from the Amazon Estuary regarding: (1) the “matrix effect” for Be analysis, (2) its extraction efficiency for pg g−1 levels Be in the presence and absence of organic matter, and (3) the data comparability with another preconcentration method. We calculated that for the determination of 9Be and 10Be in most open ocean seawater with typical 10Be concentrations of 〉 500 atoms g−1, good precisions (〈 5%) can be achieved using less than 3 liters of seawater compared to more than 20 liters routinely used previously. Even for coastal seawater with extremely low 10Be concentration (e.g., 100 atoms g−1), we estimate a maximum amount of 10 liters to be adequate.

Description: The Claromecó foreland Basin (Carboniferous–Permian; southern Buenos Aires province, Argentina) is key to understanding the paleotectonic evolution of the southwestern Gondwana margin and is relevant to energy resource exploration. This study reconstructs the thermal and burial history of the Claromecó Basin by integrating geochemical data, organic petrology, and thermal modeling techniques. Cores samples of the Tunas Formation (Pillahuincó Group, Early Permian) were studied. A 1D thermal model was constructed, calibrated with vitrinite reflectance measurements (VRo %), and corroborated with fluid inclusion and apatite fission track data from previous studies. Rock-Eval pyrolysis results show TOC% values ranging from 0.13 to 60.35 wt%. The Hydrogen index (HI 〈 50 mg HC/g TOC) and Oxygen index (OI 〈 50 mg CO2/g TOC) indicate the dominance of Type III and Type IV kerogens, most likely resulting from the thermal maturation of an original Type III kerogen. Petrologic observations confirm the presence of macerals from the inertinite group, as well as minor amounts of vitrinite and liptinite. The Tmax displays a temperature range mostly from 460 to 610 °C. The VRo % values range from 1.5 to 2.0%. Geochemical data combined with VRo % measurements confirm a late catagenesis to metagenesis stage within the wet to dry gas window for coals and organic-rich strata. In order to constrain the thermal evolution of the basin infill, different scenarios were tested by varying the heat flow and the missing section thickness associated with the uplift and erosion of the basin (Permian–Cenozoic unconformity). The best calibration results were obtained with an erosion thickness of 3000 up to 4200 m and paleo heat flow peaks of either 60 or 80 mW/m2 during the Lower Permian–Lower Cretaceous. The Tunas Formation was deposited and buried during the Permian–Triassic (Gondwanides Orogeny phase), reaching a maximum temperature of 180 °C. The results obtained by combining geochemical analysis, organic petrology, and thermal modeling techniques indicate that the coal beds of the Tunas Formation could have a current potential as gas-prone source rocks. Despite that, the hydrocarbon generation capacity of coal levels is currently low due to the high percentage of residual (Type IV) kerogen. Further research could help clarify if the hydrocarbons potentially expelled by these source rocks