ALBERT

Description: The 9th Data Science Symposium was in Bremen, 4.-5. May 2024

Description: The thermokarst lakes of permafrost regions play a major role in the global carbon cycle. These lakes are sources of methane to the atmosphere although the methane flux is restricted by an ice cover for most of the year. How methane concentrations and fluxes in these waters are affected by the presence of an ice cover is poorly understood. To relate water body morphology, ice formation and methane to each other, we studied the ice of three different water bodies in locations typical of the transition of permafrost from land to ocean in a continuous permafrost coastal region in Siberia. In total, 11 ice cores were analyzed as records of the freezing process and methane composition during the winter season. The three water bodies differed in terms of connectivity to the sea, which affected fall freezing. The first was a bay underlain by submarine permafrost (Tiksi Bay, BY), the second a shallow thermokarst lagoon cut off from the sea in winter (Polar Fox Lagoon, LG) and the third a land-locked freshwater thermokarst lake (Goltsovoye Lake, LK). Ice on all water bodies was mostly methane-supersaturated with respect to atmospheric equilibrium concentration, except for three cores from the isolated lake. In the isolated thermokarst lake, ebullition from actively thawing basin slopes resulted in the localized integration of methane into winter ice. Stable δ13C-CH4 isotope signatures indicated that methane in the lagoon ice was oxidized to concentrations close to or below the calculated atmospheric equilibrium concentration. Increasing salinity during winter freezing led to a micro-environment on the lower ice surface where methane oxidation occurred and the lagoon ice functioned as a methane sink. In contrast, the ice of the coastal marine environment was slightly supersaturated with methane, consistent with the brackish water below. Our interdisciplinary process study shows how water body morphology affects ice formation which mitigates methane fluxes to the atmosphere.

Description: The subglacial landscape of Antarctica records and influences the behaviour of its overlying ice sheet. However, in many places, the evolution of the landscape and its control on ice sheet behaviour have not been investigated in detail. Using recently released radio-echo sounding data, we investigate the subglacial landscape of the Evans–Rutford region of West Antarctica. Following quantitative analysis of the landscape morphology under ice-loaded and ice-unloaded conditions, we identify 10 flat surfaces distributed across the region. Across these 10 surfaces, we identify two distinct populations based on clustering of elevations, which potentially represent remnants of regionally coherent pre-glacial surfaces underlying the West Antarctic Ice Sheet (WAIS). The surfaces are bounded by deeply incised glacial troughs, some of which have potential tectonic controls. We assess two hypotheses for the evolution of the regional landscape: (1) passive-margin evolution associated with the break-up of the Gondwana supercontinent or (2) an extensive planation surface that may have been uplifted in association with either the West Antarctic Rift System or cessation of subduction at the base of the Antarctic Peninsula. We suggest that passive-margin evolution is the most likely of these two mechanisms, with the erosion of glacial troughs adjacent to, and incising, the flat surfaces likely having coincided with the growth of the WAIS. These flat surfaces also demonstrate similarities to other identified surfaces, indicating that a similar formational process may have been acting more widely around the Weddell Sea embayment. The subsequent fluctuations of ice flow, basal thermal regime, and erosion patterns of the WAIS are therefore controlled by the regional tectonic structures.

Description: New records of rare water mites mostly from interstitial habitat of New Zealand are presented. One new genus, Zelandostygolimnochares n. gen. (Piersigiidae) with the type species Z. curtipalpis n. sp. is described. Zelandobates occidentalis Smit & Pešić, 2020 is synonymized with Z. tongariro Smit & Pešić, 2020. The male is for the first time described for Zelandotonia orion Cook, 1992, and the female is for the first time described for Aciculacarus amalis Cook, 1983. New records of a number of rare species are presented.

Description: Type material of some South American species in the Museum für Naturkunde Berlin, belonging to different families is documented. It includes 16 species: Bulimus cuneus L. Pfeiffer, 1854, Bulimus proteus Broderip, 1832, Bulimus scalarioides L. Pfeiffer, 1867, Bulimus similaris J. Moricand, 1856, Andinia (Ehrmanniella) dedicata Weyrauch & Zilch, 1954, Helix hettneriana E. von Martens, 1897, Eurycampta hidalgonis Döring, 1877, Helix aequatoris L. Pfeiffer, 1860, Helix bituberculata L. Pfeiffer, 1853, Helix bourcieri L. Pfeiffer, 1853, Helix neogranadensis L. Pfeiffer, 1845, Cyclostoma (Cyclophorus) bourcieri L. Pfeiffer, 1854, and Helix platygyra Albers, 1857. The concept of the “salvation” of type material is explained.

