Refine
Has Fulltext
- no (216) (remove)
Year of publication
- 2019 (216) (remove)
Document Type
- Article (216) (remove)
Language
- English (216)
Is part of the Bibliography
- yes (216)
Keywords
- climate (4)
- Holocene (3)
- paleoclimate (3)
- remote sensing (3)
- Asia (2)
- Central Andes (2)
- Climate (2)
- Deformation (2)
- Earthquake (2)
- Flood risk (2)
- InSAR (2)
- Magnetotellurics (2)
- NW Argentina (2)
- Paleoclimate (2)
- Paleogene (2)
- Paleolimnology (2)
- Pamir (2)
- Permafrost (2)
- Permian (2)
- Pleistocene (2)
- Pollen (2)
- Pyrenees (2)
- Siberia (2)
- Site effects (2)
- Time-series analysis (2)
- Validation (2)
- flood risk (2)
- geochronology (2)
- geomagnetic observatories (2)
- hydrology (2)
- induced seismicity (2)
- magnetostratigraphy (2)
- melt inclusions (2)
- modelling (2)
- photogrammetry (2)
- sedimentology (2)
- silicon (2)
- tectonic geomorphology (2)
- wave-particle interactions (2)
- westerlies (2)
- 2-D numerical model (1)
- 2D numerical experiments (1)
- AMOC (1)
- Abrupt events (1)
- Absorption feature parameters (1)
- Accuracy Asseessment (1)
- Adaptation process (1)
- Agent-based models (1)
- Amplification (1)
- Anatolia (1)
- Ancient Gneiss Complex (1)
- Andean Plateau (1)
- Andean retroarc (1)
- Andes (1)
- Annual laminae (1)
- Apennine Carbonate Platform (ACP) (1)
- Archaeolithoporella (1)
- Argentinean broken foreland (1)
- Array seismology (1)
- Asian monsoon (1)
- Assimilation-fractional crystallization (1)
- Asymmetric rifting (1)
- BNDVI (1)
- Back-arc basin (1)
- Baiu (1)
- Baldeggersee (1)
- Basin analysis (1)
- Bayesian hierarchical model (1)
- Bayesian modeling (1)
- Bayesianism (1)
- Be-10 (1)
- Beacon fire (1)
- Beringia (1)
- Binalud Mountains (1)
- Biomisation (1)
- Blueschist (1)
- Blueschist metamorphism (1)
- Brazil (1)
- Britholite (1)
- British Columbia (1)
- Burned areas (1)
- Bykovsky Peninsula (1)
- C isotopes (1)
- C-14 bomb peak (1)
- CU (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)
- Chemeron Formation (1)
- Chile (1)
- City population (1)
- Clay mineralogy (1)
- Climate change (1)
- Climate change adaptation (1)
- Climate prediction (1)
- Climate warming (1)
- Common garden experiments (1)
- Complex networks (1)
- Connectivity (1)
- Conservation planning (1)
- Controlled source electromagnetics (CSEM) (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)
- Digital Elevation Model (1)
- Dryland (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)
- Electrical conductivity (1)
- Electromagnetic theory (1)
- Emission scenarios (1)
- Environmental (1)
- Environmental impact (1)
- Environmental magnetism (1)
- Eocene (1)
- Erodium cicutarium (1)
- Error (1)
- Error propagation (1)
- Europe (1)
- European air quality (1)
- Evaporitic varves (1)
- Event synchronization (1)
- Expression (1)
- Extraction (1)
- Extreme habitats (1)
- Extreme rainfall (1)
- Extremophiles (1)
- FAIR (1)
- Feature selection (1)
- Federal Water Act (1)
- Field analogue (1)
- Flood (1)
- Fluid processes (1)
- Food security (1)
- Forecasting Framework (1)
- Foreland basin fragmentation (1)
- GNSS Reflectometry (1)
- Generalized additive models (1)
- Geochemistry (1)
- Geodynamic modelling (1)
- Geomechanical Model (1)
- Geomorphology (1)
- Geothermobarometry (1)
- Germany (1)
- Glacial lake outbursts (1)
- Glacial refugia (1)
- Glacierized basins (1)
- Gravity modelling (1)
- Ground-penetrating radar (1)
- Groundwater quality (1)
- Gypsum quantification (1)
- Gyttja (1)
- Haida Gwaii (1)
- Heat flow (1)
- Herschel Island Qikiqtaruk (1)
- High pressure (1)
- Hillslope thermokarst (1)
- Himalayan hydroclimate (1)
- Hindu Kush-Karakoram (1)
- Hochwasser (1)
- Hochwasserrisikomanagementrichtlinie (1)
- Hutchinson niche (1)
- Hydraulic fracturing (1)
- Hydro-mechanical coupling (1)
- Hydrometric networks (1)
- Hypersaline lake (1)
- Hyperspectral (1)
- Hypolimnetic oxygen (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)
- International unions (1)
- Inverse theory (1)
- Irrawaddy River (1)
- Irrigation (1)
- Isotope-hydrological integrated modeling (1)
- Italy (1)
- Ivrea Zone (1)
- Jaguaribe Basin (1)
- Kameng River section (1)
- Kenya Rift (1)
- KiK-net (1)
- Kp index (1)
- Lacustrine sediments (1)
- Lake sediments (1)
- Lamina thickness (1)
- Landsat (1)
- Landsat-8 (1)
- Landscape Response (1)
- Landslide (1)
- Larix larch (1)
- Last Interglacial (1)
- Last interglacial-glacial transition (1)
- Lawsonite (1)
- Lena Delta (1)
- LiDAR (1)
- Lithospheric structure (1)
- Local knowledge (1)
- Local-to-regional scale (1)
- Long-Term (1)
- Lu/Hf dating of garnet (1)
- MO (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)
- Methane (1)
- Methane leakage (1)
- Methanogenic archaea (1)
- Microcrack (1)
- Microseismic monitoring (1)
- Middle Holocene (1)
- Middle Strand (1)
- Millennial-scale variability (1)
- Minimum mortality temperature (1)
- Model development (1)
- Model-data comparison (1)
- Mountain building (1)
- Multidimensionality (1)
- N isotopes (1)
- NW Turkey (1)
- Namibia (1)
- Namibian passive margin (1)
- Neotectonics (1)
- Neotethys (1)
- Nordeste (1)
- North America (1)
- North Anatolian Fault (1)
- North Atlantic Oscillation (1)
- Northeast Iran (1)
- Northern Andes (1)
- Northern Asia (1)
- Northwestern Anatolia (1)
- Numerical model (1)
- Numerical modelling (1)
- OAE-2 (1)
- OLCI (1)
- OSL (1)
- OSL and C-14 geochronology (1)
- OpenForecast (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)
- Particle flow code (1)
- Patterns (1)
- Peat (1)
- Per capita growth rate (1)
- Perchlorate (1)
- Permafrost degradation (1)
- Pharmaceuticals (1)
- Phosphates (1)
- Phytolith dissolution (1)
- Phytolith morphotypes (1)
- Pliocene (1)
- Pohang (Korea) (1)
- Precession (1)
- Precipitation (1)
- Predictive models (1)
- Preparedness (1)
- Principal component analysis (1)
- Proportional population growth rate (1)
- Provenance (1)
- Provenance studies (1)
- Pu239+240 (1)
- Public-private partnerships (1)
- Pull-Apart (1)
- Puna (1)
- Quantification of runoff components (1)
- Quaternary (1)
- REE (1)
- Radiogenic isotopes (1)
- Radioisotope disequilibria dating (1)
- Raman (1)
- Recurrence plot (1)
- Red River (1)
- Reefs (1)
- Regional cooperation (1)
- Remote sensing (1)
- Reservoir Networks (1)
- Review (1)
- Risk perception (1)
- Risk reduction (1)
- River Incision Model (1)
- River-groundwater-interaction (1)
- Rotifers (1)
- Russia (1)
- Rutile mineral chemistry (1)
- SIO₂ (1)
- Salt pan (1)
- Sandstone detrital modes (1)
- Sea of Marmara (1)
- Seasonal forecasting (1)
- Sediment Transport (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)
- Shear Zone (1)
- Si cycling (1)
- Si extraction (1)
- Sigmoid model (1)
- Simulated Mars-like conditions (1)
- Site proxy (1)
- Siwalik Group (1)
- Slip Distribution (1)
- Slope (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)
- 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)
- Subduction (1)
- Subduction Zone (1)
- Subfreezing temperatures (1)
- Subjective well-being (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)
- Tarim Basin (1)
- Tarim basin (1)
- Tectonic Evolution (1)
- Tectonic reconstruction (1)
- Temperature sensitivity (1)
- Tephrostratigraphy (1)
- TerraSAR-X (1)
- TerraceM (1)
- The Geysers (1)
- Theoretical seismology (1)
- Thermal evolution (1)
- Tibetan plateau (1)
- Tien Shan Mountains (1)
- Topography (1)
- Transfer functions (1)
- Transpression (1)
- Tree-ring software (1)
- Triassic (1)
- Triphylite (1)
- Tropical storms (1)
- U-Pb geochronology (1)
- UAV (1)
- UAVs (1)
- UNFCCC (1)
- Uncertainties (1)
- Uncertainty (1)
- Uncertainty Processor (1)
- Urban water cycle (1)
- Varve chronologies (1)
- Varved sediments (1)
- Varves (1)
- Vegetation driver (1)
- Viscous blocking (1)
- Volcan de Colima (1)
- Volcano monitoring (1)
- Volcano seismology (1)
- WRF-Chem (1)
- Wasserhaushaltsgesetz (1)
- Water Availability (1)
- Water stable isotope (1)
- Wave propagation (1)
- Waveform inversion (1)
- West Africa (1)
- Willingness to pay (1)
- Wuchiapingian (1)
- Xining Basin (1)
- Yarlung Tsangpo (1)
- Zooplankton (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)
- alps (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)
- assessment (1)
- atmosphere (1)
- attitudes (1)
- attraction (1)
- automatic lineament extraction (1)
- bank infiltration (1)
- best practices (1)
- biomarker (1)
- biostratigraphy (1)
- boreal forests (1)
- boundary element modeling (1)
- braided alluvial rivers (1)
- caldera collapse (1)
- carbon (1)
- carbon burial (1)
- carbon capture and storage (CCS) (1)
- carbon capture and utilization (CCU) (1)
- carbon sequestration (1)
- catchment (1)
- channel migration (1)
- circulation (1)
- climate change adaptation (1)
- climate change co-operation (1)
- climate clubs (1)
- clinopyroxenite (1)
- club goods (1)
- clustering (1)
- co-benefits (1)
- coastal erosion (1)
- coastal geomorphology (1)
- collision (1)
- complexes (1)
- compositional gap (1)
- compound flood (1)
- conditional commitments (1)
- connectivity (1)
- continuous modelling (1)
- convective available potential energy (1)
- correlation (1)
- correlation skill (1)
- corridors (1)
- cosmic-ray neutron sensing (1)
- cosmogenic nuclides (1)
- cost of electricity (COE) (1)
- costs of infrastructure (1)
- crustal fault (1)
- crystalline (1)
- data (1)
- data science (1)
- decomposition (1)
- deep convection (1)
- deforestation (1)
- delta drift (1)
- delta-c-13 (1)
- dendroclimatology (1)
- denudation (1)
- dependence (1)
- depositional environment (1)
- derived flood frequency (1)
- detrital carbonate (1)
- detrital zircons (1)
- dew-point temperature (1)
- dietary patterns (1)
- disaster risk reduction (1)
- disinfection (1)
- dislocation creep (1)
- dissolution properties (1)
- drainage networks (1)
- dynamical seasonal prediction (1)
- earthquake (1)
- earthquake nucleation (1)
- earthquake sequences (1)
- eastern south–central Andes (1)
- eclogite (1)
- economics (1)
- ecosystem services (1)
- efficient method (1)
- electrical resistivity (1)
- electricity generation (1)
- electron lifetime (1)
- electron transport (1)
- enhanced geothermal systems (EGS) (1)
- ensemble modeling (1)
- equatorial ionosphere (1)
- erosion (1)
- erosion rates (1)
- evaporites (1)
- exsolution textures (1)
- extinction (1)
- extreme rainfall (1)
- fenite (1)
- finite element modeling (1)
- flood hazard (1)
- flood loss model transfer (1)
- flowering time (1)
- fluid inclusions (1)
- fluid regime (1)
- fluids (1)
- food systems (1)
- forecasting (1)
- forest structure (1)
- fractional snow cover (1)
- functional trait correlations (1)
- gender (1)
- geochemistry (1)
- geomagnetic indices (1)
- geomagnetic observatory data (1)
- geomagnetic storm drivers (1)
- geomagnetism (1)
- geomorphic markers (1)
- geomorphology (1)
- glass (1)
- granite (1)
- gravity observations (1)
- gravity-driven slope deformation (1)
- green infrastructure (1)
- greenhouse gas emissions (1)
- groundwater (1)
- heavy rainfall (1)
- herd composition (1)
- high mountain environments (1)
- high-density reservoir network (1)
- hiss waves (1)
- historical geomagnetic storms (1)
- human activity (1)
- hummocky (1)
- hydro-meterological hazards (1)
- hydrogeology (1)
- hydrothermal alteration (1)
- image enhancement (1)
- incision (1)
- individual-based model (1)
- inner magnetosphere (1)
- intangible impacts (1)
- interferometric SAR (InSAR) (1)
- intermontane basins (1)
- intrabasinal faulting (1)
- intraplate deformation (1)
- invasive species (1)
- inverse analysis (1)
- irrigation (1)
- kettle hole (1)
- land cover change (1)
- landscape evolution (1)
- landscape evolution modeling (1)
- landslide (1)
- landslides (1)
- late Holocene (1)
- lava dome (1)
- lidar (1)
- life-history strategy (1)
- lithosphere dynamics (1)
- livestock (1)
- load stress (1)
- local adaptation (1)
- local food (1)
- lumped parameter model (1)
- magnetospheric convection (1)
- marine terraces (1)
- mean transit time (1)
- megathrust earthquake (1)
- melting (1)
- metasomatism (1)
- methane (1)
- methanol (1)
- microbialite (1)
- microfossil (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)
- nanocomposite material (1)
- nearshore zone (1)
- neotectonics (1)
- neutron simulations (1)
- nontidal ocean loading (1)
- normal faults (1)
- northern high latitudes (1)
- northwestern Europe (1)
- object-based damage modeling (1)
- ocean color remote sensing (1)
- open (1)
- operational service (1)
- opinion (1)
- optical data (1)
- optical remote sensing (1)
- orogenic peridotite (1)
- orogeny (1)
- oxygen (1)
- oxygen fractionation model (1)
- paleoceanography (1)
- paleomagnetism (1)
- paleosol (1)
- parallel evolution (1)
- part I (1)
- participation (1)
- peatlands (1)
- peri-urban agriculture (1)
- permafrost (1)
- permafrost disturbances (1)
- permafrost ecosystems (1)
- phenotypic plasticity (1)
- physical experiments (1)
- phytolith (1)
- pitch angle diffusion coefficient (1)
- plate boundary fault (1)
- plate coupling force (1)
- point cloud (1)
- point clouds (1)
- pollution (1)
- pre-monsoon (1)
- precipitation (1)
- preferences (1)
- pressures (1)
- prior derivation (1)
- probabilistic approaches (1)
- process simulation (1)
- prograding lobes (1)
- project database (1)
- proto-Paratethys (1)
- pseudomorphism (1)
- radiation belts (1)
- rain attenuation (1)
- rain effect (1)
- rainy-season (1)
- rangeland management (1)
- reactive oxygen species (1)
- record (1)
- refuge (1)
- regression (1)
- repulsion (1)
- reservoir volume (1)
- resilience (1)
- rheological weakening (1)
- ring current electrons (1)
- ring faulting (1)
- risk assessment (1)
- risk mitigation (1)
- risk modeling (1)
- river erosion (1)
- river incision (1)
- river networks (1)
- river plume (1)
- river restoration (1)
- rock glacier (1)
- runoff (1)
- rutile-type (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)
- seismotectonic segmentation (1)
- semi-arid (1)
- sentinel-1 (1)
- shrub encroachment (1)
- shrubline (1)
- siderite (1)
- silvopastoral systems (1)
- simulation model (1)
- slow slip (1)
- small baseline subset (SBAS) (1)
- snow characteristics (1)
- snow grain size (1)
- soil erosion (1)
- spatial heterogeneity (1)
- spatial sampling (1)
- spatial scales (1)
- specific surface area (1)
- spectral adjustment (1)
- speleothems (1)
- spin transition (1)
- stable isotopes (1)
- standards (1)
- strain localization (1)
- strain transient (1)
- strainmeter data (1)
- structural mapping (1)
- structure from motion (1)
- subsea permafrost (1)
- summer (1)
- supervised machine learning (1)
- surface processes (1)
- surface reflectance (1)
- surface runoff (1)
- surveys (1)
- suspended particulate matter (1)
- suspended sediment (1)
- synthetic flood events (1)
- systems (1)
- techno-economic model (1)
- tectonism (1)
- temporal disaggregation (1)
- temporal evolution (1)
- the Ogooue Delta (1)
- tillage (1)
- time series (1)
- time series analysis (1)
- tin (1)
- topographic analysis (1)
- torsion (1)
- transform faults (1)
- transient response (1)
- tundra-taiga ecotone (1)
- turbidity (1)
- two phase aggregates (1)
- uncertainty (1)
- uncertainty analysis (1)
- underground coal gasification (UCG) (1)
- urbanization (1)
- vadose zone (1)
- varved sediments (1)
- vegetation states (1)
- vegetation trajectories (1)
- vegetation-climate-fire-soil feedbacks (1)
- volcanic deformation (1)
- water (1)
- water height-area-volume curve (1)
- wave number (1)
- weathering (1)
- welfare (1)
- wind speed (1)
Institute
- Institut für Geowissenschaften (216) (remove)
The thermal structure of subduction zones exerts a major influence on deep-seated mechanical and chemical processes controlling arc magmatism, seismicity, and global element cycles. Accretionary complexes exposed inland may comprise tectonic blocks with contrasting pressure-temperature (P-T) histories, making it possible to investigate the dynamics and thermal evolution of former subduction interfaces. With this aim, we present new Lu-Hf geochronological results for mafic rocks of the Halilbagi Complex (Anatolia) that evolved along different thermal gradients. Samples include a lawsonite-epidote blueschist, a lawsonite-epidote eclogite, and an epidote eclogite (all with counter-clockwise P-T paths), a prograde lawsonite blueschist with a "hairpin"-type P-T path, and a garnet amphibolite from the overlying sub-ophiolitic metamorphic sole. Equilibrium phase diagrams suggest that the garnet amphibolite formed at similar to 0.6-0.7 GPa and 800-850 degrees C, whereas the prograde lawsonite blueschist records burial from 2.1 GPa and 420 degrees C to 2.6 GPa and 520 degrees C. Well-defined Lu-Hf isochrons were obtained for the epidote eclogite (92.38 +/- 0.22 Ma) and the lawsonite-epidote blueschist (90.19 +/- 0.54 Ma), suggesting rapid garnet growth. The lawsonite-epidote eclogite (87.30 +/- 0.39 Ma) and the prograde lawsonite blueschist (ca. 86 Ma) are younger, whereas the garnet amphibolite (104.5 +/- 3.5 Ma) is older. Our data reveal a consistent trend of progressively decreasing geothermal gradient from granulite-facies conditions at similar to 104 Ma to the epidote-eclogite facies around 92 Ma, and the lawsonite blueschist-facies between 90 Ma and 86 Ma. Three Lu-Hf garnet dates (between 92 Ma and 87 Ma) weighted toward the growth of post-peak rims (as indicated by Lu distribution in garnet) suggest that the HP/LT rocks were exhumed continuously and not episodically. We infer that HP/LT metamorphic rocks within the Halilbagi Complex were subjected to continuous return flow, with "warm" rocks being exhumed during the tectonic burial of "cold" ones. Our results, combined with regional geological constraints, allow us to speculate that subduction started at a transform fault near a mid-oceanic spreading centre. Following its formation, this ancient subduction interface evolved thermally over more than 15 Myr, most likely as a result of heat dissipation rather than crustal underplating. (C) 2018, China University of Geosciences (Beijing) and Peking University. Production and hosting by Elsevier B.V.
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.
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.
The Raman spectra of pure N-2, CO2, and CH4 were analyzed over the range 10 to 500 bars and from -160 degrees C to 200 degrees C (N-2), 22 degrees C to 350 degrees C (CO2), and -100 degrees C to 450 degrees C (CH4). At constant temperature, Raman peak position, including the more intense CO2 peak (nu+), decreases (shifts to lower wave number) with increasing pressure for all three gases over the entire pressure and temperature (PT) range studied. At constant pressure, the peak position for CO2 and CH4 increases (shifts to higher wave number) with increasing temperature over the entire PT range studied. In contrast, N-2 first shows an increase in peak position with increasing temperature at constant pressure, followed by a decrease in peak position with increasing temperature. The inflection temperature at which the trend reverses for N-2 is located between 0 degrees C and 50 degrees C at pressures above similar to 50 bars and is pressure dependent. Below similar to 50 bars, the inflection temperature was observed as low as -120 degrees C. The shifts in Raman peak positions with PT are related to relative density changes, which reflect changes in intermolecular attraction and repulsion. A conceptual model relating the Raman spectral properties of N-2, CO2, and CH4 to relative density (volume) changes and attractive and repulsive forces is presented here. Additionally, reduced temperature-dependent densimeters and barometers are presented for each pure component over the respective PT ranges. The Raman spectral behavior of the pure gases as a function of temperature and pressure is assessed to provide a framework for understanding the behavior of each component in multicomponent N-2-CO2-CH4 gas systems in a future study.
We address the question of whether all large-magnitude earthquakes produce an erosion peak in the subaerial components of fluvial catchments. We evaluate the sediment flux response to the Maule earthquake in the Chilean Andes (Mw 8.8) using daily suspended sediment records from 31 river gauges. The catchments cover drainage areas of 350 to around 10,000 km(2), including a wide range of topographic slopes and vegetation cover of the Andean western flank. We compare the 3- to 8-year postseismic record of sediment flux to each of the following preseismic periods: (1) all preseismic data, (2) a 3-year period prior to the seismic event, and (3) the driest preseismic periods, as drought conditions prevailed in the postseismic period. Following the earthquake, no increases in suspended sediment flux were observed for moderate to high percentiles of the streamflow distribution (mean, median, and >= 75th percentile). However, more than half of the examined stations showed increased sediment flux during baseflow. By using a Random Forest approach, we evaluate the contributions of seismic intensities, peak ground accelerations, co-seismic landslides, hydroclimatic conditions, topography, lithology, and land cover to explain the observed changes in suspended sediment concentration and fluxes. We find that the best predictors are hillslope gradient, low-vegetation cover, and changes in streamflow discharge. This finding suggests a combined first-order control of topography, land cover, and hydrology on the catchment-wide erosion response. We infer a reduced sediment connectivity due to the postseismic drought, which increased the residence time of sediment detached and remobilized following the Maule earthquake.
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.
Alluvial and transport-limited bedrock rivers constitute the majority of fluvial systems on Earth. Their long profiles hold clues to their present state and past evolution. We currently possess first-principles-based governing equations for flow, sediment transport, and channel morphodynamics in these systems, which we lack for detachment-limited bedrock rivers. Here we formally couple these equations for transport-limited gravel-bed river long-profile evolution. The result is a new predictive relationship whose functional form and parameters are grounded in theory and defined through experimental data. From this, we produce a power-law analytical solution and a finite-difference numerical solution to long-profile evolution. Steady-state channel concavity and steepness are diagnostic of external drivers: concavity decreases with increasing uplift rate, and steepness increases with an increasing sediment-to-water supply ratio. Constraining free parameters explains common observations of river form: to match observed channel concavities, gravel-sized sediments must weather and fine - typically rapidly - and valleys typically should widen gradually. To match the empirical square-root width-discharge scaling in equilibrium-width gravel-bed rivers, downstream fining must occur. The ability to assign a cause to such observations is the direct result of a deductive approach to developing equations for landscape evolution.
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.
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.
Debate persists concerning the timing and geodynamics of intercontinental collision, style of syncollisional deformation, and development of topography and fold-and-thrust belts along the >1,700-km-long Izmir-Ankara-Erzincan suture zone (IAESZ) in Turkey. Resolving this debate is a necessary precursor to evaluating the integrity of convergent margin models and kinematic, topographic, and biogeographic reconstructions of the Mediterranean domain. Geodynamic models argue either for a synchronous or diachronous collision during either the Late Cretaceous and/or Eocene, followed by Eocene slab breakoff and postcollisional magmatism. We investigate the collision chronology in western Anatolia as recorded in the sedimentary archives of the 90-km-long Saricakaya Basin perched at shallow structural levels along the IAESZ. Based on new zircon U-Pb geochronology and depositional environment and sedimentary provenance results, we demonstrate that the Saricakaya Basin is an Eocene sedimentary basin with sediment sourced from both the IAESZ and Sogut Thrust fault to the south and north, respectively, and formed primarily by flexural loading from north-south shortening along the syncollisional Sogut Thrust. Our results refine the timing of collision between the Anatolides and Pontide terranes in western Anatolia to Maastrichtian-Middle Paleocene and Early Eocene crustal shortening and basin formation. Furthermore, we demonstrate contemporaneous collision, deformation, and magmatism across the IAESZ, supporting synchronous collision models. We show that regional postcollisional magmatism can be explained by renewed underthrusting instead of slab breakoff. This new IAESZ chronology provides additional constraints for kinematic, geodynamic, and biogeographic reconstructions of the Mediterranean domain.
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.
The ‘bomb-pulse’ method is a chronological approach to further constrain the age of speleothems that grew between 1950 CE – present. Establishing dependable chronological constraints is crucial for modern calibration studies of speleothems to instrumental climate records, which provides the basis for paleoclimate interpretations. However, a large unknown is how 14C is transferred from the atmosphere to any individual speleothem owing to the site-specific residence times of organic matter above cave systems. Here, we employ the bomb-pulse method to build chronologies from 14C measurements in combination with a new unsaturated zone C model which considers C decomposition as a continuum, to better understand unsaturated zone 14C dynamics. The bomb-pulse curves of eight speleothems from southern Australia in three contrasting climatic regions; the semi-arid Wellington Caves site, the mediterranean Golgotha Cave site and the montane Yarrangobilly Caves site, are investigated. Overall, the modelled 14C bomb-pulse curves produce excellent fits with measured 14C speleothem data (r2 = 0.82–0.99). The C modelling reveals that unsaturated zone C is predominately young at the semi-arid site, with a weighted-mean residence time of 32 years and that tree root respiration is likely an important source of vadose CO2. At the montane site, ∼39% of C is young (<1 years), but the weighted-mean C ages are older (145–220 years). The mediterranean site has very little contribution from young C (<12%: 0–1 years), with weighted-mean ages between 157 and 245 years, likely due to greater adsorption of organic matter in the upper vadose zone during matrix flow, and remobilisation of C from young syngenetic karst. New end members for low speleothem Dead Carbon Proportion (DCP) are identified (2.19% and 1.65%, respectively) for Australian montane and semi-arid zone speleothems, where oversupply of modern CO2 in the vadose zone leads to lower DCP. It is also demonstrated that DCP can be quite variable over small time scales, that processes may be difficult to untangle and a constant DCP assumption is likely invalid. DCP variability over time is mainly controlled by the changes vadose zone CO2, where vegetation regeneration, wild-fires and karst hydrology play an important role.
An Overview of Using Weather Radar for Climatological Studies: Successes, Challenges, and Potential
(2019)
Weather radars have been widely used to detect and quantify precipitation and nowcast severe weather for more than 50 years. Operational weather radars generate huge three-dimensional datasets that can accumulate to terabytes per day. So it is essential to review what can be done with existing vast amounts of data, and how we should manage the present datasets for the future climatologists. All weather radars provide the reflectivity factor, and this is the main parameter to be archived. Saving reflectivity as volumetric data in the original spherical coordinates allows for studies of the three-dimensional structure of precipitation, which can be applied to understand a number of processes, for example, analyzing hail or thunderstorm modes. Doppler velocity and polarimetric moments also have numerous applications for climate studies, for example, quality improvement of reflectivity and rain rate retrievals, and for interrogating microphysical and dynamical processes. However, observational data alone are not useful if they are not accompanied by sufficient metadata. Since the lifetime of a radar ranges between 10 and 20 years, instruments are typically replaced or upgraded during climatologically relevant time periods. As a result, present metadata often do not apply to past data. This paper outlines the work of the Radar Task Team set by the Atmospheric Observation Panel for Climate (AOPC) and summarizes results from a recent survey on the existence and availability of long time series. We also provide recommendations for archiving current and future data and examples of climatological studies in which radar data have already been used.
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.
With the onset of the global food crisis, the discussion about the use and misuse of agricultural market interventions regained academic attention. As a result of economies of scale, centralized policy implementation at the regional level has the potential to reduce the budgetary costs of policies. Borrowing from the literature on international unions and international policy coordination, we develop a conceptual framework to analyze when regional policy implementation makes sense. This is the case whenever spill-overs from centralization are large and policy preferences, driven by country-specific characteristics, are homogeneous. Subsequently, we examine the advantageousness of centralized policy implementation for the West African region regarding the most common food security policies. We show that centralization of trade policies and emergency food reserves is beneficial, while buffer stocks, safety net policies, and producer support policies should be implemented at the national level.
The processes that control long term landscape evolution in continental interiors and, in particular, along passive margins such as in southern Africa, are still the subject of much debate (e.g. Braun, 2018). Although today the Namibian margin is characterized by an arid climate, it has experienced climatic fluctuations during the Cenozoic and, yet, to date no study has documented the potential role of climate on its erosion history. In western Namibia, the Brandberg Massif, an erosional remnant or inselberg, provides a good opportunity to document the Cenozoic denudation history of the margin using the relationship between rock cooling or exhumation ages and their elevation. Here we provide new apatite (UThSm)/He dates on the Brandberg Inselberg that range from 151 +/- 12 to 30 +/- 2 Ma. Combined with existing apatite fission track data, they yield new constraints on the denudation history of the margin. These data document two main cooling phases since continental break-up 130 Myr ago, a rapid one (similar to 10 degrees C/Myr) following break-up and a slower one (similar to 12 degrees C/Myr) between 65 and 35 Ma. We interpret them respectively to be related to escarpment erosion following rifting and continental break-up and as a phase of enhanced denudation during the Early Eocene Climatic Optimum. We propose that during the Early Eocene Climatic Optimum chemical weathering was important and contributed significantly to the denudation of the Namibian margin and the formation of a pediplain around the Brandberg and enhanced valley incision within the massif. Additionally, aridification of the region since 35 Ma has resulted in negligible denudation rates since that time. (C) 2019 Elsevier B.V. All rights reserved.
The surface deformation associated with the 2010 M-w 8.8 Maule earthquake in Chile was recorded in great detail before, during and after the event. The high data quality of the continuous GPS (cGPS) observations has facilitated a number of studies that model the postseismic deformation signal with a combination of relocking, afterslip and viscoelastic relaxation using linear rheology for the upper mantle. Here, we investigate the impact of using linear Maxwell or power-law rheology with a 2D geomechanical-numerical model to better understand the relative importance of the different processes that control the postseismic deformation signal. Our model results reveal that, in particular, the modeled cumulative vertical postseismic deformation pattern in the near field (< 300 km from the trench) is very sensitive to the location of maximum afterslip and choice of rheology. In the model with power-law rheology, the afterslip maximum is located at 20-35 km rather than > 50 km depth as suggested in previous studies. The explanation for this difference is that in the model with power-law rheology the relaxation of coseismically imposed differential stresses occurs mainly in the lower crust. However, even though the model with power-law rheology probably has more potential to explain the vertical postseismic signal in the near field, the uncertainty of the applied temperature field is substantial, and this needs further investigations and improvements.
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.
Aquifer thermal energy storage (ATES) as a complement to fluctuating renewable energy systems is a reliable technology to guarantee continuous energy supply for heating and air conditioning. We investigated a high-temperature (HT) mono-well system (c. 100 degrees C), where the well screens are separated vertically within the aquifer, as an alternative to conventional doublet ATES systems for an underground storage in northern Oman. We analysed the impact of thermal inference between injection and extraction well screens on the heat recovery factor (HRF) in order to define the optimal screento-screen distance for best possible systems efficiency. Two controlling interference parameters were considered: the vertical screen-to-screen distance and aquifer heterogeneities. The sensitivity study shows that with decreasing screen-to-screen distances, thermal interference increases storage performance. A turning point is reached if the screen distance is too close, causing either water breakthrough or negative thermal interference between the screens. Our simulations show that a combined heat plume with spherical geometry results in the highest heat recovery factors due to the lowest surface area to volume ratios. Thick aquifers for mono-well HT-ATES are thus not mandatory Our study shows that a HT-ATES mono-well system is a feasible storage design with high heat recovery factors for continuous cooling or heating purposes.
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.
Identifying the provenance signature and geodynamic setting on which sedimentary basins at convergent margins grow is challenging since they result from coupled erosional and tectonic processes, which shape the evolution of source areas and the stress regime. The Early Cretaceous evolution of the northern Andes of Colombia is characterized by extensional tectonics and the subsequent formation of a marginal basin. The Abejorral Formation and coeval volcano-sedimentary rocks are exposed along the western flank and axis of the Central Cordillera. They comprise an Early Cretaceous transgressive sequence initially accumulated in fluvial deltaic environments, which switched towards a deep-marine setting, and are interpreted as the infilling record of a marginal back-arc basin. Available provenance data suggest that Permian-Triassic metamorphic and less abundant Jurassic magmatic rocks forming the basement of the Central Cordillera sourced the Abejorral Formation. New detailed volcanic and metamorphic lithics analyses, conventional and varietal study of heavy minerals, detrital rutile mineral chemistry, allowed us to document changes in the source areas defined by the progressive appearance of both higher-grade and more distal low-grade metamorphic sources, which switched from pelitic to dominantly mafic in composition. Crystallochemical indexes of clay minerals of fine-grained rocks of the Abejorral Formation suggest that samples located close to the Romeral Fault System show characteristics of low-medium P-T low-grade metamorphism, whereas rocks located farther to the northeast preserve primary diagenetic features, which suggest a high heat-flow accumulation setting. We interpret that the Abejorral Formation records the progressive unroofing of the Central Cordillera basement that was being rapidly exhumed, as well as the incorporation of distal subduction-related metamorphic complexes to the west in response either to the widening of extensional front or the reactivation of fault structures on the oceanward margin of the basin. Although the deformational record of the Abejorral Formation would have resulted from over-imposed episodes, our new geochronological constraints suggest that this sedimentary sequence must have been deformed before the Paleocene due to the presence of arc-related intrusive non-deformed magmatic rocks with a crystallization age of ca. 60 Ma.
The paper presents a historical long-distance communication system based on beacon fires in one of the most dynamic and rugged mountain ranges of the world, the Hindu Kush-Karakoram-Himalayas. It was deployed as an early warning system for glacial lake outburst floods, which caused devastating impacts on settlement zones and infrastructure until the middle of the twentieth century. The study revealed that the beacon fire systems were operated in distinct valleys spread over the entire Hindu Kush-Karakoram Region. The remarkable fact is the establishment of fire posts in highly difficult accessible mountain environments with communication distances of several hundred kilometres for individual beacon lines. The warning system was a cooperative natural hazard management, which was operated even across distinct ethnic groups. Distant societies, formerly perceived as isolated villages by physical barriers of the high mountain relief, were in historical times connected not only by challenging trade routes but also by a fast working optical communication system. The findings are discussed in the context of a future sustainable natural hazard management.
Main group pallasite meteorites are samples of a single early magmatic planetesimal, dominated by metal and olivine but containing accessory chromite, sulfide, phosphide, phosphates, and rare phosphoran olivine. They represent mixtures of core and mantle materials, but the environment of formation is poorly understood, with a quiescent core-mantle boundary, violent core-mantle mixture, or surface mixture all recently suggested. Here, we review main group pallasite data sets and petrologic characteristics, and present new observations on the low-MnO pallasite Brahin that contains abundant fragmental olivine, but also rounded and angular olivine and potential evidence of sulfide-phosphide liquid immiscibility. A reassessment of the literature shows that low-MnO and high-FeO subgroups preferentially host rounded olivine and low-temperature P2O5-rich phases such as the Mg-phosphate farringtonite and phosphoran olivine. These phases form after metal and silicate reservoirs back-react during decreasing temperature after initial separation, resulting in oxidation of phosphorus and chromium. Farringtonite and phosphoran olivine have not been found in the common subgroup PMG, which are mechanical mixtures of olivine, chromite with moderate Al2O3 contents, primitive solid metal, and evolved liquid metal. Lower concentrations of Mn in olivine of the low-MnO PMG subgroup, and high concentrations of Mn in low-Al2O3 chromites, trace the development and escape of sulfide-rich melt in pallasites and the partially chalcophile behavior for Mn in this environment. Pallasites with rounded olivine indicate that the core-mantle boundary of their planetesimal may not be a simple interface but rather a volume in which interactions between metal, silicate, and other components occur.
We investigate the thermal and structural evolution of asymmetric rifted margin using numerical modeling and geological observations derived from the Western Pyrenees. Our numerical model provides a self-consistent physical evolution of the top basement heat flow during asymmetric rifting. The model shows a pronounced thermal asymmetry that is caused by migration of the rift center toward the upper plate. The same process creates a diachronism for the record of maximum heat flow and maximum temperatures (T-max) in basal rift sequences. The Mauleon-Arzacq basin (W-Pyrenees) corresponds to a former mid-Cretaceous asymmetric hyperextended rift basin. New vitrinite reflectance data in addition to existing data sets from this basin reveal an asymmetry in the distribution of peak heat (T-max) with respect to the rift shoulders, where highest values are located at the former upper- to lower-plate transition. This data set from the Arzacq-Mauleon field study confirms for the first time the thermal asymmetry predicted by numerical models. Numerical modeling results also suggest that complexities in synrift thermal architecture could arise when hanging-wall-derived extensional allochthons and related T-max become part of the lower plate and are transported away from the upper- to lower-plate transition. This study emphasizes the limitations of the common approach to integrate punctual thermal data from pre-rift to synrift sedimentary sequences in order to describe the rift-related thermal evolution and paleothermal gradients at the scale of a rift basin or a rifted margin.
The application of electrical resistivity tomography to peatlands supports conventional coring by providing data on the current condition of peatlands, including data on stratigraphy, peat properties and thickness of organic deposits. Data on the current condition of drained peatlands are particularly required to improve estimates of carbon storage as well as losses and emissions from agriculturally used peatlands. However, most of the studies focusing on electrical resistivity tomography surveys have been conducted on natural peatlands with higher groundwater levels. Peatlands drained for agriculture have not often been studied using geophysical techniques. Drained sites are characterized by low groundwater levels and high groundwater fluctuations during the year, which lead to varying levels of water saturation. To validate better electrical resistivity tomography surveys of drained peatlands, the aim of this laboratory study is to investigate the influence of varying water saturation levels on electrical conductivity (reciprocal of resistivity) for a variety of peat and gyttja types, as well as for different degrees of peat decomposition. Results show that different levels of water saturation strongly influence bulk electrical conductivity. Distinct differences in this relationship exist between peat and gyttja substrates and between different degrees of peat decomposition. Peat shows an exponential relationship for all degrees of decomposition, whereas gyttja, in particular organic-rich gyttja, is characterized by a rather unimodal relationship. The slopes for the relationship between electrical conductivity and water content are steeper at high degrees of decomposition than for peat of low degrees of decomposition. These results have direct implications for field electrical resistivity tomography surveys. In drained peatlands that are strongly susceptible to drying, electrical resistivity tomography surveys have a high potential to monitor the actual field water content. In addition, at comparable water saturations, high or low degrees of decomposition can be inferred from electrical conductivity.
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.
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.
Rapid local adaptation frequently occurs during the spread of invading species. It remains unclear, however, how consistent, and therefore potentially predictable, such patterns of local adaptation are. One approach to this question is to measure patterns of local differentiation in functional traits and plasticity levels in invasive species in multiple regions. Finding consistent patterns of local differentiation in replicate regions suggests that these patterns are adaptive. Further, this outcome indicates that the invading species likely responds predictably to selection along environmental gradients, even though standing genetic variation is likely to have been reduced during introduction. We studied local differentiation in the invasive annual plant Erodium cicutarium in two invaded regions, California and Chile. We collected seeds from across strong gradients in precipitation and temperature in Mediterranean-climate parts of the two regions (10 populations per region). We grew seeds from maternal families from these populations through two generations and exposed the second generation to contrasting levels of water and nutrient availability. We measured growth, flowering time and leaf functional traits across these treatments to obtain trait means and plasticity measures. We found strong differentiation among populations in all traits. Plants from drier environments flowered earlier, were less plastic in flowering time and reached greater size in all treatments. Correlations among traits within regions suggested a coordinated evolutionary response along environmental gradients associated with growing season length. There was little divergence in traits and trait intercorrelations between regions, but strongly parallel divergence in traits within regions. Similar, statistically consistent patterns of local trait differentiation across two regions suggest that local adaptation to environmental gradients has aided the spread of this invasive species, and that the formation of ecotypes in newly invaded environments has been relatively consistent and predictable.
Ground and space-based geomagnetic data were used in the investigation of the longitudinal, seasonal and lunar phase dependence of the equatorial counter electrojet (CEJ) occurrence in the Peruvian, Brazilian, African, Indian and Philippine sectors during geomagnetically quiet days from the solar cycle 24 (2008 to 2018). We found that CEJ events occur more frequently during the morning (MCEJ) than in the afternoon (ACEJ). The highest MCEJ and ACEJ occurrence rates were observed for the Brazilian sector. Distinct seasonal dependence was found for each longitudinal sector under investigation. The lunar phase dependence was determined for the first time for the Philippine sector (longitude 125 degrees E), and it was shown to be less pronounced than in the Peruvian, Brazilian and African sectors. We demonstrate that differences in CEJ rates derived from ground-based and satellite data can arise from the longitudinal separation between low-latitude and equatorial stations that are used to determine the signal and its consequent time delay in their sunrise/sunset times at ionospheric heights.
Lineament mapping, which is an important part of any structural geological investigation, is made more efficient and easier by the availability of optical as well as radar remote sensing data, such as Landsat and Sentinel with medium and high spatial resolutions. However, the results from these multi-resolution data vary due to their difference in spatial resolution and sensitivity to soil occupation. The accuracy and quality of extracted lineaments depend strongly on the spatial resolution of the imagery. Therefore, the aim of this study was to compare the optical Landsat-8, Sentinel-2A, and radar Sentinel-1A satellite data for automatic lineament extraction. The framework of automatic approach includes defining the optimal parameters for automatic lineament extraction with a combination of edge detection and line-linking algorithms and determining suitable bands from optical data suited for lineament mapping in the study area. For the result validation, the extracted lineaments are compared against the manually obtained lineaments through the application of directional filtering and edge enhancement as well as to the lineaments digitized from the existing geological maps of the study area. In addition, a digital elevation model (DEM) has been utilized for an accuracy assessment followed by the field verification. The obtained results show that the best correlation between automatically extracted lineaments, manual interpretation, and the preexisting lineament map is achieved from the radar Sentinel-1A images. The tests indicate that the radar data used in this study, with 5872 and 5865 lineaments extracted from VH and VV polarizations respectively, is more efficient for structural lineament mapping than the Landsat-8 and Sentinel-2A optical imagery, from which 2338 and 4745 lineaments were extracted respectively.
Ongoing development of dynamical atmosphere-ocean general circulation models keep expectations high regarding seasonal predictions of Indian monsoon rainfall. This study compares past and present skill of four currently operating forecasting systems, CFSv2 from NCEP, ENSEMBLES, System 4 and the newest SEAS5 from ECMWF, by analysing correlations of respective hindcasts with observed all-India summer rainfall. For the common time period 1982-2005, only ENSEMBLES and CFSv2 give significantly skilful forecasts. It is shown that skill is highly dependent on the chosen time period. Especially the intense El Nino of 1997 seems to degrade the predictions, most notably for SEAS4 and SEAS5 which seem to be linked to El Nino too strongly. We show that by discarding that year, a regime shift in the 1990s is no longer visible. Overall, we observe a convergence of skill towards the present with correlations of about 0.4 for CFSv2 and of 0.6 for System 4 and SEAS5.
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.
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.
Three small stalagmites from Zoolithencave (southern Germany) show visible laminae, which consist of a clear and a brownish, pigmented layer pair. This potentially provides the opportunity to construct precise chronologies by counting annual laminae. The growth period of the three stalagmites was constrained by the C-14 bomb peak in the youngest part of all three stalagmites and C-14-dating of a piece of charcoal in the consolidated base part of stalagmite Zoo-rez-2. These data suggest an age of AD 1970 for the top laminae and a lower age limit of AD 1973-1682 or AD 1735-1778. Laminae were counted and their thickness determined on scanned thin sections of all stalagmites. On stalagmites Zoo-rez-1 and -2, three tracks were measured near the growth axes, each separated into three sections at prominent anchor laminae (I, II, III). Each section was replicated three times (a, b, c). For Zoo-rez-3, only one track was measured. The total number of laminae counted for Zoo-rez-1 ranges from 138 to 177, for Zoo-rez-2 from 119 to 145, and for Zoo-rez-3 from 159 to 166. The numbers agree well with the range constrained by the bomb peak and the age of the charcoal, which supports the annual origin of the laminae. The replicated measurements of the different tracks as well as the three different tracks on the stalagmites Zoo-rez-1 and-2 were cross-dated using the TSAP-Win (R) tree-ring software. This software is very useful for cross-dating because it enables to insert or delete missing or false laminae as well as identifying common pattern by shifting the series back and forth in time. However, visual inspection of the thin sections was necessary to confirm detection of missing or false laminae by TSAP-Win (R). For all three Zoo-rez speleothems, crossdating of the mean lamina thickness series was not possible due to a missing common pattern. The cross-dating procedure results in three refined chronologies for the three Zoo-rez stalagmites of ranging from AD 1821-1970 (Zoo-rez-1), AD 1835-1970 (Zoo-rez-2), and AD 1808-1970 (Zoo-rez-3).
Sociocultural valuation (SCV) of ecosystem services (ES) discloses the principles, importance or preferences expressed by people towards nature. Although ES research has increasingly addressed sociocultural values in past years, little effort has been made to systematically review the components of sociocultural valuation applications for different decision contexts (i.e. awareness raising, accounting, priority setting, litigation and instrument design). In this analysis, we investigate the characteristics of 48 different sociocultural valuation applications—characterised by unique combinations of decision context, methods, data collection formats and participants—across ten European case studies. Our findings show that raising awareness for the sociocultural value of ES by capturing people’s perspective and establishing the status quo, was found the most frequent decision context in case studies, followed by priority setting and instrument development. Accounting and litigation issues were not addressed in any of the applications. We reveal that applications for particular decision contexts are methodologically similar, and that decision contexts determine the choice of methods, data collection formats and participants involved. Therefore, we conclude that understanding the decision context is a critical first step to designing and carrying out fit-for-purpose sociocultural valuation of ES in operational ecosystem management.
Niche-based species distribution models (SDMs) have become an essential tool in conservation and restoration planning. Given the current threats to freshwater biodiversity, it is of fundamental importance to address scale effects on the performance of niche-based SDMs of freshwater species’ distributions. The scale effects are addressed here in the context of hierarchical catchment ordering, considered as counterpart to coarsening grain-size by increasing grid-cell size. We combine fish occurrence data from the Danube River Basin, the hierarchical catchment ordering and multiple environmental factors representing topographic, climatic and anthropogenic effects to model fish occurrence probability across multiple scales. We focus on 1st to 5th order catchments. The spatial scale (hierarchical catchment order) only marginally influences the mean performance of SDMs, however the uncertainty of the estimates increases with scale. Key predictors and their relative importance are scale and species dependent. Our findings have useful implications for choosing proper species dependent spatial scales for river rehabilitation measures, and for conservation planning in areas where fine grain species data are unavailable.
Different tectonic episodes from Late Triassic to recent times in the eastern Binalud Mountains have resulted from convergence and transpression between the Turan and Central Iran plates. Heterogeneous deformation and variable portions of pure and simple shear, demonstrated by finite strain and vorticity analysis in the Mashhad metamorphic rocks, indicate strain partitioning during the first tectonic episode. Modern strain partitioning is characterized by reverse and strike-slip faulting along the Neyshabur fault system and Shandiz fault zone in the southern and northern flanks of the eastern Binalud, respectively. Time-transgressive regional deformation migrated from the hinterland of the belt into the foreland basin, from northeast to the southwest of the mountains. Different generations of deformation resulted in obliteration of the subduction-related accretionary wedge, and growth of an orogenic wedge resulted from collision between the Central Iran and Turan plates.
Due to increasing demands for irrigation using groundwater as a source there is an urgent need for efficient methods that shed light on the resulting anthropogenic impacts on the connected aquifers. Thus an innovative approach is introduced, that aims to identify predominant geochemical changes in the groundwater system. The approach involves a principal component analysis as a promising tool to disentangle the effects of different impacts and even to give a quantitative assessment of the respective effect strength at each site. The study was applied in an irrigation region of the Nuthe River Basin, State Brandenburg, Northeast Germany. The results identify the negative impacts on the groundwater quality in the aquifer used for irrigation. A decrease of shallow groundwater quality under irrigation due to contamination with fertilizers (NO3, Cl, K, Na) and a slight shift in the redox system is indicated. Beside this direct impact on the shallow groundwater a long-term impact on a deeper groundwater resource could be identified. There is clear evidence, that the contamination is not restricted to the shallow groundwater but that extraction from deeper wells increasingly includes deeper, uncontaminated groundwater resources into the local irrigation cycle. The approach can be used as a basic tool for the adaptation of sustainable agricultural irrigation management strategies.
As a consequence of the rapid growing worldwide seismic data set, a huge variety of automatized data-processing methods have been developed. To perform automatized waveform-based seismological studies aiming for magnitudes or source process inversion, it is crucial to identify network stations with erroneous transfer functions, gain factors, or component orientations. We developed a new tool dedicated to automated station quality control of dense seismic networks and arrays. The python-based AutoStatsQ toolbox uses the pyrocko seismic data-processing environment. The toolbox automatically downloads data and metadata for selected teleseismic events and performs different tests. As a result, relative gain factors, sensor orientation corrections, and reliable frequency bands are computed for all stations in a chosen time period. Relative gain factors are calculated for all stations and events in a time domain based on maximum P-phase amplitudes. A Rayleigh-wave polarization analysis is used to identify deviating sensor orientations. The power spectra of all stations in a given frequency range are compared with synthetic ones, accessing Global Centroid Moment Tensor (CMT) solutions. Frequency ranges of coinciding synthetic and recorded power spectral densities (PSDs) may serve as guidelines for choosing band-pass filters for moment tensor (MT) inversion and help confirm the corner frequency of the instrument. The toolbox was applied to the permanent and temporary AlpArray networks as well as to the denser SWATH-D network, a total of over 750 stations. Stations with significantly deviating gain factors were identified, as well as stations with inverse polarity and misorientations of the horizontal components. The tool can be used to quickly access network quality and to omit or correct stations before MT inversion. Electronic Supplement: List of teleseismic events and tables of median, mean, and standard deviation of relative gain factors, and figures of relative gain factors of all event-station pairs, waveform example showing inverse polarity of horizontal components on ZS.D125, histograms of median, mean, and standard deviation of the correction angles, examples of synthetic and recorded frequency spectra of ZS.D046 and NI.VINO.
Small-scale variations in mineral chemistry, textures, and platinum group element (PGE) mineralization were investigated in the Lower and Middle Group chromitite layers LG6, LG6a, MG1, MG2, and MG2 II from vertical drill core profiles at the Thaba mine in the northwestern limb of the Bushveld Complex. We present detailed geochemical profiles of chromite composition and chromite crystal size distribution curves to shed light on the processes of chromite accumulation and textural modification as well as mineralization. Multiple samples within each layer were assayed for PGE concentrations, and the respective platinum group mineral association was determined by mineral liberation analysis (MLA). There is strong evidence for postcumulus changes in the chromitites. The crystal size distribution curves suggest that the primary chromite texture was coarsened by a combination of adcumulus growth and textural equilibration, while compaction of the crystal mush played only a minor role. Mineral compositions were also modified by postcumulus processes, but because of the very high modal amount of chromite and its local preservation in orthopyroxene oikocrysts, that phase retained much primary information. Vertical variations of chromite composition within chromitite layers and from one layer to another do not support the idea of chromite accumulation from crystal-rich slurries or crystal settling from a large magma chamber. Instead, we favor a successive buildup of chromitite layers by repeated injections of relatively thin layers of chromite-saturated magmas, with in situ crystallization occurring at the crystal mush-magma interface. The adcumulus growth of chromite grains to form massive chromitite required addition of Cr to the layers, which we attribute to downward percolation from the overlying magma. The PGE concentrations are elevated in all chromitite layers compared to adjacent silicate rocks and show a systematic increase upward from LG6 (avg 807 ppb Ir + Ru + Rh + Pt + Pd + Au) to MG2 II (avg 2,062 ppb). There are also significant internal variations in all layers, with enrichments at hanging and/or footwalls. The enriched nature of chromitites in PGEs compared to host pyroxenites is a general feature, independent of the layer thickness. The MLA results distinguish two principal groups of PGE mineral associations: the LG6, LG6, and MG1 are dominated by the malanite series, laurite, and PGE sulfarsenides, while the MG2 and MG2 II layers are characterized by laurite and PGE sulfides as well as Pt-Fe-Sn and PGE-Sb-Bi-Pb alloys. Differences in the PGE associations are attributed to postcumulus alteration of the MG2 and MG2 II layer, while the chromitites below, particularly LG6 and LG6a, contain a more pristine association.
Power-to-gas (PtG) stores chemical energy by converting excess electrical energy from renewable sources into an energy-dense gas. Due to its higher available capacity compared to surface-based storage technologies, subsurface storage in geological systems is the most promising approach for efficient and economic realization of the PtG system’s storage component. For this purpose, methane (CH4) produced by methanation by means of hydrogen (H2) and carbon dioxide (CO2) is stored in a geological reservoir until required for further use. In this context, CO2 is used as the cushion gas to maintain reservoir pressure and limiting working gas, i.e., (CH4) losses during withdrawal periods. Consequently, mixing of both gases in the reservoir is inevitable. Therefore, it is necessary to minimize the gas mixing region to optimize the efficiency of the PtG system’s storage component. In the present study, the physical properties of CH4, CO2 and their mixtures are reviewed. Then, a multicomponent flow model is implemented and validated against published data. Next, a hydromechanically coupled model is established, considering fluid flow through porous media and effective stresses to investigate the mixing behavior of both gases and the mechanical reservoir stability. The simulation results show that, with increasing reservoir thickness and dip angle, the mixing region is reduced during gas injection if CO2 is employed as the cushion gas. In addition, the degree of mixing is lower at higher temperatures. Feasible injection rates and injection schedules can be derived from the integrated reservoir stability analysis. The methodology developed in the present study allows the determination of optimum strategies for storage reservoir selection and gas injection scheduling by minimizing the gas mixing region.
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.
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.
Ancient evaporite deposits are geological archives of depositional environments characterized by a long‐term negative precipitation balance and bear evidence for global ocean element mass balance calculations. Here, Cretaceous selenite pseudomorphs from western Anatolia (‘Rosetta Marble’) — characterized by their exceptional morphological preservation — and their ‘marine’ geochemical signatures are described and interpreted in a process‐oriented context. These rocks recorded Late Cretaceous high‐pressure/low‐temperature, subduction‐related metamorphism with peak conditions of 1·0 to 1·2 GPa and 300 to 400°C. Metre‐scale, rock‐forming radiating rods, now present as fibrous calcite marble, clearly point to selenitic gypsum as the precursor mineral. Stratigraphic successions are recorded along a reconstructed proximal to distal transect. The cyclical alternation of selenite beds and radiolarian ribbon‐bedded cherts in the distal portions are interpreted as a two type of seawater system. During arid intervals, shallow marine brines cascaded downward into basinal settings and induced precipitation. During more humid times, upwelling‐induced radiolarian blooms caused the deposition of radiolarite facies. Interestingly, there is no comparable depositional setting known from the Cenozoic world. Meta‐selenite geochemical data (δ13C, δ18O and 87Sr/86Sr) plot within the range of reconstructed middle Cretaceous seawater signatures. Possible sources for the 13C‐enriched (mean 2·2‰) values include methanogenesis, gas hydrates and cold seep fluid exhalation. Spatially resolved component‐specific analysis of a rock slab displays isotopic variances between meta‐selenite crystals (mean δ13C 2·2‰) and host matrix (mean δ13C 1·3‰). The Cretaceous evaporite‐pseudomorphs of Anatolia represent a basin wide event coeval with the Aptian evaporites of the Proto‐Atlantic and the pseudomorphs share many attributes, including lateral distribution of 600 km and stratigraphic thickness of 1·5 to 2·0 km, with the evaporites formed during the younger Messinian salinity crisis. The Rosetta Marble of Anatolia may represent the best‐preserved selenite pseudomorphs worldwide and have a clear potential to act as a template for the study of meta‐selenite in deep time.
Millennial-scale Dansgaard Oeschger (DO) variability at northern high latitudes has influenced climatic and environmental conditions in the Mediterranean during the last glacial period. There is evidence that the hemispheric transmission of the DO variability occurred at the end of DO event 25; however, the exact timing and the trigger that activated the environmental response in the Mediterranean remains incompletely understood. Here, we provide evidence that the clear millennial-scale teleconnection between Greenland and the Mediterranean started at similar to 111.4 ka BP and was initiated by a sub-millennial scale cooling in Greenland (GI-25b). High-resolution sediment proxies and the pollen record of Lago Grande di Monticchio (MON), Italy, reflect climatic instability during the last millennium of the last interglacial, which was characterised by a first and short cooling episode (MON 1) at 111.44 +/- 0.69 ka BP, coinciding with the Greenland cold sub-event GI-25b in duration and timing (within dating uncertainties). MON and Greenland (NorthGRIP ice core) also agree in recording a subsequent warm rebound phase that abruptly culminated in the stadial MON 2/GS-25, marking the transition into the last glacial period. Our results show that the GI-25b triggered an early environmental response at MON to centennial-scale climate change in Greenland as a prelude to the millennial-scale teleconnection that was maintained during the glacial period.
Flood risk perceptions and the willingness to pay for flood insurance in the Veneto region of Italy
(2019)
The floodplain of the Veneto region (north-east Italy) is one of the most inhabited and economically competitive urban landscapes in Europe. Moreover, recent flood events have caused millions of Euros in damage across the region. Due to the combined influence of climate change and socio-economic development, flood impacts are expected to grow. Therefore, it is important for all flood-prone individuals to actively manage and limit flood risk through property-level flood risk management as part of an integrated flood risk management strategy. This is in line with the calls for wider community engagement in risk management in the Sendai Framework for Disaster Risk Reduction and the Sustainable Development Goals. Therefore, an online-survey of Veneto region residents was conducted asking questions regarding flood risk perceptions, preparedness, and preferences towards flood insurance via self-stated willingness to pay (WTP). Our analysis provides an initial indication that while flood risk knowledge is high, it may not be sufficient to encourage proactive risk management. From the WTP values provided people seem reluctant to buy insurance. However, many respondents expressed that a compulsory insurance system may be acceptable. In such a scheme the estimated insurance premium could fall to between (sic)26 and (sic)42 per year, as compared to, potentially, (sic)800 under risk-based premiums, which falls within the majority of WTP estimates provided ((sic)0-(sic)250). Overall, we identify areas of future research that are critical for the better design of risk management policies, supporting the insurance companies in risk management and for recommendations regarding property-level risk management.
Plate tectonics is among the most important geological processes on Earth, but its emergence and evolution remain unclear. Here we extrapolate models of present-day plate tectonics to the past and propose that since about three billion years ago the rise of continents and the accumulation of sediments at continental edges and in trenches has provided lubrication for the stabilization of subduction and has been crucial in the development of plate tectonics on Earth. We conclude that the two largest surface erosion and subduction lubrication events occurred after the Palaeoproterozoic Huronian global glaciations (2.45 to 2.2 billion years ago), leading to the formation of the Columbia supercontinent, and after the Neoproterozoic ‘snowball’ Earth glaciations (0.75 to 0.63 billion years ago). The snowball Earth event followed the ‘boring billion’—a period of reduced plate tectonic activity about 1.75 to 0.75 billion years ago that was probably caused by a shortfall of sediments in trenches—and it kick-started the modern episode of active plate tectonics.
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.
Isolated extreme habitats are ideally suited to investigate pivotal ecological processes such as niche use, local adaptation and dispersal. Extremophilic animals living in isolated habitats face the problem that dispersal is limited through the absence of suitable dispersal corridors, which in turn facilitates local adaptation. We used five rotifer isolates from extremely acidic mining lakes with a pH of below 3 as model organisms to test whether these isolates are acidotolerant or acidophilic, whether they survive and reproduce at their niche edges (here pH 2 and circum-neutral pH) and whether local adaptation has evolved. To evaluate potential dispersal limitation, we tested whether animals and their parthenogenetic eggs survive and remain reproductive or viable at unfavourable pH-conditions. All five isolates were acidophilic with a pH-optimum in the range of 4-6, which is well above the pH (< 3) of their lakes of origin. At unfavourable high pH, in four out of the five isolates parthenogenetic females produced a high number of non-viable eggs. Females and eggs produced at favourable pH (4) remained vital at an otherwise unfavourable pH of 7, indicating that for dispersal no acidic dispersal corridors are necessary. Common garden experiments revealed no clear evidence for local adaptation in any of the five isolates. Despite their acidophilic nature, all five isolates can potentially disperse via circum-neutral water bodies as long as their residence time is short, suggesting a broader dispersal niche than their realized niche. Local adaptation might have been hampered by the low population sizes of the rotifers in their isolated habitat and the short time span the mining lakes have existed.
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.
Empirical evidence of the relationship between social support and post-disaster mental health provides support for a general beneficial effect of social support (main-effect model; Wheaton, 1985). From a theoretical perspective, a buffering effect of social support on the negative relationship between disaster-related stress and mental health also seems plausible (stress-buffering model; Wheaton, 1985). Previous studies, however, (a) have paid less attention to the buffering effect of social support and (b) have mainly relied on interpersonal support (but not collective-level support such as community resilience) when investigating this issue. This previous work might have underestimated the effect of support on post-disaster mental health. Building on a sample of residents in Germany recently affected by flooding (N = 118), we show that community resilience to flooding (but not general interpersonal social support) buffered against the negative effects of flooding on post-disaster mental health. The results support the stress-buffering model and call for a more detailed look at the relationship between support and resilience and post-disaster adjustment, including collective-level variables.
The marine sedimentary record contains unique information about the history of erosion, uplift and climate of the adjacent continent. Inverting this record has been the purpose of many numerical studies. However, limited attention has been given to linking continental erosion to marine sediment transport and deposition in large-scale surface process evolution models. Here we present a new numerical method for marine sediment transport and deposition that is directly coupled to a landscape evolution algorithm solving for the continental fluvial and hillslope erosion equations using implicit and O(N) algorithms. The new method takes into account the sorting of grain sizes (e.g., silt and sand) in the marine domain using a non-linear multiple grain-size diffusion equation and assumes that the sediment flux exported from the continental domain is proportional to the bathymetric slope. Specific transport coefficients and compaction factors are assumed for the two different grain sizes to simulate the stratigraphic architecture. The resulting set of equations is solved using an efficient (O(N) and implicit) algorithm. It can thus be used to invert stratigraphic geometries using a Bayesian approach that requires a large number of simulations. This new method is used to invert the sedimentary geometry of a natural example, the Ogooue Delta (Gabon), over the last similar to 5 Myr. The objective is to unravel the set of erosional histories of the adjacent continental domain compatible with the observed geometry of the offshore delta. For this, we use a Bayesian inversion scheme in which the misfit function is constructed by comparing four geometrical parameters between the natural and the simulated delta: the volume of sediments stored in the delta, the surface slope, the initial and the final shelf lengths. We find that the best-fit values of the transport coefficients for silt in the marine domain are in the range of 300 - 500 m(2)/yr, in agreement with previous studies on offshore diffusion. We also show that, in order to fit the sedimentary geometry, erosion rate on the continental domain must have increased by a factor of 6 to 8 since 5.3 Ma. (C) 2019 Elsevier B.V. All rights reserved.
Climate and weather conditions in the mid-latitudes are strongly driven by the large-scale atmosphere circulation. Observational data indicate that important components of the large-scale circulation have changed in recent decades, including the strength and the width of the Hadley cell, jets, storm tracks and planetary waves. Here, we use a new statistical-dynamical atmosphere model (SDAM) to test the individual sensitivities of the large-scale atmospheric circulation to changes in the zonal temperature gradient, meridional temperature gradient and global-mean temperature. We analyze the Northern Hemisphere Hadley circulation, jet streams, storm tracks and planetary waves by systematically altering the zonal temperature asymmetry, the meridional temperature gradient and the global-mean temperature. Our results show that the strength of the Hadley cell, storm tracks and jet streams depend, in terms of relative changes, almost linearly on both the global-mean temperature and the meridional temperature gradient, whereas the zonal temperature asymmetry has little or no influence. The magnitude of planetary waves is affected by all three temperature components, as expected from theoretical dynamical considerations. The width of the Hadley cell behaves nonlinearly with respect to all three temperature components in the SDAM. Moreover, some of these observed large-scale atmospheric changes are expected from dynamical equations and are therefore an important part of model validation.
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.
Terrestrial gravimetry is increasingly used to monitor mass transport processes in geophysics boosted by the ongoing technological development of instruments. Resolving a particular phenomenon of interest, however, requires a set of gravity corrections of which the uncertainties have not been addressed up to now. In this study, we quantify the time domain uncertainty of tide, global atmospheric, large-scale hydrological, and nontidal ocean loading corrections. The uncertainty is assessed by comparing the majority of available global models for a suite of sites worldwide. The average uncertainty expressed as root-mean-square error equals 5.1nm/s(2), discounting local hydrology or air pressure. The correction-induced uncertainty of gravity changes over various time periods of interest ranges from 0.6nm/s(2) for hours up to a maximum of 6.7nm/s(2) for 6months. The corrections are shown to be significant and should be applied for most geophysical applications of terrestrial gravimetry. From a statistical point of view, however, resolving subtle gravity effects in the order of few nanometers per square second is challenged by the uncertainty of the corrections. Plain Language Summary Many scientists are exploring ways to benefit from gravity measurements in fields of high societal relevance such as monitoring of volcanoes or measuring the amount of water in underground. Any application of such new methods, however, requires careful preparation of the gravity measurements. The intention of the preparation process is to ensure that the measurements do not contain information about processes that are not of interest. For that reason, the influence of atmosphere, ocean, tides, and hydrology needs to be reduced from the gravity. In this study, we investigate how this reduction process influences the quality of the measurement. We found that the precision degrades especially owing to the hydrology. The ocean plays an important role at sites close to the coast and the atmosphere at sites located in mountains. The overall errors of the reductions may complicate a reliable use of gravity measurements in certain studies focusing on very small signals. Nevertheless, the precision of gravity reductions alone does not obstruct a meaningful use of gravity measurements in most research fields. Details specifying the reduction precision are provided in this study allowing scientist dealing with gravity measurements to decide if their signal of interest can be reliably resolved.
A comprehensive description of electromagnetic processes related to equatorial plasma depletions (EPDs) is essential for understanding their evolution and day-to-day variability. Recently, field-aligned currents (FACs) flowing at both western and eastern edges of EPDs were observed to be interhemispheric rather than anti-parallel about the dip equator, as suggested by previous theoretical studies. In this paper, we investigate the spatial and temporal variability of the FACs orientation using simultaneous measurements of electron density and magnetic field gathered by ESA’s Swarm constellation mission. By using empirical models, we assess the role of the Pedersen conductance in the preference of the FACs to close either in the northern or southern magnetic hemisphere. Here we show that the closure of the FACs agrees with an electrostatic regime determined by a hemispherical asymmetry of the Pedersen conductance. That is, the EPD-related FACs close at lower altitudes in the hemisphere with the highest conductivity. The evidence of this conclusion stands on the general agreement between the longitudinal and seasonal variability of both the conductivity and the FACs orientation.
Breakup Without Borders
(2019)
Relative plate motions during continental rifting result from the interplay of local with far-field forces. Here we study the dynamics of rifting and breakup using large-scale numerical simulations of mantle convection with self-consistent evolution of plate boundaries. We show that continental separation follows a characteristic evolution with four distinctive phases: (1) an initial slow rifting phase with low divergence velocities and maximum tensional stresses, (2) a synrift speed-up phase featuring an abrupt increase of extension rate with a simultaneous drop of tensional stress, (3) the breakup phase with inception of fast sea-floor spreading, and (4) a deceleration phase occurring in most but not all models where extensional velocities decrease. We find that the speed-up during rifting is compensated by subduction acceleration or subduction initiation even in distant localities. Our study illustrates new links between local rift dynamics, plate motions, and subduction kinematics during times of continental separation.
Interannual variations in ecosystem primary productivity are dominated by water availability. Until recently, characterizing the photosynthetic response of different ecosystems to soil moisture anomalies was hampered by observational limitations. Here, we use a number of satellite-based proxies for productivity, including spectral indices, sun-induced chlorophyll fluorescence, and data-driven estimates of gross primary production, to reevaluate the relationship between terrestrial photosynthesis and water. In contrast to nonwoody vegetation, we find a resilience of forested ecosystems to reduced soil moisture. Sun-induced chlorophyll fluorescence and data-driven gross primary production indicate an increase in photosynthesis as a result of the accompanying higher amounts of light and temperature despite lowered light-use-efficiency. Conversely, remote sensing indicators of greenness reach their detection limit and largely remain stable. Our study thus highlights the differential responses of ecosystems along a tree cover gradient and illustrates the importance of differentiating photosynthesis indicators from those of greenness for the monitoring and understanding of ecosystems. Plain Language Summary The capacity of vegetation to thrive and to sequester carbon depends on how much water they can have access to. In this work, we evaluate how different types of satellite observations can describe the response of vegetation to changes in soil moisture over the entire planet. The first source of observation measures only the greenness of the land surface, the second measures light that is emitted by pigments in plants which are photosynthetically active (chlorophyll fluorescence), and the third are simulations of gross carbon uptake derived from machine learning techniques. For periods of water shortage all three indicate a reduction of growth in ecosystems with few trees. However, in cold boreal forests, when soil moisture is particularly low, we still detect an increase in photosynthesis due to higher light and temperature conditions, but this is not reflected in the greenness indicator. This work illustrates how lack of water is not necessarily harmful for catching carbon through photosynthesis, but to monitor this effect, we need remote sensing indicators that measure more than just how green the plants are, and fluorescence is likely a good candidate.
Using optical absorption and Raman spectroscopic measurements, in conjunction with the first-principles calculations, a pressure-induced high-spin (HS)-to-low-spin (LS) state electronic transition of Fe2+ (M2-octahedral site) was resolved around 76-80GPa in a natural triphylite-lithiophilite sample with chemical composition (LiFe0.7082+Mn0.292PO4)-Li-M1-Fe-M2 (theoretical composition (LiFe0.52+Mn0.5PO4)-Li-M1-Fe-M2). The optical absorption spectra at ambient conditions consist of a broad doublet band with two constituents (1) (similar to 9330cm(-1)) and (2) (similar to 7110cm(-1)), resulting from the electronic spin-allowed transition (T2gEg)-T-5-E-5 of octahedral (HSFe2+)-Fe-M2. Both (1) and (2) bands shift non-linearly with pressure to higher energies up to similar to 55GPa. In the optical absorption spectrum measured at similar to 81GPa, the aforementioned HS-related bands disappear, whereas a new broadband with an intensity maximum close to 16,360cm(-1) appears, superimposed on the tail of the high-energy ligand-to-metal O2-Fe2+ charge-transfer absorption edge. We assign this new band to the electronic spin-allowed dd-transition (1)A(1g)(1)T(1g) of LS Fe2+ in octahedral coordination. The high-pressure Raman spectra evidence the Fe2+ HS-to-LS transition mainly from the abrupt shift of the P-O symmetric stretching modes to lower frequencies at similar to 76GPa, the highest pressure achieved in the Raman spectroscopic experiments. Calculations indicated that the presence of Mn-M2(2+) simply shifts the isostructural HS-to-LS transition to higher pressures compared to the triphylite Fe-M2(2+) end-member, in qualitative agreement with our experimental observations.
Proxy records from climate archives provide evidence about past climate changes, but the recorded signal is affected by non-climate-related effects as well as time uncertainty. As proxy-based climate reconstructions are frequently used to test climate models and to quantitatively infer past climate, we need to improve our understanding of the proxy record signal content as well as the uncertainties involved. In this study, we empirically estimate signal-to-noise ratios (SNRs) of temperature proxy records used in global compilations of the middle to late Holocene (last 6000 years). This is achieved through a comparison of the correlation of proxy time series from nearby sites of three compilations and model time series extracted at the proxy sites from two transient climate model simulations: a Holocene simulation of the ECHAM5/MPI-OM model and the Holocene part of the TraCE-21ka simulation. In all comparisons, we found the mean correlations of the proxy time series on centennial to millennial timescales to be low (R < 0.2), even for nearby sites, which resulted in low SNR estimates. The estimated SNRs depend on the assumed time uncertainty of the proxy records, the timescale analysed, and the model simulation used. Using the spatial correlation structure of the ECHAM5/MPI-OM simulation, the estimated SNRs on centennial timescales ranged from 0.05 - assuming no time uncertainty - to 0.5 for a time uncertainty of 400 years. On millennial timescales, the estimated SNRs were generally higher. Use of the TraCE-21ka correlation structure generally resulted in lower SNR estimates than for ECHAM5/MPI-OM. As the number of available high-resolution proxy records continues to grow, a more detailed analysis of the signal content of specific proxy types should become feasible in the near future. The estimated low signal content of Holocene temperature compilations should caution against over-interpretation of these multi-proxy and multisite syntheses until further studies are able to facilitate a better characterisation of the signal content in paleoclimate records.
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.
The Last Interglacial (Eemian, MIS 5e) can be considered a test-bed for climate dynamics under a warmer-than-present climate. In this study we present a chironomid record from the high latitude Sokli site (N Finland), where a long continuous sediment sequence from the last interglacial has been preserved from glacial erosion. The chironomid-analysis shows a diverse fauna, with dominance of warm-water indicators and shifts in assemblage composition that can be attributed to temperature, lake depth, productivity and habitat availability. Quantitative mean July paleotemperature estimates based on the chironomid data indicate overall mean July air temperatures up to 1 degrees C warmer than present. Two cooling events can be discerned, the Tunturi event, dated to about 127.5kaBP, in the lower part of the sequence, and the Varrio event, dated to about 119kaBP, associated with the beginning of a cooling trend in the upper part of the record. Warm conditions already at the onset of the interglacial contrast with a recent chironomid-based last interglacial temperature reconstruction from Denmark, which suggests a late onset of Eemian warming. The relatively small increase in inferred temperatures compared to present day temperatures at Sokli differs from other high latitude Eemian sites, and likely reflects the influence of the Atlantic Meridional Overturning Circulation in maintaining already elevated temperatures in Fennoscandia during interglacials.
Quantifying erosion rates, and how they compare to rock uplift rates, is fundamental for understanding landscape response to tectonics and associated sediment fluxes from upland areas. The erosional response to uplift is well-represented by river incision and the associated landslide activity. However, characterising the relationship between these processes remains a major challenge in tectonically active areas, in some cases because landslides can preclude obtaining reliable erosion rates from cosmogenic radionuclide (CRN) concentrations. Here, we quantify the control of tectonics and its coupled geomorphic response on the erosion rates of catchments in southern Italy that are experiencing a transient response to normal faulting. We analyse in-situ Be-10 concentrations for detrital sediment samples, collected along the strike of faults with excellent tectonic constraints and landslide inventories. We demonstrate that Be-10-derived erosion rates are controlled by fault throw rates and the extent of transient incision and associated landsliding in the catchments. We show that the low-relief sub-catchments above knickpoints erode at uniform background rates of similar to 0.10 mm/yr, while downstream of knickpoints, erosion removes similar to 50% of the rock uplifted by the faults, at rates of 0.10-0.64 mm/yr. Despite widespread landsliding, CRN samples provide relatively consistent and accurate erosion rates, most likely because landslides are frequent, small, and shallow, and represent the integrated record of landsliding over several seismic cycles. Consequently, we combine these validated Be-10 erosion rates and data from a geomorphological landslide inventory in a published numerical model, to gain further insight into the long-term landslide rates and sediment mixing, highlighting the potential of CRN data to study landslide dynamics. (C) 2018 Elsevier B.V. All rights reserved.
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.
SiO(2 )is the main component of silicate melts and thus controls their network structure and physical properties. The compressibility and viscosities of melts at depth are governed by their short range atomic and electronic structure. We measured the O K-edge and the Si L-2,L-3-edge in silica up to 110 GPa using X-ray Raman scattering spectroscopy, and found a striking match to calculated spectra based on structures from molecular dynamic simulations. Between 20 and 27 GPa, Si-[4] species are converted into a mixture of Si-[5] and Si-[6] species and between 60 and 70 GPa, Si-[6] becomes dominant at the expense of Si-[5] with no further increase up to at least 110 GPa. Coordination higher than 6 is only reached beyond 140 GPa, corroborating results from Brillouin scattering. Network modifying elements in silicate melts may shift this change in coordination to lower pressures and thus magmas could be denser than residual solids at the depth of the core-mantle boundary.
Metamorphic geology
(2019)
Results of a paleolimnological investigation of a well-dated lake sediment section from Shikotan Island (Southern Kurils) showed that from ca 8.0 to 5.8 cal ka BP a warm and humid period corresponding to middle Holocene optimum took place. Cooling thereafter corresponds to Neoglacial. A reconstructed from ca 0.9 to ca 0.58 cal ka BP warm period can correspond to a Medieval Warm Period. Cooling after 0.58cal ka BP can be correlated with the LIA. Marine regression stages were identified at ca 6.2-5.9, 5.5-5.1 and 1.07-0.36 cal ka BP. The general chronology of major climatic events of Holocene in the island is in accordance with the climate records from the North Pacific region. Revealed spatial differences in timing and magnitude of the Late Holocene climatic episodes (LIA, MWP) in the region needs further investigations.
During eruptive activity of andesitic stratovolcanoes, the extrusion of lava domes, their collapse and intermittent explosions are common volcanic hazards. Many lava domes grow in a preferred direction, in turn affecting the direction of lava flows and pyroclastic density currents. Access to active lava domes is difficult and hazardous, so detailed data characterizing lava dome growth are typically limited, keeping the processes controlling the directionality of extrusions unclear. Here we combine TerraSAR-X satellite radar observations with high-resolution airborne photogrammetry to assess morphological changes, and perform finite element modeling to investigate the impact of loading stress on shallow magma ascent directions associated with lava dome extrusion and crater formation at Volcan de Colima, Mexico. The TerraSAR-X data, acquired in similar to 1-m resolution spotlight mode, enable us to derive a chronology of the eruptive processes from intensity-based time-lapse observations of the general crater and dome evolution. The satellite images are complemented by close-range airborne photos, processed by the Structure-from-Motion workflow. This allows the derivation of high-resolution digital elevation models, providing insight into detailed loading and unloading features. During the observation period from Jan-2013 to Feb-2016, we identify a dominantly W-directed dome growth and lava flow production until Jan-2015. In Feb-2015, following the removal of the active summit dome, the surface crater widened and elongated along a NE-SW axis. Later in May-2015, a new dome grew toward the SW of the crater while a separate vent developed in the NE of the crater, reflecting a change in the direction of magma ascent and possible conduit bifurcation. Finite element models show a significant stress change in agreement with the observed magma ascent direction changes in response to the changing surface loads, both for loading (dome growth) and unloading (crater forming excavation) cases. These results allow insight into shallow dome growth dynamics and the migration of magma ascent in response to changing volcano summit morphology. They further highlight the importance of detailed volcano summit morphology surveillance, as changes in direction or location of dome extrusion may have major implications regarding the directions of potential volcanic hazards, such as pyroclastic density currents generated by dome collapse.
Glacial-interglacial variations in CO2 and methane in polar ice cores have been attributed, in part, to changes in global wetland extent, but the wetland distribution before the Last Glacial Maximum (LGM, 21 ka to 18 ka) remains virtually unknown. We present a study of global peatland extent and carbon (C) stocks through the last glacial cycle (130 ka to present) using a newly compiled database of 1,063 detailed stratigraphic records of peat deposits buried by mineral sediments, as well as a global peatland model. Quantitative agreement between modeling and observations shows extensive peat accumulation before the LGM in northern latitudes (> 40 degrees N), particularly during warmer periods including the last interglacial (130 ka to 116 ka, MIS 5e) and the interstadial (57 ka to 29 ka, MIS 3). During cooling periods of glacial advance and permafrost formation, the burial of northern peatlands by glaciers and mineral sediments decreased active peatland extent, thickness, and modeled C stocks by 70 to 90% from warmer times. Tropical peatland extent and C stocks show little temporal variation throughout the study period. While the increased burial of northern peats was correlated with cooling periods, the burial of tropical peat was predominately driven by changes in sea level and regional hydrology. Peat burial by mineral sediments represents a mechanism for long-term terrestrial C storage in the Earth system. These results show that northern peatlands accumulate significant C stocks during warmer times, indicating their potential for C sequestration during the warming Anthropocene.
The near-polar orbit satellites of Swarm mission provide a good opportunity to investigate the conjugacy of equatorial plasma irregularities (EPIs) since their trajectories at low latitudes are basically aligned with fixed geographical longitude. However, the Swarm in situ electron density occasionally shows EPIs at only one hemisphere at this longitude. In this study, we provide detailed analysis of such EPI events from the in situ electron densities and onboard global positioning system (GPS) measurements of Swarm low pair satellites, and simultaneous GPS data from two geomagnetically conjugate ground stations at the Africa longitudes. The result indicates that when Swam in situ electron density sometime shows EPIs at only one hemisphere, the GPS scintillations are still observed from the Swarm onboard receiver and by the two conjugate ground stations. It implies that the EPIs should generally elongate along the geomagnetic flux tube. More than two-year statistic results show that the onset time of scintillation in the northern station is on average 16 and 18 min earlier than that in the southern station for September equinox and December solstice in 2015, while for March equinox in 2016 the onset time of scintillation of northern station is about 11 min later than that of southern station, which indicates the asymmetry features of EPIs along the flux tube. Further analysis of nearly three-year GPS data from two conjugate stations at the Asia longitudes, we find that during solar maximum years the local sunset time plays an important role for causing the difference of onset time of scintillation between two conjugate stations.
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.
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.
Following the Middle Permian (Capitanian) mass extinction there was a global ‘reef eclipse’, and this event had an important role in the Paleozoic-Mesozoic transition of reef ecosystems. Furthermore, the recovery pattern of reef ecosystems in the Wuchiapingian of South China, before the radiation of Changhsingian reefs, is poorly understood. Here, we present a detailed sedimentological account of the Tieqiao section, South China, which records the only known Wuchiapingian reef setting from South China. Six reef growing phases were identified within six transgressive-regressive cycles. The cycles represent changes of deposition in a shallow basin to a subtidal outer platform setting, and the reefal build-ups are recorded in the shallowest part of the cycles above wave base in the euphotic zone. Our results show that the initial reef recovery started from the shallowing up part of the 1st cycle, within the Clarkina leveni conodont zone, which is two conodont zones earlier than previously recognized. In addition, even though metazoans, such as sponges, do become important in the development of the reef bodies, they are not a major component until later in the Wuchiapingian in the 5th and 6th transgressive-regressive cycles. This suggests a delayed recovery of metazoan reef ecosystems following the Middle Permian extinction. Furthermore, even though sponges do become abundant within the reefs, it is the presence and growth of the encrusters Archaeolithoporella and Tubiphytes and abundance of microbial micrites that play an important role in stabilizing the reef structures that form topographic highs.
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).
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.
A well-preserved sequence, by Archean standards, of mantle-derived metabasalts and metakomatiites forms large parts of the lower Onverwacht Group of the Barberton Greenstone Belt (South Africa). To elucidate the origin of mafic and ultramafic rocks from this 3.55 to 3.45 Ga sequence, we present a comprehensive geochemical dataset including major and trace elements as well as Lu-Hf and Sm-Nd isotope compositions for a variety of metavolcanic rocks. These include metabasalts of the amphibolite-facies Sandspruit and Theespruit Formations as well as metabasalts and metakomatiites of the lower greenschist-facies Komati Formation. Based on their incompatible trace element patterns, the basalts of the Sandspruit and Theespruit Formations can be subdivided into a light rare earth element (LREE) depleted group, a LREE-undepleted group, and a LREE-enriched group. Positive epsilon Hf-(t) and epsilon Nd-(t) values of ca. +3 to +4 and 0 to +2, respectively, together with depletions in Th and La-CN/Yb-CN indicate derivation of the LREE-depleted basalts from a depleted mantle source. However, chondritic epsilon Hf-(t) and epsilon Nd-(t) values combined with positive Th and La-CN/Yb-CN of the LREE-enriched samples indicate a contribution from older granitoid crust in the petrogenesis of these samples. Trace element patterns of komatiites and basalts of the Komati Formation are generally flat relative to primitive mantle with slight depletions in heavy rare earth elements and Th and overall positive epsilon Hf-(t) of + 2.5 +/- 3.5 (2 s.d.) and epsilon Nd-(t) of + 0.5 +/- 2.2 (2 s. d.). The coherence in trace element characteristics suggests a common magmatic origin for basalts and komatiites. This study reveals that the two lavas were derived from the same mantle plume, i. e. komatiites were formed by high degrees of melting of a depleted mantle source containing residual garnet and the basalts were formed by moderate degrees of partial melting in shallower regions of the mantle. Based on the current dataset, combined with published data, we propose a geodynamic model for the oldest units of the Barberton Greenstone Belt that describes the development from a submerged continental setting (for the Sandspruit and Theespruit Formations) to a submarine plateau setting (for the Komati Formation) as a consequence of continental rifting.
Evolution of Large-Scale Magnetic Fields From Near-Earth Space During the Last 11 Solar Cycles
(2019)
We use hourly mean magnetic field measurements from 34 midlatitude geomagnetic observatories between 1900 and 2015 to investigate the long-term evolution and driving mechanism of the large-scale external magnetic field at ground. The Hourly Magnetospheric Currents index (HMC) is derived as a refinement of the Annual Magnetospheric Currents index (HMC, Pick & Korte, 2017, https://doi.org/10.1093/gji/ggx367). HMC requires an extensive revision of the observatory hourly means. It depends on three third party geomagnetic field models used to eliminate the core, the crustal, and the ionospheric solar-quiet field contributions. We mitigate the dependency of HMC on the core field model by subtracting only nondipolar components of the model from the data. The separation of the residual (dipolar) signal into internal and external (HMC) parts is the main methodological challenge. Observatory crustal biases are updated with respect to AMC, and the solar-quiet field estimation is extended to the past based on a reconstruction of solar radio flux (F10.7). We find that HMC has more power at low frequencies (periods = 1 year) than the Dcx index, especially at periods relevant to the solar cycle. Most of the slow variations in HMC can be explained by the open solar magnetic flux. There is a weakly decreasing linear trend in absolute HMC from 1900 to present, which depends sensitively on the data rejection criteria at early years. HMC is well suited for studying long-term variations of the geomagnetic field.
Underground coal gasification (UCG) enables utilization of coal reserves, currently not economically exploitable due to complex geological boundary conditions. Hereby, UCG produces a high-calorific synthesis gas that can be used for generation of electricity, fuels, and chemical feedstock. The present study aims to identify economically-competitive, site-specific end-use options for onshore- and offshore-produced UCG synthesis gas, taking into account the capture and storage (CCS) and/or utilization (CCU) of produced CO2. Modeling results show that boundary conditions favoring electricity, methanol, and ammonia production expose low costs for air separation, low compression power requirements, and appropriate shares of H-2/N-2. Hereby, a gasification agent ratio of more than 30% oxygen by volume is not favorable from the economic and CO2 mitigation viewpoints. Compared to the costs of an offshore platform with its technical equipment, offshore drilling costs are marginal. Thus, uncertainties related to parameters influenced by drilling costs are negligible. In summary, techno-economic process modeling results reveal that air-blown gasification scenarios are the most cost-effective ones, while offshore UCG-CCS/CCU scenarios are up to 1.7 times more expensive than the related onshore processes. Hereby, all investigated onshore scenarios except from ammonia production under the assumed worst-case conditions are competitive on the European market.
Ferns are among the most popular groups of plants in the tropics and subtropics, and their role as carbon sequestrators has been widely recognized. However, there is little understanding of the silicaceous structures (phytoliths) of ferns, rate of phytolith turnover, the consequences for organic matter sequestered in phytoliths and consequences for other soil properties. In the study reported here, high-resolution X-ray tomographic microscopy and chemical characterization were applied to examine the traits of phytoliths of the fern Dicranopteris linearis (Burm.f.) Underw. (D. linearis), with a focus on their dissolution properties and accumulation in northern Vietnamese soils in relation to soil properties. Tomographic images revealed an inter-embedding structure of silica and organic matter, especially in leaf-derived material. We propose that organic matter and silica can preserve each other against decomposition. In batch experiments, there was a relatively small rate of dissolution of phytoliths with dry ashing and subsequent H2O2 treatment. Silicon (Si) dissolution for D. linearis phytolith samples was much less than that for rice phytoliths. Despite the fact that the aluminum (Al) content was large in D. linearis leaves, batch dissolution data did not confirm a relation between Al and the slow rate of phytolith dissolution. The soil phytolith content varied from 0.9 to 7.5 g kg(-1) in the topsoil across the mountainous areas in northern Vietnam, whereas it tended to be smaller in the subsoil. The data indicate a relation between phytolith and soil organic matter, clay content, oxalate-soluble Al and electrical conductivity, suggesting that these soil properties are among the important factors affecting the size of the soil phytolith Si pool. Highlights
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.
Increased landslide activity on forested hillslopes following two recent volcanic eruptions in Chile
(2019)
Large explosive eruptions can bury landscapes beneath thick layers of tephra. Rivers subsequently overloaded with excess pyroclastic sediments have some of the highest reported specific sediment yields. Much less is known about how hillslopes respond to tephra loads. Here, we report a pulsed and distinctly delayed increase in landslide activity following the eruptions of the Chaiten (2008) and Puyehue-Cordon Caulle (2011) volcanoes in southern Chile. Remote-sensing data reveal that land-slides clustered in densely forested hillslopes mostly two to six years after being covered by tephra. This lagged instability is consistent with a gradual loss of shear strength of decaying tree roots in areas of high tephra loads. Surrounding areas with comparable topography, forest cover, rainfall and lithology maintained landslide rates roughly ten times lower. The landslides eroded the landscape by up to 4.8 mm on average within 30 km of both volcanoes, mobilizing up to 1.6 MtC at rates of about 265 tC km(-2) yr(-1). We suggest that these yields may reinforce the elevated river loads of sediment and organic carbon in the decade after the eruptions. We recommend that studies of post-eruptive mass fluxes and hazards include lagged landslide responses of tephra-covered forested hillslopes, to avoid substantial underestimates.
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 show that near-real-time seismic monitoring of fluid injection allowed control of induced earthquakes during the stimulation of a 6.1-km-deep geothermal well near Helsinki, Finland. A total of 18,160 m(3) of fresh water was pumped into crystalline rocks over 49 days in June to July 2018. Seismic monitoring was performed with a 24-station borehole seismometer network. Using near-real-time information on induced-earthquake rates, locations, magnitudes, and evolution of seismic and hydraulic energy, pumping was either stopped or varied-in the latter case, between well-head pressures of 60 and 90 MPa and flow rates of 400 and 800 liters/min. This procedure avoided the nucleation of a project-stopping magnitude M-W 2.0 induced earthquake, a limit set by local authorities. Our results suggest a possible physics-based approach to controlling stimulation-induced seismicity in geothermal projects.
We analyze a large transient strainmeter signal recorded at 62.5 m depth along the southern shore of the eastern Sea of Marmara region in northwestern Turkey. This region represents a passage of stress transfer from the Izmit rupture to the Marmara seismic gap. The strain signal was recorded at the Esenkoy site by one of the ICDP-GONAF (International Continental Drilling Programme - Geophysical Observatory at the North Anatolian Fault) strainmeters on the Armutlu peninsula with a maximum amplitude of 5 microstrain and lasting about 50 days. The onset of the strain signal coincided with the origin time of a M-w 4.4 earthquake offshore Yalova, which occurred as part of a seismic sequence including eight M-w >= 3.5 earthquakes. The Mw 4.4 event occurred at a distance of about 30 km from Esenkoy on June 25th 2016 representing the largest earthquake in this region since 2008. Before the event, the maximum horizontal strain was subparallel to the regional maximum horizontal stress derived from stress inversion of local seismicity. During the strain transient, we observe a clockwise rotation in the local horizontal strain field of about 20 degrees. The strain signal does not correlate with known environmental parameters such as annual changes of sea level, rainfall or temperature. The strain signal could indicate local slow slip on the Cinarcik fault and thus a transfer of stress to the eastern Marmara seismic gap.
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 intangible impacts of floods on welfare are not well investigated, even though they are important aspects of welfare. Moreover, flooding has gender based impacts on welfare. These differing impacts create a gender based flood risk resilience gap. We study the intangible impacts of flood risk on the subjective well-being of residents in central Vietnam. The measurement of intangible impacts through subjective well-being is a growing field within flood risk research. We find an initial drop in welfare through subjective well-being across genders when a flood is experienced. Male respondents tended to recover their welfare losses by around 80% within 5 years while female respondents were associated with a welfare recovery of around 70%. A monetization of the impacts floods have on an individual’s subjective well-being shows that for the average female respondent, between 41% to 86% of annual income would be required to compensate subjective well-being losses after 5 years of experiencing a flood. The corresponding value for males is 30% to 57% of annual income. This shows that the intangible impacts of flood risk are important (across genders) and need to be integrated into flood (or climate) risk assessments to develop more socially appropriate risk management strategies.
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.
The Arctic is directly impacted by climate change. The increase in air temperature drives the thawing of permafrost and an increase in coastal erosion and river discharge. This leads to a greater input of sediment and organic matter into coastal waters, which substantially impacts the ecosystems, the subsistence economy of the local population, and the climate because of the transformation of organic matter into greenhouse gases. Yet, the patterns of sediment dispersal in the nearshore zone are not well known, because ships do not often reach shallow waters and satellite remote sensing is traditionally focused on less dynamic environments. The goal of this study is to use the extensive Landsat archive to investigate sediment dispersal patterns specifically on an exemplary Arctic nearshore environment, where field measurements are often scarce. Multiple Landsat scenes were combined to calculate means of sediment dispersal and sea surface temperature under changing seasonal wind conditions in the nearshore zone of Herschel Island Qikiqtaruk in the western Canadian Arctic since 1982. We use observations in the Landsat red and thermal wavebands, as well as a recently published water turbidity algorithm to relate archive wind data to turbidity and sea surface temperature. We map the spatial patterns of turbidity and water temperature at high spatial resolution in order to resolve transport pathways of water and sediment at the water surface. Our results show that these pathways are clearly related to the prevailing wind conditions, being ESE and NW. During easterly wind conditions, both turbidity and water temperature are significantly higher in the nearshore area. The extent of the Mackenzie River plume and coastal erosion are the main explanatory variables for sediment dispersal and sea surface temperature distributions in the study area. During northwesterly wind conditions, the influence of the Mackenzie River plume is negligible. Our results highlight the potential of high spatial resolution Landsat imagery to detect small-scale hydrodynamic processes, but also show the need to specifically tune optical models for Arctic nearshore environments.
Thermokarst results from the thawing of ice-rich permafrost and alters the biogeochemical cycling in the Arctic by reworking soil material and redistributing soil organic carbon (SOC) and total nitrogen (TN) along uplands, hillslopes, and lowlands. Understanding the impact of this redistribution is key to better estimating the storage of SOC in permafrost terrains. However, there are insufficient studies quantifying long-term impacts of thaw processes on the distribution of SOC and TN along hillslopes. We address this issue by providing estimates of SOC and TN stocks along the hillslopes of three valleys located on Herschel Island (Yukon, Canada), and by discussing the impact of hillslope thermokarst on the variability of SOC and TN stocks. We found that the average SOC and TN 0-100 cm stocks in the valleys were 26.4 +/- 8.9 kg C m(-2) and 2.1 +/- 0.6 kg N m(-2). We highlight the strong variability in the soils physical and geochemical properties within hillslope positions. High SOC stocks were found at the summits, essentially due to burial of organic matter by cryoturbation, and at the toeslopes due to impeded drainage which favored peat formation and SOC accumulation. The average carbon-to-nitrogen ratio in the valleys was 12.9, ranging from 9.7 to 18.9, and was significantly higher at the summits compared to the backslopes and footslopes (p < 0.05), suggesting a degradation of SOC downhill. Carbon and nitrogen contents and stocks were significantly lower on 16% of the sites that were previously affected by hillslope thermokarst (p < 0.05). Our results showed that lateral redistribution of SOC and TN due to hillslope thermokarst has a strong impact on the SOC storage in ice-rich permafrost terrains.
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.
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%.
When dealing with issues that are of high societal relevance, Earth sciences still face a lack of acceptance, which is partly rooted in insufficient communication strategies on the individual and local community level. To increase the efficiency of communication routines, science has to transform its outreach concepts to become more aware of individual needs and demands. The “encoding/decoding” concept as well as critical intercultural communication studies can offer pivotal approaches for this transformation.
The Arctic-Boreal regions experience strong changes of air temperature and precipitation regimes, which affect the thermal state of the permafrost. This results in widespread permafrost-thaw disturbances, some unfolding slowly and over long periods, others occurring rapidly and abruptly. Despite optical remote sensing offering a variety of techniques to assess and monitor landscape changes, a persistent cloud cover decreases the amount of usable images considerably. However, combining data from multiple platforms promises to increase the number of images drastically. We therefore assess the comparability of Landsat-8 and Sentinel-2 imagery and the possibility to use both Landsat and Sentinel-2 images together in time series analyses, achieving a temporally-dense data coverage in Arctic-Boreal regions. We determined overlapping same-day acquisitions of Landsat-8 and Sentinel-2 images for three representative study sites in Eastern Siberia. We then compared the Landsat-8 and Sentinel-2 pixel-pairs, downscaled to 60 m, of corresponding bands and derived the ordinary least squares regression for every band combination. The acquired coefficients were used for spectral bandpass adjustment between the two sensors. The spectral band comparisons showed an overall good fit between Landsat-8 and Sentinel-2 images already. The ordinary least squares regression analyses underline the generally good spectral fit with intercept values between 0.0031 and 0.056 and slope values between 0.531 and 0.877. A spectral comparison after spectral bandpass adjustment of Sentinel-2 values to Landsat-8 shows a nearly perfect alignment between the same-day images. The spectral band adjustment succeeds in adjusting Sentinel-2 spectral values to Landsat-8 very well in Eastern Siberian Arctic-Boreal landscapes. After spectral adjustment, Landsat and Sentinel-2 data can be used to create temporally-dense time series and be applied to assess permafrost landscape changes in Eastern Siberia. Remaining differences between the sensors can be attributed to several factors including heterogeneous terrain, poor cloud and cloud shadow masking, and mixed pixels.
High spectral resolution (hyperspectral) remote sensing has already demonstrated its capabilities for soil constituent mapping based on absorption feature parameters. This paper tests different parametrizations of the 1.75 μm gypsum feature for the determination of gypsum abundances, from the laboratory to remote sensing applications of recent as well as upcoming hyperspectral sensors. In particular, this study focuses on remote sensing imagery over the large body of the Omongwa pan located in the Namibian Kalahari. Four common absorption feature parameters are compared: band ratio through the introduction of the Normalized Differenced Gypsum Index (NDGI), the shape-based parameters Slope, and Half-Area, and the Continuum Removed Absorption Depth (CRAD). On laboratory soil samples from the pan, CRAD and NDGI approaches perform best to determine gypsum content tested in cross validated regression models with XRD mineralogical data (R² = 0.84 for NDGI and R² = 0.86 for CRAD). Subsequently the laboratory prediction functions are transferred to remote sensing imagery of spaceborne Hyperion, airborne HySpex and simulated spaceborne EnMAP sensor. Variable results were obtained depending on sensor characteristics, data quality, preprocessing and spectral parameters. Overall, the CRAD parameter in this wavelength region proved not to be robust for remote sensing applications, and the simple band ratio based parameter, the NDGI, proved robust and is recommended for future use for the determination of gypsum content in bare soils based on remote sensing hyperspectral imagery.
With the advent of the two Sentinel-1 (S1) satellites, Synthetic Aperture Radar (SAR) data with high temporal and spatial resolution are freely available. This provides a promising framework to facilitate detailed investigations of surface instabilities and movements on large scales with high temporal resolution, but also poses substantial processing challenges because of storage and computation requirements. Methods are needed to efficiently detect short term changes in dynamic environments. Approaches considering pair-wise processing of a series of consecutive scenes to retain maximum temporal resolution in conjunction with time series analyses are required. Here we present OSARIS, the “Open Source SAR Investigation System,” as a framework to process large stacks of S1 data on high-performance computing clusters. Based on Generic Mapping Tools SAR, shell scripts, and the workload manager Slurm, OSARIS provides an open and modular framework combining parallelization of high-performance C programs, flexible processing schemes, convenient configuration, and generation of geocoded stacks of analysis-ready base data, including amplitude, phase, coherence, and unwrapped interferograms. Time series analyses can be conducted by applying automated modules to the data stacks. The capabilities of OSARIS are demonstrated in a case study from the northwestern Tien Shan, Central Asia. After merging of slices, a total of 80 scene pairs were processed from 174 total input scenes. The coherence time series exhibits pronounced seasonal variability, with relatively high coherence values prevailing during the summer months in the nival zone. As an example of a time series analysis module, we present OSARIS' “Unstable Coherence Metric” which identifies pixels affected by significant drops from high to low coherence values. Measurements of motion provided by LOSD measurements require careful evaluation because interferometric phase unwrapping is prone to errors. Here, OSARIS provides a series of modules to detect and mask unwrapping errors, correct for atmospheric disturbances, and remove large-scale trends. Wall clock processing time for the case study (area ~9,000 km2) was ~12 h 4 min on a machine with 400 cores and 2 TB RAM. In total, ~12 d 10 h 44 min (~96%) were saved through parallelization. A comparison of selected OSARIS datasets to results from two state-of-the-art SAR processing suites, ISCE and SNAP, shows that OSARIS provides products of competitive quality despite its high level of automatization. OSARIS thus facilitates efficient S1-based region-wide investigations of surface movement events over multiple years.
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.
Flow accumulation algorithms estimate the steady state of flow on real or modeled topographic surfaces and are crucial for hydrological and geomorphological assessments, including delineation of river networks, drainage basins, and sediment transport processes. Existing flow accumulation algorithms are typically designed to compute flows on regular grids and are not directly applicable to arbitrarily sampled topographic data such as lidar point clouds. In this study we present a random sampling scheme that generates homogeneous point densities, in combination with a novel flow path tracing approach-the Facet-Flow Network (FFN)-that estimates flow accumulation in terms of specific catchment area (SCA) on triangulated surfaces. The random sampling minimizes biases due to spatial sampling and the FFN allows for direct flow estimation from point clouds. We validate our approach on a Gaussian hill surface and study the convergence of its SCA compared to the analytical solution. Here, our algorithm outperforms the multiple flow direction algorithm, which is optimized for divergent surfaces. We also compute the SCA of a 6-km(2)-steep, vegetated catchment on Santa Cruz Island, California, based on airborne lidar point-cloud data. Point-cloud-based SCA values estimated by our method compare well with those estimated by the D-infinity or multiple flow direction algorithm on gridded data. The advantage of computing SCA from point clouds becomes relevant especially for divergent topography and for small drainage areas: These are depicted with much more detail due to the higher sampling density of point clouds.
Caldera unrest can lead to major volcanic eruptions. Analysis of subtle subsidence or inflation at calderas helps understanding of their subsurface volcanic processes and related hazards. Several subsiding calderas have shown similar patterns of ground deformation composed of broad subsidence affecting the entire volcanic edifice and stronger localized subsidence focused inside the caldera. Physical models of internal deformation sources used to explain these observations typically consist of two magma reservoirs at different depths in an elastic half-space. However, such models ignore important subsurface structures, such as ring faults, that may influence the deformation pattern. Here we use both analog subsidence experiments and boundary element modeling to study the three-dimensional geometry and kinematics of caldera subsidence processes, evolving from an initial downsag to a later collapse stage. We propose that broad subsidence is mainly caused by volume decrease within a single magma reservoir, whereas buried ring-fault activity localizes the deformation within the caldera. Omitting ring faulting in physical models of subsiding calderas and using multiple point/sill-like sources instead can result in erroneous estimates of magma reservoir depths and volume changes. (C) 2019 Elsevier B.V. All rights reserved.
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.
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.