ALBERT

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

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

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

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

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

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

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

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

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

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

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

Beschreibung: This article summarizes the ground-motion characterization (GMC) model component of the 2022 New Zealand National Seismic Hazard Model (2022 NZ NSHM). The model development process included establishing a NZ-specific context through the creation of a new ground-motion database, and consideration of alternative ground-motion models (GMMs) that have been historically used in NZ or have been recently developed for global application with or without NZ-specific regionalizations. Explicit attention was given to models employing state-of-the-art approaches in terms of their ability to provide robust predictions when extrapolated beyond the predictor variable scenarios that are well constrained by empirical data alone. We adopted a “hybrid” logic tree that combined both a “weightson- models” approach along with backbone models (i.e., metamodels), the former being the conventional approach to GMC logic tree modeling for NSHM applications using published models, and the latter being increasingly used in research literature and site-specific studies. In this vein, two NZ-specific GMMs were developed employing the backbone model construct. All of the adopted subduction GMMs in the logic tree were further modified from their published versions to include the effects of increased attenuation in the back-arc region; and, all but one model was modified to account for the reduction in ground-motion standard deviations as a result of nonlinear surficial site response. As well as being based on theoretical arguments, these adjustments were implemented as a result of hazard sensitivity analyses using models without these effects, which we consider gave unrealistically high hazard estimates.

Beschreibung: Seismicity usually exhibits a non-Poisson spatiotemporal distribution and could undergo nonstationary processes. However, the Poisson assumption is still deeply rooted in current probabilistic seismic hazard analysis models, especially when input catalogs must be declustered to obtain a Poisson background rate. In addition, nonstationary behavior and scarce earthquake records in regions of low seismicity can bias hazard estimates that use stationary or spatially precise forecasts. In this work, we implement hazard formulations using forecasts that trade-off spatial precision to account for overdispersion and nonstationarity of seismicity in the form of uniform rate zones (URZs), which describe rate variability using non-Poisson probabilistic distributions of earthquake numbers. The impact of these forecasts in the hazard space is investigated by implementing a negative- binomial formulation in the OpenQuake hazard software suite, which is adopted by the 2022 Aotearoa New Zealand National Seismic Hazard Model. For a 10% exceedance probability of peak ground acceleration (PGA) in 50 yr, forecasts that only reduce the spatial precision, that is, stationary Poisson URZ models, cause up to a twofold increase in hazard for low-seismicity regions compared to spatially precise forecasts. Furthermore, the inclusion of non-Poisson temporal processes in URZ models increases the expected PGA by up to three times in low-seismicity regions, whereas the effect on high-seismicity is minimal (∼5%). The hazard estimates presented here highlight the relevance, as well as the feasibility, of incorporating analytical formulations of seismicity that go beyond the inadequate stationary Poisson description of seismicity.

Beschreibung: The distribution of earthquakes in time and space is seldom stationary, which could hinder a robust statistical analysis, particularly in low-seismicity regions with limited data. This work investigates the performance of stationary Poisson and spatially precise forecasts, such as smoothed seismicity models (SSMs), in terms of the available training data. Catalog bootstrap experiments are conducted to: (1) identify the number of training data necessary for SSMs to perform spatially better than the least-informative Uniform Rate Zone (URZ) models; and (2) describe the rate temporal variability accounting for the overdispersion and nonstationarity of seismicity. Formally, the strict-stationarity assumption used in traditional forecasts is relaxed into local and incremental stationarity (i.e., a catalog is only stationary in the vicinity of a given time point t) along with self-similar behavior described by a power law. The results reveal rate dispersion up to 10 times higher than predicted by Poisson models and highlight the impact of nonstationarity in assuming a constant mean rate within training-forecast intervals. The temporal rate variability is translated into a reduction of spatial precision by means of URZmodels. First, counting processes are devised to capture rate distributions, considering the rate as a random variable. Second, we devise a data-driven method based on geodetic strain rate to spatially delimit the precision of URZs, assuming that strain/stress rate is related to the timescales of earthquake interactions. Finally, rate distributions are inferred from the available data within each URZ. We provide forecasts for the New Zealand National Seismic Hazard Model update,which can exhibit rates up to ten times higher in low-seismicity regions compared with SSMs. This study highlights the need to consider nonstationarity in seismicity models and underscores the importance of appropriate statistical descriptions of rate variability in probabilistic seismic hazard analysis.

Beschreibung: National-scale seismic hazard models with large logic trees can be difficult to calculate using traditional seismic hazard software. To calculate the complete 2022 revision of the New Zealand National Seismic Hazard Model—Te Tauira Matapae Pūmate Rū i Aotearoa, including epistemic uncertainty, we have developed a method in which the calculation is broken into two separate stages. This method takes advantage of logic tree structures that comprise multiple, independent logic trees from which complete realizations are formed by combination. In the first stage, we precalculate the independent realizations of the logic trees. In the second stage, we assemble the full ensemble of logic tree realizations by combining components from the first stage. Once all realizations of the full logic tree have been calculated, we can compute aggregate statistics for the model. This method benefits both from the reduction in the amount of computation necessary and its parallelism. In addition to facilitating the computation of a large seismic hazard model, the method described can also be used for sensitivity testing of model components and to speed up experimentation with logic tree structure and weights.