Description: Here we report a well-preserved isolated physeteroid tooth of Late Miocene age from Liessel, the Netherlands. The presence of several morphological features allows attribution to the macroraptorial physeteroids. Size and morphology are to some extent comparable to Zygophyseter and almost identical to the primarily tooth-based Tortonian taxon Scaldicetus caretti. However, the genus Scaldicetus was declared unutilizable, which is supported here with an overview of modern classifications of Scaldicetus species and specimens. Despite the restrictions, the type species S. caretti is still valid, although the name is to be restricted to the type material. Based on its morphological resemblance, the tooth is identified as Physeteroidea indet. cf. Scaldicetus caretti.

Description: Currently there are several problems with the genus-level taxa within the Scolodontidae. Multiple type-species designations have been proposed for some genera, type species are poorly described, and in some cases no clear decisions have been made in cases of homonymy or synonymy. This has resulted in wrongly identified species and genera within this group, which, among other problems, hinders the discovery and description of new species as well as the identification of known species. This paper is the first in a series in which all scolodontid genera will be redescribed based on type materials, starting with Happia Bourguignat, 1890 and its allies. Nomenclatural issues are resolved where possible. One new genus and a new species are described: Luteostriatella gen. nov. and Austroselenites pichinchense sp. nov. The following new combinations are made: Happia andia (Pilsbry, 1932) comb. nov., Systrophiella altivaga (Crawford, 1939) comb. nov., Systrophiella cayennensis (L. Pfeiffer, 1842) comb. nov., Systrophiella pygmea (Spix in Spix & Wagner, 1827) comb. nov., Systrophiella snethlagei (F. Baker, 1913) comb. nov., Systrophiella vitrina ( J.A. Wagner in Spix & Wagner, 1827) comb. nov., and Luteostriatella variegata (F. Haas, 1949) comb. nov.

Description: Heterozygous mutations in COL10A1 lead to metaphyseal chondrodysplasia type Schmid (MCDS), a skeletal disorder characterized by epiphyseal abnormalities. Prior analysis revealed impaired trimerization and intracellular retention of mutant collagen type X alpha 1 chains as cause for elevated endoplasmic reticulum (ER) stress. However, how ER stress translates into structural defects remained unclear. We generated a medaka (Oryzias latipes) MCDS model harboring a 5 base pair deletion in col10a1, which led to a frameshift and disruption of 11 amino acids in the conserved trimerization domain. col10a1D633a heterozygotes recapitulated key features of MCDS and revealed early cell polarity defects as cause for dysregulated matrix secretion and deformed skeletal structures. Carbamazepine, an ER stress -reducing drug, rescued this polarity impairment and alleviated skeletal defects in col10a1D633a heterozygotes. Our data imply cell polarity dysregulation as a potential contributor to MCDS and suggest the col10a1D633a medaka mutant as an attractive MCDS animal model for drug screening.

Description: The Mekong River Basin: Ecohydrological Complexity from Catchment to Coast, Volume Three presents real facts, data and predictions for quantifying human-induced changes throughout the Mekong watershed, including its estuaries and coasts, and proposes solutions to decrease or mitigate the negative effect and enable sustainable development. This is the first work to link socio–ecological interaction study over the whole Mekong River basin through the lens of ecohydrology. Each chapter is written by a leading expert, with coverage on climate change, groundwater, land use, flooding drought, biodiversity and anthropological issues. Human activities are enormous in the whole watershed and are still increasing throughout the catchment, with severe negative impacts on natural resources are emerging. Among these activities, hydropower dams, especially a series of 11 dams in China, are the most critical as they generate massive changes throughout the system, including in the delta and to the livelihoods of millions of people and they threaten sustainability.

