Gold Open-Access
Refine
Document Type
- Article (239)
- Review (4)
- Other (3)
- Doctoral Thesis (1)
Is part of the Bibliography
- yes (247)
Keywords
- climate change (9)
- permafrost (9)
- reactive transport (5)
- PHREEQC (4)
- digital rock physics (3)
- micro-CT (3)
- model (3)
- precipitation (3)
- remote sensing (3)
- time series analysis (3)
- Arctic Ocean (2)
- Chinese loess (2)
- Germany (2)
- InSAR (2)
- Mediterranean (2)
- Russia (2)
- Sentinel-1 (2)
- Siberia (2)
- VIIRS DNB (2)
- adaptation (2)
- arctic (2)
- black carbon (2)
- change detection (2)
- climate change adaptation (2)
- climate-change (2)
- deep biosphere (2)
- dynamics (2)
- elastic properties (2)
- elastic rock properties (2)
- flooding (2)
- formation damage (2)
- geothermal energy (2)
- heterogeneity (2)
- insurance (2)
- lake monitoring (2)
- lakes (2)
- migration (2)
- nitrogen (2)
- numerical simulation (2)
- patterns (2)
- permeability (2)
- photogrammetry (2)
- phreeqc (2)
- playa (2)
- pollen (2)
- radial flow (2)
- radiocarbon (2)
- repository far-field (2)
- salt pan (2)
- sediment (2)
- sedimentary ancient DNA (2)
- sorption (2)
- subduction zone (2)
- tectonics (2)
- treeline (2)
- vulnerability (2)
- AHP (1)
- ALOS World 3D (1)
- ASTER GDEM (1)
- AgI (1)
- Andes (1)
- Antarctic ice (1)
- Apatite (U-TH)/HE (1)
- Aphanizomenon (1)
- Arabian carbonate platform (1)
- Arctic (1)
- Arctic Russia (1)
- Arctic lakes (1)
- Arctic ocean (1)
- Asia (1)
- Atmospheric chemistry (1)
- Atmospheric dynamics (1)
- Atmospheric science (1)
- Australia (1)
- Barremian-Aptian (1)
- Bayesian model (1)
- Beaufort Sea (1)
- Bentheim sandstone (1)
- Biomedical engineering (1)
- Blueschist (1)
- Bohemian Massif (1)
- Brain injuries (1)
- Brazilian Earth System Model (1)
- CDOM (1)
- CMIP5 (1)
- CMIP5 models (1)
- CO2 (1)
- CO2 emissions (1)
- Cantera (1)
- Central Asia (1)
- Central-Asia (1)
- Chilean Coastal Cordillera (1)
- China (1)
- Climate change (1)
- Code_Aster (1)
- Computational models (1)
- Computational neuroscience (1)
- Copernicus DEM (1)
- Cyanobium (1)
- DDM simulation (1)
- DEM analysis (1)
- DEM noise (1)
- DNA preservation (1)
- DOC (1)
- Data acquisition (1)
- Data processing (1)
- Digitale Gesteinsphysik (1)
- ERA5 (1)
- East Asian summer monsoon (1)
- East Siberia (1)
- Eastern Alps (1)
- Eclogite (1)
- Eger Rift (1)
- Eisenbuhl (1)
- Elastische Gesteinseigenschaften (1)
- Electrical and electronic engineering (1)
- Emission scenarios (1)
- Energy and society (1)
- Environmental impact (1)
- Environmental sciences (1)
- Erodium cicutarium (1)
- Europe (1)
- European air quality (1)
- Extreme (1)
- FTIR (1)
- Firmicutes (1)
- Fission-track thermochronology (1)
- Fluid Flow (1)
- Forschungsmethodik (1)
- Fourier analysis (1)
- GEDI (1)
- GIS (1)
- GNSS Reflectometry (1)
- GNSS-integrated water vapour (1)
- Gemmatimonadetes (1)
- Generalized Extreme Value distribution (1)
- Geochemistry (1)
- Geodynamics (1)
- Geology (1)
- Geomorphology (1)
- Geophysics (1)
- Global Navigation Satellite System (GNSS) (1)
- Greek Islands (1)
- Greenland ice cores (1)
- Guaymas Basin (1)
- HBV-light model (1)
- Halobacteria (1)
- Himalaya (1)
- Himalayas (1)
- ICESat-2 (1)
- IGP (1)
- IPCC AR5 scenarios (1)
- Indian summer monsoon (1)
- Ion-Beam (FIB) (1)
- Ionic liquids (1)
- Ionosphere (1)
- Kalahari (1)
- Kardia mine (1)
- Kenya Rift (1)
- L-Asterisk (1)
- L-moments estimation (1)
- LM (1)
- Lake Issyk‐ Kul (1)
- Landsat (1)
- Landsat-8 (1)
- Larix (1)
- Lawsonite (1)
- LiDAR (1)
- Lower crust (1)
- Lu/Hf dating of garnet (1)
- MCDM (1)
- Mallik (1)
- Metamorphic sole (1)
- Methane leakage (1)
- Middle Eocene deformation (1)
- Mikro-CT (1)
- Mont Terri (1)
- Nanorocks (1)
- Natural hazard (1)
- Nd isotopic composition (1)
- Nearshore zone (1)
- Neogene (1)
- Nepal (1)
- Neural circuits (1)
- North (1)
- Northern Alpine Foreland Basin (1)
- Northern Sierras Pampeanas (1)
- OLCI (1)
- Ocean color remote sensing (1)
- Ocean sciences (1)
- OpenForecast (1)
- OpenStreetMap (1)
- Ozone (1)
- Palynology (1)
- Paratethys (1)
- Permeabilitätsentwicklung (1)
- Pitzer (1)
- Planktothrix (1)
- Plate tectonics (1)
- Precipitation (1)
- Puna Plateau (1)
- Qilian mountains (1)
- RCP scenarios (1)
- RCP4,5 and 8,5 (1)
- RCP4.5 (1)
- RCP8.5 (1)
- Random Forest (1)
- SEM (1)
- SRTM (1)
- Santa-Barbara system (1)
- Sava River (1)
- Sentinel 2 (1)
- Sentinel 3 (1)
- Sentinel-2 (1)
- Shale gas (1)
- Shan (1)
- Slope (1)
- Snowella (1)
- Solid Earth sciences (1)
- South Asia (1)
- South China Sea (1)
- South Korea (1)
- Southern Norway (1)
- Soziale Medien (1)
- Sporopollenin (1)
- Steppe-desert (1)
- Subduction (1)
- Synthetische Sandsteine (1)
- TRMM (1)
- TanDEM-X (1)
- TanDEM-X DEM (1)
- Tauern Window; (1)
- Tectonics (1)
- Temperature (1)
- TerraSAR-X (1)
- TerraceM (1)
- Tian (1)
- Tibet (1)
- Turbidity retrieval (1)
- Twitter (1)
- UAV (1)
- Unterrichtsmethoden (1)
- Vinoren (1)
- Visual system (1)
- Volatiles (1)
- Volcan de Colima (1)
- Volcanology (1)
- WRF-Chem (1)
- Western Interior Basin (1)
- Western Qaidam Basin (1)
- WorldDEM (1)
- Yamal (1)
- Zelezna hurka (1)
- absorption (1)
- acceptability (1)
- acoustic emissions (1)
- aegean (1)
- age-heterogeneity (1)
- agemodeling (1)
- agricultural productivity (1)
- agroecosystem (1)
- aillikite (1)
- air pollution (1)
- airglow (1)
- albedo (1)
- ammonia (1)
- amphibole (1)
- amplicon sequencing (1)
- ancient DNA (1)
- apatite geochemistry and U-Pb dating (1)
- aquatic carbon cycle (1)
- architecture (1)
- arctic water bodies (1)
- argon dating (1)
- artificial light (1)
- artificial light at night (1)
- assessment (1)
- astronomical calibration (1)
- atmospheric aerosols (1)
- attitudes (1)
- automated monitoring (1)
- automatic lineament extraction (1)
- automatic processing (1)
- bacterial diversity (1)
- bacterial-community structure (1)
- bat fatalities (1)
- biodiversity decline (1)
- biogeochemistry (1)
- biomarker (1)
- bioturbation (1)
- blueSeis-3A sensors (1)
- boreal (1)
- bromeliad (1)
- brown (1)
- bryophytes (1)
- buildings (1)
- calibration (1)
- canopy height (1)
- carbon (1)
- carbon capture and storage (CCS) (1)
- carbon capture and utilization (CCU) (1)
- carbon cycle (1)
- carbon dioxide (1)
- carbon sequestration (1)
- carbon stabilization (1)
- carbon-cycle (1)
- carbonate (1)
- carbonate assimilation (1)
- carnallite (1)
- cementation (1)
- central south island (1)
- change (1)
- characteristics (1)
- charcoal (1)
- chemical-mechanical interaction (1)
- chemistry (1)
- chronology (AICC2012) (1)
- city district (1)
- clay (1)
- cleaner energy transitions (1)
- climat change (1)
- climate (1)
- climate adaptation (1)
- climate benefits (1)
- climate gradient (1)
- clinopyroxenite (1)
- clonal growth (1)
- clumped isotopes (1)
- co-benefits (1)
- coastal geomorphology (1)
- code (1)
- coexistence (1)
- coexisting phases (1)
- coherency (1)
- coherent noise (1)
- competition-defense trade-off (1)
- composite properties (1)
- composition (1)
- compound flood (1)
- connectivity (1)
- continental lithosphere (1)
- convection (1)
- convective available potential energy (CAPE) (1)
- core-shell (1)
- corridors (1)
- cost of electricity (COE) (1)
- crust (1)
- cyanobacteria sedimentation (1)
- damage (1)
- damage model (1)
- data fusion (1)
- data-mining (1)
- deep carbon cycle (1)
- defense against predation (1)
- deforestation (1)
- deformation monitoring (1)
- density-driven (1)
- density-driven flow (1)
- deterministic (1)
- diamond (1)
- diatom indicator species (1)
- dietary patterns (1)
- diffusion (1)
- digital core reconstruction (1)
- digital twin (1)
- directivity (1)
- disaster risk (1)
- disaster risk reduction (1)
- discharge (1)
- discharge time series (1)
- dissolved silica concentration (1)
- disturbance tracking (1)
- diversity (1)
- downscaling (1)
- dryland (1)
- dust (1)
- dust sources (1)
- early eocene (1)
- earth system model (1)
- earthquake (1)
- earthquake source (1)
- earthquake swarms (1)
- eclogite (1)
- ecological succession (1)
- ecology (1)
- economics (1)
- ecosystem service cascade (1)
- ecosystem services (1)
- ecosystems (1)
- electrical resistivity tomography (1)
- electricity generation (1)
- electrons (1)
- emergency management (1)
- empirical (1)
- energy transfer (1)
- environmental (1)
- environmental reconstruction (1)
- eocene thermal maximum (1)
- equations (1)
- equatorial ionosphere (1)
- erosion (1)
- evaporite minerals (1)
- evaporites (1)
- exposure (1)
- extension (1)
- extensional tectonics (1)
- extinction event (1)
- extracellular DNA (1)
- extreme rainfall (1)
- facies (1)
- fault bend (1)
- faults (1)
- finite element modeling (1)
- fire (1)
- flat-slab (1)
- flat-slab subduction (1)
- flood risk analysis (1)
- flood risk management (1)
- floodplains (1)
- floods (1)
- flow-through experiment (1)
- flowering time (1)
- flux tracking (1)
- food systems (1)
- food web (1)
- forecasting (1)
- forecasting (1922, 4315, 7924, 7964) (1)
- forest (1)
- forest structure (1)
- fracture aperture (1)
- fracturing (1)
- frictional properties (1)
- functional response (1)
- functional trait correlations (1)
- functional traits (1)
- gas flux (1)
- gas hydrate (1)
- gas supply conditions (1)
- geo-bio interaction (1)
- geochronology (1)
- geogenic degassing (1)
- geomagnetic observatory data (1)
- geomagnetic storm drivers (1)
- geomorphic markers (1)
- glacial (1)
- glacial geomorphology (1)
- glacier (1)
- glacier chemistry (1)
- global (1)
- global change (1)
- global environmental change (1)
- global monsoon (1)
- governance (1)
- grain size (1)
- grain-size distribution (1)
- green infrastructure (1)
- green-green dilemma (1)
- greenhouse gas (1)
- greenhouse gas emissions (1)
- ground motion fields (1)
- groundwater vulnerability (1)
- growth rate (1)
- gulf coastal plain (1)
- hanging wall (1)
- hazards (1)
- heat-stress (1)
- hellenic volcanic arc (1)
- high mountain environments (1)
- high spatial resolution (1)
- high-resolution record (1)
- historical geomagnetic storms (1)
- holocene (1)
- huddle test (1)
- hybrid (1)
- hybrid pumped hydro power storage (1)
- hydrate formation (1)
- hydrate formation process (1)
- hydraulic (1)
- hydro power (1)
- hydro storage (1)
- hydro-geological system (1)
- hydrological drought (1)
- hydrological modeling (1)
- hydrology (1)
- hyperspectral (1)
- hysteresis (1)
- ice complex (1)
- ice sheet (1)
- ilmenite (1)
- image enhancement (1)
- impact (1)
- in (1)
- indirect facilitation (1)
- individual-based model (1)
- induced seismicity (1)
- injectivity (1)
- insights (1)
- interception (1)
- internal structure (1)
- intracellular DNA (1)
- invasive species (1)
- ionic conductivity (1)
- ionospheric irregularities (1)
- lake Donggi Cona (1)
- lake catchments (1)
- lake stratification (1)
- landscape evolution (1)
- landscape transience (1)
- lanthanides (1)
- lanthanoid migration (1)
- larch (1)
- last glacial period (1)
- lateral fluxes (1)
- lava dome (1)
- leading edge (1)
- legume-grass mixture (1)
- length measurements (1)
- level (1)
- lichens (1)
- lidar (1)
- life-history strategy (1)
- light pollution (1)
- lithobiont (1)
- livestock (1)
- load stress (1)
- local adaptation (1)
- local food (1)
- low-relief (1)
- machine learning (1)
- machine learning (0555) (1)
- magma-poor (1)
- magnetosphere (1)
- magnetostratigraphy (1)
- mapping (1)
- margins (1)
- marine sediments (1)
- marine terraces (1)
- mass (1)
- melilitite (1)
- melt inclusions (1)
- metabarcoding (1)
- metagenomic analysis (1)
- metal recycling plants (1)
- metasomatism (1)
- meteoric 10Be (1)
- methane (1)
- methane hydrate (1)
- methanol (1)
- mica (1)
- micro-CT scan (1)
- microbial diversity (1)
- microbial lipid (1)
- microbiology (1)
- mixed cropping (1)
- mixed gas hydrates (1)
- modelling (1)
- models (1)
- modern pollen assemblages (1)
- mofette systems (1)
- moment tensor (1)
- monitoring (1)
- monsoon (1)
- morphological traits (1)
- morphometry (1)
- movement (1)
- multi-stage diagenesis (1)
- multi-year insurance (1)
- multisensor system (1)
- multitemporal (1)
- multivariate statistics (1)
- nature-based solutions (NBS) (1)
- nearshore zone (1)
- neotectonics (1)
- night lights (1)
- nightlights (1)
- non-vascular epiphyte (1)
- normal faults (1)
- numerical (1)
- numerical simulations (1)
- ocean acidification (1)
- ocean color remote sensing (1)
- oligocene climate (1)
- online (1)
- open (1)
- open data (1)
- operational service (1)
- optical remote sensing (1)
- origination (1)
- orogenic peridotite (1)
- oxidizer (1)
- pH (1)
- paleoceanography (1)
- paleosol sequence (1)
- parallel evolution (1)
- part 1 (1)
- participation (1)
- passive microwave (1)
- peatlands (1)
- pedogenic carbonate (1)
- perennial crop (1)
- peri-urban agriculture (1)
- permafost (1)
- permafrost lowlands (1)
- permafrost thaw (1)
- permeability evolution (1)
- permeability-porosity relationship (1)
- phenotypic plasticity (1)
- phlogopite (1)
- phylogenetic and taxonomic plant diversity (1)
- phytoplankton (1)
- picocyanobacteria diversity (1)
- planktonic foraminifera (1)
- planning indicators (1)
- plate reconstructions (1)
- plateau (1)
- point cloud (1)
- pollen-vegetation-climate relationships (1)
- pollution (1)
- pollution indices (1)
- pore-scale (1)
- porous media (1)
- potentially toxic elements (1)
- precaution (1)
- predator-prey cycles (1)
- preferences (1)
- preparedness (1)
- pressure solution (1)
- process simulation (1)
- processes and landforms of glacial erosion (1)
- protists (1)
- provenance (1)
- psychrotolerant (1)
- quantitative vegetation reconstruction (1)
- quartz sandstone (1)
- radiation belts (1)
- radiation-damage (1)
- radio science (1)
- radio wave propagation (1)
- radiocarbon dating (1)
- rain attenuation (1)
- rain effect (1)
- reactivation (1)
- record (1)
- recovery (1)
- redox processes (1)
- reduction (1)
- regional flood frequency analysis (1)
- relative pollen productivity (1)
- remote (1)
- remote sensing data (1)
- renewable energy (1)
- residential buildings (1)
- response time (1)
- rheology (1)
- risk (1)
- risk management (1)
- risk reduction (1)
- river (1)
- river dynamics (1)
- river floods (1)
- river restoration (1)
- river sediment (1)
- rock (1)
- rock glacier (1)
- rotational seismology (1)
- runoff (1)
- saline groundwater (1)
- saline pan (1)
- saline pan cycle (1)
- salt dissolution (1)
- san andreas fault (1)
- santorini (1)
- savanna (1)
- scale (1)
- scaling (1)
- sea (1)
- sea level rise (1)
- sea level trends (1)
- sea surface height (1)
- sea turtle (1)
- sea-level projections (1)
- security and economy (1)
- sediment dynamics (1)
- sediment routing (1)
- sediment routing system (1)
- sediment transit time (1)
- sediment traps (1)
- sedimentary ancient DNA (sedaDNA) (1)
- seismicity (1)
- seismicity model (1)
- sensing (1)
- sensitivity (1)
- sentinel-1 (1)
- shallow earthquakes (1)
- shallow-water chemostratigraphy (1)
- shoreline movement direction (1)
- shoreline movement rate (1)
- shotgun sequencing (1)
- signal propagation (1)
- silica concentration (1)
- silvopastoral systems (1)
- simulation (1)
- situ Raman spectroscopy (1)
- slags (1)
- slip partitioning (1)
- slope aspect (1)
- snow (1)
- snow characteristics (1)
- snow grain size (1)
- software (1)
- soil formation (1)
- soils (1)
- solid Earth degassing (1)
- solid earth (1)
- source direction (1)
- source time function (1)
- source-to-sink (1)
- southern hemisphere oceans (1)
- special sensor microwave imager (1)
- special sensor microwave imager/sounder (1)
- specific surface area (1)
- spectral analysis (1)
- spinel (1)
- stable carbon isotopes (1)
- storm surge (1)
- strain localization (1)
- stratigraphy (1)
- structural mapping (1)
- structure from motion (1)
- subarctic (1)
- submarine gas vents (1)
- submarine volcanism (1)
- subsidence (1)
- subsurface life (1)
- sulfate reduction (1)
- summer monsoon (1)
- supervised machine learning (1)
- surface type-function-concept (1)
- suspended sediment (1)
- synthesis gas (1)
- synthetic sandstone (1)
- systems (1)
- taiga (1)
- techno-economic model (1)
- tectonic geomorphology (1)
- tectonic inversion (1)
- temperature sensor (1)
- thaw (1)
- thermal properties (1)
- thermodynamic equilibrium (1)
- thermokarst (1)
- thermokarst lakes (1)
- thermophiles; (1)
- thickness-displacement relationships (1)
- tibetan plateau (1)
- tide gauge observations (1)
- time-scale (1)
- time-series mapping (1)
- titanite (1)
- trait-based approaches (1)
- transfer (1)
- transferability (1)
- transit-time (1)
- transport processes (1)
- travel time distributions (1)
- tropical biodiversity (1)
- tropical swamp (1)
- tundra (1)
- tundra-taiga ecotone (1)
- underground coal gasification (1)
- underground coal gasification (UCG) (1)
- upconversion nanoparticles (1)
- uranium speciation (1)
- urban green (1)
- urban planning practice (1)
- urbanization (1)
- vascular epiphyte (1)
- very high resolution imagery (1)
- water (1)
- water age (1)
- water quality (1)
- water rock interaction (1)
- water storage capacity (1)
- weakening mechanisms (1)
- weathering (1)
- wildflower mixture (1)
- wind energy production (1)
- wind energy-biodiversity conflict (1)
- wind speed (1)
- yedoma (1)
- yedoma ice complex (1)
Institute
- Institut für Geowissenschaften (247) (remove)
Understanding eastern African paleoclimate is critical for contextualizing early human evolution, adaptation, and dispersal, yet Pleistocene climate of this region and its governing mechanisms remain poorly understood due to the lack of long, orbitally-resolved, terrestrial paleoclimate records. Here we present leaf wax hydrogen isotope records of rainfall from paleolake sediment cores from key time windows that resolve long-term trends, variations, and high-latitude effects on tropical African precipitation. Eastern African rainfall was dominantly controlled by variations in low-latitude summer insolation during most of the early and middle Pleistocene, with little evidence that glacial-interglacial cycles impacted rainfall until the late Pleistocene. We observe the influence of high-latitude-driven climate processes emerging from the last interglacial (Marine Isotope Stage 5) to the present, an interval when glacial-interglacial cycles were strong and insolation forcing was weak. Our results demonstrate a variable response of eastern African rainfall to low-latitude insolation forcing and high-latitude-driven climate change, likely related to the relative strengths of these forcings through time and a threshold in monsoon sensitivity. We observe little difference in mean rainfall between the early, middle, and late Pleistocene, which suggests that orbitally-driven climate variations likely played a more significant role than gradual change in the relationship between early humans and their environment.
The immense advances in computer power achieved in the last decades have had a significant impact in Earth science, providing valuable research outputs that allow the simulation of complex natural processes and systems, and generating improved forecasts. The development and implementation of innovative geoscientific software is currently evolving towards a sustainable and efficient development by integrating models of different aspects of the Earth system. This will set the foundation for a future digital twin of the Earth. The codification and update of this software require great effort from research groups and therefore, it needs to be preserved for its reuse by future generations of geoscientists. Here, we report on Geo-Soft-CoRe, a Geoscientific Software & Code Repository, hosted at the archive DIGITAL.CSIC. This is an open source, multidisciplinary and multiscale collection of software and code developed to analyze different aspects of the Earth system, encompassing tools to: 1) analyze climate variability; 2) assess hazards, and 3) characterize the structure and dynamics of the solid Earth. Due to the broad range of applications of these software packages, this collection is useful not only for basic research in Earth science, but also for applied research and educational purposes, reducing the gap between the geosciences and the society. By providing each software and code with a permanent identifier (DOI), we ensure its self-sustainability and accomplish the FAIR (Findable, Accessible, Interoperable and Reusable) principles. Therefore, we aim for a more transparent science, transferring knowledge in an easier way to the geoscience community, and encouraging an integrated use of computational infrastructure.
Pedogenic carbonate is widespread at mid latitudes where warm and dry conditions favor soil carbonate growth from spring to fall. The mechanisms and timing of pedogenic carbonate formation are more ambiguous in the tropical domain, where long periods of soil water saturation and high soil respiration enhance calcite dissolution. This paper provides stable carbon, oxygen and clumped isotope values from Quaternary and Miocene pedogenic carbonates in the tropical domain of Myanmar, in areas characterized by warm (>18°C) winters and annual rainfall up to 1,700 mm. We show that carbonate growth in Myanmar is delayed to the driest and coldest months of the year by sustained monsoonal rainfall from mid spring to late fall. The range of isotopic variability in Quaternary pedogenic carbonates can be solely explained by temporal changes of carbonate growth within the dry season, from winter to early spring. We propose that high soil moisture year-round in the tropical domain narrows carbonate growth to the driest months and makes it particularly sensitive to the seasonal distribution of rainfall. This sensitivity is also enabled by high winter temperatures, allowing carbonate growth to occur outside the warmest months of the year. This high sensitivity is expected to be more prominent in the geological record during times with higher temperatures and greater expansion of the tropical realm. Clumped isotope temperatures, δ13C and δ18O values of tropical pedogenic carbonates are impacted by changes of both rainfall seasonality and surface temperatures; this sensitivity can potentially be used to track past tropical rainfall distribution.
At the junction of greenhouse and icehouse climate states, the Eocene-Oligocene Transition (EOT) is a key moment in Cenozoic climate history. While it is associated with severe extinctions and biodiversity turnovers on land, the role of terrestrial climate evolution remains poorly resolved, especially the associated changes in seasonality. Some paleobotanical and geochemical continental records in parts of the Northern Hemisphere suggest the EOT is associated with a marked cooling in winter, leading to the development of more pronounced seasons (i.e., an increase in the mean annual range of temperature, MATR). However, the MATR increase has been barely studied by climate models and large uncertainties remain on its origin, geographical extent and impact. In order to better understand and describe temperature seasonality changes between the middle Eocene and the early Oligocene, we use the Earth system model IPSL-CM5A2 and a set of simulations reconstructing the EOT through three major climate forcings: pCO(2) decrease (1120, 840 and 560 ppm), the Antarctic ice-sheet (AIS) formation and the associated sea-level decrease. Our simulations suggest that pCO(2) lowering alone is not sufficient to explain the seasonality evolution described by the data through the EOT but rather that the combined effects of pCO(2) , AIS formation and increased continentality provide the best data-model agreement.pCO(2) decrease induces a zonal pattern with alternating increasing and decreasing seasonality bands particularly strong in the northern high latitudes (up to 8 degrees C MATR increase) due to sea-ice and surface albedo feedback. Conversely, the onset of the AIS is responsible for a more constant surface albedo yearly, which leads to a strong decrease in seasonality in the southern midlatitudes to high latitudes (> 40 degrees S). Finally, continental areas that emerged due to the sea-level lowering cause the largest increase in seasonality and explain most of the global heterogeneity in MATR changes (1MATR) patterns. The Delta MATR patterns we reconstruct are generally consistent with the variability of the EOT biotic crisis intensity across the Northern Hemisphere and provide insights on their underlying mechanisms.
Quantifying the resilience of vegetated ecosystems is key to constraining both present-day and future global impacts of anthropogenic climate change. Here we apply both empirical and theoretical resilience metrics to remotely-sensed vegetation data in order to examine the role of water availability and variability in controlling vegetation resilience at the global scale. We find a concise global relationship where vegetation resilience is greater in regions with higher water availability. We also reveal that resilience is lower in regions with more pronounced inter-annual precipitation variability, but find less concise relationships between vegetation resilience and intra-annual precipitation variability. Our results thus imply that the resilience of vegetation responds differently to water deficits at varying time scales. In view of projected increases in precipitation variability, our findings highlight the risk of ecosystem degradation under ongoing climate change.
Vegetation dynamics depend on both the amount of precipitation and its variability over time. Here, the authors show that vegetation resilience is greater where water availability is higher and where precipitation is more stable from year to year.
Changes in snowpack associated with climatic warming has drastic impacts on surface energy balance in the cryosphere. Yet, traditional monitoring techniques, such as punctual measurements in the field, do not cover the full snowpack spatial and temporal variability, which hampers efforts to upscale measurements to the global scale. This variability is one of the primary constraints in model development. In terms of spatial resolution, active microwaves (synthetic aperture radar - SAR) can address the issue and outperform methods based on passive microwaves. Thus, high-spatial-resolution monitoring of snow depth (SD) would allow for better parameterization of local processes that drive the spatial variability of snow. The overall objective of this study is to evaluate the potential of the TerraSAR-X (TSX) SAR sensor and the wave co-polar phase difference (CPD) method for characterizing snow cover at high spatial resolution. Consequently, we first (1) investigate SD and depth hoar fraction (DHF) variability between different vegetation classes in the Ice Creek catchment (Qikiqtaruk/Herschel Island, Yukon, Canada) using in situ measurements collected over the course of a field campaign in 2019; (2) evaluate linkages between snow characteristics and CPD distribution over the 2019 dataset; and (3) determine CPD seasonality considering meteorological data over the 2015-2019 period. SD could be extracted using the CPD when certain conditions are met. A high incidence angle (>30 circle) with a high topographic wetness index (TWI) (>7.0) showed correlation between SD and CPD (R2 up to 0.72). Further, future work should address a threshold of sensitivity to TWI and incidence angle to map snow depth in such environments and assess the potential of using interpolation tools to fill in gaps in SD information on drier vegetation types.
Organic carbon (OC) stored in Arctic permafrost represents one of Earth's largest and most vulnerable terrestrial carbon pools. Amplified climate warming across the Arctic results in widespread permafrost thaw. Permafrost deposits exposed at river cliffs and coasts are particularly susceptible to thawing processes. Accelerating erosion of terrestrial permafrost along shorelines leads to increased transfer of organic matter (OM) to nearshore waters. However, the amount of terrestrial permafrost carbon and nitrogen as well as the OM quality in these deposits is still poorly quantified. We define the OM quality as the intrinsic potential for further transformation, decomposition and mineralisation. Here, we characterise the sources and the quality of OM supplied to the Lena River at a rapidly eroding permafrost river shoreline cliff in the eastern part of the delta (Sobo-Sise Island). Our multi-proxy approach captures bulk elemental, molecu- lar geochemical and carbon isotopic analyses of Late Pleistocene Yedoma permafrost and Holocene cover deposits, discontinuously spanning the last similar to 52 kyr. We showed that the ancient permafrost exposed in the Sobo-Sise cliff has a high organic carbon content (mean of about 5 wt %). The oldest sediments stem from Marine Isotope Stage (MIS) 3 interstadial deposits (dated to 52 to 28 cal ka BP) and are overlaid by last glacial MIS 2 (dated to 28 to 15 cal ka BP) and Holocene MIS 1 (dated to 7-0 cal ka BP) deposits. The relatively high average chain length (ACL) index of n-alkanes along the cliff profile indicates a predominant contribution of vascular plants to the OM composition. The elevated ratio of isoand anteiso-branched fatty acids (FAs) relative to mid- and long-chain (C >= 20) n-FAs in the interstadial MIS 3 and the interglacial MIS 1 deposits suggests stronger microbial activity and consequently higher input of bacterial biomass during these climatically warmer periods. The overall high carbon preference index (CPI) and higher plant fatty acid (HPFA) values as well as high C/N ratios point to a good quality of the preserved OM and thus to a high potential of the OM for decomposition upon thaw. A decrease in HPFA values downwards along the profile probably indicates stronger OM decomposition in the oldest (MIS 3) deposits of the cliff. The characterisation of OM from eroding permafrost leads to a better assessment of the greenhouse gas potential of the OC released into river and nearshore waters in the future.
Arctic river deltas and deltaic near-shore zones represent important land-ocean transition zones influencing sediment dynamics and nutrient fluxes from permafrost-affected terrestrial ecosystems into the coastal Arctic Ocean. To accurately model fluvial carbon and freshwater export from rapidly changing river catchments as well as assess impacts of future change on the Arctic shelf and coastal ecosystems, we need to understand the sea floor characteristics and topographic variety of the coastal zones. To date, digital bathymetrical data from the poorly accessible, shallow, and large areas of the eastern Siberian Arctic shelves are sparse. We have digitized bathymetrical information for nearly 75 000 locations from large-scale (1 V 25000-1 V 500000) current and historical nautical maps of the Lena Delta and the Kolyma Gulf region in northeastern Siberia. We present the first detailed and seamless digital models of coastal zone bathymetry for both delta and gulf regions in 50 and 200m spatial resolution. We validated the resulting bathymetry layers using a combination of our own water depth measurements and a collection of available depth measurements, which showed a strong correlation (r>0.9). Our bathymetrical models will serve as an input for a high-resolution coupled hydrodynamic-ecosystem model to better quantify fluvial and coastal carbon fluxes to the Arctic Ocean, but they may be useful for a range of other studies related to Arctic delta and near-shore dynamics such as modeling of submarine permafrost, near-shore sea ice, or shelf sediment transport. The new digital high-resolution bathymetry products are available on the PANGAEA data set repository for the Lena Delta (https://doi.org/10.1594/PANGAEA.934045; Fuchs et al., 2021a) and Kolyma Gulf region (https://doi.org/10.1594/PANGAEA.934049; Fuchs et al., 2021b), respectively. Likewise, the depth validation data are available on PANGAEA as well (https://doi.org/10.1594/PANGAEA.933187; Fuchs et al., 2021c).
The Palaeocene-Eocene Thermal Maximum (ca. 56 million years ago) offers a primary analogue for future global warming and carbon cycle recovery. Yet, where and how massive carbon emissions were mitigated during this climate warming event remains largely unknown. Here we show that organic carbon burial in the vast epicontinental seaways that extended over Eurasia provided a major carbon sink during the Palaeocene-Eocene Thermal Maximum. We coupled new and existing stratigraphic analyses to a detailed paleogeographic framework and using spatiotemporal interpolation calculated ca. 720–1300 Gt organic carbon excess burial, focused in the eastern parts of the Eurasian epicontinental seaways. A much larger amount (2160–3900 Gt C, and when accounting for the increase in inundated shelf area 7400–10300 Gt C) could have been sequestered in similar environments globally. With the disappearance of most epicontinental seas since the Oligocene-Miocene, an effective negative carbon cycle feedback also disappeared making the modern carbon cycle critically dependent on the slower silicate weathering feedback.
An earthquake swarm affected the Bransfield Strait, Antarctica, a unique rift basin in transition from intra-arc rifting to ocean spreading. The swarm, counting similar to 85,000 volcano-tectonic earthquakes since August 2020, is located close to the Orca submarine volcano, previously considered inactive. Simultaneously, geodetic data reported up to similar to 11 cm north-westward displacement over King George Island. We use a broad variety of geophysical data and methods to reveal the complex migration of seismicity, accompanying the intrusion of 0.26-0.56 km(3) of magma. Strike-slip earthquakes mark the intrusion at depth, while shallower normal faulting the similar to 20 km long lateral growth of a dike. Seismicity abruptly decreased after a Mw 6.0 earthquake, suggesting the magmatic dike lost pressure with the slipping of a large fault. A seafloor eruption is likely, but not confirmed by sea surface temperature anomalies. The unrest documents episodic magmatic intrusion in the Bransfield Strait, providing unique insights into active continental rifting.
Oxygen (O-2) availability in soils is vital for plant growth and productivity. The transport and consumption of O-2 in the root zone is closely linked to soil moisture content, the spatial distribution of roots, as well as structure and heterogeneity of the surrounding soil. In this study, we measure three-dimensional root system architecture and the spatiotemporal dynamics of soil moisture (& theta;) and O-2 concentrations in the root zone of maize (Zea mays) via non-invasive imaging, and then construct and parameterize a reactive transport model based on the experimental data. The combination of three non-invasive imaging methods allowed for a direct comparison of simulation results with observations at high spatial and temporal resolution. In three different modeling scenarios, we investigated how the results obtained for different levels of conceptual complexity in the model were able to match measured & theta; and O-2 concentration patterns. We found that the modeling scenario that considers heterogeneous soil structure and spatial variability of hydraulic parameters (permeability, porosity, and van Genuchten & alpha; and n), better reproduced the measured & theta; and O-2 patterns relative to a simple model with a homogenous soil domain. The results from our combined imaging and modeling analysis reveal that experimental O-2 and water dynamics can be reproduced quantitatively in a reactive transport model, and that O-2 and water dynamics are best characterized when conditions unique to the specific system beyond the distribution of roots, such as soil structure and its effect on water saturation and macroscopic gas transport pathways, are considered.
Ore precipitation in porphyry copper systems is generally characterized by metal zoning (Cu-Mo to Zn-Pb-Ag), which is suggested to be variably related to solubility decreases during fluid cooling, fluid-rock interactions, partitioning during fluid phase separation and mixing with external fluids. Here, we present new advances of a numerical process model by considering published constraints on the temperature- and salinity-dependent solubility of Cu, Pb and Zn in the ore fluid. We quantitatively investigate the roles of vapor-brine separation, halite saturation, initial metal contents, fluid mixing and remobilization as first-order controls of the physical hydrology on ore formation. The results show that the magmatic vapor and brine phases ascend with different residence times but as miscible fluid mixtures, with salinity increases generating metal-undersaturated bulk fluids. The release rates of magmatic fluids affect the location of the thermohaline fronts, leading to contrasting mechanisms for ore precipitation: higher rates result in halite saturation without significant metal zoning, lower rates produce zoned ore shells due to mixing with meteoric water. Varying metal contents can affect the order of the final metal precipitation sequence. Redissolution of precipitated metals results in zoned ore shell patterns in more peripheral locations and also decouples halite saturation from ore precipitation.
Knowledge on the response of sediment export to recent climate change in glacierized areas in the European Alps is limited, primarily because long-term records of suspended sediment concentrations (SSCs) are scarce. Here we tested the estimation of sediment export of the past five decades using quantile regression forest (QRF), a nonparametric, multivariate regression based on random forest. The regression builds on short-term records of SSCs and long records of the most important hydroclimatic drivers (discharge, precipitation and air temperature - QPT). We trained independent models for two nested and partially glacier-covered catchments, Vent (98 km(2)) and Vernagt (11.4 km(2)), in the upper otztal in Tyrol, Austria (1891 to 3772 m a.s.l.), where available QPT records start in 1967 and 1975. To assess temporal extrapolation ability, we used two 2-year SSC datasets at gauge Vernagt, which are almost 20 years apart, for a validation. For Vent, we performed a five-fold cross-validation on the 15 years of SSC measurements. Further, we quantified the number of days where predictors exceeded the range represented in the training dataset, as the inability to extrapolate beyond this range is a known limitation of QRF. Finally, we compared QRF performance to sediment rating curves (SRCs). We analyzed the modeled sediment export time series, the predictors and glacier mass balance data for trends (Mann-Kendall test and Sen's slope estimator) and step-like changes (using the widely applied Pettitt test and a complementary Bayesian approach).Our validation at gauge Vernagt demonstrated that QRF performs well in estimating past daily sediment export (Nash-Sutcliffe efficiency (NSE) of 0.73) and satisfactorily for SSCs (NSE of 0.51), despite the small training dataset. The temporal extrapolation ability of QRF was superior to SRCs, especially in periods with high-SSC events, which demonstrated the ability of QRF to model threshold effects. Days with high SSCs tended to be underestimated, but the effect on annual yields was small. Days with predictor exceedances were rare, indicating a good representativity of the training dataset. Finally, the QRF reconstruction models outperformed SRCs by about 20 percent points of the explained variance.Significant positive trends in the reconstructed annual suspended sediment yields were found at both gauges, with distinct step-like increases around 1981. This was linked to increased glacier melt, which became apparent through step-like increases in discharge at both gauges as well as change points in mass balances of the two largest glaciers in the Vent catchment. We identified exceptionally high July temperatures in 1982 and 1983 as a likely cause. In contrast, we did not find coinciding change points in precipitation. Opposing trends at the two gauges after 1981 suggest different timings of "peak sediment". We conclude that, given large-enough training datasets, the presented QRF approach is a promising tool with the ability to deepen our understanding of the response of high-alpine areas to decadal climate change.
The Kolumbo submarine volcano in the southern Aegean (Greece) is associated with repeated seismic unrest since at least two decades and the causes of this unrest are poorly understood. We present a ten-month long microseismicity data set for the period 2006-2007. The majority of earthquakes cluster in a cone-shaped portion of the crust below Kolumbo. The tip of this cone coincides with a low Vp-anomaly at 2-4 km depth, which is interpreted as a crustal melt reservoir. Our data set includes several earthquake swarms, of which we analyze the four with the highest events numbers in detail. Together the swarms form a zone of fracturing elongated in the SW-NE direction, parallel to major regional faults. All four swarms show a general upward migration of hypocenters and the cracking front propagates unusually fast, compared to swarms in other volcanic areas. We conclude that the swarm seismicity is most likely triggered by a combination of pore-pressure perturbations and the re-distribution of elastic stresses. Fluid pressure perturbations are induced likely by obstructions in the melt conduits in a rheologically strong layer between 6 and 9 km depth. We conclude that the zone of fractures below Kolumbo is exploited by melts ascending from the mantle and filling the crustal melt reservoir. Together with the recurring seismic unrest, our study suggests that a future eruption is probable and monitoring of the Kolumbo volcanic system is highly advisable.
Wildfires play an essential role in the ecology of boreal forests.
In eastern Siberia, fire activity has been increasing in recent years, challenging the livelihoods of local communities. Intensifying fire regimes also increase disturbance pressure on the boreal forests, which currently protect the permafrost beneath from accelerated degradation.
However, long-term relationships between changes in fire regime and forest structure remain largely unknown.
We assess past fire-vegetation feedbacks using sedimentary proxy records from Lake Satagay, Central Yakutia, Siberia, covering the past c. 10,800 years.
Results from macroscopic and microscopic charcoal analyses indicate high amounts of burnt biomass during the Early Holocene, and that the present-day, low-severity surface fire regime has been in place since c. 4,500 years before present.
A pollen-based quantitative reconstruction of vegetation cover and a terrestrial plant record based on sedimentary ancient DNA metabarcoding suggest a pronounced shift in forest structure toward the Late Holocene.
Whereas the Early Holocene was characterized by postglacial open larch-birch woodlands, forest structure changed toward the modern, mixed larch-dominated closed-canopy forest during the Mid-Holocene.
We propose a potential relationship between open woodlands and high amounts of burnt biomass, as well as a mediating effect of dense larch forest on the climate-driven intensification of fire regimes.
Considering the anticipated increase in forest disturbances (droughts, insect invasions, and wildfires), higher tree mortality may force the modern state of the forest to shift toward an open woodland state comparable to the Early Holocene.
Such a shift in forest structure may result in a positive feedback on currently intensifying wildfires.
These new long-term data improve our understanding of millennial-scale fire regime changes and their relationships to changes of vegetation in Central Yakutia, where the local population is already being confronted with intensifying wildfire seasons.
Through its magnetic activity, the Sun governs the conditions in Earth's vicinity, creating space weather events, which have drastic effects on our space- and ground-based technology.
One of the most important solar magnetic features creating the space weather is the solar wind that originates from the coronal holes (CHs).
The identification of the CHs on the Sun as one of the source regions of the solar wind is therefore crucial to achieve predictive capabilities.
In this study, we used an unsupervised machine-learning method, k-means, to pixel-wise cluster the passband images of the Sun taken by the Atmospheric Imaging Assembly on the Solar Dynamics Observatory in 171, 193, and 211 angstrom in different combinations.
Our results show that the pixel-wise k-means clustering together with systematic pre- and postprocessing steps provides compatible results with those from complex methods, such as convolutional neural networks.
More importantly, our study shows that there is a need for a CH database where a consensus about the CH boundaries is reached by observers independently.
This database then can be used as the "ground truth," when using a supervised method or just to evaluate the goodness of the models.
Modern pollen-vegetation-climate relationships underpin palaeovegetation and palaeoclimate reconstructions from fossil pollen records. East Siberia is an ideal area for investigating the relationships between modern pollen assemblages and near natural vegetation under cold continental climate conditions. Reliable pollen-based quantitative vegetation and climate reconstructions are still scarce due to the limited number of modern pollen datasets. Furthermore, differences in pollen representation of samples from lake sediments and soils are not well understood. Here, we present a new pollen dataset of 48 moss/soil and 24 lake surface-sediment samples collected in Chukotka and central Yakutia in East Siberia. The pollen-vegetation-climate relationships were investigated by ordination analyses. Generally, tundra and taiga vegetation types can be well distinguished in the surface pollen assemblages. Moss/soil and lake samples contain generally similar pollen assemblages as revealed by a Procrustes comparison with some exceptions. Overall, modern pollen assemblages reflect the temperature and precipitation gradients in the study areas as revealed by constrained ordination analysis. We estimate the relative pollen productivity (RPP) of major taxa and the relevant source area of pollen (RSAP) for moss/soil samples from Chukotka and central Yakutia using Extended R-Value (ERV) analysis. The RSAP of the tundra-forest transition area in Chukotka and taiga area in central Yakutia are ca. 1300 and 360 m, respectively. For Chukotka, RPPs relative to both Poaceae and Ericaceae were estimated while RPPs for central Yakutia were relative only to Ericaceae. Relative to Ericaceae (reference taxon, RPP = 1), Larix, Betula, Picea, and Pinus are overrepresented while Alnus, Cyperaceae, Poaceae, and Salix are underrepresented in the pollen spectra. Our estimates are in general agreement with previously published values and provide the basis for reliable quantitative reconstructions of East Siberian vegetation.
As the recent permafrost thawing of northern Asia proceeds due to anthropogenic climate change, precise and detailed palaeoecological records from past warm periods are essential to anticipate the extent of future permafrost variations. Here, based on the modern relationship between permafrost and vegetation (represented by pollen assemblages), we trained a Random Forest model using pollen and permafrost data and verified its reliability to reconstruct the history of permafrost in northern Asia during the Holocene. An early Holocene (12-8 cal ka BP) strong thawing trend, a middle-to-late Holocene (8-2 cal ka BP) relatively slow thawing trend, and a late Holocene freezing trend of permafrost in northern Asia are consistent with climatic proxies such as summer solar radiation and Northern Hemisphere temperature. The extensive distribution of permafrost in northern Asia inhibited the spread of evergreen coniferous trees during the early Holocene warming and might have decelerated the enhancement of the East Asian summer monsoon (EASM) by altering hydrological processes and albedo. Based on these findings, we suggest that studies of the EASM should consider more the state of permafrost and vegetation in northern Asia, which are often overlooked and may have a profound impact on climate change in this region.
We investigate whether the distribution of maximum seasonal streamflow is significantly affected by catchment or climate state of the season/month ahead. We fit the Generalized Extreme Value (GEV) distribution to extreme seasonal streamflow for around 600 stations across Europe by conditioning the GEV location and scale parameters on 14 indices, which represent the season-ahead climate or catchment state. The comparison of these climate-informed models with the classical GEV distribution, with time-constant parameters, suggests that there is a substantial potential for seasonal forecasting of flood probabilities. The potential varies between seasons and regions. Overall, the season-ahead catchment wetness shows the highest potential, although climate indices based on large-scale atmospheric circulation, sea surface temperature or sea ice concentration also show some skill for certain regions and seasons. Spatially coherent patterns and a substantial fraction of climate-informed models are promising signs towards early alerts to increase flood preparedness already a season ahead.
Fast Holocene slip and localized strain along the Liquiñe-Ofqui strike-slip fault system, Chile
(2021)
In active tectonic settings dominated by strike-slip kinematics, slip partitioning across subparallel faults is a common feature; therefore, assessing the degree of partitioning and strain localization is paramount for seismic hazard assessments. Here, we estimate a slip rate of 18.8 +/- 2.0 mm/year over the past 9.0 +/- 0.1 ka for a single strand of the Liquirie-Ofqui Fault System, which straddles the Main Cordillera in Southern Chile. This Holocene rate accounts for similar to 82% of the trench-parallel component of oblique plate convergence and is similar to million-year estimates integrated over the entire fault system. Our results imply that strain localizes on a single fault at millennial time scale but over longer time scales strain localization is not sustained. The fast millennial slip rate in the absence of historical Mw> 6.5 earthquakes along the Liquine-Ofqui Fault System implies either a component of aseismic slip or Mw similar to 7 earthquakes involving multi-trace ruptures and > 150-year repeat times. Our results have implications for the understanding of strike-slip fault system dynamics within volcanic arcs and seismic hazard assessments.
Rupture directivity, implying a predominant earthquake rupture propagation direction, is typically inferred upon the identification of 2D azimuthal patterns of seismic observations for weak to large earthquakes using surface-monitoring networks. However, the recent increase of 3D monitoring networks deployed in the shallow subsurface and underground laboratories toward the monitoring of microseismicity allows to extend the directivity analysis to 3D modeling, beyond the usual range of magnitudes. The high-quality full waveforms recorded for the largest, decimeter-scale acoustic emission (AE) events during a meter-scale hydraulic fracturing experiment in granites at similar to 410 m depth allow us to resolve the apparent durations observed at each AE sensor to analyze 3D-directivity effects. Unilateral and (asymmetric) bilateral ruptures are then characterized by the introduction of a parameter kappa, representing the angle between the directivity vector and the station vector. While the cloud of AE activity indicates the planes of the hydrofractures, the resolved directivity vectors show off-plane orientations, indicating that rupture planes of microfractures on a scale of centimeters have different geometries. Our results reveal a general alignment of the rupture directivity with the orientation of the minimum horizontal stress, implying that not only the slip direction but also the fracture growth produced by the fluid injections is controlled by the local stress conditions.
Data driven high resolution modeling and spatial analyses of the COVID-19 pandemic in Germany
(2021)
The SARS-CoV-2 virus has spread around the world with over 100 million infections to date, and currently many countries are fighting the second wave of infections. With neither sufficient vaccination capacity nor effective medication, non-pharmaceutical interventions (NPIs) remain the measure of choice.
However, NPIs place a great burden on society, the mental health of individuals, and economics. Therefore the cost/benefit ratio must be carefully balanced and a target-oriented small-scale implementation of these NPIs could help achieve this balance.
To this end, we introduce a modified SEIRD-class compartment model and parametrize it locally for all 412 districts of Germany. The NPIs are modeled at district level by time varying contact rates. This high spatial resolution makes it possible to apply geostatistical methods to analyse the spatial patterns of the pandemic in Germany and to compare the results of different spatial resolutions.
We find that the modified SEIRD model can successfully be fitted to the COVID-19 cases in German districts, states, and also nationwide. We propose the correlation length as a further measure, besides the weekly incidence rates, to describe the current situation of the epidemic.
Woody plants are expanding into the Arctic in response to the warming climate. The impact on arctic plant communities is not well understood due to the limited knowledge about plant assembly rules.
Records of past plant diversity over long time series are rare. Here, we applied sedimentary ancient DNA metabarcoding targeting the P6 loop of the chloroplast trnL gene to a sediment record from Lake Ilirney (central Chukotka, Far Eastern Russia) covering the last 28 thousand years.
Our results show that forb-rich steppe-tundra and dwarf-shrub tundra dominated during the cold climate before 14 ka, while deciduous erect-shrub tundra was abundant during the warm period since 14 ka. Larix invasion during the late Holocene substantially lagged behind the likely warmest period between 10 and 6 ka, where the vegetation biomass could be highest.
We reveal highest richness during 28-23 ka and a second richness peak during 13-9 ka, with both periods being accompanied by low relative abundance of shrubs. During the cold period before 14 ka, rich plant assemblages were phylogenetically clustered, suggesting low genetic divergence in the assemblages despite the great number of species. This probably originates from environmental filtering along with niche differentiation due to limited resources under harsh environmental conditions. In contrast, during the warmer period after 14 ka, rich plant assemblages were phylogenetically overdispersed.
This results from a high number of species which were found to harbor high genetic divergence, likely originating from an erratic recruitment process in the course of warming. Some of our evidence may be of relevance for inferring future arctic plant assembly rules and diversity changes. By analogy to the past, we expect a lagged response of tree invasion. Plant richness might overshoot in the short term; in the long-term, however, the ongoing expansion of deciduous shrubs will eventually result in a phylogenetically more diverse community.
Efforts have been made in the past to enhance building exposure models on a regional scale with increasing spatial resolutions by integrating different data sources. This work follows a similar path and focuses on the downscaling of the existing SARA exposure model that was proposed for the residential building stock of the communes of Valparaiso and Vina del Mar (Chile). Although this model allowed great progress in harmonising building classes and characterising their differential physical vulnerabilities, it is now outdated, and in any case, it is spatially aggregated over large administrative units. Hence, to more accurately consider the impact of future earthquakes on these cities, it is necessary to employ more reliable exposure models. For such a purpose, we propose updating this existing model through a Bayesian approach by integrating ancillary data that has been made increasingly available from Volunteering Geo-Information (VGI) activities. Its spatial representation is also optimised in higher resolution aggregation units that avoid the inconvenience of having incomplete building-by-building footprints. A worst-case earthquake scenario is presented to calculate direct economic losses and highlight the degree of uncertainty imposed by exposure models in comparison with other parameters used to generate the seismic ground motions within a sensitivity analysis. This example study shows the great potential of using increasingly available VGI to update worldwide building exposure models as well as its importance in scenario-based seismic risk assessment.
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.
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.
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.
Sedimentary ancient DNA-based studies have been used to probe centuries of climate and environmental changes and how they affected cyanobacterial assemblages in temperate lakes. Due to cyanobacteria containing potential bloom-forming and toxin-producing taxa, their approximate reconstruction from sediments is crucial, especially in lakes lacking long-term monitoring data. To extend the resolution of sediment record interpretation, we used high-throughput sequencing, amplicon sequence variant (ASV) analysis, and quantitative PCR to compare pelagic cyanobacterial composition to that in sediment traps (collected monthly) and surface sediments in Lake Tiefer See. Cyanobacterial composition, species richness, and evenness was not significantly different among the pelagic depths, sediment traps and surface sediments (p > 0.05), indicating that the cyanobacteria in the sediments reflected the cyanobacterial assemblage in the water column. However, total cyanobacterial abundances (qPCR) decreased from the metalimnion down the water column. The aggregate-forming (Aphanizomenon) and colony-forming taxa (Snowella) showed pronounced sedimentation. In contrast, Planktothrix was only very poorly represented in sediment traps (meta- and hypolimnion) and surface sediments, despite its highest relative abundance at the thermocline (10 m water depth) during periods of lake stratification (May-October). We conclude that this skewed representation in taxonomic abundances reflects taphonomic processes, which should be considered in future DNA-based paleolimnological investigations.
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.
Planned decommissioning of coal-fired plants in Europe requires innovative technical and economic strategies to support coal regions on their path towards a climate-resilient future. The repurposing of open pit mines into hybrid pumped hydro power storage (HPHS) of excess energy from the electric grid, and renewable sources will contribute to the EU Green Deal, increase the economic value, stabilize the regional job market and contribute to the EU energy supply security. This study aims to present a preliminary phase of a geospatial workflow used to evaluate land suitability by implementing a multi-criteria decision making (MCDM) technique with an advanced geographic information system (GIS) in the context of an interdisciplinary feasibility study on HPHS in the Kardia lignite open pit mine (Western Macedonia, Greece). The introduced geospatial analysis is based on the utilization of the constraints and ranking criteria within the boundaries of the abandoned mine regarding specific topographic and proximity criteria. The applied criteria were selected from the literature, while for their weights, the experts' judgement was introduced by implementing the analytic hierarchy process (AHP), in the framework of the ATLANTIS research program. According to the results, seven regions were recognized as suitable, with a potential energy storage capacity from 1.09 to 5.16 GWh. Particularly, the present study's results reveal that 9.27% (212,884 m(2)) of the area had a very low suitability, 15.83% (363,599 m(2)) had a low suitability, 23.99% (550,998 m(2)) had a moderate suitability, 24.99% (573,813 m(2)) had a high suitability, and 25.92% (595,125 m(2)) had a very high suitability for the construction of the upper reservoir. The proposed semi-automatic geospatial workflow introduces an innovative tool that can be applied to open pit mines globally to identify the optimum design for an HPHS system depending on the existing lower reservoir.
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.
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.
The concepts of CO2 emission, global warming, climate change, and their environmental impacts are of utmost importance for the understanding and protection of the ecosystems.
Among the natural sources of gases into the atmosphere, the contribution of geogenic sources plays a crucial role. However, while subaerial emissions are widely studied, submarine outgassing is not yet well understood.
In this study, we review and catalog 122 literature and unpublished data of submarine emissions distributed in ten coastal areas of the Aegean Sea. This catalog includes descriptions of the degassing vents through in situ observations, their chemical and isotopic compositions, and flux estimations.
Temperatures and pH data of surface seawaters in four areas affected by submarine degassing are also presented.
This overview provides useful information to researchers studying the impact of enhanced seawater CO2 concentrations related either to increasing CO2 levels in the atmosphere or leaking carbon capture and storage systems.
In-depth understanding of the reorganization of the hydrological cycle in response to global climate change is crucial in highly sensitive regions like the eastern Mediterranean, where water availability is a major factor for socioeconomic and political development.
The sediments of Lake Lisan provide a unique record of hydroclimatic change during the last glacial to Holocene transition (ca. 24-11 ka) with its tremendous water level drop of similar to 240 m that finally led to its transition into the present hypersaline water body-the Dead Sea.
Here we utilize high-resolution sedimentological analyses from the marginal terraces and deep lake to reconstruct an unprecedented seasonal record of the last millennia of Lake Lisan. Aragonite varve formation in intercalated intervals of our record demonstrates that a stepwise long-term lake level decline was interrupted by almost one millennium of rising or stable water level.
Even periods of pronounced water level drops indicated by gypsum deposition were interrupted by decades of positive water budgets.
Our results thus highlight that even during major climate change at the end of the last glacial, decadal to millennial periods of relatively stable or positive moisture supply occurred which could have been an important premise for human sedentism.
The Arctic is rich in aquatic systems and experiences rapid warming due to climate change. The accelerated warming causes permafrost thaw and the mobilization of organic carbon.
When dissolved organic carbon is mobilized, this DOC can be transported to aquatic systems and degraded in the water bodies and further downstream. Here, we analyze the influence of different landscape components on DOC concentrations and export in a small (6.45 km(2)) stream catchment in the Lena River Delta.
The catchment includes lakes and ponds, with the flow path from Pleistocene yedoma deposits across Holocene non-yedoma deposits to the river outlet. In addition to DOC concentrations, we use radiocarbon dating of DOC as well as stable oxygen and hydrogen isotopes (delta O-18 and delta D) to assess the origin of DOC.
We find significantly higher DOC concentrations in the Pleistocene yedoma area of the catchment compared to the Holocene non-yedoma area with medians of 5 and 4.5 mg L-1 (p < 0.05), respectively. When yedoma thaw streams with high DOC concentration reach a large yedoma thermokarst lake, we observe an abrupt decrease in DOC concentration, which we attribute to dilution and lake processes such as mineralization. The DOC ages in the large thermokarst lake (between 3,428 and 3,637 C-14 y BP) can be attributed to a mixing of mobilized old yedoma and Holocene carbon.
Further downstream after the large thermokarst lake, we find progressively younger DOC ages in the stream water to its mouth, paired with decreasing DOC concentrations. This process could result from dilution with leaching water from Holocene deposits and/or emission of ancient yedoma carbon to the atmosphere. Our study shows that thermokarst lakes and ponds may act as DOC filters, predominantly by diluting incoming waters of higher DOC concentrations or by re-mineralizing DOC to CO2 and CH4.
Nevertheless, our results also confirm that the small catchment still contributes DOC on the order of 1.2 kg km(-2) per day from a permafrost landscape with ice-rich yedoma deposits to the Lena River.
In this study, we analyzed a large seismological dataset from temporary and permanent networks in the southern and eastern Alps to establish high-precision hypocenters and 1-D V-P and V-P/V-S models. The waveform data of a subset of local earthquakes with magnitudes in the range of 1-4.2 M-L were recorded by the dense, temporary SWATH-D network and selected stations of the AlpArray network between September 2017 and the end of 2018. The first arrival times of P and S waves of earthquakes are determined by a semi-automatic procedure. We applied a Markov chain Monte Carlo inversion method to simultaneously calculate robust hypocenters, a 1-D velocity model, and station corrections without prior assumptions, such as initial velocity models or earthquake locations. A further advantage of this method is the derivation of the model parameter uncertainties and noise levels of the data. The precision estimates of the localization procedure is checked by inverting a synthetic travel time dataset from a complex 3-D velocity model and by using the real stations and earthquakes geometry. The location accuracy is further investigated by a quarry blast test. The average uncertainties of the locations of the earthquakes are below 500m in their epicenter and similar to 1.7 km in depth. The earthquake distribution reveals seismicity in the upper crust (0-20 km), which is characterized by pronounced clusters along the Alpine frontal thrust, e.g., the Friuli-Venetia (FV) region, the Giudicarie-Lessini (GL) and Schio-Vicenza domains, the Austroalpine nappes, and the Inntal area. Some seismicity also occurs along the Periadriatic Fault. The general pattern of seismicity reflects head-on convergence of the Adriatic indenter with the Alpine orogenic crust. The seismicity in the FV and GL regions is deeper than the modeled frontal thrusts, which we interpret as indication for southward propagation of the southern Alpine deformation front (blind thrusts).
Dryland xeric conditions exert a deterministic effect on microbial communities, forcing life into refuge niches. Deposited rocks can form a lithic niche for microorganisms in desert regions. Mineral weathering is a key process in soil formation and the importance of microbial-driven mineral weathering for nutrient extraction is increasingly accepted. Advances in geobiology provide insight into the interactions between microorganisms and minerals that play an important role in weathering processes. In this study, we present the examination of the microbial diversity in dryland rocks from the Tsauchab River banks in Namibia. We paired culture-independent 16S rRNA gene amplicon sequencing with culture-dependent (isolation of bacteria) techniques to assess the community structure and diversity patterns. Bacteria isolated from dryland rocks are typical of xeric environments and are described as being involved in rock weathering processes. For the first time, we extracted extra- and intracellular DNA from rocks to enhance our understanding of potentially rock-weathering microorganisms. We compared the microbial community structure in different rock types (limestone, quartz-rich sandstone and quartz-rich shale) with adjacent soils below the rocks. Our results indicate differences in the living lithic and sublithic microbial communities.
The release of greenhouse gases from the large organic carbon stock in permafrost deposits in the circumarctic regions may accelerate global warming upon thaw. The extent of this positive climate feedback is thought to be largely controlled by the microbial degradability of the organic matter preserved in these sediments. In addition, weathering and oxidation processes may release inorganic carbon preserved in permafrost sediments as CO2, which is generally not accounted for. We used C-13 and C-14 analysis and isotopic mass balances to differentiate and quantify organic and inorganic carbon released as CO2 in the field from an active retrogressive thaw slump of Pleistocene-age Yedoma and during a 1.5-years incubation experiment. The results reveal that the dominant source of the CO2 released from freshly thawed Yedoma exposed as thaw mound is Pleistocene-age organic matter (48-80%) and to a lesser extent modern organic substrate (3-34%). A significant portion of the CO2 originated from inorganic carbon in the Yedoma (17-26%). The mixing of young, active layer material with Yedoma at a site on the slump floor led to the preferential mineralization of this young organic carbon source. Admixtures of younger organic substrates in the Yedoma thaw mound were small and thus rapidly consumed as shown by lower contributions to the CO2 produced during few weeks of aerobic incubation at 4 degrees C corresponding to approximately one thaw season. Future CO2 fluxes from the freshly thawed Yedoma will contain higher proportions of ancient inorganic (22%) and organic carbon (61-78%) as suggested by the results at the end, after 1.5 years of incubation. The increasing contribution of inorganic carbon during the incubation is favored by the accumulation of organic acids from microbial organic matter degradation resulting in lower pH values and, in consequence, in inorganic carbon dissolution. Because part of the inorganic carbon pool is assumed to be of pedogenic origin, these emissions would ultimately not alter carbon budgets. The results of this study highlight the preferential degradation of younger organic substrates in freshly thawed Yedoma, if available, and a substantial release of CO2 from inorganic sources.
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.
Regional flood frequency analysis (RFFA) is a powerful method for interrogating hydrological series since it combines observational time series from several sites within a region to estimate risk-relevant statistical parameters with higher accuracy than from single-site series. Since RFFA extreme value estimates depend on the shape of the selected distribution of the data-generating stochastic process, there is need for a suitable goodness-of-distributional-fit measure in order to optimally utilize given data. Here we present a novel, least-squares-based measure to select the optimal fit from a set of five distributions, namely Generalized Extreme Value (GEV), Generalized Logistic, Gumbel, Log-Normal Type III and Log-Pearson Type III. The fit metric is applied to annual maximum discharge series from six hydrological stations along the Sava River in South-eastern Europe, spanning the years 1961 to 2020. Results reveal that (1) the Sava River basin can be assessed as hydrologically homogeneous and (2) the GEV distribution provides typically the best fit. We offer hydrological-meteorological insights into the differences among the six stations. For the period studied, almost all stations exhibit statistically insignificant trends, which renders the conclusions about flood risk as relevant for hydrological sciences and the design of regional flood protection infrastructure.
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.
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.
Ice-rich permafrost has been subject to abrupt thaw and thermokarst formation in the past and is vulnerable to current global warming. The ice-rich permafrost domain includes Yedoma sediments that have never thawed since deposition during the late Pleistocene and Alas sediments that were formed by previous thermokarst processes during the Lateglacial and Holocene warming. Permafrost thaw unlocks organic carbon (OC) and minerals from these deposits and exposes OC to mineralization. A portion of the OC can be associated with iron (Fe), a redox-sensitive element acting as a trap for OC. Post-depositional thaw processes may have induced changes in redox conditions in these deposits and thereby affected Fe distribution and interactions between OC and Fe, with knock-on effects on the role that Fe plays in mediating present day OC mineralization. To test this hypothesis, we measured Fe concentrations and proportion of Fe oxides and Fe complexed with OC in unthawed Yedoma and previously thawed Alas deposits. Total Fe concentrations were determined on 1,292 sediment samples from the Yedoma domain using portable X-ray fluorescence; these concentrations were corrected for trueness using a calibration based on a subset of 144 samples measured by inductively coupled plasma optical emission spectrometry after alkaline fusion (R (2) = 0.95). The total Fe concentration is stable with depth in Yedoma deposits, but we observe a depletion or accumulation of total Fe in Alas deposits, which experienced previous thaw and/or flooding events. Selective Fe extractions targeting reactive forms of Fe on unthawed and previously thawed deposits highlight that about 25% of the total Fe is present as reactive species, either as crystalline or amorphous oxides, or complexed with OC, with no significant difference in proportions of reactive Fe between Yedoma and Alas deposits. These results suggest that redox driven processes during past thermokarst formation impact the present-day distribution of total Fe, and thereby the total amount of reactive Fe in Alas versus Yedoma deposits. This study highlights that ongoing thermokarst lake formation and drainage dynamics in the Arctic influences reactive Fe distribution and thereby interactions between Fe and OC, OC mineralization rates, and greenhouse gas emissions.
How biased are our models?
(2021)
Geophysical process simulations play a crucial role in the understanding of the subsurface. This understanding is required to provide, for instance, clean energy sources such as geothermal energy. However, the calibration and validation of the physical models heavily rely on state measurements such as temperature. In this work, we demonstrate that focusing analyses purely on measurements introduces a high bias. This is illustrated through global sensitivity studies. The extensive exploration of the parameter space becomes feasible through the construction of suitable surrogate models via the reduced basis method, where the bias is found to result from very unequal data distribution. We propose schemes to compensate for parts of this bias. However, the bias cannot be entirely compensated. Therefore, we demonstrate the consequences of this bias with the example of a model calibration.
Cyanobacteria are important primary producers in temperate freshwater ecosystems. However, studies on the seasonal and spatial distribution of cyanobacteria in deep lakes based on high-throughput DNA sequencing are still rare. In this study, we combined monthly water sampling and monitoring in 2019, amplicon sequence variants analysis (ASVs; a proxy for different species) and quantitative PCR targeting overall cyanobacteria abundance to describe the seasonal and spatial dynamics of cyanobacteria in the deep hard-water oligo-mesotrophic Lake Tiefer See, NE Germany. We observed significant seasonal variation in the cyanobacterial community composition (p < 0.05) in the epi- and metalimnion layers, but not in the hypolimnion. In winter-when the water column is mixed-picocyanobacteria (Synechococcus and Cyanobium) were dominant. With the onset of stratification in late spring, we observed potential niche specialization and coexistence among the cyanobacteria taxa driven mainly by light and nutrient dynamics. Specifically, ASVs assigned to picocyanobacteria and the genus Planktothrix were the main contributors to the formation of deep chlorophyll maxima along a light gradient. While Synechococcus and different Cyanobium ASVs were abundant in the epilimnion up to the base of the euphotic zone from spring to fall, Planktothrix mainly occurred in the metalimnetic layer below the euphotic zone where also overall cyanobacteria abundance was highest in summer. Our data revealed two potentially psychrotolerant (cold-adapted) Cyanobium species that appear to cope well under conditions of lower hypolimnetic water temperature and light as well as increasing sediment-released phosphate in the deeper waters in summer. The potential cold-adapted Cyanobium species were also dominant throughout the water column in fall and winter. Furthermore, Snowella and Microcystis-related ASVs were abundant in the water column during the onset of fall turnover. Altogether, these findings suggest previously unascertained and considerable spatiotemporal changes in the community of cyanobacteria on the species level especially within the genus Cyanobium in deep hard-water temperate lakes.
Strong hydroclimatic controls on vulnerability to subsurface nitrate contamination across Europe
(2020)
Subsurface contamination due to excessive nutrient surpluses is a persistent and widespread problem in agricultural areas across Europe. The vulnerability of a particular location to pollution from reactive solutes, such as nitrate, is determined by the interplay between hydrologic transport and biogeochemical transformations. Current studies on the controls of subsurface vulnerability do not consider the transient behaviour of transport dynamics in the root zone. Here, using state-of-the-art hydrologic simulations driven by observed hydroclimatic forcing, we demonstrate the strong spatiotemporal heterogeneity of hydrologic transport dynamics and reveal that these dynamics are primarily controlled by the hydroclimatic gradient of the aridity index across Europe. Contrasting the space-time dynamics of transport times with reactive timescales of denitrification in soil indicate that similar to 75% of the cultivated areas across Europe are potentially vulnerable to nitrate leaching for at least onethird of the year. We find that neglecting the transient nature of transport and reaction timescale results in a great underestimation of the extent of vulnerable regions by almost 50%. Therefore, future vulnerability and risk assessment studies must account for the transient behaviour of transport and biogeochemical transformation processes.
Water bodies are a highly abundant feature of Arctic permafrost ecosystems and strongly influence their hydrology, ecology and biogeochemical cycling. While very high resolution satellite images enable detailed mapping of these water bodies, the increasing availability and abundance of this imagery calls for fast, reliable and automatized monitoring. This technical work presents a largely automated and scalable workflow that removes image noise, detects water bodies, removes potential misclassifications from infrastructural features, derives lake shoreline geometries and retrieves their movement rate and direction on the basis of ortho-ready very high resolution satellite imagery from Arctic permafrost lowlands. We applied this workflow to typical Arctic lake areas on the Alaska North Slope and achieved a successful and fast detection of water bodies. We derived representative values for shoreline movement rates ranging from 0.40-0.56 m yr(-1) for lake sizes of 0.10 ha-23.04 ha. The approach also gives an insight into seasonal water level changes. Based on an extensive quantification of error sources, we discuss how the results of the automated workflow can be further enhanced by incorporating additional information on weather conditions and image metadata and by improving the input database. The workflow is suitable for the seasonal to annual monitoring of lake changes on a sub-meter scale in the study areas in northern Alaska and can readily be scaled for application across larger regions within certain accuracy limitations.
The subsurface is a temporally dynamic and spatially heterogeneous compartment of the Earth's critical zone, and biogeochemical transformations taking place in this compartment are crucial for the cycling of nutrients.
The impact of spatial heterogeneity on such microbially mediated nutrient cycling is not well known, which imposes a severe challenge in the prediction of in situ biogeochemical transformation rates and further of nutrient loading contributed by the groundwater to the surface water bodies.
Therefore, we used a numerical modelling approach to evaluate the sensitivity of groundwater microbial biomass distribution and nutrient cycling to spatial heterogeneity in different scenarios accounting for various residence times.
The model results gave us an insight into domain characteristics with respect to the presence of oxic niches in predominantly anoxic zones and vice versa depending on the extent of spatial heterogeneity and the flow regime.
The obtained results show that microbial abundance, distribution, and activity are sensitive to the applied flow regime and that the mobile (i.e. observable by groundwater sampling) fraction of microbial biomass is a varying, yet only a small, fraction of the total biomass in a domain. Furthermore, spatial heterogeneity resulted in anaerobic niches in the domain and shifts in microbial biomass between active and inactive states. The lack of consideration of spatial heterogeneity, thus, can result in inaccurate estimation of microbial activity. In most cases this leads to an overestimation of nutrient removal (up to twice the actual amount) along a flow path.
We conclude that the governing factors for evaluating this are the residence time of solutes and the Damkohler number (Da) of the biogeochemical reactions in the domain. We propose a relationship to scale the impact of spatial heterogeneity on nutrient removal governed by the logioDa.
This relationship may be applied in upscaled descriptions of microbially mediated nutrient cycling dynamics in the subsurface thereby resulting in more accurate predictions of, for example, carbon and nitrogen cycling in groundwater over long periods at the catchment scale.
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.
We present a new set of global and local sea‐level projections at example tide gauge locations under the RCP2.6, RCP4.5, and RCP8.5 emissions scenarios. Compared to the CMIP5‐based sea‐level projections presented in IPCC AR5, we introduce a number of methodological innovations, including (i) more comprehensive treatment of uncertainties, (ii) direct traceability between global and local projections, and (iii) exploratory extended projections to 2300 based on emulation of individual CMIP5 models. Combining the projections with observed tide gauge records, we explore the contribution to total variance that arises from sea‐level variability, different emissions scenarios, and model uncertainty. For the period out to 2300 we further breakdown the model uncertainty by sea‐level component and consider the dependence on geographic location, time horizon, and emissions scenario. Our analysis highlights the importance of local variability for sea‐level change in the coming decades and the potential value of annual‐to‐decadal predictions of local sea‐level change. Projections to 2300 show a substantial degree of committed sea‐level rise under all emissions scenarios considered and highlight the reduced future risk associated with RCP2.6 and RCP4.5 compared to RCP8.5. Tide gauge locations can show large ( > 50%) departures from the global average, in some cases even reversing the sign of the change. While uncertainty in projections of the future Antarctic ice dynamic response tends to dominate post‐2100, we see substantial differences in the breakdown of model variance as a function of location, time scale, and emissions scenario.
At the interface between the lithosphere and the atmosphere, the critical zone records the complex interactions between erosion, climate, geologic substrate, and life and can be directly monitored. Long data records (30 consecutive years for sediment yields) collected in the sparsely vegetated, steep, and small marly badland catchments of the Draix-Bleone Critical Zone Observatory (CZO), SE France, allow analyzing potential climatic controls on regolith dynamics and sediment export. Although widely accepted as a first-order control, rainfall variability does not fully explain the observed interannual variability in sediment export. Previous studies in this area have suggested that frost-weathering processes could drive regolith production and potentially modulate the observed pattern of sediment export. Here, we define sediment export anomalies as the residuals from a predictive model with annual rainfall intensity above a threshold as the control. We then use continuous soil temperature data recorded at different locations over multiple years to highlight the role of different frost-weathering processes (i.e., ice segregation versus volumetric expansion) in regolith production. Several proxies for different frost-weathering processes have been calculated from these data and compared to the sediment export anomalies, with careful consideration of field data quality. Our results suggest that frost-cracking intensity (linked to ice segregation) can explain about half (47 %-64 %) of the sediment export anomalies. In contrast, the number of freeze-thaw cycles (linked to volumetric expansion) has only a minor impact on catchment sediment response. The time spent below 0 degrees C also correlates well with the sediment export anomalies and requires fewer field data to be calculated than the frost-cracking intensity. Thus, frost-weathering processes modulate sediment export by controlling regolith production in these catchments and should be taken into account when building predictive models of sediment export from these badlands under a changing climate.
The hyperthermal events of the Cenozoic, including the Paleocene-Eocene Thermal Maximum, provide an opportunity to investigate the potential effects of climate warming on marine ecosystems. Here, we examine the shallow benthic marine communities preserved in the late Cretaceous to Eocene strata on the Gulf Coastal Plain (United States). In stark contrast to the ecological shifts following the end-Cretaceous mass extinction, our data show that the early Cenozoic hyperthermals did not have a long-term impact on the generic diversity nor composition of the Gulf Coastal Plain molluscan communities. We propose that these communities were resilient to climate change because molluscs are better adapted to high temperatures than other taxa, as demonstrated by their physiology and evolutionary history. In terms of resilience, these communities differ from other shallow-water carbonate ecosystems, such as reef communities, which record significant changes during the early Cenozoic hyperthermals. These data highlight the strikingly different responses of community types, i.e., the almost imperceptible response of molluscs versus the marked turnover of foraminifera and reef faunas. The impact on molluscan communities may have been low because detrimental conditions did not devastate the entire Gulf Coastal Plain, allowing molluscs to rapidly recolonise vacated areas once harsh environmental conditions ameliorated.
The origin of Asian monsoons
(2020)
The Cenozoic inception and development of the Asian monsoons remain unclear and have generated much debate, as several hypotheses regarding circulation patterns at work in Asia during the Eocene have been proposed in the few last decades. These include (a) the existence of modern-like monsoons since the early Eocene; (b) that of a weak South Asian monsoon (SAM) and little to no East Asian monsoon (EAM); or (c) a prevalence of the Intertropical Convergence Zone (ITCZ) migrations, also referred to as Indonesian-Australian monsoon (I-AM). As SAM and EAM are supposed to have been triggered or enhanced primarily by Asian palaeogeographic changes, their possible inception in the very dynamic Eocene palaeogeographic context remains an open question, both in the modelling and field-based communities. We investigate here Eocene Asian climate conditions using the IPSL-CM5A2 (Sepulchre et al., 2019) earth system model and revised palaeogeographies. Our Eocene climate simulation yields atmospheric circulation patterns in Asia substantially different from modern conditions. A large high-pressure area is simulated over the Tethys ocean, which generates intense low tropospheric winds blowing southward along the western flank of the proto-Himalayan-Tibetan plateau (HTP) system. This low-level wind system blocks, to latitudes lower than 10 degrees N, the migration of humid and warm air masses coming from the Indian Ocean. This strongly contrasts with the modern SAM, during which equatorial air masses reach a latitude of 20-25 degrees N over India and southeastern China. Another specific feature of our Eocene simulation is the widespread subsidence taking place over northern India in the midtroposphere (around 5000 m), preventing deep convective updraught that would transport water vapour up to the condensation level. Both processes lead to the onset of a broad arid region located over northern India and over the HTP. More humid regions of high seasonality in precipitation encircle this arid area, due to the prevalence of the Intertropical Convergence Zone (ITCZ) migrations (or Indonesian-Australian monsoon, I-AM) rather than monsoons. Although the existence of this central arid region may partly result from the specifics of our simulation (model dependence and palaeogeographic uncertainties) and has yet to be confirmed by proxy records, most of the observational evidence for Eocene monsoons are located in the highly seasonal transition zone between the arid area and the more humid surroundings. We thus suggest that a zonal arid climate prevailed over Asia before the initiation of monsoons that most likely occurred following Eocene palaeogeographic changes. Our results also show that precipitation seasonality should be used with caution to infer the presence of a monsoonal circulation and that the collection of new data in this arid area is of paramount importance to allow the debate to move forward.
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.
New Zealand's Alpine Fault is a large, platebounding strike-slip fault, which ruptures in large (M-w > 8) earthquakes. We conducted field and laboratory analyses of fault rocks to assess its fault zone architecture. Results reveal that the Alpine Fault Zone has a complex geometry, comprising an anastomosing network of multiple slip planes that have accommodated different amounts of displacement. This contrasts with the previous perception of the Alpine Fault Zone, which assumes a single principal slip zone accommodated all displacement. This interpretation is supported by results of drilling projects and geophysical investigations. Furthermore, observations presented here show that the young, largely unconsolidated sediments that constitute the footwall at shallow depths have a significant influence on fault gouge rheological properties and structure.
We produce climate projections through the 21st century using the fractional energy balance equation (FEBE): a generalization of the standard energy balance equation (EBE). The FEBE can be derived from Budyko-Sellers models or phenomenologically through the application of the scaling symmetry to energy storage processes, easily implemented by changing the integer order of the storage (derivative) term in the EBE to a fractional value.
The FEBE is defined by three parameters: a fundamental shape parameter, a timescale and an amplitude, corresponding to, respectively, the scaling exponent h, the relaxation time tau and the equilibrium climate sensitivity (ECS). Two additional parameters were needed for the forcing: an aerosol recalibration factor alpha to account for the large aerosol uncertainty and a volcanic intermittency correction exponent upsilon. A Bayesian framework based on historical temperatures and natural and anthropogenic forcing series was used for parameter estimation. Significantly, the error model was not ad hoc but rather predicted by the model itself: the internal variability response to white noise internal forcing.
The 90 % credible interval (CI) of the exponent and relaxation time were h = [0.33, 0.44] (median = 0.38) and tau = [2.4, 7.0] (median = 4.7) years compared to the usual EBE h = 1, and literature values of tau typically in the range 2-8 years. Aerosol forcings were too strong, requiring a decrease by an average factor alpha = [0.2, 1.0] (median = 0.6); the volcanic intermittency correction exponent was upsilon = [0.15, 0.41] (median = 0.28) compared to standard values alpha = upsilon = 1. The overpowered aerosols support a revision of the global modern (2005) aerosol forcing 90 % CI to a narrower range [ -1.0, -0.2] W m(-2). The key parameter ECS in comparison to IPCC AR5 (and to the CMIP6 MME), the 90 % CI range is reduced from [1.5, 4.5] K ([2.0, 5.5] K) to [1.6, 2.4] K ([1.5, 2.2] K), with median value lowered from 3.0 K (3.7 K) to 2.0 K (1.8 K) Similarly we found for the transient climate response (TCR), the 90 % CI range shrinks from [1.0, 2.5] K ([1.2, 2.8] K) to [1.2, 1.8] K ([1.1, 1.6] K) and the median estimate decreases from 1.8 K (2.0 K) to 1.5 K (1.4 K). As often seen in other observational-based studies, the FEBE values for climate sensitivities are therefore somewhat lower but still consistent with those in IPCC AR5 and the CMIP6 MME. <br /> Using these parameters, we made projections to 2100 using both the Representative Concentration Pathway (RCP) and Shared Socioeconomic Pathway (SSP) scenarios, and compared them to the corresponding CMIP5 and CMIP6 multi-model ensembles (MMEs). The FEBE historical reconstructions (1880-2020) closely follow observations, notably during the 1998-2014 slowdown ("hiatus"). We also reproduce the internal variability with the FEBE and statistically validate this against centennial-scale temperature observations. Overall, the FEBE projections were 10 %-15 % lower but due to their smaller uncertainties, their 90 % CIs lie completely within the GCM 90 % CIs. This agreement means that the FEBE validates the MME, and vice versa.
The effect of lithology on the relationship between denudation rate and chemical weathering pathways
(2022)
The denudation of rocks in mountain belts exposes a range of fresh minerals to the surface of the Earth that are chemically weathered by acidic and oxygenated fluids. The impact of the resulting coupling between denudation and weathering rates fundamentally depends on the types of minerals that are weathering. Whereas silicate weathering sequesters CO2, the combination of sulfide oxidation and carbonate dissolution emits CO2 to the atmosphere. Here, we combine the concentrations of dissolved major elements in stream waters with Be-10 basin-wide denudation rates from 35 small catchments in eastern Tibet to elucidate the importance of lithology in modulating the relationships between denudation rate, chemical weathering pathways, and CO2 consumption or release. Our catchments span 3 orders of magnitude in denudation rate in low-grade flysch, high-grade metapelites, and granitoid rocks. For each stream, we estimate the concentrations of solutes sourced from silicate weathering, carbonate dissolution, and sulfide oxidation using a mixing model. We find that for all lithologies, cation concentrations from silicate weathering are largely independent of denudation rate, but solute concentrations from carbonates and, where present, sulfides increase with increasing denudation rate. With increasing denudation rates, weathering may therefore shift from consuming to releasing CO2 in both (meta)sedimentary and granitoid lithologies. For a given denudation rate, we report dissolved solid concentrations and inferred weathering fluxes in catchments underlain by (meta)sedimentary rock that are 2-10 times higher compared to catchments containing granitoid lithologies, even though climatic and topographic parameters do not vary systematically between these catchments. Thus, varying proportions of exposed (meta)sedimentary and igneous rocks during orogenesis could lead to changes in the sequestration and release of CO2 that are independent of denudation rate.
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.
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.
Leaching zones within potash seams generally represent a significant risk to subsurface mining operations and the construction of technical caverns in salt rocks, but their temporal and spatial formation has been investigated only rudimentarily to date. To the knowledge of the authors, current reactive transport simulation implementations are not capable to address hydraulic-chemical interactions within potash salt. For this reason, a reactive transport model has been developed and complemented by an innovative approach to calculate the interchange of minerals and solution at the water-rock interface. Using this model, a scenario analysis was carried out based on a carnallite-bearing potash seam. The results show that the evolution of leaching zones depends on the mineral composition and dissolution rate of the original salt rock, and that the formation can be classified by the dimensionless parameters of Peclet (Pe) and Damkohler (Da). For Pe > 2 and Da > 1, a funnel-shaped leaching zone is formed, otherwise the dissolution front is planar. Additionally, Da > 1 results in the formation of a sylvinitic zone and a flow barrier. Most scenarios represent hybrid forms of these cases. The simulated shapes and mineralogies are confirmed by literature data and can be used to assess the hazard potential.
Barite scalings are a common cause of permanent formation damage to deep geothermal reservoirs. Well injectivity can be impaired because the ooling of saline fluids reduces the solubility of barite, and the continuous re-injection of supersaturated fluids forces barite to precipitate in the host rock. Stimulated reservoirs in the Upper Rhine Graben often have multiple relevant flow paths in the porous matrix and fracture zones, sometimes spanning multiple stratigraphical units to achieve the economically necessary injectivity. While the influence of barite scaling on injectivity has been investigated for purely porous media, the role of fractures within reservoirs consisting of both fractured and porous sections is still not well understood. Here, we present hydro-chemical simulations of a dual-layer geothermal reservoir to study the long-term impact of barite scale formation on well injectivity. Our results show that, compared to purely porous reservoirs, fractured porous reservoirs have a significantly reduced scaling risk by up to 50%, depending on the flow rate ratio of fractures. Injectivity loss is doubled, however, if the amount of active fractures is increased by one order of magnitude, while the mean fracture aperture is decreased, provided the fractured aquifer dictates the injection rate. We conclude that fractured, and especially hydraulically stimulated, reservoirs are generally less affected by barite scaling and that large, but few, fractures are favourable. We present a scaling score for fractured-porous reservoirs, which is composed of easily derivable quantities such as the radial equilibrium length and precipitation potential. This score is suggested for use approximating the scaling potential and its impact on injectivity of a fractured-porous reservoir for geothermal exploitation.
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.
The spatial pattern of extreme precipitation from 40 years of gauge data in the central Himalaya
(2022)
The topography of the Himalaya exerts a substantial control on the spatial distribution of monsoonal rainfall, which is a vital water source for the regional economy and population. But the occurrence of short-lived and high-intensity precipitation results in socio-economic losses. This study relies on 40 years of daily data from 204 ground stations in Nepal to derive extreme precipitation thresholds, amounts, and days at the 95th percentile. We additionally determine the precipitation magnitude-frequency relation. We observe that extreme precipitation amounts follow an almost uniform band parallel to topographic contour lines in the southern Himalaya mountains in central and eastern Nepal but not in western Nepal. The relationship of extreme precipitation indices with topographic relief shows that extreme precipitation thresholds decrease with increasing elevation, but extreme precipitation days increase in higher elevation areas. Furthermore, stations above 1 km elevation exhibit a power-law relation in the rainfall magnitude-frequency framework. Stations at higher elevations generally have lower values of power-law exponents than low elevation areas. This suggests a fundamentally different behaviour of the rainfall distribution and an increased occurrence of extreme rainfall storms in the high elevation areas of Nepal.
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 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.
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.
Barite scales in geothermal installations are a highly unwanted effect of circulating deep saline fluids. They build up in the reservoir if supersaturated fluids are re-injected, leading to irreversible loss of injectivity. A model is presented for calculating the total expected barite precipitation. To determine the related injectivity decline over time, the spatial precipitation distribution in the subsurface near the injection well is assessed by modelling barite growth kinetics in a radially diverging Darcy flow domain. Flow and reservoir properties as well as fluid chemistry are chosen to represent reservoirs subject to geothermal exploration located in the North German Basin (NGB) and the Upper Rhine Graben (URG) in Germany. Fluids encountered at similar depths are hotter in the URG, while they are more saline in the NGB. The associated scaling amount normalised to flow rate is similar for both regions. The predicted injectivity decline after 10 years, on the other hand, is far greater for the NGB (64%) compared to the URG (24%), due to the temperature- and salinity-dependent precipitation rate. The systems in the NGB are at higher risk. Finally, a lightweight score is developed for approximating the injectivity loss using the Damkohler number, flow rate and total barite scaling potential. This formula can be easily applied to geothermal installations without running complex reactive transport simulations.
This study, based on new and high quality in situ observations, quantifies for the first time, the individual contributions of light-absorbing aerosols (black carbon (BC), brown carbon (BrC) and dust) to aerosol absorption over the Indo-Gangetic Plain (IGP) and the Himalayan foothill region, a relatively poorly studied region with several sensitive ecosystems of global importance, as well as highly vulnerable populations. The annual and seasonal average single scattering albedo (SSA) over Kathmandu is the lowest of all the locations. The SSA over Kathmandu is < 0.89 during all seasons, which confirms the dominance of light-absorbing carbonaceous aerosols from local and regional sources over Kathmandu. It is observed here that the SSA decreases with increasing elevation, confirming the dominance of light absorbing carbonaceous aerosols at higher elevations. In contrast, the SSA over the IGP does not exhibit a pronounced spatial variation. BC dominates (>= 75%) the aerosol absorption over the IGP and the Himalayan foothills throughout the year. Higher BC concentration at elevated locations in the Himalayas leads to lower SSA at elevated locations in the Himalayas. The contribution of dust to aerosol absorption is higher throughout the year over the IGP than over the Himalayan foothills. The aerosol absorption over South Asia is very high, exceeding available observations over East Asia, and also exceeds previous model estimates. This quantification will be valuable as observational constraints to help improve regional simulations of climate change, impacts on the glaciers and the hydrological cycle, and will help to direct the focus towards BC as the main contributor to aerosol-induced warming in the region.
The SiDroForest (Siberian drone-mapped forest inventory) data collection is an attempt to remedy the scarcity of forest structure data in the circumboreal region by providing adjusted and labeled tree-level and vegetation plot-level data for machine learning and upscaling purposes. We present datasets of vegetation composition and tree and plot level forest structure for two important vegetation transition zones in Siberia, Russia; the summergreen-evergreen transition zone in Central Yakutia and the tundra-taiga transition zone in Chukotka (NE Siberia). The SiDroForest data collection consists of four datasets that contain different complementary data types that together support in-depth analyses from different perspectives of Siberian Forest plot data for multi-purpose applications. i. Dataset 1 provides unmanned aerial vehicle (UAV)-borne data products covering the vegetation plots surveyed during fieldwork (Kruse et al., 2021, ). The dataset includes structure-from-motion (SfM) point clouds and red-green-blue (RGB) and red-green-near-infrared (RGN) orthomosaics. From the orthomosaics, point-cloud products were created such as the digital elevation model (DEM), canopy height model (CHM), digital surface model (DSM) and the digital terrain model (DTM). The point-cloud products provide information on the three-dimensional (3D) structure of the forest at each plot. Dataset 2 contains spatial data in the form of point and polygon shapefiles of 872 individually labeled trees and shrubs that were recorded during fieldwork at the same vegetation plots (van Geffen et al., 2021c, ). The dataset contains information on tree height, crown diameter, and species type. These tree and shrub individually labeled point and polygon shapefiles were generated on top of the RGB UVA orthoimages. The individual tree information collected during the expedition such as tree height, crown diameter, and vitality are provided in table format. This dataset can be used to link individual information on trees to the location of the specific tree in the SfM point clouds, providing for example, opportunity to validate the extracted tree height from the first dataset. The dataset provides unique insights into the current state of individual trees and shrubs and allows for monitoring the effects of climate change on these individuals in the future. Dataset 3 contains a synthesis of 10 000 generated images and masks that have the tree crowns of two species of larch ( and ) automatically extracted from the RGB UAV images in the common objects in context (COCO) format (van Geffen et al., 2021a, ). As machine-learning algorithms need a large dataset to train on, the synthetic dataset was specifically created to be used for machine-learning algorithms to detect Siberian larch species. Larix gmeliniiLarix cajanderiDataset 4 contains Sentinel-2 (S-2) Level-2 bottom-of-atmosphere processed labeled image patches with seasonal information and annotated vegetation categories covering the vegetation plots (van Geffen et al., 2021b, ). The dataset is created with the aim of providing a small ready-to-use validation and training dataset to be used in various vegetation-related machine-learning tasks. It enhances the data collection as it allows classification of a larger area with the provided vegetation classes. The SiDroForest data collection serves a variety of user communities. <br /> The detailed vegetation cover and structure information in the first two datasets are of use for ecological applications, on one hand for summergreen and evergreen needle-leaf forests and also for tundra-taiga ecotones. Datasets 1 and 2 further support the generation and validation of land cover remote-sensing products in radar and optical remote sensing. In addition to providing information on forest structure and vegetation composition of the vegetation plots, the third and fourth datasets are prepared as training and validation data for machine-learning purposes. For example, the synthetic tree-crown dataset is generated from the raw UAV images and optimized to be used in neural networks. Furthermore, the fourth SiDroForest dataset contains S-2 labeled image patches processed to a high standard that provide training data on vegetation class categories for machine-learning classification with JavaScript Object Notation (JSON) labels provided. The SiDroForest data collection adds unique insights into remote hard-to-reach circumboreal forest regions.
Mineral resource exploration and mining is an essential part of today's high-tech industry. Elements such as rare-earth elements (REEs) and copper are, therefore, in high demand. Modern exploration techniques from multiple platforms (e.g., spaceborne and airborne), to detect and map the spectral characteristics of the materials of interest, require spectral libraries as an essential reference. They include field and laboratory spectral information in combination with geochemical analyses for validation. Here, we present a collection of REE- and copper-related hyperspectral spectra with associated geochemical information. The libraries contain reflectance spectra from rare-earth element oxides, REE-bearing minerals, copper-bearing minerals and mine surface samples from the Apliki copper-gold-pyrite mine in the Republic of Cyprus. The samples were measured with the HySpex imaging spectrometers in the visible and near infrared (VNIR) and shortwave infrared (SWIR) range (400-2500 nm). The geochemical validation of each sample is provided with the reflectance spectra. The spectral libraries are openly available to assist future mineral mapping campaigns and laboratory spectroscopic analyses. The spectral libraries and corresponding geochemistry are published via GFZ Data Services with the following DOIs: https://doi.org/10.5880/GFZ.1.4.2019.004 (13 REE-bearing minerals and 16 oxide powders, Koerting et al., 2019a), https://doi.org/10.5880/GFZ.1.4.2019.003 (20 copper-bearing minerals, Koellner et al., 2019), and https://doi.org/10.5880/GFZ.1.4.2019.005 (37 copper-bearing surface material samples from the Apliki coppergold-pyrite mine in Cyprus, Koerting et al., 2019b). All spectral libraries are united and comparable by the internally consistent method of hyperspectral data acquisition in the laboratory.
Tectonically active coasts are dynamic environments characterized by the presence of multiple marine terraces formed by the combined effects of wave erosion, tectonic uplift, and sea-level oscillations at glacialcycle timescales. Well-preserved erosional terraces from the last interglacial sea-level highstand are ideal marker horizons for reconstructing past sea-level positions and calculating vertical displacement rates. We carried out an almost continuous mapping of the last interglacial marine terrace along similar to 5000 km of the western coast of South America between 1 degrees N and 40 degrees S. We used quantitatively replicable approaches constrained by published terrace-age estimates to ultimately compare elevations and patterns of uplifted terraces with tectonic and climatic parameters in order to evaluate the controlling mechanisms for the formation and preservation of marine terraces and crustal deformation. Uncertainties were estimated on the basis of measurement errors and the distance from referencing points. Overall, our results indicate a median elevation of 30.1 m, which would imply a median uplift rate of 0.22 m kyr(-1) averaged over the past similar to 125 kyr. The patterns of terrace elevation and uplift rate display high-amplitude (similar to 100-200 m) and long-wavelength (similar to 10(2) km) structures at the Manta Peninsula (Ecuador), the San Juan de Marcona area (central Peru), and the Arauco Peninsula (south-central Chile). Medium-wavelength structures occur at the Mejillones Peninsula and Topocalma in Chile, while short-wavelength (< 10 km) features are for instance located near Los Vilos, Valparaiso, and Carranza, Chile. We interpret the long-wavelength deformation to be controlled by deep-seated processes at the plate interface such as the subduction of major bathymetric anomalies like the Nazca and Carnegie ridges. In contrast, short-wavelength deformation may be primarily controlled by sources in the upper plate such as crustal faulting, which, however, may also be associated with the subduction of topographically less pronounced bathymetric anomalies. Latitudinal differences in climate additionally control the formation and preservation of marine terraces. Based on our synopsis we propose that increasing wave height and tidal range result in enhanced erosion and morphologically well-defined marine terraces in south-central Chile. Our study emphasizes the importance of using systematic measurements and uniform, quantitative methodologies to characterize and correctly interpret marine terraces at regional scales, especially if they are used to unravel the tectonic and climatic forcing mechanisms of their formation. This database is an integral part of the World Atlas of Last Interglacial Shorelines (WALIS), published online at https://doi.org/10.5281/zenodo.4309748 (Freisleben et al., 2020).
Nature-based solutions (NBS) are seen as a promising adaptation measure that sustainably deals with diverse societal challenges, while simultaneously delivering multiple benefits. Nature-based solutions have been highlighted as a resilient and sustainable means of mitigating floods and other hazards globally. This study examined diverging conceptualizations of NBS, as well as the attitudinal (for example, emotions and beliefs) and contextual (for example, legal and political aspects) barriers and drivers of NBS for flood risks in South Korea. Semistructured interviews were conducted with 11 experts and focused on the topic of flood risk measures and NBS case studies. The analysis found 11 barriers and five drivers in the attitudinal domain, and 13 barriers and two drivers in the contextual domain. Most experts see direct monetary benefits as an important attitudinal factor for the public. Meanwhile, the cost-effectiveness of NBS and their capacity to cope with flood risks were deemed influential factors that could lead decision makers to opt for NBS. Among the contextual factors, insufficient systems to integrate NBS in practice and the ideologicalization of NBS policy were found to be peculiar barriers, which hinder consistent realization of initiatives and a long-term national plan for NBS. Understanding the barriers and drivers related to the mainstreaming of NBS is critical if we are to make the most of such solutions for society and nature. It is also essential that we have a shared definition, expectation, and vision of NBS.
The LArge-scale Reservoir Simulator (LARS) has been previously developed to study hydrate dissociation in hydrate-bearing systems under in-situ conditions. In the present study, a numerical framework of equations of state describing hydrate formation at equilibrium conditions has been elaborated and integrated with a numerical flow and transport simulator to investigate a multi-stage hydrate formation experiment undertaken in LARS. A verification of the implemented modeling framework has been carried out by benchmarking it against another established numerical code. Three-dimensional (3D) model calibration has been performed based on laboratory data available from temperature sensors, fluid sampling, and electrical resistivity tomography. The simulation results demonstrate that temperature profiles, spatial hydrate distribution, and bulk hydrate saturation are consistent with the observations. Furthermore, our numerical framework can be applied to calibrate geophysical measurements, optimize post-processing workflows for monitoring data, improve the design of hydrate formation experiments, and investigate the temporal evolution of sub-permafrost methane hydrate reservoirs.
Transport properties of potential host rocks for nuclear waste disposal are typically determined in laboratory or in-situ experiments under geochemically controlled and constant conditions. Such a homogeneous assumption is no longer applicable on the host rock scale as can be seen from the pore water profiles of the potential host rock Opalinus Clay at Mont Terri (Switzerland). The embedding aquifers are the hydro-geological boundaries, that established gradients in the 210 m thick low permeable section through diffusive exchange over millions of years. Present-day pore water profiles were confirmed by a data-driven as well as by a conceptual scenario. Based on the modelled profiles, the influence of the geochemical gradient on uranium migration was quantified by comparing the distances after one million years with results of common homogeneous models. Considering the heterogeneous system, uranium migrated up to 24 m farther through the formation depending on the source term position within the gradient and on the partial pressure of carbon dioxide pCO2 of the system. Migration lengths were almost equal for single- and multicomponent diffusion. Differences can predominantly be attributed to changes in the sorption capacity, whereby pCO2 governs how strong uranium migration is affected by the geochemical gradient. Thus, the governing parameters for uranium migration in the Opalinus Clay can be ordered in descending priority: pCO2, geochemical gradients, mineralogical heterogeneity.</p>
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.
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.
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.
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.
During the period 750-600 Ma ago, prior to the final break-up of the supercontinent Rodinia, the crust of both the North American Craton and Baltica was intruded by significant amounts of rift-related magmas originating from the mantle. In the Proterozoic crust of Southern Norway, the 580 Ma old Fen carbonatite-ultramafic complex is a representative of this type of rocks. In this paper, we report the occurrence of an ultramafic lamprophyre dyke which possibly is linked to the Fen complex, although Ar-40/Ar-39 data from phenocrystic phlogopite from the dyke gave an age of 686 +/- 9 Ma. The lamprophyre dyke was recently discovered in one of the Kongsberg silver mines at Vinoren, Norway. Whole rock geochemistry, geochronological and mineralogical data from the ultramafic lamprophyre dyke are presented aiming to elucidate its origin and possible geodynamic setting. From the whole-rock composition of the Vinoren dyke, the rock could be recognized as transitional between carbonatite and kimberlite-II (orangeite). From its diagnostic mineralogy, the rock is classified as aillikite. The compositions and xenocrystic nature of several of the major and accessory minerals from the Vinoren aillikite are characteristic for diamondiferous rocks (kimberlites/lamproites/UML): Phlogopite with kinoshitalite-rich rims, chromite-spinel-ulvospinel series, Mg- and Mn-rich ilmenites, rutile and lucasite-(Ce). We suggest that the aillikite melt formed during partial melting of a MARID (mica-amphibole-rutile-ilmenite-diopside)-like source under CO2 fluxing. The pre-rifting geodynamic setting of the Vinoren aillikite before the Rodinia supercontinent breakup suggests a relatively thick SCLM (Subcontinental Lithospheric Mantle) during this stage and might indicate a diamond-bearing source for the parental melt. This is in contrast to the about 100 Ma younger Fen complex, which were derived from a thin SCLM.
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.
How insoluble inclusions and intersecting layers affect the leaching process within potash seams
(2021)
Potash seams are a valuable resource containing several economically interesting, but also highly soluble minerals. In the presence of water, uncontrolled leaching can occur, endangering subsurface mining operations. In the present study, the influence of insoluble inclusions and intersecting layers on leaching zone evolution was examined by means of a reactive transport model. For that purpose, a scenario analysis was carried out, considering different rock distributions within a carnallite-bearing potash seam. The results show that reaction-dominated systems are not affected by heterogeneities at all, whereas transport-dominated systems exhibit a faster advance in homogeneous rock compositions. In return, the ratio of permeated rock in vertical direction is higher in heterogeneous systems. Literature data indicate that most natural potash systems are transport-dominated. Accordingly, insoluble inclusions and intersecting layers can usually be seen as beneficial with regard to reducing hazard potential as long as the mechanical stability of leaching zones is maintained. Thereby, the distribution of insoluble areas is of minor impact unless an inclined, intersecting layer occurs that accelerates leaching zone growth in one direction. Moreover, it is found that the saturation dependency of dissolution rates increases the growth rate in the long term, and therefore must be considered in risk assessments.
Monitoring the response of volcanic CO2 emissions to changes in the Los Humeros hydrothermal system
(2021)
Carbon dioxide is the most abundant, non-condensable gas in volcanic systems, released into the atmosphere through either diffuse or advective fluid flow. The emission of substantial amounts of CO2 at Earth's surface is not only controlled by volcanic plumes during periods of eruptive activity or fumaroles, but also by soil degassing along permeable structures in the subsurface. Monitoring of these processes is of utmost importance for volcanic hazard analyses, and is also relevant for managing geothermal resources. Fluid-bearing faults are key elements of economic value for geothermal power generation. Here, we describe for the first time how sensitively and quickly natural gas emissions react to changes within a deep hydrothermal system due to geothermal fluid reinjection. For this purpose, we deployed an automated, multi-chamber CO2 flux monitoring system within the damage zone of a deep-rooted major normal fault in the Los Humeros Volcanic Complex (LHVC) in Mexico and recorded data over a period of five months. After removing the atmospheric effects on variations in CO2 flux, we calculated correlation coefficients between residual CO2 emissions and reinjection rates, identifying an inverse correlation of rho = - 0.51 to - 0.66. Our results indicate that gas emissions respond to changes in reinjection rates within 24 h, proving an active hydraulic communication between the hydrothermal system and Earth's surface. This finding is a promising indication not only for geothermal reservoir monitoring but also for advanced long-term volcanic risk analysis. Response times allow for estimation of fluid migration velocities, which is a key constraint for conceptual and numerical modelling of fluid flow in fracture-dominated systems.
POET (v0.1): speedup of many-core parallel reactive transport simulations with fast DHT lookups
(2021)
Coupled reactive transport simulations are extremely demanding in terms of required computational power, which hampers their application and leads to coarsened and oversimplified domains. The chemical sub-process represents the major bottleneck: its acceleration is an urgent challenge which gathers increasing interdisciplinary interest along with pressing requirements for subsurface utilization such as spent nuclear fuel storage, geothermal energy and CO2 storage. In this context we developed POET (POtsdam rEactive Transport), a research parallel reactive transport simulator integrating algorithmic improvements which decisively speed up coupled simulations. In particular, POET is designed with a master/worker architecture, which ensures computational efficiency in both multicore and cluster compute environments. POET does not rely on contiguous grid partitions for the parallelization of chemistry but forms work packages composed of grid cells distant from each other. Such scattering prevents particularly expensive geochemical simulations, usually concentrated in the vicinity of a reactive front, from generating load imbalance between the available CPUs (central processing units), as is often the case with classical partitions. Furthermore, POET leverages an original implementation of the distributed hash table (DHT) mechanism to cache the results of geochemical simulations for further reuse in subsequent time steps during the coupled simulation. The caching is hence particularly advantageous for initially chemically homogeneous simulations and for smooth reaction fronts. We tune the rounding employed in the DHT on a 2D benchmark to validate the caching approach, and we evaluate the performance gain of POET's master/worker architecture and the DHT speedup on a 3D benchmark comprising around 650 000 grid elements. The runtime for 200 coupling iterations, corresponding to 960 simulation days, reduced from about 24 h on 11 workers to 29 min on 719 workers. Activating the DHT reduces the runtime further to 2 h and 8 min respectively. Only with these kinds of reduced hardware requirements and computational costs is it possible to realistically perform the longterm complex reactive transport simulations, as well as perform the uncertainty analyses required by pressing societal challenges connected with subsurface utilization.
Compound natural hazards likeEl Ninoevents cause high damage to society, which to manage requires reliable risk assessments. Damage modelling is a prerequisite for quantitative risk estimations, yet many procedures still rely on expert knowledge, and empirical studies investigating damage from compound natural hazards hardly exist. A nationwide building survey in Peru after theEl Ninoevent 2017 - which caused intense rainfall, ponding water, flash floods and landslides - enables us to apply data-mining methods for statistical groundwork, using explanatory features generated from remote sensing products and open data. We separate regions of different dominant characteristics through unsupervised clustering, and investigate feature importance rankings for classifying damage via supervised machine learning. Besides the expected effect of precipitation, the classification algorithms select the topographic wetness index as most important feature, especially in low elevation areas. The slope length and steepness factor ranks high for mountains and canyons. Partial dependence plots further hint at amplified vulnerability in rural areas. An example of an empirical damage probability map, developed with a random forest model, is provided to demonstrate the technical feasibility.
The global geological volatile cycle (H, C, N) plays an important role in the long term self-regulation of the Earth system. However, the complex interaction between its deep, solid Earth components (i.e. crust and mantle), Earth's fluid envelopes (i.e. atmosphere and hydrosphere) and plate tectonic processes is a subject of ongoing debate. In this study we want to draw attention to how the presence of primary melt (MI) and fluid (FI) inclusions in high-grade metamorphic minerals could help constrain the crustal component of the volatile cycle. To that end, we review the distribution of MI and FI throughout Earth's history, from ca. 3.0 Ga ago up to the present day. We argue that the lower crust might constitute an important, long-term, volatile storage unit, capable to influence the composition of the surface envelopes through the mean of weathering, crustal thickening, partial melting and crustal assimilation during volcanic activity. Combined with thermodynamic modelling, our compilation indicates that periods of well-established plate tectonic regimes at <0.85 Ga and 1.7-2.1 Ga, might be more prone to the reworking of supracrustal lithologies and the storage of volatiles in the lower crust. Such hypothesis has implication beyond the scope of metamorphic petrology as it potentially links geodynamic mechanisms to habitable surface conditions. MI and FI in metamorphic crustal rocks then represent an invaluable archive to assess and quantify the co-joint evolution of plate tectonics and Earth's external processes. (C) 2021 China University of Geosciences (Beijing) and Peking University. Production and hosting by Elsevier B.V. This is an open access article under the CC BY-NC-ND license (http://creativecommons.org/licenses/by-nc-nd/4.0/).
Multi-component (MC) diffusion simulations enable a process based and more precise approach to calculate transport and sorption compared to the commonly used single-component (SC) models following Fick's law. The MC approach takes into account the interaction of chemical species in the porewater with the diffuse double layer (DDL) adhering clay mineral surfaces. We studied the shaly, sandy and carbonate-rich facies of the Opalinus Clay. High clay contents dominate diffusion and sorption of uranium. The MC simulations show shorter diffusion lengths than the SC models due to anion exclusion from the DDL. This hampers diffusion of the predominant species CaUO2(CO3)32-. On the one side, species concentrations and ionic strengths of the porewater and on the other side surface charge of the clay minerals control the composition and behaviour of the DDL. For some instances, it amplifies the diffusion of uranium. We developed a workflow to transfer computationally intensive MC simulations to SC models via calibrated effective diffusion and distribution coefficients. Simulations for one million years depict maximum uranium diffusion lengths between 10 m and 35 m. With respect to the minimum requirement of a thickness of 100 m, the Opalinus Clay seems to be a suitable host rock for nuclear waste repositories.
The occurrence of refugia beyond the arctic treeline and genetic adaptation therein play a crucial role of largely unknown effect size. While refugia have potential for rapidly colonizing the tundra under global warming, the taxa may be maladapted to the new environmental conditions. Understanding the genetic composition and age of refugia is thus crucial for predicting any migration response.
Here, we genotype 194 larch individuals from an similar to 1.8 km(2)area in northcentral Siberia on the southern Taimyr Peninsula by applying an assay of 16 nuclear microsatellite markers. For estimating the age of clonal individuals, we counted tree rings at sections along branches to establish a lateral growth rate that was then combined with geographic distance.
Findings reveal that the predominant reproduction type is clonal (58.76%) by short distance spreading of ramets. One outlier of clones 1 km apart could have been dispersed by reindeer. In clonal groups and within individuals, we find that somatic mutations accumulate with geographic distance. Clonal groups of two or more individuals are observed. Clonal age estimates regularly suggest individuals as old as 2,200 years, which coincides with a major environmental change that forced a treeline retreat in the region.
We conclude that individuals with clonal growth mode were naturally selected as it lowers the likely risk of extinction under a harsh environment. We discuss this legacy from the past that might now be a maladaptation and hinder expansion under currently strongly increasing temperatures.
The Issyk-Kul Basin (Kyrgyzstan), situated in the central Tian Shan Mountains, hosts the largest and deepest mountain lake in Central Asia. Erosion of the surrounding Terskey and Kungey ranges led to the accumulation of up to 4 km of sediment in the adjacent depression. Creation of the basin from regional shortening and uplift likely initiated around the Oligocene-Miocene, yet precise age control is sparse. To better understand the timing of these processes, we obtained magnetostratigraphic age constraints on fossil-poor, fluvio-lacustrine sediments exposed south of Lake Issyk-Kul, that agree well with previous age constraints of the equivalent strata outside the Issyk-Kul Basin. Two 500-650 m thick sections comprised mainly of Chu Group sediments were dated at 6.3-2.8 Ma and 7.0-2.4 Ma (late Miocene to early Pleistocene). Together with reinterpreted magnetostratigraphic constraints from underlying strata, we find that syn-tectonic deposition commenced at similar to 22 Ma with average sedimentation rates <10 cm/ka. Sedimentation rates increased to 10-30 cm/ka at 7 Ma, concurrent with accelerated uplift in the Terskey Range to the south. A deformation event in one section (Kaji-Say) between 5 and 3 Ma together with concurrent shifts of depositional centers throughout the basin signal the onset of substantial uplift of the Kungey Range to the north at similar to 5 Ma. This uplift and deformation transformed the Issyk-Kul area into a closed basin that facilitated the formation of a deep lake. Lacustrine facies deposited around 3 Ma mark the existence of Lake Issyk-Kul by that time.
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.
Diagenetic trends of synthetic reservoir sandstone properties assessed by digital rock physics
(2021)
Quantifying interactions and dependencies among geometric, hydraulic and mechanical properties of reservoir sandstones is of particular importance for the exploration and utilisation of the geological subsurface and can be assessed by synthetic sandstones comprising the microstructural complexity of natural rocks. In the present study, three highly resolved samples of the Fontainebleau, Berea and Bentheim sandstones are generated by means of a process-based approach, which combines the gravity-driven deposition of irregularly shaped grains and their diagenetic cementation by three different schemes. The resulting evolution in porosity, permeability and rock stiffness is examined and compared to the respective micro-computer tomographic (micro-CT) scans. The grain contact-preferential scheme implies a progressive clogging of small throats and consequently produces considerably less connected and stiffer samples than the two other schemes. By contrast, uniform quartz overgrowth continuously alters the pore space and leads to the lowest elastic properties. The proposed stress-dependent cementation scheme combines both approaches of contact-cement and quartz overgrowth, resulting in granulometric, hydraulic and elastic properties equivalent to those of the respective micro-CT scans, where bulk moduli slightly deviate by 0.8%, 4.9% and 2.5% for the Fontainebleau, Berea and Bentheim sandstone, respectively. The synthetic samples can be further altered to examine the impact of mineral dissolution or precipitation as well as fracturing on various petrophysical correlations, which is of particular relevance for numerous aspects of a sustainable subsurface utilisation.
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.
The imagination of clearly separated core-shell structures is already outdated by the fact, that the nanoparticle core-shell structures remain in terms of efficiency behind their respective bulk material due to intermixing between core and shell dopant ions. In order to optimize the photoluminescence of core-shell UCNP the intermixing should be as small as possible and therefore, key parameters of this process need to be identified. In the present work the Ln(III) ion migration in the host lattices NaYF4 and NaGdF4 was monitored. These investigations have been performed by laser spectroscopy with help of lanthanide resonance energy transfer (LRET) between Eu(III) as donor and Pr(III) or Nd(III) as acceptor. The LRET is evaluated based on the Forster theory. The findings corroborate the literature and point out the migration of ions in the host lattices. Based on the introduced LRET model, the acceptor concentration in the surrounding of one donor depends clearly on the design of the applied core-shell-shell nanoparticles. In general, thinner intermediate insulating shells lead to higher acceptor concentration, stronger quenching of the Eu(III) donor and subsequently stronger sensitization of the Pr(III) or the Nd(III) acceptors. The choice of the host lattice as well as of the synthesis temperature are parameters to be considered for the intermixing process.
The computational costs associated with coupled reactive transport simulations are mostly due to the chemical subsystem: replacing it with a pre-trained statistical surrogate is a promising strategy to achieve decisive speedups at the price of small accuracy losses and thus to extend the scale of problems which can be handled. We introduce a hierarchical coupling scheme in which "full-physics" equation-based geochemical simulations are partially replaced by surrogates. Errors in mass balance resulting from multivariate surrogate predictions effectively assess the accuracy of multivariate regressions at runtime: inaccurate surrogate predictions are rejected and the more expensive equation-based simulations are run instead. Gradient boosting regressors such as XGBoost, not requiring data standardization and being able to handle Tweedie distributions, proved to be a suitable emulator. Finally, we devise a surrogate approach based on geochemical knowledge, which overcomes the issue of robustness when encountering previously unseen data and which can serve as a basis for further development of hybrid physics-AI modelling.
Projected changes in compound flood hazard from riverine and coastal floods in northwestern Europe
(2020)
Compound flooding in coastal regions, that is, the simultaneous or successive occurrence of high sea levels and high river flows, is expected to increase in a warmer world. To date, however, there is no robust evidence on projected changes in compound flooding for northwestern Europe. We combine projected storm surges and river floods with probabilistic, localized relative sea-level rise (SLR) scenarios to assess the future compound flood hazard over northwestern coastal Europe in the high (RCP8.5) emission scenario. We use high-resolution, dynamically downscaled regional climate models (RCM) to drive a storm surge model and a hydrological model, and analyze the joint occurrence of high coastal water levels and associated river peaks in a multivariate copula-based approach. The RCM-forced multimodel mean reasonably represents the observed spatial pattern of the dependence strength between annual maxima surge and peak river discharge, although substantial discrepancies exist between observed and simulated dependence strength. All models overestimate the dependence strength, possibly due to limitations in model parameterizations. This bias affects compound flood hazard estimates and requires further investigation. While our results suggest decreasing compound flood hazard over the majority of sites by 2050s (2040-2069) compared to the reference period (1985-2005), an increase in projected compound flood hazard is limited to around 34% of the sites. Further, we show the substantial role of SLR, a driver of compound floods, which has frequently been neglected. Our findings highlight the need to be aware of the limitations of the current generation of Earth system models in simulating coastal compound floods.
Natural gas hydrate occurrences contain predominantly methane; however, there are increasing reports of complex mixed gas hydrates and coexisting hydrate phases. Changes in the feed gas composition due to the preferred incorporation of certain components into the hydrate phase and an inadequate gas supply is often assumed to be the cause of coexisting hydrate phases. This could also be the case for the gas hydrate system in Qilian Mountain permafrost (QMP), which is mainly controlled by pores and fractures with complex gas compositions. This study is dedicated to the experimental investigations on the formation process of mixed gas hydrates based on the reservoir conditions in QMP. Hydrates were synthesized from water and a gas mixture under different gas supply conditions to study the effects on the hydrate formation process. In situ Raman spectroscopic measurements and microscopic observations were applied to record changes in both gas and hydrate phase over the whole formation process. The results demonstrated the effects of gas flow on the composition of the resulting hydrate phase, indicating a competitive enclathration of guest molecules into the hydrate lattice depending on their properties. Another observation was that despite significant changes in the gas composition, no coexisting hydrate phases were formed.
Multi-year index-based insurance for adapting Water Utility Companies to hydrological drought
(2020)
The sustainability of water utility companies is threatened by non-stationary drivers, such as climate and anthropogenic changes. To cope with potential economic losses, instruments such as insurance are useful for planning scenarios and mitigating impacts, but data limitations and risk uncertainties affect premium estimation and, consequently, business sustainability. This research estimated the possible economic impacts of business interruption to the Sao Paulo Water Utility Company derived from hydrological drought and how this could be mitigated with an insurance scheme. Multi-year insurance (MYI) was proposed through a set of "change" drivers: the climate driver, through forcing the water evaluation and planning system (WEAP) hydrological tool; the anthropogenic driver, through water demand projections; and the economic driver, associated with recent water price policies adopted by the utility company during water scarcity periods. In our study case, the evaluated indices showed that MYI contracts that cover only longer droughts, regardless of the magnitude, offer better financial performance than contracts that cover all events (in terms of drought duration). Moreover, through MYI contracts, we demonstrate solvency for the insurance fund in the long term and an annual average actuarially fair premium close to the total expected revenue reduction.
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.
The Fram Strait is an area with a relatively low and irregular distribution of diatom microfossils in surface sediments, and thus microfossil records are scarce, rarely exceed the Holocene, and contain sparse information about past richness and taxonomic composition. These attributes make the Fram Strait an ideal study site to test the utility of sedimentary ancient DNA (sedaDNA) metabarcoding. Amplifying a short, partial rbcL marker from samples of sediment core MSM05/5-712-2 resulted in 95.7% of our sequences being assigned to diatoms across 18 different families, with 38.6% of them being resolved to species and 25.8% to genus level. Independent replicates show a high similarity of PCR products, especially in the oldest samples. Diatom sedaDNA richness is highest in the Late Weichselian and lowest in Mid- and Late Holocene samples. Taxonomic composition is dominated by cold-water and sea-ice-associated diatoms and suggests several reorganisations - after the Last Glacial Maximum, after the Younger Dryas, and after the Early and after the Mid-Holocene. Different sequences assigned to, amongst others, Chaetoceros socialis indicate the detectability of intra-specific diversity using sedaDNA. We detect no clear pattern between our diatom sedaDNA record and the previously published IP25 record of this core, although proportions of pennate diatoms increase with higher IP25 concentrations and proportions of Nitzschia cf. frigida exceeding 2% of the assemblage point towards past sea-ice presence.
The Victoria microplate between the Eastern and Western Branches of the East African Rift System is one of the largest continental microplates on Earth. In striking contrast to its neighboring plates, Victoria rotates counterclockwise with respect to Nubia. The underlying cause of this distinctive rotation has remained elusive so far. Using 3D numerical models, we investigate the role of pre-existing lithospheric heterogeneities in continental microplate rotation. We find that Victoria's rotation is primarily controlled by the distribution of rheologically stronger zones that transmit the drag of the major plates to the microplate and of the mechanically weaker mobile belts surrounding Victoria that facilitate rotation. Our models reproduce Victoria's GPS-derived counterclockwise rotation as well as key complexities of the regional tectonic stress field. These results reconcile competing ideas on the opening of the rift system by highlighting differences in orientation of the far-field divergence, local extension, and the minimum horizontal stress. One of the largest continental microplates on Earth is situated in the center of the East African Rift System, and oddly, the Victoria microplate rotates counterclockwise with respect to the neighboring African tectonic plate. Here, the authors' modelling results suggest that Victoria microplate rotation is caused by edge-driven lithospheric processes related to the specific geometry of rheologically weak and strong regions.