Beschreibung: Central America is a seismically active region where six tectonic plates (North America, Caribbean, Cocos, Nazca, Panama, and South America) interact in a subduction zone with transform faults and two triple points. This complex tectonic setting makes the maximum magnitude—Mmax—estimation a challenging task, with the crustal fault earthquakes being the most damaging in the seismic history of Central America. The empirical source scaling relations (ESSR) allow the Mmax of faults to be determined from rupture parameters. In this study, we use a dataset of well-characterized earthquakes in the region, comprising 64 events from 1972 to 2021 with magnitudes between Mw 4.1 and 7.7. The dataset incorporates records of rupture parameters (length, width, area, slip, and magnitude) and information on the faults and aftershocks associated. This database is an important product in itself, and through its use we determine which global relations fit best to our data via a residual analysis. Moreover, based on the best-quality records, we develop scaling relations for Central America (CA-ESSR) for rupture length, width, and area. These new relations were tested and compared with recent earthquakes, and logic trees are proposed to combine the CA-ESSR and the best-fit global relations. Therefore, we estimate the Mmax for 30 faults using the logic tree for rupture length, considering a total rupture of the fault andmultifault scenarios. Our results suggest that in CentralAmerica rupture areas larger than other regions are required to generate the samemagnitudes.We associate this with the shear modulus (μ), which seems to be lower (∼ 30% less) than the global mean values for crustal rocks. Furthermore, considering multifault ruptures, we found several fault systems with potential Mmax ≥Mw 7.0. These findings contribute to a better understanding of regional seismotectonics and to the efficient characterization of fault rupture models for seismic hazards.

Beschreibung: In near-Earth space, a large population of high-energy electrons are trapped by Earth’s magnetic field. These energetic electrons are trapped in the regions called Earth’s ring current and radiation belts. They are very dynamic and show a very strong dependence on solar wind and geomagnetic conditions. These energetic electrons can be dangerous to satellites in the near-Earth space. Therefore, it is very important to understand the mechanisms which drive the dynamics of these energetic electrons. Wave-particle interaction is one of the most important mechanisms. Among the waves that can be encountered by the energetic electrons when they move around our Earth, whistler mode chorus waves can cause both acceleration and the loss of energetic electrons in the Earth's radiation belts and ring current. Using more than 5 years of wave measurements from NASA’s Van Allen Probe mission, Wang et al (2019) developed chorus wave models which depend on magnetic local time (MLT), Magnetic Latitude (MLat), L-shell, and geomagnetic condition index Kp. To quantify the effect of chorus waves on energetic electrons, we calculated the bounce-averaged quasi-linear diffusion coefficients using the chorus wave model developed by Wang et al (2019) and extended to higher latitudes according to Wang and Shprits (2019). Using these diffusion coefficients, we calculated the lifetime of the electrons with an energy range from 1 keV to 2 MeV. In each MLT, we calculate the lifetime for each energy and L-shell using two different methods according to Shprits et al (2007) and Albert and Shprits (2009). We make the calculated electron lifetime database available here. Please notice that the chorus wave model by Wang et al (2019) is valid when Kp 〈= 6. If the user wants to use this lifetime database for Kp 〉6, please be careful and contact the authors.

Beschreibung: Dynamic rupture simulations generate synthetic waveforms that account for nonlinear source and path complexity. Here, we analyze millions of spatially dense waveforms from 3D dynamic rupture simulations in a novel way to illuminate the spectral fingerprints of earthquake physics. We define a Brune-type equivalent near-field corner frequency (f c ) to analyze the spatial variability of ground-motion spectra and unravel their link to source complexity. We first investigate a simple 3D strike-slip setup, including an asper- ity and a barrier, and illustrate basic relations between source properties and f c varia- tions. Next, we analyze 〉 13,000,000 synthetic near-field strong-motion waveforms generated in three high-resolution dynamic rupture simulations of real earthquakes, the 2019 Mw 7.1 Ridgecrest mainshock, the Mw 6.4 Searles Valley foreshock, and the 1992 Mw 7.3 Landers earthquake. All scenarios consider 3D fault geometries, topography, off-fault plasticity, viscoelastic attenuation, and 3D velocity structure and resolve frequencies up to 1–2 Hz. Our analysis reveals pronounced and localized patterns of elevated f c , specifically in the vertical components. We validate such f c variability with observed near-fault spectra. Using isochrone analysis, we identify the complex dynamic mechanisms that explain rays of elevated f c and cause unexpectedly impulsive, localized, vertical ground motions. Although the high vertical frequencies are also associated with path effects, rupture directivity, and coalescence of multiple rupture fronts, we show that they are dominantly caused by rake-rotated surface-breaking rupture fronts that decel- erate due to fault heterogeneities or geometric complexity. Our findings highlight the potential of spatially dense ground-motion observations to further our understanding of earthquake physics directly from near-field data. Observed near-field f c variability may inform on directivity, surface rupture, and slip segmentation. Physics-based models can identify “what to look for,” for example, in the potentially vast amount of near-field large array or distributed acoustic sensing data.

Beschreibung: Garnet is a prominent mineral in skarn deposits and its rare earth elements (REE) geochemistry is pivotal for understanding skarn mineralization and fluid evolution. In contrast to magmatic and metamorphic garnets, skarn garnets are mainly grossular-andradite in composition. They exhibit variable REE patterns, spanning from notable heavy (H)-REE enrichment to significant light (L)-REE enrichment, accompanied by negative to positive europium (Eu) anomalies. However, the key factors governing REE fractionation in skarn garnets remain uncertain. This study applies the lattice-strain theory (LST) to investigate the influence of crystal chemistry and structure on REE fractionation in garnets from the Lazhushan Fe skarn deposit in eastern China. Our results demonstrate that the garnet-liquid partition coefficient ratios of DLa/DYb significantly increase (up to 5–7 orders of magnitude) with rising andradite content in garnet. This variation underscores the pivotal role of garnet structure in controlling LREE/HREE fractionation. The results further show that partition coefficient ratios of DLa/DSm are strongly dependent on andradite content in garnets, whereas the DGd/DYb ratios only show a weak correlation to the garnet composition. This contrast suggests that fractionation of LREE in garnet is more sensitive to variations of andradite content than HREE. Data compilation of major elements and REE for garnet from the Lazhushan Fe skarn deposit and other skarn deposits worldwide shows that the garnet REE patterns vary from positive through concave to negative shapes with the garnet ranging from grossularitic to andraditic compositions. Such variations in garnet REE patterns are consistent with the results of geochemical modeling based on the LST. This study demonstrates that, through LST equations, the shape of fluid REE patterns can be predicted from garnet REE patterns, and vice versa. Furthermore, the Eu anomaly (Eu/Eu*Grt) in skarn garnet depends mainly on fluid Eu anomaly (Eu/Eu*fluid) and garnet-fluid partition coefficient ratio of D(Eu2+)/D(Eu3+) with the latter being influenced by garnet composition. These findings highlight the critical role of crystal chemistry and structure in garnet REE fractionation, enhancing our ability to utilize garnet REE in tracing the origin and evolution of skarn-forming fluids.