Description: Lower crustal flow in regions of post-orogenic extension has been inferred to explain the exhumation of metamorphic core complexes and associated low-angle normal (detachment) fault systems. However, the origin of detachment faults, whether initially formed as high-angle or low-angle shear zones, and the extension is symmetric or asymmetric remains enigmatic. Here, we use numerical modeling constrained by geophysical and geological data to show that symmetric extension in the central Menderes Massif of western Anatolia is accommodated by the crustal flow. Our geodynamic model explains how opposite dipping Gediz and Büyük Menderes detachment faults are formed by ∼40° footwall rotation. Model predictions agree with seismic tomography data that suggests updoming of lower crust beneath the exhumed massifs, represented as “twin domes” and a flat Moho. Our work helps to account for the genetic relation between the exhumation of metamorphic core complexes and low-angle normal faulting in both Cordillera and Aegean orogenic regions and has important implications on crustal dynamics in extensional provinces.

Description: 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.

Description: Kimberlites are volatile-rich, occasionally diamond-bearing magmas that have erupted explosively at Earth’s surface in the geologic past1,2,3. These enigmatic magmas, originating from depths exceeding 150 km in Earth’s mantle1, occur in stable cratons and in pulses broadly synchronous with supercontinent cyclicity4. Whether their mobilization is driven by mantle plumes5 or by mechanical weakening of cratonic lithosphere4,6 remains unclear. Here we show that most kimberlites spanning the past billion years erupted about 30 million years (Myr) after continental breakup, suggesting an association with rifting processes. Our dynamical and analytical models show that physically steep lithosphere–asthenosphere boundaries (LABs) formed during rifting generate convective instabilities in the asthenosphere that slowly migrate many hundreds to thousands of kilometres inboard of rift zones. These instabilities endure many tens of millions of years after continental breakup and destabilize the basal tens of kilometres of the cratonic lithosphere, or keel. Displaced keel is replaced by a hot, upwelling mixture of asthenosphere and recycled volatile-rich keel in the return flow, causing decompressional partial melting. Our calculations show that this process can generate small-volume, low-degree, volatile-rich melts, closely matching the characteristics expected of kimberlites1,2,3. Together, these results provide a quantitative and mechanistic link between kimberlite episodicity and supercontinent cycles through progressive disruption of cratonic keels.

Description: Soil organic carbon (SOC) distribution and interaction with light is influenced by soil texture parameters (clay, silt and sand), which makes SOC prediction complicated, especially in samples with considerable pedological variability. Hence, understanding the relationship between SOC and soil texture is important within the context of SOC prediction using remote sensing data. The main objective of this study was to find the impact of soil texture on the performance of local SOC prediction models that were developed on Sentinel-2 (S2) multispectral and CASI/SASI (CS) hyperspectral airborne data as the main predictor variables. One approach to that objective was to lowering the texture variance by stratification of the samples. Therefore, soil samples collected from four agricultural sites in the Czech Republic were segregated based on the i) site-based and ii) texture-based stratification strategies. Random forest (RF) models were then developed on all stratified classes with and without considering the soil texture parameters as predictor variables and results were compared with those obtained by the RF models developed on the non-stratified (NS) samples. Both stratification strategies provided more homogeneous classes, which enhanced the accuracy of SOC prediction, compared to using the NS samples. In addition, the texture-based RF models yielded higher accuracy predictions than the site-based ones. Except sand, adding texture parameters to the main predictors improved accuracy of the models, so that the highest prediction performance was obtained by a texture-based model developed on clay added CS data. Overall, texture-based stratification could significantly enhance the SOC prediction, when the texture parameters were added to the S2 and CS data as the main predictor variables.

Description: Sahara sands have been proposed to result from the extensive and repetitive recycling of much older sedimentary rocks -- a necessary mechanism to explain their petrographic maturity and the similarity of their detrital zircon populations at continental scale. Where and how this recycling occurs today remain poorly understood. This study investigates the source of modern sands from the southeastern Sahara by leveraging on a large (n 〉 7800) new dataset of detrital zircon ages from source rocks, modern and ancient dune fields in Chad and Cameroon. We show that zircon age populations show noticeable regional differences when analyzing a large n amount of ages, questioning the similarity of detrital zircon populations in Saharan sands. Dune fields from the driest parts of our sampling area have distinct age distributions that imply discrete sources with differences in bedrock zircon age populations at regional (several 100 km) scale. In the wetter, Sahelian part of our sampling area, the zircon age distribution of dune fields is best-explained by a significant contribution of recent alluvium from local wadis and rivers to the aeolian sedimentary budget. The origin of aeolian sands in the southeastern Sahara is thus local and polygenetic. Recycling of older sedimentary rocks via physical abrasion is only prominent in the driest parts of our sampling area and does not result in the homogenization of Saharan sands.

