Refine
Year of publication
- 2022 (157) (remove)
Document Type
- Article (125)
- Doctoral Thesis (19)
- Postprint (9)
- Monograph/Edited Volume (2)
- Part of a Book (1)
- Review (1)
Is part of the Bibliography
- yes (157)
Keywords
- climate change (6)
- permafrost (6)
- Andes (3)
- machine learning (3)
- ocean color remote sensing (3)
- radiation belts (3)
- Arctic ocean (2)
- Argentina (2)
- Central Andes (2)
- Diffusion (2)
- ERA5 (2)
- GEDI (2)
- ICESat-2 (2)
- Landsat (2)
- Magnetotellurics (2)
- Opalinus Clay (2)
- Opalinuston (2)
- PHREEQC (2)
- Precipitation (2)
- Sentinel 2 (2)
- Tibetan Plateau (2)
- Vertical flux (2)
- analysis (2)
- bacteria (2)
- bat fatalities (2)
- biodiversity decline (2)
- canopy height (2)
- diffusion (2)
- earthquake (2)
- extreme events (2)
- flood risk (2)
- floods (2)
- food web (2)
- geodynamics (2)
- green-green dilemma (2)
- hydraulic fracturing (2)
- landscape evolution (2)
- lidar (2)
- mapping (2)
- modelling (2)
- monitoring (2)
- nearshore zone (2)
- numerical modelling (2)
- numerische Modellierung (2)
- radiocarbon (2)
- reactive transport (2)
- reaktiver Transport (2)
- renewable energy (2)
- savanna (2)
- sensitivity (2)
- subduction (2)
- surface processes (2)
- suspended sediment (2)
- wind energy production (2)
- %Ro (1)
- 16S rRNA genes (1)
- 2D Numerical Modelling (1)
- ASPECT (1)
- Active tectonics (1)
- Africa (1)
- Ahr (1)
- Ahr River (1)
- Air pollution (1)
- Alexander von Humboldt (1)
- Amazon region (1)
- Analog historical seismograms (1)
- Andean back-arc; (1)
- Anden (1)
- Antarctica (1)
- Arctic (1)
- Arctic Ocean (1)
- Arctic lakes (1)
- Arctic nearhore zone (1)
- Argentinien (1)
- Asia (1)
- Assapaat landslide (1)
- B-isotopes (1)
- Baryt (1)
- Bayesian inference (1)
- Bayesian model (1)
- Beaufort Sea (1)
- Beni Suef Basin (1)
- Best management practice (1)
- Biodiversity (1)
- Blattverschiebung (1)
- Bocono (1)
- Brazil (1)
- Bruchausbreitung (1)
- Building exposure modelling (1)
- CMIP6 (1)
- CO2 emissions (1)
- COH-fluid (1)
- CRNS (1)
- Carbonation (1)
- Causality (1)
- Cenozoic flat-slab (1)
- Chain structure (1)
- Change-point regression model (1)
- Chilean Coastal Cordillera (1)
- Chilean subduction zone (1)
- Climate (1)
- Climate change (1)
- Climate governance (1)
- Cluster analysis moment (1)
- Cologne (1)
- Colorado mineral belt (1)
- Comparative hydrology (1)
- Complex terrain (1)
- Conductive channel (1)
- Continental Rifts (1)
- Coordination number (1)
- Cost-benefit (1)
- DEM analysis (1)
- DOC (1)
- Data collection (1)
- Data-Mining (1)
- Deep carbon cycle (1)
- Deformation (1)
- Deformationsmechanismen (1)
- Densification (1)
- Dissolved organic matter (DOM) (1)
- Distributed acoustic sensing (1)
- EMIC (1)
- EMIC waves (1)
- Early pleistocene; (1)
- Earthquake loss modelling (1)
- Earthquake scenario (1)
- Earthquake source observations (1)
- Earthquake swarm (1)
- Eastern Alps (1)
- Eastern Dharwar craton (1)
- Egypt (1)
- El Nino Southern Oscillation (1)
- Electrical resistivity tomography (1)
- Embedding (1)
- Endlagerung nuklearer Abfälle (1)
- Ensemble (1)
- Epistemic uncertainty (1)
- Erdbeben (1)
- Erosion (1)
- Ethiopia rift (1)
- Etna (1)
- European Alps (1)
- Evolutionary ecology (1)
- Extreme (1)
- Extreme value statistics (1)
- Extremereignisse (1)
- Extremwertstatistik (1)
- Faceted taxonomy (1)
- Fagradalsfjall (1)
- FastScape (1)
- Field measurements (1)
- Fluid inclusion geochemistry (1)
- Fluid mixing (1)
- Fluid/magma reservoir (1)
- Forearc (1)
- Formationsschaden (1)
- Fourier analysis (1)
- Functional connectivity (1)
- GNSS (1)
- Ganga River (1)
- Garnet-orthopyroxene-cordierite gneiss (1)
- Gashydrate (1)
- Gefahrenkarten (1)
- Generalized Extreme Value distribution (1)
- Genetic algorithm (1)
- Geochronology (1)
- Geodynamic Modelling (1)
- Geodynamics (1)
- Geodynamik (1)
- Geodynamische Modellierung (1)
- Geomorphic systems (1)
- Geothermie (1)
- Germany (1)
- Gewässerfernerkundung (1)
- Global inversion (1)
- Global warming (1)
- Global warming potential (1)
- Gorkha earthquake (1)
- Grabenbrüche (1)
- Granulite enclaves (1)
- Graph theory (1)
- Greek Islands (1)
- Greenhouse gas emissions (1)
- Ground-based estimates (1)
- Ground-motion model (1)
- Groundwater recharge (1)
- Grundgestein (1)
- Guaymas Basin (1)
- HBV-light model (1)
- HP-experiments (1)
- Hadley-Walker Circulation (1)
- Heteroscedastic (1)
- High pressure and high temperature; (1)
- Himalaya (1)
- Hochwasserrisiko (1)
- Homo sapiens (1)
- Horizontal (1)
- Horizontal flux (1)
- Hydrology (1)
- Hydrothermal veins (1)
- Iceland (1)
- Image processing (1)
- InSAR (1)
- India (1)
- Indian monsoon (1)
- Indian summer monsoon (1)
- Induzierte Seismizität (1)
- Injektion (1)
- Injektionsschema (1)
- Intermediate-deep earthquakes (1)
- International policy (1)
- Inverse methods (1)
- Ischigualasto-Villa Union Basin (1)
- Isotopenfraktionierung (1)
- James Ross archipelago (1)
- Joint Inversion (1)
- Kenya Rift (1)
- Kerogen petrography (1)
- L-moments estimation (1)
- LP (1)
- Lacustrine sediment (1)
- Landscapes (1)
- Landslide chronosequences (1)
- Landslides (1)
- Las Chacras-Potrerillos (1)
- Limnology (1)
- Lithosphäre (1)
- Lokalisierung von Deformation (1)
- Low impact development (1)
- Lu-Hf (1)
- Machine learning (1)
- Magma-Entgasung (1)
- Magnitude scalin (1)
- Mallik (1)
- Mediterranean (1)
- Merida Andes (1)
- Methane (1)
- Microbial processes (1)
- Micropollutants (1)
- Mikroplatte (1)
- Mikrostrukturelle (1)
- Modellierung (1)
- Modellierung der Wassertrübung (1)
- Molybdenum mineralization (1)
- Mountain meteorology (1)
- Multivariate (1)
- N-butylpyridinium bromide (1)
- NSGA-II (1)
- Nd isotopic composition (1)
- Near-source saturation (1)
- Near-surface geophysics (1)
- Neotektonik (1)
- Nepal (1)
- Nocardioides alcanivorans (1)
- Non-uniqueness (1)
- Northern Alpine Foreland Basin (1)
- Numerische 2D Modellierung (1)
- Nutrients (1)
- ODP 659 (1)
- ODP 721/722 (1)
- ODP 967 (1)
- Oberflächenprozesse (1)
- Ohre Rift (1)
- OpenStreetMap (1)
- Optimization (1)
- Ore (1)
- Orinoco River system (1)
- PM balances (1)
- PM10, PM2, PM1 (1)
- PM10, PM2.5 and PM1.0 concentrations (1)
- PSD (1)
- Palaeoclimate (1)
- Paleoseismologie (1)
- Paläoklimatologie (1)
- Particle swarm optimization (1)
- Permafrostsedimente (1)
- Phytoplankton (1)
- Plio-Pleistocene (1)
- Plio-pleistocene transition (1)
- Prediction (1)
- Probabilistic exposure modelling (1)
- RCP4,5 and 8,5 (1)
- Raman spectroscopy (1)
- Raman-Spektroskopie (1)
- Reaktivierung von Störungszonen (1)
- Recharge controls (1)
- Recurrence analysis (1)
- Respiration (1)
- Reykjanes (1)
- Rhine basin (1)
- Rift (1)
- Risstransmissivität (1)
- Russian Arctic (1)
- SIMS (1)
- Salzgestein (1)
- Sava River (1)
- Scenario (1)
- Schadensmodellierung (1)
- Schallemissionen (1)
- Scheme (1)
- Sediment characterization (1)
- Sediment reuse (1)
- Seesystemreaktionen (1)
- Seismic damage assessment (1)
- Seismic nest (1)
- Seismic noise (1)
- Seismic risk (1)
- Seismicity and tectonics (1)
- Semiarid (1)
- Sensitivity analysis (1)
- Sentinel-1 InSAR (1)
- Short-lived climate pollutants (1)
- Siberia (1)
- Sierra de San Luis (1)
- Silicate melts (1)
- Slope temperature (1)
- Sm-Nd (1)
- Soil (1)
- Sorption (1)
- South America (1)
- South China Sea (1)
- Spannungsmessung (1)
- Spatially cross-correlated ground motion (1)
- Spectral accelerations (1)
- Spectroscopy (1)
- State space reconstruction (1)
- Statistical analysis (1)
- StorAge Selection function (1)
- Storm water management model (1)
- Strain Localisation (1)
- Streamflow (1)
- Strike-slip faults (1)
- Strukturgeologie (1)
- Subduktion (1)
- Surface (1)
- TanDEM-X DEM (1)
- Tauern Window; (1)
- Thick-skinned tectonics (1)
- Time series analysis (1)
- Time-series analysis (1)
- Torsion Experiments (1)
- Torsionsexperimente (1)
- Triassic rifting (1)
- Tropical lake (1)
- Trujillo Block (1)
- Ultra-low velocity zones (1)
- Uran (1)
- VERB (1)
- VLP (1)
- VT events and tremor (1)
- Walker circulation (1)
- Wasser-Gesteins-Wechselwirkungen (1)
- Water budget / balance (1)
- Waveform inversion (1)
- Windböen (1)
- Winderosion (1)
- Wirtsgesteinsskala (1)
- Within-model uncertainty (1)
- XANES (1)
- XRD in the clay fraction (1)
- XRF analysis (1)
- Yangtze River (1)
- Yedoma (1)
- Zentral Asien (1)
- acoustic emissions (1)
- afterslip inversion (1)
- age-heterogeneity (1)
- agemodeling (1)
- air pollution (1)
- ambient noise (1)
- anaerobe Inkubationensexperimente (1)
- anaerobic incubation experiments (1)
- andean geology (1)
- apatite geochemistry and U-Pb dating (1)
- aquatic carbon cycle (1)
- arctic (1)
- arktischer Nahküstenbereich (1)
- bacterial diversity (1)
- bacterial-community structure (1)
- barite (1)
- basement rock (1)
- basin modeling (1)
- biogeography (1)
- bioturbation (1)
- black carbon (1)
- body waves (1)
- bond strength (1)
- boreal forest (1)
- buildings (1)
- burial (1)
- carbon cycle (1)
- catchment (1)
- catchment nitrate export (1)
- catchment response (1)
- change (1)
- characteristics (1)
- chemistry (1)
- chorus waves (1)
- city district (1)
- clay (1)
- cleaner energy transitions (1)
- climat change (1)
- climate adaptation (1)
- climate benefits (1)
- climate change adaptation (1)
- climate gradient (1)
- climate modeling (1)
- climate time series analysis (1)
- climate transition (1)
- computational geosciences (1)
- conflict (1)
- cosmic-rays (1)
- cosmogenic nuclides (1)
- covalency (1)
- creep and unmapped faults; (1)
- crustal structure (1)
- damage modelling (1)
- data-mining (1)
- deep biosphere (1)
- deep long-period earthquakes (1)
- deformation (1)
- degradation (1)
- density-driven flow (1)
- deposition (1)
- derived flood risk analysis (1)
- deterministic (1)
- diamond anvil cell (1)
- dichtegetriebene Strömung (1)
- digital elevation models (1)
- directivity analysis (1)
- disaster risk (1)
- discharge (1)
- discharge time series (1)
- disturbance (1)
- downscaling (1)
- drainage-divide migration (1)
- earth system sciences (1)
- earthquake damage (1)
- earthquake hydrology (1)
- earthquakes (1)
- ecosystem service cascade (1)
- electrical resistivity tomography (1)
- electrochemistry (1)
- emergency management (1)
- empirical (1)
- enhanced geothermal system (1)
- enhanced geothermal systems (EGS) (1)
- environmental (1)
- environmentalism (1)
- equifinality (1)
- erosion (1)
- eruption (1)
- estimation (1)
- event characteristics (1)
- expedition (1)
- exposure (1)
- extensional tectonics (1)
- extreme seasons (1)
- fault network (1)
- fault reactivation (1)
- faults (1)
- fields (1)
- flat-slab (1)
- flood frequency (1)
- flood hazard mapping; (1)
- flood modelling; (1)
- flood risk analysis (1)
- floodplain wetlands (1)
- fluctuation (1)
- fluvial erosion (1)
- flux (1)
- forecasting (1)
- forecasting (1922, 4315, 7924, 7964) (1)
- formation damage (1)
- fracture growth (1)
- fracture transmissivity (1)
- frequency-domain electromagnetics (1)
- gas flux (1)
- gas hydrate (1)
- gas hydrates (1)
- generative model (1)
- geochronology (1)
- geogenic degassing (1)
- geoinformatics (1)
- geomodeling (1)
- geomorphic connectivity (1)
- geostatistics (1)
- geothermal energy (1)
- geyser (1)
- grain size (1)
- greenhouse gas (1)
- ground motion fields (1)
- ground motion modeling (1)
- groundwater recharge (1)
- habitat (1)
- heat flow (1)
- heavy tail behavior (1)
- hexadecane (1)
- high spatial resolution (1)
- high-pressure (1)
- historical floods (1)
- historical geography (1)
- historische Hochwasser (1)
- hominid (1)
- horizontaler Fluss (1)
- host rock scale (1)
- human evolution (1)
- hybrid (1)
- hydrate formation (1)
- hydraulische Risserzeugung (1)
- hydraulisches Aufbrechen (1)
- hydrodynamic interactions (1)
- hydrogen bond (1)
- hydrological modeling (1)
- hydrological uplift and subsidence (1)
- hydrothermal system; (1)
- ice complex (1)
- imaging (1)
- impact (1)
- incommensurate modulation (1)
- induced seismicity (1)
- injection (1)
- injection scheme (1)
- intercultural understanding (1)
- interdisciplinarity (1)
- interseismic strain rate (1)
- inundation simulation (1)
- inverse analysis (1)
- inverse modeling (1)
- ionic liquids (1)
- iron (1)
- isotopic fractionation (1)
- knickpoints (1)
- lake system responses (1)
- larch forest (1)
- laterally constrained inversion (1)
- lethal violence (1)
- level (1)
- lithosphere (1)
- low-frequency (1)
- machine learning (0555) (1)
- magma degassing (1)
- magma plumbing system (1)
- marine sediments (1)
- maschinelles Lernen (1)
- mathematical geology (1)
- mesoscale heterogeneous catchment (1)
- metal-containing ionic liquids; (1)
- methane (1)
- methane hydrate (1)
- microbial community (1)
- microbial diversity (1)
- microbial mats (1)
- microplate (1)
- microstructural deformation mechanisms (1)
- minimum gradient support regularization (1)
- mining sciences (1)
- model (1)
- modeling; (1)
- molecular biomarkers (1)
- molekulare Biomarker (1)
- morphometric indices (1)
- mountain growth (1)
- multiperspective reception of AvH (1)
- multisensor system (1)
- multivariate statistics (1)
- muons (1)
- near-surface monitoring (1)
- neutrons (1)
- nitrogen (1)
- nitrogen legacies (1)
- normal faults (1)
- nuclear waste disposal (1)
- numerical simulation (1)
- numerical simulations (1)
- offene Daten (1)
- online (1)
- open data (1)
- organic carbon stabilization (1)
- organic matter (1)
- organisches Material (1)
- orographic rainfall (1)
- paleoceanography (1)
- paleoclimatology (1)
- paleoseismology (1)
- parameter (1)
- peat characterization (1)
- pegmatite (1)
- periglacial process (1)
- permafost (1)
- permafrost sediments (1)
- permutation entropy (1)
- petrothermales System (EGS) (1)
- phase transition (1)
- phytoplankton (1)
- planktonic foraminifera (1)
- planning indicators (1)
- plasma density (1)
- plastic degradation (1)
- plastisphere (1)
- pollen (1)
- population pressure (1)
- poroelasticity (1)
- power-law rheology (1)
- precaution (1)
- precipitation (1)
- propagating uplift (1)
- provenance (1)
- quartäre Geochronologie (1)
- radial flow (1)
- radiale Strömung (1)
- radiocarbon dating (1)
- reactive transport simulation (1)
- reaktive Transportsimulation (1)
- redox processes (1)
- reduction (1)
- regional flood frequency analysis (1)
- relaxation (1)
- remote sensing (1)
- renya Rift (1)
- reservoirs (1)
- rheology (1)
- rift (1)
- rift transfer zone (1)
- rifts (1)
- risk (1)
- risk management (1)
- river dynamics (1)
- river incision (1)
- rotational sensor (1)
- russische Arktis (1)
- salt rock (1)
- sediment (1)
- sediment dynamics (1)
- sediment routing system (1)
- sedimentary aquifer (1)
- seismic monitoring (1)
- seismic tomography (1)
- seismicity (1)
- seismicity model (1)
- sensitivity analysis (1)
- slope aspect (1)
- soil formation (1)
- soil-moisture (1)
- sorption (1)
- source rock evaluation (1)
- spin (1)
- stabile Isotope (1)
- stable carbon isotopes (1)
- stable isotopes (1)
- statistics (1)
- steric factors (1)
- stress measurement (1)
- strike-slip (1)
- subcatchment response (1)
- subduction zone (1)
- submarine gas vents (1)
- subsurface life (1)
- sulfate reduction (1)
- surface type-function-concept (1)
- surface velocity (1)
- surface wave (1)
- surface wave, inversion, near surface (1)
- tectonic geodesy (1)
- temperature sensor (1)
- tensor inversion (1)
- terrestrial (1)
- tetrahalidometallates (1)
- thaw (1)
- thermochronological modelling (AFT, AHe and ZHe) (1)
- thermochronology (1)
- thermodynamic and kinetic properties (1)
- thermodynamische und kinetische Eigenschaften (1)
- thermokarst (1)
- thermokarst lakes (1)
- thermophiles; (1)
- time series analysis (1)
- tourmaline (1)
- transition (1)
- travel time distribution (1)
- triaxial deformation experiments (1)
- triaxiale Deformationsexperimente (1)
- tropical climate (1)
- turbidity modelling (1)
- twin domains (1)
- ultra-relativistic electrons (1)
- upper tail (1)
- uranium (1)
- urban green (1)
- urban planning practice (1)
- variability (1)
- varved lake sediments (1)
- verbesserte geothermische Systeme (1)
- vertikaler Fluss (1)
- volcano-seismology (1)
- volcanoseismology (1)
- vulnerability (1)
- warfare (1)
- warvierte Seesedimente (1)
- water balance (1)
- water quality (1)
- water quality modeling (1)
- water rock interactions (1)
- water table (1)
- wetland (1)
- wetland hydrology (1)
- wind energy-biodiversity (1)
- wind energy-biodiversity conflict (1)
- wind gusts (1)
- winderosion (1)
- yedoma (1)
- yedoma ice complex (1)
- Überflutungssimulation (1)
Institute
- Institut für Geowissenschaften (157) (remove)
Ground-penetrating radar (GPR) is a method that can provide detailed information about the near subsurface in sedimentary and carbonate environments.
The classical interpretation of GPR data (e.g., based on manual feature selection) often is labor-intensive and limited by the experience of the intercally used for seismic interpretation, can provide faster, more repeatable, and less biased interpretations. We have recorded a 3D GPD data set collected across a paleokarst breccia pipe in the Billefjorden area on Spitsbergen, Svalbard. After performing advanced processing, we compare the results of a classical GPR interpretation to the results of an attribute-based classification.
Our attribute classification incorporates a selection of dip and textural attributes as the input for a k-means clustering approach. Similar to the results of the classical interpretation, the resulting classes differentiate between undisturbed strata and breccias or fault zones.
The classes also reveal details inside the breccia pipe that are not discerned in the classical fer that the intrapipe GPR facies result from subtle differences, such as breccia lithology, clast size, or pore-space filling.
The region of West Bohemia and Upper Palatinate belongs to the West Bohemian Massif. The study area is situated at the junction of three different Variscan tectonic units and hosts the ENE-WSW trending Ohre Rift as well as many different fault systems. The entire region is characterized by ongoing magmatic processes in the intra-continental lithospheric mantle expressed by a series of phenomena, including e.g. the occurrence of repeated earthquake swarms and massive degassing of mantle derived CO2 in form of mineral springs and mofettes. Ongoing active tectonics is mainly manifested by Cenozoic volcanism represented by different Quaternary volcanic structures. All these phenomena make the Ohre Rift a unique target area for European intra-continental geo-scientific research. With magnetotelluric (MT) measurements we image the subsurface distribution of the electrical resistivity and map possible fluid pathways. Two-dimensional (2D) inversion results by Munoz et al. (2018) reveal a conductive channel in the vicinity of the earthquake swarm region that extends from the lower crust to the surface forming a pathway for fluids into the region of the mofettes. A second conductive channel is present in the south of their model; however, their 2D inversions allow ambiguous interpretations of this feature. Therefore, we conducted a large 3D MT field experiment extending the study area towards the south. The 3D inversion result matches well with the known geology imaging different fluid/magma reservoirs at crust-mantle depth and mapping possible fluid pathways from the reservoirs to the surface feeding known mofettes and spas. A comparison of 3D and 2D inversion results suggests that the 2D inversion results are considerably characterized by 3D and off-profile structures. In this context, the new results advocate for the swarm earthquakes being located in the resistive host rock surrounding the conductive channels; a finding in line with observations e.g. at the San Andreas Fault, California.
40Ar/39Ar dating of a hydrothermal pegmatitic buddingtonite–muscovite assemblage from Volyn, Ukraine
(2022)
We determined Ar-40/Ar-39 ages of buddingtonite, occurring together with muscovite, with the laser-ablation method. This is the first attempt to date the NH4-feldspar buddingtonite, which is typical for sedimentary-diagenetic environments of sediments, rich in organic matter, or in hydrothermal environments, associated with volcanic geyser systems. The sample is a hydrothermal breccia, coming from the Paleoproterozoic pegmatite field of the Korosten Plutonic Complex, Volyn, Ukraine. A detailed characterization by optical methods, electron microprobe analyses, backscattered electron imaging, and IR analyses showed that the buddingtonite consists of euhedral-appearing platy crystals of tens of micrometers wide, 100 or more micrometers in length, which consist of fine-grained fibers of <= 1 mu m thickness. The crystals are sector and growth zoned in terms of K-NH4-H3O content. The content of K allows for an age determination with the Ar-40/Ar-39 method, as well as in the accompanying muscovite, intimately intergrown with the buddingtonite. The determinations on muscovite yielded an age of 1491 +/- 9 Ma, interpreted as the hydrothermal event forming the breccia. However, buddingtonite apparent ages yielded a range of 563 +/- 14 Ma down to 383 +/- 12 Ma, which are interpreted as reset ages due to Ar loss of the fibrous buddingtonite crystals during later heating. We conclude that buddingtonite is suited for Ar-40/Ar-39 age determinations as a supplementary method, together with other methods and minerals; however, it requires a detailed mineralogical characterization, and the ages will likely represent minimum ages.
The quantification of spatial propagation of extreme precipitation events is vital in water resources planning and disaster mitigation. However, quantifying these extreme events has always been challenging as many traditional methods are insufficient to capture the nonlinear interrelationships between extreme event time series. Therefore, it is crucial to develop suitable methods for analyzing the dynamics of extreme events over a river basin with a diverse climate and complicated topography. Over the last decade, complex network analysis emerged as a powerful tool to study the intricate spatiotemporal relationship between many variables in a compact way. In this study, we employ two nonlinear concepts of event synchronization and edit distance to investigate the extreme precipitation pattern in the Ganga river basin. We use the network degree to understand the spatial synchronization pattern of extreme rainfall and identify essential sites in the river basin with respect to potential prediction skills. The study also attempts to quantify the influence of precipitation seasonality and topography on extreme events. The findings of the study reveal that (1) the network degree is decreased in the southwest to northwest direction, (2) the timing of 50th percentile precipitation within a year influences the spatial distribution of degree, (3) the timing is inversely related to elevation, and (4) the lower elevation greatly influences connectivity of the sites. The study highlights that edit distance could be a promising alternative to analyze event-like data by incorporating event time and amplitude and constructing complex networks of climate extremes.
A comprehensive workflow to analyze ensembles of globally inverted 2D electrical resistivity models
(2022)
Electrical resistivity tomography (ERT) aims at imaging the subsurface resistivity distribution and provides valuable information for different geological, engineering, and hydrological applications. To obtain a subsurface resistivity model from measured apparent resistivities, stochastic or deterministic inversion procedures may be employed. Typically, the inversion of ERT data results in non-unique solutions; i.e., an ensemble of different models explains the measured data equally well. In this study, we perform inference analysis of model ensembles generated using a well-established global inversion approach to assess uncertainties related to the nonuniqueness of the inverse problem. Our interpretation strategy starts by establishing model selection criteria based on different statistical descriptors calculated from the data residuals. Then, we perform cluster analysis considering the inverted resistivity models and the corresponding data residuals. Finally, we evaluate model uncertainties and residual distributions for each cluster. To illustrate the potential of our approach, we use a particle swarm optimization (PSO) algorithm to obtain an ensemble of 2D layer-based resistivity models from a synthetic data example and a field data set collected in Loon-Plage, France. Our strategy performs well for both synthetic and field data and allows us to extract different plausible model scenarios with their associated uncertainties and data residual distributions. Although we demonstrate our workflow using 2D ERT data and a PSObased inversion approach, the proposed strategy is general and can be adapted to analyze model ensembles generated from other kinds of geophysical data and using different global inversion approaches.
Deep hydrothermal Mo, W, and base metal mineralization at the Sweet Home mine (Detroit City portal) formed in response to magmatic activity during the Oligocene. Microthermometric data of fluid inclusions trapped in greisen quartz and fluorite suggest that the early-stage mineralization at the Sweet Home mine precipitated from low- to medium-salinity (1.5-11.5 wt% equiv. NaCl), CO2-bearing fluids at temperatures between 360 and 415 degrees C and at depths of at least 3.5 km. Stable isotope and noble gas isotope data indicate that greisen formation and base metal mineralization at the Sweet Home mine was related to fluids of different origins. Early magmatic fluids were the principal source for mantle-derived volatiles (CO2, H2S/SO2, noble gases), which subsequently mixed with significant amounts of heated meteoric water. Mixing of magmatic fluids with meteoric water is constrained by delta H-2(w)-delta O-18(w) relationships of fluid inclusions. The deep hydrothermal mineralization at the Sweet Home mine shows features similar to deep hydrothermal vein mineralization at Climax-type Mo deposits or on their periphery. This suggests that fluid migration and the deposition of ore and gangue minerals in the Sweet Home mine was triggered by a deep-seated magmatic intrusion. The findings of this study are in good agreement with the results of previous fluid inclusion studies of the mineralization of the Sweet Home mine and from Climax-type Mo porphyry deposits in the Colorado Mineral Belt.
A large landslide (frozen debris avalanche) occurred at Assapaat on the south coast of the Nuussuaq Peninsula in Central West Greenland on June 13, 2021, at 04:04 local time. We present a compilation of available data from field observations, photos, remote sensing, and seismic monitoring to describe the event. Analysis of these data in combination with an analysis of pre- and post-failure digital elevation models results in the first description of this type of landslide. The frozen debris avalanche initiated as a 6.9 * 10(6) m(3) failure of permafrozen talus slope and underlying colluvium and till at 600-880 m elevation. It entrained a large volume of permafrozen colluvium along its 2.4 km path in two subsequent entrainment phases accumulating a total volume between 18.3 * 10(6) and 25.9 * 10(6) m(3). About 3.9 * 10(6) m(3) is estimated to have entered the Vaigat strait; however, no tsunami was reported, or is evident in the field. This is probably because the second stage of entrainment along with a flattening of slope angle reduced the mobility of the frozen debris avalanche. We hypothesise that the initial talus slope failure is dynamically conditioned by warming of the ice matrix that binds the permafrozen talus slope. When the slope ice temperature rises to a critical level, its shear resistance is reduced, resulting in an unstable talus slope prone to failure. Likewise, we attribute the large-scale entrainment to increasing slope temperature and take the frozen debris avalanche as a strong sign that the permafrost in this region is increasingly at a critical state. Global warming is enhanced in the Arctic and frequent landslide events in the past decade in Western Greenland let us hypothesise that continued warming will lead to an increase in the frequency and magnitude of these types of landslides. Essential data for critical arctic slopes such as precipitation, snowmelt, and ground and surface temperature are still missing to further test this hypothesis. It is thus strongly required that research funds are made available to better predict the change of landslide threat in the Arctic.
Van Allen Probes measurements revealed the presence of the most unusual structures in the ultra-relativistic radiation belts. Detailed modeling, analysis of pitch angle distributions, analysis of the difference between relativistic and ultra-realistic electron evolution, along with theoretical studies of the scattering and wave growth, all indicate that electromagnetic ion cyclotron (EMIC) waves can produce a very efficient loss of the ultra-relativistic electrons in the heart of the radiation belts. Moreover, a detailed analysis of the profiles of phase space densities provides direct evidence for localized loss by EMIC waves. The evolution of multi-MeV fluxes shows dramatic and very sudden enhancements of electrons for selected storms. Analysis of phase space density profiles reveals that growing peaks at different values of the first invariant are formed at approximately the same radial distance from the Earth and show the sequential formation of the peaks from lower to higher energies, indicating that local energy diffusion is the dominant source of the acceleration from MeV to multi-MeV energies. Further simultaneous analysis of the background density and ultra-relativistic electron fluxes shows that the acceleration to multi-MeV energies only occurs when plasma density is significantly depleted outside of the plasmasphere, which is consistent with the modeling of acceleration due to chorus waves.
A novel approach for estimating precipitation patterns is developed here and applied to generate a new hydrologically corrected daily precipitation dataset, called RAIN4PE (Rain for Peru and Ecuador), at 0.1 degrees spatial resolution for the period 1981-2015 covering Peru and Ecuador. It is based on the application of 1) the random forest method to merge multisource precipitation estimates (gauge, satellite, and reanalysis) with terrain elevation, and 2) observed and modeled streamflow data to first detect biases and second further adjust gridded precipitation by inversely applying the simulated results of the ecohydrological model SWAT (Soil and Water Assessment Tool). Hydrological results using RAIN4PE as input for the Peruvian and Ecuadorian catchments were compared against the ones when feeding other uncorrected (CHIRP and ERA5) and gauge-corrected (CHIRPS, MSWEP, and PISCO) precipitation datasets into the model. For that, SWAT was calibrated and validated at 72 river sections for each dataset using a range of performance metrics, including hydrograph goodness of fit and flow duration curve signatures. Results showed that gauge-corrected precipitation datasets outperformed uncorrected ones for streamflow simulation. However, CHIRPS, MSWEP, and PISCO showed limitations for streamflow simulation in several catchments draining into the Pacific Ocean and the Amazon River. RAIN4PE provided the best overall performance for streamflow simulation, including flow variability (low, high, and peak flows) and water budget closure. The overall good performance of RAIN4PE as input for hydrological modeling provides a valuable criterion of its applicability for robust countrywide hydrometeorological applications, including hydroclimatic extremes such as droughts and floods. Significance StatementWe developed a novel precipitation dataset RAIN4PE for Peru and Ecuador by merging multisource precipitation data (satellite, reanalysis, and ground-based precipitation) with terrain elevation using the random forest method. Furthermore, RAIN4PE was hydrologically corrected using streamflow data in watersheds with precipitation underestimation through reverse hydrology. The results of a comprehensive hydrological evaluation showed that RAIN4PE outperformed state-of-the-art precipitation datasets such as CHIRP, ERA5, CHIRPS, MSWEP, and PISCO in terms of daily and monthly streamflow simulations, including extremely low and high flows in almost all Peruvian and Ecuadorian catchments. This underlines the suitability of RAIN4PE for hydrometeorological applications in this region. Furthermore, our approach for the generation of RAIN4PE can be used in other data-scarce regions.
A review of source models to further the understanding of the seismicity of the Groningen field
(2022)
The occurrence of felt earthquakes due to gas production in Groningen has initiated numerous studies and model attempts to understand and quantify induced seismicity in this region. The whole bandwidth of available models spans the range from fully deterministic models to purely empirical and stochastic models. In this article, we summarise the most important model approaches, describing their main achievements and limitations. In addition, we discuss remaining open questions and potential future directions of development.
Groundwater recharge (GWR) is one of the most challenging water fluxes to estimate, as it relies on observed data that are often limited in many developing countries.
This study developed an innovative water budget method using satellite products for estimating the spatially distributed GWR at monthly and annual scales in tropical wet sedimentary regions despite cloudy conditions.
The distinctive features proposed in this study include the capacity to address 1) evapotranspiration estimations in tropical wet regions frequently overlaid by substantial cloud cover; and 2) seasonal root-zone water storage estimations in sedimentary regions prone to monthly variations.
The method also utilises satellite-based information of the precipitation and surface runoff. The GWR was estimated and validated for the hydrologically contrasting years 2016 and 2017 over a tropical wet sedimentary region located in North-eastern Brazil, which has substantial potential for groundwater abstraction.
This study showed that applying a cloud-cleaning procedure based on monthly compositions of biophysical data enables the production of a reasonable proxy for evapotranspiration able to estimate groundwater by the water budget method.
The resulting GWR rates were 219 (2016) and 302 (2017) mm yr(-1), showing good correlations (CC = 0.68 to 0.83) and slight underestimations (PBIAS =-13 to-9%) when compared with the referenced estimates obtained by the water table fluctuation method for 23 monitoring wells. Sensitivity analysis shows that water storage changes account for +19% to-22% of our monthly evaluation.
The satellite-based approach consistently demonstrated that the consideration of cloud-cleaned evapotranspiration and root-zone soil water storage changes are essential for a proper estimation of spatially distributed GWR in tropical wet sedimentary regions because of their weather seasonality and cloudy conditions.
Different lake systems might reflect different climate elements of climate changes, while the responses of lake systems are also divers, and are not completely understood so far. Therefore, a comparison of lakes in different climate zones, during the high-amplitude and abrupt climate fluctuations of the Last Glacial to Holocene transition provides an exceptional opportunity to investigate distinct natural lake system responses to different abrupt climate changes. The aim of this doctoral thesis was to reconstruct climatic and environmental fluctuations down to (sub-) annual resolution from two different lake systems during the Last Glacial-Interglacial transition (~17 and 11 ka). Lake Gościąż, situated in the temperate central Poland, developed in the Allerød after recession of the Last Glacial ice sheets. The Dead Sea is located in the Levant (eastern Mediterranean) within a steep gradient from sub-humid to hyper-arid climate, and formed in the mid-Miocene. Despite their differences in sedimentation processes, both lakes form annual laminations (varves), which are crucial for studies of abrupt climate fluctuations. This doctoral thesis was carried out within the DFG project PALEX-II (Paleohydrology and Extreme Floods from the Dead Sea ICDP Core) that investigates extreme hydro-meteorological events in the ICDP core in relation to climate changes, and ICLEA (Virtual Institute of Integrated Climate and Landscape Evolution Analyses) that intends to better the understanding of climate dynamics and landscape evolutions in north-central Europe since the Last Glacial. Further, it contributes to the Helmholtz Climate Initiative REKLIM (Regional Climate Change and Humans) Research Theme 3 “Extreme events across temporal and spatial scales” that investigates extreme events using climate data, paleo-records and model-based simulations. The three main aims were to (1) establish robust chronologies of the lakes, (2) investigate how major and abrupt climate changes affect the lake systems, and (3) to compare the responses of the two varved lakes to these hemispheric-scale climate changes.
Robust chronologies are a prerequisite for high-resolved climate and environmental reconstructions, as well as for archive comparisons. Thus, addressing the first aim, the novel chronology of Lake Gościąż was established by microscopic varve counting and Bayesian age-depth modelling in Bacon for a non-varved section, and was corroborated by independent age constrains from 137Cs activity concentration measurements, AMS radiocarbon dating and pollen analysis. The varve chronology reaches from the late Allerød until AD 2015, revealing more Holocene varves than a previous study of Lake Gościąż suggested. Varve formation throughout the complete Younger Dryas (YD) even allowed the identification of annually- to decadal-resolved leads and lags in proxy responses at the YD transitions.
The lateglacial chronology of the Dead Sea (DS) was thus far mainly based on radiocarbon and U/Th-dating. In the unique ICDP core from the deep lake centre, continuous search for cryptotephra has been carried out in lateglacial sediments between two prominent gypsum deposits – the Upper and Additional Gypsum Units (UGU and AGU, respectively). Two cryptotephras were identified with glass analyses that correlate with tephra deposits from the Süphan and Nemrut volcanoes indicating that the AGU is ~1000 years younger than previously assumed, shifting it into the YD, and the underlying varved interval into the Bølling/Allerød, contradicting previous assumptions.
Using microfacies analyses, stable isotopes and temperature reconstructions, the second aim was achieved at Lake Gościąż. The YD lake system was dynamic, characterized by higher aquatic bioproductivity, more re-suspended material and less anoxia than during the Allerød and Early Holocene, mainly influenced by stronger water circulation and catchment erosion due to stronger westerly winds and less lake sheltering. Cooling at the YD onset was ~100 years longer than the final warming, while environmental proxies lagged the onset of cooling by ~90 years, but occurred contemporaneously during the termination of the YD. Chironomid-based temperature reconstructions support recent studies indicating mild YD summer temperatures. Such a comparison of annually-resolved proxy responses to both abrupt YD transitions is rare, because most European lake archives do not preserve varves during the YD.
To accomplish the second aim at the DS, microfacies analyses were performed between the UGU (~17 ka) and Holocene onset (~11 ka) in shallow- (Masada) and deep-water (ICDP core) environments. This time interval is marked by a huge but fluctuating lake level drop and therefore the complete transition into the Holocene is only recorded in the deep-basin ICDP core. In this thesis, this transition was investigated for the first time continuously and in detail. The final two pronounced lake level drops recorded by deposition of the UGU and AGU, were interrupted by one millennium of relative depositional stability and a positive water budget as recorded by aragonite varve deposition interrupted by only a few event layers. Further, intercalation of aragonite varves between the gypsum beds of the UGU and AGU shows that these generally dry intervals were also marked by decadal- to centennial-long rises in lake level. While continuous aragonite varves indicate decadal-long stable phases, the occurrence of thicker and more frequent event layers suggests general more instability during the gypsum units. These results suggest a pattern of complex and variable hydroclimate at different time scales during the Lateglacial at the DS.
The third aim was accomplished based on the individual studies above that jointly provide an integrated picture of different lake responses to different climate elements of hemispheric-scale abrupt climate changes during the Last Glacial-Interglacial transition. In general, climatically-driven facies changes are more dramatic in the DS than at Lake Gościąż. Further, Lake Gościąż is characterized by continuous varve formation nearly throughout the complete profile, whereas the DS record is widely characterized by extreme event layers, hampering the establishment of a continuous varve chronology. The lateglacial sedimentation in Lake Gościąż is mainly influenced by westerly winds and minor by changes in catchment vegetation, whereas the DS is primarily influenced by changes in winter precipitation, which are caused by temperature variations in the Mediterranean. Interestingly, sedimentation in both archives is more stable during the Bølling/Allerød and more dynamic during the YD, even when sedimentation processes are different.
In summary, this doctoral thesis presents seasonally-resolved records from two lake archives during the Lateglacial (ca 17-11 ka) to investigate the impact of abrupt climate changes in different lake systems. New age constrains from the identification of volcanic glass shards in the lateglacial sediments of the DS allowed the first lithology-based interpretation of the YD in the DS record and its comparison to Lake Gościąż. This highlights the importance of the construction of a robust chronology, and provides a first step for synchronization of the DS with other eastern Mediterranean archives. Further, climate reconstructions from the lake sediments showed variability on different time scales in the different archives, i.e. decadal- to millennial fluctuations in the lateglacial DS, and even annual variations and sub-decadal leads and lags in proxy responses during the rapid YD transitions in Lake Gościąż. This showed the importance of a comparison of different lake archives to better understand the regional and local impacts of hemispheric-scale climate variability. An unprecedented example is demonstrated here of how different lake systems show different lake responses and also react to different climate elements of abrupt climate changes. This further highlights the importance of the understanding of the respective lake system for climate reconstructions.
Flood risk management in Germany follows an integrative approach in which both private households and businesses can make an important contribution to reducing flood damage by implementing property-level adaptation measures. While the flood adaptation behavior of private households has already been widely researched, comparatively less attention has been paid to the adaptation strategies of businesses. However, their ability to cope with flood risk plays an important role in the social and economic development of a flood-prone region. Therefore, using quantitative survey data, this study aims to identify different strategies and adaptation drivers of 557 businesses damaged by a riverine flood in 2013 and 104 businesses damaged by pluvial or flash floods between 2014 and 2017. Our results indicate that a low perceived self-efficacy may be an important factor that can reduce the motivation of businesses to adapt to flood risk. Furthermore, property-owners tended to act more proactively than tenants. In addition, high experience with previous flood events and low perceived response costs could strengthen proactive adaptation behavior. These findings should be considered in business-tailored risk communication.
Reducing greenhouse gas emissions in food systems is becoming more challenging as food is increasingly consumed away from producer regions, highlighting the need to consider emissions embodied in trade in agricultural emissions accounting.
To address this, our study explores recent trends in trade-adjusted agricultural emissions of food items at the global, regional, and national levels.
We find that emissions are largely dependent on a country’s consumption patterns and their agricultural emission intensities relative to their trading partners’.
The absolute differences between the production-based and trade-adjusted emissions accounting approaches are especially apparent for major agricultural exporters and importers and where large shares of emission-intensive items such as ruminant meat, milk products and rice are involved.
In relative terms, some low-income and emerging and developing economies with consumption of high emission intensity food products show large differences between approaches.
Similar trends are also found under various specifications that account for trade and re-exports differently.
These findings could serve as an important element towards constructing national emissions reduction targets that consider trading partners, leading to more effective emissions reductions overall.
Marine sedimentary archives are routinely used to reconstruct past environmental changes. In many cases, bioturbation and sedimentary mixing affect the proxy time-series and the age-depth relationship. While idealized models of bioturbation exist, they usually assume homogeneous mixing, thus that a single sample is representative for the sediment layer it is sampled from.
However, it is largely unknown to which extent this assumption holds for sediments used for paleoclimate reconstructions.
To shed light on
1) the age-depth relationship and its full uncertainty,
2) the magnitude of mixing processes affecting the downcore proxy variations, and
3) the representativity of the discrete sample for the sediment layer, we designed and performed a case study on South China Sea sediment material which was collected using a box corer and which covers the last glacial cycle.
Using the radiocarbon content of foraminiferal tests as a tracer of time, we characterize the spatial age-heterogeneity of sediments in a three-dimensional setup. In total, 118 radiocarbon measurements were performed on defined small- and large-volume bulk samples ( similar to 200 specimens each) to investigate the horizontal heterogeneity of the sediment. Additionally, replicated measurements on small numbers of specimens (10 x 5 specimens) were performed to assess the heterogeneity within a sample volume. Visual assessment of X-ray images and a quantitative assessment of the mixing strength show typical mixing from bioturbation corresponding to around 10 cm mixing depth.
Notably, our 3D radiocarbon distribution reveals that the horizontal heterogeneity (up to 1,250 years), contributing to the age uncertainty, is several times larger than the typically assumed radiocarbon based age-model error (single errors up to 250 years). Furthermore, the assumption of a perfectly bioturbated layer with no mixing underneath is not met.
Our analysis further demonstrates that the age-heterogeneity might be a function of sample size; smaller samples might contain single features from the incomplete mixing and are thus less representative than larger samples.
We provide suggestions for future studies, optimal sampling strategies for quantitative paleoclimate reconstructions and realistic uncertainty in age models, as well as discuss possible implications for the interpretation of paleoclimate records.
Alexander von Humboldt
(2022)
This book aims to view and to understand Alexander von Humboldt from different perspectives and in varying disciplinary contexts. His contributions addressed numerous topics in the earth but also life sciences—spanning from geo-botany, climatology, paleontology, oceanography, mineralogy, resources, and hydrogeology to links between the environmental impact of humans, erosion, and climate change. From the very beginning, he paved the way for a modern, integrated earth system science approach to decipher, characterize, and model the different forcing factors and their feedback mechanisms. It becomes obvious that Humboldt’s holistic approach is far beyond simple description and empiric data collection. As documented and analyzed in the different texts of this volume, he combines observation and analysis with emotions and subjective perceptions in a very affectionate way. However, this publication does not intend to add another encyclopedic text compilation but to observe and critically analyze this unique personality´s relevance in a modern context, particularly in discussing environmental and social key issues in the twenty-first century.
Measuring the variability of incoming neutrons locally would be usefull for the cosmic-ray neutron sensing (CRNS) method. As the measurement of high energy neutrons is not so easy, alternative particles can be considered for such purpose. Among them, muons are particles created from the same cascade of primary cosmic-ray fluxes that generate neutrons at the ground. In addition, they can be easily detected by small and relatively inexpensive detectors. For these reasons they could provide a suitable local alternative to incoming corrections based on remote neutron monitor data. The reported measurements demonstrated that muon detection system can detect incoming cosmic-ray variations locally. Furthermore the precision of this measurement technique is considered adequate for many CRNS applications.
Die Hochwasserkatastrophe im Juli 2021 in Westdeutschland erfordert eine kritische Diskussion über die Abschätzung der Hochwassergefährdung, Aktualisierung von Hochwassergefahrenkarten und Kommunikation von extremen Hochwasserszenarien. In der vorliegenden Arbeit wurde die Extremwertstatistik für die jährlichen maximalen Spitzenabflüsse am Pegel Altenahr im Ahrtal mit und ohne Berücksichtigung historischer Hochwasser berechnet und verglichen. Die Schätzung der Wiederkehrperiode für das aktuelle Hochwasser mittels Generalisierter Extremwertverteilung (GEV) unter Berücksichtigung historischer Hochwasser schwankt zwischen etwa 2.600 und über 58.700 Jahren (90%-Konfidenzintervall) mit einem Median bei etwa 8.600 Jahren, wogegen die Schätzung, die nur auf der systematisch gemessenen Abflusszeitreihe von 74 Jahren basiert, theoretisch eine Wiederkehrperiode von über 100 Millionen Jahren ergeben würde. Die Berücksichtigung der historischen Hochwasser führt zu einer dramatischen Änderung der Hochwasserquan-
tile, die für eine Gefahrenkartierung zugrunde gelegt werden. Die Anpassung der GEV an die Zeitreihe mit historischen Hochwassern zeigt dennoch, dass das GEV-Modell möglicherweise die Grundgesamtheit der Hochwasser im Ahrtal nicht adäquat abbilden kann. Es könnte sich im vorliegenden Fall um eine gemischte Stichprobe handeln, in der die extremen Hochwasser im Vergleich zu kleineren Ereignissen durch besondere Prozesse hervorgerufen werden. Somit könnten die Wahrscheinlichkeiten von extremen Hochwassern deutlich größer sein, als aus dem GEV-Modell hervorgeht. Hier sollte in Zukunft die Anwendung einer prozessbasierten Mischverteilung
untersucht werden. Der Vergleich von amtlichen Gefahrenkarten zu Extremhochwassern (HQextrem) im Ahrtal mit den Überflutungsflächen vom Juli 2021
zeigt eine deutliche Diskrepanz in den betroffenen Gebieten und die Notwendigkeit, die Grundlagen zur Erstellung der Extremszenarien zu überdenken. Die hydrodynamisch-numerischen Simulationen von 1.000-jährlichen Hochwassern (HQ1000) unter Berücksichtigung historischer Ereignisse und des größten historischen Hochwassers 1804 können die Gefährdung des Juli-Hochwassers 2021 deutlich besser widerspiegeln, wenngleich auch diese beiden Szenarien die Überflutungsflächen unterschätzen. Besondere Effekte wie die Verklausung von Brücken und die geomorphologischen Änderungen im Flussschlauch führten zu noch größeren Überflutungs- flächen im Juli 2021, als die Simulationsergebnisse zeigten. Basierend auf dieser Analyse wird eine einheitliche Festlegung von HQextrem bei Hochwassergefahrenkartierungen in Deutschland vorgeschlagen, die sich an höheren Hochwasserquantilen im Bereich von HQ1000 orientiert. Zusätzlich sollen simulationsbasierte Rekonstruktionen von den größten verlässlich dokumentierten historischen Hochwassern und/oder synthetische Worst-Case-Szenarien in den Hochwassergefahrenkarten gesondert dargestellt werden. Damit wird ein wichtiger Beitrag geleistet, um die potenziell betroffene Bevölkerung und das Katastrophenmanagement vor Überraschungen durch sehr seltene und extreme Hochwasser in Zukunft besser zu schützen.
The main Marmara fault (MMF) extends for 150 km through the Sea of Marmara and forms the only portion of the North Anatolian fault zone that has not ruptured in a large event (Mw >7) for the last 250 yr. Accordingly, this portion is potentially a major source contributing to the seismic hazard of the Istanbul region. On 26 September 2019, a sequence of moderate-sized events started along the MMF only 20 km south of Istanbul and were widely felt by the population. The largest three events, 26 September Mw 5.8 (10:59 UTC), 26 September 2019 Mw 4.1 (11:26 UTC), and 20 January 2020 Mw 4.7 were recorded by numerous strong-motion seismic stations and the resulting ground motions were compared to the predicted means resulting from a set of the most recent ground-motion prediction equations (GMPEs). The estimated residuals were used to investigate the spatial variation of ground motion across the Marmara region. Our results show a strong azimuthal trend in ground-motion residuals, which might indicate systematically repeating directivity effects toward the eastern Marmara region.
Extracting information about past tectonic or climatic environmental changes from sedimentary records is a key objective of provenance research. Interpreting the imprint of such changes remains challenging as signals might be altered in the sediment-routing system.
We investigate the sedimentary provenance of the Oligocene/Miocene Upper Austrian Northern Alpine Foreland Basin and its response to the tectonically driven exhumation of the Tauern Window metamorphic dome (28 +/- 1 Ma) in the Eastern European Alps by using the unprecedented combination of Nd isotopic composition of bulk-rock clay-sized samples and partly previously published multi-proxy (Nd isotopic composition, trace-element geochemistry, U-Pb dating) sand-sized apatite single-grain analysis.
The basin offers an excellent opportunity to investigate environmental signal propagation into the sedimentary record because comprehensive stratigraphic and seismic datasets can be combined with present research results. The bulk-rock clay-sized fraction epsilon Nd values of well-cutting samples from one well on the northern basin slope remained stable at similar to-9.7 from 27 to 19 Ma but increased after 19 Ma to similar to-9.1. In contrast, apatite single-grain distributions, which were extracted from 22 drill-core samples, changed significantly around 23.3 Ma from apatites dominantly from low-grade (<upper amphibolite-facies) metamorphic sources with Permo-Mesozoic and late Variscan U-Pb ages and epsilon Nd values of -4.4 to dominantly high-grade metamorphic apatites with late Variscan U-Pb ages and epsilon Nd values of -2.2.
The change in apatite single-grain distributions at 23.3 Ma is interpreted to result from the exposure of a new Upper Austroalpine source nappe with less negative epsilon Nd values triggered by the ongoing Tauern Window exhumation. Combining these data with the clay-sized bulk-rock epsilon Nd values reveals that the provenance changed 4-5 Myrs later at 19 Ma in the clay-sized fraction.
Reasons for the delayed provenance-change recording are rooted in the characteristics of the applied methods.
Whereas single-grain distributions of orogen-wide sediment-routing systems can be dominated by geographically small areas with high erosion and mineral fertility rates, bulk-rock methods integrate over the entire drainage basin, thus diminishing extreme values. Hence, by combining these two methods, spatial information are uncovered, enabling a previously unattained understanding of the underlying environmental change.
Technological progress allows for producing ever more complex predictive models on the basis of increasingly big datasets. For risk management of natural hazards, a multitude of models is needed as basis for decision-making, e.g. in the evaluation of observational data, for the prediction of hazard scenarios, or for statistical estimates of expected damage. The question arises, how modern modelling approaches like machine learning or data-mining can be meaningfully deployed in this thematic field. In addition, with respect to data availability and accessibility, the trend is towards open data. Topic of this thesis is therefore to investigate the possibilities and limitations of machine learning and open geospatial data in the field of flood risk modelling in the broad sense. As this overarching topic is broad in scope, individual relevant aspects are identified and inspected in detail.
A prominent data source in the flood context is satellite-based mapping of inundated areas, for example made openly available by the Copernicus service of the European Union. Great expectations are directed towards these products in scientific literature, both for acute support of relief forces during emergency response action, and for modelling via hydrodynamic models or for damage estimation. Therefore, a focus of this work was set on evaluating these flood masks. From the observation that the quality of these products is insufficient in forested and built-up areas, a procedure for subsequent improvement via machine learning was developed. This procedure is based on a classification algorithm that only requires training data from a particular class to be predicted, in this specific case data of flooded areas, but not of the negative class (dry areas). The application for hurricane Harvey in Houston shows the high potential of this method, which depends on the quality of the initial flood mask.
Next, it is investigated how much the predicted statistical risk from a process-based model chain is dependent on implemented physical process details. Thereby it is demonstrated what a risk study based on established models can deliver. Even for fluvial flooding, such model chains are already quite complex, though, and are hardly available for compound or cascading events comprising torrential rainfall, flash floods, and other processes. In the fourth chapter of this thesis it is therefore tested whether machine learning based on comprehensive damage data can offer a more direct path towards damage modelling, that avoids explicit conception of such a model chain. For that purpose, a state-collected dataset of damaged buildings from the severe El Niño event 2017 in Peru is used. In this context, the possibilities of data-mining for extracting process knowledge are explored as well. It can be shown that various openly available geodata sources contain useful information for flood hazard and damage modelling for complex events, e.g. satellite-based rainfall measurements, topographic and hydrographic information, mapped settlement areas, as well as indicators from spectral data. Further, insights on damaging processes are discovered, which mainly are in line with prior expectations. The maximum intensity of rainfall, for example, acts stronger in cities and steep canyons, while the sum of rain was found more informative in low-lying river catchments and forested areas. Rural areas of Peru exhibited higher vulnerability in the presented study compared to urban areas. However, the general limitations of the methods and the dependence on specific datasets and algorithms also become obvious.
In the overarching discussion, the different methods – process-based modelling, predictive machine learning, and data-mining – are evaluated with respect to the overall research questions. In the case of hazard observation it seems that a focus on novel algorithms makes sense for future research. In the subtopic of hazard modelling, especially for river floods, the improvement of physical models and the integration of process-based and statistical procedures is suggested. For damage modelling the large and representative datasets necessary for the broad application of machine learning are still lacking. Therefore, the improvement of the data basis in the field of damage is currently regarded as more important than the selection of algorithms.
Increasing arctic coastal erosion rates imply a greater release of sediments and organic matter into the coastal zone. With 213 sediment samples taken around Herschel Island-Qikiqtaruk, Canadian Beaufort Sea, we aimed to gain new insights on sediment dynamics and geochemical properties of a shallow arctic nearshore zone. Spatial characteristics of nearshore sediment texture (moderately to poorly sorted silt) are dictated by hydrodynamic processes, but ice-related processes also play a role. We determined organic matter (OM) distribution and inferred the origin and quality of organic carbon by C/N ratios and stable carbon isotopes delta C-13. The carbon content was higher offshore and in sheltered areas (mean: 1.0 wt.%., S.D.: 0.9) and the C/N ratios also showed a similar spatial pattern (mean: 11.1, S.D.: 3.1), while the delta C-13 (mean: -26.4 parts per thousand VPDB, S.D.: 0.4) distribution was more complex. We compared the geochemical parameters of our study with terrestrial and marine samples from other studies using a bootstrap approach. Sediments of the current study contained 6.5 times and 1.8 times less total organic carbon than undisturbed and disturbed terrestrial sediments, respectively. Therefore, degradation of OM and separation of carbon pools take place on land and continue in the nearshore zone, where OM is leached, mineralized, or transported beyond the study area.
ArcticBeach v1.0
(2022)
In the Arctic, air temperatures are increasing and sea ice is declining, resulting in larger waves and a longer open water season, all of which intensify the thaw and erosion of ice-rich coasts. Climate change has been shown to increase the rate of Arctic coastal erosion, causing problems for Arctic cultural heritage, existing industrial, military, and civil infrastructure, as well as changes in nearshore biogeochemistry. Numerical models that reproduce historical and project future Arctic erosion rates are necessary to understand how further climate change will affect these problems, and no such model yet exists to simulate the physics of erosion on a pan-Arctic scale. We have coupled a bathystrophic storm surge model to a simplified physical erosion model of a permafrost coastline. This Arctic erosion model, called ArcticBeach v1.0, is a first step toward a physical parameterization of Arctic shoreline erosion for larger-scale models. It is forced by wind speed and direction, wave period and height, sea surface temperature, all of which are masked during times of sea ice cover near the coastline. Model tuning requires observed historical retreat rates (at least one value), as well as rough nearshore bathymetry. These parameters are already available on a pan-Arctic scale. The model is validated at three study sites at 1) Drew Point (DP), Alaska, 2) Mamontovy Khayata (MK), Siberia, and 3) Veslebogen Cliffs, Svalbard. Simulated cumulative retreat rates for DP and MK respectively (169 and 170 m) over the time periods studied at each site (2007-2016, and 1995-2018) are found to the same order of magnitude as observed cumulative retreat (172 and 120 m). The rocky Veslebogen cliffs have small observed cumulative retreat rates (0.05 m over 2014-2016), and our model was also able to reproduce this same order of magnitude of retreat (0.08 m). Given the large differences in geomorphology between the study sites, this study provides a proof-of-concept that ArcticBeach v1.0 can be applied on very different permafrost coastlines. ArcticBeach v1.0 provides a promising starting point to project retreat of Arctic shorelines, or to evaluate historical retreat in places that have had few observations.
Assessment of climate change impact on discharge of the lakhmass catchment (Northwest Tunisia)
(2022)
The Mediterranean region is increasingly recognized as a climate change hotspot but is highly underrepresented in hydrological climate change studies. This study aims to investigate the climate change effects on the hydrology of Lakhmass catchment in Tunisia. Lakhmass catchment is a part of the Medium Valley of Medjerda in northwestern Tunisia that drains an area of 126 km(2). First, the Hydrologiska Byrans Vattenbalansavdelning light (HBV-light) model was calibrated and validated successfully at a daily time step to simulate discharge during the 1981-1986 period. The Nash Sutcliffe Efficiency and Percent bias (NSE, PBIAS) were (0.80, +2.0%) and (0.53, -9.5%) for calibration (September 1982-August 1984) and validation (September 1984-August 1986) periods, respectively. Second, HBV-light model was considered as a predictive tool to simulate discharge in a baseline period (1981-2009) and future projections using data (precipitation and temperature) from thirteen combinations of General Circulation Models (GCMs) and Regional Climatic Models (RCMs). We used two trajectories of Representative Concentration Pathways, RCP4.5 and RCP8.5, suggested by the Intergovernmental Panel on Climate Change (IPCC). Each RCP is divided into three projection periods: near-term (2010-2039), mid-term (2040-2069) and long-term (2070-2099). For both scenarios, a decrease in precipitation and discharge will be expected with an increase in air temperature and a reduction in precipitation with almost 5% for every +1 degrees C of global warming. By long-term (2070-2099) projection period, results suggested an increase in temperature with about 2.7 degrees C and 4 degrees C, and a decrease in precipitation of approximately 7.5% and 15% under RCP4.5 and RCP8.5, respectively. This will likely result in a reduction of discharge of 12.5% and 36.6% under RCP4.5 and RCP8.5, respectively. This situation calls for early climate change adaptation measures under a participatory approach, including multiple stakeholders and water users.
In this study we present a novel method for the automatic detection of minerals and elements using hyperspectral transmittance imaging microscopy measurements of complete thin sections (HyperTIM).
This is accomplished by using a hyperspectral camera system that operates in the visible and near-infrared (VNIR) range with a specifically designed sample holder, scanning setup, and a microscope lens.
We utilize this method on a monazite ore thin section from Steenkampskraal (South Africa), which we analyzed for the rare earth element (REE)-bearing mineral monazite ((Ce,Nd,La)PO4), with high concentrations of Nd. The transmittance analyses with the hyperspectral VNIR camera can be used to identify REE minerals and Nd in thin sections.
We propose a three-point band depth index, the Nd feature depth index (NdFD), and its related product the Nd band depth index (NdBDI), which enables automatic mineral detection and classification for the Nd-bearing monazites in thin sections. In combination with the average concentration of the relative Nd content, it permits a destruction-free, total concentration calculation for Nd across the entire thin section.
During the last 5 Ma the Earth's ocean-atmosphere system passed through several major transitions, many of which are discussed as possible triggers for human evolution. A classic in this context is the possible influence of the closure of the Panama Strait, the intensification of Northern Hemisphere Glaciation, a stepwise increase in aridity in Africa, and the first appearance of the genus Homo about 2.5 - 2.7 Ma ago. Apart from the fact that the correlation between these events does not necessarily imply causality, many attempts to establish a relationship between climate and evolution fail due to the challenge of precisely localizing an a priori unknown number of changes potentially underlying complex climate records. The kernel-based Bayesian inference approach applied here allows inferring the location, generic shape, and temporal scale of multiple transitions in established records of Plio-Pleistocene African climate. By defining a transparent probabilistic analysis strategy, we are able to identify conjoint changes occurring across the investigated terrigenous dust records from Ocean Drilling Programme (ODP) sites in the Atlantic Ocean (ODP 659), Arabian (ODP 721/722) and Mediterranean Sea (ODP 967). The study indicates a two-step transition in the African climate proxy records at (2.35-2.10) Ma and (1.70 - 1.50) Ma, that may be associated with the reorganization of the Hadley-Walker Circulation. .
Beyond CO2 equivalence
(2022)
In this article we review the physical and chemical properties of methane (CH4) relevant to impacts on climate, ecosystems, and air pollution, and examine the extent to which this is reflected in climate and air pollution governance. Although CH4 is governed under the UNFCCC climate regime, its treatment there is limited to the ways in which it acts as a "CO2 equivalent" climate forcer on a 100-year time frame. The UNFCCC framework neglects the impacts that CH4 has on near-term climate, as well its impacts on human health and ecosystems, which are primarily mediated by methane's role as a precursor to tropospheric ozone. Frameworks for air quality governance generally address tropospheric ozone as a pollutant, but do not regulate CH4 itself. Methane's climate and air quality impacts, together with its alarming rise in atmospheric concentrations in recent years, make it clear that mitigation of CH4 emissions needs to be accelerated globally. We examine challenges and opportunities for further progress on CH4 mitigation within the international governance landscapes for climate change and air pollution.
Sulfate reduction is the quantitatively most important process to degrade organic matter in anoxic marine sediment and has been studied intensively in a variety of settings. Guaymas Basin, a young marginal ocean basin, offers the unique opportunity to study sulfate reduction in an environment characterized by organic-rich sediment, high sedimentation rates, and high geothermal gradients (100-958 degrees C km(-1)). We measured sulfate reduction rates (SRR) in samples taken during the International Ocean Discovery Program (IODP) Expedition 385 using incubation experiments with radiolabeled (SO42-)-S-35 carried out at in situ pressure and temperature. The highest SRR (387 nmol cm(-3) d(-1)) was recorded in near-surface sediments from Site U1548C, which had the steepest geothermal gradient (958 degrees C km(-1)). At this site, SRR were generally over an order of magnitude higher than at similar depths at other sites (e.g., 387-157 nmol cm(-3) d(-1) at 1.9 mbsf from Site U1548C vs. 46-1.0 nmol cm(-3) d(-1) at 2.1 mbsf from Site U1552B). Site U1546D is characterized by a sill intrusion, but it had already reached thermal equilibrium and SRR were in the same range as nearby Site U1545C, which is minimally affected by sills. The wide temperature range observed at each drill site suggests major shifts in microbial community composition with very different temperature optima but awaits confirmation by molecular biological analyses. At the transition between the mesophilic and thermophilic range around 40 degrees C-60 degrees C, sulfate-reducing activity appears to be decreased, particularly in more oligotrophic settings, but shows a slight recovery at higher temperatures.
The creation of building exposure models for seismic risk assessment is frequently challenging due to the lack of availability of detailed information on building structures. Different strategies have been developed in recent years to overcome this, including the use of census data, remote sensing imagery and volunteered graphic information (VGI). This paper presents the development of a building-by-building exposure model based exclusively on openly available datasets, including both VGI and census statistics, which are defined at different levels of spatial resolution and for different moments in time. The initial model stemming purely from building-level data is enriched with statistics aggregated at the neighbourhood and city level by means of a Monte Carlo simulation that enables the generation of full realisations of damage estimates when using the exposure model in the context of an earthquake scenario calculation. Though applicable to any other region of interest where analogous datasets are available, the workflow and approach followed are explained by focusing on the case of the German city of Cologne, for which a scenario earthquake is defined and the potential damage is calculated. The resulting exposure model and damage estimates are presented, and it is shown that the latter are broadly consistent with damage data from the 1978 Albstadt earthquake, notwithstanding the differences in the scenario. Through this real-world application we demonstrate the potential of VGI and open data to be used for exposure modelling for natural risk assessment, when combined with suitable knowledge on building fragility and accounting for the inherent uncertainties.
Improving nitrogen (N) status in European water bodies is a pressing issue. N levels depend not only on current but also past N inputs to the landscape, that have accumulated through time in legacy stores (e.g., soil, groundwater).
Catchment-scale N models, that are commonly used to investigate in-stream N levels, rarely examine the magnitude and dynamics of legacy components.
This study aims to gain a better understanding of the long-term fate of the N inputs and its uncertainties, using a legacy-driven N model (ELEMeNT) in Germany's largest national river basin (Weser; 38,450 km(2)) over the period 1960-2015.
We estimate the nine model parameters based on a progressive constraining strategy, to assess the value of different observational data sets.
We demonstrate that beyond in-stream N loading, soil N content and in-stream N concentration allow to reduce the equifinality in model parameterizations.
We find that more than 50% of the N surplus denitrifies (1480-2210 kg ha(-1)) and the stream export amounts to around 18% (410-640 kg ha(-1)), leaving behind as much as around 230-780 kg ha(-1) of N in the (soil) source zone and 10-105 kg ha(-1) in the subsurface.
A sensitivity analysis reveals the importance of different factors affecting the residual uncertainties in simulated N legacies, namely hydrologic travel time, denitrification rates, a coefficient characterizing the protection of organic N in source zone and N surplus input.
Our study calls for proper consideration of uncertainties in N legacy characterization, and discusses possible avenues to further reduce the equifinality in water quality modeling.
Climate Benefits of Cleaner Energy Transitions in East and South Asia Through Black Carbon Reduction
(2022)
The state of air pollution has historically been tightly linked to how we produce and use energy. Air pollutant emissions over Asia are now changing rapidly due to cleaner energy transitions; however, magnitudes of benefits for climate and air quality remain poorly quantified. The associated risks involve adverse health impacts, reduced agricultural yields, reduced freshwater availability, contributions to climate change, and economic costs. We focus particularly on climate benefits of energy transitions by making first-time use of two decades of high quality observations of atmospheric loading of light-absorbing black carbon (BC) over Kanpur (South Asia) and Beijing (East Asia) and relating these observations to changing energy, emissions, and economic trends in India and China. Our analysis reveals that absorption aerosol optical depth (AAOD) due to BC has decreased substantially, by 40% over Kanpur and 60% over Beijing between 2001 and 2017, and thus became decoupled from regional economic growth. Furthermore, the resultant decrease in BC emissions and BC AAOD over Asia is regionally coherent and occurs primarily due to transitions into cleaner energies (both renewables and fossil fuels) and not due to the decrease in primary energy supply or decrease in use of fossil use and biofuels and waste. Model simulations show that BC aerosols alone contribute about half of the surface temperature change (warming) of the total forcing due to greenhouse gases, natural and internal variability, and aerosols, thus clearly revealing the climate benefits due to a reduction in BC emissions, which would significantly reduce global warming. However, this modeling study excludes responses from natural variability, circulation, and sea ice responses, which cause relatively strong temperature fluctuations that may mask signals from BC aerosols. Our findings show additional benefits for climate (beyond benefits of CO2 reduction) and for several other issues of sustainability over South and East Asia, provide motivation for ongoing cleaner energy production, and consumption transitions, especially when they are associated with reduced emissions of air pollutants. Such an analysis connecting the trends in energy transitions and aerosol absorption loading, unavailable so far, is crucial for simulating the aerosol climate impacts over Asia which is quite uncertain.
Understanding the influence of climate change and population pressure on human conflict remains a critically important topic in the social sciences. Long-term records that evaluate these dynamics across multiple centuries and outside the range of modern climatic variation are especially capable of elucidating the relative effect of-and the interaction between-climate and demography. This is crucial given that climate change may structure population growth and carrying capacity, while both climate and population influence per capita resource availability. This study couples paleoclimatic and demographic data with osteological evaluations of lethal trauma from 149 directly accelerator mass spectrometry C-14-dated individuals from the Nasca highland region of Peru. Multiple local and supraregional precipitation proxies are combined with a summed probability distribution of 149 C-14 dates to estimate population dynamics during a 700-y study window. Counter to previous findings, our analysis reveals a precipitous increase in violent deaths associated with a period of productive and stable climate, but volatile population dynamics. We conclude that favorable local climate conditions fostered population growth that put pressure on the marginal and highly circumscribed resource base, resulting in violent resource competition that manifested in over 450 y of internecine warfare. These findings help support a general theory of intergroup violence, indicating that relative resource scarcity-whether driven by reduced resource abundance or increased competition-can lead to violence in subsistence societies when the outcome is lower per capita resource availability.
Devolatilization of subducting lithologies liberates COH-fluids. These may become partially sequestered in peridotites in the slab and the overlying forearc mantle, affecting the cycling of volatiles and fluid mobile elements in subduction zones. Here we assess the magnitudes, timescales and mechanism of channelized injection of COH-fluids doped with Ca-aq(2+), Sr-aq(2+) and Ba-aq(2+) into the dry forearc mantle by performing piston cylinder experiments between 1-2.5 GPa and 600-700 degrees C. Cylindrical cores of natural spinel-bearing harzburgites were used as starting materials. Based on mineral assemblage and composition three reaction zones are distinguishable from the rim towards the core of primary olivine and orthopyroxene grains. Zone 1 contains carbonates + quartz +/- kyanite and zone 2 contains carbonates + talc +/- chlorite. Olivine is further replaced in zone 3 by either antigorite+ magnesite or magnesite +talc within or above antigorite stability, respectively. Orthopyroxene is replaced in zone 3 by talc + chlorite. Mineral assemblages and the compositions of secondary minerals depend on fluid composition and the replaced primary silicate. The extent of alteration depends on fluid CO2 content and fluid/rock-ratio, and is further promoted by fluid permeable reaction zones and reaction driven cracking. Our results show that COH-fluid induced metasomatism of the forearc mantle is self-perpetuating and efficient at sequestering Ca-aq(2+), Sr-aq(2+), Ba-aq(2+) and CO2aq into newly formed carbonates. This process is fast with 90% of the available C sequestered and nearly 50% of the initial minerals altered at 650 degrees C, 2 GPa within 55 h. The dissolution of primary silicates under high COH-fluid/rock-ratios, as in channelized fluid flow, enriches SiO2aq in the fluid, while CO2aq is sequestered into carbonates. In an open system, the remaining CO2-depleted, Si-enriched aqueous fluid may cause Si-metasomatism in the forearc further away from the injection of the COH-fluid into peridotite.
Here I present a comparison between two of the most widely used reduced-complexity models for the representation of sediment transport and deposition processes, namely the transport-limited (or TL) model and the under-capacity (or xi-q) model more recently developed by Davy and Lague (2009). Using both models, I investigate the behavior of a sedimentary continental system of length L fed by a fixed sedimentary flux from a catchment of size A(0) in a nearby active orogen through which sediments transit to a fixed base level representing a large river, a lake or an ocean. This comparison shows that the two models share the same steady-state solution, for which I derive a simple 1D analytical expression that reproduces the major features of such sedimentary systems: a steep fan that connects to a shallower alluvial plain. The resulting fan geometry obeys basic observational constraints on fan size and slope with respect to the upstream drainage area, A(0). The solution is strongly dependent on the size of the system, L, in comparison to a distance L-0, which is determined by the size of A(0), and gives rise to two fundamentally different types of sedimentary systems: a constrained system where L < L-0 and open systems where L > L-0. I derive simple expressions that show the dependence of the system response time on the system characteristics, such as its length, the size of the upstream catchment area, the amplitude of the incoming sedimentary flux and the respective rate parameters (diffusivity or erodibility) for each of the two models. I show that the xi-q model predicts longer response times. I demonstrate that although the manner in which signals propagates through the sedimentary system differs greatly between the two models, they both predict that perturbations that last longer than the response time of the system can be recorded in the stratigraphy of the sedimentary system and in particular of the fan. Interestingly, the xi-q model predicts that all perturbations in the incoming sedimentary flux will be transmitted through the system, whereas the TL model predicts that rapid perturbations cannot. I finally discuss why and under which conditions these differences are important and propose observational ways to determine which of the two models is most appropriate to represent natural systems.
In this study, 3-D models of P-wave velocity (Vp) and P-wave and S-wave ratio (Vp/Vs) of the crust and upper mantle in the Eastern and eastern Southern Alps (northern Italy and southern Austria) were calculated using local earthquake tomography (LET). The data set includes high-quality arrival times from well-constrained hypocenters observed by the dense, temporary seismic networks of the AlpArray AASN and SWATH-D. The resolution of the LET was checked by synthetic tests and analysis of the model resolution matrix. The small inter-station spacing (average of similar to 15 km within the SWATH-D network) allowed us to image crustal structure at unprecedented resolution across a key part of the Alps. The derived P velocity model revealed a highly heterogeneous crustal structure in the target area. One of the main findings is that the lower crust is thickened, forming a bulge at 30-50 km depth just south of and beneath the Periadriatic Fault and the Tauern Window. This indicates that the lower crust decoupled both from its mantle substratum as well as from its upper crust. The Moho, taken to be the iso-velocity contour of Vp = 7.25 km/s, agrees with the Moho depth from previous studies in the European and Adriatic forelands. It is shallower on the Adriatic side than on the European side. This is interpreted to indicate that the European Plate is subducted beneath the Adriatic Plate in the Eastern and eastern Southern Alps.
Localisation of deformation is a ubiquitous feature in continental rift dynamics and observed across drastically different time and length scales. This thesis comprises one experimental and two numerical modelling studies investigating strain localisation in (1) a ductile shear zone induced by a material heterogeneity and (2) in an active continental rift setting. The studies are related by the fact that the weakening mechanisms on the crystallographic and grain size scale enable bulk rock weakening, which fundamentally enables the formation of shear zones, continental rifts and hence plate tectonics. Aiming to investigate the controlling mechanisms on initiation and evolution of a shear zone, the torsion experiments of the experimental study were conducted in a Patterson type apparatus with strong Carrara marble cylinders with a weak, planar Solnhofen limestone inclusion. Using state-of-the-art numerical modelling software, the torsion experiments were simulated to answer questions regarding localisation procedure like stress distribution or the impact of rheological weakening. 2D numerical models were also employed to integrate geophysical and geological data to explain characteristic tectonic evolution of the Southern and Central Kenya Rift. Key elements of the numerical tools are a randomized initial strain distribution and the usage of strain softening. During the torsion experiments, deformation begins to localise at the limestone inclusion tips in a process zone, which propagates into the marble matrix with increasing deformation until a ductile shear zone is established. Minor indicators for coexisting brittle deformation are found close to the inclusion tip and presumed to slightly facilitate strain localisation besides the dominant ductile deformation processes. The 2D numerical model of the torsion experiment successfully predicts local stress concentration and strain rate amplification ahead of the inclusion in first order agreement with the experimental results. A simple linear parametrization of strain weaking enables high accuracy reproduction of phenomenological aspects of the observed weakening. The torsion experiments suggest that loading conditions do not affect strain localisation during high temperature deformation of multiphase material with high viscosity contrasts. A numerical simulation can provide a way of analysing the process zone evolution virtually and extend the examinable frame. Furthermore, the nested structure and anastomosing shape of an ultramylonite band was mimicked with an additional second softening step. Rheological weakening is necessary to establish a shear zone in a strong matrix around a weak inclusion and for ultramylonite formation.
Such strain weakening laws are also incorporated into the numerical models of the
Southern and Central Kenya Rift that capture the characteristic tectonic evolution. A three-stage early rift evolution is suggested that starts with (1) the accommodation of strain by a single border fault and flexure of the hanging-wall crust, after which (2) faulting in the hanging-wall and the basin centre increases before (3) the early-stage asymmetry is lost and basinward localisation of deformation occurs. Along-strike variability of rifts can be produced by modifying the initial random noise distribution. In summary, the three studies address selected aspects of the broad range of mechanisms and processes that fundamentally enable the deformation of rock and govern the localisation patterns across the scales. In addition to the aforementioned results, the first and second manuscripts combined, demonstrate a procedure to find new or improve on existing numerical formulations for specific rheologies and their dynamic weakening. These formulations are essential in addressing rock deformation from the grain to the global scale. As within the third study of this thesis, where geodynamic controls on the evolution of a rift were examined and acquired by the integration of geological and geophysical data into a numerical model.
In an ocean-continent subduction zone, the assessment of the lithospheric thermal state is essential to determine the controls of the deformation within the upper plate and the dip angle of the subducting lithosphere. In this study, we evaluate the degree of influence of both the configuration of the upper plate (i.e., thickness and composition of the rock units) and variations of the subduction angle on the lithospheric thermal field of the southern Central Andes (29 degrees-39 degrees S). Here, the subduction angle increases from subhorizontal (5 degrees) north of 33 degrees S to steep (similar to 30 degrees) in the south. We derived the 3D temperature and heat flow distribution of the lithosphere in the southern Central Andes considering conversion of S wave tomography to temperatures together with steady-state conductive thermal modeling. We found that the orogen is overall warmer than the forearc and the foreland and that the lithosphere of the northern part of the foreland appears colder than its southern counterpart. Sedimentary blanketing and the thickness of the radiogenic crust exert the main control on the shallow thermal field (<50km depth). Specific conditions are present where the oceanic slab is relatively shallow (<85 km depth) and the radiogenic crust is thin. This configuration results in relatively colder temperatures compared to regions where the radiogenic crust is thick and the slab is steep. At depths >50km, the temperatures of the overriding plate are mainly controlled by the mantle heat input and the subduction angle. The thermal field of the upper plate likely preserves the flat subduction angle and influences the spatial distribution of shortening.
Data recorded by distributed acoustic sensing (DAS) along an optical fibre sample the spatial and temporal properties of seismic wavefields at high spatial density. Often leading to massive amount of data when collected for seismic monitoring along many kilometre long cables. The spatially coherent signals from weak seismic arrivals within the data are often obscured by incoherent noise. We present a flexible and computationally efficient filtering technique, which makes use of the dense spatial and temporal sampling of the data and that can handle the large amount of data. The presented adaptive frequency-wavenumber filter suppresses the incoherent seismic noise while amplifying the coherent wavefield. We analyse the response of the filter in time and spectral domain, and we demonstrate its performance on a noisy data set that was recorded in a vertical borehole observatory showing active and passive seismic phase arrivals. Lastly, we present a performant open-source software implementation enabling real-time filtering of large DAS data sets.
We use a dense seismic network on the Reykjanes Peninsula, Iceland, to image a group of earthquakes at 10-12 km depth, 2 km north-east of 2021 Fagradalsfjall eruption site. These deep earthquakes have a lower frequency content compared to earthquakes located in the upper, brittle crust and are similar to deep long period (DLP) seismicity observed at other volcanoes in Iceland and around the world. We observed several swarms of DLP earthquakes between the start of the study period (June 2020) and the initiation of the 3-week-long dyke intrusion that preceded the eruption in March 2021. During the eruption, DLP earthquake swarms returned 1 km SW of their original location during periods when the discharge rate or fountaining style of the eruption changed. The DLP seismicity is therefore likely to be linked to the magma plumbing system beneath Fagradalsfjall. However, the DLP seismicity occurred similar to 5 km shallower than where petrological modelling places the near-Moho magma storage region in which the Fagradalsfjall lava was stored. We suggest that the DLP seismicity was triggered by the exsolution of CO2-rich fluids or the movement of magma at a barrier to the transport of melt in the lower crust. Increased flux through the magma plumbing system during the eruption likely adds to the complexity of the melt migration process, thus causing further DLP seismicity, despite a contemporaneous magma channel to the surface.
Fast-localized electron loss, resulting from interactions with electromagnetic ion cyclotron (EMIC) waves, can produce deepening minima in phase space density (PSD) radial profiles. Here, we perform a statistical analysis of local PSD minima to quantify how readily these are associated with radiation belt depletions. The statistics of PSD minima observed over a year are compared to the Versatile Electron Radiation Belts (VERB) simulations, both including and excluding EMIC waves. The observed minima distribution can only be achieved in the simulation including EMIC waves, indicating their importance in the dynamics of the radiation belts. By analyzing electron flux depletions in conjunction with the observed PSD minima, we show that, in the heart of the outer radiation belt (L* < 5), on average, 53% of multi-MeV electron depletions are associated with PSD minima, demonstrating that fast localized loss by interactions with EMIC waves are a common and crucial process for ultra-relativistic electron populations.
Diet analysis of bats killed at wind turbines suggests large-scale losses of trophic interactions
(2022)
Agricultural practice has led to landscape simplification and biodiversity decline, yet recently, energy-producing infrastructures, such as wind turbines, have been added to these simplified agroecosystems, turning them into multi-functional energy-agroecosystems. Here, we studied the trophic interactions of bats killed at wind turbines using a DNA metabarcoding approach to shed light on how turbine-related bat fatalities may possibly affect local habitats. Specifically, we identified insect DNA in the stomachs of common noctule bats (Nyctalus noctula) killed by wind turbines in Germany to infer in which habitats these bats hunted. Common noctule bats consumed a wide variety of insects from different habitats, ranging from aquatic to terrestrial ecosystems (e.g., wetlands, farmland, forests, and grasslands). Agricultural and silvicultural pest insects made up about 20% of insect species consumed by the studied bats. Our study suggests that the potential damage of wind energy production goes beyond the loss of bats and the decline of bat populations. Bat fatalities at wind turbines may lead to the loss of trophic interactions and ecosystem services provided by bats, which may add to the functional simplification and impaired crop production, respectively, in multi-functional ecosystems.
Diet analysis of bats killed at wind turbines suggests large-scale losses of trophic interactions
(2022)
Agricultural practice has led to landscape simplification and biodiversity decline, yet recently, energy-producing infrastructures, such as wind turbines, have been added to these simplified agroecosystems, turning them into multi-functional energy-agroecosystems. Here, we studied the trophic interactions of bats killed at wind turbines using a DNA metabarcoding approach to shed light on how turbine-related bat fatalities may possibly affect local habitats. Specifically, we identified insect DNA in the stomachs of common noctule bats (Nyctalus noctula) killed by wind turbines in Germany to infer in which habitats these bats hunted. Common noctule bats consumed a wide variety of insects from different habitats, ranging from aquatic to terrestrial ecosystems (e.g., wetlands, farmland, forests, and grasslands). Agricultural and silvicultural pest insects made up about 20% of insect species consumed by the studied bats. Our study suggests that the potential damage of wind energy production goes beyond the loss of bats and the decline of bat populations. Bat fatalities at wind turbines may lead to the loss of trophic interactions and ecosystem services provided by bats, which may add to the functional simplification and impaired crop production, respectively, in multi-functional ecosystems.
Understanding catchment controls on catchment solute export is a prerequisite for water quality management. StorAge Selection (SAS) functions encapsulate essential information about catchment functioning in terms of discharge selection preference and solute export dynamics. However, they lack information on the spatial origin of solutes when applied at the catchment scale, thereby limiting our understanding of the internal (subcatchment) functioning. Here, we parameterized SAS functions in a spatially explicit way to understand the internal catchment responses and transport dynamics of reactive dissolved nitrate (N-NO3). The model was applied in a nested mesoscale catchment (457 km(2)), consisting of a mountainous partly forested, partly agricultural subcatchment, a middle-reach forested subcatchment, and a lowland agricultural subcatchment. The model captured flow and nitrate concentration dynamics not only at the catchment outlet but also at internal gauging stations. Results reveal disparate subsurface mixing dynamics and nitrate export among headwater and lowland subcatchments. The headwater subcatchment has high seasonal variation in subsurface mixing schemes and younger water in discharge, while the lowland subcatchment has less pronounced seasonality in subsurface mixing and much older water in discharge. Consequently, nitrate concentration in discharge from the headwater subcatchment shows strong seasonality, whereas that from the lowland subcatchment is stable in time. The temporally varying responses of headwater and lowland subcatchments alternate the dominant contribution to nitrate export in high and low-flow periods between subcatchments. Overall, our results demonstrate that the spatially explicit SAS modeling provides useful information about internal catchment functioning, helping to develop or evaluate spatial management practices.
Draft Genome Sequence of Nocardioides alcanivorans NGK65(T), a Hexadecane-Degrading Bacterium
(2022)
The Gram-positive bacterium Nocardioides alcanivorans NGK65(T) was isolated from plastic-polluted soil and cultivated on medium with polyethylene as the single carbon source. Nanopore sequencing revealed the presence of candidate enzymes for the biodegradation of polyethylene. Here, we report the draft genome of this newly described member of the terrestrial plastisphere.
The Arctic is greatly affected by climate change. Increasing air temperatures drive permafrost thaw and an increase in coastal erosion and river discharge. This results in a greater input of sediment and organic matter into nearshore waters, impacting ecosystems by reducing light transmission through the water column and altering biogeochemistry. This potentially results in impacts on the subsistence economy of local people as well as the climate due to the transformation of suspended organic matter into greenhouse gases. Even though the impacts of increased suspended sediment concentrations and turbidity in the Arctic nearshore zone are well-studied, the mechanisms underpinning this increase are largely unknown. Wave energy and tides drive the level of turbidity in the temperate and tropical parts of the world, and this is generally assumed to also be the case in the Arctic. However, the tidal range is considerably lower in the Arctic, and processes related to the occurrence of permafrost have the potential to greatly contribute to nearshore turbidity. In this study, we use high-resolution satellite imagery alongside in situ and ERA5 reanalysis data of ocean and climate variables in order to identify the drivers of nearshore turbidity, along with its seasonality in the nearshore waters of Herschel Island Qikiqtaruk, in the western Canadian Arctic. Nearshore turbidity correlates well to wind direction, wind speed, significant wave height, and wave period. Nearshore turbidity is superiorly correlated to wind speed at the Beaufort Shelf compared to in situ measurements at Herschel Island Qikiqtaruk, showing that nearshore turbidity, albeit being of limited spatial extent, is influenced by large-scale weather and ocean phenomenons. We show that, in contrast to the temperate and tropical ocean, freshly eroded material is the predominant driver of nearshore turbidity in the Arctic, rather than resuspension, which is caused by the vulnerability of permafrost coasts to thermo-erosion.
The Arctic is greatly affected by climate change. Increasing air temperatures drive permafrost thaw and an increase in coastal erosion and river discharge. This results in a greater input of sediment and organic matter into nearshore waters, impacting ecosystems by reducing light transmission through the water column and altering biogeochemistry. This potentially results in impacts on the subsistence economy of local people as well as the climate due to the transformation of suspended organic matter into greenhouse gases. Even though the impacts of increased suspended sediment concentrations and turbidity in the Arctic nearshore zone are well-studied, the mechanisms underpinning this increase are largely unknown. Wave energy and tides drive the level of turbidity in the temperate and tropical parts of the world, and this is generally assumed to also be the case in the Arctic. However, the tidal range is considerably lower in the Arctic, and processes related to the occurrence of permafrost have the potential to greatly contribute to nearshore turbidity. In this study, we use high-resolution satellite imagery alongside in situ and ERA5 reanalysis data of ocean and climate variables in order to identify the drivers of nearshore turbidity, along with its seasonality in the nearshore waters of Herschel Island Qikiqtaruk, in the western Canadian Arctic. Nearshore turbidity correlates well to wind direction, wind speed, significant wave height, and wave period. Nearshore turbidity is superiorly correlated to wind speed at the Beaufort Shelf compared to in situ measurements at Herschel Island Qikiqtaruk, showing that nearshore turbidity, albeit being of limited spatial extent, is influenced by large-scale weather and ocean phenomenons. We show that, in contrast to the temperate and tropical ocean, freshly eroded material is the predominant driver of nearshore turbidity in the Arctic, rather than resuspension, which is caused by the vulnerability of permafrost coasts to thermo-erosion.
Effect of temperature on the densification of silicate melts to lower earth's mantle conditions
(2022)
Physical properties of silicate melts play a key role for global planetary dynamics, controlling for example volcanic eruption styles, mantle convection and elemental cycling in the deep Earth. They are significantly modified by structural changes at the atomic scale due to external parameters such as pressure and temperature or due to chemistry. Structural rearrangements such as 4- to 6-fold coordination change of Si with increasing depth may profoundly influence melt properties, but have so far mostly been studied at ambient temperature due to experimental difficulties. In order to investigate the structural properties of silicate melts and their densification mechanisms at conditions relevant to the deep Earth's interior, we studied haplo basaltic glasses and melts (albite-diopside composition) at high pressure and temperature conditions in resistively and laser-heated diamond anvil cells using X-ray absorption near edge structure spectroscopy. Samples were doped with 10 wt% of Ge, which is accessible with this experimental technique and which commonly serves as a structural analogue for the network forming cation Si. We acquired spectra on the Ge K edge up to 48 GPa and 5000 K and derived the average Ge-O coordination number NGe-O, and bond distance RGe-O as functions of pressure. Our results demonstrate a continuous transformation from tetrahedral to octahedral coordination between ca. 5 and 30 GPa at ambient temperature. Above 1600 K the data reveal a reduction of the pressure needed to complete conversion to octahedral coordination by ca. 30 %. The results allow us to determine the influence of temperature on the Si coordination number changes in natural melts in the Earth's interior. We propose that the complete transition to octahedral coordination in basaltic melts is reached at about 40 GPa, corresponding to a depth of ca. 1200 km in the uppermost lower mantle. At the core-mantle boundary (2900 km, 130 GPa, 3000 K) the existence of non-buoyant melts has been proposed to explain observed low seismic wave velocity features. Our results highlight that the melt composition can affect the melt density at such extreme conditions and may strongly influence the structural response.
In the past decade, sediment connectivity has become a widely recognized characteristic of a geomorphic system. However, the quantification of functional connectivity (i.e. connectivity which arises due to the actual occurrence of sediment transport processes) and its variation over space and time is still a challenge. In this context, this study assesses the effects of expected future phenomena in the context of climate change (i.e. glacier retreat, permafrost degradation or meteorological extreme events) on sediment transport dynamics in a glacierised Alpine basin. The study area is the Sulden river basin (drainage area 130 km(2)) in the Italian Alps, which is composed of two geomorphologically diverse sub-basins. Based on graph theory, we evaluated the spatio-temporal variations in functional connectivity in these two sub-basins. The graph-object, obtained by manually mapping sediment transport processes between landforms, was adapted to 6 different hydro-meteorological scenarios, which derive from combining base, heatwave and rainstorm conditions with snowmelt and glacier-melt periods. For each scenario and each sub-basin, the sediment transport network and related catchment characteristics were analysed. To compare the effects of the scenarios on functional connectivity, we introduced a connectivity degree, calculated based on the area of the landforms involved in sediment cascades. Results indicate that the area of the basin connected to its outlet in terms of sediment transport might feature a six-fold increase in case of rainstorm conditions compared to "average " meteorological conditions assumed for the base scenario. Furthermore, markedly different effects of climate change on sediment connectivity are expected between the two sub-catchments due to their contrasting morphological and lithological characteristics, in terms of relative importance of rainfall triggered colluvial processes vs temperature-driven proglacial fluvial dynamics.
The authors demonstrate that a vegetation system's ability to recover from disturbances-its resilience-can be estimated from its natural variability. Global patterns of resilience loss and gains since the early 1990s reveal shifts towards widespread resilience loss since the early 2000s.
The character and health of ecosystems worldwide is tightly coupled to changes in Earth's climate. Theory suggests that ecosystem resilience-the ability of ecosystems to resist and recover from external shocks such as droughts and fires-can be inferred from their natural variability. Here, we quantify vegetation resilience globally with complementary metrics based on two independent long-term satellite records. We first empirically confirm that the recovery rates from large perturbations can be closely approximated from internal vegetation variability across vegetation types and climate zones. On the basis of this empirical relationship, we quantify vegetation resilience continuously and globally from 1992 to 2017. Long-term vegetation resilience trends are spatially heterogeneous, with overall increasing resilience in the tropics and decreasing resilience at higher latitudes. Shorter-term trends, however, reveal a marked shift towards a global decline in vegetation resilience since the early 2000s, particularly in the equatorial rainforest belt.
In seismic risk assessment, the sources of uncertainty associated with building exposure modelling have not received as much attention as other components related to hazard and vulnerability. Conventional practices such as assuming absolute portfolio compositions (i.e., proportions per building class) from expert-based assumptions over aggregated data crudely disregard the contribution of uncertainty of the exposure upon earthquake loss models. In this work, we introduce the concept that the degree of knowledge of a building stock can be described within a Bayesian probabilistic approach that integrates both expert-based prior distributions and data collection on individual buildings. We investigate the impact of the epistemic uncertainty in the portfolio composition on scenario-based earthquake loss models through an exposure-oriented logic tree arrangement based on synthetic building portfolios. For illustrative purposes, we consider the residential building stock of Valparaiso (Chile) subjected to seismic ground-shaking from one subduction earthquake. We have found that building class reconnaissance, either from prior assumptions by desktop studies with aggregated data (top-down approach), or from building-by-building data collection (bottom-up approach), plays a fundamental role in the statistical modelling of exposure. To model the vulnerability of such a heterogeneous building stock, we require that their associated set of structural fragility functions handle multiple spectral periods. Thereby, we also discuss the relevance and specific uncertainty upon generating either uncorrelated or spatially cross-correlated ground motion fields within this framework. We successively show how various epistemic uncertainties embedded within these probabilistic exposure models are differently propagated throughout the computed direct financial losses. This work calls for further efforts to redesign desktop exposure studies, while also highlighting the importance of exposure data collection with standardized and iterative approaches.
A volcanic eruption is usually preceded by seismic precursors, but their interpretation and use for forecasting the eruption onset time remain a challenge. A part of the eruptive processes in open conduits of volcanoes may be similar to those encountered in geysers. Since geysers erupt more often, they are useful sites for testing new forecasting methods. We tested the application of Permutation Entropy (PE) as a robust method to assess the complexity in seismic recordings of the Strokkur geyser, Iceland. Strokkur features several minute-long eruptive cycles, enabling us to verify in 63 recorded cycles whether PE behaves consistently from one eruption to the next one. We performed synthetic tests to understand the effect of different parameter settings in the PE calculation. Our application to Strokkur shows a distinct, repeating PE pattern consistent with previously identified phases in the eruptive cycle. We find a systematic increase in PE within the last 15 s before the eruption, indicating that an eruption will occur. We quantified the predictive power of PE, showing that PE performs better than seismic signal strength or quiescence when it comes to forecasting eruptions.
In some catchments, the distribution of annual maximum streamflow shows heavy tail behavior, meaning the occurrence probability of extreme events is higher than if the upper tail decayed exponentially. Neglecting heavy tail behavior can lead to an underestimation of the likelihood of extreme floods and the associated risk. Partly contradictory results regarding the controls of heavy tail behavior exist in the literature and the knowledge is still very dispersed and limited. To better understand the drivers, we analyze the upper tail behavior and its controls for 480 catchments in Germany and Austria over a period of more than 50 years. The catchments span from quickly reacting mountain catchments to large lowland catchments, allowing for general conclusions. We compile a wide range of event and catchment characteristics and investigate their association with an indicator of the tail heaviness of flood distributions, namely the shape parameter of the GEV distribution. Following univariate analyses of these characteristics, along with an evaluation of different aggregations of event characteristics, multiple linear regression models, as well as random forests, are constructed. A novel slope indicator, which represents the relation between the return period of flood peaks and event characteristics, captures the controls of heavy tails best. Variables describing the catchment response are found to dominate the heavy tail behavior, followed by event precipitation, flood seasonality, and catchment size. The pre-event moisture state in a catchment has no relevant impact on the tail heaviness even though it does influence flood magnitudes.
Evolution of chemical bonding and spin-pairing energy in ferropericlase across Its spin transition
(2022)
The evolution of chemical bonding in ferropericlase, (Mg,Fe)O, with pressure may affect the physical and chemical properties of the Earth's lower mantle. Here, we report high-pressure optical absorption spectra of single-crystalline ferropericlase ((Mg0.87Fe0.13)O) up to 135 GPa. Combined with a re-evaluation of published partial fluorescence yield X-ray absorption spectroscopy data, we show that the covalency of the Fe-O bond increases with pressure, but the iron spin transition at 57-76.5 GPa reverses this trend. The qualitative crossover in chemical bonding suggests that the spin-pairing transition weakens the Fe-O bond in ferropericlase. We find, that the spin transition in ferropericlase is caused by both the increase of the ligand field-splitting energy and the decrease in the spin-pairing energy of high-spin Fe2+.
Plate tectonics describes the movement of rigid plates at the surface of the Earth as well as their complex deformation at three types of plate boundaries: 1) divergent boundaries such as rift zones and mid-ocean ridges, 2) strike-slip boundaries where plates grind past each other, such as the San Andreas Fault, and 3) convergent boundaries that form large mountain ranges like the Andes. The generally narrow deformation zones that bound the plates exhibit complex strain patterns that evolve through time. During this evolution, plate boundary deformation is driven by tectonic forces arising from Earth’s deep interior and from within the lithosphere, but also by surface processes, which erode topographic highs and deposit the resulting sediment into regions of low elevation. Through the combination of these factors, the surface of the Earth evolves in a highly dynamic way with several feedback mechanisms. At divergent boundaries, for example, tensional stresses thin the lithosphere, forcing uplift and subsequent erosion of rift flanks, which creates a sediment source. Meanwhile, the rift center subsides and becomes a topographic low where sediments accumulate. This mass transfer from foot- to hanging wall plays an important role during rifting, as it prolongs the activity of individual normal faults. When rifting continues, continents are eventually split apart, exhuming Earth’s mantle and creating new oceanic crust. Because of the complex interplay between deep tectonic forces that shape plate boundaries and mass redistribution at the Earth’s surface, it is vital to understand feedbacks between the two domains and how they shape our planet.
In this study I aim to provide insight on two primary questions: 1) How do divergent and strike-slip plate boundaries evolve? 2) How is this evolution, on a large temporal scale and a smaller structural scale, affected by the alteration of the surface through erosion and deposition? This is done in three chapters that examine the evolution of divergent and strike-slip plate boundaries using numerical models. Chapter 2 takes a detailed look at the evolution of rift systems using two-dimensional models. Specifically, I extract faults from a range of rift models and correlate them through time to examine how fault networks evolve in space and time. By implementing a two-way coupling between the geodynamic code ASPECT and landscape evolution code FastScape, I investigate how the fault network and rift evolution are influenced by the system’s erosional efficiency, which represents many factors like lithology or climate. In Chapter 3, I examine rift evolution from a three-dimensional perspective. In this chapter I study linkage modes for offset rifts to determine when fast-rotating plate-boundary structures known as continental microplates form. Chapter 4 uses the two-way numerical coupling between tectonics and landscape evolution to investigate how a strike-slip boundary responds to large sediment loads, and whether this is sufficient to form an entirely new type of flexural strike-slip basin.
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 (SP) 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: (a) distributed deformation and coalescence, (b) fault system growth, (c) fault system decline and basinward localization, (d) rift migration, and (e) 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 SP 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.
Sea level rise and coastal erosion have inundated large areas of Arctic permafrost. Submergence by warm and saline waters increases the rate of inundated permafrost thaw compared to sub-aerial thawing on land. Studying the contact between the unfrozen and frozen sediments below the seabed, also known as the ice-bearing permafrost table (IBPT), provides valuable information to understand the evolution of sub-aquatic permafrost, which is key to improving and understanding coastal erosion prediction models and potential greenhouse gas emissions. In this study, we use data from 2D electrical resistivity tomography (ERT) collected in the nearshore coastal zone of two Arctic regions that differ in their environmental conditions (e.g., seawater depth and resistivity) to image and study the subsea permafrost. The inversion of 2D ERT data sets is commonly performed using deterministic approaches that favor smoothed solutions, which are typically interpreted using a user-specified resistivity threshold to identify the IBPT position. In contrast, to target the IBPT position directly during inversion, we use a layer-based model parameterization and a global optimization approach to invert our ERT data. This approach results in ensembles of layered 2D model solutions, which we use to identify the IBPT and estimate the resistivity of the unfrozen and frozen sediments, including estimates of uncertainties. Additionally, we globally invert 1D synthetic resistivity data and perform sensitivity analyses to study, in a simpler way, the correlations and influences of our model parameters. The set of methods provided in this study may help to further exploit ERT data collected in such permafrost environments as well as for the design of future field experiments.
Sea level rise and coastal erosion have inundated large areas of Arctic permafrost. Submergence by warm and saline waters increases the rate of inundated permafrost thaw compared to sub-aerial thawing on land. Studying the contact between the unfrozen and frozen sediments below the seabed, also known as the ice-bearing permafrost table (IBPT), provides valuable information to understand the evolution of sub-aquatic permafrost, which is key to improving and understanding coastal erosion prediction models and potential greenhouse gas emissions. In this study, we use data from 2D electrical resistivity tomography (ERT) collected in the nearshore coastal zone of two Arctic regions that differ in their environmental conditions (e.g., seawater depth and resistivity) to image and study the subsea permafrost. The inversion of 2D ERT data sets is commonly performed using deterministic approaches that favor smoothed solutions, which are typically interpreted using a user-specified resistivity threshold to identify the IBPT position. In contrast, to target the IBPT position directly during inversion, we use a layer-based model parameterization and a global optimization approach to invert our ERT data. This approach results in ensembles of layered 2D model solutions, which we use to identify the IBPT and estimate the resistivity of the unfrozen and frozen sediments, including estimates of uncertainties. Additionally, we globally invert 1D synthetic resistivity data and perform sensitivity analyses to study, in a simpler way, the correlations and influences of our model parameters. The set of methods provided in this study may help to further exploit ERT data collected in such permafrost environments as well as for the design of future field experiments.
Salt deposits offer a variety of usage types. These include the mining of rock salt and potash salt as important raw materials, the storage of energy in man-made underground caverns, and the disposal of hazardous substances in former mines. The most serious risk with any of these usage types comes from the contact with groundwater or surface water. It causes an uncontrolled dissolution of salt rock, which in the worst case can result in the flooding or collapse of underground facilities. Especially along potash seams, cavernous structures can spread quickly, because potash salts show a much higher solubility than rock salt. However, as their chemical behavior is quite complex, previous models do not account for these highly soluble interlayers. Therefore, the objective of the present thesis is to describe the evolution of cavernous structures along potash seams in space and time in order to improve hazard mitigation during the utilization of salt deposits.
The formation of cavernous structures represents an interplay of chemical and hydraulic processes. Hence, the first step is to systematically investigate the dissolution and precipitation reactions that occur when water and potash salt come into contact. For this purpose, a geochemical reaction model is used. The results show that the minerals are only partially dissolved, resulting in a porous sponge like structure. With the saturation of the solution increasing, various secondary minerals are formed, whose number and type depend on the original rock composition. Field data confirm a correlation between the degree of saturation and the distance from the center of the cavern, where solution is entering. Subsequently, the reaction model is coupled with a flow and transport code and supplemented by a novel approach called ‘interchange’. The latter enables the exchange of solution and rock between areas of different porosity and mineralogy, and thus ultimately the growth of the cavernous structure. By means of several scenario analyses, cavern shape, growth rate and mineralogy are systematically investigated, taking also heterogeneous potash seams into account. The results show that basically four different cases can be distinguished, with mixed forms being a frequent occurrence in nature. The classification scheme is based on the dimensionless numbers Péclet and Damköhler, and allows for a first assessment of the hazard potential. In future, the model can be applied to any field case, using measurement data for calibration.
The presented research work provides a reactive transport model that is able to spatially and temporally characterize the propagation of cavernous structures along potash seams for the first time. Furthermore, it allows to determine thickness and composition of transition zones between cavern center and unaffected salt rock. The latter is particularly important in potash mining, so that natural cavernous structures can be located at an early stage and the risk of mine flooding can thus be reduced. The models may also contribute to an improved hazard prevention in the construction of storage caverns and the disposal of hazardous waste in salt deposits. Predictions regarding the characteristics and evolution of cavernous structures enable a better assessment of potential hazards, such as integrity or stability loss, as well as of suitable mitigation measures.
The Permo-Triassic period marks the time interval between Hercynian (Variscan) orogenic events in the Tien Shan and the North Pamir, and the Cimmerian accretion of the Gondwana-derived Central and South Pamir to the southern margin of the Paleo-Asian continent. A well-preserved Permo-Triassic volcano-sedimentary sequence from the Chinese North Pamir yields important information on the geodynamic evolution of Asia’s pre-Cimmerian southern margin. The oldest volcanic rocks from that section are dated to the late Guadalupian epoch by a rhyolite and a dacitic dike that gave zircon U-Pb ages of ~260 Ma. Permian volcanism was largely pyroclastic and mafic to intermediate. Upsection, a massive ignimbritic crystal tuff in the Chinese Qimgan valley was dated to 244.1 +/- 1.1 Ma, a similar unit in the nearby Gez valley to 245 +/- 11 Ma, and an associated rhyolite to 233.4 +/- 1.1 Ma. Deposition of the locally ~200 m thick crystal tuff unit follows an unconformity and marks the onset of intense, mainly mafic to intermediate, calc-alkaline magmatic activity. Triassic volcanic activity in the North Pamir was coeval with the major phase of Cimmerian intrusive activity in the Karakul-Mazar arc-accretionary complex to the south, caused by northward subduction of the Paleo-Tethys. It also coincided with the emplacement of basanitic and carbonatitic dikes and a thermal event in the South Tien Shan, to the north of our study area. Evidence for arc-related magmatic activity in a back-arc position provides strong arguments for back-arc extension or transtension and basin formation. This puts the Qimgan succession in line with a more than 1000 km long realm of extensional Triassic back-arc basins known from the North Pamir in the Kyrgyz Altyn Darya valley (Myntekin formation), the North Pamir of Tajikistan and Afghanistan, and the Afghan Hindukush (Doab formation) and further west from the Paropamisus and Kopet Dag (Aghdarband, NE Iran).
The Permo-Triassic period marks the time interval between Hercynian (Variscan) orogenic events in the Tien Shan and the North Pamir, and the Cimmerian accretion of the Gondwana-derived Central and South Pamir to the southern margin of the Paleo-Asian continent. A well-preserved Permo-Triassic volcano-sedimentary sequence from the Chinese North Pamir yields important information on the geodynamic evolution of Asia’s pre-Cimmerian southern margin. The oldest volcanic rocks from that section are dated to the late Guadalupian epoch by a rhyolite and a dacitic dike that gave zircon U-Pb ages of ~260 Ma. Permian volcanism was largely pyroclastic and mafic to intermediate. Upsection, a massive ignimbritic crystal tuff in the Chinese Qimgan valley was dated to 244.1 +/- 1.1 Ma, a similar unit in the nearby Gez valley to 245 +/- 11 Ma, and an associated rhyolite to 233.4 +/- 1.1 Ma. Deposition of the locally ~200 m thick crystal tuff unit follows an unconformity and marks the onset of intense, mainly mafic to intermediate, calc-alkaline magmatic activity. Triassic volcanic activity in the North Pamir was coeval with the major phase of Cimmerian intrusive activity in the Karakul-Mazar arc-accretionary complex to the south, caused by northward subduction of the Paleo-Tethys. It also coincided with the emplacement of basanitic and carbonatitic dikes and a thermal event in the South Tien Shan, to the north of our study area. Evidence for arc-related magmatic activity in a back-arc position provides strong arguments for back-arc extension or transtension and basin formation. This puts the Qimgan succession in line with a more than 1000 km long realm of extensional Triassic back-arc basins known from the North Pamir in the Kyrgyz Altyn Darya valley (Myntekin formation), the North Pamir of Tajikistan and Afghanistan, and the Afghan Hindukush (Doab formation) and further west from the Paropamisus and Kopet Dag (Aghdarband, NE Iran).
The Gofa Province and the Chew Bahir Basin of southern Ethiopia constitute tectonically active regions, where the Southern Main Ethiopian Rift converges with the Northern Kenya Rift through a wide zone of extensional deformation with several north to northeast-trending, left-stepping en-e & PRIME;chelon basins. This sector of the Southern Main Ethiopian Rift is characterized by a semi-arid climate and a largely uniform lithology, and thus provides ideal conditions for studying the different parameters that define the tectonic and geomorphic features of this complex kinematic transfer zone. In this study, the degree of tectonic activity, spatiotemporal variations in extension, and the nature of kinematic linkage between different fault systems of the transfer zone are constrained by detailed quantitative geomorphic analysis of river catchments and focused field work. We analyzed fluvial and landscape morphometric characteristics in combination with structural, seismicity, and climatic data to better evaluate the tectono-geomorphic history of this transfer zone. Our data reveal significant north-south variations in the degree of extension from the Sawula Basin in the north (mature) to the Chew Bahir Basin in the south (juvenile). First, normalized channel-steepness indices and the spatial arrangement of knickpoints in footwall-draining streams suggest a gradual, southward shift in extensional deformation and recent tectonic activity. Second, based on 1-k(m) radius local relief and mean-hillslope maximum values that are consistent with ksn anomalies, we confirm strain localization within zones of fault interaction. Third, morphometric indices such as hypsometry, basin asymmetry factor, and valley floor width to valley height ratio also indicate a north to south gradient in tectonic activity, highlighting the importance of such a wide transfer zone with diffuse extension linking different rift segments during the break-up of continental crust.
We construct and examine the prototype of a deep learning-based ground-motion model (GMM) that is both fully data driven and nonergodic. We formulate ground-motion modeling as an image processing task, in which a specific type of neural network, the U-Net, relates continuous, horizontal maps of earthquake predictive parameters to sparse observations of a ground-motion intensity measure (IM). The processing of map-shaped data allows the natural incorporation of absolute earthquake source and observation site coordinates, and is, therefore, well suited to include site-, source-, and path-specific amplification effects in a nonergodic GMM. Data-driven interpolation of the IM between observation points is an inherent feature of the U-Net and requires no a priori assumptions. We evaluate our model using both a synthetic dataset and a subset of observations from the KiK-net strong motion network in the Kanto basin in Japan. We find that the U-Net model is capable of learning the magnitude???distance scaling, as well as site-, source-, and path-specific amplification effects from a strong motion dataset. The interpolation scheme is evaluated using a fivefold cross validation and is found to provide on average unbiased predictions. The magnitude???distance scaling as well as the site amplification of response spectral acceleration at a period of 1 s obtained for the Kanto basin are comparable to previous regional studies.
The response of rapidly compressed highly oriented pyrolytic graphite (HOPG) normal to its basal plane was investigated at a pressure of & SIM;80 GPa. Ultrafast x-ray diffraction using & SIM;100 fs pulses at the Materials Under Extreme Conditions sector of the Linac Coherent Light Source was used to probe the changes in crystal structure resulting from picosecond timescale compression at laser drive energies ranging from 2.5 to 250 mJ. A phase transformation from HOPG to a highly textured hexagonal diamond structure is observed at the highest energy, followed by relaxation to a still highly oriented, but distorted graphite structure following release. We observe the formation of a highly oriented lonsdaleite within 20 ps, subsequent to compression. This suggests that a diffusionless martensitic mechanism may play a fundamental role in phase transition, as speculated in an early work on this system, and more recent static studies of diamonds formed in impact events. Published by AIP Publishing.
River-valley morphology preserves information on tectonic and climatic conditions that shape landscapes. Observations suggest that river discharge and valley-wall lithology are the main controls on valley width. Yet, current models based on these observations fail to explain the full range of cross-sectional valley shapes in nature, suggesting hitherto unquantified controls on valley width. In particular, current models cannot explain the existence of paired terrace sequences that form under cyclic climate forcing. Paired river terraces are staircases of abandoned floodplains on both valley sides, and hence preserve past valley widths. Their formation requires alternating phases of predominantly river incision and predominantly lateral planation, plus progressive valley narrowing. While cyclic Quaternary climate changes can explain shifts between incision and lateral erosion, the driving mechanism of valley narrowing is unknown. Here, we extract valley geometries from climatically formed, alluvial river-terrace sequences and show that across our dataset, the total cumulative terrace height (here: total valley height) explains 90%–99% of the variance in valley width at the terrace sites. This finding suggests that valley height, or a parameter that scales linearly with valley height, controls valley width in addition to river discharge and lithology. To explain this valley-width-height relationship, we reformulate existing valley-width models and suggest that, when adjusting to new boundary conditions, alluvial valleys evolve to a width at which sediment removal from valley walls matches lateral sediment supply from hillslope erosion. Such a hillslope-channel coupling is not captured in current valley-evolution models. Our model can explain the existence of paired terrace sequences under cyclic climate forcing and relates valley width to measurable field parameters. Therefore, it facilitates the reconstruction of past climatic and tectonic conditions from valley topography.
River-valley morphology preserves information on tectonic and climatic conditions that shape landscapes. Observations suggest that river discharge and valley-wall lithology are the main controls on valley width. Yet, current models based on these observations fail to explain the full range of cross-sectional valley shapes in nature, suggesting hitherto unquantified controls on valley width. In particular, current models cannot explain the existence of paired terrace sequences that form under cyclic climate forcing. Paired river terraces are staircases of abandoned floodplains on both valley sides, and hence preserve past valley widths. Their formation requires alternating phases of predominantly river incision and predominantly lateral planation, plus progressive valley narrowing. While cyclic Quaternary climate changes can explain shifts between incision and lateral erosion, the driving mechanism of valley narrowing is unknown. Here, we extract valley geometries from climatically formed, alluvial river-terrace sequences and show that across our dataset, the total cumulative terrace height (here: total valley height) explains 90%–99% of the variance in valley width at the terrace sites. This finding suggests that valley height, or a parameter that scales linearly with valley height, controls valley width in addition to river discharge and lithology. To explain this valley-width-height relationship, we reformulate existing valley-width models and suggest that, when adjusting to new boundary conditions, alluvial valleys evolve to a width at which sediment removal from valley walls matches lateral sediment supply from hillslope erosion. Such a hillslope-channel coupling is not captured in current valley-evolution models. Our model can explain the existence of paired terrace sequences under cyclic climate forcing and relates valley width to measurable field parameters. Therefore, it facilitates the reconstruction of past climatic and tectonic conditions from valley topography.
A detailed analysis of horizontal and vertical particulate matter (PM) fluxes during wind erosion has been done, based on measurements of PM smaller than 10, 2.5, and 1.0 mu mm, at windward and leeward positions on a measuring field. The three fractions of PM measurement are differently influenced by the increasing wind and shear velocities of the wind. The measured concentrations of the coarser fractions of the fine dust, PM10, and PM2.5, increase with wind and shear velocity, whereas the PM1.0 concentrations show no clear correlation to the shear velocity. The share of PM2.5 on PM10 depends on the measurement height and wind speed and varies between 4 and 12 m/s at the 1 m height ranging from 25% to 7% (average 10%), and at the 4 m height from 39% to 23% (average 30%). Although general relationships between wind speed, PM concentration, and horizontal and vertical fluxes could be found, the contribution of the measuring field was very low, as balances of incoming and outgoing fluxes show. Consequently, the measured PM concentrations are determined from a variety of sources, such as traffic on unpaved roads, cattle drives, tillage operations, and wind erosion, and thus, represent all components of land use and landscape structure in the near and far surroundings of the measuring field. The current results may reflect factors from the landscape scale rather than the influence of field-related variables. The measuring devices used to monitor PM concentrations showed differences of up to 20%, which led to considerable deviations when determining total balances. Differences up to 67% between the calculated fluxes prove the necessity of a previous calibration of the devices used. (c) 2022 International Research and Training Centre on Erosion and Sedimentation/the World Association for Sedimentation and Erosion Research.
Hydraulic-driven fractures play a key role in subsurface energy technologies across several scales. By injecting fluid at high hydraulic pressure into rock with intrinsic low permeability, in-situ stress field and fracture development pattern can be characterised as well as rock permeability can be enhanced. Hydraulic fracturing is a commercial standard procedure for enhanced oil and gas production of rock reservoirs with low permeability in petroleum industry. However, in EGS utilization, a major geological concern is the unsolicited generation of earthquakes due to fault reactivation, referred to as induced seismicity, with a magnitude large enough to be felt on the surface or to damage facilities and buildings. Furthermore, reliable interpretation of hydraulic fracturing tests for stress measurement is a great challenge for the energy technologies. Therefore, in this cumulative doctoral thesis the following research questions are investigated. (1): How do hydraulic fractures grow in hard rock at various scales?; (2): Which parameters control hydraulic fracturing and hydro-mechanical coupling?; and (3): How can hydraulic fracturing in hard rock be modelled?
In the laboratory scale study, several laboratory hydraulic fracturing experiments are investigated numerically using Irazu2D that were performed on intact cubic Pocheon granite samples from South Korea applying different injection protocols. The goal of the laboratory experiments is to test the concept of cyclic soft stimulation which may enable sustainable permeability enhancement (Publication 1).
In the borehole scale study, hydraulic fracturing tests are reported that were performed in boreholes located in central Hungary to determine the in-situ stress for a geological site investigation. At depth of about 540 m, the recorded pressure versus time curves in mica schist with low dip angle foliation show atypical evolution. In order to provide explanation for this observation, a series of discrete element computations using Particle Flow Code 2D are performed (Publication 2).
In the reservoir scale study, the hydro-mechanical behaviour of fractured crystalline rock due to one of the five hydraulic stimulations at the Pohang Enhanced Geothermal site in South Korea is studied. Fluid pressure perturbation at faults of several hundred-meter lengths during hydraulic stimulation is simulated using FracMan (Publication 3).
The doctoral research shows that the resulting hydraulic fracturing geometry will depend “locally”, i.e. at the length scale of representative elementary volume (REV) and below that (sub-REV), on the geometry and strength of natural fractures, and “globally”, i.e. at super-REV domain volume, on far-field stresses. Regarding hydro-mechanical coupling, it is suggested to define separate coupling relationship for intact rock mass and natural fractures. Furthermore, the relative importance of parameters affecting the magnitude of formation breakdown pressure, a parameter characterising hydro-mechanical coupling, is defined. It can be also concluded that there is a clear gap between the capacity of the simulation software and the complexity of the studied problems. Therefore, the computational time of the simulation of complex hydraulic fracture geometries must be reduced while maintaining high fidelity simulation results. This can be achieved either by extending the computational resources via parallelization techniques or using time scaling techniques. The ongoing development of used numerical models focuses on tackling these methodological challenges.
The new in situ geodynamic laboratory established in the framework of the ICDP Eger project aims to develop the most modern, comprehensive, multiparameter laboratory at depth for studying earthquake swarms, crustal fluid flow, mantle-derived CO2 and helium degassing, and processes of the deep biosphere. In order to reach a new level of high-frequency, near-source and multiparameter observation of earthquake swarms and related phenomena, such a laboratory comprises a set of shallow boreholes with high-frequency 3-D seismic arrays as well as modern continuous real-time fluid monitoring at depth and the study of the deep biosphere.
This laboratory is located in the western part of the Eger Rift at the border of the Czech Republic and Germany (in the West Bohemia-Vogtland geodynamic region) and comprises a set of five boreholes around the seismoactive zone. To date, all monitoring boreholes have been drilled. This includes the seismic monitoring boreholes S1, S2 and S3 in the crystalline units north and east of the major Novy Kostel seismogenic zone, borehole F3 in the Hartousov mofette field and borehole S4 in the newly discovered Bazina maar near Liba. Supplementary borehole P1 is being prepared in the Neualbenreuth maar for paleoclimate and biological research. At each of these sites, a borehole broadband seismometer will be installed, and sites S1, S2 and S3 will also host a 3-D seismic array composed of a vertical geophone chain and surface seismic array. Seismic instrumenting has been completed in the S1 borehole and is in preparation in the remaining four monitoring boreholes. The continuous fluid monitoring site of Hartousov includes three boreholes, F1, F2 and F3, and a pilot monitoring phase is underway. The laboratory also enables one to analyze microbial activity at CO2 mofettes and maar structures in the context of changes in habitats. The drillings into the maar volcanoes contribute to a better understanding of the Quaternary paleoclimate and volcanic activity.
The new in situ geodynamic laboratory established in the framework of the ICDP Eger project aims to develop the most modern, comprehensive, multiparameter laboratory at depth for studying earthquake swarms, crustal fluid flow, mantle-derived CO2 and helium degassing, and processes of the deep biosphere. In order to reach a new level of high-frequency, near-source and multiparameter observation of earthquake swarms and related phenomena, such a laboratory comprises a set of shallow boreholes with high-frequency 3-D seismic arrays as well as modern continuous real-time fluid monitoring at depth and the study of the deep biosphere.
This laboratory is located in the western part of the Eger Rift at the border of the Czech Republic and Germany (in the West Bohemia–Vogtland geodynamic region) and comprises a set of five boreholes around the seismoactive zone. To date, all monitoring boreholes have been drilled. This includes the seismic monitoring boreholes S1, S2 and S3 in the crystalline units north and east of the major Nový Kostel seismogenic zone, borehole F3 in the Hartoušov mofette field and borehole S4 in the newly discovered Bažina maar near Libá. Supplementary borehole P1 is being prepared in the Neualbenreuth maar for paleoclimate and biological research. At each of these sites, a borehole broadband seismometer will be installed, and sites S1, S2 and S3 will also host a 3-D seismic array composed of a vertical geophone chain and surface seismic array. Seismic instrumenting has been completed in the S1 borehole and is in preparation in the remaining four monitoring boreholes. The continuous fluid monitoring site of Hartoušov includes three boreholes, F1, F2 and F3, and a pilot monitoring phase is underway. The laboratory also enables one to analyze microbial activity at CO2 mofettes and maar structures in the context of changes in habitats. The drillings into the maar volcanoes contribute to a better understanding of the Quaternary paleoclimate and volcanic activity.
Enhanced geothermal systems (EGS) are considered a cornerstone of future sustainable energy production. In such systems, high-pressure fluid injections break the rock to provide pathways for water to circulate in and heat up. This approach inherently induces small seismic events that, in rare cases, are felt or can even cause damage. Controlling and reducing the seismic impact of EGS is crucial for a broader public acceptance. To evaluate the applicability of hydraulic fracturing (HF) in EGS and to improve the understanding of fracturing processes and the hydromechanical relation to induced seismicity, six in-situ, meter-scale HF experiments with different injection schemes were performed under controlled conditions in crystalline rock in a depth of 410 m at the Äspö Hard Rock Laboratory (Sweden).
I developed a semi-automated, full-waveform-based detection, classification, and location workflow to extract and characterize the acoustic emission (AE) activity from the continuous recordings of 11 piezoelectric AE sensors. Based on the resulting catalog of 20,000 AEs, with rupture sizes of cm to dm, I mapped and characterized the fracture growth in great detail. The injection using a novel cyclic injection scheme (HF3) had a lower seismic impact than the conventional injections. HF3 induced fewer AEs with a reduced maximum magnitude and significantly larger b-values, implying a decreased number of large events relative to the number of small ones. Furthermore, HF3 showed an increased fracture complexity with multiple fractures or a fracture network. In contrast, the conventional injections developed single, planar fracture zones (Publication 1).
An independent, complementary approach based on a comparison of modeled and observed tilt exploits transient long-period signals recorded at the horizontal components of two broad-band seismometers a few tens of meters apart from the injections. It validated the efficient creation of hydraulic fractures and verified the AE-based fracture geometries. The innovative joint analysis of AEs and tilt signals revealed different phases of the fracturing process, including the (re-)opening, growth, and aftergrowth of fractures, and provided evidence for the reactivation of a preexisting fault in one of the experiments (Publication 2). A newly developed network-based waveform-similarity analysis applied to the massive AE activity supports the latter finding.
To validate whether the reduction of the seismic impact as observed for the cyclic injection schemes during the Äspö mine-scale experiments is transferable to other scales, I additionally calculated energy budgets for injection experiments from previously conducted laboratory tests and from a field application. Across all three scales, the cyclic injections reduce the seismic impact, as depicted by smaller maximum magnitudes, larger b-values, and decreased injection efficiencies (Publication 3).
Soil bacteria play a fundamental role in pedogenesis. However, knowledge about both the impact of climate and slope aspects on microbial communities and the consequences of these items in pedogenesis is lacking. Therefore, soil-bacterial communities from four sites and two different aspects along the climate gradient of the Chilean Coastal Cordillera were investigated. Using a combination of microbiological and physicochemical methods, soils that developed in arid, semi-arid, mediterranean, and humid climates were analyzed. Proteobacteria, Acidobacteria, Chloroflexi, Verrucomicrobia, and Planctomycetes were found to increase in abundance from arid to humid climates, while Actinobacteria and Gemmatimonadetes decreased along the transect. Bacterial-community structure varied with climate and aspect and was influenced by pH, bulk density, plant-available phosphorus, clay, and total organic-matter content. Higher bacterial specialization was found in arid and humid climates and on the south-facing slope and was likely promoted by stable microclimatic conditions. The presence of specialists was associated with ecosystem-functional traits, which shifted from pioneers that accumulated organic matter in arid climates to organic decomposers in humid climates. These findings provide new perspectives on how climate and slope aspects influence the composition and functional capabilities of bacteria, with most of these capabilities being involved in pedogenetic processes.
Phase transitions in molecular crystals are often determined by intermolecular interactions. The cage complex of [Co(C12H30N8)](3+) . 3 NO3- is reported to undergo a disorder-order phase transition at T-c1 approximate to 133 K upon cooling. Temperature-dependent neutron and synchrotron diffraction experiments revealed satellite reflections in addition to main reflections in the diffraction patterns below T-c1. The modulation wave vector varies as function of temperature and locks in at T-c3 approximate to 98 K. Here, we demonstrate that the crystal symmetry lowers from hexagonal to monoclinic in the incommensurately modulated phases in T-c1<T<T-c3. Distinctive levels of competitions: trade-off between longer N-H...O and shorter C-H...O hydrogen bonds; steric constraints to dense C-H...O bonds give rise to pronounced modulation of the basic structure. Severely frustrated crystal packing in the incommensurate phase is precursor to optimal balance of intermolecular interactions in the lock-in phase.
Despite its high-seismogenic potential, the details of the seismogenic processes of Zagros Simply Folded Belt (SFB) remains debated. Three large earthquakes (M-w 7.3, 5.9 and 6.3) struck in the Lurestan arc of the Zagros SFB in 2017 and 2018. The sequence was recorded by seismic stations at regional, and teleseismic distances. Coseismic surface displacements, measured by Sentinel-1A/B satellites, provide additional data and a unique opportunity to study these earthquakes in detail. Here, we complement previous studies of the coseismic slip distribution of the 12 November 2017 M-w 7.3 Ezgeleh earthquake by a detailed analysis of its aftershocks, and we analysed the rupture process of the two interrelated earthquakes (25 August 2018 M-w 5.9 Tazehabad and the 25 November 2018 M-w 6.3 Sarpol-e Zahab earthquakes). We model the surface displacements obtained from Interferometric Synthetic Aperture Radar (InSAR) measurements and seismic records. We conduct non-linear probabilistic optimizations based on joint InSAR and seismic data to obtain finite-fault rupture of these earthquakes. The Lurestan arc earthquakes were followed by a sustained aftershock activity, with 133 aftershocks exceeding M-n 4.0 until 30 December 2019. We rely on the permanent seismic networks of Iran and Iraq to relocate similar to 700 M-n 3 + events and estimate moment tensor solutions for 85 aftershocks down to M-w 4.0. The 2017 Ezgeleh earthquake has been considered to activate a low-angle (similar to 17 degrees) dextral-thrust fault at the depth of 10-20 km. However, most of its aftershocks have shallow centroid depths (8-12 km). The joint interpretation of finite source models, moment tensor and hypocentral location indicate that the 2018 Tazehabad and Sarpol-e Zahab earthquakes ruptured different strike-slip structures, providing evidence for the activation of the sinistral and dextral strike-slip faults, respectively. The deformation in the Lurestan arc is seismically accommodated by a complex fault system involving both thrust and strike-slip faults. Knowledge about the deformation characteristics is important for the understanding of crustal shortening, faulting and hazard and risk assessment in this region.
Understanding the key factors influencing the water quality of large river systems forms an important basis for the assessment and protection of cross-regional ecosystems and the implementation of adapted water management concepts. However, identifying these factors requires in-depth comprehension of the unique environmental systems, which can only be achieved by detailed water quality monitoring.
Within the scope of the joint science and sports event "Elbschwimmstaffel" (swimming relay on the river Elbe) in June/July 2017 organized by the German Ministry of Education and Research, water quality data were acquired along a 550 km long stretch of the Elbe River in Germany. During the survey, eight physiochemical water quality parameters were recorded in high spatial and temporal resolution with the BIOFISH multisensor system. Multivariate statistical methods were applied to identify and delineate processes influencing the water quality.
The BIOFISH dataset revealed that phytoplankton activity has a major impact on the water quality of the Elbe River in the summer months. The results suggest that phytoplankton biomass constitutes a substantial proportion of the suspended particles and that photosynthetic activity of phytoplankton is closely related to significant temporal changes in pH and oxygen saturation.
An evaluation of the BIOFISH data based on the combination of statistical analysis with weather and discharge data shows that the hydrological and meteorological history of the sampled water body was the main driver of phytoplankton dynamics. This study demonstrates the capacity of longitudinal river surveys with the BIOFISH or similar systems for water quality assessment, the identification of pollution sources and their utilization for online in situ monitoring of rivers.
Rainfall-intense summer monsoon seasons on the Indian subcontinent that are exceeding long-term averages cause widespread floods and landslides.
Here we show that the latest generation of coupled climate models robustly project an intensification of very rainfall-intense seasons (June-September).
Under the shared socioeconomic pathway SSP5-8.5, very wet monsoon seasons as observed in only 5 years in the period 1965-2015 are projected to occur 8 times more often in 2050-2100 in the multi-model average.
Under SSP2-4.5, these seasons become only a factor of 6 times more frequent, showing that even modest efforts to mitigate climate change can have a strong impact on the frequency of very strong rainfall seasons.
Besides, we find that the increasing risk of extreme seasonal rainfall is accompanied by a shift from days with light rainfall to days with moderate or heavy rainfall. Additionally, the number of wet days is projected to increase.
The geometry of carbonate platforms reflects the interaction of several factors. However, the impact of carbonate-producing organisms has been poorly investigated so far. This study applies stratigraphic forward modelling (SFM) and sensitivity analysis to examine, referenced to the Miocene Llucmajor Platform, the effect of changes of dominant biotic production in the oligophotic and euphotic zones on platform geometry. Our results show that the complex interplay of carbonate production rates, bathymetry and variations in accommodation space control the platform geometry. The main driver of progradation is the oligophotic production of rhodalgal sediments during the lowstands. This study demonstrates that platform geometry and internal architecture varies significantly according to the interaction of the predominant carbonate-producing biotas. The input parameters for this study are based on well-understood Miocene carbonate biotas with characteristic euphotic, oligophotic and photo-independent carbonate production in which it is crucial that each carbonate-producing class is modelled explicitly within the simulation run and not averaged with a single carbonate production-depth profile. This is important in subsurface exploration studies based on stratigraphic forward models where the overall platform geometry may be approximated through calibration runs, and constrained by seismic surveys and wellbores. However, the internal architecture is likely to be oversimplified without an in-depth understanding of the target carbonate system and a transfer to forward modelling parameters.
ABSTRACT: Structural evolution of cesium triiodide at high pressures has been revealed by synchrotron single-crystal X-ray diffraction. Cesium triiodide undergoes a first-order phase transition above 1.24(3) GPa from an orthorhombic to a trigonal system. This transition is coupled with severe reorganization of the polyiodide network from a layered to three-dimensional architecture. Quantum chemical calculations show that even though the two polymorphic phases are nearly isoenergetic under ambient conditions, the PV term is decisive in stabilizing the trigonal polymorph above the transition point. Phonon calculations using a non-local correlation functional that accounts for dispersion interactions confirm that this polymorph is dynamically unstable under ambient conditions. The high-pressure behavior of crystalline CsI3 can be correlated with other alkali metal trihalides, which undergo a similar sequence of structural changes upon load.
On their way from inland to the ocean, flowing water bodies, their constituents and their biotic communities are ex-posed to complex transport and transformation processes. However, detailed process knowledge as revealed by La-grangian measurements adjusted to travel time is rare in large rivers, in particular at hydrological extremes. To fill this gap, we investigated autotrophic processes, heterotrophic carbon utilization, and micropollutant concentrations applying a Lagrangian sampling design in a 600 km section of the River Elbe (Germany) at historically low discharge. Under base flow conditions, we expect the maximum intensity of instream processes and of point source impacts. Phy-toplankton biomass and photosynthesis increased from upstream to downstream sites but maximum chlorophyll con-centration was lower than at mean discharge. Concentrations of dissolved macronutrients decreased to almost complete phosphate depletion and low nitrate values. The longitudinal increase of bacterial abundance and production was less pronounced than in wetter years and bacterial community composition changed downstream. Molecular analyses revealed a longitudinal increase of many DOM components due to microbial production, whereas saturated lipid-like DOM, unsaturated aromatics and polyphenols, and some CHOS surfactants declined. In decomposition exper-iments, DOM components with high O/C ratios and high masses decreased whereas those with low O/C ratios, low masses, and high nitrogen content increased at all sites. Radiocarbon age analyses showed that DOC was relatively old (890-1870 years B.P.), whereas the mineralized fraction was much younger suggesting predominant oxidation of algal lysis products and exudates particularly at downstream sites. Micropollutants determining toxicity for algae (terbuthylazine, terbutryn, isoproturon and lenacil), hexachlorocyclohexanes and DDTs showed higher concentrations from the middle towards the downstream part but calculated toxicity was not negatively correlated to phytoplankton. Overall, autotrophic and heterotrophic process rates and micropollutant concentrations increased from up-to down-stream reaches, but their magnitudes were not distinctly different to conditions at medium discharges.
Landslides
(2022)
Erosion by landslides is a common phenomenon in mountain regions around the globe, affecting all climatic zones. Landslides facilitate bedrock weathering, pedogenesis and ecological succession, being key drivers of biodiversity. Landslide chronosequences have long been used for studies of vegetation succession in initial ecosystems, but they further offer ideal model systems for studies of soil development and microbial community succession. In this review we synthesize the state of knowledge on the role of landslides in ecosystems, their influence on element cycles and interactions with biota. Further, we discuss feedback mechanisms between global warming, landslide activity and greenhouse gas emissions. In the view of increasing anthropogenic influence and climate change, soils are becoming a critical resource. Due to their ubiquity, landslide chronosequences have the potential to provide critical insights into soil development under different climates and thereby contribute to future soil restoration efforts.
Large-Scale interseismic strain mapping of the NE Tibetan Plateau from Sentinel-1 Interferometry
(2022)
The launches of the Sentinel-1 synthetic aperture radar satellites in 2014 and 2016 started a new era of high-resolution velocity and strain rate mapping for the continents. However, multiple challenges exist in tying independently processed velocity data sets to a common reference frame and producing high-resolution strain rate fields. We analyze Sentinel-1 data acquired between 2014 and 2019 over the northeast Tibetan Plateau, and develop new methods to derive east and vertical velocities with similar to 100 m resolution and similar to 1 mm/yr accuracy across an area of 440,000 km(2). By implementing a new method of combining horizontal gradients of filtered east and interpolated north velocities, we derive the first similar to 1 km resolution strain rate field for this tectonically active region. The strain rate fields show concentrated shear strain along the Haiyuan and East Kunlun Faults, and local contractional strain on fault junctions, within the Qilianshan thrusts, and around the Longyangxia Reservoir. The Laohushan-Jingtai creeping section of the Haiyuan Fault is highlighted in our data set by extremely rapid strain rates. Strain across unknown portions of the Haiyuan Fault system, including shear on the eastern extension of the Dabanshan Fault and contraction at the western flank of the Quwushan, highlight unmapped tectonic structures. In addition to the uplift across most of the lowlands, the vertical velocities also contain climatic, hydrological or anthropogenic-related deformation signals. We demonstrate the enhanced view of large-scale active tectonic processes provided by high-resolution velocities and strain rates derived from Sentinel-1 data and highlight associated wide-ranging research applications.
The origin of the First Bend of the Yangtze River is key to understanding the birth of the modern Yangtze River. Despite considerable efforts, the timing and mechanism of formation of the First Bend remain highly debated. Inverse river-profile modeling of three tributaries (Chongjiang, Lima, and Gudu) of the Jinsha River, integrated with regional tectonic and geomorphic interpretations, allows the onset of incision at the First Bend to be constrained to 28-20 Ma. The spatio-temporal coincidence of initial river incision and activity of Yulong Thrust Belt in southeastern Tibet highlights thrusting to be fundamental in reshaping the pre-existing stream network at the First Bend. These results enable us to reinterpret a change in sedimentary environment from a braided river to a swamp-like lake in the Jianchuan Basin south of the First Bend, recording the destruction of the hypothesized southwards-flowing paleo-Jinsha and Shuiluo Rivers at ~36-35 Ma by magmatism. During the late Oligoceneearly Miocene, the paleo-Shuiluo River was diverted to the north by focused rock uplift due to thrusting along the Yulong Thrust Belt, which also led to exhumation of the Jianchuan Basin. Diversion of the paleo-Shuiluo River can be explained by capture from a downstream river in the footwall of the Yulong Thrust Belt. Subsequent rapid headward erosion, that was caused by thrusting-induced drop of local base level, is recorded by upstream younging ages for the onset of incision and led to the formation of the First Bend. The combination of new ages for the onset of incision at 28-20 Ma at the First Bend and younger ages upstream indicates northwards expansion of the Jinsha River at a rate of 62 +/- 18 mm/yr. Our results suggest that the origin of the First Bend was likely triggered by thrusting at 28-20 Ma, after which the Yangtze River formed.
Frequency-domain electromagnetic (FDEM) data are commonly inverted to characterize subsurface geoelectrical properties using smoothness constraints in 1D inversion schemes assuming a layered medium.
Smoothness constraints are suitable for imaging gradual transitions of subsurface geoelectrical properties caused, for example, by varying sand, clay, or fluid content. However, such inversion approaches are limited in characterizing sharp interfaces. Alternative regularizations based on the minimum gradient support (MGS) stabilizers can, instead, be used to promote results with different levels of smoothness/sharpness selected by simply acting on the so-called focusing parameter.
The MGS regularization has been implemented for different kinds of geophysical data inversion strategies. However, concerning FDEM data, the MGS regularization has only been implemented for vertically constrained inversion (VCI) approaches but not for laterally constrained inversion (LCI) approaches.
We present a novel LCI approach for FDEM data using the MGS regularization for the vertical and lateral direction. Using synthetic and field data examples, we demonstrate that our approach can efficiently and automatically provide a set of model solutions characterized by different levels of sharpness and variable lateral consistencies.
In terms of data misfit, the obtained set of solutions contains equivalent models allowing us also to investigate the non-uniqueness of FDEM data inversion.
LegacyPollen 1.0
(2022)
Here we describe the LegacyPollen 1.0, a dataset of 2831 fossil pollen records with metadata, a harmonized taxonomy, and standardized chronologies.
A total of 1032 records originate from North America, 1075 from Europe, 488 from Asia, 150 from Latin America, 54 from Africa, and 32 from the Indo-Pacific.
The pollen data cover the late Quaternary (mostly the Holocene). The original 10 110 pollen taxa names (including variations in the notations) were harmonized to 1002 terrestrial taxa (including Cyperaceae), with woody taxa and major herbaceous taxa harmonized to genus level and other herbaceous taxa to family level.
The dataset is valuable for synthesis studies of, for example, taxa areal changes, vegetation dynamics, human impacts (e.g., deforestation), and climate change at global or continental scales.
The harmonized pollen and metadata as well as the harmonization table are available from PANGAEA (https://doi.org/10.1594/PANGAEA.929773; Herzschuh et al., 2021). R code for the harmonization is provided at Zenodo (https://doi.org/10.5281/zenodo.5910972; Herzschuh et al., 2022) so that datasets at a customized harmonization level can be easily established.
Li and B in ascending magmas: an experimental study on their mobility and isotopic fractionation
(2022)
This research study focuses on the behaviour of Li and B during magmatic ascent, and decompression-driven degassing related to volcanic systems. The main objective of this dissertation is to determine whether it is possible to use the diffusion properties of the two trace elements as a tool to trace magmatic ascent rate. With this objective, diffusion-couple and decompression experiments have been performed in order to study Li and B mobility in intra-melt conditions first, and then in an evolving system during decompression-driven degassing.
Synthetic glasses were prepared with rhyolitic composition and an initial water content of 4.2 wt%, and all the experiments were performed using an internally heated pressure vessel, in order to ensure a precise control on the experimental parameters such as temperature and pressure.
Diffusion-couple experiments were performed with a fix pressure 300 MPa. The temperature was varied in the range of 700-1250 °C with durations between 0 seconds and 24 hours. The diffusion-couple results show that Li diffusivity is very fast and starts already at very low temperature. Significant isotopic fractionation occurs due to the faster mobility of 6Li compared to 7Li. Boron diffusion is also accelerated by the presence of water, but the results of the isotopic ratios are unclear, and further investigation would be necessary to well constrain the isotopic fractionation process of boron in hydrous silicate melts. The isotopic ratios results show that boron isotopic fractionation might be affected by the speciation of boron in the silicate melt structure, as 10B and 11B tend to have tetrahedral and trigonal coordination, respectively.
Several decompression experiments were performed at 900 °C and 1000 °C, with pressures going from 300 MPa to 71-77 MPa and durations of 30 minutes, two, five and ten hours, in order to trigger water exsolution and the formation of vesicles in the sample. Textural observations and the calculation of the bubble number density confirmed that the bubble size and distribution after decompression is directly proportional to the decompression rate.
The overall SIMS results of Li and B show that the two trace elements tend to progressively decrease their concentration with decreasing decompression rates. This is explained because for longer decompression times, the diffusion of Li and B into the bubbles has more time to progress and the melt continuously loses volatiles as the bubbles expand their volumes.
For fast decompression, Li and B results show a concentration increase with a δ7Li and δ11B decrease close to the bubble interface, related to the sudden formation of the gas bubble, and the occurrence of a diffusion process in the opposite direction, from the bubble meniscus to the unaltered melt. When the bubble growth becomes dominant and Li and B start to exsolve into the gas phase, the silicate melt close to the bubble gets depleted in Li and B, because of a stronger diffusion of the trace elements into the bubble.
Our data are being applied to different models, aiming to combine the dynamics of bubble nucleation and growth with the evolution of trace elements concentration and isotopic ratios. Here, first considerations on these models will be presented, giving concluding remarks on this research study. All in all, the final remarks constitute a good starting point for further investigations. These results are a promising base to continue to study this process, and Li and B can indeed show clear dependences on decompression-related magma ascent rates in volcanic systems.
We examined the relationship between the mechanical strength of the lithosphere and the distribution of seismicity within the overriding continental plate of the southern Central Andes (SCA, 29 degrees-39 degrees S), where the oceanic Nazca Plate changes its subduction angle between 33 degrees S and 35 degrees S, from subhorizontal in the north (<5 degrees) to steep in the south (similar to 30 degrees). We computed the long-term lithospheric strength based on an existing 3D model describing variations in thickness, density, and temperature of the main geological units forming the lithosphere of the SCA and adjacent forearc and foreland regions. The comparison between our results and seismicity within the overriding plate (upper-plate seismicity) shows that most of the events occur within the modeled brittle domain of the lithosphere. The depth where the deformation mode switches from brittle frictional to thermally activated ductile creep provides a conservative lower bound to the seismogenic zone in the overriding plate of the study area. We also found that the majority of upper-plate earthquakes occurs within the realm of first-order contrasts in integrated strength (12.7-13.3 log Pam in the Andean orogen vs. 13.5-13.9 log Pam in the forearc and the foreland). Specific conditions characterize the mechanically strong northern foreland of the Andes, where seismicity is likely explained by the effects of slab steepening.
The current awareness of the high importance of urban green leads to a stronger need for tools to comprehensively represent urban green and its benefits. A common scientific approach is the development of urban ecosystem services (UES) based on remote sensing methods at the city or district level. Urban planning, however, requires fine-grained data that match local management practices. Hence, this study linked local biotope and tree mapping methods to the concept of ecosystem services. The methodology was tested in an inner-city district in SW Germany, comparing publicly accessible areas and non-accessible courtyards. The results provide area-specific [m(2)] information on the green inventory at the microscale, whereas derived stock and UES indicators form the basis for comparative analyses regarding climate adaptation and biodiversity. In the case study, there are ten times more micro-scale green spaces in private courtyards than in the public space, as well as twice as many trees. The approach transfers a scientific concept into municipal planning practice, enables the quantitative assessment of urban green at the microscale and illustrates the importance for green stock data in private areas to enhance decision support in urban development. Different aspects concerning data collection and data availability are critically discussed.
The Ice, Cloud, and Land Elevation Satellite-2 (ICESat-2) with its land and vegetation height data product (ATL08), and Global Ecosystem Dynamics Investigation (GEDI) with its terrain elevation and height metrics data product (GEDI Level 2A) missions have great potential to globally map ground and canopy heights. Canopy height is a key factor in estimating above-ground biomass and its seasonal changes; these satellite missions can also improve estimated above-ground carbon stocks. This study presents a novel Sparse Vegetation Detection Algorithm (SVDA) which uses ICESat-2 (ATL03, geolocated photons) data to map tree and vegetation heights in a sparsely vegetated savanna ecosystem. The SVDA consists of three main steps: First, noise photons are filtered using the signal confidence flag from ATL03 data and local point statistics. Second, we classify ground photons based on photon height percentiles. Third, tree and grass photons are classified based on the number of neighbors. We validated tree heights with field measurements (n = 55), finding a root-mean-square error (RMSE) of 1.82 m using SVDA, GEDI Level 2A (Geolocated Elevation and Height Metrics product): 1.33 m, and ATL08: 5.59 m. Our results indicate that the SVDA is effective in identifying canopy photons in savanna ecosystems, where ATL08 performs poorly. We further identify seasonal vegetation height changes with an emphasis on vegetation below 3 m; widespread height changes in this class from two wet-dry cycles show maximum seasonal changes of 1 m, possibly related to seasonal grass-height differences. Our study shows the difficulties of vegetation measurements in savanna ecosystems but provides the first estimates of seasonal biomass changes.
The Ice, Cloud, and Land Elevation Satellite-2 (ICESat-2) with its land and vegetation height data product (ATL08), and Global Ecosystem Dynamics Investigation (GEDI) with its terrain elevation and height metrics data product (GEDI Level 2A) missions have great potential to globally map ground and canopy heights. Canopy height is a key factor in estimating above-ground biomass and its seasonal changes; these satellite missions can also improve estimated above-ground carbon stocks. This study presents a novel Sparse Vegetation Detection Algorithm (SVDA) which uses ICESat-2 (ATL03, geolocated photons) data to map tree and vegetation heights in a sparsely vegetated savanna ecosystem. The SVDA consists of three main steps: First, noise photons are filtered using the signal confidence flag from ATL03 data and local point statistics. Second, we classify ground photons based on photon height percentiles. Third, tree and grass photons are classified based on the number of neighbors. We validated tree heights with field measurements (n = 55), finding a root-mean-square error (RMSE) of 1.82 m using SVDA, GEDI Level 2A (Geolocated Elevation and Height Metrics product): 1.33 m, and ATL08: 5.59 m. Our results indicate that the SVDA is effective in identifying canopy photons in savanna ecosystems, where ATL08 performs poorly. We further identify seasonal vegetation height changes with an emphasis on vegetation below 3 m; widespread height changes in this class from two wet-dry cycles show maximum seasonal changes of 1 m, possibly related to seasonal grass-height differences. Our study shows the difficulties of vegetation measurements in savanna ecosystems but provides the first estimates of seasonal biomass changes.
Deep geological repositories represent a promising solution for the final disposal of nuclear waste. Due to its low permeability, high sorption capacity and self-sealing potential, Opalinus Clay (OPA) is considered a suitable host rock formation for the long-term storage of nuclear waste in Switzerland and Germany. However, the clay formation is characterized by compositional and structural variabilities including the occurrence of carbonate- and quartz-rich layers, pronounced bedding planes as well as tectonic elements such as pre-existing fault zones and fractures, suggesting heterogeneous rock mass properties.
Characterizing the heterogeneity of host rock properties is therefore essential for safety predictions of future repositories. This includes a detailed understanding of the mechanical and hydraulic properties, deformation behavior and the underlying deformation processes for an improved assessment of the sealing integrity and long-term safety of a deep repository in OPA. Against this background, this thesis presents the results of deformation experiments performed on intact and artificially fractured specimens of the quartz-rich, sandy and clay-rich, shaly facies of OPA. The experiments focus on the influence of mineralogical composition on the deformation behavior as well as the reactivation and sealing properties of pre-existing faults and fractures at different boundary conditions (e.g., pressure, temperature, strain rate).
The anisotropic mechanical properties of the sandy facies of OPA are presented in the first section, which were determined from triaxial deformation experiments using dried and resaturated samples loaded at 0°, 45° and 90° to the bedding plane orientation. A Paterson-type deformation apparatus was used that allowed to investigate how the deformation behavior is influenced by the variation of confining pressure (50 – 100 MPa), temperature (25 – 200 °C), and strain rate (1 × 10-3 – 5 × 10-6 s-1). Constant strain rate experiments revealed brittle to semi-brittle deformation behavior of the sandy facies at the applied conditions. Deformation behavior showed a strong dependence on confining pressure, degree of water saturation as well as bedding orientation, whereas the variation of temperature and strain rate had no significant effect on deformation. Furthermore, the sandy facies displays higher strength and stiffness compared to the clay-rich shaly facies deformed at similar conditions by Nüesch (1991). From the obtained results it can be concluded that cataclastic mechanisms dominate the short-term deformation behavior of dried samples from both facies up to elevated pressure (<200 MPa) and temperature (<200 °C) conditions.
The second part presents triaxial deformation tests that were performed to investigate how structural discontinuities affect the deformation behavior of OPA and how the reactivation of preexisting faults is influenced by mineral composition and confining pressure. To this end, dried cylindrical samples of the sandy and shaly facies of OPA were used, which contained a saw-cut fracture oriented at 30° to the long axis. After hydrostatic pre-compaction at 50 MPa, constant strain rate deformation tests were performed at confining pressures of 5, 20 or 35 MPa. With increasing confinement, a gradual transition from brittle, highly localized fault slip including a stress drop at fault reactivation to semi-brittle deformation behavior, characterized by increasing delocalization and non-linear strain hardening without dynamic fault reactivation, can be observed. Brittle localization was limited by the confining pressure at which the fault strength exceeded the matrix yield strength, above which strain partitioning between localized fault slip and distributed matrix deformation occurred. The sandy facies displayed a slightly higher friction coefficient (≈0.48) compared to the shaly facies (≈0.4). In addition, slide-hold-slide tests were conducted, revealing negative or negligible frictional strengthening, which suggests stable creep and long-term weakness of faults in both facies of OPA. The conducted experiments demonstrate that dilatant brittle fault reactivation in OPA may be favored at high overconsolidation ratios and shallow depths, increasing the risk of seismic hazard and the creation of fluid pathways.
The final section illustrates how the sealing capacity of fractures in OPA is affected by mineral composition. Triaxial flow-through experiments using Argon-gas were performed with dried samples from the sandy and shaly facies of OPA containing a roughened, artificial fracture. Slate, graywacke, quartzite, natural fault gouge, and granite samples were also tested to highlight the influence of normal stress, mineralogy and diagenesis on the sustainability of fracture transmissivity. With increasing normal stress, a non-linear decrease of fracture transmissivity can be observed that resulted in a permanent reduction of transmissivity after stress release. The transmissivity of rocks with a high portion of strong minerals (e.g., quartz) and high unconfined compressive strength was less sensitive to stress changes. In accordance with this, the sandy facies of OPA displayed a higher initial transmissivity that was less sensitive to stress changes compared to the shaly facies. However, transmissivity of rigid slate was less sensitive to stress changes than the sandy facies of OPA, although the slate is characterized by a higher phyllosilicate content. This demonstrates that in addition to mineral composition, other factors such as the degree of metamorphism, cementation and consolidation have to be considered when evaluating the sealing capacity of phyllosilicate-rich rocks.
The results of this thesis highlighted the role of confining pressure on the failure behavior of intact and artificially fractured OPA. Although the quartz-rich sandy facies may be considered as being more favorable for underground constructions due to its higher shear strength and stiffness than the shaly facies, the results indicate that when fractures develop in the sandy facies, they are more conductive and remain more permeable compared to fractures in the clay-dominated shaly facies at a given stress. The results may provide the basis for constitutive models to predict the integrity and evolution of a future repository. Clearly, the influence of composition and consolidation, e.g., by geological burial and uplift, on the mechanical sealing behavior of OPA highlights the need for a detailed site-specific material characterization for a future repository.
Magmatic continental rifts often constitute nascent plate boundaries, yet long-term extension rates and transient rate changes associated with these early stages of continental breakup remain difficult to determine. Here, we derive a time-averaged minimum extension rate for the inner graben of the Northern Kenya Rift (NKR) of the East African Rift System for the last 0.5 m.y. We use the TanDEM-X science digital elevation model to evaluate fault-scarp geometries and determine fault throws across the volcano-tectonic axis of the inner graben of the NKR. Along rift-perpendicular profiles, amounts of cumulative extension are determined, and by integrating four new Ar-40/Ar-39 radiometric dates for the Silali volcano into the existing geochronology of the faulted volcanic units, time-averaged extension rates are calculated. This study reveals that in the inner graben of the NKR, the long-term extension rate based on mid-Pleistocene to recent brittle deformation has minimum values of 1.0-1.6 mm yr(-1), locally with values up to 2.0 mm yr(-1). A comparison with the decadal, geodetically determined extension rate reveals that at least 65% of the extension must be accommodated within a narrow, 20-km-wide zone of the inner rift. In light of virtually inactive border faults of the NKR, we show that extension is focused in the region of the active volcano-tectonic axis in the inner graben, thus highlighting the maturing of continental rifting in the NKR.
Eight d-metal-containing N-butylpyridinium ionic liquids (ILs) with the nominal composition (C4Py)2[Ni0.5M0.5Cl4] or (C4Py)2[Zn0.5M0.5Cl4] (M = Cu, Co, Mn, Ni, Zn; C4Py = N-butylpyridinium) were synthesized, characterized, and investigated for their optical properties. Single crystal and powder X-ray analysis shows that the compounds are isostructural to existing examples based on other d-metal ions. Inductively coupled plasma optical emission spectroscopy measurements confirm that the metal/metal ratio is around 50 : 50. UV-Vis spectroscopy shows that the optical absorption can be tuned by selection of the constituent metals. Moreover, the compounds can act as an optical sensor for the detection of gases such as ammonia as demonstrated via a simple prototype setup.
Urban surface runoff management via best management practices (BMP) and low impact development (LID) has earned significant recognition owing to positive environmental and ecological impacts. However, due to the complexity of the parameters involved, the estimation of LID efficiency in attenuating the urban surface runoff at the watershed scale is challenging. A planning analysis of employing Green Roofs and Infiltration Trenches as BMPs/LIDs practices for urban surface runoff control is presented in this study. A multi-objective optimization decision-making framework is established by coupling SWMM (Storm Water Management Model) with NSGA-II models to check the performance of BMPs/LIDs concerning the cost-benefit analysis of LID at the watershed scale. Two urbanized areas belonging to Central Delhi in India were used as case studies. The results showed that the SWMM model is useful in simulating optimization problems for managing urban surface runoff. The optimum scenarios efficiently minimized the urban runoff volume while maintaining the BMPs/LIDs implementation costs and size. With BMPs/LIDs implementation, the reduction in runoff volume increases as expenses increase initially; however, there is no noticeable reduction in flood volume after a certain threshold. Contrasted with the haphazard arrangement of BMPs/LIDs, the proposed approach demonstrates 22%-24% runoff reductions for the same expenditures in watershed 1 and 23%-26% in watershed 2. The result of the study provides insights into planning and management of the urban surface runoff control with LID practices. The proposed framework assists the hydrologists in optimum selection and placements of BMPs/LIDs practices to acquire the most extreme ecological advantages with the least expenses.
We present two new empirical models of radiation belt electron flux at geostationary orbit. GOES-15 measurements of 0.8 MeV electrons were used to train a Nonlinear Autoregressive with Exogenous input (NARX) neural network for both modeling GOES-15 flux values and an upper boundary condition scaling factor (BF). The GOES-15 flux model utilizes an input and feedback delay of 2 and 2 time steps (i.e., 5 min time steps) with the most efficient number of hidden layers set to 10. Magnetic local time, Dst, Kp, solar wind dynamic pressure, AE, and solar wind velocity were found to perform as predicative indicators of GOES-15 flux and therefore were used as the exogenous inputs. The NARX-derived upper boundary condition scaling factor was used in conjunction with the Versatile Electron Radiation Belt (VERB) code to produce reconstructions of the radiation belts during the period of July-November 1990, independent of in-situ observations. Here, Kp was chosen as the sole exogenous input to be more compatible with the VERB code. This Combined Release and Radiation Effects Satellite-era reconstruction showcases the potential to use these neural network-derived boundary conditions as a method of hindcasting the historical radiation belts. This study serves as a companion paper to another recently published study on reconstructing the radiation belts during Solar Cycles 17-24 (Saikin et al., 2021, ), for which the results featured in this paper were used.
Ground-motion models (GMMs) are often used to predict the random distribution of Spectral accelerations (SAs) at a site due to a nearby earthquake. In probabilistic seismic hazard and risk assessment, large earthquakes occurring close to a site are considered as critical scenarios. GMMs are expected to predict realistic SAs with low within-model uncertainty (sigma(mu)) for such rare scenarios. However, the datasets used to regress GMMs are usually deficient of data from critical scenarios. The (Kotha et al., A Regionally Adaptable Ground-Motion Model for Shallow Crustal Earthquakes in Europe Bulletin of Earthquake Engineering 18:4091-4125, 2020) GMM developed from the Engineering strong motion (ESM) dataset was found to predict decreasing short-period SAs with increasing M-W >= M-h = 6.2, and with large sigma(mu) at near-source distances <= 30km. In this study, we updated the parametrisation of the GMM based on analyses of ESM and the Near source strong motion (NESS) datasets. With M-h = 5.7, we could rectify the M-W scaling issue, while also reducing sigma(mu). at M-W >= M-h. We then evaluated the GMM against NESS data, and found that the SAs from a few large, thrust-faulting events in California, New Zealand, Japan, and Mexico are significantly higher than GMM median predictions. However, recordings from these events were mostly made on soft-soil geology, and contain anisotropic pulse-like effects. A more thorough non-ergodic treatment of NESS was not possible because most sites sampled unique events in very diverse tectonic environments. We provide an updated set of GMM coefficients,sigma(mu), and heteroscedastic variance models; while also cautioning against its application for M-W <= 4 in low-moderate seismicity regions without evaluating the homogeneity of M-W estimates between pan-European ESM and regional datasets.
Neoarchean (similar to 2.73-2.70 Ga) accretionary history of the eastern Dharwar Craton, India
(2022)
Cratonic mid-crustal plutons may contain supracrustal enclaves that preserve evidence of an earlier growth history. The Eastern Dharwar craton records Neoarchean two-stage accretionary sequential growth (2.70 and 2.55 Ga) and a chronology of their enclaves could refine orogenic models. To test whether the metamorphic history of their enclaves was related to any of these stages, phase equilibria modelling and combined Lu-Hf and Sm-Nd geochronology on garnet were conducted on metapsammite, now preserved as garnet-orthopyroxene-cordierite gneiss. Phase equilibria modelling indicates peak metamorphic conditions, similar to 850 degrees C and similar to 8.5 kbar (M1a), were followed by near isothermal decompression to 5-6 kbar (M1b) and isobaric cooling to similar to 800 degrees C (M1c). The thermobaric gradient related to peak metamorphic conditions, similar to 30 degrees C kbar(-1), is typical of collisional orogens. Regression of the whole-rock and garnet, for sample S17b, yield Lu-Hf isochron ages of 2733 +/- 29 Ma, and for sample S18, 2724 +/- 13 Ma. A Lu-Hf weighted mean age for the porphyroblastic garnet suggests growth at 2725.5 +/- 11.9 Ma during the M1a-M1b stages. In contrast, the whole-rock sample S17b and the garnet fractions yield a Sm-Nd isochron age of 2696 +/- 10 Ma. From sample S18 the whole rock, garnet fractions, and orthopyroxene yield an isochron age of 2683 +/- 15 Ma. The garnet Sm-Nd weighted mean age at 2692.0 +/- 8.3 Ma constrains the M1b-M1c stages. We suggest that the protoliths to these supracrustal enclaves were deposited in an arc tectonic setting and underwent thickening followed by heating during peeled-back lithospheric convergence. Therefore, the earliest of the craton-forming accretionary stages is preserved as the similar to 2.73 Ga granulite-facies enclaves, marginally older than the 2.70-2.65 Ga cratonic greenstone volcanism. Tectonic exhumation of these mid-crustal granulite enclaves was in response to the late-Proterozoic (similar to 1.7 Ga) Bhopalpatnam orogeny.
Nocardioides alcanivorans sp. nov., a novel hexadecane-degrading species isolated from plastic waste
(2022)
Strain NGK65(T), a novel hexadecane degrading, non-motile, Gram-positive, rod-to-coccus shaped, aerobic bacterium, was isolated from plastic polluted soil sampled at a landfill.
Strain NGK65(T) hydrolysed casein, gelatin, urea and was catalase-positive. It optimally grew at 28 degrees C. in 0-1% NaCl and at pH 7.5-8.0. Glycerol, D-glucose, arbutin, aesculin, salicin, potassium 5-ketogluconate. sucrose, acetate, pyruvate and hexadecane were used as sole carbon sources.
The predominant membrane fatty acids were iso-C-16:0 followed by iso-C(17:)0 and C-18:1 omega 9c. The major polar lipids were phosphatidylglycerol, phosphatidylethanolamine, phosphatidylinositol and hydroxyphosphatidylinositol.
The cell-wall peptidoglycan type was A3 gamma, with LL-diaminopimelic acid and glycine as the diagnostic amino acids. MK 8 (H-4) was the predominant menaquinone. Phylogenetic analysis based on 16S rRNA gene sequences indicated that strain NGK65(T) belongs to the genus Nocardioides (phylum Actinobacteria). appearing most closely related to Nocardioides daejeonensis MJ31(T) (98.6%) and Nocardioides dubius KSL-104(T) (98.3%).
The genomic DNA G+C content of strain NGK65(T) was 68.2%.
Strain NGK65(T) and the type strains of species involved in the analysis had average nucleotide identity values of 78.3-71.9% as well as digital DNA-DNA hybridization values between 22.5 and 19.7%, which clearly indicated that the isolate represents a novel species within the genus Nocardioides.
Based on phenotypic and molecular characterization, strain NGK65(T) can clearly be differentiated from its phylogenetic neighbours to establish a novel species, for which the name Nocardioides alcanivorans sp. nov. is proposed.
The type strain is NGK65(T) (=DSM 113112(T)=NCCB 100846(T)).
Due to the major role of greenhouse gas emissions in global climate change, the development of non-fossil energy technologies is essential. Deep geothermal energy represents such an alternative, which offers promising properties such as a high base load capability and a large untapped potential. The present work addresses barite precipitation within geothermal systems and the associated reduction in rock permeability, which is a major obstacle to maintaining high efficiency. In this context, hydro-geochemical models are essential to quantify and predict the effects of precipitation on the efficiency of a system.
The objective of the present work is to quantify the induced injectivity loss using numerical and analytical reactive transport simulations. For the calculations, the fractured-porous reservoirs of the German geothermal regions North German Basin (NGB) and Upper Rhine Graben (URG) are considered.
Similar depth-dependent precipitation potentials could be determined for both investigated regions (2.8-20.2 g/m3 fluid). However, the reservoir simulations indicate that the injectivity loss due to barite deposition in the NGB is significant (1.8%-6.4% per year) and the longevity of the system is affected as a result; this is especially true for deeper reservoirs (3000 m). In contrast, simulations of URG sites indicate a minor role of barite (< 0.1%-1.2% injectivity loss per year). The key differences between the investigated regions are reservoir thicknesses and the presence of fractures in the rock, as well as the ionic strength of the fluids. The URG generally has fractured-porous reservoirs with much higher thicknesses, resulting in a greater distribution of precipitates in the subsurface. Furthermore, ionic strengths are higher in the NGB, which accelerates barite precipitation, causing it to occur more concentrated around the wellbore. The more concentrated the precipitates occur around the wellbore, the higher the injectivity loss.
In this work, a workflow was developed within which numerical and analytical models can be used to estimate and quantify the risk of barite precipitation within the reservoir of geothermal systems. A key element is a newly developed analytical scaling score that provides a reliable estimate of induced injectivity loss. The key advantage of the presented approach compared to fully coupled reservoir simulations is its simplicity, which makes it more accessible to plant operators and decision makers. Thus, in particular, the scaling score can find wide application within geothermal energy, e.g., in the search for potential plant sites and the estimation of long-term efficiency.
The Mackenzie Delta (MD) is a permafrost-bearing region along the coasts of the Canadian Arctic which exhibits high sub-permafrost gas hydrate (GH) reserves. The GH occurring at the Mallik site in the MD is dominated by thermogenic methane (CH4), which migrated from deep conventional hydrocarbon reservoirs, very likely through the present fault systems. Therefore, it is assumed that fluid flow transports dissolved CH4 upward and out of the deeper overpressurized reservoirs via the existing polygonal fault system and then forms the GH accumulations in the Kugmallit-Mackenzie Bay Sequences. We investigate the feasibility of this mechanism with a thermo-hydraulic-chemical numerical model, representing a cross section of the Mallik site. We present the first simulations that consider permafrost formation and thawing, as well as the formation of GH accumulations sourced from the upward migrating CH4-rich formation fluid. The simulation results show that temperature distribution, as well as the thickness and base of the ice-bearing permafrost are consistent with corresponding field observations. The primary driver for the spatial GH distribution is the permeability of the host sediments. Thus, the hypothesis on GH formation by dissolved CH4 originating from deeper geological reservoirs is successfully validated. Furthermore, our results demonstrate that the permafrost has been substantially heated to 0.8-1.3 degrees C, triggered by the global temperature increase of about 0.44 degrees C and further enhanced by the Arctic Amplification effect at the Mallik site from the early 1970s to the mid-2000s.