Beschreibung: Tsunamigenic earthquakes pose considerable risks, both economically and socially, yet earthquake and tsunami hazard assessments are typically conducted separately. Earthquakes associated with unexpected tsunamis, such as the 2018 Mw  7.5 strike-slip Sulawesi earthquake, emphasize the need to study the tsunami potential of active submarine faults in different tectonic settings. Here, we investigate physics-based scenarios combining simulations of 3D earthquake dynamic rupture and seismic wave propagation with tsunami generation and propagation. We present time-dependent modeling of one-way linked and 3D fully coupled earthquakes and tsunamis for the ∼ 100 km long Húsavík–Flatey Fault Zone (HFFZ) in North Iceland. Our analysis shows that the HFFZ has the potential to generate sizable tsunamis. The six dynamic rupture models sourcing our tsunami scenarios vary regarding hypocenter location, spatiotemporal evolution, fault slip, and fault structure complexity but coincide with historical earthquake magnitudes. Earthquake dynamic rupture scenarios on a less segmented fault system, particularly with a hypocenter location in the eastern part of the fault system, have a larger potential for local tsunami generation. Here, dynamically evolving large shallow fault slip (∼ 8 m), near-surface rake rotation (± 20∘), and significant coseismic vertical displacements of the local bathymetry (± 1 m) facilitate strike-slip faulting tsunami generation. We model tsunami crest to trough differences (total wave heights) of up to ∼ 0.9 m near the town Ólafsfjörður. In contrast, none of our scenarios endanger the town of Akureyri, which is shielded by multiple reflections within the narrow Eyjafjörður bay and by Hrísey island. We compare the modeled one-way linked tsunami waveforms with simulation results using a 3D fully coupled approach. We find good agreement in the tsunami arrival times and location of maximum tsunami heights. While seismic waves result in transient motions of the sea surface and affect the ocean response, they do not appear to contribute to tsunami generation. However, complex source effects arise in the fully coupled simulations, such as tsunami dispersion effects and the complex superposition of seismic and acoustic waves within the shallow continental shelf of North Iceland. We find that the vertical velocity amplitudes of near-source acoustic waves are unexpectedly high – larger than those corresponding to the actual tsunami – which may serve as a rapid indicator of surface dynamic rupture. Our results have important implications for understanding the tsunamigenic potential of strike-slip fault systems worldwide and the coseismic acoustic wave excitation during tsunami generation and may help to inform future tsunami early warning systems.

Beschreibung: The ocean basins contain numerous volcanic ridges, seamounts and large igneous provinces (LIPs). Numerous studies have focused on the origin of seamount chains and LIPs but much less focus has been applied to understanding the genesis of large volcanic structures formed from a combination or series of volcanic drivers. Here we propose the term Oceanic Mid-plate Superstructures (OMS) to describe independent bathymetric swells or volcanic structures that are constructed through superimposing pulses of volcanism, over long time periods and from multiple sources. These sources can represent periods when the lithosphere drifted over different mantle plumes and/or experienced pulses of volcanism associated with shallow tectonic drivers (e.g. plate flexure; lithospheric extension). Here we focus on the Melanesian Border Plateau (MBP), one example of an OMS that has a complex and enigmatic origin. The MBP is a region of shallow Pacific lithosphere consisting of high volumes of volcanic guyots, ridges and seamounts that resides on the northern edge of the Vitiaz Lineament. Here we reconcile recently published constraints to build a comprehensive volcanic history of the MBP. The MBP was built through four distinct episodes: (1) Volcanism associated with the Louisville hotspot likely generating Robbie Ridge and some Cretaceous seamounts near the MBP. (2) Construction of oceanic islands and seamounts during the Eocene when the lithosphere passed over the Rurutu-Arago hotspot. (3) Reactivation of previous oceanic islands/seamounts and construction of new volcanos in the Miocene when the lithosphere passed over the Samoa hotspot. (4) Miocene to modern volcanism driven by lithospheric deformation and/or westward entrainment of enriched plume mantle due to toroidal mantle flow driven by the rollback of the Pacific plate beneath the Tonga trench. The combination of these processes is responsible for ∼222,000 km2 of intraplate volcanism in the MBP and indicates that this OMS was constructed from multiple volcanic drivers.