Description: The present work proposes a simulation-based Bayesian method for parameter estimation and fragility model selection for mutually exclusive and collectively exhaustive (MECE) damage states. This method uses an adaptive Markov chain Monte Carlo simulation (MCMC) based on likelihood estimation using point-wise intensity values. It identifies the simplest model that fits the data best, among the set of viable fragility models considered. The proposed methodology is demonstrated for empirical fragility assessments for two different tsunami events and different classes of buildings with varying numbers of observed damage and flow depth data pairs. As case studies, observed pairs of data for flow depth and the corresponding damage level from the South Pacific tsunami on 29 September 2009 and the Sulawesi–Palu tsunami on 28 September 2018 are used. Damage data related to a total of five different building classes are analysed. It is shown that the proposed methodology is stable and efficient for data sets with a very low number of damage versus intensity data pairs and cases in which observed data are missing for some of the damage levels.

Description: The Mekong Delta is affected by annual seasonal salinity intrusion in December–May. In this season, the discharge of the Mekong is low, and high tidal water levels cause a deep intrusion of saline water. Critical salt concentrations can be observed 50–80 km deep during high tides. These high salt concentrations cause severe problems in terms of fresh irrigation and drinking water. A quantitative knowledge of the current and likely future salinity intrusion is thus urgently required to develop adaptation plans and mitigation measures. This chapter quantifies the current salinity intrusion in the Mekong Delta, as well as a model-based estimation of the possible future salinity intrusion. The different factors affecting salinity intrusion—changing river discharge by climate change and dam development, changes in effective tidal water levels, changing ocean salinity—and their impacts are quantified by hydraulic modeling. Based on the modeling results, recommendations for future planning are derived.

Description: Research data and software are widely accepted as an outcome of scientific work. However, in comparison to text-based publications, there is not yet an established process to assess and evaluate quality of research data and research software publications. This paper presents an attempt to fill this gap. Initiated by the Working Group Open Science of the Helmholtz Association, the Task Group Helmholtz Quality Indicators for Data and Software Publications currently develops a quality indicator for research data and research software publications to be used within the Association. This report summarizes the vision of the group of what all contributes to such an indicator. The proposed approach relies on generic well-established concepts for quality criteria, such as the FAIR Principles and the COBIT Maturity Model. It does – on purpose – not limit itself to technical implementation possibilities to avoid using an existing metric for a new purpose. The intention of this paper is to share the current state for further discussion with all stakeholders, particularly with other groups also working on similar metrics but also with entities that use the metrics.

Description: Observations of rift and rifted margin architecture suggest that significant spatial and temporal structural heterogeneity develops during the multiphase evolution of continental rifting. Inheritance is often invoked to explain this heterogeneity, such as pre‐existing anisotropies in rock composition, rheology, and deformation. Here, we use high‐resolution 3D thermal‐mechanical numerical models of continental extension to demonstrate that rift‐parallel heterogeneity may develop solely through fault network evolution during the transition from distributed to localized deformation. In our models, the initial phase of distributed normal faulting is seeded through randomized initial strength perturbations in an otherwise laterally homogeneous lithosphere extending at a constant rate. Continued extension localizes deformation onto lithosphere‐scale faults, which are laterally offset by 10’s of km and discontinuous along‐strike. These results demonstrate that rift‐ and margin‐parallel heterogeneity of large‐scale fault patterns may in‐part be a natural byproduct of fault network coalescence.

Description: Continental microplates are enigmatic plate boundary features, which can occur in extensional and compressional regimes. Here we focus on microplate formation and their temporal evolution in continental rift settings. To this aim, we employ the geodynamic finite element software ASPECT to conduct 3D lithospheric‐scale numerical models from rift inception to continental breakup. We find that depending on the strike‐perpendicular offset and crustal strength, rift segments connect or interact through one of four regimes: (1) an oblique rift, (2) a transform fault, (3) a rotating continental microplate or (4) a rift jump. We highlight that rotating microplates form at offsets 〉200 km in weak to moderately strong crustal setups. We describe the dynamics of microplate evolution from initial rift propagation, to segment overlap, vertical‐axis rotation, and eventually continental breakup. These models may explain microplate size and kinematics of the Flemish Cap, the Sao Paulo Plateau and other continental microplates that formed during continental rifting worldwide.

