Refine
Year of publication
- 2019 (276) (remove)
Document Type
- Article (216)
- Doctoral Thesis (24)
- Postprint (15)
- Other (9)
- Review (7)
- Habilitation Thesis (3)
- Conference Proceeding (1)
- Master's Thesis (1)
Is part of the Bibliography
- yes (276)
Keywords
- climate (5)
- remote sensing (4)
- Andes (3)
- Climate (3)
- Holocene (3)
- geochronology (3)
- melt inclusions (3)
- paleoclimate (3)
- photogrammetry (3)
- silicon (3)
- AMOC (2)
- Accuracy Asseessment (2)
- Anden (2)
- Asia (2)
- Baldeggersee (2)
- Bayesianism (2)
- CU (2)
- Central Andes (2)
- Climate change (2)
- Climate change adaptation (2)
- Connectivity (2)
- Deformation (2)
- Digital Elevation Model (2)
- Earthquake (2)
- Error (2)
- Expression (2)
- Extraction (2)
- Flood risk (2)
- Forecasting Framework (2)
- Geomechanical Model (2)
- Hochwasser (2)
- InSAR (2)
- Jaguaribe Basin (2)
- Landscape Response (2)
- Long-Term (2)
- MO (2)
- Magnetotellurics (2)
- Middle Strand (2)
- NW Argentina (2)
- NW Turkey (2)
- Nordeste (2)
- North Anatolian Fault (2)
- Northwestern Anatolia (2)
- OpenForecast (2)
- Paleoclimate (2)
- Paleogene (2)
- Paleolimnology (2)
- Pamir (2)
- Patterns (2)
- Permafrost (2)
- Permian (2)
- Pleistocene (2)
- Pollen (2)
- Precipitation (2)
- Pull-Apart (2)
- Pyrenees (2)
- Recurrence plot (2)
- Reservoir Networks (2)
- River Incision Model (2)
- Russia (2)
- SIO₂ (2)
- Sediment Transport (2)
- Shear Zone (2)
- Siberia (2)
- Site effects (2)
- Slip Distribution (2)
- Slope (2)
- Subjective well-being (2)
- Tectonic Evolution (2)
- Time-series analysis (2)
- UAV (2)
- Uncertainties (2)
- Uncertainty Processor (2)
- Validation (2)
- Water Availability (2)
- alps (2)
- carbon dioxide (2)
- climate change adaptation (2)
- clinopyroxenite (2)
- complexes (2)
- convective available potential energy (2)
- correlation (2)
- crystalline (2)
- data science (2)
- deep convection (2)
- delta-c-13 (2)
- detrital carbonate (2)
- dew-point temperature (2)
- eastern south–central Andes (2)
- eclogite (2)
- extreme rainfall (2)
- flood risk (2)
- fluid inclusions (2)
- forecasting (2)
- forest structure (2)
- geochemistry (2)
- geomagnetic observatories (2)
- geomagnetic observatory data (2)
- geomagnetic storm drivers (2)
- glass (2)
- historical geomagnetic storms (2)
- hydro-meterological hazards (2)
- hydrogeology (2)
- hydrology (2)
- induced seismicity (2)
- land cover change (2)
- magnetostratigraphy (2)
- metasomatism (2)
- methane (2)
- modelling (2)
- northern high latitudes (2)
- numerische Modellierung (2)
- object-based damage modeling (2)
- open (2)
- operational service (2)
- opinion (2)
- optical data (2)
- orogenic peridotite (2)
- oxygen (2)
- paleomagnetism (2)
- participation (2)
- permafrost (2)
- permafrost disturbances (2)
- point cloud (2)
- precipitation (2)
- pressures (2)
- prior derivation (2)
- probabilistic approaches (2)
- record (2)
- risk assessment (2)
- runoff (2)
- rutile-type (2)
- sedimentology (2)
- spatial scales (2)
- spectral adjustment (2)
- stable isotopes (2)
- structure from motion (2)
- supervised machine learning (2)
- surface reflectance (2)
- tectonic geomorphology (2)
- time series (2)
- tin (2)
- uncertainty (2)
- uncertainty analysis (2)
- varved sediments (2)
- wave-particle interactions (2)
- westerlies (2)
- 2-D numerical model (1)
- 2D numerical experiments (1)
- Abrupt events (1)
- Abschiebungshorizonte (1)
- Absorption feature parameters (1)
- Adaptation (1)
- Adaptation process (1)
- Agent-based models (1)
- Alterationsgeochemie (1)
- Ambient noise (1)
- Amplification (1)
- Anatolia (1)
- Ancient Gneiss Complex (1)
- Andean Plateau (1)
- Andean retroarc (1)
- Angewandte Geophysik (1)
- Annual laminae (1)
- Apennine Carbonate Platform (ACP) (1)
- Applied Geophysics (1)
- Archaeolithoporella (1)
- Arctic tundra (1)
- Argentina (1)
- Argentinean broken foreland (1)
- Argentinien (1)
- Array seismology (1)
- Asian monsoon (1)
- Assimilation-fractional crystallization (1)
- Asymmetric rifting (1)
- BNDVI (1)
- Back-arc basin (1)
- Baiu (1)
- Basin analysis (1)
- Bayesian hierarchical model (1)
- Bayesian modeling (1)
- Be-10 (1)
- Beacon fire (1)
- Beringia (1)
- Binalud Mountains (1)
- Biomisation (1)
- Black Sea (1)
- Blueschist (1)
- Blueschist metamorphism (1)
- Bohrloch-Rekonstruktion (1)
- Bolivian tin belt (1)
- Bor (1)
- Bor-isotopen (1)
- Borehole reconstruction (1)
- Boron isotopes (1)
- Brazil (1)
- Britholite (1)
- British Columbia (1)
- Burned areas (1)
- Bykovsky Peninsula (1)
- C isotopes (1)
- C-14 bomb peak (1)
- Carbon (1)
- Carbon 14 (1)
- Carbon stable isotopes (1)
- Carbonate contourite drift (1)
- Catchment geomorphology (1)
- Catchment order (1)
- Cenozoic (1)
- Cenozoic basins (1)
- Central Asia (1)
- Central Chile (1)
- Central Myanmar Basin (1)
- Chaco-Paraná Becken (1)
- Chaco-Paraná basin (1)
- Chemeron Formation (1)
- Chile (1)
- City population (1)
- Clay mineralogy (1)
- Climate Change (1)
- Climate prediction (1)
- Climate warming (1)
- Common garden experiments (1)
- Complex networks (1)
- Compound dislocation models (CDMs) (1)
- Conservation planning (1)
- Controlled source electromagnetics (CSEM) (1)
- Coupled heat and mass transport (1)
- Cross-dating (1)
- CryoGRID (1)
- Cryospheric science (1)
- Cryostratigraphy (1)
- Cumbres Calchaquies (1)
- Cyanobacteria (1)
- Cyclic soft stimulation (CSS) (1)
- DDM simulation (1)
- DFT (1)
- Danube (1)
- Data assimilation (1)
- Data fusion (1)
- Databases (1)
- Deep carbon cycle (1)
- Denitrification (1)
- Density modelling (1)
- Dichtemodellierung (1)
- Diffusion (1)
- Discrete Element Method (1)
- Diskrete-Elemente-Methode (1)
- Dryland (1)
- Dämpfungstomographie (1)
- EU Floods Directive (1)
- Early Eocene Climatic Optimum (1)
- Early cretaceous (1)
- Early warning system (1)
- Earth in Five Reactions: A Deep Carbon Perspective (1)
- Earth tides (1)
- Earthquake hazards (1)
- East Africa (1)
- Eastern Africa (1)
- Eastern Himalaya (1)
- Eastern Sierras Pampeanas (1)
- Eastern Tibet (1)
- Eastern Tibetan Plateau (1)
- Eccentricity (1)
- Eclogite (1)
- Ecogeomorphology (1)
- Ecosystem dissimilarity (1)
- Ecosystem dynamics (1)
- Ecosystem services (1)
- Elastizitätsmodul (1)
- Electrical conductivity (1)
- Electromagnetic theory (1)
- Emission scenarios (1)
- Environmental (1)
- Environmental impact (1)
- Environmental magnetism (1)
- Eocene (1)
- Epithermal Ag-Sn deposits (1)
- Erdbebenmodelierung (1)
- Erdbebenquellinversion (1)
- Erdbebenschwärme (1)
- Erdfälle (1)
- Erodium cicutarium (1)
- Error propagation (1)
- Europe (1)
- European air quality (1)
- Evaporitic varves (1)
- Event synchronization (1)
- Exhumierung (1)
- Extreme habitats (1)
- Extreme rainfall (1)
- Extremophiles (1)
- FAIR (1)
- FISH (1)
- Falten (1)
- Feature selection (1)
- Federal Water Act (1)
- Festigkeit (1)
- Field analogue (1)
- Flood (1)
- Flooding (1)
- Fluid inclusions (1)
- Fluid processes (1)
- Fluid-Gesteins Wechselwirkung (1)
- Fluid-strömungen (1)
- Flüssigkeitseinschlüsse (1)
- Flüssigkeitsinklusionen (1)
- Food security (1)
- Foreland basin (1)
- Foreland basin fragmentation (1)
- Foreland basins (1)
- GNSS Reflectometry (1)
- GST (1)
- Gekoppelter Wärme- und Massetransport (1)
- Generalized additive models (1)
- Geochemistry (1)
- Geochronologie (1)
- Geodynamic modelling (1)
- Geomechanical Modelling (1)
- Geomechanische Modellierung (1)
- Geomorphology (1)
- Geothermobarometry (1)
- Germany (1)
- Glacial lake outbursts (1)
- Glacial refugia (1)
- Glacierized basins (1)
- Glasstruktur (1)
- Gravity modelling (1)
- Ground-penetrating radar (1)
- Groundwater quality (1)
- Gypsum quantification (1)
- Gyttja (1)
- H/V Verhältnis (1)
- H/V ratio (1)
- Haida Gwaii (1)
- Heat flow (1)
- Herschel Island Qikiqtaruk (1)
- High pressure (1)
- Hillslope thermokarst (1)
- Himalaja (1)
- Himalaya (1)
- Himalayan hydroclimate (1)
- Hindu Kush-Karakoram (1)
- Hochwasserrisikomanagementrichtlinie (1)
- Hutchinson niche (1)
- Hydraulic fracturing (1)
- Hydro-mechanical coupling (1)
- Hydrometric networks (1)
- Hypersaline lake (1)
- Hyperspectral (1)
- Hypolimnetic oxygen (1)
- Hämmerlein (1)
- ICDP Dead Sea deep drilling (1)
- Ice wedges (1)
- In-situ degradation rates (1)
- India-Asia collision (1)
- Indian summer monsoon (1)
- Induced seismicity (1)
- Informative prior (1)
- Infrared (1)
- Insurance (1)
- Intangible losses (1)
- International unions (1)
- Inundation (1)
- Inverse modelling (1)
- Inverse theory (1)
- Irrawaddy River (1)
- Irrigation (1)
- Isotope-hydrological integrated modeling (1)
- Isotopen (1)
- Isotopes (1)
- Isótopos de Boro (1)
- Italy (1)
- Ivrea Zone (1)
- Kameng River section (1)
- Karst (1)
- Kenya Rift (1)
- KiK-net (1)
- Klima (1)
- Kohlenstoffdioxid (1)
- Kp index (1)
- Kriecheigenschaften (1)
- Krustenverformungen (1)
- LA-ICP-MS (1)
- Lacustrine sediments (1)
- Lagerstätte (1)
- Lake sediments (1)
- Lamina thickness (1)
- Landsat (1)
- Landsat-8 (1)
- Landslide (1)
- Larix larch (1)
- Last Interglacial (1)
- Last interglacial-glacial transition (1)
- Lawsonite (1)
- Lena Delta (1)
- LiDAR (1)
- Lithium (1)
- Lithospheric structure (1)
- Local knowledge (1)
- Local-to-regional scale (1)
- Lu/Hf dating of garnet (1)
- Machine learning (1)
- Magma migration and fragmentation (1)
- Magnetic properties (1)
- Magnetostratigraphy (1)
- Maiella Mountains (1)
- Mantle source (1)
- Martian regolith analogs (1)
- Mass extinction (1)
- Maule megathrust (1)
- Mediterranean (1)
- Mediterranean shrubland (1)
- Mekong Delta (1)
- Metamorphic evolution (1)
- Metamorphic sole (1)
- Metering (1)
- Methan (1)
- Methane (1)
- Methane leakage (1)
- Methanogenic archaea (1)
- Microcrack (1)
- Microseismic monitoring (1)
- Middle Holocene (1)
- Millennial-scale variability (1)
- Minimum mortality temperature (1)
- Model development (1)
- Model-data comparison (1)
- Mountain building (1)
- Multidimensionality (1)
- Muscovit (1)
- N isotopes (1)
- Namibia (1)
- Namibian passive margin (1)
- Natural Hazards (1)
- Naturgefahren (1)
- Neotectonics (1)
- Neotethys (1)
- North America (1)
- North Atlantic Oscillation (1)
- Northeast Iran (1)
- Northern Andes (1)
- Northern Asia (1)
- Numerical model (1)
- Numerical modelling (1)
- OAE-2 (1)
- OLCI (1)
- OSL (1)
- OSL and C-14 geochronology (1)
- Ocean Circulation (1)
- Operational use (1)
- Optical absorption spectroscopy (1)
- Optimal city size distribution (1)
- Orfento Formation (1)
- Organic matter degradation (1)
- Orogenic wedge (1)
- Ozone (1)
- Palaeoclimate reconstruction (1)
- Palaeoclimatology (1)
- Palaeodrainage (1)
- Palaeoecology (1)
- Paleoarchean (1)
- Paleoclimate dynamics (1)
- Paleocurrent (1)
- Paläoklima (1)
- Paläomagnetik (1)
- Paläosekularvariation (1)
- Panasqueira (1)
- Paris Agreement (1)
- Particle flow code (1)
- Passiv-seismische Interferometrie (1)
- Passive seismic interferometry (1)
- Peat (1)
- Per capita growth rate (1)
- Perchlorate (1)
- Permafrost degradation (1)
- Pharmaceuticals (1)
- Phosphates (1)
- Photogrammetrie (1)
- Photogrammetry (1)
- Phytolith dissolution (1)
- Phytolith morphotypes (1)
- Pirquitas (1)
- Pliocene (1)
- Pohang (Korea) (1)
- Porentypen in Karbonate (1)
- Precession (1)
- Predictive models (1)
- Preparedness (1)
- Principal component analysis (1)
- Proportional population growth rate (1)
- Provenance (1)
- Provenance studies (1)
- Psychologie (1)
- Pu239+240 (1)
- Public-private partnerships (1)
- Puna (1)
- Quantification of runoff components (1)
- Quaternary (1)
- REE (1)
- Radiogenic isotopes (1)
- Radioisotope disequilibria dating (1)
- Raman (1)
- Randelementmethode (1)
- Recurrence network (1)
- Red River (1)
- Reefs (1)
- Regional cooperation (1)
- Remote sensing (1)
- Review (1)
- Rheologie (1)
- Rheology (1)
- Risk perception (1)
- Risk reduction (1)
- River-groundwater-interaction (1)
- Rotifers (1)
- Rutile mineral chemistry (1)
- Röntgenabsorptionsspektroskopie (1)
- SIMS (1)
- Salt pan (1)
- Salt transport (1)
- Salztransport (1)
- Sandstone detrital modes (1)
- Schmelzstruktur (1)
- Schwarzes Meer (1)
- Sea of Marmara (1)
- Seasonal forecasting (1)
- Sediment connectivity (1)
- Sediment load (1)
- Sedimentary facies (1)
- Sedimentology (1)
- Seismic interferometry (1)
- Seismogenic sources (1)
- Seismotectonic segmentation (1)
- Sentinel 3 (1)
- Sentinel-1 (1)
- Sentinel-2 (1)
- Shale gas (1)
- Shallow borehole sensors (1)
- Shallow marine deposits (1)
- Shallow subsurface (1)
- Si cycling (1)
- Si extraction (1)
- Sigmoid model (1)
- Simulated Mars-like conditions (1)
- Sinkholes (1)
- Site proxy (1)
- Siwalik Group (1)
- Skarn (1)
- Sn (1)
- Social-ecological systems (1)
- Socio-economy (1)
- Sociocultural valuation (1)
- Soil organic carbon storage (1)
- South-America (1)
- Southern Oscillation (1)
- Species distribution modelling (1)
- Spektroskopie (1)
- Speleothems (1)
- Spin transition (1)
- Stability properties (1)
- Stable isotopes (1)
- Stage-discharge (1)
- Starkregen (1)
- Statistical methods (1)
- Stress measurement (1)
- Structural geology (1)
- Sturzfluten (1)
- Störungen (1)
- Störungszone (1)
- Subduction (1)
- Subduction Zone (1)
- Subfreezing temperatures (1)
- Submarine permafrost (1)
- Submariner Permafrost (1)
- Subsea permafrost (1)
- Surface roughness parameters (1)
- Surface-to-borehole spectral ratios (1)
- Sustainable natural hazard management (1)
- Syntectonic sedimentation (1)
- TRMM (1)
- Tajik Basin (1)
- Talik (1)
- Tangible losses (1)
- Tarim Basin (1)
- Tarim basin (1)
- Tectonic reconstruction (1)
- Tektonik (1)
- Temperature reconstruction (1)
- Temperature sensitivity (1)
- Temperaturerekonstruktion (1)
- Tephrostratigraphy (1)
- TerraSAR-X (1)
- TerraceM (1)
- The Geysers (1)
- Theoretical seismology (1)
- Thermal evolution (1)
- Tibetan plateau (1)
- Tien Shan Mountains (1)
- Tonminerale (1)
- Topography (1)
- Transfer functions (1)
- Transpression (1)
- Tree-ring software (1)
- Treibhausgase (1)
- Triangular dislocations (TDs) (1)
- Triassic (1)
- Triphylite (1)
- Tropical storms (1)
- Turmalin (1)
- U-Pb geochronology (1)
- UAVs (1)
- UNFCCC (1)
- Uncertainty (1)
- Urban water cycle (1)
- Varve chronologies (1)
- Varved sediments (1)
- Varves (1)
- Vegetation (1)
- Vegetation driver (1)
- Viscous blocking (1)
- Volcan de Colima (1)
- Volcano deformation modelling (1)
- Volcano monitoring (1)
- Volcano seismology (1)
- Vorlandbecken (1)
- Vulkanverformungsmodellierung (1)
- WRF-Chem (1)
- Wasserhaushaltsgesetz (1)
- Water stable isotope (1)
- Wave propagation (1)
- Waveform inversion (1)
- West Africa (1)
- Wiederkehr (1)
- Willingness to pay (1)
- Wofram-Zinn Lagerstätte (1)
- Wuchiapingian (1)
- X-ray absorption spectroscopy (1)
- Xining Basin (1)
- Yarlung Tsangpo (1)
- Young's modulus (1)
- Zinn (1)
- Zooplankton (1)
- Zusammengesetztes Dislokationsmodel (CDM) (1)
- acceptability (1)
- accumulation in soils (1)
- administrative units (1)
- agent-based models (1)
- agglomeration (1)
- agricultural productivity (1)
- agricultural soils (1)
- albedo (1)
- alkaline granites (1)
- alpine grassland degradation (1)
- alteration geochemistry (1)
- ambition (1)
- ammonia (1)
- analog models (1)
- anatexis (1)
- anoxia (1)
- apatite (1)
- apatite (U-Th-Sm)/He thermochronology (1)
- area-average snow monitoring (1)
- aridification (1)
- arktische Tundra (1)
- assessment (1)
- atmosphere (1)
- attenuation tomography (1)
- attitudes (1)
- attraction (1)
- automatic lineament extraction (1)
- bank infiltration (1)
- best practices (1)
- biomarker (1)
- biostratigraphy (1)
- bolivianischen Zinngürtels (1)
- boreal forests (1)
- boron (1)
- boundary element method (1)
- boundary element modeling (1)
- braided alluvial rivers (1)
- bromeliad (1)
- bryophytes (1)
- caldera collapse (1)
- carbon (1)
- carbon burial (1)
- carbon capture and storage (CCS) (1)
- carbon capture and utilization (CCU) (1)
- carbon flows (1)
- carbon sequestration (1)
- carbonate assimilation (1)
- carbonate pore types (1)
- catchment (1)
- channel migration (1)
- circulation (1)
- clay mineral (1)
- climate change co-operation (1)
- climate change cooperation (1)
- climate clubs (1)
- closed chamber method (1)
- club goods (1)
- clustering (1)
- co-benefits (1)
- coastal erosion (1)
- coastal geomorphology (1)
- collision (1)
- complex systems (1)
- compliance (1)
- compositional gap (1)
- compound flood (1)
- conditional commitments (1)
- connectivity (1)
- contamination processes (1)
- continental break-up (1)
- continuous modelling (1)
- controlo estrutural (1)
- correlation skill (1)
- corridors (1)
- cosmic-ray neutron sensing (1)
- cosmogenic nuclides (1)
- cost of electricity (COE) (1)
- cost-benefit analysis (1)
- costs of infrastructure (1)
- creep properties (1)
- crustal deformations (1)
- crustal fault (1)
- data (1)
- decomposition (1)
- deep carbon cycle (1)
- deep learning (1)
- deforestation (1)
- delta drift (1)
- dendroclimatology (1)
- denudation (1)
- dependence (1)
- depositional environment (1)
- depósitos de volfrâmio-estanho (1)
- derived flood frequency (1)
- detrital zircons (1)
- dietary patterns (1)
- digital image analysis (1)
- digitale Bildanalyse (1)
- disaster risk reduction (1)
- disinfection (1)
- dislocation creep (1)
- dissolution properties (1)
- drainage networks (1)
- dreieckige Dislokationen (TDs) (1)
- dynamical seasonal prediction (1)
- décollement (1)
- earthquake (1)
- earthquake modeling (1)
- earthquake nucleation (1)
- earthquake sequences (1)
- earthquake swarms (1)
- economics (1)
- ecosystem services (1)
- effective porosity (1)
- effektive Porosität (1)
- efficient method (1)
- electrical resistivity (1)
- electricity generation (1)
- electron lifetime (1)
- electron transport (1)
- enhanced geothermal systems (EGS) (1)
- ensemble modeling (1)
- epithermale Zinn-Silber-lagerstätte (1)
- equatorial ionosphere (1)
- ererbte Strukturen (1)
- erosion (1)
- erosion rates (1)
- evaporites (1)
- exhumation (1)
- exsolution textures (1)
- extinction (1)
- fault zone (1)
- faults (1)
- fenite (1)
- finite element modeling (1)
- flash floods (1)
- flood (1)
- flood genesis (1)
- flood hazard (1)
- flood loss model transfer (1)
- flood mechanisms (1)
- flood typology (1)
- floods (1)
- flow cytometry (1)
- flowering time (1)
- fluid flow (1)
- fluid regime (1)
- fluid-rock interaction (1)
- fluids (1)
- fluorescence (1)
- fluxo de fluidos (1)
- folds (1)
- food systems (1)
- fractional snow cover (1)
- fracture-controlled (1)
- functional trait correlations (1)
- gender (1)
- gene sequencing (1)
- geomagnetic excursions (1)
- geomagnetic indices (1)
- geomagnetische Exkursionen (1)
- geomagnetism (1)
- geomorphic markers (1)
- geomorphology (1)
- geoquímica da alteração hidrotermal (1)
- geschlossene Haubenmessmethode (1)
- glass structure (1)
- granite (1)
- gravity observations (1)
- gravity-driven slope deformation (1)
- green infrastructure (1)
- greenhouse gas emissions (1)
- greenhouse gases (1)
- groundwater (1)
- heavy rainfall (1)
- herd composition (1)
- high mountain environments (1)
- high-density reservoir network (1)
- hiss waves (1)
- historical floods (1)
- human activity (1)
- hummocky (1)
- hydroclimatology of floods (1)
- hydrothermal alteration (1)
- hyperspectral remote sensing (1)
- hyperspektral Fernerkundung (1)
- image enhancement (1)
- imaging spectroscopy (1)
- incision (1)
- individual-based model (1)
- inherited structures (1)
- inner magnetosphere (1)
- intangible impacts (1)
- interception (1)
- interferometric SAR (InSAR) (1)
- intermontane basins (1)
- intrabasinal faulting (1)
- intraplate deformation (1)
- invasive species (1)
- inverse Modellirung (1)
- inverse analysis (1)
- irrigation (1)
- isotopic methods (1)
- kettle hole (1)
- komplexe Systeme (1)
- landscape evolution (1)
- landscape evolution modeling (1)
- landslide (1)
- landslides (1)
- late Holocene (1)
- lava dome (1)
- lichens (1)
- lidar (1)
- life-history strategy (1)
- lithium (1)
- lithosphere dynamics (1)
- livestock (1)
- load stress (1)
- local adaptation (1)
- local food (1)
- lumped parameter model (1)
- magmatic processes (1)
- magmatic-hydrothermal systems (1)
- magmatisch-hydrothermale Systeme (1)
- magnetospheric convection (1)
- mantle evolution (1)
- mantle formation (1)
- marine terraces (1)
- mean transit time (1)
- megathrust earthquake (1)
- melt structure (1)
- melting (1)
- methanol (1)
- methods (1)
- microbialite (1)
- microfossil (1)
- microscopy (1)
- migration (1)
- miscibility gap (1)
- modelling systems (1)
- monsoon (1)
- monsoon onset (1)
- moraine landscape (1)
- morphometry (1)
- movement (1)
- multidrug-resistant Escherichia coli (1)
- multilevel probabilistic flood loss model (1)
- muscovite (1)
- nanocomposite material (1)
- nearshore zone (1)
- neotectonics (1)
- neutron simulations (1)
- noble gases (1)
- non-vascular epiphyte (1)
- nontidal ocean loading (1)
- normal faults (1)
- northwestern Europe (1)
- numerical modeling (1)
- numerical simulation (1)
- ocean color remote sensing (1)
- ocean-crust formation (1)
- optical remote sensing (1)
- ore deposit (1)
- orogeny (1)
- oxygen fractionation model (1)
- palaeoclimate (1)
- paleoceanography (1)
- paleosecular variations (1)
- paleosol (1)
- parallel evolution (1)
- part I (1)
- peatlands (1)
- peri-urban agriculture (1)
- permafrost ecosystems (1)
- phagotrophy (1)
- phenotypic plasticity (1)
- phylotypes (1)
- physical experiments (1)
- phytolith (1)
- phytoplankton (1)
- pitch angle diffusion coefficient (1)
- plate boundary fault (1)
- plate coupling force (1)
- plate reconstructions (1)
- plume-ridge interaction (1)
- point clouds (1)
- pollution (1)
- pre-monsoon (1)
- preferences (1)
- process simulation (1)
- prograding lobes (1)
- project database (1)
- proto-Paratethys (1)
- pseudomorphism (1)
- psychology (1)
- radiation belts (1)
- radiogenic isotopes (1)
- rain attenuation (1)
- rain effect (1)
- rainy-season (1)
- rangeland management (1)
- reactive oxygen species (1)
- recurrence (1)
- recurrence network (1)
- recurrence plot (1)
- refuge (1)
- regression (1)
- repulsion (1)
- reservoir volume (1)
- resilience (1)
- resources (1)
- rheological weakening (1)
- ring current electrons (1)
- ring faulting (1)
- risk mitigation (1)
- risk modeling (1)
- risk reduction (1)
- river erosion (1)
- river incision (1)
- river networks (1)
- river plume (1)
- river restoration (1)
- rock glacier (1)
- räumliche Wiederkehr (1)
- salt diffusion (1)
- satellite data (1)
- sea level change (1)
- sea turtle (1)
- seasonal precipitation (1)
- sediment routing (1)
- sediment-transport model (1)
- sedimentary loading and unloading cycles (1)
- seismic hazard (1)
- seismisches Hintergrundrauschen (1)
- seismotectonic segmentation (1)
- semi-arid (1)
- sentinel-1 (1)
- shrub encroachment (1)
- shrubline (1)
- siderite (1)
- silvopastoral systems (1)
- simulation model (1)
- simulações numéricas (1)
- sistemas magmático-hidrotermais (1)
- skarn (1)
- slow slip (1)
- small baseline subset (SBAS) (1)
- snow characteristics (1)
- snow grain size (1)
- soil erosion (1)
- solid Earth degassing (1)
- source inversion (1)
- spatial heterogeneity (1)
- spatial planning (1)
- spatial recurrence (1)
- spatial sampling (1)
- specific surface area (1)
- speleothems (1)
- spin transition (1)
- standards (1)
- strain localization (1)
- strain transient (1)
- strainmeter data (1)
- strength (1)
- structural mapping (1)
- strukturelle Kontrolle (1)
- subduction zone (1)
- subsea permafrost (1)
- summer (1)
- surface processes (1)
- surface runoff (1)
- surveys (1)
- suspended particulate matter (1)
- suspended sediment (1)
- synthetic flood events (1)
- systems (1)
- techno-economic model (1)
- tectonics (1)
- tectonism (1)
- temporal disaggregation (1)
- temporal evolution (1)
- the Ogooue Delta (1)
- tiefes Lernen (1)
- tillage (1)
- time series analysis (1)
- topographic analysis (1)
- torsion (1)
- tourmaline (1)
- trace elements (1)
- transform faults (1)
- transient response (1)
- tundra-taiga ecotone (1)
- tundra–taiga ecotone (1)
- tungsten-tin deposits (1)
- turbidity (1)
- turmalina (1)
- two phase aggregates (1)
- underground coal gasification (UCG) (1)
- urbanization (1)
- vadose zone (1)
- vascular epiphyte (1)
- vegetation (1)
- vegetation states (1)
- vegetation trajectories (1)
- vegetation-climate-fire-soil feedbacks (1)
- volcanic deformation (1)
- water (1)
- water height-area-volume curve (1)
- water storage capacity (1)
- wave number (1)
- weathering (1)
- welfare (1)
- wind speed (1)
- zoning (1)
- Überflutung (1)
Institute
- Institut für Geowissenschaften (276) (remove)
Sea surface temperature (SST) patterns can – as surface climate forcing – affect weather and climate at large distances. One example is El Niño-Southern Oscillation (ENSO) that causes climate anomalies around the globe via teleconnections. Although several studies identified and characterized these teleconnections, our understanding of climate processes remains incomplete, since interactions and feedbacks are typically exhibited at unique or multiple temporal and spatial scales. This study characterizes the interactions between the cells of a global SST data set at different temporal and spatial scales using climate networks. These networks are constructed using wavelet multi-scale correlation that investigate the correlation between the SST time series at a range of scales allowing instantaneously deeper insights into the correlation patterns compared to traditional methods like empirical orthogonal functions or classical correlation analysis. This allows us to identify and visualise regions of – at a certain timescale – similarly evolving SSTs and distinguish them from those with long-range teleconnections to other ocean regions. Our findings re-confirm accepted knowledge about known highly linked SST patterns like ENSO and the Pacific Decadal Oscillation, but also suggest new insights into the characteristics and origins of long-range teleconnections like the connection between ENSO and Indian Ocean Dipole.
Sea surface temperature (SST) patterns can – as surface climate forcing – affect weather and climate at large distances. One example is El Niño-Southern Oscillation (ENSO) that causes climate anomalies around the globe via teleconnections. Although several studies identified and characterized these teleconnections, our understanding of climate processes remains incomplete, since interactions and feedbacks are typically exhibited at unique or multiple temporal and spatial scales. This study characterizes the interactions between the cells of a global SST data set at different temporal and spatial scales using climate networks. These networks are constructed using wavelet multi-scale correlation that investigate the correlation between the SST time series at a range of scales allowing instantaneously deeper insights into the correlation patterns compared to traditional methods like empirical orthogonal functions or classical correlation analysis. This allows us to identify and visualise regions of – at a certain timescale – similarly evolving SSTs and distinguish them from those with long-range teleconnections to other ocean regions. Our findings re-confirm accepted knowledge about known highly linked SST patterns like ENSO and the Pacific Decadal Oscillation, but also suggest new insights into the characteristics and origins of long-range teleconnections like the connection between ENSO and Indian Ocean Dipole.
The Central Asian Pamir Mountains (Pamirs) are a high-altitude region sensitive to climatic change, with only few paleoclimatic records available. To examine the glacial-interglacial hydrological changes in the region, we analyzed the geochemical parameters of a 31-kyr record from Lake Karakul and performed a set of experiments with climate models to interpret the results. delta D values of terrestrial biomarkers showed insolation-driven trends reflecting major shifts of water vapor sources. For aquatic biomarkers, positive delta D shifts driven by changes in precipitation seasonality were observed at ca. 31-30, 28-26, and 17-14 kyr BP. Multiproxy paleoecological data and modelling results suggest that increased water availability, induced by decreased summer evaporation, triggered higher lake levels during those episodes, possibly synchronous to northern hemispheric rapid climate events. We conclude that seasonal changes in precipitation-evaporation balance significantly influenced the hydrological state of a large waterbody such as Lake Karakul, while annual precipitation amount and inflows remained fairly constant.
Sinkholes and depressions are typical landforms of karst regions. They pose a considerable natural hazard to infrastructure, agriculture, economy and human life in affected areas worldwide. The physio-chemical processes of sinkholes and depression formation are manifold, ranging from dissolution and material erosion in the subsurface to mechanical subsidence/failure of the overburden. This thesis addresses the mechanisms leading to the development of sinkholes and depressions by using complementary methods: remote sensing, distinct element modelling and near-surface geophysics.
In the first part, detailed information about the (hydro)-geological background, ground structures, morphologies and spatio-temporal development of sinkholes and depressions at a very active karst area at the Dead Sea are derived from satellite image analysis, photogrammetry and geologic field surveys. There, clusters of an increasing number of sinkholes have been developing since the 1980s within large-scale depressions and are distributed over different kinds of surface materials: clayey mud, sandy-gravel alluvium and lacustrine evaporites (salt). The morphology of sinkholes differs depending in which material they form: Sinkholes in sandy-gravel alluvium and salt are generally deeper and narrower than sinkholes in the interbedded evaporite and mud deposits. From repeated aerial surveys, collapse precursory features like small-scale subsidence, individual holes and cracks are identified in all materials. The analysis sheds light on the ongoing hazardous subsidence process, which is driven by the base-level fall of the Dead Sea and by the dynamic formation of subsurface water channels.
In the second part of this thesis, a novel, 2D distinct element geomechanical modelling approach with the software PFC2D-V5 to simulating individual and multiple cavity growth and sinkhole and large-scale depression development is presented. The approach involves a stepwise material removal technique in void spaces of arbitrarily shaped geometries and is benchmarked by analytical and boundary element method solutions for circular cavities. Simulated compression and tension tests are used to calibrate model parameters with bulk rock properties for the materials of the field site. The simulations show that cavity and sinkhole evolution is controlled by material strength of both overburden and cavity host material, the depth and relative speed of the cavity growth and the developed stress pattern in the subsurface. Major findings are: (1) A progressively deepening differential subrosion with variable growth speed yields a more fragmented stress pattern with stress interaction between the cavities. It favours multiple sinkhole collapses and nesting within large-scale depressions. (2) Low-strength materials do not support large cavities in the material removal zone, and subsidence is mainly characterised by gradual sagging into the material removal zone with synclinal bending. (3) High-strength materials support large cavity formation, leading to sinkhole formation by sudden collapse of the overburden. (4) Large-scale depression formation happens either by coalescence of collapsing holes, block-wise brittle failure, or gradual sagging and lateral widening.
The distinct element based approach is compared to results from remote sensing and geophysics at the field site. The numerical simulation outcomes are generally in good agreement with derived morphometrics, documented surface and subsurface structures as well as seismic velocities. Complementary findings on the subrosion process are provided from electric and seismic measurements in the area.
Based on the novel combination of methods in this thesis, a generic model of karst landform evolution with focus on sinkhole and depression formation is developed. A deepening subrosion system related to preferential flow paths evolves and creates void spaces and subsurface conduits. This subsequently leads to hazardous subsidence, and the formation of sinkholes within large-scale depressions. Finally, a monitoring system for shallow natural hazard phenomena consisting of geodetic and geophysical observations is proposed for similarly affected areas.
The 2-D distinct element method (DEM) code (PFC2D_V5) is used here to simulate the evolution of subsidence-related karst landforms, such as single and clustered sinkholes, and associated larger-scale depressions. Subsurface material in the DEM model is removed progressively to produce an array of cavities; this simulates a network of subsurface groundwater conduits growing by chemical/mechanical erosion. The growth of the cavity array is coupled mechanically to the gravitationally loaded surroundings, such that cavities can grow also in part by material failure at their margins, which in the limit can produce individual collapse sinkholes. Two end-member growth scenarios of the cavity array and their impact on surface subsidence were examined in the models: (1) cavity growth at the same depth level and growth rate; (2) cavity growth at progressively deepening levels with varying growth rates. These growth scenarios are characterised by differing stress patterns across the cavity array and its overburden, which are in turn an important factor for the formation of sinkholes and uvalalike depressions. For growth scenario (1), a stable compression arch is established around the entire cavity array, hindering sinkhole collapse into individual cavities and favouring block-wise, relatively even subsidence across the whole cavity array. In contrast, for growth scenario (2), the stress system is more heterogeneous, such that local stress concentrations exist around individual cavities, leading to stress interactions and local wall/overburden fractures. Consequently, sinkhole collapses occur in individual cavities, which results in uneven, differential subsidence within a larger-scale depression. Depending on material properties of the cavity-hosting material and the overburden, the larger-scale depression forms either by sinkhole coalescence or by widespread subsidence linked geometrically to the entire cavity array. The results from models with growth scenario (2) are in close agreement with surface morphological and subsurface geophysical observations from an evaporite karst area on the eastern shore of the Dead Sea.
Ground-penetrating radar is widely used to provide highly resolved images of subsurface sedimentary structures, with implications for processes active in the vadose zone. Frequently overlooked among these structures are tunnels excavated by fossorial animals (e.g., moles). We present two repeated ground-penetrating radar surveys performed a year apart in 2016 and 2017. Careful three-dimensional data processing reveals, in each data set, a pattern of elongated structures that are interpreted as a subsurface mole tunnel network. Our data demonstrate the ability of three-dimensional ground-penetrating radar imaging to non-invasively delineate the small animal tunnels (similar to 5 cm diameter) at a higher spatial and geolocation resolution than has previously been achieved. In turn, this makes repeated surveys and, therefore, long-term monitoring possible. Our results offer valuable insight into the understanding of the near-surface and showcase a potential new application for a geophysical method as well as a non-invasive method of ecological surveying.
Thawing of subsea permafrost can impact offshore infrastructure, affect coastal erosion, and release permafrost organic matter. Thawing is usually modeled as the result of heat transfer, although salt diffusion may play an important role in marine settings. To better quantify nearshore subsea permafrost thawing, we applied the CryoGRID2 heat diffusion model and coupled it to a salt diffusion model. We simulated coastline retreat and subsea permafrost evolution as it develops through successive stages of a thawing sequence at the Bykovsky Peninsula, Siberia. Sensitivity analyses for seawater salinity were performed to compare the results for the Bykovsky Peninsula with those of typical Arctic seawater. For the Bykovsky Peninsula, the modeled ice-bearing permafrost table (IBPT) for ice-rich sand and an erosion rate of 0.25m/year was 16.7 m below the seabed 350m offshore. The model outputs were compared to the IBPT depth estimated from coastline retreat and electrical resistivity surveys perpendicular to and crossing the shoreline of the Bykovsky Peninsula. The interpreted geoelectric data suggest that the IBPT dipped to 15-20m below the seabed at 350m offshore. Both results suggest that cold saline water forms beneath grounded ice and floating sea ice in shallow water, causing cryotic benthic temperatures. The freezing point depression produced by salt diffusion can delay or prevent ice formation in the sediment and enhance the IBPT degradation rate. Therefore, salt diffusion may facilitate the release of greenhouse gasses to the atmosphere and considerably affect the design of offshore and coastal infrastructure in subsea permafrost areas.
We present a model of the electrical resistivity structure of the lithosphere in the Central Andes between 20 degrees and 24 degrees S from 3-D inversion of 56 long-period magnetotelluric sites. Our model shows a complex resistivity structure with significant variability parallel and perpendicular to the trench direction. The continental forearc is characterized mainly by high electrical resistivity (>1,000m), suggesting overall low volumes of fluids. However, low resistivity zones (LRZs, <5m) were found in the continental forearc below areas where major trench-parallel faults systems intersect NW-SE transverse faults. Forearc LRZs indicate circulation and accumulation of fluids in highly permeable fault zones. The continental crust along the arc shows three distinctive resistivity domains, which coincide with segmentation in the distribution of volcanoes. The northern domain (20 degrees-20.5 degrees S) is characterized by resistivities >1,000m and the absence of active volcanism, suggesting the presence of a low-permeability block in the continental crust. The central domain (20.5 degrees-23 degrees S) exhibits a number of LRZs at varying depths, indicating different levels of a magmatic plumbing system. The southern domain (23 degrees-24 degrees S) is characterized by resistivities >1,000m, suggesting the absence of large magma reservoirs below the volcanic chain at crustal depths. Magma reservoirs located below the base of the crust or in the backarc may fed active volcanism in the southern domain. In the subcontinental mantle, the model exhibits LRZs in the forearc mantle wedge and above clusters of intermediate-depth seismicity, likely related to fluids produced by serpentinization of the mantle and eclogitization of the slab, respectively.
Ring current electrons (1–100 keV) have received significant attention in recent decades, but many questions regarding their major transport and loss mechanisms remain open. In this study, we use the four‐dimensional Versatile Electron Radiation Belt code to model the enhancement of phase space density that occurred during the 17 March 2013 storm. Our model includes global convection, radial diffusion, and scattering into the Earth's atmosphere driven by whistler‐mode hiss and chorus waves. We study the sensitivity of the model to the boundary conditions, global electric field, the electric field associated with subauroral polarization streams, electron loss rates, and radial diffusion coefficients. The results of the code are almost insensitive to the model parameters above 4.5 RERE, which indicates that the general dynamics of the electrons between 4.5 RE and the geostationary orbit can be explained by global convection. We found that the major discrepancies between the model and data can stem from the inaccurate electric field model and uncertainties in lifetimes. We show that additional mechanisms that are responsible for radial transport are required to explain the dynamics of ≥40‐keV electrons, and the inclusion of the radial diffusion rates that are typically assumed in radiation belt studies leads to a better agreement with the data. The overall effect of subauroral polarization streams on the electron phase space density profiles seems to be smaller than the uncertainties in other input parameters. This study is an initial step toward understanding the dynamics of these particles inside the geostationary orbit.
The novel space-borne Global Navigation Satellite System Reflectometry (GNSS-R) technique has recently shown promise in monitoring the ocean state and surface wind speed with high spatial coverage and unprecedented sampling rate. The L-band signals of GNSS are structurally able to provide a higher quality of observations from areas covered by dense clouds and under intense precipitation, compared to those signals at higher frequencies from conventional ocean scatterometers. As a result, studying the inner core of cyclones and improvement of severe weather forecasting and cyclone tracking have turned into the main objectives of GNSS-R satellite missions such as Cyclone Global Navigation Satellite System (CYGNSS). Nevertheless, the rain attenuation impact on GNSS-R wind speed products is not yet well documented. Evaluating the rain attenuation effects on this technique is significant since a small change in the GNSS-R can potentially cause a considerable bias in the resultant wind products at intense wind speeds. Based on both empirical evidence and theory, wind speed is inversely proportional to derived bistatic radar cross section with a natural logarithmic relation, which introduces high condition numbers (similar to ill-posed conditions) at the inversions to high wind speeds. This paper presents an evaluation of the rain signal attenuation impact on the bistatic radar cross section and the derived wind speed. This study is conducted simulating GNSS-R delay-Doppler maps at different rain rates and reflection geometries, considering that an empirical data analysis at extreme wind intensities and rain rates is impossible due to the insufficient number of observations from these severe conditions. Finally, the study demonstrates that at a wind speed of 30 m/s and incidence angle of 30 degrees, rain at rates of 10, 15, and 20 mm/h might cause overestimation as large as approximate to 0.65 m/s (2%), 1.00 m/s (3%), and 1.3 m/s (4%), respectively, which are still smaller than the CYGNSS required uncertainty threshold. The simulations are conducted in a pessimistic condition (severe continuous rainfall below the freezing height and over the entire glistening zone) and the bias is expected to be smaller in size in real environments.
OpenForecast
(2019)
The development and deployment of new operational runoff forecasting systems are a strong focus of the scientific community due to the crucial importance of reliable and timely runoff predictions for early warnings of floods and flashfloods for local businesses and communities. OpenForecast, the first operational runoff forecasting system in Russia, open for public use, is presented in this study. We developed OpenForecast based only on open-source software and data-GR4J hydrological model, ERA-Interim meteorological reanalysis, and ICON deterministic short-range meteorological forecasts. Daily forecasts were generated for two basins in the European part of Russia. Simulation results showed a limited efficiency in reproducing the spring flood of 2019. Although the simulations managed to capture the timing of flood peaks, they failed in estimating flood volume. However, further implementation of the parsimonious data assimilation technique significantly alleviates simulation errors. The revealed limitations of the proposed operational runoff forecasting system provided a foundation to outline its further development and improvement.
OpenForecast
(2019)
The development and deployment of new operational runoff forecasting systems are a strong focus of the scientific community due to the crucial importance of reliable and timely runoff predictions for early warnings of floods and flashfloods for local businesses and communities. OpenForecast, the first operational runoff forecasting system in Russia, open for public use, is presented in this study. We developed OpenForecast based only on open-source software and data-GR4J hydrological model, ERA-Interim meteorological reanalysis, and ICON deterministic short-range meteorological forecasts. Daily forecasts were generated for two basins in the European part of Russia. Simulation results showed a limited efficiency in reproducing the spring flood of 2019. Although the simulations managed to capture the timing of flood peaks, they failed in estimating flood volume. However, further implementation of the parsimonious data assimilation technique significantly alleviates simulation errors. The revealed limitations of the proposed operational runoff forecasting system provided a foundation to outline its further development and improvement.
The removal, redistribution, and transient storage of sediments in tectonically active mountain belts is thought to exert a first-order control on shallow crustal stresses, fault activity, and hence on the spatiotemporal pattern of regional deformation processes. Accordingly, sediment loading and unloading cycles in intermontane sedimentary basins may inhibit or promote intrabasinal faulting, respectively, but unambiguous evidence for this potential link has been elusive so far. Here we combine 2D numerical experiments that simulate contractional deformation in a broken-foreland setting (i.e., a foreland where shortening is diachronously absorbed by spatially disparate, reverse faults uplifting basement blocks) with field data from intermontane basins in the NW Argentine Andes. Our modeling results suggest that thicker sedimentary fills (>0.7-1.0 km) may suppress basinal faulting processes, while thinner fills (<0.7 km) tend to delay faulting. Conversely, the removal of sedimentary loads via fluvial incision and basin excavation promotes renewed intrabasinal faulting. These results help to better understand the tectono-sedimentary history of intermontane basins that straddle the eastern border of the Andean Plateau in northwestern Argentina. For example, the Santa Maria and the Humahuaca basins record intrabasinal deformation during or after sediment unloading, while the Quebrada del Toro Basin reflects the suppression of intrabasinal faulting due to loading by coarse conglomerates. We conclude that sedimentary loading and unloading cycles may exert a fundamental control on spatiotemporal deformation patterns in intermontane basins of tectonically active broken forelands. (C) 2018 Elsevier B.V. All rights reserved.
In Germany, the irrigation sector accounts for only 1% of water use. In recent years, however, this sector has attracted more attention due to the occurrence of severe drought periods. Irrigation scheduling systems could support adaptation strategies but little is known about current providers, performance and users. In this study we aimed to depict the current situation of the existence and functioning of irrigation scheduling systems available in Germany. Six methods were identified and assessed based on direct interviews with end-users and a comparative analysis. The results showed a positive feedback from the users. However, the recommendations were rarely implemented, while only the seasonal irrigation requirement was considered to support actual water abstraction. These results were corroborated by the comparative analysis. Five of the six irrigation scheduling systems estimated the seasonal irrigation amount consistently, while wider differences were found by looking at the irrigation season and at the number of irrigations. Overall, it is found that irrigation support systems are valuable tools for supporting adaptation strategies to fast changes in agro-environmental conditions. However, specific assessments based on real measurements should be considered in order to improve the performance of the systems and provide more consistent support to end-users. (c) 2019 John Wiley & Sons, Ltd.
In this paper, we examine the influence of the 27 October 2012, M-w 7.8 earthquake on landslide occurrence in the southern half of Haida Gwaii (formerly Queen Charlotte Islands), British Columbia, Canada. Our 1350 km(2) study area is undisturbed, primarily forested terrain that has not experienced road building or timber harvesting. Our inventory of landslide polygons is based on optical airborne and spaceborne images acquired between 2007 and 2018, from which we extracted and mapped 446 individual landslides (an average of 33 landslides per 100 km(2)). The landslide rate in years without major earthquakes averages 19.4 per year, or 1.4/100 km(2)/year, and the annual average area covered by non-seismically triggered landslides is 35 ha/year. The number of landslides identified in imagery closely following the 2012 earthquake, and probably triggered by it, is 244 or an average of about 18 landslides per 100 km(2). These landslides cover a total area of 461 ha. In the following years-2013-2016 and 2016-2018-the number of landslides fell, respectively, to 26 and 13.5 landslides per year. In non-earthquake years, most landslides happen on south-facing slopes, facing the prevailing winds. In contrast, during or immediately after the earthquake, up to 32% of the landslides occurred on north and northwest-facing slopes. Although we could not find imagery from the day after the earthquake, overview reconnaissance flights 10 and 16 days later showed that most of the landslides were recent, suggesting they were co-seismic.
The Seismic Hazard Inferred from Tectonics based on the Global Strain Rate Map (SHIFT_GSRM) earthquake forecast was designed to provide high-resolution estimates of global shallow seismicity to be used in seismic hazard assessment. This model combines geodetic strain rates with global earthquake parameters to characterize long-term rates of seismic moment and earthquake activity. Although SHIFT_GSRM properly computes seismicity rates in seismically active continental regions, it underestimates earthquake rates in subduction zones by an average factor of approximately 3. We present a complementary method to SHIFT_GSRM to more accurately forecast earthquake rates in 37 subduction segments, based on the conservation of moment principle and the use of regional interface seismicity parameters, such as subduction dip angles, corner magnitudes, and coupled seismogenic thicknesses. In seven progressive steps, we find that SHIFT_GSRM earthquake-rate underpredictions are mainly due to the utilization of a global probability function of seismic moment release that poorly captures the great variability among subduction megathrust interfaces. Retrospective test results show that the forecast is consistent with the observations during the 1 January 1977 to 31 December 2014 period. Moreover, successful pseudoprospective evaluations for the 1 January 2015 to 31 December 2018 period demonstrate the power of the regionalized earthquake model to properly estimate subduction-zone seismicity.
Hyperspectral remote sensing of the spatial and temporal heterogeneity of low Arctic vegetation
(2019)
Arctic tundra ecosystems are experiencing warming twice the global average and Arctic vegetation is responding in complex and heterogeneous ways. Shifting productivity, growth, species composition, and phenology at local and regional scales have implications for ecosystem functioning as well as the global carbon and energy balance. Optical remote sensing is an effective tool for monitoring ecosystem functioning in this remote biome. However, limited field-based spectral characterization of the spatial and temporal heterogeneity limits the accuracy of quantitative optical remote sensing at landscape scales. To address this research gap and support current and future satellite missions, three central research questions were posed:
• Does canopy-level spectral variability differ between dominant low Arctic vegetation communities and does this variability change between major phenological phases?
• How does canopy-level vegetation colour images recorded with high and low spectral resolution devices relate to phenological changes in leaf-level photosynthetic pigment concentrations?
• How does spatial aggregation of high spectral resolution data from the ground to satellite scale influence low Arctic tundra vegetation signatures and thereby what is the potential of upcoming hyperspectral spaceborne systems for low Arctic vegetation characterization?
To answer these questions a unique and detailed database was assembled. Field-based canopy-level spectral reflectance measurements, nadir digital photographs, and photosynthetic pigment concentrations of dominant low Arctic vegetation communities were acquired at three major phenological phases representing early, peak and late season. Data were collected in 2015 and 2016 in the Toolik Lake Research Natural Area located in north central Alaska on the North Slope of the Brooks Range. In addition to field data an aerial AISA hyperspectral image was acquired in the late season of 2016. Simulations of broadband Sentinel-2 and hyperspectral Environmental and Mapping Analysis Program (EnMAP) satellite reflectance spectra from ground-based reflectance spectra as well as simulations of EnMAP imagery from aerial hyperspectral imagery were also obtained.
Results showed that canopy-level spectral variability within and between vegetation communities differed by phenological phase. The late season was identified as the most discriminative for identifying many dominant vegetation communities using both ground-based and simulated hyperspectral reflectance spectra. This was due to an overall reduction in spectral variability and comparable or greater differences in spectral reflectance between vegetation communities in the visible near infrared spectrum.
Red, green, and blue (RGB) indices extracted from nadir digital photographs and pigment-driven vegetation indices extracted from ground-based spectral measurements showed strong significant relationships. RGB indices also showed moderate relationships with chlorophyll and carotenoid pigment concentrations. The observed relationships with the broadband RGB channels of the digital camera indicate that vegetation colour strongly influences the response of pigment-driven spectral indices and digital cameras can track the seasonal development and degradation of photosynthetic pigments.
Spatial aggregation of hyperspectral data from the ground to airborne, to simulated satel-lite scale was influenced by non-photosynthetic components as demonstrated by the distinct shift of the red edge to shorter wavelengths. Correspondence between spectral reflectance at the three scales was highest in the red spectrum and lowest in the near infra-red. By artificially mixing litter spectra at different proportions to ground-based spectra, correspondence with aerial and satellite spectra increased. Greater proportions of litter were required to achieve correspondence at the satellite scale.
Overall this thesis found that integrating multiple temporal, spectral, and spatial data is necessary to monitor the complexity and heterogeneity of Arctic tundra ecosystems. The identification of spectrally similar vegetation communities can be optimized using non-peak season hyperspectral data leading to more detailed identification of vegetation communities. The results also highlight the power of vegetation colour to link ground-based and satellite data. Finally, a detailed characterization non-photosynthetic ecosystem components is crucial for accurate interpretation of vegetation signals at landscape scales.
Growing attention to phytoplankton mixotrophy as a trophic strategy has led to significant revisions of traditional pelagic food web models and ecosystem functioning. Although some empirical estimates of mixotrophy do exist, a much broader set of in situ measurements are required to (i) identify which organisms are acting as mixotrophs in real time and to (ii) assess the contribution of their heterotrophy to biogeochemical cycling. Estimates are needed through time and across space to evaluate which environmental conditions or habitats favour mixotrophy: conditions still largely unknown. We review methodologies currently available to plankton ecologists to undertake estimates of plankton mixotrophy, in particular nanophytoplankton phago-mixotrophy. Methods are based largely on fluorescent or isotopic tracers, but also take advantage of genomics to identify phylotypes and function. We also suggest novel methods on the cusp of use for phago-mixotrophy assessment, including single-cell measurements improving our capacity to estimate mixotrophic activity and rates in wild plankton communities down to the single-cell level. Future methods will benefit from advances in nanotechnology, micromanipulation and microscopy combined with stable isotope and genomic methodologies. Improved estimates of mixotrophy will enable more reliable models to predict changes in food web structure and biogeochemical flows in a rapidly changing world.
The sedimentary record of the Dead Sea provides an exceptional high-resolution archive of past climate changes in the drought-sensitive eastern Mediterranean-Levant, a key region for the development of humankind at the boundary of global climate belts. Moreover, it is the only deep hypersaline lake known to have deposited long sequences of finely laminated, annually deposited sediments (i.e. varves) of varied compositions, including aragonite, gypsum, halite and clastic sediments. Vast efforts have been made over the years to decipher the environmental information stored in these evaporitic-clastic sequences spanning from the Pleistocene Lake Amora to the Holocene Dead Sea. A general characterisation of sediment facies has been derived from exposed sediment sections, as well as from shallow- and deep-water sediment cores. During high lake stands and episodes of positive water budget, mostly during glacial times, alternating aragonite and detritus laminae (‘aad’ facies) were accumulated, whereas during low lake stands and droughts, prevailing during interglacials, laminated detritus (‘ld’ facies) and laminated halite (‘lh’ facies) dominate the sequence. In this paper, we (i) review the three types of laminated sediments of the Dead Sea sedimentary record (‘aad’, ‘ld’ and ‘lh’ facies), (ii) discuss their modes of formation, deposition and accumulation, and their interpretation as varves, and (iii) illustrate how Dead Sea varves are utilized for palaeoclimate reconstructions and for establishing floating chronologies.
Preparatory mechanisms accompanying or leading to nucleation of larger earthquakes have been observed at both laboratory and field scales, but conditions favoring the occurrence of observable preparatory processes are still largely unknown. In particular, it remains a matter of debate why some earthquakes occur spontaneously without noticeable precursors as opposed to events that are preceded by an extended failure process. In this study, we have generated new high-resolution seismicity catalogs framing the occurrence of 20 M-L > 2.5 earthquakes at The Geysers geothermal field in California. To this end, a seismicity catalog of the 11 days framing each large event was created. We selected 20 sequences sampling different hypocentral depths and hydraulic conditions within the field. Seismic activity and magnitude frequency distributions displayed by the different earthquake sequences are correlated with their location within the reservoir. Sequences located in the northwestern part of the reservoir show overall increased seismic activity and low b values, while the southeastern part is dominated by decreased seismic activity and higher b values. Periods of high injection coincide with high b values and vice versa. These observations potentially reflect varying differential and mean stresses and damage of the reservoir rocks across the field. About 50% of analyzed sequences exhibit no change in seismicity rate in response to the large main event. However, we find complex waveforms at the onset of the main earthquake, suggesting that small ruptures spontaneously grow into or trigger larger events.
We derive a set of regional ground-motion prediction equations (GMPEs) in the Fourier amplitude spectra (FAS-GMPE) and in the spectral acceleration (SA-GMPE) domains for the purpose of interpreting the between-event residuals in terms of source parameter variability. We analyze a dataset of about 65,000 recordings generated by 1400 earthquakes (moment magnitude 2: 5 <= M-w <= 6: 5, hypocentral distance R-hypo <= 150 km) that occurred in central Italy between January 2008 and October 2017. In a companion article (Bindi, Spallarossa, et al., 2018), the nonparametric acceleration source spectra were interpreted in terms of omega-square models modified to account for deviations from a high-frequency flat plateau through a parameter named k(source). Here, the GMPEs are derived considering the moment (M-w), the local (M-L), and the energy (M-e) magnitude scales, and the between-event residuals are computed as random effects. We show that the between-event residuals for the FAS-GMPE implementing M-w are correlated with stress drop, with correlation coefficients increasing with increasing frequency up to about 10 Hz. Contrariwise, the correlation is weak for the FAS-GMPEs implementing M-L and M-e, in particular between 2 and 5 Hz, where most of the corner frequencies lie. At higher frequencies, all models show a strong correlation with k(source). The correlation with the source parameters reflects in a different behavior of the standard deviation tau of the between-event residuals with frequency. Although tau is smaller for the FAS-GMPE using M-w below 1.5 Hz, at higher frequencies, the model implementing either M-L or M-e shows smaller values, with a reduction of about 30% at 3 Hz (i.e., from 0.3 for M-w to 0.1 for M-L). We conclude that considering magnitude scales informative for the stress-drop variability allows to reduce the between-event variability with a significant impact on the hazard assessment, in particular for studies in which the ergodic assumption on site is removed.
Rapidly changing climate in the Northern Hemisphere and associated socio-economic impacts require reliable understanding of lake systems as important freshwater resources and sensitive sentinels of environmental change. To better understand time-series data in lake sediment cores, it is necessary to gain information on within-lake spatial variabilities of environmental indicator data. Therefore, we retrieved a set of 38 samples from the sediment surface along spatial habitat gradients in the boreal, deep, and yet pristine Lake Bolshoe Toko in southern Yakutia, Russia. Our methods comprise laboratory analyses of the sediments for multiple proxy parameters, including diatom and chironomid taxonomy, oxygen isotopes from diatom silica, grain-size distributions, elemental compositions (XRF), organic carbon content, and mineralogy (XRD). We analysed the lake water for cations, anions, and isotopes. Our results show that the diatom assemblages are strongly influenced by water depth and dominated by planktonic species, i.e. Pliocaenicus bolshetokoensis. Species richness and diversity are higher in the northern part of the lake basin, associated with the availability of benthic, i.e. periphytic, niches in shallower waters. delta O-18(diatom) values are higher in the deeper south-western part of the lake, probably related to water temperature differences. The highest amount of the chironomid taxa underrepresented in the training set used for palaeoclimate inference was found close to the Utuk River and at southern littoral and profundal sites. Abiotic sediment components are not symmetrically distributed in the lake basin, but vary along restricted areas of differential environmental forcing. Grain size and organic matter are mainly controlled by both river input and water depth. Mineral (XRD) data distributions are influenced by the methamorphic lithology of the Stanovoy mountain range, while elements (XRF) are intermingled due to catchment and diagenetic differences. We conclude that the lake represents a valuable archive for multiproxy environmental reconstruction based on diatoms (including oxygen isotopes), chironomids, and sediment-geochemical parameters. Our analyses suggest multiple coring locations preferably at intermediate depth in the northern basin and the deep part in the central basin, to account for representative bioindicator distributions and higher temporal resolution, respectively.
The northward indentation of the Pamir salient into the Tarim basin at the western syntaxis of the India-Asia collision zone is the focus of controversial models linking lithospheric to surface and atmospheric processes. Here we report on tectonic events recorded in the most complete and best-dated sedimentary sequences from the western Tarim basin flanking the eastern Pamir (the Aertashi section), based on sedimentologic, provenance, and magnetostratigraphic analyses. Increased tectonic subsidence and a shift from marine to continental fluvio-deltaic deposition at 41Ma indicate that far-field deformation from the south started to affect the Tarim region. A sediment accumulation hiatus from 24.3 to 21.6Ma followed by deposition of proximal conglomerates is linked to fault propagation into the Tarim basin. From 21.6 to 15.0Ma, increasing accumulation rates of fining upward clastics is interpreted as the expression of a major dextral transtensional system linking the Kunlun to the Tian Shan ahead of the northward Pamir indentation. At 15.0Ma, the appearance of North Pamir-sourced conglomerates followed at 11Ma by Central Pamir-sourced volcanics coincides with a shift to E-W compression, clockwise vertical-axis rotations and the onset of growth strata associated with the activation of the local east vergent Qimugen thrust wedge. Together, this enables us to interpret that Pamir indentation into Tarim had started by 24.3Ma, reached the study location by 15.0Ma and had passed it by 11Ma, providing kinematic constraints on proposed tectonic models involving intracontinental subduction and delamination.
Secondary mica minerals collected from the Santa Helena (W- (Cu) mineralization) and Venise (W-Mo mineralization) endogenic breccia structures were Ar-40/Ar-39 dated. The muscovite Ar-40/Ar-39 data yielded 286.8 +/- 1.2 (+/- 1 sigma) Ma (samples 6Ha and 11Ha) which reflect the age of secondary muscovite formation probably from magmatic biotite or feldspar alteration. Sericite Ar-40/Ar-39 data yielded 280.9 +/- 1.2 (+/- 1 sigma) Ma to 279.0 +/- 1.1 (+/- 1 sigma) Ma (samples 6Hb and 11Hb) reflecting the age of greisen alteration (T similar to 300 degrees C) where the W- disseminated mineralization occurs. The muscovite 40Ar/39Ar data of 277.3 +/- 1.3 (+/- 1 sigma) Ma and 281.3 +/- 1.2 (+/- 1 sigma) Ma (samples 5 and 6) also reflect the age of muscovite (selvage) crystallized adjacent to molybdenite veins within the Venise breccia. Geochronological data obtained confirmed that the W mineralization at Santa Helena breccia is older than Mo-mineralization at Venise breccia. Also, the timing of hydrothermal circulation and the cooling history for the W-stage deposition was no longer than 7 Ma and 4 Ma for Mo-deposition.
Climatic observables are often correlated across long spatial distances, and extreme events, such as heatwaves or floods, are typically assumed to be related to such teleconnections(1,2). Revealing atmospheric teleconnection patterns and understanding their underlying mechanisms is of great importance for weather forecasting in general and extreme-event prediction in particular(3,4), especially considering that the characteristics of extreme events have been suggested to change under ongoing anthropogenic climate change(5-8). Here we reveal the global coupling pattern of extreme-rainfall events by applying complex-network methodology to high-resolution satellite data and introducing a technique that corrects for multiple-comparison bias in functional networks. We find that the distance distribution of significant connections (P < 0.005) around the globe decays according to a power law up to distances of about 2,500 kilometres. For longer distances, the probability of significant connections is much higher than expected from the scaling of the power law. We attribute the shorter, power-law-distributed connections to regional weather systems. The longer, super-power-law-distributed connections form a global rainfall teleconnection pattern that is probably controlled by upper-level Rossby waves. We show that extreme-rainfall events in the monsoon systems of south-central Asia, east Asia and Africa are significantly synchronized. Moreover, we uncover concise links between south-central Asia and the European and North American extratropics, as well as the Southern Hemisphere extratropics. Analysis of the atmospheric conditions that lead to these teleconnections confirms Rossby waves as the physical mechanism underlying these global teleconnection patterns and emphasizes their crucial role in dynamical tropical-extratropical couplings. Our results provide insights into the function of Rossby waves in creating stable, global-scale dependencies of extreme-rainfall events, and into the potential predictability of associated natural hazards.
Measuring (dis)similarity between ecosystem states is a key theme in ecology. Much of community and ecosystem ecology is devoted to searching for patterns in ecosystem similarity from an external observer's viewpoint, using variables such as species abundances, measures of diversity and complexity. However, from the point of view of organisms in the ecosystem, proportional population growth rates are the only relevant aspect of ecosystem state, because natural selection acts on groups of organisms with different proportional population growth rates. We therefore argue that two ecosystem states are equivalent if and only if, for each species they contain, the proportional population growth rate does not differ between the states. Based on this result, we develop species-level and aggregated summary measures of ecosystem state and discuss their ecological meaning. We illustrate our approach using a long-term dataset on the plankton community from the Central European Lake Constance. We show that the first three principal components of proportional population growth rates describe most of the variation in ecosystem state in Lake Constance. We strongly recommend using proportional population growth rates and the derived equivalence classes for comparative ecosystem studies. This opens up new perspectives on important existing topics such as alternative stable ecosystem states, community assembly, and the processes generating regularities in ecosystems.
Adjustment of median ground motion prediction equations (GMPEs) from one region to another region is one of the major challenges within the current practice of seismic hazard analysis. In our approach of generating response spectra, we derive two separate empirical models for a) Fourier amplitude spectrum (FAS) and b) duration of ground motion. To calculate response spectra, the two models are combined within the random vibration theory (RVT) framework. The models are calibrated on recordings obtained from shallow crustal earthquakes in active tectonic regions. We use a subset of NGA-West2 database with M3.2-7.9 earthquakes at distances 0-300 km. The NGA-West2 database expanded over a wide magnitude and distance range facilitates a better constraint over derived models. A frequency-dependent duration model is derived to obtain adjustable response spectral ordinates. Excellent comparison of our approach with other NGA-West2 models implies that it can also be used as a stand-alone model.
We study optimal and equilibrium sizes of cities in a city system model with pollution. Pollution is a function of population size. If pollution is local or per-capita pollution increases with population, equilibrium cities are too large under symmetry; with asymmetric cities, the largest cities are too large and the smallest too small. When pollution is global and per-capita pollution declines with city size, cities may be too small under symmetry; with asymmetric cities, the largest cities are too small and the smallest too large if the marginal damage of pollution is large enough. We calibrate the model to US cities and find that the largest cities may be undersized by 3-4%.
Garnet of eclogite (formerly termed garnet clinopyroxenite) hosted in lenses of orogenic garnet peridotite from the Granulitgebirge, NW Bohemian Massif, contains unique inclusions of granitic melt, now either glassy or crystallized. Analysed glasses and re‐homogenized inclusions are hydrous, peraluminous, and enriched in highly incompatible elements characteristic of the continental crust such as Cs, Li, B, Pb, Rb, Th, and U. The original melt thus represents a pristine, chemically evolved metasomatic agent, which infiltrated the mantle via deep continental subduction during the Variscan orogeny. The bulk chemical composition of the studied eclogites is similar to that of Fe‐rich basalt and the enrichment in LILE and U suggest a subduction‐related component. All these geochemical features confirm metasomatism. In comparison with many other garnet+clinopyroxene‐bearing lenses in peridotites of the Bohemian Massif, the studied samples from Rubinberg and Klatschmühle are more akin to eclogite than pyroxenites, as reflected in high jadeite content in clinopyroxene, relatively low Mg, Cr, and Ni but relatively high Ti. However, trace elements of both bulk rock and individual mineral phases show also important differences making these samples rather unique. Metasomatism involving a melt requiring a trace element pattern very similar to the composition reported here has been suggested for the source region of rocks of the so‐called durbachite suite, that is, ultrapotassic melanosyenites, which are found throughout the high‐grade Variscan basement. Moreover, the Th, U, Pb, Nb, Ta, and Ti patterns of these newly studied melt inclusions (MI) strongly resemble those observed for peridotite and its enclosed pyroxenite from the T‐7 borehole (Staré, České Středhoři Mountains) in N Bohemia. This suggests that a similar kind of crustal‐derived melt also occurred here. This study of granitic MI in eclogites from peridotites has provided the first direct characterization of a preserved metasomatic melt, possibly responsible for the metasomatism of several parts of the mantle in the Variscides.
Garnet of eclogite (formerly termed garnet clinopyroxenite) hosted in lenses of orogenic garnet peridotite from the Granulitgebirge, NW Bohemian Massif, contains unique inclusions of granitic melt, now either glassy or crystallized. Analysed glasses and re‐homogenized inclusions are hydrous, peraluminous, and enriched in highly incompatible elements characteristic of the continental crust such as Cs, Li, B, Pb, Rb, Th, and U. The original melt thus represents a pristine, chemically evolved metasomatic agent, which infiltrated the mantle via deep continental subduction during the Variscan orogeny. The bulk chemical composition of the studied eclogites is similar to that of Fe‐rich basalt and the enrichment in LILE and U suggest a subduction‐related component. All these geochemical features confirm metasomatism. In comparison with many other garnet+clinopyroxene‐bearing lenses in peridotites of the Bohemian Massif, the studied samples from Rubinberg and Klatschmühle are more akin to eclogite than pyroxenites, as reflected in high jadeite content in clinopyroxene, relatively low Mg, Cr, and Ni but relatively high Ti. However, trace elements of both bulk rock and individual mineral phases show also important differences making these samples rather unique. Metasomatism involving a melt requiring a trace element pattern very similar to the composition reported here has been suggested for the source region of rocks of the so‐called durbachite suite, that is, ultrapotassic melanosyenites, which are found throughout the high‐grade Variscan basement. Moreover, the Th, U, Pb, Nb, Ta, and Ti patterns of these newly studied melt inclusions (MI) strongly resemble those observed for peridotite and its enclosed pyroxenite from the T‐7 borehole (Staré, České Středhoři Mountains) in N Bohemia. This suggests that a similar kind of crustal‐derived melt also occurred here. This study of granitic MI in eclogites from peridotites has provided the first direct characterization of a preserved metasomatic melt, possibly responsible for the metasomatism of several parts of the mantle in the Variscides.
Forest structure is a crucial component in the assessment of whether a forest is likely to act as a carbon sink under changing climate. Detailed 3D structural information about the tundra–taiga ecotone of Siberia is mostly missing and still underrepresented in current research due to the remoteness and restricted accessibility. Field based, high-resolution remote sensing can provide important knowledge for the understanding of vegetation properties and dynamics. In this study, we test the applicability of consumer-grade Unmanned Aerial Vehicles (UAVs) for rapid calculation of stand metrics in treeline forests. We reconstructed high-resolution photogrammetric point clouds and derived canopy height models for 10 study sites from NE Chukotka and SW Yakutia. Subsequently, we detected individual tree tops using a variable-window size local maximum filter and applied a marker-controlled watershed segmentation for the delineation of tree crowns. With this, we successfully detected 67.1% of the validation individuals. Simple linear regressions of observed and detected metrics show a better correlation (R2) and lower relative root mean square percentage error (RMSE%) for tree heights (mean R2 = 0.77, mean RMSE% = 18.46%) than for crown diameters (mean R2 = 0.46, mean RMSE% = 24.9%). The comparison between detected and observed tree height distributions revealed that our tree detection method was unable to representatively identify trees <2 m. Our results show that plot sizes for vegetation surveys in the tundra–taiga ecotone should be adapted to the forest structure and have a radius of >15–20 m to capture homogeneous and representative forest stands. Additionally, we identify sources of omission and commission errors and give recommendations for their mitigation. In summary, the efficiency of the used method depends on the complexity of the forest’s stand structure.
Forest structure is a crucial component in the assessment of whether a forest is likely to act as a carbon sink under changing climate. Detailed 3D structural information about the tundra–taiga ecotone of Siberia is mostly missing and still underrepresented in current research due to the remoteness and restricted accessibility. Field based, high-resolution remote sensing can provide important knowledge for the understanding of vegetation properties and dynamics. In this study, we test the applicability of consumer-grade Unmanned Aerial Vehicles (UAVs) for rapid calculation of stand metrics in treeline forests. We reconstructed high-resolution photogrammetric point clouds and derived canopy height models for 10 study sites from NE Chukotka and SW Yakutia. Subsequently, we detected individual tree tops using a variable-window size local maximum filter and applied a marker-controlled watershed segmentation for the delineation of tree crowns. With this, we successfully detected 67.1% of the validation individuals. Simple linear regressions of observed and detected metrics show a better correlation (R2) and lower relative root mean square percentage error (RMSE%) for tree heights (mean R2 = 0.77, mean RMSE% = 18.46%) than for crown diameters (mean R2 = 0.46, mean RMSE% = 24.9%). The comparison between detected and observed tree height distributions revealed that our tree detection method was unable to representatively identify trees <2 m. Our results show that plot sizes for vegetation surveys in the tundra–taiga ecotone should be adapted to the forest structure and have a radius of >15–20 m to capture homogeneous and representative forest stands. Additionally, we identify sources of omission and commission errors and give recommendations for their mitigation. In summary, the efficiency of the used method depends on the complexity of the forest’s stand structure.
Water stable isotope ratios of tropical precipitation predominantly reflect moisture source and precipitation intensity. Trees can incorporate the isotopic signals into annual tree-ring cellulose records, permitting reconstruction of the temporal changes of hydroclimate over decades to millennia. This is especially valuable in the Himalayas where the understanding of monsoon dynamics is limited by the lack of a dense and representative observational network. We have analyzed tree ring delta O-18 records from two distinct physiographic sites along the upper Kali Gandaki valley in the central Nepal Himalayas, representing the wet High-Himalayas and the Trans-Himalayan dryland to the north. Empirical correlations and regression analyses were compared to an in-situ calibrated oxygen isotope fractionation model, exploring the relationships between tree ring delta O-18 and seasonal-mean variability of hydroclimatic forcing at the different locations. For this purpose, gridded precipitation data from the Asian rain gauge dataset APHRODITE, as well as high resolution onsite observations (relative humidity, air temperature, delta O-18 of precipitation and radial tree growth) were used. We found that two distinct sets of meteorological values, reflecting pre-monsoon and monsoon conditions, are needed to reproduce the measured tree ring delta O-18 values from the High-Himalayan site, but that a single set of monsoonal values performs best for the Trans-Himalayan site. We conclude that Trans-Himalayan trees capture long-term changes in strength of the Indian summer monsoon. In contrast, High-Himalayan tree ring delta(18)Orecords a more complex hydro-climatic signal reflecting both pre-monsoon and monsoon seasons with very contrasting isotopic signatures of precipitation. This difference in the two hydroclimatic proxy records offers an opportunity to reconstruct first-order hydroclimate conditions, such as local precipitation rates, and to gain new insights into monsoon timing and seasonal water source determination across the Himalayan orographic region. (C) 2019 Elsevier B.V. All rights reserved.
Lateral movements of alluvial river channels control the extent and reworking rates of alluvial fans, floodplains, deltas, and alluvial sections of bedrock rivers. These lateral movements can occur by gradual channel migration or by sudden changes in channel position (avulsions). Whereas models exist for rates of river avulsion, we lack a detailed understanding of the rates of lateral channel migration on the scale of a channel belt. In a two-step process, we develop here an expression for the lateral migration rate of braided channel systems in coarse, non-cohesive sediment. On the basis of photographic and topographic data from laboratory experiments of braided channels performed under constant external boundary conditions, we first explore the impact of autogenic variations of the channel-system geometry (i.e. channel-bank heights, water depths, channel-system width, and channel slope) on channel-migration rates. In agreement with theoretical expectations, we find that, under such constant boundary conditions, the laterally reworked volume of sediment is constant and lateral channel-migration rates scale inversely with the channel-bank height. Furthermore, when channel-bank heights are accounted for, lateral migration rates are independent of the remaining channel geometry parameters. These constraints allow us, in a second step, to derive two alternative expressions for lateral channel-migration rates under different boundary conditions using dimensional analysis. Fits of a compilation of laboratory experiments to these expressions suggest that, for a given channel bank-height, migration rates are strongly sensitive to water discharges and more weakly sensitive to sediment discharges. In addition, external perturbations, such as changes in sediment and water discharges or base level fall, can indirectly affect lateral channel-migration rates by modulating channel-bank heights.
The Atlantic Meridional Overturning Circulation (AMOC) is likely the most well-known system of ocean currents on Earth, redistributing heat, nutrients and carbon over a large part of the Earth’s surface and affecting global climate as a result. Due to enhanced freshwater fluxes into the subpolar North Atlantic as a response to global warming, the AMOC is expected, and may have already started, to weaken and these changes will likely have global impacts. It is therefore of considerable relevance to improve our understanding of past and future AMOC changes. My thesis tries to answer some of the open questions in this field by giving strong evidence that the AMOC has already weakened over the last century, by narrowing future projections of this slowdown and
by studying the impacts on global surface warming.
While there have been various studies trying to reconstruct the strength of the overturning circulation in the past, often based on model simulations in combination with observations (Jackson et al., 2016, Kanzow et al., 2010) or proxies (Frajka-Williams, 2015, Latif et al., 2006), the results so far, due to lack of direct measurements, have been inconclusive. In the first paper I build on previous work that links the anomalously low sea surface temperatures (SSTs) in the North Atlantic with the reduced meridional heat transport due to a weaker AMOC. Using the output of a high-resolution global climate model, I derive a characteristic spatial and seasonal SST fingerprint of an AMOC slowdown and an improved SST-based AMOC index. The same fingerprint is seen in
the observational SSTs since the late 19th Century, giving strong evidence that since then the AMOC has slowed down. In addition, the reconstruction of the historical overturning strength with the new AMOC index agrees well with and extends the results of earlier studies as well as the direct measurements from the RAPID project and shows a strong decline of the AMOC by about 15% (3±1 Sv) since the mid-20th Century (Caesar et al., 2018).
The reconstruction of the historical overturning strength with the AMOC index enables us to weight future AMOC projections based on their skill in modeling the historical AMOC as described in the second paper of this thesis (Olson et al., 2018). Using Bayesian model averaging we considerably narrow the projections of the CMIP5 ensemble to a decrease of -4.0 Sv and -6.8 Sv between the years 1960-1999 and 2060-2099 for the RCP4.5 and RCP8.5 emission scenarios, respectively. These values fit to, yet are at the lower end of, previously published estimates.
In the third paper I examine how the AMOC slowdown affects the global mean surface temperature (GMST) with a focus on how it will change the ocean heat uptake (OHC). Accounting for the effect of changes in the radiative forcing on the GMST, I test how AMOC variations correlate with the residual part of surface temperature changes in the past. I find that the correlation is positive which fits the understanding that the deep-water formation that is important in driving the AMOC cools the deep ocean and therefore warms the surface (Caesar et al., 2019). The future weakening of the overturning circulation could therefore delay global surface warming.
Due to nonlinear behavior and scale specific changes it can be difficult to study the dominant processes and modes that drive climate variability. In the fourth paper we develop and test a new technique based on the wavelet multiscale correlation (WMC) similarity measure to study climate variability on different temporal and spatial scales (Agarwal et al., 2018). In a fifth contribution to my thesis this method is applied to the observed sea surface temperatures. The results reconfirm well-known relations between SST anomalies such as the El Niño-Southern Oscillation (ENSO) and the Pacific Decadal Oscillation (PDO) on inter-annual and decadal timescales, respectively. They
furthermore give new insights into the characteristics and origins of long-range teleconnections, for example, that the teleconnection between ENSO and Indian Ocean dipole exist mainly between the northern part of the ENSO tongue and the equatorial Indian Ocean, and provides therefore valuable knowledge about the regions that are necessary to include when modeling regional climate variability at a certain scale (Agarwal et al., 2019).
In summary, my PhD thesis investigates past and future AMOC variability and its effects on global mean surface temperature by utilizing a combination of observational sea surface data and the output of historical and future climate model simulations from both the high-resolution CM2.6 model as well as the CMIP5 ensemble. It further includes the development and validation of a new method to study climate variability, that, applied to the observed sea surface temperatures, gives new insight about teleconnections in the Earth System. My findings provide evidence that the AMOC has already slowed down, will continue to do so in the future, and will impact the global mean temperature. Further impacts of an AMOC slowdown may include increased sea-level rise at the U.S. east coast (Ezer, 2015), heat extremes in Europe (Duchez et al., 2016) and increased storm activity in the North Atlantic region (Jackson et al., 2015), all of which have significant socio-economic implications.
Two principal groups of processes shape mass fluxes from and into a soil: vertical profile development and lateral soil redistribution. Periods having predominantly progressive soil forming processes (soil profile development) alternate with periods having predominantly regressive processes (erosion). As a result, short‐term soil redistribution – years to decades – can differ substantially from long‐term soil redistribution; i.e. centuries to millennia. However, the quantification of these processes is difficult and consequently their rates are poorly understood. To assess the competing roles of erosion and deposition we determined short‐ and long‐term soil redistribution rates in a formerly glaciated area of the Uckermark, northeast Germany. We compared short‐term erosion or accumulation rates using plutonium‐239 and ‐240 (239+240Pu) and long‐term rates using both in situ and meteoric cosmogenic beryllium‐10 (10Be). Three characteristic process domains have been analysed in detail: a flat landscape position having no erosion/deposition, an erosion‐dominated mid‐slope, and a deposition‐dominated lower‐slope site. We show that the short‐term mass erosion and accumulation rates are about one order of magnitude higher than long‐term redistribution rates. Both, in situ and meteoric 10Be provide comparable results. Depth functions, and therefore not only an average value of the topsoil, give the most meaningful rates. The long‐term soil redistribution rates were in the range of −2.1 t ha‐1 yr‐1 (erosion) and +0.26 t ha‐1 yr‐1 (accumulation) whereas the short‐term erosion rates indicated strong erosion of up to 25 t ha‐1 yr‐1 and accumulation of 7.6 t ha‐1 yr‐1. Our multi‐isotope method identifies periods of erosion and deposition, confirming the ‘time‐split approach’ of distinct different phases (progressive/regressive) in soil evolution. With such an approach, temporally‐changing processes can be disentangled, which allows the identification of both the dimensions of and the increase in soil erosion due to human influence
Ongoing and past biome transitions are generally assigned to climate and atmospheric changes (e.g. temperature, precipitation, CO2), but the major regional factors or factor combinations that drive vegetation change often remain unknown. Modelling studies applying ensemble runs can help to partition the effects of the different drivers. Such studies require careful validation with observational data. In this study, fossil pollen records from 741 sites in Europe, 728 sites in North America, and 418 sites in Asia (extracted from terrestrial archives including lake sediments) are used to reconstruct biomes at selected time slices between 40 cal ka BP (calibrated thousand years before present) and today. These results are used to validate Northern Hemisphere biome distributions (>30 degrees N) simulated by the biome model BIOME4 that has been forced with climate data simulated by a General Circulation model. Quantitative comparisons between pollen- and model-based results show a generally good fit at a broad spatial scale. Mismatches occur in central-arid Asia with a broader extent of grassland throughout the last 40 ka (likely due to the over-representation of Artemisia and Chenopodiaceae pollen) and in Europe with over-estimation of tundra at 0 cal ka BP (likely due to human impacts to some extent). Sensitivity analysis reveals that broad-scale biome changes follow the global signal of major postglacial temperature change, although the climatic variables vary in their regional and temporal importance. Temperature is the dominant variable in Europe and other rather maritime areas for biome changes between 21 and 14 ka, while precipitation is highly important in the arid inland regions of Asia and North America. The ecophysiological effect of changes in the atmospheric CO2-concentration has the highest impact during this transition than in other intervals. With respect to modern vegetation in the course of global warming, our findings imply that vegetation change in the Northern Hemisphere may be strongly limited by effective moisture changes, i.e. the combined effect of temperature and precipitation, particularly in inland areas. (C) 2019 Elsevier Ltd. All rights reserved.
Pollen-based quantitative land-cover reconstruction for northern Asia covering the last 40 ka cal BP
(2019)
We collected the available relative pollen productivity estimates (PPEs) for 27 major pollen taxa from Eurasia and applied them to estimate plant abundances during the last 40 ka cal BP (calibrated thousand years before present) using pollen counts from 203 fossil pollen records in northern Asia (north of 40 degrees N). These pollen records were organized into 42 site groups and regional mean plant abundances calculated using the REVEALS (Regional Estimates of Vegetation Abundance from Large Sites) model. Time-series clustering, constrained hierarchical clustering, and detrended canonical correspondence analysis were performed to investigate the regional pattern, time, and strength of vegetation changes, respectively. Reconstructed regional plant functional type (PFT) components for each site group are generally consistent with modern vegetation in that vegetation changes within the regions are characterized by minor changes in the abundance of PFTs rather than by an increase in new PFTs, particularly during the Holocene. We argue that pollen-based REVEALS estimates of plant abundances should be a more reliable reflection of the vegetation as pollen may overestimate the turnover, particularly when a high pollen producer invades areas dominated by low pollen producers. Comparisons with vegetation-independent climate records show that climate change is the primary factor driving land-cover changes at broad spatial and temporal scales. Vegetation changes in certain regions or periods, however, could not be explained by direct climate change, e.g. inland Siberia, where a sharp increase in evergreen conifer tree abundance occurred at ca. 7-8 ka cal BP despite an unchanging climate, potentially reflecting their response to complex climate-permafrost-fire-vegetation interactions and thus a possible long-term lagged climate response.
We investigate the inclusions hosted in peritectic garnet from metapelitic migmatites of the Kinzigite Formation (Ivrea Zone, NW Italy) to evaluate the starting composition of the anatectic melt and fluid regime during anatexis throughout the upper amphibolite facies, transition, and granulite facies zones. Inclusions have negative crystal shapes, sizes from 2 to 10 mu m and are regularly distributed in the core of the garnet. Microstructural and micro-Raman investigations indicate the presence of two types of inclusions: crystallized silicate melt inclusions (i.e., nanogranitoids, NI), and fluid inclusions (FI). Microstructural evidence suggests that FI and NI coexist in the same cluster and are primary (i.e., were trapped simultaneously during garnet growth). FI have similar compositions in the three zones and comprise variable proportions of CO2, CH4, and N-2, commonly with siderite, pyrophyllite, and kaolinite, suggesting a COHN composition of the trapped fluid. The mineral assemblage in the NI contains K-feldspar, plagioclase, quartz, biotite, muscovite, chlorite, graphite and, rarely, calcite. Polymorphs such as kumdykolite, cristobalite, tridymite, and less commonly kokchetavite, were also found. Rehomogenized NI from the different zones show that all the melts are leucogranitic but have slightly different compositions. In samples from the upper amphibolite facies, melts are less mafic (FeO + MgO = 2.0-3.4 wt%), contain 860-1700 ppm CO2 and reach the highest H2O contents (6.5-10 wt%). In the transition zone melts have intermediate H2O (4.8-8.5 wt%), CO2 (457-1534 ppm) and maficity (FeO + MgO = 2.3-3.9 wt%). In contrast, melts at granulite facies reach highest CaO, FeO + MgO (3.2-4.7 wt%), and CO2 (up to 2,400 ppm), with H2O contents comparable (5.4-8.3 wt%) to the other two zones. Our results represent the first clear evidence for carbonic fluid-present melting in the Ivrea Zone. Anatexis of metapelites occurred through muscovite and biotite breakdown melting in the presence of a COH fluid, in a situation of fluid-melt immiscibility. The fluid is assumed to have been internally derived, produced initially by devolatilization of hydrous silicates in the graphitic protolith, then as a result of oxidation of carbon by consumption of Fe3+-bearing biotite during melting. Variations in the compositions of the melts are interpreted to result from higher T of melting. The H2O contents of the melts throughout the three zones are higher than usually assumed for initial H2O contents of anatectic melts. The CO2 contents are highest at granulite facies, and show that carbon-contents of crustal magmas are not negligible at high T. The activity of H2O of the fluid dissolved in granitic melts decreases with increasing metamorphic grade. Carbonic fluid-present melting of the deep continental crust represents, together with hydrate-breakdown melting reactions, an important process in the origin of crustal anatectic granitoids.
We studied FeCO3 using Fe K-edge X-ray absorption near-edge structure (XANES) spectroscopy at pressures up to 54 GPa and temperatures above 2000 K. First-principles calculations of Fe at the K-edge in FeCO3 were performed to support the interpretation of the XANES spectra. The variation of iron absorption edge features with pressure and temperature in FeCO3 matches well with recently reported observations on FeCO3 at extreme conditions, and provides new insight into the stability of Fe-carbonates in Earth's mantle. Here we show that at conditions of the mid-lower mantle, ~50 GPa and ~2200 K, FeCO3 melts and partially decomposes to high-pressure Fe3O4. Carbon (diamond) and oxygen are also inferred products of the reaction. We constrained the thermodynamic phase boundary between crystalline FeCO3 and melt to be at 51(1) GPa and ~1850 K. We observe that at 54(1) GPa, temperature-induced spin crossover of Fe2+ takes place from low to high spin such that at 1735(100) K, all iron in FeCO3 is in the high-spin state. A comparison between experiment and theory provides a more detailed understanding of FeCO3 decomposition observed in X-ray absorption spectra and helps to explain spectral changes due to pressure-induced spin crossover in FeCO3 at ambient temperature.
In response to collision and convergence between India and Asia during the Cenozoic, convergence took place between the Pamir and South Tian Shan. Here we present new detrital zircon U-Pb ages coupled with conglomerate clast counting and sedimentary data from the late Cenozoic Wuheshalu section in the convergence zone, to shed light on the convergence process of the Pamir and South Tian Shan. Large Triassic zircon U-Pb age populations in all seven samples suggest that Triassic igneous rocks from the North Pamir were the major source area for the late Cenozoic Wuheshalu section. In the Miocene, large populations of the North Pamir component supports rapid exhumation in the North Pamir and suggest that topography already existed there since the early Miocene. Exhumation of the South Tian Shan was relatively less important in the Miocene and its detritus could only reach a limited area in the foreland area. Gradually increasing sediment loading and convergence of the Pamir and South Tian Shan caused rapid subsidence in the convergence area. Since ca. 6-5.3 Ma, the combination of a major North Pamir component and a minor South Tian Shan component at the Wuheshalu section is consistent with active deformation of the South Tian Shan and the North Pamir. During deposition of the upper Atushi Formation, a larger proportion of North Pamir-derived sediments was deposited in the Wuheshalu section, maybe because faulting and northward propagation of the North Pamir caused northward displacement of the depocenter to north of the Wuheshalu section.
The structure and organization of river networks has been used for decades to investigate the influence of climate and tectonics on landscapes. The majority of these studies either analyze rivers in profile view by extracting channel steepness or calculate planform metrics such as drainage density. However, these techniques rely on the assumption of homogeneity: that intrinsic and external factors are spatially or temporally invariant over the measured profile. This assumption is violated for the majority of Earth's landscapes, where variations in uplift rate, rock strength, climate, and geomorphic process are almost ubiquitous. We propose a method for classifying river profiles to identify landscape regions with similar characteristics by adapting hierarchical clustering algorithms developed for time series data. We first test our clustering on two landscape evolution scenarios and find that we can successfully cluster regions with different erodibility and detect the transient response to sudden base level fall. We then test our method in two real landscapes: first in Bitterroot National Forest, Idaho, where we demonstrate that our method can detect transient incision waves and the topographic signature of fluvial and debris flow process regimes; and second, on Santa Cruz Island, California, where our technique identifies spatial patterns in lithology not detectable through normalized channel steepness analysis. By calculating channel steepness separately for each cluster, our method allows the extraction of more reliable steepness metrics than if calculated for the landscape as a whole. These examples demonstrate the method's ability to disentangle fluvial morphology in complex lithological and tectonic settings.
Climate change is affecting the rate of carbon cycling, particularly in the Arctic. Permafrost degradation through deeper thaw and physical disturbances results in the release of carbon dioxide and methane to the atmosphere and to an increase in lateral dissolved organic matter (DOM) fluxes. Whereas riverine DOM fluxes of the large Arctic rivers are well assessed, knowledge is limited with regard to small catchments that cover more than 40% of the Arctic drainage basin. Here, we use absorption measurements to characterize changes in DOM quantity and quality in a low Arctic (Herschel Island, Yukon, Canada) and a high Arctic (Cape Bounty, Melville Island, Nunavut, Canada) setting with regard to geographical differences, impacts of permafrost degradation, and rainfall events. We find that DOM quantity and quality is controlled by differences in vegetation cover and soil organic carbon content (SOCC). The low Arctic site has higher SOCC and greater abundance of plant material resulting in higher chromophoric dissolved organic matter (cDOM) and dissolved organic carbon (DOC) than in the high Arctic. DOC concentration and cDOM in surface waters at both sites show strong linear relationships similar to the one for the great Arctic rivers. We used the optical characteristics of DOM such as cDOM absorption, specific ultraviolet absorbance (SUVA), ultraviolet (UV) spectral slopes (S275-295), and slope ratio (SR) for assessing quality changes downstream, at base flow and storm flow conditions, and in relation to permafrost disturbance. DOM in streams at both sites demonstrated optical signatures indicative of photodegradation downstream processes, even over short distances of 2000 m. Flow pathways and the connected hydrological residence time control DOM quality. Deeper flow pathways allow the export of permafrost-derived DOM (i.e. from deeper in the active layer), whereas shallow pathways with shorter residence times lead to the export of fresh surface- and near-surface-derived DOM. Compared to the large Arctic rivers, DOM quality exported from the small catchments studied here is much fresher and therefore prone to degradation. Assessing optical properties of DOM and linking them to catchment properties will be a useful tool for understanding changing DOM fluxes and quality at a pan-Arctic scale.
This dissertation combines field and geochemical observations and analyses with numerical modeling to understand the formation of vein-hosted Sn-W ore in the Panasqueira deposit of Portugal, which is among the ten largest worldwide. The deposit is located above a granite body that is altered by magmatic-hydrothermal fluids in its upper part (greisen). These fluids are thought to be the source of metals, but that was still under debate. The goal of this study is to determine the composition and temperature of hydrothermal fluids at Panasqueira, and with that information to construct a numerical model of the hydrothermal system. The focus is on analysis of the minerals tourmaline and white mica, which formed during mineralization and are widespread throughout the deposit. Tourmaline occurs mainly in alteration zones around mineralized veins and is less abundant in the vein margins. White mica is more widespread. It is abundant in vein margins as well as alteration zones, and also occurs in the granite greisen. The laboratory work involved in-situ microanalysis of major- and trace elements in tourmaline and white mica, and boron-isotope analysis in both minerals by secondary ion mass spectrometry (SIMS).
The boron-isotope composition of tourmaline and white mica suggests a magmatic source. Comparison of hydrothermally-altered and unaltered rocks from drill cores shows that the ore metals (W, Sn, Cu, and Zn) and As, F, Li, Rb, and Cs were introduced during the alteration. Most of these elements are also enriched in tourmaline and mica, which confirms their potential value as exploration guides to Sn-W ores elsewhere.
The thermal evolution of the hydrothermal system was estimated by B-isotope exchange thermometry and the Ti-in-quartz method. Both methods yielded similar temperatures for the early hydrothermal phase: 430° to 460°C for B-isotopes and 503° ± 24°C for Ti-in-quartz. Mineral pairs from a late fault zone yield significantly lower median temperatures of 250°C. The combined results of thermometry with variations in chemical and B-isotope composition of tourmaline and mica suggest that a similar magmatic-hydrothermal fluid was active at all stages of mineralization. Mineralization in the late stage shows the same B-isotope composition as in the main stage despite a ca. 250°C cooling, which supports a multiple injection model of magmatic-hydrothermal fluids.
Two-dimensional numerical simulations of convection in a multiphase NaCl hydrothermal system were conducted: (a) in order to test a new approach (lower dimensional elements) for flow through fractures and faults and (b) in order to identify conditions for horizontal fluid flow as observed in the flat-lying veins at Panasqueira. The results show that fluid flow over an intrusion (heat and fluid source) develops a horizontal component if there is sufficient fracture connectivity. Late, steep fault zones have been identified in the deposit area, which locally contain low-temperature Zn-Pb mineralization. The model results confirm that the presence of subvertical faults with enhanced permeability play a crucial role in the ascent of magmatic fluids to the surface and the recharge of meteoric waters. Finally, our model results suggest that recharge of meteoric fluids and mixing processes may be important at later stages, while flow of magmatic fluids dominate the early stages of the hydrothermal fluid circulation.
Although quantitative isotope data from speleothems has been used to evaluate isotope-enabled model simulations, currently no consensus exists regarding the most appropriate methodology through which to achieve this. A number of modelling groups will be running isotope-enabled palaeoclimate simulations in the framework of the Coupled Model Intercomparison Project Phase 6, so it is timely to evaluate different approaches to using the speleothem data for data–model comparisons. Here, we illustrate this using 456 globally distributed speleothem δ18O records from an updated version of the Speleothem Isotopes Synthesis and Analysis (SISAL) database and palaeoclimate simulations generated using the ECHAM5-wiso isotope-enabled atmospheric circulation model. We show that the SISAL records reproduce the first-order spatial patterns of isotopic variability in the modern day, strongly supporting the application of this dataset for evaluating model-derived isotope variability into the past. However, the discontinuous nature of many speleothem records complicates the process of procuring large numbers of records if data–model comparisons are made using the traditional approach of comparing anomalies between a control period and a given palaeoclimate experiment. To circumvent this issue, we illustrate techniques through which the absolute isotope values during any time period could be used for model evaluation. Specifically, we show that speleothem isotope records allow an assessment of a model's ability to simulate spatial isotopic trends. Our analyses provide a protocol for using speleothem isotope data for model evaluation, including screening the observations to take into account the impact of speleothem mineralogy on δ18O values, the optimum period for the modern observational baseline and the selection of an appropriate time window for creating means of the isotope data for palaeo-time-slices.
We present a new algorithm for solving the common problem of flow trapped in closed depressions within digital elevation models, as encountered in many applications relying on flow routing. Unlike other approaches (e.g., the Priority-Flood depression filling algorithm), this solution is based on the explicit computation of the flow paths both within and across the depressions through the construction of a graph connecting together all adjacent drainage basins. Although this represents many operations, a linear time complexity can be reached for the whole computation, making it very efficient. Compared to the most optimized solutions proposed so far, we show that this algorithm of flow path enforcement yields the best performance when used in landscape evolution models. In addition to its efficiency, our proposed method also has the advantage of letting the user choose among different strategies of flow path enforcement within the depressions (i.e., filling vs. carving). Furthermore, the computed graph of basins is a generic structure that has the potential to be reused for solving other problems as well, such as the simulation of erosion. This sequential algorithm may be helpful for those who need to, e.g., process digital elevation models of moderate size on single computers or run batches of simulations as part of an inference study.
Costs of adaptation in the developing world have been mostly equated to those of climate proofing infrastructure under the assumption of unconstrained knowledge and planning capacities. To correct this, we introduce a cost-scaling methodology estimating sectoral investments to enhance the knowledge and planning capacities of countries based on an empirical collection of 385 climate-related projects. We estimate that circa 9.2 billion USD are required for financing knowledge and planning activities in developing countries in 2015. The agricultural and water sectors demand the higher investments ? 3.8 and 3.5 billion USD, respectively. Average investments between 2015 and 2050 are projected at 7 billion USD per year ? the largest fraction of which (4 billion) in Africa. Investments in this study were found to constitute approximately 40%, 20?60% and 5?15% of previous cost estimates to climate-proof infrastructure in the agricultural, water, and coastal sectors, respectively. The effort to finance the knowledge and planning capacities in developing countries is therefore not marginal relative to the costs of adapting infrastructure. The findings underline the potential of using empirical collections of climate-related projects for adaptation cost assessments as complementary to process and economic models.
Stochastic modeling is a common practice for modeling uncertainty in hydrogeology. In stochastic modeling, aquifer properties are characterized by their probability density functions (PDFs). The Bayesian approach for inverse modeling is often used to assimilate information from field measurements collected at a site into properties’ posterior PDFs. This necessitates the definition of a prior PDF, characterizing the knowledge of hydrological properties before undertaking any investigation at the site, and usually coming from previous studies at similar sites. In this paper, we introduce a Bayesian hierarchical algorithm capable of assimilating various information–like point measurements, bounds and moments–into a single, informative PDF that we call ex-situ prior. This informative PDF summarizes the ex-situ information available about a hydrogeological parameter at a site of interest, which can then be used as a prior PDF in future studies at the site. We demonstrate the behavior of the algorithm on several synthetic case studies, compare it to other methods described in the literature, and illustrate the approach by applying it to a public open-access hydrogeological dataset.
Environmental factors shape the spatial distribution and dynamics of populations. Understanding how these factors interact with movement behavior is critical for efficient conservation, in particular for migratory species. Adult female green sea turtles, Chelonia mydas, migrate between foraging and nesting sites that are generally separated by thousands of kilometers. As an emblematic endangered species, green turtles have been intensively studied, with a focus on nesting, migration, and foraging. Nevertheless, few attempts integrated these behaviors and their trade‐offs by considering the spatial configurations of foraging and nesting grounds as well as environmental heterogeneity like oceanic currents and food distribution. We developed an individual‐based model to investigate the impact of local environmental conditions on emerging migratory corridors and reproductive output and to thereby identify conservation priority sites. The model integrates movement, nesting, and foraging behavior. Despite being largely conceptual, the model captured realistic movement patterns which confirm field studies. The spatial distribution of migratory corridors and foraging hot spots was mostly constrained by features of the regional landscape, such as nesting site locations, distribution of feeding patches, and oceanic currents. These constraints also explained the mixing patterns in regional forager communities. By implementing alternative decision strategies of the turtles, we found that foraging site fidelity and nesting investment, two characteristics of green turtles' biology, are favorable strategies under unpredictable environmental conditions affecting their habitats. Based on our results, we propose specific guidelines for the regional conservation of green turtles as well as future research suggestions advancing spatial ecology of sea turtles. Being implemented in an easy to learn open‐source software, our model can coevolve with the collection and analysis of new data on energy budget and movement into a generic tool for sea turtle research and conservation. Our modeling approach could also be useful for supporting the conservation of other migratory marine animals.
Geomorphic strain markers accumulating the effects of many earthquake cycles help to constrain the mechanical behaviour of continental rift systems as well as the related seismic hazards. In the Corinth Rift (Greece), the unique record of onshore and offshore markers of Pleistocene similar to 100-ka climate cycles provides an outstanding possibility to constrain rift mechanics over a range of timescales. Here we use high-resolution topography to analyse the 3D geometry of a sequence of Pleistocene emerged marine terraces associated with flexural rift-flank uplift. We integrate this onshore dataset with offshore seismic data to provide a synoptic view of the flexural deformation across the rift. This allows us to derive an average slip rate of 4.5-9.0 mm.yr(-1) on the master fault over the past similar to 610 ka and an uplift/ subsidence ratio of 1:1.1-2.4. We reproduce the observed flexure patterns, using 3 and 5-layered lithospheric scale finite element models. Modelling results imply that the observed elastic flexure is produced by coseismic slip along 40-60 degrees planar normal faults in the elastic upper crust, followed by postseismic viscous relaxation occurring within the basal lower crust or upper mantle. We suggest that such a mechanism may typify rapid localised extension of continental lithosphere.
The Olorgesailie Drilling Project and the related Hominin Sites and Paleolakes Drilling Project in East Africa were initiated to test hypotheses and models linking environmental change to hominin evolution by drilling lake basin sediments adjacent to important archeological and paleoanthropological sites. Drill core OL012-1A recovered 139 m of sedimentary and volcaniclastic strata from the Koora paleolake basin, southern Kenya Rift, providing the opportunity to compare paleoenvironmental influences over the past million years with the parallel record exposed at the nearby Olorgesailie archeological site. To refine our ability to link core-to-outcrop paleoenvironmental records, we institute here a methodological framework for deriving a robust age model for the complex lithostratigraphy of OL012-1A. Firstly, chronostratigraphic control points for the core were established based on 4 Ar/39Ar ages from intercalated tephra deposits and a basal trachyte flow, as well as the stratigraphic position of the Brunhes-Matuyama geomagnetic reversal. This dataset was combined with the position and duration of paleosols, and analyzed using a new Bayesian algorithm for high-resolution age-depth modeling of hiatus-bearing stratigraphic sections. This model addresses three important aspects relevant to highly dynamic, nonlinear depositional environments: 1) correcting for variable rates of deposition, 2) accommodating hiatuses, and 3) quantifying realistic age uncertainty with centimetric resolution. Our method is applicable to typical depositional systems in extensional rifts as well as to drill cores from other dynamic terrestrial or aquatic environments. We use the core age model and lithostratigraphy to examine the inter connectivity of the Koora Basin to adjacent areas and sources of volcanism. (C) 2019 Elsevier Ltd. All rights reserved.
The Baringo-Tugen-Barsemoi 2013 drillcore (BTB13), acquired as part of the Hominin Sites and Paleolakes Drilling Project, recovered 228 m of fluviolacustrine sedimentary rocks and tuffs spanning a similar to 3.29-2.56 Ma interval of the highly fossiliferous and hominin-bearing Chemeron Formation, Tugen Hills, Kenya. Here we present a Bayesian stratigraphic age model for the core employing chronostratigraphic control points derived from Ar-40/Ar-39 dating of tuffs from core and outcrop, Ar-40/Ar-39 age calibration of related outcrop diatomaceous units, and core magnetostratigraphy. The age model reveals three main intervals with distinct sediment accumulation rates: an early rapid phase from 3.2 to 2.9 Ma; a relatively slow phase from 2.9 to 2.7 Ma; and the highest rate of accumulation from 2.7 to 2.6 Ma. The intervals of rapid accumulation correspond to periods of high Earth orbital eccentricity, whereas the slow accumulation interval corresponds to low eccentricity at 2.9-2.7 Ma, suggesting that astronomically mediated climate processes may be responsible for the observed changes in sediment accumulation rate. Lacustrine transgression-regression events, as delineated using sequence stratigraphy, dominantly operate on precession scale, particularly within the high eccentricity periods. A set of erosively based fluvial conglomerates correspond to the 2.9-2.7 Ma interval, which could be related to either the depositional response to low eccentricity or to the development of unconformities due to local tectonic activity. Age calibration of core magnetic susceptibility and gamma density logs indicates a close temporal correspondence between a shift from high- to low-frequency signal variability at similar to 3 Ma, approximately coincident the end of the mid-Piacenzian Warm Period, and the beginning of the cooling of world climate leading to the initiation of Northern Hemispheric glaciation c. 2.7 Ma. BTB13 and the Baringo Basin records may thus provide evidence of a connection between high-latitude glaciation and equatorial terrestrial climate toward the end of the Pliocene.
The Central Andes host large reserves of base and precious metals. The region represented, in 2017, an important part of the worldwide mining activity. Three principal types of deposits have been identified and studied: 1) porphyry type deposits extending from central Chile and Argentina to Bolivia, and Northern Peru, 2) iron oxide-copper-gold (IOCG) deposits, extending from central Peru to central Chile, and 3) epithermal tin polymetallic deposits extending from Southern Peru to Northern Argentina, which compose a large part of the deposits of the Bolivian Tin Belt (BTB). Deposits in the BTB can be divided into two major types: (1) tin-tungsten-zinc pluton-related polymetallic deposits, and (2) tin-silver-lead-zinc epithermal polymetallic vein deposits.
Mina Pirquitas is a tin-silver-lead-zinc epithermal polymetallic vein deposit, located in north-west Argentina, that used to be one of the most important tin-silver producing mine of the country. It was interpreted to be part of the BTB and it shares similar mineral associations with southern pluton related BTB epithermal deposits. Two major mineralization events related to three pulses of magmatic fluids mixed with meteoric water have been identified. The first event can be divided in two stages: 1) stage I-1 with quartz, pyrite, and cassiterite precipitating from fluids between 233 and 370 °C and salinity between 0 and 7.5 wt%, corresponding to a first pulse of fluids, and 2) stage I-2 with sphalerite and tin-silver-lead-antimony sulfosalts precipitating from fluids between 213 and 274 °C with salinity up to 10.6 wt%, corresponding to a new pulse of magmatic fluids in the hydrothermal system. The mineralization event II deposited the richest silver ores at Pirquitas. Event II fluids temperatures and salinities range between 190 and 252 °C and between 0.9 and 4.3 wt% respectively. This corresponds to the waning supply of magmatic fluids. Noble gas isotopic compositions and concentrations in ore-hosted fluid inclusions demonstrate a significant contribution of magmatic fluids to the Pirquitas mineralization although no intrusive rocks are exposed in the mine area.
Lead and sulfur isotopic measurements on ore minerals show that Pirquitas shares a similar signature with southern pluton related polymetallic deposits in the BTB. Furthermore, the major part of the sulfur isotopic values of sulfide and sulfosalt minerals from Pirquitas ranges in the field for sulfur derived from igneous rocks. This suggests that the main contribution of sulfur to the hydrothermal system at Pirquitas is likely to be magma-derived. The precise age of the deposit is still unknown but the results of wolframite dating of 2.9 ± 9.1 Ma and local structural observations suggest that the late mineralization event is younger than 12 Ma.
With the Late Cretaceous onset of Africa-Iberia-Europe convergence Central Europe experienced a pulse of intraplate shortening lasting some 15-20 Myr. This deformation event documents area-wide deviatoric compression of Europe and has been interpreted as a far-field response to Africa-Iberia-Europe convergence. However, the factors that governed the compression of Europe and conditioned the transient character of the deformation event have remained unclear. Based on mechanical considerations, numerical simulations, and geological reconstructions, we examine how the dynamics of intraplate deformation were governed by the formation of a convergent plate boundary fault between Iberia and Europe. During the Late Cretaceous, plate convergence was accommodated by the inversion of a young hyperextended rift system separating Iberia from Europe. Our analysis shows that the strength of the lithosphere beneath this rift was initially sufficient to transmit large compressive stresses far into Europe, though the lithosphere beneath the rift was thinned and thermally weakened. Continued convergence forced the formation of the plate boundary fault between Iberia and Europe. The fault evolved progressively and constituted a lithospheric-scale structure at the southern margin of Europe that weakened rheologically. This development caused a decrease in mechanical coupling between Iberia and Europe and a reduction of compressional far field stresses, which eventually terminated intraplate deformation in Central Europe. Taken together, our findings suggest that the Late Cretaceous intraplate deformation event records a high force transient that relates to the earliest strength evolution of a lithospheric-scale plate boundary fault.
The 2015 Paris Agreement (PA) has been widely hailed as a diplomatic triumph and a breakthrough in global climate cooperation. However, it is commonly accepted that the PA's collective goal—keeping global warming “well below” 2°C above preindustrial levels—remains ambitious. Making matters even more challenging, in 2017, global CO2 emissions resumed growth after 3 years of near standstill. In 2018, this growth accelerated. It is therefore extremely important that the PA's institutional architecture meet expectations concerning its ability to induce member countries to promise and deliver emissions reductions. This study offers a review of the rapidly growing literature on the PA, to assess its strengths and weaknesses, its significance, and its prospects. We focus on evaluations of its institutional structure and its ability to induce member countries to implement policies. We frame the issues as a trilemma: the challenge of simultaneously satisfying all three main conditions for effectiveness—broad participation, deep commitments, and satisfactory compliance rates. Based on our review, we conclude that the key challenge for the PA will likely be to facilitate sufficiently fast ratcheting‐up of nationally determined contributions, while keeping compliance rates high.
Stable carbon isotopes of sediment organic matter (delta C-13(OM)) are widely applied in paleoenvironmental studies. Interpretations of delta C-13(OM), however, remain challenging and factors that influence delta C-13(OM) may not apply across all lakes. Common explanations for stratigraphic shifts in delta C-13(OM) include changes in lake productivity or changes in inputs of allochthonous OM. We investigated the influence of different oxygen conditions (oxic versus anoxic) on the delta C-13(OM) values in the sediments of Lake Tiefer See. We analysed (1) a long sediment core from the deepest part of the lake, (2) two short, sediment-water interface cores from shallower water depths, and (3) OM in the water column, i.e. from sediment traps. Fresh OM throughout the entire water column showed a relatively constant delta C-13(OM) value of approximately -30.5 parts per thousand. Similar values, about -31 parts per thousand, were obtained for well-varved sediments in both the long and short, sediment-water interface cores. In contrast, delta C-13(OM) values from non-varved sediments in all cores were significantly less negative (-29 parts per thousand). The delta C-13(OM) values in the sediment-water interface cores from different water depths differ for sediments of the same age, if oxygen conditions at the time of deposition were different at these sites, as suggested by the state of varve preservation. Sediments deposited from AD 1924 to 1980 at 62m water depth are varved and exhibit delta C-13(OM) values around -31 parts per thousand, whereas sediments of the same age in the core from 35m water depth are not varved and show less negative delta C-13(OM) values of about -29 parts per thousand. The relation between varve occurrence and delta C-13(OM) values suggests that delta C-13(OM) is associated with oxygen conditions because varve preservation depends on hypolimnetic anoxia. A mechanism that likely influences delta C-13(OM) is selective degradation of OM under oxic conditions, such that organic components with more negative delta C-13(OM) are preferably decomposed, leading to less negative delta C-13(OM) values in the remaining, undegraded OM pool. Greater decomposition of OM in non-varved sediments is supported by lower TOC concentrations in these deposits (similar to 5%) compared to well-varved sediments (similar to 15%). Even in lakes that display small variations in productivity and terrestrial OM input through time, large spatial and temporal differences in hypolimnetic oxygen concentrations may be an important factor controlling sediment delta C-13(OM).
Localization processes in the viscous lower crust generate ductile shear zones over a broad range of scales affecting long‐term lithosphere deformation and the mechanical response of faults during the seismic cycle. Here we use centimeter‐scale numerical models in order to gain detailed insight into the processes involved in strain localization and rheological weakening in viscously deforming rocks. Our 2‐D Cartesian models are benchmarked to high‐temperature and high‐pressure torsion experiments on Carrara marble samples containing a single weak Solnhofen limestone inclusion. The models successfully reproduce bulk stress‐strain transients and final strain distributions observed in the experiments by applying a simple, first‐order softening law that mimics rheological weakening. We find that local stress concentrations forming at the inclusion tips initiate strain localization inside the host matrix. At the tip of the propagating shear zone, weakening occurs within a process zone, which expands with time from the inclusion tips toward the matrix. Rheological weakening is a precondition for shear zone localization, and the width of this shear zone is found to be controlled by the degree of softening. Introducing a second softening step at elevated strain, a high strain layer develops inside the localized shear zone, analogous to the formation of ultramylonite bands in mylonites. These results elucidate the transient evolution of stress and strain rate during inception and maturation of ductile shear zones.
The Upper Cretaceous (Campanian-Maastrichtian) bioclastic wedge of the Orfento Formation in the Montagna della Maiella, Italy, is compared to newly discovered contourite drifts in the Maldives. Like the drift deposits in the Maldives, the Orfento Formation fills a channel and builds a Miocene delta-shaped and mounded sedimentary body in the basin that is similar in size to the approximately 350 km(2) large coarse-grained bioclastic Miocene delta drifts in the Maldives. The composition of the bioclastic wedge of the Orfento Formation is also exclusively bioclastic debris sourced from the shallow-water areas and reworked clasts of the Orfento Formation itself. In the near mud-free succession, age-diagnostic fossils are sparse. The depositional textures vary from wackestone to float-rudstone and breccia/conglomerates, but rocks with grainstone and rudstone textures are the most common facies. In the channel, lensoid convex-upward breccias, cross-cutting channelized beds and thick grainstone lobes with abundant scours indicate alternating erosion and deposition from a high-energy current. In the basin, the mounded sedimentary body contains lobes with a divergent progradational geometry. The lobes are built by decametre thick composite megabeds consisting of sigmoidal clinoforms that typically have a channelized topset, a grainy foreset and a fine-grained bottomset with abundant irregular angular clasts. Up to 30 m thick channels filled with intraformational breccias and coarse grainstones pinch out downslope between the megabeds. In the distal portion of the wedge, stacked grainstone beds with foresets and reworked intraclasts document continuous sediment reworking and migration. The bioclastic wedge of the Orfento Formation has been variously interpreted as a succession of sea-level controlled slope deposits, a shoaling shoreface complex, or a carbonate tidal delta. Current-controlled delta drifts in the Maldives, however, offer a new interpretation because of their similarity in architecture and composition. These similarities include: (i) a feeder channel opening into the basin; (ii) an excavation moat at the exit of the channel; (iii) an overall mounded geometry with an apex that is in shallower water depth than the source channel; (iv) progradation of stacked lobes; (v) channels that pinch out in a basinward direction; and (vi) smaller channelized intervals that are arranged in a radial pattern. As a result, the Upper Cretaceous (Campanian-Maastrichtian) bioclastic wedge of the Orfento Formation in the Montagna della Maiella, Italy, is here interpreted as a carbonate delta drift.
Anthropogenic pressures increasingly alter natural systems. Therefore, understanding the resilience of agent-based complex systems such as ecosystems, i.e. their ability to absorb these pressures and sustain their functioning and services, is a major challenge. However, the mechanisms underlying resilience are still poorly understood. A main reason for this is the multidimensionality of both resilience, embracing the three fundamental stability properties recovery, resistance and persistence, and of the specific situations for which stability properties can be assessed. Agent-based models (ABM) complement empirical research which is, for logistic reasons, limited in coping with these multiple dimensions. Besides their ability to integrate multidimensionality through extensive manipulation in a fully controlled system, ABMs can capture the emergence of system resilience from individual interactions and feedbacks across different levels of organization. To assess the extent to which this potential of ABMs has already been exploited, we reviewed the state of the art in exploring resilience and its multidimensionality in ecological and socio-ecological systems with ABMs. We found that the potential of ABMs is not utilized in most models, as they typically focus on a single dimension of resilience by using variability as a proxy for persistence, and are limited to one reference state, disturbance type and scale. Moreover, only few studies explicitly test the ability of different mechanisms to support resilience. To overcome these limitations, we recommend to simultaneously assess multiple stability properties for different situations and under consideration of the mechanisms that are hypothesised to render a system resilient. This will help us to better exploit the potential of ABMs to understand and quantify resilience mechanisms, and hence support solving real-world problems related to the resilience of agent-based complex systems.
Geysers are hot springs whose frequency of water eruptions remain poorly understood. We set up a local broadband seismic network for 1 year at Strokkur geyser, Iceland, and developed an unprecedented catalog of 73,466 eruptions. We detected 50,135 single eruptions but find that the geyser is also characterized by sets of up to six eruptions in quick succession. The number of single to sextuple eruptions exponentially decreased, while the mean waiting time after an eruption linearly increased (3.7 to 16.4 min). While secondary eruptions within double to sextuple eruptions have a smaller mean seismic amplitude, the amplitude of the first eruption is comparable for all eruption types. We statistically model the eruption frequency assuming discharges proportional to the eruption multiplicity and a constant probability for subsequent events within a multituple eruption. The waiting time after an eruption is predictable but not the type or amplitude of the next one. <br /> Plain Language Summary Geysers are springs that often erupt in hot water fountains. They erupt more often than volcanoes but are quite similar. Nevertheless, it is poorly understood how often volcanoes and also geysers erupt. We created a list of 73,466 eruption times of Strokkur geyser, Iceland, from 1 year of seismic data. The geyser erupted one to six times in quick succession. We found 50,135 single eruptions but only 1 sextuple eruption, while the mean waiting time increased from 3.7 min after single eruptions to 16.4 min after sextuple eruptions. Mean amplitudes of each eruption type were higher for single eruptions, but all first eruptions in a succession were similar in height. Assuming a constant heat inflow at depth, we can predict the waiting time after an eruption but not the type or amplitude of the next one.
In this study, a low-cost unmanned aerial vehicle was used to obtain multi-spectral high-resolution imagery (1.4 cmpx(-1)) from2 microcatchments (3.3 ha) with burned Mediterranean shrubland and pine forests. This imagery was used to calculate the blue normalized differential vegetation index and to generate digital elevation models for estimating the sediment connectivity index. Both indices enabled an integrated approach for deciphering how hydrological and sediment connectivity interact with vegetation as well as soil conservation structures. The application of spatial analysis improves our understanding of the feedback between biological and geomorphological processes. Local spatial data analysis established a significant link between local geomorphological and biological factors, enabling a precise identification of homogeneous areas at micro-catchment scale and the minimal size of vegetation units reacting to geomorphology as natural groups at plot-scale where management strategies and efforts should be applied. Establishing this local relationship between sediment connectivity and vegetation patterns through new and interdisciplinary methodologies represents a new strategy for the assessment of ecosystem dynamics and management.
Several hydraulic fracturing tests were performed in boreholes located in central Hungary in order to determine the in-situ stress for a geological site investigation. At a depth of about 540m, the observed pressure versus time curves in mica schist with low dip angle foliation shows atypical pressure versus time results. After each pressurization cycle, the fracture breakdown pressure in the first fracturing cycle is lower than the refracturing or reopening pressure in the subsequent pressurizations. It is assumed that the viscosity of the drilling mud and observed foliation of the mica schist have a significant influence on the pressure values. In order to study this problem, numerical modeling was performed using the distinct element code particle flow code, which has been proven to be a valuable tool to investigate rock engineering problems such as hydraulic fracturing. The two-dimensional version of the code applied in this study can simulate hydro-mechanically coupled fluid flow in crystalline rock with low porosity and pre-existing fractures. In this study, the effect of foliation angle and fluid viscosity on the peak pressure is tested. The atypical characteristics of the pressure behaviour are interpreted so that mud with higher viscosity penetrates the sub-horizontal foliation plane, blocks the plane of weakness and makes the partly opened fracture tight and increase the pore pressure which decreases slowly with time. We see this viscous blocking effect as one explanation for the observed increase in fracture reopening pressure in subsequent pressurization cycles.
Obtaining representative hydrometric values is essential for characterizing extreme events, hydrological dynamics and detecting possible changes on the long-term hydrology. Reliability of streamflow data requires a temporal continuity and a maintenance of the gauging stations, which data are affected by epistemic and random sources of error. An assessment of discharge meterings' and stage-discharge rating curves' uncertainties were carried out by comparing the accuracy of the measuring instruments of two different hydrometric networks (i.e., one analogical and one digital) established in the same river location at the Mediterranean island of Mallorca. Furthermore, the effects of such uncertainties were assessed on the hydrological dynamics, considering the significant global change impacts beset this island. Evaluation was developed at four representative gauging stations of the hydrographic network with analogic (≈40 years) and digital (≈10 years) data series. The study revealed that the largest source of uncertainty in the analogical (28 to 274%) and in the digital (17–37%) networks were the stage-discharge rating curves. Their impact on the water resources was also evaluated at the event and annual scales, resulting in an average difference of water yields of 183% and 142% respectively. Such improvement on the comprehension of hydrometric networks uncertainties will dramatically benefit the interpretation of the long-term streamflow by providing better insights into the hydrologic and flood hazard planning, management and modelling.
Thrombolite and stromatolite habitats are becoming increasingly recognized as important refuges for invertebrates during Phanerozoic Oceanic Anoxic Events (OAEs); it is posited that oxygenic photosynthesis by cyanobacteria in these microbialites provided a refuge from anoxic conditions (i.e., the "microbialite refuge" hypothesis). Here, we test this hypothesis by investigating the distribution of ~34, 500 benthic invertebrate fossils found in ~100 samples from a microbialite succession that developed following the latest Permian mass extinction event on the Great Bank of Guizhou (South China), representing microbial (stromatolites and thrombolites) and non-microbial facies. The stromatolites were the least taxonomically diverse facies, and the thrombolites also recorded significantly lower diversities when compared to the non-microbial facies. Based on the distribution and ornamentation of the bioclasts within the thrombolites and stromatolites, the bioclasts are inferred to have been transported and concentrated in the non-microbial fabrics, that is, cavities around the microbial framework. Therefore, many of the identified metazoans from the post-extinction microbialites are not observed to have been living within a microbial mat. Furthermore, the lifestyle of many of the taxa identified from the microbialites was not suited for, or even amenable to, life within a benthic microbial mat. The high diversity of oxygen-dependent metazoans in the non-microbial facies on the Great Bank of Guizhou, and inferences from geochemical records, suggests that the microbialites and benthic communities developed in oxygenated environments, which disproves that the microbes were the source of the oxygenation. Instead, we posit that microbialite successions represent a taphonomic window for exceptional preservation of the biota, similar to a Konzentrat-Lagerstatte, which has allowed for diverse fossil assemblages to be preserved during intervals of poor preservation.
Oceanic Anoxic Event-2 (OAE-2; Cenomanian-Turonian) is characterized by extensive deposition of organic carbon-rich deposits (black shales) in ocean basins worldwide as result of a major perturbation of the global carbon cycle. While the sedimentological, geochemical, and paleontological aspects of deep water expressions of OAE-2 have been intensively studied in the last few decades, much less attention has been given to the coeval shallow water deposits. In this study, we present the results of a detailed facies and petrographic (optical microscope and scanning electron microscopy) and geochemical studies (delta C-13(carb), delta C-13 (org), delta N-15(bulk), TOC, and Rock-Eval pyrolysis) on two key shallow marine sections from the Apennine Carbonate Platform (ACP; Italy). Here a continuous record of shallow water carbonates through the OAE-2 interval is preserved, offering the unique opportunity to document the archive of paleoenvironmental changes in a neritic setting, at a tropical latitude and far from the influence of a large continental block. Two conspicuous intervals are characterized by abundant and closely spaced dark microbial laminites found at correlative stratigraphic horizons in the two stratigraphic sections. These laminites contain elevated concentrations of TOC (up to 1%) relative to microbial capping cycles laminites stratigraphically above and below. The organic matter preserved in these fine-grained laminites is dominated by cyanobacteria remains, which accounted for most of the organic matter produced on the ACP in these intervals. Our study suggests that Tethyan carbonate platforms experienced significant biological changes during OAE-2, alternating, in few kiloyears, between eutrophic phases dominated by microbial communities and mesotrophic/oligotrophic conditions favoring normal carbonate production/sedimentation. The synchronous occurrence of microbialite facies at different locations across the ACP and on other platforms worldwide (e.g., Mexico and Croatia) suggests a causal link between Large Igneous Province volcanism and the environmental conditions necessary to trigger cyanobacterial proliferation on shallow carbonate platforms.
Thermokarst lake landscapes are permafrost regions, which are prone to rapid (on seasonal to decadal time scales) changes, affecting carbon and nitrogen cycles. However, there is a high degree of uncertainty related to the balance between carbon and nitrogen cycling and storage. We collected 12 permafrost soil cores from six drained thermokarst lake basins (DTLBs) along a chronosequence north of Teshekpuk Lake in northern Alaska and analyzed them for carbon and nitrogen contents. For comparison, we included three lacustrine cores from an adjacent thermokarst lake and one soil core from a non thermokarst affected remnant upland. This allowed to calculate the carbon and nitrogen stocks of the three primary landscape units (DTLB, lake, and upland), to reconstruct the landscape history, and to analyze the effect of thermokarst lake formation and drainage on carbon and nitrogen stocks. We show that carbon and nitrogen contents and the carbon-nitrogen ratio are considerably lower in sediments of extant lakes than in the DTLB or upland cores indicating degradation of carbon during thermokarst lake formation. However, we found similar amounts of total carbon and nitrogen stocks due to the higher density of lacustrine sediments caused by the lack of ground ice compared to DTLB sediments. In addition, the radiocarbon-based landscape chronology for the past 7,000years reveals five successive lake stages of partially, spatially overlapping DTLBs in the study region, reflecting the dynamic nature of ice-rich permafrost deposits. With this study, we highlight the importance to include these dynamic landscapes in future permafrost carbon feedback models. Plain Language Summary When permanently frozen soils (permafrost) contain ice-rich sediments, the thawing of this permafrost causes the surface to sink, which may result in lake formation. This process, the thaw of ice-rich permafrost and melting of ground ice leads to characteristic landforms-known as thermokarst. Once such a thaw process is initiated in ice-rich sediments, a thaw lake forms and grows by shoreline erosion, eventually expanding until a drainage pathway is encountered and the lake eventually drains, resulting in a drained thermokarst lake basin. In our study, we show that such a thermokarst-affected landscape north of Teshekpuk Lake in northern Alaska is shaped by repeated thaw lake formation and lake drainage events during the past 7,000years, highlighting the dynamic nature of these landscapes. These landscape-scale processes have a big effect on the carbon and nitrogen stored in permafrost soils. We show that large amounts of carbon (>45kg C/m(2)) and nitrogen (>2.6kg N/m(2)) are stored in unfrozen lake sediments and in frozen soil sediments. The findings are important when considering the potential effect that permafrost thaw has for the global climate through releasing carbon and nitrogen, which was frozen and therefore locked away for millennia, from the active carbon cycle.
Changes in the steepness of river profiles or abrupt vertical steps (i.e. waterfalls) are thought to be indicative of changes in erosion rates, lithology or other factors that affect landscape evolution. These changes are referred to as knickpoints or knickzones and are pervasive in bedrock river systems. Such features are thought to reveal information about landscape evolution and patterns of erosion, and therefore their locations are often reported in the geomorphic literature. It is imperative that studies reporting knickpoints and knickzones use a reproducible method of quantifying their locations, as their number and spatial distribution play an important role in interpreting tectonically active landscapes. In this contribution we introduce a reproducible knickpoint and knickzone extraction algorithm that uses river profiles transformed by integrating drainage area along channel length (the so-called integral or chi method). The profile is then statistically segmented and the differing slopes and step changes in the elevations of these segments are used to identify knickpoints, knickzones and their relative magnitudes. The output locations of identified knickpoints and knickzones compare favourably with human mapping: we test the method on Santa Cruz Island, CA, using previously reported knickzones and also test the method against a new dataset from the Quadrilatero Ferrifero in Brazil. The algorithm allows for the extraction of varying knickpoint morphologies, including stepped, positive slope-break (concave upward) and negative slope-break knickpoints. We identify parameters that most affect the resulting knickpoint and knickzone locations and provide guidance for both usage and outputs of the method to produce reproducible knickpoint datasets.
We analyze trends in compound flooding resulting from high coastal water levels (HCWLs) and peak river discharge over northwestern Europe during 1901-2014. Compound peak discharge associated with 37 stream gauges with at least 70 years of record availability near the North and Baltic Sea coasts is used. Compound flooding is assessed using a newly developed index, compound hazard ratio, that compares the severity of river flooding associated with HCWL with the at-site, T-year (a flood with 1/T chance of being exceeded in any given year) fluvial peak discharge. Our findings suggest a spatially coherent pattern in the dependence between HCWL and river peaks and in compound flood magnitudes and frequency. For higher return levels, we find upward trends in compound hazard ratio frequency at midlatitudes (gauges from 47 degrees N to 60 degrees N) and downward trends along the high latitude (>60 degrees N) regions of northwestern Europe. Plain Language Summary Compound floods in delta areas, that is, the co-occurrence of high coastal water levels (HCWLs) and high river discharge, are a particular challenge for disaster management. Such events are caused by two distinct mechanisms: (1) HCWLs may affect river flows and water levels by backwater effects or by reversing the seaward flow of rivers, particularly in regions with elevation less than 10 m in northwestern Europe. (2) The correlation between HCWL and river flow peaks may also stem from a common meteorological driver. Severe storm periods may be associated with high winds leading to storm surges, and at the same time with high precipitation followed by inland flooding. Understanding the historical trends in compound flooding, owing to changes in relative sea levels, in river flooding and in the dependence between these two drivers, is essential for projecting future changes and disaster management. The risk assessment frameworks are often limited to assessing flood risk from a single driver only. We present a new approach to assess compound flood severity resulting from extreme coastal water level and peak river discharge. We find upward trends in compound flooding for midlatitude regions and downward trends for high latitudes in northwestern Europe.
Compound flooding, such as the co-occurrence of fluvial floods and extreme coastal water levels (CWL), may lead to significant impacts in densely-populated Low Elevation Coastal Zones. They may overstrain disaster management owing to the co-occurrence of inundation from rivers and the sea. Recent studies are limited by analyzing joint dependence between river discharge and either CWL or storm surges, and little is known about return levels of compound flooding, accounting for the covariance between drivers. Here, we assess the compound flood severity and identify hotspots for northwestern Europe during 1970–2014, using a newly developed Compound Hazard Ratio (CHR) that compares the severity of compound flooding associated with extreme CWL with the unconditional T-year fluvial peak discharge. We show that extreme CWL and stronger storms greatly amplify fluvial flood hazards. Our results, based on frequency analyses of observational records during 2013/2014’s winter storm Xaver, reveal that the river discharge of the 50-year compound flood is up to 70% larger, conditioned on the occurrence of extreme CWL, than that of the at-site peak discharge. For this event, nearly half of the stream gauges show increased flood hazards, demonstrating the importance of including the compounding effect of extreme CWL in river flood risk management.
We analyzed the Lomas de Carabajal area in the intermontane Lerma valley of the Cordillera Oriental to assess the level of neotectonic activity in a densely populated region of northwestern Argentina. In this region, Plio-Pleistocene synorogenic conglomerates are deformed, locally associated with high-angle faults, and NNW-SSE oriented en-echelon folds characterized by wavelengths of < 1 km. The deformed Quaternary units follow the same pattern of deformation as observed in the underlying Neogene deposits; growth-strata geometries are observed near faults. This configuration is compatible with local left-lateral transpressional tectonism driven by ENE-WSW buttressing against the NW-oriented border of a Cretaceous extensional basin (Alemania sub-basin). Optically Stimulated Luminescence analysis of sandy-silty layers interbedded within the folded late Pleistocene conglomeratic sequence helps to determine uplift rates of 0.83-0.87 mm/a during the last 30-40 ka. Nearby the Lomas de Carabajal, a WNW-striking, 3-m-high fault scarp disrupts radiocarbon dated, 10-ka-old loessic deposits providing a Holocene mean uplift rate of 0.30 mm/a. Our data unambiguously show that shallow crustal deformation in the intermontane Lerma valley is ongoing; some of this deformation may be associated with seismicity. Our findings support the notion of temporally and spatially disparate deformation processes in the broken foreland of the northwestern Argentinean Andes.
Vermetid reefs and rocky shores are hot spots of biodiversity, often referred to as the subtropical equivalent of coral reefs. The development of the ecosystem depends on the activity of several reef builders, including red crustose coralline algae (CCA) such as Neogoniolithon brassica-florida. Despite its importance, little is known about Neogoniolithon sp. acclimation to rapid changes in light intensity and corresponding photosynthetic activity. To overcome the large spatial variability in the light field (due to location and the porous nature of the rocks) we grew Neogoniolithon sp. on glass slides and characterized its photosynthetic performance in response to various light intensities by following O-2 exchange and fluorescence parameters. This was also performed on rock-inhabiting thalli collected from the east Mediterranean basin. Generally, maximal photosynthetic rate was reached when Neogoniolithon sp. thalli grown under low illumination (such as in protected niches where the light intensity can be as low as 1% of surface illumination) were examined. When exposed to light intensities higher than those experienced during growth, Neogoniolithon sp. activates adaptive/protective mechanisms such as state transition and nonphotochemical fluorescence quenching and increases the dark respiration thereafter. We find that the Fv/Fm parameter (variable/maximal fluorescence) is not suitable to assess photosynthetic performance in Neogoniolithon sp. and propose using instead an alternative parameter recently developed. Our findings help to clarify why Neogoniolithon sp. is usually observed in shaded niches along the reef surfaces.
Encroachment of shrubs into the unique pastoral grassland ecosystems of the Tibetan Plateau has significant impact on ecosystem services, especially forage production. We developed a process-based ecohydrological model to identify the relative importance of the main drivers of shrub encroachment for the alpine meadows within the Qinghai province. Specifically, we explored the effects of summer livestock grazing (intensity and type of livestock) together with the effects of climate warming, including interactions between herbaceous and woody vegetation and feedback loops between soil, water and vegetation. Under current climatic conditions and a traditional herd composition, an increasing grazing intensity above a threshold value of 0.32 +/- 0.10 large stock units (LSU) ha(-1) day(-1) changes the vegetation composition from herbaceous towards a woody and bare soil dominated system. Very high grazing intensity (above 0.8 LSU ha(-1) day(-1)) leads to a complete loss of any vegetation. Under warmer conditions, the vegetation showed a higher resilience against livestock farming. This resilience is enhanced when the herd has a higher browser : grazer ratio. A cooler climate has a shrub encroaching effect, whereas warmer conditions increase the cover of the herbaceous vegetation. This effect was primarily due to season length and an accompanied competitive loss of slower growing shrubs, rather than evaporative water loss leading to less soil water in deeper soil layers for deeper rooting shrubs. If climate warming is driving current shrub encroachment, we conclude it is only indirectly so. It would be manifest by an advancing shrubline and could be regarded as a climatic escape of specific shrub species such as Potentilla fruticosa. Under the recent high intensity of grazing, only herding by more browsing animals can potentially prevent both shrub encroachment and the complete loss of herbaceous vegetation.
An essential, respected, and critical aspect of the modern practice of science and scientific publishing is peer review. The process of peer review facilitates best practices in scientific conduct and communication, ensuring that manuscripts published are as accurate, valuable, and clearly communicated. The over 216 papers published in Tectonics in 2018 benefit from the time, effort, and expertise of our reviewers who have provided thoughtfully considered advice on each manuscript. This role is critical to advancing our understanding of the evolution of the continents and their margins, as these reviews lead to even clearer and higher-quality papers. In 2018, the over 443 papers submitted to Tectonics were the beneficiaries of more than 1,010 reviews provided by 668 members of the tectonics community and related disciplines. To everyone who has volunteered their time and intellect to peer reviewing, thank you for helping Tectonics and all other AGU Publications provide the best science possible.
Variability of the Cold Season Climate in Central Asia. Part II: Hydroclimatic Predictability
(2019)
Central Asia (CA) is subjected to a large variability of precipitation. This study presents a statistical model, relating precipitation anomalies in three subregions of CA in the cold season (November-March) with various predictors in the preceding October. Promising forecast skill is achieved for two subregions covering 1) Uzbekistan, Turkmenistan, Kyrgyzstan, Tajikistan, and southern Kazakhstan and 2) Iran, Afghanistan, and Pakistan. ENSO in October is identified as the major predictor. Eurasian snow cover and the quasi-biennial oscillation further improve the forecast performance. To understand the physical mechanisms, an analysis of teleconnections between these predictors and the wintertime circulation over CA is conducted. The correlation analysis of predictors and large-scale circulation indices suggests a seasonal persistence of tropical circulation modes and a dynamical forcing of the westerly circulation by snow cover variations over Eurasia. An EOF analysis of pressure and humidity patterns allows separating the circulation variability over CA into westerly and tropical modes and confirms that the identified predictors affect the respective circulation characteristics. Based on the previously established weather type classification for CA, the predictors are investigated with regard to their effect on the regional circulation. The results suggest a modification of the Hadley cell due to ENSO variations, with enhanced moisture supply from the Arabian Gulf during El Nino. They further indicate an influence of Eurasian snow cover on the wintertime Arctic Oscillation (AO) and Northern Hemispheric Rossby wave tracks. Positive anomalies favor weather types associated with dry conditions, while negative anomalies promote the formation of a quasi-stationary trough over CA, which typically occurs during positive AO conditions.
Fold and thrust belts are characteristic features of collisional orogen that grow laterally through time by deforming the upper crust in response to stresses caused by convergence. The deformation propagation in the upper crust is accommodated by shortening along major folds and thrusts. The formation of these structures is influenced by the mechanical strength of décollements, basement architecture, presence of preexisting structures and taper of the wedge. These factors control not only the sequence of deformation but also cause differences in the structural style.
The Himalayan fold and thrust belt exhibits significant differences in the structural style from east to west. The external zone of the Himalayan fold and thrust belt, also called the Subhimalaya, has been extensively studied to understand the temporal development and differences in the structural style in Bhutan, Nepal and India; however, the Subhimalaya in Pakistan remains poorly studied. The Kohat and Potwar fold and thrust belts (herein called Kohat and Potwar) represent the Subhimalaya in Pakistan. The Main Boundary Thrust (MBT) marks the northern boundary of both Kohat and Potwar, showing that these belts are genetically linked to foreland-vergent deformation within the Himalayan orogen, despite the pronounced contrast in structural style. This contrast becomes more pronounced toward south, where the active strike-slip Kalabagh Fault Zone links with the Kohat and Potwar range fronts, known as the Surghar Range and the Salt Range, respectively. The Surghar and Salt Ranges developed above the Surghar Thrust (SGT) and Main Frontal Thrust (MFT). In order to understand the structural style and spatiotemporal development of the major structures in Kohat and Potwar, I have used structural modeling and low temperature thermochronolgy methods in this study. The structural modeling is based on construction of balanced cross-sections by integrating surface geology, seismic reflection profiles and well data. In order to constrain the timing and magnitude of exhumation, I used apatite (U-Th-Sm)/He (AHe) and apatite fission track (AFT) dating. The results obtained from both methods are combined to document the Paleozoic to Recent history of Kohat and Potwar.
The results of this research suggest two major events in the deformation history. The first major deformation event is related to Late Paleozoic rifting associated with the development of the Neo-Tethys Ocean. The second major deformation event is related to the Late Miocene to Pliocene development of the Himalayan fold and thrust belt in the Kohat and Potwar. The Late Paleozoic rifting is deciphered by inverse thermal modelling of detrital AFT and AHe ages from the Salt Range. The process of rifting in this area created normal faulting that resulted in the exhumation/erosion of Early to Middle Paleozoic strata, forming a major unconformity between Cambrian and Permian strata that is exposed today in the Salt Range. The normal faults formed in Late Paleozoic time played an important role in localizing the Miocene-Pliocene deformation in this area. The combination of structural reconstructions and thermochronologic data suggest that deformation initiated at 15±2 Ma on the SGT ramp in the southern part of Kohat. The early movement on the SGT accreted the foreland into the Kohat deforming wedge, forming the range front. The development of the MBT at 12±2 Ma formed the northern boundary of Kohat and Potwar. Deformation propagated south of the MBT in the Kohat on double décollements and in the Potwar on a single basal décollement. The double décollement in the Kohat adopted an active roof-thrust deformation style that resulted in the disharmonic structural style in the upper and lower parts of the stratigraphic section. Incremental shortening resulted in the development of duplexes in the subsurface between two décollements and imbrication above the roof thrust. Tectonic thickening caused by duplexes resulted in cooling and exhumation above the roof thrust by removal of a thick sequence of molasse strata. The structural modelling shows that the ramps on which duplexes formed in Kohat continue as tip lines of fault propagation folds in the Potwar. The absence of a double décollement in the Potwar resulted in the preservation of a thick sequence of molasse strata there. The temporal data suggest that deformation propagated in-sequence from ~ 8 to 3 Ma in the northern part of Kohat and Potwar; however, internal deformation in the Kohat was more intense, probably required for maintaining a critical taper after a significant load was removed above the upper décollement. In the southern part of Potwar, a steeper basement slope (β≥3°) and the presence of salt at the base of the stratigraphic section allowed for the complete preservation of the stratigraphic wedge, showcased by very little internal deformation. Activation of the MFT at ~4 Ma allowed the Salt Range to become the range front of the Potwar. The removal of a large amount of molasse strata above the MFT ramp enhanced the role of salt in shaping the structural style of the Salt Range and Kalabagh Fault Zone. Salt accumulation and migration resulted in the formation of normal faults in both areas. Salt migration in the Kalabagh fault zone has triggered out-of-sequence movement on ramps in the Kohat.
The amount of shortening calculated between the MBT and the SGT in Kohat is 75±5 km and between the MBT and the MFT in Potwar is 65±5 km. A comparable amount of shortening is accommodated in the Kohat and Potwar despite their different widths: 70 km Kohat and 150 km Potwar. In summary, this research suggests that deformation switched between different structures during the last ~15 Ma through different modes of fault propagation, resulting in different structural styles and the out-of-sequence development of Kohat and Potwar.
Abstract. The Sea of Marmara, in northwestern Turkey, is a transition zone where the dextral North Anatolian Fault zone (NAFZ) propagates westward from the Anatolian Plate to the Aegean Sea Plate. The area is of interest in the context of seismic hazard of Istanbul, a metropolitan area with about 15 million inhabitants. Geophysical observations indicate that the crust is heterogeneous beneath the Marmara basin, but a detailed characterization of the crustal heterogeneities is still missing. To assess if and how crustal heterogeneities are related to the NAFZ segmentation below the Sea of Marmara, we develop new crustal-scale 3-D density models which integrate geological and seismological data and that are additionally constrained by 3-D gravity modeling. For the latter, we use two different gravity datasets including global satellite data and local marine gravity observation. Considering the two different datasets and the general non-uniqueness in potential field modeling, we suggest three possible “end-member” solutions that are all consistent with the observed gravity field and illustrate the spectrum of possible solutions. These models indicate that the observed gravitational anomalies originate from significant density heterogeneities within the crust. Two layers of sediments, one syn-kinematic and one pre-kinematic with respect to the Sea of Marmara formation are underlain by a heterogeneous crystalline crust. A felsic upper crystalline crust (average density of 2720 kgm⁻³) and an intermediate to mafic lower crystalline crust (average density of 2890 kgm⁻³) appear to be cross-cut by two large, dome-shaped mafic highdensity bodies (density of 2890 to 3150 kgm⁻³) of considerable thickness above a rather uniform lithospheric mantle (3300 kgm⁻³). The spatial correlation between two major bends of the main Marmara fault and the location of the highdensity bodies suggests that the distribution of lithological heterogeneities within the crust controls the rheological behavior along the NAFZ and, consequently, maybe influences fault segmentation and thus the seismic hazard assessment in the region.
Abstract. The Sea of Marmara, in northwestern Turkey, is a transition zone where the dextral North Anatolian Fault zone (NAFZ) propagates westward from the Anatolian Plate to the Aegean Sea Plate. The area is of interest in the context of seismic hazard of Istanbul, a metropolitan area with about 15 million inhabitants. Geophysical observations indicate that the crust is heterogeneous beneath the Marmara basin, but a detailed characterization of the crustal heterogeneities is still missing. To assess if and how crustal heterogeneities are related to the NAFZ segmentation below the Sea of Marmara, we develop new crustal-scale 3-D density models which integrate geological and seismological data and that are additionally constrained by 3-D gravity modeling. For the latter, we use two different gravity datasets including global satellite data and local marine gravity observation. Considering the two different datasets and the general non-uniqueness in potential field modeling, we suggest three possible “end-member” solutions that are all consistent with the observed gravity field and illustrate the spectrum of possible solutions. These models indicate that the observed gravitational anomalies originate from significant density heterogeneities within the crust. Two layers of sediments, one syn-kinematic and one pre-kinematic with respect to the Sea of Marmara formation are underlain by a heterogeneous crystalline crust. A felsic upper crystalline crust (average density of 2720 kgm⁻³) and an intermediate to mafic lower crystalline crust (average density of 2890 kgm⁻³) appear to be cross-cut by two large, dome-shaped mafic highdensity bodies (density of 2890 to 3150 kgm⁻³) of considerable thickness above a rather uniform lithospheric mantle (3300 kgm⁻³). The spatial correlation between two major bends of the main Marmara fault and the location of the highdensity bodies suggests that the distribution of lithological heterogeneities within the crust controls the rheological behavior along the NAFZ and, consequently, maybe influences fault segmentation and thus the seismic hazard assessment in the region.
Climate science is highly interdisciplinary by nature, so understanding interactions between Earth processes inherently warrants the use of analytical software that can operate across the disciplines of Earth science. Toward this end, we present the Climate Data Toolbox for MATLAB, which contains more than 100 functions that span the major climate-related disciplines of Earth science. The toolbox enables streamlined, entirely scriptable workflows that are intuitive to write and easy to share. Included are functions to evaluate uncertainty, perform matrix operations, calculate climate indices, and generate common data displays. Documentation is presented pedagogically, with thorough explanations of how each function works and tutorials showing how the toolbox can be used to replicate results of published studies. As a well-tested, well-documented platform for interdisciplinary collaborations, the Climate Data Toolbox for MATLAB aims to reduce time spent writing low-level code, let researchers focus on physics rather than coding and encourage more efficacious code sharing. Plain Language Summary This article describes a collection of computer code that has recently been released to help scientists analyze many types of Earth science data. The code in this toolbox makes it easy to investigate things like global warming, El Nino, or other major climate-related processes such as how winds affect ocean circulation. Although the toolbox was designed to be used by expert climate scientists, its instruction manual is well written, and beginners may be able to learn a great deal about coding and Earth science, simply by following along with the provided examples. The toolbox is intended to help scientists save time, help them ensure their analysis is accurate, and make it easy for other scientists to repeat the results of previous studies.
The in-phase response collected by portable loop-loop electromagnetic induction (EMI) sensors operating at low and moderate induction numbers (<= 1) is typically used for sensing the magnetic permeability (or susceptibility) of the subsurface. This is due to the fact that the in-phase response contains a small induction fraction and a preponderant induced magnetization fraction. The magnetization fraction follows the magneto-static equations similarly to the magnetic method but with an active magnetic source. The use of an active source offers the possibility to collect data with several loop-loop configurations, which illuminate the subsurface with different sensitivity patterns. Such multiconfiguration soundings thereby allows the imaging of subsurface magnetic permeability/susceptibility variations through an inversion procedure. This method is not affected by the remnant magnetization and theoretically overcomes the classical depth ambiguity generally encountered with passive geomagnetic data. To invert multiconfiguration in-phase data sets, we propose a novel methodology based on a full-grid 3-D multichannel deconvolution (MCD) procedure. This method allows us to invert large data sets (e.g. consisting of more than a hundred thousand of data points) for a dense voxel-based 3-D model of magnetic susceptibility subject to smoothness constraints. In this study, we first present and discuss synthetic examples of our imaging procedure, which aim at simulating realistic conditions. Finally, we demonstrate the applicability of our method to field data collected across an archaeological site in Auvergne (France) to image the foundations of a Gallo-Roman villa built with basalt rock material. Our synthetic and field data examples demonstrate the potential of the proposed inversion procedure offering new and complementary ways to interpret data sets collected with modern EMI instruments.
Nature-based solutions (NBS) have recently received attention due to their potential ability to sustainably reduce hydro-meteorological risks, providing co-benefits for both ecosystems and affected people. Therefore, pioneering research has dedicated efforts to optimize the design of NBS, to evaluate their wider co-benefits and to understand promoting and/or hampering governance conditions for the uptake of NBS. In this article, we aim to complement this research by conducting a comprehensive literature review of factors shaping people’s perceptions of NBS as a means to reduce hydro-meteorological risks. Based on 102 studies, we identified six topics shaping the current discussion in this field of research: (1) valuation of the co-benefits (including those related to ecosystems and society); (2) evaluation of risk reduction efficacy; (3) stakeholder participation; (4) socio-economic and location-specific conditions; (5) environmental attitude, and (6) uncertainty. Our analysis reveals that concerned empirical insights are diverse and even contradictory, they vary in the depth of the insights generated and are often not comparable for a lack of a sound theoretical-methodological grounding. We, therefore, propose a conceptual model outlining avenues for future research by indicating potential inter-linkages between constructs underlying perceptions of NBS to hydro-meteorological risks.
Arboreal epiphytes (plants residing in forest canopies) are present across all major climate zones and play important roles in forest biogeochemistry. The substantial water storage capacity per unit area of the epiphyte "bucket" is a key attribute underlying their capability to influence forest hydrological processes and their related mass and energy flows. It is commonly assumed that the epiphyte bucket remains saturated, or near-saturated, most of the time; thus, epiphytes (particularly vascular epiphytes) can store little precipitation, limiting their impact on the forest canopy water budget. We present evidence that contradicts this common assumption from (i) an examination of past research; (ii) new datasets on vascular epiphyte and epi-soil water relations at a tropical montane cloud forest (Monteverde, Costa Rica); and (iii) a global evaluation of non-vascular epiphyte saturation state using a process-based vegetation model, LiBry. All analyses found that the external and internal water storage capacity of epiphyte communities is highly dynamic and frequently available to intercept precipitation. Globally, non-vascular epiphytes spend <20% of their time near saturation and regionally, including the humid tropics, model results found that non-vascular epiphytes spend similar to 1/3 of their time in the dry state (0-10% of water storage capacity). Even data from Costa Rican cloud forest sites found the epiphyte community was saturated only 1/3 of the time and that internal leaf water storage was temporally dynamic enough to aid in precipitation interception. Analysis of the epi-soils associated with epiphytes further revealed the extent to which the epiphyte bucket emptied-as even the canopy soils were often <50% saturated (29-53% of all days observed). Results clearly show that the epiphyte bucket is more dynamic than currently assumed, meriting further research on epiphyte roles in precipitation interception, redistribution to the surface and chemical composition of "net" precipitation waters reaching the surface.
Water stable isotope signatures can provide valuable insights into the catchment internal runoff processes. However, the ability of the water isotope data to constrain the internal apportionments of runoff components in hydrological models for glacierized basins is not well understood. This study developed an approach to simultaneously model the water stable isotopic compositions and runoff processes in a glacierized basin in Central Asia. The fractionation and mixing processes of water stable isotopes in and from the various water sources were integrated into a glacio-hydrological model. The model parameters were calibrated on discharge, snow cover and glacier mass balance data, and additionally isotopic composition of streamflow. We investigated the value of water isotopic compositions for the calibration of model parameters, in comparison to calibration methods without using such measurements. Results indicate that: (1) The proposed isotope-hydrological integrated modeling approach was able to reproduce the isotopic composition of streamflow, and improved the model performance in the evaluation period; (2) Involving water isotopic composition for model calibration reduced the model parameter uncertainty, and helped to reduce the uncertainty in the quantification of runoff components; (3) The isotope-hydrological integrated modeling approach quantified the contributions of runoff components comparably to a three-component tracer-based end-member mixing analysis method for summer peak flows, and required less water tracer data. Our findings demonstrate the value of water isotopic compositions to improve the quantification of runoff components using hydrological models in glacierized basins.
The computation of such synthetic GFs is computationally and operationally demanding. As a consequence, the onthe-fly recalculation of synthetic GFs in each iteration of an optimisation is time-consuming and impractical. Therefore, the pre-calculation and efficient storage of synthetic GFs on a dense grid of source to receiver combinations enables the efficient lookup and utilisation of GFs in time-critical scenarios. We present a Python-based framework and toolkit - Pyrocko-GF - that enables the pre-calculation of synthetic GF stores, which are independent of their numerical calculation method and GF transfer function. The framework aids in the creation of such GF stores by interfacing a suite of established numerical forward modelling codes in seismology (computational back ends). So far, interfaces to back ends for layered Earth model cases have been provided; however, the architecture of Pyrocko-GF is designed to cover back ends for other geometries (e.g. full 3-D heterogeneous media) and other physical quantities (e.g. gravity, pressure, tilt). Therefore, Pyrocko-GF defines an extensible GF storage format suitable for a wide range of GF types, especially handling elasticity and wave propagation problems. The framework assists with visualisations, quality control, and the exchange of GF stores, which is supported through an online platform that provides many pre-calculated GF stores for local, regional, and global studies. The Pyrocko-GF toolkit comes with a well-documented application programming interface (API) for the Python programming language to efficiently facilitate forward modelling of geophysical processes, e.g. synthetic waveforms or static displacements for a wide range of source models.
The occurrence of deep low-frequency (DLF) microearthquakes beneath volcanoes is commonly attributed to mass transport in the volcanic plumbing system and used to infer feeding channels from and into magma reservoirs. The key question is how magmas migrate from depth to the shallow crust and whether magma reservoirs are currently being recharged. For the first time since the improvement of the local seismic networks in the East Eifel region (Rhineland-Palatinate, Germany), we detect and locate recurrent DLF earthquakes in the lower crust and upper mantle beneath the Laacher See Volcano (LSV), using a joint data set of permanent sensors and a temporary deployment. So far, eight DLF earthquake sequences were observed in four distinct clusters between 10 and 40 km depth. These clusters of weak events (M-L< 2) align along an approximately 80. southeast dipping line south of the LSV. Moment tensor solutions of these events have large shear components, and the irregular dispersion and long coda of body waves indicate interaction processes between shear cracks and fluids. We find a rotation of P-axes orientation for shallow tectonic earthquakes compared to DLF events, indicating that the stress field in the depth interval of DLF events might favour a vertical migration of magma or magmatic fluids. The caldera of the LSV was formed by the last major eruption of the East Eifel Volcanic Field only 12.9 kyr ago, fed by a shallow magma chamber at 5-8 km depth and erupting a total magma volume of 6.7 km(3). The observed DLF earthquake activity and continuous volcanic gas emissions around the LSV indicate an active magmatic system, possibly connected with an upper mantle melt zone.
The thesis comprises three experimental studies, which were carried out to unravel the short- as well as the long-term mechanical properties of shale rocks. Short-term mechanical properties such as compressive strength and Young’s modulus were taken from recorded stress-strain curves of constant strain rate tests. Long-term mechanical properties are represented by the time– dependent creep behavior of shales. This was obtained from constant stress experiments, where the test duration ranged from a couple minutes up to two weeks. A profound knowledge of the mechanical behavior of shales is crucial to reliably estimate the potential of a shale reservoir for an economical and sustainable extraction of hydrocarbons (HC). In addition, healing of clay-rich forming cap rocks involving creep and compaction is important for underground storage of carbon dioxide and nuclear waste.
Chapter 1 introduces general aspects of the research topic at hand and highlights the motivation for conducting this study. At present, a shift from energy recovered from conventional resources e.g., coal towards energy provided by renewable resources such as wind or water is a big challenge. Gas recovered from unconventional reservoirs (shale plays) is considered a potential bridge technology.
In Chapter 2, short-term mechanical properties of two European mature shale rocks are presented, which were determined from constant strain rate experiments performed at ambient and in situ deformation conditions (confining pressure, pc ≤ 100 MPa, temperature, T ≤ 125 °C, representing pc, T - conditions at < 4 km depth) using a Paterson– type gas deformation apparatus. The investigated shales were mainly from drill core material of Posidonia (Germany) shale and weathered material of Bowland (United Kingdom) shale. The results are compared with mechanical properties of North American shales. Triaxial compression tests performed perpendicular to bedding revealed semibrittle deformation behavior of Posidonia shale with pronounced inelastic deformation. This is in contrast to Bowland shale samples that deformed brittle and displayed predominantly elastic deformation. The static Young’s modulus, E, and triaxial compressive strength, σTCS, determined from recorded stress-strain curves strongly depended on the applied confining pressure and sample composition, whereas the influence of temperature and strain rate on E and σTCS was minor. Shales with larger amounts of weak minerals (clay, mica, total organic carbon) yielded decreasing E and σTCS. This may be related to a shift from deformation supported by a load-bearing framework of hard phases (e.g., quartz) towards deformation of interconnected weak minerals, particularly for higher fractions of about 25 – 30 vol% weak phases. Comparing mechanical properties determined at reservoir conditions with mechanical data applying effective medium theories revealed that E and σTCS of Posidonia and Bowland shale are close to the lower (Reuss) bound. Brittleness B is often quoted as a measure indicating the response of a shale formation to stimulation and economic production. The brittleness, B, of Posidonia and Bowland shale, estimated from E, is in good agreement with the experimental results. This correlation may be useful to predict B from sonic logs, from which the (dynamic) Young’s modulus can be retrieved.
Chapter 3 presents a study of the long-term creep properties of an immature Posidonia shale. Constant stress experiments (σ = const.) were performed at elevated confining pressures (pc = 50 – 200 MPa) and temperatures (T = 50 – 200 °C) to simulate reservoir pc, T - conditions. The Posidonia shale samples were acquired from a quarry in South Germany. At stresses below ≈ 84 % compressive strength of Posidonia shale, at high temperature and low confining pressure, samples showed pronounced transient (primary) creep with high deformation rates in the semibrittle regime. Sample deformation was mainly accommodated by creep of weak sample constituents and pore space reduction. An empirical power law relation between strain and time, which also accounts for the influence of pc, T and σ on creep strain was formulated to describe the primary creep phase. Extrapolation of the results to a creep period of several years, which is the typical time interval for a large production decline, suggest that fracture closure is unlikely at low stresses. At high stresses as expected for example at the contact between the fracture surfaces and proppants added during stimulation measures, subcritical crack growth may lead to secondary and tertiary creep. An empirical power law is suggested to describe secondary creep of shale rocks as a function of stress, pressure and temperature. The predicted closure rates agree with typical production decline curves recorded during the extraction of hydrocarbons. At the investigated conditions, the creep behavior of Posidonia shale was found to correlate with brittleness, calculated from sample composition.
In Chapter 4 the creep properties of mature Posidonia and Bowland shales are presented. The observed long-term creep behavior is compared to the short-term behavior determined in Chapter 2. Creep experiments were performed at simulated reservoir conditions of pc = 50 – 115 MPa and T = 75 – 150 °C. Similar to the mechanical response of immature Posidonia shale samples investigated in Chapter 3, creep strain rates of mature Bowland and Posidonia shales were enhanced with increasing stress and temperature and decreasing confining pressures. Depending on applied deformation conditions, samples displayed either only a primary (decelerating) or in addition also a secondary (quasi-steady state) and subsequently a tertiary (accelerating) creep phase before failure. At the same deformation conditions, creep strain of Posidonia shale, which is rich in weak constituents, is tremendously higher than of quartz-rich Bowland shale. Typically, primary creep strain is again mostly accommodated by deformation of weak minerals and local pore space reduction. At the onset of tertiary creep most of the deformation was accommodated by micro crack growth. A power law was used to characterize the primary creep phase of Posidonia and Bowland shale. Primary creep strain of shale rocks is inversely correlated to triaxial compressive strength and brittleness, as described in Chapter 2.
Chapter 5 provides a synthesis of the experimental findings and summarizes the major results of the studies presented in Chapters 2 – 4 and potential applications in the Exploration & Production industry.
Chapter 6 gives a brief outlook on potential future experimental research that would help to further improve our understanding of processes leading to fracture closure involving proppant embedment in unconventional shale gas reservoirs. Such insights may allow to improve stimulation techniques aimed at maintaining economical extraction of hydrocarbons over several years.
Issue Despite their rather similar climatic conditions, eastern Eurasia and northern North America are largely covered by different plant functional types (deciduous or evergreen boreal forest) composed of larch or pine, spruce and fir, respectively. I propose that these deciduous and evergreen boreal forests represent alternative quasi-stable states, triggered by their different northern tree refugia that reflect the different environmental conditions experienced during the Last Glacial. Evidence This view is supported by palaeoecological and environmental evidence. Once established, Asian larch forests are likely to have stabilized through a complex vegetation-fire-permafrost soil-climate feedback system. Conclusion With respect to future forest developments, this implies that Asian larch forests are likely to be governed by long-term trajectories and are therefore largely resistant to natural climate variability on time-scales shorter than millennia. The effects of regional human impact and anthropogenic global warming might, however, cause certain stability thresholds to be crossed, meaning that irreversible transitions occur and resulting in marked consequences for ecosystem services on these human-relevant time-scales.
Permafrost thaw subjects previously frozen soil organic carbon (SOC) to microbial degradation to the greenhouse gases carbon dioxide (CO2) and methane (CH4). Emission of these gases constitutes a positive feedback to climate warming. Among numerous uncertainties in estimating the strength of this permafrost carbon feedback (PCF), two are: (i) how mineralization of permafrost SOC thawed in saturated anaerobic conditions responds to changes in temperature and (ii) how microbial communities and temperature sensitivities change over time since thaw. To address these uncertainties, we utilized a thermokarst-lake sediment core as a natural chronosequence where SOC thawed and incubated in situ under saturated anaerobic conditions for up to 400 years following permafrost thaw. Initial microbial communities were characterized, and sediments were anaerobically incubated in the lab at four temperatures (0 °C, 3 °C, 10 °C, and 25 °C) bracketing those observed in the lake's talik. Net CH4 production in freshly-thawed sediments near the downward-expanding thaw boundary at the base of the talik were most sensitive to warming at the lower incubation temperatures (0 °C to 3 °C), while the overlying sediments which had been thawed for centuries had initial low abundant methanogenic communities (< 0.02%) and did not experience statistically significant increases in net CH4 production potentials until higher incubation temperatures (10 °C to 25 °C). We propose these observed differences in temperature sensitivities are due to differences in SOM quality and functional microbial community composition that evolve over time; however further research is necessary to better constrain the roles of these factors in determining temperature controls on anaerobic C mineralization.
In this study, we detect high percentile rainfall events in the eastern central Andes, based on Tropical Rainfall Measuring Mission (TRMM) with a spatial resolution of 0.25 × 0.25°, a temporal resolution of 3 h, and for the duration from 2001 to 2018. We identify three areas with high mean accumulated rainfall and analyze their atmospheric behaviour and rainfall characteristics with specific focus on extreme events. Extreme events are defined by events above the 95th percentile of their daily mean accumulated rainfall. Austral summer (DJF) is the period of the year presenting the most frequent extreme events over these three regions. Daily statistics show that the spatial maxima, as well as their associated extreme events, are produced during the night. For the considered period, ERA-Interim reanalysis data, provided by the European Centre for Medium-Range Weather Forecasts (ECMWF) with 0.75° x0.75° spatial and 6-hourly temporal resolutions, were used for the analysis of the meso- and synoptic-scale atmospheric patterns. Night- and day-time differences indicate a nocturnal overload of northerly and northeasterly low-level humidity flows arriving from tropical South America. Under these conditions, cooling descending air from the mountains may find unstable air at the surface, giving place to the development of strong local convection. Another possible mechanism is presented here: a forced ascent of the low-level flow due to the mountains, disrupting the atmospheric stratification and generating vertical displacement of air trajectories. A Principal Component Analysis (PCA) in T-mode is applied to day- and night-time data during the maximum and extreme events. The results show strong correlation areas over each subregion under study during night-time, whereas during day-time no defined patterns are found. This confirms the observed nocturnal behavior of rainfall within these three hotspots.
The trace gases CO2 and CH4 pertain to the most relevant greenhouse gases and are important exchange fluxes of the global carbon (C) cycle. Their atmospheric quantity increased significantly as a result of the intensification of anthropogenic activities, such as especially land-use and land-use change, since the mid of the 18th century. To mitigate global climate change and ensure food security, land-use systems need to be developed, which favor reduced trace gas emissions and a sustainable soil carbon management. This requires the accurate and precise quantification of the influence of land-use and land-use change on CO2 and CH4 emissions. A common method to determine the trace gas dynamics and C sink or source function of a particular ecosystem is the closed chamber method. This method is often used assuming that accuracy and precision are high enough to determine differences in C gas emissions for e.g., treatment comparisons or different ecosystem components.
However, the broad range of different chamber designs, related operational procedures and data-processing strategies which are described in the scientific literature contribute to the overall uncertainty of closed chamber-based emission estimates. Hence, the outcomes of meta-analyses are limited, since these methodical differences hamper the comparability between studies. Thus, a standardization of closed chamber data acquisition and processing is much-needed.
Within this thesis, a set of case studies were performed to: (I) develop standardized routines for an unbiased data acquisition and processing, with the aim of providing traceable, reproducible and comparable closed chamber based C emission estimates; (II) validate those routines by comparing C emissions derived using closed chambers with independent C emission estimates; and (III) reveal processes driving the spatio-temporal dynamics of C emissions by developing (data processing based) flux separation approaches.
The case studies showed: (I) the importance to test chamber designs under field conditions for an appropriate sealing integrity and to ensure an unbiased flux measurement. Compared to the sealing integrity, the use of a pressure vent and fan was of minor importance, affecting mainly measurement precision; (II) that the developed standardized data processing routines proved to be a powerful and flexible tool to estimate C gas emissions and that this tool can be successfully applied on a broad range of flux data sets from very different ecosystem; (III) that automatic chamber measurements display temporal dynamics of CO2 and CH4 fluxes very well and most importantly, that they accurately detect small-scale spatial differences in the development of soil C when validated against repeated soil inventories; and (IV) that a simple algorithm to separate CH4 fluxes into ebullition and diffusion improves the identification of environmental drivers, which allows for an accurate gap-filling of measured CH4 fluxes.
Overall, the proposed standardized data acquisition and processing routines strongly improved the detection accuracy and precision of source/sink patterns of gaseous C emissions. Hence, future studies, which consider the recommended improvements, will deliver valuable new data and insights to broaden our understanding of spatio-temporal C gas dynamics, their particular environmental drivers and underlying processes.
Large-magnitude fluid-injection induced seismic events are a potential risk for geothermal energy developments worldwide. One potential risk mitigation measure is the application of cyclic injection schemes. After validation at small (laboratory) and meso (mine) scale, the concept has now been applied for the first time at field scale at the Pohang Enhanced Geothermal System (EGS) site in Korea. From 7 August until 14 August 2017 a total of 1756 m(3) of surface water was injected into Pohang well PX-1 at flow rates between 1 and 10 l s(-1), with a maximum wellhead pressure (WHP) of 22.8 MPa, according to a site-specific cyclic soft stimulation schedule and traffic light system. A total of 52 induced microearthquakes were detected in real-time during and shortly after the injection, the largest of M-w 1.9. After that event a total of 1771 m(3) of water was produced back from the well over roughly 1 month, during which time no larger-magnitude seismic event was observed. The hydraulic data set exhibits pressure-dependent injectivity increase with fracture opening between 15 and 17 MPa WHP, but no significant permanent transmissivity increase was observed. The maximum magnitude of the induced seismicity during the stimulation period was below the target threshold of M-w 2.0 and additional knowledge about the stimulated reservoir was gained. Additionally, the technical feasibility of cyclic injection at field scale was evaluated. The major factors that limited the maximum earthquake magnitude are believed to be: limiting the injected net fluid volume, flowback after the occurrence of the largest induced seismic event, using a cyclic injection scheme, the application of a traffic light system, and including a priori information from previous investigations and operations in the treatment design.
Although the Paris Agreement arguably made some progress, interest in supplementary approaches to climate change co-operation persist. This article examines the conditions under which a climate club might emerge and grow. Using agent-based simulations, it shows that even with less than a handful of major actors as initial members, a club can eventually reduce global emissions effectively. To succeed, a club must be initiated by the ‘right’ constellation of enthusiastic actors, offer sufficiently large incentives for reluctant countries and be reasonably unconstrained by conflicts between members over issues beyond climate change. A climate club is particularly likely to persist and grow if initiated by the United States and the European Union. The combination of club-good benefits and conditional commitments can produce broad participation under many conditions.
One commonly proposed method to limit flood risk is land-use or zoning policies which regulates construction in high-risk areas, in order to reduce economic exposure and its vulnerability to flood events. Although such zoning regulations can be effective in limiting trends in flood risk, they also have adverse impacts on society, for instance by limiting local development of areas near the water. In order to judge whether proposed land-use or zoning policies are a net benefit to society, they should be accepted or rejected based on a societal cost-benefit analysis (CBA). However, conducting a CBA of zoning regulation is complex and comprehensive guidelines of how to do such an analysis are lacking. We offer guidelines for good practice. In order to assess the costs and benefits of zoning as a climate change adaption strategy, they should be assessed at a societal level in order to account for public good features of flood risk reduction strategies, and because costs in one area can be benefits in another region. We propose a multistep process: first, determine the spatial extent of the zoning policy and how interconnected the zoned area is to other locations; second, conduct a CBA using monetary costs and benefits estimated from an integrated hydro-economic model to investigate if total benefits exceed total costs; third, conduct a sensitivity analysis regarding the main assumptions; fourth, conduct a multicriteria analysis (MCA) of the normative outcomes of a zoning policy. A desirable policy is preferred in both the CBA and MCA. This article is categorized under: Engineering Water > Planning Water Human Water > Value of Water Science of Water > Water Extremes Human Water > Methods
Flood risk will increase in many areas around the world due to climate change and increase in economic exposure. This implies that adequate flood insurance schemes are needed to adapt to increasing flood risk and to minimise welfare losses for households in flood-prone areas. Flood insurance markets may need reform to offer sufficient and affordable financial protection and incentives for risk reduction. Here, we present the results of a study that aims to evaluate the ability of flood insurance arrangements in Europe to cope with trends in flood risk, using criteria that encompass common elements of the policy debate on flood insurance reform. We show that the average risk-based flood insurance premium could double between 2015 and 2055 in the absence of more risk reduction by households exposed to flooding. We show that part of the expected future increase in flood risk could be limited by flood insurance mechanisms that better incentivise risk reduction by policyholders, which lowers vulnerability. The affordability of flood insurance can be improved by introducing the key features of public-private partnerships (PPPs), which include public reinsurance, limited premium cross-subsidisation between low- and high-risk households, and incentives for policyholder-level risk reduction. These findings were evaluated in a comprehensive sensitivity analysis and support ongoing reforms in Europe and abroad that move towards risk-based premiums and link insurance with risk reduction, strengthen purchase requirements, and engage in multi-stakeholder partnerships.
Flood disasters severely impact human subjective well-being (SWB). Nevertheless, few studies have examined the influence of flood events on individual well-being and how such impacts may be limited by flood protection measures. This study estimates the long term impacts on individual subjective well-being of flood experiences, individual subjective flood risk perceptions, and household flood preparedness decisions. These effects are monetised and placed in context through a comparison with impacts of other adverse events on well-being. We collected data from households in flood-prone areas in France. The results indicate that experiencing a flood has a large negative impact on subjective well-being that is incompletely attenuated over time. Moreover, individuals do not need to be directly affected by floods to suffer SWB losses since subjective well-being is lower for those who expect their flood risk to increase or who have seen a neighbour being flooded. Floodplain inhabitants who prepared for flooding by elevating their home have a higher subjective well-being. A monetisation of the aforementioned well-being impacts shows that a flood requires Euro150,000 in immediate compensation to attenuate SWB losses. The decomposition of the monetised impacts of flood experience into tangible losses and intangible effects on SWB shows that intangible effects are about twice as large as the tangible direct monetary flood losses. Investments in flood protection infrastructure may be under funded if the intangible SWB benefits of flood protection are not taken into account.