Beschreibung: We compared the performance of DREAM3D simulations in reproducing the long‐term radiation belt dynamics observed by Van Allen Probes over the entire year of 2017 with various boundary conditions (BCs) and model inputs. Specifically, we investigated the effects of three different outer boundary conditions, two different low‐energy boundary conditions for seed electrons, four different radial diffusion (RD) coefficients (DLL), four hiss wave models, and two chorus wave models from the literature. Using the outer boundary condition driven by GOES data, our benchmark simulation generally well reproduces the observed radiation belt dynamics inside L* = 6, with a better model performance at lower μ than higher μ, where μ is the first adiabatic invariant. By varying the boundary conditions and inputs, we find that: (a) The data‐driven outer boundary condition is critical to the model performance, while adding in the data‐driven seed population doesn't further improve the performance. (b) The model shows comparable performance with DLL from Brautigam and Albert (2000, https://doi.org/10.1029/1999ja900344), Ozeke et al. (2014, https://doi.org/10.1002/2013ja019204), and Liu et al. (2016, https://doi.org/10.1002/2015gl067398), while with DLL from Ali et al. (2016, https://doi.org/10.1002/2016ja023002) the model shows less RD compared to data. (c) The model performance is similar with data‐based hiss models, but the results show faster loss is still needed inside the plasmasphere. (d) The model performs similarly with the two different chorus models, but better capturing the electron enhancement at higher μ using the Wang et al. (2019, https://doi.org/10.1029/2018ja026183) model due to its stronger wave power, since local heating for higher energy electrons is under‐reproduced in the current model.

Beschreibung: Rare metals (Nb, Ta, Y, Zr, Sn, U, W and REE) are economically important and new supplies need to be found. In order to understand Neoproterozoic rare metal granites of the Arabian–Nubian Shield (ANS), six samples from five rare-metal mineralized alkali feldspar granites, syenogranites and granodiorite from the Central and SE Desert of Egypt were studied for zircon U–Pb ages and O-isotopic compositions as well as whole-rock Sr- and Nd- and alkali feldspar Pb-isotopic compositions. These are transitional between I-type and A-type granites, mostly high-K calc-alkaline, peraluminous granites with gullwing-shaped REE patterns and strongly negative Eu anomalies. Four granites gave mantle-like zircon 18OV-SMOW between 4.2 and 5.96‰ and yielded ages of 628–633 Ma. This is about when subduction-related magmatism began to be replaced by collision-related magmatism. Igla Ahmr granites are older, formed at 691.7–678.9 Ma with 18OV-SMOW c. 5.95‰. All have positive initial Nd values (+3.3 to +6.9) typical for mantle and juvenile crust. Pb isotopic compositions are unusually radiogenic compared with unmineralized ANS granitic rocks. The data indicate similar magmatic sources for ANS mineralized and unmineralized granites. Exploration for other rare-metal mineralized granites in the ANS should focus on bodies with similar characteristics

Beschreibung: The Tieshajie Cu deposit, located in the northeastern part of the Qin-Hang Metallogenic Belt (QHMB), South China, has long been regarded as a representative Meso-Neoproterozoic volcanogenic massive sulfide (VMS) deposit. Here we present a hydrothermal titanite U-Pb age, Re-Os and in-situ S-Cu isotope data for chalcopyrite to constrain the timing and ore genesis of the Tieshajie deposit. Laser ablation-inductively coupled plasma-mass spectrometry (LA-ICP-MS) U-Pb dating of titanite from the disseminated Cu ore yielded a weighted mean 206Pb/238U age of 160.1 ± 4.4 Ma. Chalcopyrite from different ore types has low 187Os/188Os (0.85–3.60) and 187Re/188Os (46.1–614.0) ratios, combined with initial 187Os/188Os (0.74–2.00), excluding a mantle source. A Re-Os isochron age (188 ± 30 Ma) for five chalcopyrite samples is consistent with the titanite U-Pb age within errors. Moreover, the variations in Cu isotope compositions (δ65Cu: −1.13 to +0.12 ‰) and δ34S values (+3.8 to +7.7 ‰) of chalcopyrite are inconsistent with those reported from the ancient VMS deposits in previous studies. Therefore, our results are indicative of a Late Jurassic magmatic-hydrothermal origin instead of a VMS origin for the Tieshajie deposit. In combination with previous studies, we propose that the Tieshajie Cu deposit belongs to the distal part of the Mid-Late Jurassic (170–150 Ma) porphyry-skarn Cu mineralization event in the QHMB, likely triggered by the subduction of the Paleo-Pacific plate during the Late Mesozoic. This study also has new insights into the genesis of Cu mineralization in the QHMB and further provides implications for future exploration.

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

Beschreibung: Low Earth Orbit satellites offer extensive data of the radiation belt region, but utilizing these observations is challenging due to potential contamination and difficulty of intercalibration with spacecraft measurements at Highly Elliptic Orbit that can observe all equatorial pitch-angles. This study introduces a new intercalibration method for satellite measurements of energetic electrons in the radiation belts using a Data assimilation (DA) approach. We demonstrate our technique by intercalibrating the electron flux measurements of the National Oceanic and Atmospheric Administration (NOAA) Polar-orbiting Operational Environmental Satellites (POES) NOAA-15,-16,-17,-18,-19, and MetOp-02 against Van Allen Probes observations from October 2012 to September 2013. We use a reanalysis of the radiation belts obtained by assimilating Van Allen Probes and Geostationary Operational Environmental Satellites observations into 3-D Versatile Electron Radiation Belt (VERB-3D) code simulations via a standard Kalman filter. We compare the reanalysis to the POES data set and estimate the flux ratios at each time, location, and energy. From these ratios, we derive energy and L* dependent recalibration coefficients. To validate our results, we analyze on-orbit conjunctions between POES and Van Allen Probes. The conjunction recalibration coefficients and the data-assimilative estimated coefficients show strong agreement, indicating that the differences between POES and Van Allen Probes observations remain within a factor of two. Additionally, the use of DA allows for improved statistics, as the possible comparisons are increased 10-fold. Data-assimilative intercalibration of satellite observations is an efficient approach that enables intercalibration of large data sets using short periods of data.

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

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

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

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