Description: Continental rifting is responsible for the generation of major sedimentary basins, both during rift inception and during the formation of rifted continental margins. Geophysical and field studies revealed that rifts feature complex networks of normal faults but the factors controlling fault network properties and their evolution are still matter of debate. Here, we employ high-resolution 2D geodynamic models (ASPECT) including two-way coupling to a surface processes code (FastScape) to conduct 12 models of major rift types that are exposed to various degrees of erosion and sedimentation. We further present a novel quantitative fault analysis toolbox (Fatbox), which allows us to isolate fault growth patterns, the number of faults, and their length and displacement throughout rift history. Our analysis reveals that rift fault networks may evolve through five major phases: 1) distributed deformation and coalescence, 2) fault system growth, 3) fault system decline and basinward localization, 4) rift migration, and 5) breakup. These phases can be correlated to distinct rifted margin domains. Models of asymmetric rifting suggest rift migration is facilitated through both ductile and brittle deformation within a weak exhumation channel that rotates subhorizontally and remains active at low angles. In sedimentation-starved settings, this channel satisfies the conditions for serpentinization. We find that surface processes are not only able to enhance strain localization and to increase fault longevity but that they also reduce the total length of the fault system, prolong rift phases and delay continental breakup.

Description: Over the past few decades, azimuthal seismic anisotropy measurements have been widely used proxy to study past and present‐day deformation of the lithosphere and to characterize convection in the mantle. Beneath continental regions, distinguishing between shallow and deep sources of anisotropy remains difficult due to poor depth constraints of measurements and a lack of regional‐scale geodynamic modeling. Here, we constrain the sources of seismic anisotropy beneath Madagascar where a complex pattern cannot be explained by a single process such as absolute plate motion, global mantle flow, or geology. We test the hypotheses that either Edge‐Driven Convection (EDC) or mantle flow derived from mantle wind interactions with lithospheric topography is the dominant source of anisotropy beneath Madagascar. We, therefore, simulate two sets of mantle convection models using regional‐scale 3‐D computational modeling. We then calculate Lattice Preferred Orientation that develops along pathlines of the mantle flow models and use them to calculate synthetic splitting parameters. Comparison of predicted with observed seismic anisotropy shows a good fit in northern and southern Madagascar for the EDC model, but the mantle wind case only fits well in northern Madagascar. This result suggests the dominant control of the measured anisotropy may be from EDC, but the role of localized fossil anisotropy in narrow shear zones cannot be ruled out in southern Madagascar. Our results suggest that the asthenosphere beneath northern and southern Madagascar is dominated by dislocation creep. Dislocation creep rheology may be dominant in the upper asthenosphere beneath other regions of continental lithosphere.

Description: Silicic caldera volcanoes present major volcanic and seismic hazards but also host dynamic hydrothermal and groundwater systems and a rich but largely unexplored subsurface biosphere. Many of these volcanoes are hosted in rift settings. The intricate connections and feedbacks among magmatism, rifting, hydrothermal processes, and the biosphere in these complex systems remain poorly understood, necessitating subsurface joint observations that are only enabled by scientific drilling. The CALDERA (Connections Among Life, geo-Dynamics and Eruptions in a Rifting Arc caldera) project workshop funded by the International Continental Scientific Drilling Program (ICDP) gathered multi-disciplinary international experts in January 2023 to advance planning of a scientific drilling project within one of these dynamic, rift-hosted calderas, the Okataina Volcanic Centre (OVC), Aotearoa New Zealand. The OVC's high eruption rate, frequent unrest events and earthquake swarms, location in a densely faulted rapidly extending rift, abundant groundwater–geothermal fluid circulations, and diverse surface hot spring microbiota make it an ideal location for exploring a connected geo-hydro-biosphere via scientific drilling and developing a test bed for novel volcano monitoring approaches. Drilling configurations with at least two boreholes (∼ 200 and ∼ 1000–1500 m deep) were favoured to achieve the multi-disciplinary objectives of the CALDERA project. Decadal monitoring including biosphere activity and composition has the potential to evaluate the response of the hydro-bio system to volcano-tectonic activity. In addition to the OVC caldera-scale datasets already available, site surveys will be conducted to select the best drilling locations. The CALDERA project at the OVC would provide, for the first time, an understanding of volcanic–tectonic–hydrological–biological connections in a caldera–rift system and a baseline for global comparisons with other volcanoes, rifts, and hydrothermal systems. CALDERA would serve as an unprecedented model system to understand how and how quickly the subsurface biosphere responds to geologic activities. Discoveries will improve assessment of volcanic and seismic hazards, guide the sustainable management and/or conservation of groundwater and geothermal resources and microbial ecosystems, and provide a forum for interweaving mātauranga Māori and Western knowledge systems.

Description: 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.

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

Description: The 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.

Description: GPS satellite observations indicate that in the tectonically complex eastern Mediterranean and east African regions microplates rotate counterclockwise with respect to the neighboring African plate. Using 3D numerical models, Glerum relates these observations of crustal deformation to the dynamics of the lithosphere and the underlying mantle that may cause this deformation. Glerum first describes her additions to the ASPECT software necessary for numerically modeling the upper mantle and lithosphere dynamics of convergent and divergent plate boundaries. These additions include the tracking of multiple materials with different physical properties and nonlinear viscous as well as viscoplastic rheologies. The implementations of complex, multi-material rheologies are verified with well-known 2D benchmarks and multi-material viscoplasticity is applied in 3D time-dependent thermomechanical models of oceanic subduction. Subsequently, Glerum uses ASPECT to investigate the sensitivity of horizontal surface motions to individual geodynamic processes in the eastern Mediterranean. Identification of all mantle drivers that should participate in modeling attempts to explain observations of crustal flow is essential to fully exploit the information contained by surface motions about their driving processes. Glerum therefore employs 3D data-driven instantaneous dynamics models of compressible flow including a complete set of possible mantle drivers of surface deformation. The reference instantaneous flow model results indicate that mantle processes can explain a large part of the crustal motion of the Aegean-Anatolian microplate. Subsequent systematic perturbations of model properties with respect to this reference model help estimate the individual contributions of tectonic plate motions, slab pull and trench suction, and density-induced mantle flow interacting with the slab and overlying plates while moderated by the mantle’s bulk viscosity. In order of regional importance, the predicted crustal flow of the Aegean-Anatolian region is most sensitive to slab pull, followed by slab-mantle interaction and basal drag, mantle rheology, and the absolute plate motion reference frame. Lastly, Glerum demonstrates a possible mechanism for the counterclockwise rotation of the Victoria microplate in the East African Rift System, which is in striking contrast to the clockwise motion of the surrounding plates. 3D models of the divergent system show that Victoria’s rotation can be caused by the drag of the African and Somalian plates along the strong edges of the microplate, while the rift segments along inherited lithospheric weaknesses facilitate Victoria’s rotation. The amount of rotation is therefore primarily controlled by the distribution of preexisting stronger regions and the weaker Precambrian mobile belts that surround Victoria. The induced counterclockwise rotation of the microplate leads to a clockwise shift of the local extension direction from E-W to more WNW-ESE along the overlapping rift branches. Comparison of the resulting predicted stress field and tectonic regimes to observations helps to elucidate the interpretation of local stress and strain indicators and to reconcile different opening models used to interpret the East African Rift System.

Description: 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.

Description: The lithosphere is often assumed to reside in a thermal steady‐state when quantitatively describing the temperature distribution in continental interiors and sedimentary basins, but also at active plate boundaries. Here, we investigate the applicability limit of this assumption at slowly deforming continental rifts. To this aim, we assess the tectonic thermal imprint in numerical experiments that cover a range of realistic rift configurations. For each model scenario, the deviation from thermal equilibrium is evaluated. This is done by comparing the transient temperature field of every model to a corresponding steady‐state model with identical structural configuration. We find that the validity of the thermal steady‐state assumption strongly depends on rift type, divergence velocity, sample location and depth within the rift. Maximum differences between transient and steady‐state models occur in narrow rifts, at the rift sides, and if the extension rate exceeds 0.5‐2 mm/a. Wide rifts, however, reside close to thermal steady‐state even for high extension velocities. The transient imprint of rifting appears to be overall negligible for shallow isotherms with a temperature less than 100°C. Contrarily, a steady‐state treatment of deep crustal isotherms leads to underestimation of crustal temperatures, especially for narrow rift settings. Thus, not only relatively fast rifts like the Gulf of Corinth, Red Sea, and Main Ethiopian Rift, but even slow rifts like the Kenya Rift, Rhine Graben, and Rio Grande Rift must be expected to feature a pronounced transient component in the temperature field and to therefore violate the thermal steady‐state assumption for deeper crustal isotherms.