Beschreibung: Seafloor massive sulfides are modern analogues to ancient volcanogenic massive sulfide deposits, which are particularly enriched in volatile and precious metals (e.g., Te, Au, Ag, Cu, Bi, Se) in subduction-related settings. However, the sources of metals are still poorly constrained, and it remains elusive, whether magmatic volatile influx controls their distribution in submarine hydrothermal systems on the plate tectonic-scale. Here, we demonstrate, for the first time, that Te, As, and Sb contents as well as related Te/As and Te/Sb ratios vary systematically with the δ34S composition of pyrite and native S, as reported by high-resolution coupled SIMS δ34S and trace element LA-ICP-MS micro-analysis. The better correlation of element ratios (Te/As, Te/Sb) opposed to trace element contents (e.g., Te) with δ34S in pyrite demonstrates that element ratios provide a more robust record of magmatic volatile influx than their absolute contents. On this basis, we define a quantitative threshold of high Te/As (〉0.004) and Te/Sb (〉0.6) ratios in pyrite that are indicative of magmatic volatile influx to submarine subduction-related hydrothermal systems. Two-component fluid mixing simulations further revealed that 〈5 % of magmatic volatile influx drastically changes the Te/As (and Te/Sb) ratio of the modelled fluid, but only slightly changes its δ34S composition. This suggests that Te/As and Te/Sb ratios are more sensitive to a magmatic volatile influx into seawater-dominated hydrothermal systems than δ34S signatures if the magmatic volatile influx was low. Beyond this, our results demonstrate that magma-derived fluid mixing with seawater only has a negligible effect on the magmatic volatile record of Te/As and Te/Sb, while the S isotope system is prone for seawater overprinting leading to commonly ambiguous source signatures. Thus, Te/As and Te/Sb systematics in pyrite provide a robust proxy to evaluate the contribution of magmatic volatiles to submarine hydrothermal systems from the grain- to plate tectonic-scale.

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

Beschreibung: Earlier experiments have shown that cyclic hydraulic fracturing (CHF) systematically reduces the monotonic breakdown pressure (MBP). However, cyclic injection also causes a significantly longer injection time to failure as compared to the monotonic injection tests and complex fracture propagation that is hard to predict. In this study, a different injection scheme employing rock fatigue behavior, named creep injection, was tested on granite cylinders. The creep injection creates continuous pressurization under a constant borehole pressure (CBP) with a pre-defined maximum value below the MBP. Three different pressure ratios (CBP/MBP) of 0.85, 0.9 and 0.95 were tested. We found that both the CHF and hydraulic fracturing with creep injection can reduce the breakdown pressure by ca. 15 ~ 20% without confining pressure. Two mechanisms could explain the reduction: the influence of fluid infiltration within the theory of linear poroelasticity and stress corrosion within the subcritical crack growth theory. The lifetime of the granite cores subjected to creep injection is comparable with previous CHF experiments employing the same pressure ratio. In addition, the lifetime increases logarithmically when the ratio of CBP/MBP is decreased. This relationship has a high regression coefficient of R2 = 0.97, and the lifetime can be well predicted using a stress corrosion index of 70. On the contrary, CHF shows a significantly larger variance in the lifetime with a regression coefficient of R2 = 0.19 and, therefore, is hard to predict. Our results also point out that the injection scheme can modify hydraulic fracture patterns, in terms of fracture aperture, branching, and fracture propagation.

Beschreibung: Qongjiagang pegmatite-type Li deposit in Tibet is the first discovered pegmatite-type deposit with economic value in the Himalayan region, which confirms that the Himalayan region has the potential to become a strategic base of rare metal in China, and provide indications to find pegmatite-type Li deposit in the Himalayan region. In this study, we use SEM to identify the type, frequency and occurrence (relationship with cracks) of mineral inclusions in the three main accessory minerals, monazite, apatite and zircon from granite and pegmatite of Qongjiagang Li deposit, combining with the EPMA analysis of feldspar inclusions in apatite to comprehensively trace the property and evolution of the melts and fluids. Our study indicates that: (1) the main mineral inclusions in monazite, apatite and zircon from Qongjiagang Li deposit include silicate, oxide, phosphates and a small amount of sulfide, not only the REE-rich monazite and apatite filling or intersecting cracks are formed by hydrothermal alteration, but also the uraninite and thorianite isolated from cracks occur in the self-irradiation region of zircon are related to fluids; (2) the appearances of columbite and pyrochlore inclusions in the apatite from tourmaline-muscovite granite demonstrate that the initial melt is enriched in Nb and Ta, the amount and type of rare metal mineral inclusions can be used as an indicator for rare metal mineralization in highly evolved granite and pegmatite; (3) the plagioclase inclusions with high and a large range of An values in apatite from spodumene pegmatite represent the capture of less-differentiated melt and continuously fractional crystallization. Our results indicate that the types and compositions of mineral inclusions in accessory minerals can be good tracers for the characteristics and evolution of melts and fluids in the highly evolved granite-pegmatite system