Description: We combine numerical modeling of lithospheric extension with analysis of seismic moment release and earthquake b‐value in order to elucidate the mechanism for deep crustal seismicity and seismic swarms in the Main Ethiopian Rift (MER). We run 2D numerical simulations of lithospheric deformation calibrated by appropriate rheology and extensional history of the MER to simulate migration of deformation from mid‐Miocene border faults to ~30 km wide zone of Pliocene to recent rift floor faults. While the highest strain rate is localized in a narrow zone within the rift axis, brittle strain has been accumulated in a wide region of the rift. The magnitude of deviatoric stress shows strong variation with depth. The uppermost crust deforms with maximum stress of 80 MPa, at 8‐14 km depth stress sharply decreases to 10 MPa and then increases to a maximum of 160 MPa at ~18 km depth. These 2 peaks at which the crust deforms with maximum stress of 80 MPa or above correspond to peaks in the seismic moment release. Correspondingly, the drop in stress at 8‐14 km correlates to a low in seismic moment release. At this depth range, the crust is weaker and deformation is mainly accommodated in a ductile manner. We therefore see a good correlation between depths at which the crust is strong and elevated seismic deformation, while regions where the crust is weaker deform more aseismically. Overall the bimodal depth distribution of seismic moment release is best explained by rheology of the deforming crust.

Description: Strike-slip faults are classically associated with pull-apart basins where continental crust is thinned between two laterally offset fault segments. We propose a subsidence mechanism to explain the formation of a new type of basin where no substantial segment offset or syn-strike-slip thinning is observed. Such “flexural strike-slip basins” form due to a sediment load creating accommodation space by bending the lithosphere. We use a two-way coupling between the geodynamic code ASPECT and surface-processes code FastScape to show that flexural strike-slip basins emerge if sediment is deposited on thin lithosphere close to a strike-slip fault. These conditions were met at the Andaman Basin Central fault (Andaman Sea, Indian Ocean), where seismic reflection data provide evidence of a laterally extensive flexural basin with a depocenter located parallel to the strike-slip fault trace.

Description: Due to the difficulties in estimating groundwater recharge and cross-boundary nature of many aquifers, estimating groundwater recharge at large scale has been called upon. Process-based models as well as data-driven models have been established to meet this need. Meanwhile, with the advent of explainable artificial intelligence (XAI) methods, data-driven machine learning models can take advantage of enhanced explainability while keeping the strength of high flexibility. In this study, an ensemble neural network model was built to check the suitability of the model to predict groundwater recharge and the possibility to gain new insights from large data set. Recent large inputs of groundwater recharge data and additional input for the Arabian Peninsula collated in this study were fed to the model with multiple predictors related to climatology considering seasonality, soil and plant characteristics, topography, and hydrogeology. The model showed higher performance (adjusted R2: 0.702, RMSE: 193.35 mm yr−1) than a recent global process-based model in predicting groundwater recharge. Using XAI methods as individual conditional expectations and Shapley Additive Explanation interaction values, the model behavior was analyzed and possible linear and non-linear relationships between the predictors and the groundwater recharge rate were found. Long-term averaged precipitation and enhanced vegetation index showed non-linear relationships with groundwater recharge rate, while slope, compound topographic index, and water table depth showed low importance to the model results. Most model behaviors followed the domain knowledge, while multi-correlation between predictors and data skewness hindered the model from learning.

Description: During February 2023, a total of 32 individual distributed acoustic sensing (DAS) systems acted jointly as a global seismic monitoring network. The aim of this Global DAS Month campaign was to coordinate a diverse network of organizations, instruments, and file formats to gain knowledge and move toward the next generation of earthquake monitoring networks. During this campaign, 156 earthquakes of magnitude 5 or larger were reported by the U.S. Geological Survey and contributors shared data for 60 min after each event’s origin time. Participating systems represent a variety of manufacturers, a range of recording parameters, and varying cable emplacement settings (e.g., shallow burial, borehole, subaqueous, and dark fiber). Monitored cable lengths vary between 152 and 120,129 m, with channel spacing between 1 and 49 m. The data has a total size of 6.8 TB, and are available for free download. Organizing and executing the Global DAS Month has produced a unique dataset for further exploration and highlighted areas of further development for the seismological community to address.

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