Beschreibung: This publication provides the codes produced for the article "Temporally dynamic carbon dioxide and methane emission factors for rewetted peatlands. Nature Communications Earth and Environment" by Aram Kalhori, Christian Wille, Pia Gottschalk, Zhan Li, Josh Hashemi, Karl Kemper, and Torsten Sachs (https://doi.org/10.1038/s43247-024-01226-9). In the article, the authors estimate the cumulative GHG emissions of a rewetted peatland in Germany using the long-term ecosystem flux measurements. They observe a source-to-sink transition of annual carbon dioxide (CO2) fluxes and decreasing trend of methane (CH4) emissions. This software is written in R and MATLAB. Running the codes ([R files and .m files](Code)) and loading the data files ([CSV files and .mat files](Data)) requires the pre-installation of [R and RStudio] (https://posit.co/downloads/) and ([MATLAB]. The RStudio 2022.07.2 Build 576 version has been used for the R scripts. The land cover classification work was performed in QGIS, v.3.16.11-Hannover. Data were analyzed in both MATLAB and R and plots created with R (R Core Development Team 2020) in RStudio®.

Beschreibung: The data provided here is an exemplary dataset for the flux site Zarnekow from one year (2018). The complete dataset that is needed to run the codes for all the years can be obtained from the European Fluxes Database Cluster under site ID DE-Zrk (Sachs et al., 2016) or provided upon request. This repository is intended to provide the necessary MATLAB and R code to reproduce the results by Kalhori et al. (2024). The data are provided as zip folder containing (1) a csv file with associated definition of variables and units (file: 2023-004_Kalhori-et-al_README_2018_units.txt), (2) a shapefile (file: 2023-004_Kalhori-et-al_2018_LAiV_DOP.shp) and (3) a Geotiff (file: 2023-004_Kalhori-et-al_2018_LAiV_DOP.tiff).

Beschreibung: Hydrothermal alteration is crucial in the formation of many ore deposits, with potassium (K) mobilization and cycling being prevalent. Potassic metasomatism of wall rocks generally forms K-bearing minerals, such as hydrothermal feldspar and mica. However, determining the source and redistribution of K (and other elements transported by the same fluid) in hydrothermal systems is challenging. K isotopes offer a potential solution to this problem. This study presents new K isotope data from two K-rich alteration assemblages — K-feldspar and sericite-quartz-pyrite — in the Jiaodong gold province of China. The data covers a compositional range from unaltered granites to syn-magmatic potassic alteration (formation of K-feldspar) and post-magmatic syn-mineralization phyllic alteration (formation of sericite). Potassic alteration in granite correlates with significant K addition, whereas phyllic alteration of earlier phases of magmatic and hydrothermal K-feldspar resulted in K loss. K-feldspar altered granites display similar δ41K values (–0.55 to –0.42 ‰ for whole-rocks and –0.56 to –0.48 ‰ for K-feldspar separates) as unaltered granite (–0.52 to –0.47 ‰). The narrow δ41K range suggests that magmatic fluid exsolution and magmatic-hydrothermal alteration have a minor effect on δ41K of the altered rock. Phyllic alteration of K-feldspar altered precursor rock leads to K loss and elevated δ41K values ranging from –0.36 to –0.19 ‰ for whole-rocks and –0.34 to –0.17 ‰ for sericite mineral separates. As sericite preferentially incorporates 41K, sericite will have higher δ41K values than the precursor K-feldspar, whereas the fluids will have lower δ41K values. Our study demonstrates that hydrothermal alteration may affect the K isotope composition of altered rocks in several ways, contingent on the nature of the involved phases, making K isotopes a promising tool for studying hydrothermal alteration and associated mineralization.

Beschreibung: Significant progress in permafrost carbon science made over the past decades include the identification of vast permafrost carbon stocks, the development of new pan-Arctic permafrost maps, an increase in terrestrial measurement sites for CO2 and methane fluxes, and important factors affecting carbon cycling, including vegetation changes, periods of soil freezing and thawing, wildfire, and other disturbance events. Process-based modeling studies now include key elements of permafrost carbon cycling and advances in statistical modeling and inverse modeling enhance understanding of permafrost region C budgets. By combining existing data syntheses and model outputs, the permafrost region is likely a wetland methane source and small terrestrial ecosystem CO2 sink with lower net CO2 uptake toward higher latitudes, excluding wildfire emissions. For 2002–2014, the strongest CO2 sink was located in western Canada (median: −52 g C m−2 y−1) and smallest sinks in Alaska, Canadian tundra, and Siberian tundra (medians: −5 to −9 g C m−2 y−1). Eurasian regions had the largest median wetland methane fluxes (16–18 g CH4 m−2 y−1). Quantifying the regional scale carbon balance remains challenging because of high spatial and temporal variability and relatively low density of observations. More accurate permafrost region carbon fluxes require: (a) the development of better maps characterizing wetlands and dynamics of vegetation and disturbances, including abrupt permafrost thaw; (b) the establishment of new year-round CO2 and methane flux sites in underrepresented areas; and (c) improved models that better represent important permafrost carbon cycle dynamics, including non-growing season emissions and disturbance effects.

Beschreibung: This data set includes the results of high-resolution digital elevation models (DEM) and digital image correlation (DIC) analysis applied to analogue modelling experiments. Twenty generic analogue models are extended on top of a rubber sheet. Two benchmark experiments are also reported. Detailed descriptions of the experiments can be found in Liu et al. (submitted) to which this data set is supplement. The data presented here are visualized as topography and the horizontal cumulative surface strain (principal strain and slip rake).

Beschreibung: Rapidly growing seismic and macroseismic databases and simplified access to advanced machine learning methods have in recent years opened up vast opportunities to address challenges in engineering and strong motion seismology from novel, datacentric perspectives. In this thesis, I explore the opportunities of such perspectives for the tasks of ground motion modeling and rapid earthquake impact assessment, tasks with major implications for long-term earthquake disaster mitigation. In my first study, I utilize the rich strong motion database from the Kanto basin, Japan, and apply the U-Net artificial neural network architecture to develop a deep learning based ground motion model. The operational prototype provides statistical estimates of expected ground shaking, given descriptions of a specific earthquake source, wave propagation paths, and geophysical site conditions. The U-Net interprets ground motion data in its spatial context, potentially taking into account, for example, the geological properties in the vicinity of observation sites. Predictions of ground motion intensity are thereby calibrated to individual observation sites and earthquake locations. The second study addresses the explicit incorporation of rupture forward directivity into ground motion modeling. Incorporation of this phenomenon, causing strong, pulse like ground shaking in the vicinity of earthquake sources, is usually associated with an intolerable increase in computational demand during probabilistic seismic hazard analysis (PSHA) calculations. I suggest an approach in which I utilize an artificial neural network to efficiently approximate the average, directivity-related adjustment to ground motion predictions for earthquake ruptures from the 2022 New Zealand National Seismic Hazard Model. The practical implementation in an actual PSHA calculation demonstrates the efficiency and operational readiness of my model. In a follow-up study, I present a proof of concept for an alternative strategy in which I target the generalizing applicability to ruptures other than those from the New Zealand National Seismic Hazard Model. In the third study, I address the usability of pseudo-intensity reports obtained from macroseismic observations by non-expert citizens for rapid impact assessment. I demonstrate that the statistical properties of pseudo-intensity collections describing the intensity of shaking are correlated with the societal impact of earthquakes. In a second step, I develop a probabilistic model that, within minutes of an event, quantifies the probability of an earthquake to cause considerable societal impact. Under certain conditions, such a quick and preliminary method might be useful to support decision makers in their efforts to organize auxiliary measures for earthquake disaster response while results from more elaborate impact assessment frameworks are not yet available. The application of machine learning methods to datasets that only partially reveal characteristics of Big Data, qualify the majority of results obtained in this thesis as explorative insights rather than ready-to-use solutions to real world problems. The practical usefulness of this work will be better assessed in the future by applying the approaches developed to growing and increasingly complex data sets.

Beschreibung: Das rapide Wachstum seismischer und makroseismischer Datenbanken und der vereinfachte Zugang zu fortschrittlichen Methoden aus dem Bereich des maschinellen Lernens haben in den letzen Jahren die datenfokussierte Betrachtung von Fragestellungen in der Seismologie ermöglicht. In dieser Arbeit erforsche ich das Potenzial solcher Betrachtungsweisen im Hinblick auf die Modellierung erdbebenbedingter Bodenerschütterungen und der raschen Einschätzung von gesellschaftlichen Erdbebenauswirkungen, Disziplinen von erheblicher Bedeutung für den langfristigen Erdbebenkatastrophenschutz in seismisch aktiven Regionen. In meiner ersten Studie nutze ich die Vielzahl an Bodenbewegungsdaten aus der Kanto Region in Japan, sowie eine spezielle neuronale Netzwerkarchitektur (U-Net) um ein Bodenbewegungsmodell zu entwickeln. Der einsatzbereite Prototyp liefert auf Basis der Charakterisierung von Erdbebenherden, Wellenausbreitungspfaden und Bodenbeschaffenheiten statistische Schätzungen der zu erwartenden Bodenerschütterungen. Das U-Net interpretiert Bodenbewegungsdaten im räumlichen Kontext, sodass etwa die geologischen Beschaffenheiten in der Umgebung von Messstationen mit einbezogen werden können. Auch die absoluten Koordinaten von Erdbebenherden und Messstationen werden berücksichtigt. Die zweite Studie behandelt die explizite Berücksichtigung richtungsabhängiger Verstärkungseffekte in der Bodenbewegungsmodellierung. Obwohl solche Effekte starke, impulsartige Erschütterungen in der Nähe von Erdbebenherden erzeugen, die eine erhebliche seismische Beanspruchung von Gebäuden darstellen, wird deren explizite Modellierung in der seismischen Gefährdungsabschätzung aufgrund des nicht vertretbaren Rechenaufwandes ausgelassen. Mit meinem, auf einem neuronalen Netzwerk basierenden, Ansatz schlage ich eine Methode vor, umdieses Vorhaben effizient für Erdbebenszenarien aus dem neuseeländischen seismischen Gefährdungsmodell für 2022 (NSHM) umzusetzen. Die Implementierung in einer seismischen Gefährdungsrechnung unterstreicht die Praktikabilität meines Modells. In einer anschließenden Machbarkeitsstudie untersuche ich einen alternativen Ansatz der auf die Anwendbarkeit auf beliebige Erdbebeszenarien abzielt. Die abschließende dritte Studie befasst sich mit dem potenziellen Nutzen der von makroseismischen Beobachtungen abgeleiteten pseudo-Erschütterungsintensitäten für die rasche Abschätzung von gesellschaftlichen Erdbebenauswirkungen. Ich zeige, dass sich aus den Merkmalen solcher Daten Schlussfolgerungen über die gesellschaftlichen Folgen eines Erdbebens ableiten lassen. Basierend darauf formuliere ich ein statistisches Modell, welches innerhalb weniger Minuten nach einem Erdbeben die Wahrscheinlichkeit für das Auftreten beachtlicher gesellschaftlicher Auswirkungen liefert. Ich komme zu dem Schluss, dass ein solches Modell, unter bestimmten Bedingungen, hilfreich sein könnte, um EntscheidungsträgerInnen in ihren Bestrebungen Hilfsmaßnahmen zu organisieren zu unterstützen. Die Anwendung von Methoden des maschinellen Lernens auf Datensätze die sich nur begrenzt als Big Data charakterisieren lassen, qualifizieren die Mehrheit der Ergebnisse dieser Arbeit als explorative Einblicke und weniger als einsatzbereite Lösungen für praktische Fragestellungen. Der praktische Nutzen dieser Arbeit wird sich in erst in Zukunft an der Anwendung der erarbeiteten Ansätze auf wachsende und zunehmend komplexe Datensätze final abschätzen lassen.

Beschreibung: Enhanced knowledge of the Pamir salient formation can contribute to comprehending the tectonic evolution of Himalaya-Tibetan orogen. However, whether the Pamir salient formed along a linear or a curved southern Asian margin between the Tarim and Tajik cratons remains controversial. Likewise, the role of the two craton blocks during the evolution of the Pamir salient is unclear. Here we present three sandbox experiments exploring the effect of the geometry of the southern Asian margin, as well as the presence of Tarim and Tajik cratons. The results show that the highly curved shape of the Pamir salient, transpressional faults in its wings and strike-slip faults within its interior only form along a curved southern Asian margin. A westward-deflecting arcuate thrust wedge formed along the asymmetric curved southern Asian margin. Together with the Tarim craton and the Tajik craton, this wedge facilitated the westward transfer of materials in the Pamir, and resulted in the westward deflection of the velocity field in Pamir and the formation of the Tajik fold-thrust belt. The oblique slip of arcuate thrust wedge along the western edge of the Tarim craton generated the Kongur extensional system. Moreover, the Tarim and Tajik cratons concentrated deformation mainly along the non-cratonic continental margin and promoted the formation of transpressional faults surrounding the Pamir and the strike-slip faults within the Pamir.

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

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

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

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

Beschreibung: Global coupled climate models are in continuous need for evaluation against independent observations to reveal systematic model deficits and uncertainties. Changes in terrestrial water storage (TWS) as measured by satellite gravimetry missions GRACE and GRACE-FO provide valuable information on wetting and drying trends over the continents. Challenges arising from a comparison of observed and modelled water storage trends are related to gravity observations including non-water related variations such as, for example, glacial isostatic adjustment (GIA). Therefore, correcting secular changes in the Earth's gravity field caused by ongoing GIA is important for the monitoring of long-term changes in terrestrial water from GRACE in particular in former ice-covered regions. By utilizing a new ensemble of 56 individual realizations of GIA signals based on perturbations of mantle viscosities and ice history, we find that many of those alternative GIA corrections change the direction of GRACE-derived water storage trends, for example, from gaining mass into drying conditions, in particular in Eastern Canada. The change in the sign of the TWS trends subsequently impacts the conclusions drawn from using GRACE as observational basis for the evaluation of climate models as it influences the dis-/agreement between observed and modelled wetting/drying trends. A modified GIA correction, a combined GRACE/GRACE-FO data record extending over two decades, and a new generation of climate model experiments leads to substantially larger continental areas where wetting/drying trends currently observed by satellite missions coincide with long-term predictions obtained from climate model experiments.

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

Beschreibung: The investigation of key minerals including zircon, apatite, titanite, rutile, monazite, xenotime, allanite, baddeleyite and garnet can retain critical information about petrogenetic and geodynamic processes and may be utilized to understand complex geological histories and the dynamic evolution of the continental crust. They act as small but often robust petrochronological capsules and provide information about crustal evolution, from local processes to plate tectonics and supercontinent cycles. They offer us insights into processes of magmatism, sedimentation, metamorphism and alteration, even when the original protolith is not preserved. In situ techniques have enabled a more in-depth understanding of trace element behaviour in these minerals within their textural context. This has led to more meaningful ages for many stages of geological events. New developments of analytical procedures have further allowed us to expand our petrochronological toolbox while improving precision and accuracy. Combining multiple proxies with multiple minerals has contributed to new interpretations of the crustal history of our planet.

Beschreibung: Silicon isotope fractionation during silicification is poorly understood and impedes our ability to decipher paleoenvironmental conditions from Si isotopes in ancient cherts. To investigate isotope fractionation during silica-for‑carbonate replacement we analyzed the microscale Si and O isotope composition in different silica phases in a silicified zebra dolostone as well as their bulk δ18O and Δ’17O compositions. The subsequent replacement of carbonate layers is mimicked by decreasing δ18O and δ30Si. The textural relationship and magnitude of Si and O isotope fractionation is best explained by near-quantitative silica precipitation in an open system with finite Si. A Rayleigh model for silicification suggests positive Ɛ30/28Si during silicification, conforming with predictions for isotope distribution at chemical equilibrium from ab-initio models. Application of the modelled Ɛ30Si-T relationship yields silicification temperatures of approx. 50 °C. To reconcile the δ18Ochert composition with these temperatures, the δ18O of the fluid must have been between −2.5 and − 4 ‰, compositions for which the quartz phases fall close to the oxygen equilibrium fractionation line in three-isotope space. Diagenetic silica replacement appears to occur in O and Si isotopic equilibrium allowing reconstructions of temperatures of silicification from Si isotopes and derive the δ18O composition of the fluid – a highly desired value needed for accurate reconstructions of the temperature- and δ18O histories of the oceans.

Beschreibung: In this work, IMF By effects on field-aligned currents (FACs) are examined in different local time sectors, seasons, and hemispheres. At dusk and 09–14 MLT, when the eastward polar electrojet (PEJ) prevails, the northern FACp (poleward side FACs) are stronger when IMF By 〈 0 than when IMF By 〉 0. Conversely, at dawn, 21–02 MLT, and 09–14 MLT with westward PEJ, the northern FACp are stronger with IMF By 〉 0 compared to IMF By 〈 0. The southern FACp shows a reversed relationship with IMF By direction. The dependence of FACe (equatorward side FACs) on IMF By is weaker, except for the midday FACe, which shows opposite variations with respect to IMF By when compared to FACp. Stronger IMF By effect is observed in local summer in most of local times. The northern FACs are located at higher latitude for IMF By 〉 0 than for IMF By 〈 0 in local times with eastward PEJ, while the opposite trend is observed in other local times and in the Southern Hemisphere. The hemispheric difference in the peak latitude of FACs demonstrates an inverse relationship with its intensity, with stronger FACs located at lower latitudes. Overall, the local time and hemispheric differences in FACs strength and latitude are discussed in the context of interhemispheric field-aligned currents linked to IMF By.

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

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

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

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

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

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

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