Refine
Has Fulltext
- no (237)
Year of publication
- 2018 (237) (remove)
Document Type
- Article (237) (remove)
Language
- English (237)
Is part of the Bibliography
- yes (237)
Keywords
- Himalaya (5)
- Germany (4)
- climate change (4)
- permafrost (4)
- Palynology (3)
- Quaternary (3)
- landslides (3)
- Arctic (2)
- Australia (2)
- Chile (2)
- Climate (2)
- Climate change (2)
- Earthquake hazards (2)
- Earthquake source observations (2)
- Microbial abundance (2)
- Northern Andes (2)
- Paleoclimatology (2)
- Paleogeography (2)
- Permafrost (2)
- Pollen (2)
- Seismicity and tectonics (2)
- South American Monsoon (2)
- Time-series analysis (2)
- Topography (2)
- Western Central Africa (2)
- geophysics (2)
- grain size (2)
- modeling (2)
- natural hazards (2)
- stable isotopes (2)
- vertical coupling (2)
- ABSH-system (1)
- Acadian (1)
- Accuracy Assessment (1)
- Acidithiobacillus (1)
- Acidobacteriaceae (1)
- Acidothermus (1)
- Adaptation (1)
- Advanced scatterometer (ASCAT) (1)
- Affordability (1)
- African Humid Period (1)
- Agrobacterium tumefaciens (1)
- Alkaline igneous rocks (1)
- Alpine Fault (1)
- Altai Mountains (1)
- Altyn Tagh Fault (1)
- Analog material (1)
- Anatexis (1)
- Anatolia westward motion (1)
- Andes (1)
- Antarctica (1)
- Aquifer thermal energy storage (ATES) (1)
- Ar-40/Ar-39 (1)
- Ar-40/Ar-39 dating (1)
- Arabia-Eurasia collision (1)
- Arbuscular mycorrhizal fungi (AMF) (1)
- Argentina (1)
- Argentine margine (1)
- Aridity gradient (1)
- Armenia (1)
- Asthenospheric fluid (1)
- Ataturk Dam Lake (1)
- Authigenic mineral formation (1)
- Authigenic mineral transformation (1)
- Batch experiments (1)
- Bayesian Network (1)
- Bayesian logistic regression (1)
- Benefits (1)
- Beringia (1)
- Bet-hedging germination (1)
- Bivalvia (1)
- Black shales (1)
- Body waves (1)
- Bolivia (1)
- Boreal ecosystems (1)
- Boundary line analysis (1)
- Brachypodium hybridum (1)
- Brahmaputra River (1)
- Brittle deformation (1)
- Bulfat (1)
- C and O K-edge spectra (1)
- CDOM (1)
- Canadian Beaufort Sea (1)
- Carbon isotope ratio (delta C-13(TOC)) (1)
- Catchment wetness (1)
- Central Asia (1)
- Central Europe (1)
- Central Pontides (1)
- Change detection (1)
- Change point analysis (1)
- Change points (1)
- China (1)
- Chota formation (1)
- Climate variability (1)
- Code_Aster (1)
- Complex network (1)
- Continuous temperature logging (1)
- Convergent margins (1)
- Cooling and unroofing history (1)
- Coordination failure (1)
- Coping appraisal (1)
- Coseismic rupture (1)
- Cosmic-ray neutrons (1)
- Cretaceous (1)
- Crustal stretching style (1)
- DIN EN 1998-1/NA (1)
- Dansgaard-Oeschger cycles (1)
- Death Valley (1)
- Debris avalanche deposit (1)
- Debris flow (1)
- Deformation (1)
- Density dependence (1)
- Denudation (1)
- Deposition (1)
- Digital Elevation Models (1)
- Direct economic loss (1)
- Disturbance index (1)
- Divergent mantle flow (1)
- Dora-Maira (1)
- Dormancy (1)
- Downholemethods (1)
- Drought indices (1)
- EXAFS (1)
- Ecosystem services (1)
- Ecuador (1)
- Eemian (1)
- Eger Rift (1)
- El Nino Southern Oscillation (ENSO) (1)
- Electrical resistivity tomography (ERT) (1)
- Electromagnetics (1)
- Elymus trachycaulus (1)
- Emergence time (1)
- Empirical site amplification functions (1)
- Enstatite (1)
- Euler-Liouville (1)
- Eurasian Dry Grassland Group (EDGG) (1)
- European Centre for Medium-Range Weather Forecasts (ECMWF) (1)
- European Vegetation Archive (EVA) (1)
- European basin system (1)
- European strong motion data (1)
- Event synchronization (1)
- Extension (1)
- Extreme discharge data (1)
- Extreme event (1)
- FGF23 (1)
- FTIR (1)
- Fault healing (1)
- Fault-rock microstructures (1)
- Final run (1)
- Firmicutes (1)
- Flash flood (1)
- Flash flood analysis (1)
- Flood (1)
- Flood regimes (1)
- Flood risk (1)
- Flood timing (1)
- Flowback (1)
- Fluid inclusions (1)
- Fluid-rock interaction (1)
- Forensic disaster analysis (1)
- Forest-savannah (1)
- Fourier analysis (1)
- Functions (1)
- Fuzzy classification (1)
- Fuzzy logic (1)
- GAIA (1)
- GNSS Reflectometry (1)
- GNSS forward scatterometry (1)
- GNSS reflectometry (1)
- Gehyra variegata (1)
- Gemmatimonadetes (1)
- Geochemistry (1)
- Geochronology (1)
- Geomorphology (1)
- Geophysical prospecting (1)
- Geothermal applications in Oman (1)
- Glacial geomorphology (1)
- Glacial lake (1)
- Glacier Mass Balances (1)
- Glacier foreland (1)
- Global precipitation measurement (1)
- Gondwana (1)
- Gongjue basin (1)
- GrassPlot (1)
- Gravel-bed (1)
- Gravitational potential energy (1)
- Gravity inversion (1)
- Green investment (1)
- Green surge (1)
- Gridded rainfall (1)
- Ground Motion Prediction Equation (1)
- Ground motion (1)
- Ground motion prediction equation (1)
- Ground squirrel nest (1)
- Hail (1)
- Halobacteria (1)
- Harran Plain (1)
- Hazard (1)
- He thermochronology (1)
- High Asia (1)
- High-Grade Metamorphism (1)
- Hindu Kush Himalayas (HKH) (1)
- Holocene (1)
- Hydrogeomorphic type (1)
- Hyetograph classification (1)
- IAPETUS suture (1)
- Illitization (1)
- Imaging (1)
- In situ Raman spectroscopy (1)
- Insurance (1)
- Integrated modelling (1)
- Inter-Andean valley (1)
- Intraplate earthquakes (1)
- Intraplate processes (1)
- Inverse Theory (1)
- Invertebrates (1)
- Ionosphere (1)
- Iran (1)
- Iraq (1)
- Island arc (1)
- JULIA (1)
- Japan subduction zone (1)
- Jurassic arc (1)
- Kalahari (1)
- Kanas Lake (1)
- Kaya Identity (1)
- Kazakhstan (1)
- Kettle holes (1)
- Kinematics of crustal (1)
- Klotho (1)
- LIA (1)
- LM (1)
- La Pampa (1)
- Lake Magadi (1)
- Lakes (1)
- Land and water (1)
- Land cover maps (1)
- Land use (1)
- Landform assemblages (1)
- Landsat time-series (1)
- Landslide (1)
- Landslides (1)
- Last cold stage (1)
- Late Pleistocene (1)
- Liguride Complex (1)
- Local structure Yttrium, Strontium, Lanthanum (1)
- Loess (1)
- Loss modelling (1)
- Lotus corniculatus (1)
- Lower crustal flow (1)
- MCA (1)
- MHT (1)
- MIT (1)
- MRD (1)
- MWAC (1)
- Magnetic+Susceptibility (1)
- Magnetostratigraphy (1)
- Markov Blanket (1)
- Mass transport (1)
- Measurement uncertainty (1)
- Mediterranean (1)
- Mediterranean Sea (1)
- Metapelites (1)
- Metasomatism (1)
- Method comparison (1)
- Mineralization (1)
- Mixed-effects regression (1)
- Model coupling (1)
- Model-data comparison (1)
- Modeling (1)
- Modern analogues (1)
- Monsoon (1)
- Mt. Quincan (1)
- Multifunctionality (1)
- Multiple regression (1)
- NAC transcription factor ATAF1 (1)
- NOR (1)
- NW Argentina (1)
- Namche Barwa (1)
- Namche Barwa syntaxis (1)
- Nanogranitoids (1)
- Natural hazards (1)
- Near+Surface (1)
- Neoacadian orogenic cycles (1)
- Neotethys (1)
- Neotropics (1)
- Nepal (1)
- Nevado Coropuna (1)
- Niger River Basin (1)
- Nordic catchments (1)
- North America (1)
- North Atlantic forcing (1)
- North Calabrian Unit (1)
- Northeast of Brazil (1)
- Northern Asia (1)
- Northern Tibetan Plateau (1)
- Northern Zagros Suture Zone (NZSZ) (1)
- Numerical modelling (1)
- Ocean tidal model (1)
- Oman (1)
- Ophiolite (1)
- Optical methods (1)
- Orai1 (1)
- Organic geochemistry (1)
- PT pseudosection (1)
- Pacific-Asia convergence (1)
- Palaeo-vegetation (1)
- Palaeoenvironmental reconstruction (1)
- Palaeotropics (1)
- Paleobiogeography (1)
- Paleoclimate (1)
- Paleoseismology (1)
- Pamir (1)
- Paratethys (1)
- Pargasite (1)
- Paris-Edinburgh press (1)
- Parlung River (1)
- Participation (1)
- Patagonian Ice Sheet (1)
- Pebas/Acre system (1)
- Peclet number (1)
- Pedestal moraine (1)
- Pedogenic carbonate nodules (1)
- Performance (1)
- Peri-urban agriculture (1)
- Petrography (1)
- Physical rock properties (1)
- Plant macrofossils (1)
- Plate tectonics (1)
- Pollen trap (1)
- Ponds (1)
- Potassium (1)
- Principal component analysis (1)
- Provenance (1)
- Pseudoleucite (1)
- Public policy (1)
- QBO (1)
- Quaternary climate change (1)
- Quaternary maar volcanism (1)
- Quaternary terraces (1)
- Radar rainfall data (1)
- Radioactivity (1)
- Radiogenic isotopes (1)
- Rainfall event (1)
- Rainfall network (1)
- Rainfall-triggered landslide (1)
- Raman spectroscopy (1)
- Random Forest (1)
- Rb-Sr mineral isochron (1)
- Recurrence plot (1)
- Reflectance spectroscopy (1)
- Relative dating (1)
- Remote sensing (1)
- Reservoir network (1)
- Residual analysis (1)
- Ring complex (1)
- Rio Grande (1)
- Risk dominance (1)
- Riskspreading (1)
- River (1)
- River flow (1)
- Rock-Eval pyrolysis (1)
- Rodinia (1)
- Romania (1)
- Root morphology (1)
- SEM (1)
- SSW (1)
- SWAT (1)
- Salinic (1)
- Savanna (1)
- Sediment Flux (1)
- Sediment load (1)
- Sediment retention (1)
- Sediment tails (1)
- Sedimentary basin (1)
- Seed competition (1)
- Seismic attenuation (1)
- Seismic hazard (1)
- Seismic interferometry (1)
- Seismic load parameters (1)
- Seismic noise (1)
- Seismicity (1)
- Seismicity and tectonic (1)
- Semiarid catchment (1)
- Sensitivity analysis (1)
- Sentinel-1 (1)
- Shear zones (1)
- Shore vegetationa (1)
- Siberia (1)
- Siderite (1)
- Silicate- carbonate melts/glasses (1)
- Single station sigma (1)
- Site classification (1)
- Site effects (1)
- Site response (1)
- Siwaliks (1)
- Slope exposure (1)
- Soil C- and N stocks (1)
- Soil microorganisms (1)
- Soil moisture (1)
- Soil process modelling (1)
- Soil respiration (1)
- Soil texture (1)
- Soil water balance (1)
- Soil water variability (1)
- South America (1)
- Southern Apennines (1)
- Southern Ethiopian Rift (1)
- Spatial coherence (1)
- Spectral clustering analysis (1)
- Sporopollenin (1)
- Sr-Nd-Pb isotopes (1)
- Stag hunt (1)
- Stakeholder (1)
- Statistical Methods (1)
- Steppe-desert (1)
- Storm (1)
- Strategic uncertainty (1)
- Structure-from-motion (1)
- Subduction zone processes (1)
- Supercontinent breakup (1)
- Supercontinent cycle (1)
- Supra-subduction zone (1)
- Surface waves and free oscillations (1)
- Sustainable cities (1)
- Synthetic pyroxene (1)
- TDS-1 (1)
- TL of feldspar (1)
- TL-thermochronology (1)
- TLS (1)
- Tasseled Cap transformation (1)
- TechDemoSat-1 (TDS-1) (1)
- Tectonic geomorphology (1)
- Temperature (1)
- Temporal spread (1)
- The northern margin of the South China Sea (1)
- Thermochronology (1)
- Tian Shan piedmont (1)
- Tian shan (1)
- Tibet (1)
- Tillage erosion (1)
- Total organic carbon (1)
- Tracer (1)
- Tropical high mountains (1)
- Tropical lake (1)
- U/Pb dating of titanite (1)
- UAS (1)
- Uncertainty (1)
- Uncertainty analysis (1)
- Unconventional gas production (1)
- Urban CO2 emissions (1)
- Urban Kaya relation (1)
- Urban learning lab (1)
- Urban scaling (1)
- VIIRS DNB (1)
- VLBI analysis (1)
- Van Krevelen diagram (1)
- Variability (1)
- Vegetation-climate disequilibrium (1)
- Vertisol genesis (1)
- Volcanic collapse (1)
- Volume reconstruction (1)
- Water management (1)
- Water quality (1)
- Water storage dynamic (1)
- Wave propagation (1)
- Weathering (1)
- Western Africa (1)
- White mica Ar-40/ Ar-39 spot ages (1)
- Wind erosion (1)
- Yamal (1)
- Younger Dryas (1)
- Yucatan Peninsula (1)
- Yukon (1)
- Z-P approach (1)
- Zeolites (1)
- acidophilic microorganisms (1)
- active continental margin (1)
- age class forest (1)
- age-elevation relationships (1)
- alluvial fans (1)
- analysis (1)
- and Granite Magmatism (1)
- arctic tundra (1)
- artefact avoidance (1)
- artificial light at night (1)
- assimilation (1)
- attribution (1)
- beaver (1)
- behavioral adaptation (1)
- biodiversity (1)
- biogeochemistry (1)
- body condition (1)
- body growth rate (1)
- calcium (1)
- carbon dioxide (1)
- central Andean foreland (1)
- charcoal (1)
- chemical weathering (1)
- chemical-mechanical interaction (1)
- climate (1)
- climatology (1)
- coastal erosion (1)
- coexistence (1)
- coherence (1)
- competition-defense trade-off (1)
- composite properties (1)
- concentration (1)
- core flow (1)
- cosmogenic nuclides (1)
- cross-calibration (1)
- crustal density structure (1)
- data archeology (1)
- data scarcity (1)
- decompression (1)
- deep biosphere (1)
- defense against predation (1)
- depletion amplitude (1)
- detrital geochronology and thermochronology (1)
- diatom indicator species (1)
- differential split-sample test (1)
- digital rock physics (1)
- dissolved silica concentration (1)
- downscaling (1)
- drought stress (1)
- dumortierite (1)
- duplex (1)
- early Holocene (1)
- earthquake (1)
- earthquakes (1)
- eastern Qiangtang terrane (1)
- elastic geobarometry (1)
- elastic properties (1)
- electromagnetic (1)
- electromagnetic scattering (1)
- elevated CO2 (1)
- equatorial electrojet (1)
- equatorial electrojet (EEJ) (1)
- equatorial ionosphere (1)
- equatorial plasma depletion (1)
- equifinality (1)
- erosion-climate link (1)
- evapotranspiration (1)
- exhumation (1)
- fission track thermochronology (1)
- flood generating processes (1)
- flood loss (1)
- flood seasonality (1)
- floods (1)
- fluid flow and thermal modelling (1)
- fluids (1)
- fluvial geomorphology (1)
- fluvial response (1)
- fracking (1)
- functional response (1)
- gabbros (1)
- geo–bio interaction (1)
- glacial buzzsaw (1)
- glacial equilibrium (1)
- glacial erosion (1)
- glacial geomorphology (1)
- glacier chemistry (1)
- glacierized basins (1)
- grassland vegetation (1)
- gravity modeling (1)
- ground ice (1)
- ground motion (1)
- ground motion prediction equations (1)
- hillslopes (1)
- hominins (1)
- human activity (1)
- human settlement (1)
- hydrograph partitioning curves (1)
- hydrologic modeling (1)
- hydrological modelling (1)
- hydrology (1)
- hydropower (1)
- hyperaridity (1)
- hypothetical trend (1)
- hysteresis (1)
- inclination shallowing (1)
- inclusions (1)
- indirect facilitation (1)
- induced seismicity (1)
- industrial ecology (1)
- input-output (I-O) analysis (1)
- interannual variability (1)
- inversion tectonics (1)
- ionospheric irregularities (1)
- kokchetavite (1)
- kumdykolite (1)
- lake (1)
- lake catchments (1)
- lakes (1)
- land-use intensity (1)
- landscape evolution (1)
- late Holocene (1)
- lateral fluxes (1)
- length of day (1)
- light pollution (1)
- lithium-ion battery electrolyte (1)
- lithology (1)
- loading conditions (1)
- localization (1)
- lower crustal bodies (1)
- luminescence dating (1)
- lunar tide (1)
- macroecology (1)
- magnesiodumortierite (1)
- magnetic surveying (1)
- mantle deformation (1)
- marble (1)
- maximum magnitude (1)
- meadow (1)
- micro-CT (1)
- microbes (1)
- microbiology (1)
- mineral deposits (1)
- minirhizotron (1)
- model calibration (1)
- model enhancement (1)
- moraine (1)
- multi-taxon (1)
- multi-thermochronometer (1)
- nanogranitoids (1)
- natural dams (1)
- nested plot (1)
- night lights (1)
- non-resonant inelastic X-ray scattering (1)
- north-atlantic (1)
- numerical simulation (1)
- ocean surface (1)
- oceanic transform and fracture zone processes (1)
- ophiolite (1)
- optical cell (1)
- organic carbon (1)
- pH (1)
- palaeoclimate records (1)
- paleo-sediment (1)
- paleoaltimetry (1)
- paleoclimate (1)
- paleohydrology (1)
- paleolimnology (1)
- paleomagnetism (1)
- parameter emulation (1)
- parathyroid hormone (1)
- partial melting (1)
- pasture (1)
- phase space time delay embedding reconstruction (1)
- phosphate (1)
- picoseismicity (1)
- planetary waves (1)
- plant transformation (1)
- plasma drifts (1)
- polymorphs (1)
- population dynamics (1)
- population recovery (1)
- population size (1)
- postmidnight (1)
- potentially affected areas (PAA) (1)
- pre-Maya settlement (1)
- precipitation reconstruction (1)
- predator-prey cycles (1)
- prediction (1)
- prediction uncertainty (1)
- process life cycle assessment (1)
- processes and landforms of glacial erosion (1)
- production rate (1)
- protection motivation theory (PMT) (1)
- proxy calibration (1)
- radio science (1)
- radio wave propagation (1)
- radiogenic heat-production (1)
- rain detection (1)
- rain splash (1)
- rainforest crisis (1)
- rapid damage mapping (1)
- reanalysis (1)
- redox reactions (1)
- relief development (1)
- remote sensing data (1)
- risk communication (1)
- river incision (1)
- rock exhumation (1)
- rock glaciers (1)
- rock uplift-relief scaling (1)
- root biomass (1)
- root density (1)
- root lengths density (1)
- root penetration (1)
- rotation (1)
- runoff (1)
- sPlot (1)
- saline pan (1)
- salmon (1)
- scale-dependence (1)
- scaling relation (1)
- sea level (1)
- sea level rise (1)
- sediment (1)
- sediment flux (1)
- sediment supply (1)
- sediment transport (1)
- sediment yield (1)
- sedimentary basins (1)
- seismomechanics (1)
- self-similarity (1)
- sensitivity analyses (1)
- service sectors (1)
- shaking (1)
- shrub expansion (1)
- snow melt (1)
- social vulnerability (1)
- soil erosion (1)
- soil horizonation (1)
- source parameters (1)
- species-area relationship (SAR) (1)
- speleothem (1)
- speleothems (1)
- statistical methods (1)
- steady state topography (1)
- stream (1)
- stream water (1)
- streamflow (1)
- subduction interface (1)
- submerged cave (1)
- subsidence analysis (1)
- subsurface heterogeneities (1)
- sulfur speciation (1)
- survival (1)
- swarm LP (1)
- synthesis (1)
- system boundary (1)
- table-top experiment (1)
- taiga forest (1)
- tectonic evolution (1)
- tectonic shortening (1)
- tectonics (1)
- tetrahedral boron (1)
- thermal field (1)
- thermal modeling (1)
- thermodynamic modeling (1)
- thermokarst (1)
- thermoregulation (1)
- tide gauge (1)
- tomato (1)
- torsion (1)
- trace elements (1)
- trench advance (1)
- trend analysis (1)
- truncation error estimate (1)
- ultrahigh-pressure (UHP) (1)
- ultramafic granulites (1)
- ungauged basins (1)
- unmanaged forest (1)
- uplift and erosion (1)
- vegetation dynamics (1)
- vegetation-plot database (1)
- volcanism (1)
- voring basin (1)
- vulnerability (1)
- water isotopes (1)
- weather pattern (1)
- wind speed (1)
Institute
- Institut für Geowissenschaften (237) (remove)
High-latitude treeless ecosystems represent spatially highly heterogeneous landscapes with small net carbon fluxes and a short growing season. Reliable observations and process understanding are critical for projections of the carbon balance of the climate-sensitive tundra. Space-borne remote sensing is the only tool to obtain spatially continuous and temporally resolved information on vegetation greenness and activity in remote circumpolar areas. However, confounding effects from persistent clouds, low sun elevation angles, numerous lakes, widespread surface inundation, and the sparseness of the vegetation render it highly challenging. Here, we conduct an extensive analysis of the timing of peak vegetation productivity as shown by satellite observations of complementary indicators of plant greenness and photosynthesis. We choose to focus on productivity during the peak of the growing season, as it importantly affects the total annual carbon uptake. The suite of indicators are as follows: (1) MODIS-based vegetation indices (VIs) as proxies for the fraction of incident photosynthetically active radiation (PAR) that is absorbed (fPAR), (2) VIs combined with estimates of PAR as a proxy of the total absorbed radiation (APAR), (3) sun-induced chlorophyll fluorescence (SIF) serving as a proxy for photosynthesis, (4) vegetation optical depth (VOD), indicative of total water content and (5) empirically upscaled modelled gross primary productivity (GPP). Averaged over the pan-Arctic we find a clear order of the annual peak as APAR <= GPP < SIF < VIs/VOD. SIF as an indicator of photosynthesis is maximised around the time of highest annual temperatures. The modelled GPP peaks at a similar time to APAR. The time lag of the annual peak between APAR and instantaneous SIF fluxes indicates that the SIF data do contain information on light-use efficiency of tundra vegetation, but further detailed studies are necessary to verify this. Delayed peak greenness compared to peak photosynthesis is consistently found across years and land-cover classes. A particularly late peak of the normalised difference vegetation index (NDVI) in regions with very small seasonality in greenness and a high amount of lakes probably originates from artefacts. Given the very short growing season in circumpolar areas, the average time difference in maximum annual photosynthetic activity and greenness or growth of 3 to 25 days (depending on the data sets chosen) is important and needs to be considered when using satellite observations as drivers in vegetation models.
The Gutenberg-Richter relation for earthquake magnitudes is the most famous empirical law in seismology. It states that the frequency of earthquake magnitudes follows an exponential distribution; this has been found to be a robust feature of seismicity above the completeness magnitude, and it is independent of whether global, regional, or local seismicity is analyzed. However, the exponent b of the distribution varies significantly in space and time, which is important for process understanding and seismic hazard assessment; this is particularly true because of the fact that the Gutenberg-Richter b-value acts as a proxy for the stress state and quantifies the ratio of large-to-small earthquakes. In our work, we focus on the automatic detection of statistically significant temporal changes of the b-value in seismicity data. In our approach, we use Bayes factors for model selection and estimate multiple change-points of the frequency-magnitude distribution in time. The method is first applied to synthetic data, showing its capability to detect change-points as function of the size of the sample and the b-value contrast. Finally, we apply this approach to examples of observational data sets for which b-value changes have previously been stated. Our analysis of foreshock and after-shock sequences related to mainshocks, as well as earthquake swarms, shows that only a portion of the b-value changes is statistically significant.
Tundra be dammed
(2018)
Increasing air temperatures are changing the arctic tundra biome. Permafrost is thawing, snow duration is decreasing, shrub vegetation is proliferating, and boreal wildlife is encroaching. Here we present evidence of the recent range expansion of North American beaver (Castor canadensis) into the Arctic, and consider how this ecosystem engineer might reshape the landscape, biodiversity, and ecosystem processes. We developed a remote sensing approach that maps formation and disappearance of ponds associated with beaver activity. Since 1999, 56 new beaver pond complexes were identified, indicating that beavers are colonizing a predominantly tundra region (18,293km(2)) of northwest Alaska. It is unclear how improved tundra stream habitat, population rebound following overtrapping for furs, or other factors are contributing to beaver range expansion. We discuss rates and likely routes of tundra beaver colonization, as well as effects on permafrost, stream ice regimes, and freshwater and riparian habitat. Beaver ponds and associated hydrologic changes are thawing permafrost. Pond formation increases winter water temperatures in the pond and downstream, likely creating new and more varied aquatic habitat, but specific biological implications are unknown. Beavers create dynamic wetlands and are agents of disturbance that may enhance ecosystem responses to warming in the Arctic.
Bedrock landsliding in mountain belts can elevate overall chemical weathering rates through rapid dissolution of exhumed reactive mineral phases in transiently stored deposits. This link between a key process of erosion and the resultant weathering affects the sequestering of carbon dioxide through weathering of silicate minerals and broader links between erosion in active orogens and climate change. Here we address the effect on the carbon cycle of weathering induced by bedrock landsliding in Taiwan and the Western Southern Alps of New Zealand. Using solute chemistry data from samples of seepage from landslide deposits and river discharge from catchments with variable proportions of landsliding, we model the proportion of silicate and carbonate weathering and the balance of sulfuric and carbonic acids that act as weathering agents. We correct for secondary precipitation, geothermal, and cyclic input, to find a closer approximation of the weathering explicitly occurring within landslide deposits. We find highly variable proportions of sulfuric and carbonic acids driving weathering in landslides and stable hillslopes. Despite this variability, the predominance of rapid carbonate weathering within landslides and catchments where mass wasting is prevalent results at best in limited sequestration of carbon dioxide by this process of rapid erosion. In many cases where sulfuric acid is a key weathering agent, a net release of CO2 to the atmosphere occurs. This suggests that a causal link between erosion in mountain belts and climate change through the sequestration of CO2, if it exists, must operate through a process other than chemical weathering driven by landsliding. Plain Language Summary There is a long-standing debate surrounding the link between erosion and climate. It is often suggested that as temperatures increase, rainier and stormier weather could increase erosion of rock; as that rock is exposed, silicate minerals within could break down, which, on long time scales, can remove CO2 from the atmosphere, lowering global temperatures and acting as a negative feedback. Recent studies have shown that landslide deposits are key locations for the link between chemical weathering and physical erosion in some mountain belts. To test how landslides affect the erosion-climate link, we used samples of water seeping through landslides in Taiwan and New Zealand to calculate the amount of carbon dioxide that is either absorbed or released through this chemical reaction. We find that the large amount of freshly exposed rock in Taiwanese landslide deposits contains significant carbonate rock and sulfide minerals; the net result of the weathering of these minerals is a release of carbon dioxide, which inverts the traditional perspective on the role erosion plays in controlling carbon dioxide release. In some mountain belts, it seems that increased erosion and resulting landsliding may act to increase carbon dioxide in the air, opening further questions into the nature of erosional-climatic links.
Many regions central to our understanding of tectonics and landscape evolution are active or ancient magmatic terranes, and robust interpretation of low-temperature thermochronologic ages in these settings requires careful attention to the drivers of rock heating and cooling, including magmatism. However, we currently lack a quantitative framework for evaluating the potential role of magmatic coolingthat is, post-magmatic thermal relaxationin shaping cooling age patterns in regions with a history of intrusive magmatism. Here we use analytical approximations and numerical models to characterize how low-temperature thermochronometers document cooling inside and around plutons in steadily exhuming environments. Our models predict that the thermal field a pluton intrudes into, specifically the ambient temperatures relative to the closure temperature of a given thermochronometer, is as important as the pluton size and temperature in controlling the pattern and extent of thermochronometer resetting in the country rocks around a pluton. We identify one advective and several conductive timescales that govern the relationship between the crystallization and cooling ages inside a pluton. In synthetic vertical age-elevation relationships (AERs), resetting next to plutons results in changes in AER slope that could be misinterpreted as past changes in exhumation rate if the history of magmatism is not accounted for. Finally, we find that large midcrustal plutons, such as those emplaced at similar to 10-15-km depth, can reset the low-temperature thermochronometers far above them in the upper crusta result with considerable consequences for thermochronology in arcs and regions with a history of magmatic activity that may not have a surface expression.
The variation of Rayleigh ellipticity versus frequency is gaining popularity in site characterization. It becomes a necessary observable to complement dispersion curves when inverting shear wave velocity profiles. Various methods have been proposed so far to extract polarization from ambient vibrations recorded on a single three-component station or with an array of three-component sensors. If only absolute values were recovered 10 yr ago, new array-based techniques were recently proposed with enhanced efficiencies providing also the ellipticity sign. With array processing, higher-order modes are often detected even in the ellipticity domain. We suggest to explore the properties of a high-resolution beamforming where radial and vertical components are explicitly included. If N is the number of three-component sensors, 2N x 2N cross-spectral density matrices are calculated for all presumed directions of propagation. They are built with N radial and N vertical channels. As a first approach, steering vectors are designed to fit with Rayleigh wave properties: the phase shift between radial and vertical components is either -Pi/2 or Pi/2. We show that neglecting the ellipticity tilt due to attenuation has only minor effects on the results. Additionally, we prove analytically that it is possible to retrieve the ellipticity value from the usual maximization of the high-resolution beam power. The method is tested on synthetic data sets and on experimental data. Both are reference sites already analysed by several authors. A detailed comparison with previous results on these cases is provided.
Near-surface geophysical imaging of alluvial fan settings is a challenging task but crucial for understating geological processes in such settings. The alluvial fan of Ghor Al-Haditha at the southeast shore of the Dead Sea is strongly affected by localized subsidence and destructive sinkhole collapses, with a significantly increasing sinkhole formation rate since ca. 1983. A similar increase is observed also on the western shore of the Dead Sea, in correlation with an ongoing decline in the Dead Sea level. Since different structural models of the upper 50 m of the alluvial fan and varying hypothetical sinkhole processes have been suggested for the Ghor Al-Haditha area in the past, this study aimed to clarify the subsurface characteristics responsible for sinkhole development. For this purpose, high-frequency shear wave reflection vibratory seismic surveys were carried out in the Ghor Al-Haditha area along several crossing and parallel profiles with a total length of 1.8 and 2.1 km in 2013 and 2014, respectively. The sedimentary architecture of the alluvial fan at Ghor Al-Haditha is resolved down to a depth of nearly 200 m at a high resolution and is calibrated with the stratigraphic profiles of two boreholes located inside the survey area. The most surprising result of the survey is the absence of evidence of a thick (>2-10 m) compacted salt layer formerly suggested to lie at ca. 35-40 m depth. Instead, seismic reflection amplitudes and velocities image with good continuity a complex interlocking of alluvial fan deposits and lacustrine sediments of the Dead Sea between 0 and 200 m depth. Furthermore, the underground section of areas affected by sinkholes is characterized by highly scattering wave fields and reduced seismic interval velocities. We propose that the Dead Sea mud layers, which comprise distributed inclusions or lenses of evaporitic chloride, sulfate, and carbonate minerals as well as clay silicates, become increasingly exposed to unsaturated water as the sea level declines and are consequently destabilized and mobilized by both dissolution and physical erosion in the subsurface. This new interpretation of the underlying cause of sinkhole development is supported by surface observations in nearby channel systems. Overall, this study shows that shear wave seismic reflection technique is a promising method for enhanced near-surface imaging in such challenging alluvial fan settings.
Woodlands and steppes
(2018)
Based on fossil organism remains including plant macrofossils, charcoal, pollen, and invertebrates preserved in syngenetic deposits of the Batagay permafrost sequence in the Siberian Yana Highlands, we reconstructed the environmental history during marine isotope stages (MIS) 6 to 2. Two fossil assemblages, exceptionally rich in plant remains, allowed for a detailed description of the palaeo-vegetation during two climate extremes of the Late Pleistocene, the onset of the last glacial maximum (LGM) and the last interglacial. In addition, altogether 41 assemblages were used to outline the vegetation history since the penultimate cold stage of MIS 6. Accordingly, meadow steppes analogue to modern communities of the phytosociological order Festucetalia lenensis formed the primary vegetation during the Saalian and Weichselian cold stages. Cold-resistant tundra-steppe communities (Carici rupestris-Kobresietea bellardii) as they occur above the treeline today were, in contrast to more northern locations, mostly lacking. During the last interglacial, open coniferous woodland similar to modern larch taiga was the primary vegetation at the site. Abundant charcoal indicates wildfire events during the last interglacial. Zoogenic disturbances of the local vegetation were indicated by the presence of ruderal plants, especially by abundant Urtica dioica, suggesting that the area was an interglacial refugium for large herbivores. Meadow steppes, which formed the primary vegetation during cold stages and provided potentially suitable pastures for herbivores, were a significant constituent of the plant cover in the Yana Highlands also under the full warm stage conditions of the last interglacial. Consequently, meadow steppes occurred in the Yana Highlands during the entire investigated timespan from MIS 6 to MIS 2 documenting a remarkable environmental stability. Thus, the proportion of meadow steppe vegetation merely shifted in response to the respectively prevailing climatic conditions. Their persistence indicates low precipitation and a relatively warm growing season throughout and beyond the late Pleistocene. The studied fossil record also proves that modern steppe occurrences in the Yana Highlands did not establish as late as in the Holocene but instead are relicts of a formerly continuous steppe belt extending from Central Siberia to Northeast Yakutia during the Pleistocene. The persistence of plants and invertebrates characteristic of meadow steppe vegetation in interior Yakutia throughout the late Quaternary indicates climatic continuity and documents the suitability of this region as a refugium also for other organisms of the Pleistocene mammoth steppe including the iconic large herbivores. (C)2018 Elsevier Ltd. All rights reserved.
Both enstatite spectra are very bright in the VIS and NIR and show almost neutral to slightly bluish spectral slopes with a steep absorption in the UV. Very low iron in the enstatites (below similar to 0.04 wt% FeO) already results in weak albeit noticeable absorptions in the VNIR between 0.4 and 0.9 mu m. Orthoenstatite and clinoenstatite are not distinguishable based only on their spectra in the VIS and NIR. At the Reststrahlen bands in the MIR a systematic difference in the number and exact position of local minima at similar to 10 mu m between clinoenstatite and orthoenstatite is evident. This can be used to discern between the polymorphs in this wavelength range. Additionally, we can distinguish between Fe-free low- and high-Ca pyroxenes in the MIR.
We present the first high-resolution trace element (Mg/Ca, Sr/Ca, Ba/Ca) record from a stalagmite in southwestern Romania covering the last 3.6 ka, which provides the potential for quantitative climate reconstruction. Precise age control is based on three independent dating methods, in particular for the last 250 yr, where chemical lamina counting is combined with the identification of the 20th century radiocarbon bomb peak and Th-230/U dating. Long-term cave monitoring and model simulations of drip water and speleothem elemental variability indicate that precipitation-related processes are the main drivers of speleothem Mg/Ca ratios. Calibration against instrumental climate data shows a significant anti-correlation of speleothem Mg/Ca ratios with autumn/winter (October to March) precipitation (r = -0.61, p < 0.01), which is statistically robust when considering age uncertainties and auto-correlation. This relationship is used to develop a quantitative reconstruction of autumn/winter precipitation. During the late Holocene, our data suggest a heterogeneous pattern of past regional winter hydroclimate in the Carpathian/Balkan realm, along with intermittent weakening of the dominant influence of North Atlantic forcing. In agreement with other regional paleo-hydrological reconstructions, the observed variability reveals periodically occurring strong NW-SE hydro-climate gradients. We hypothesize, that this pattern is caused by shifts of the eastern edge of the area of influence of the NAO across central eastern Europe due to the confluence of North Atlantic forcing, and other climatic features such as the East Atlantic/Western Russia (EAWR) pattern. (C) 2018 Elsevier B.V. All rights reserved.
The basic seismic load parameters for the upcoming national design regulation for DIN EN 1998-1/NA result from the reassessment of the seismic hazard supported by the German Institution for Civil Engineering (DIBt). This 2016 version of the national seismic hazard assessment for Germany is based on a comprehensive involvement of all accessible uncertainties in models and parameters and includes the provision of a rational framework for integrating ranges of epistemic uncertainties and aleatory variabilities in a comprehensive and transparent way. The developed seismic hazard model incorporates significant improvements over previous versions. It is based on updated and extended databases, it includes robust methods to evolve sets of models representing epistemic uncertainties, and a selection of the latest generation of ground motion prediction equations. The new earthquake model is presented here, which consists of a logic tree with 4040 end branches and essential innovations employed for a realistic approach. The output specifications were designed according to the user oriented needs as suggested by two review teams supervising the entire project. Seismic load parameters, for rock conditions of nu(S30) = 800 m/s, are calculated for three hazard levels (10, 5 and 2% probability of occurrence or exceedance within 50 years) and delivered in the form of uniform hazard spectra, within the spectral period range 0.02-3 s, and seismic hazard maps for peak ground acceleration, spectral response accelerations and for macroseismic intensities. Results are supplied as the mean, the median and the 84th percentile. A broad analysis of resulting uncertainties of calculated seismic load parameters is included. The stability of the hazard maps with respect to previous versions and the cross-border comparison is emphasized.
Nearly 13,000 years ago, the warming trend into the Holocene was sharply interrupted by a reversal to near glacial conditions. Climatic causes and ecological consequences of the Younger Dryas (YD) have been extensively studied, however proxy archives from the Mediterranean basin capturing this period are scarce and do not provide annual resolution. Here, we report a hydroclimatic reconstruction from stable isotopes (delta O-18, delta C-13) in subfossil pines from southern France. Growing before and during the transition period into the YD (12 900-12 600 cal BP), the trees provide an annually resolved, continuous sequence of atmospheric change. Isotopic signature of tree sourcewater (delta O-18(sw)) and estimates of relative air humidity were reconstructed as a proxy for variations in air mass origin and precipitation regime. We find a distinct increase in inter-annual variability of sourcewater isotopes (delta O-18(sw)), with three major downturn phases of increasing magnitude beginning at 12 740 cal BP. The observed variation most likely results from an amplified intensity of North Atlantic (low delta O-18(sw)) versus Mediterranean (high delta O-18(sw)) precipitation. This marked pattern of climate variability is not seen in records from higher latitudes and is likely a consequence of atmospheric circulation oscillations at the margin of the southward moving polar front.
The Sierra Nevada de Santa Marta in NW Colombia is an isolated massif at the northernmost end of the Andes chain near the boundary with the Caribbean plate. Previous geologic mapping and K-Ar dating have shown that Jurassic plutonic and volcanic units make up a large part of the Santa Marta Massif (SMM). These rocks have been considered to be part of a Jurassic magmatic arc extending from NW Colombia to northern Chile, but without any geochemical basis for comparison. This paper reports on a geochemical and Sr-Nd-Pb isotope study of the Jurassic rocks in the SMM and provides 12 new U-Pb zircon ages from in-situ laser ICP-MS dating. The plutonic and volcanic units span a range from 45 to 78 wt.% SiO2, with a dominance of intermediate to felsic compositions with SiO2 > 57 wt.%. They classify as calc-alkaline, medium to high-K, metaluminous rocks with trace-element features typical for arc-derived magma series. In terms of their major and trace-element compositions, the SMM Jurassic units overlap with contemporary plutonic and volcanic rocks from other regions of the Central and Eastern Cordilleras of Colombia, and confirm an arc affinity. The new U-Pb ages range from 176 +/- 1 Ma to 192 +/- 2 Ma (n = 12), with most between 180 and 188 Ma (n = 7). The initial Sr isotope ratios (at 180 Ma) are between 0.7012 and 0.7071 (n = 29), with 3 outliers attributed to mobilization of Rb and/or Sr, Nearly all samples have negative( )epsilon Nd-(180) values between - 10.3 and 0.0 (n = 30), the two exceptions being only slightly positive (1.1 and 1.9). Measured Pb isotope ratios fall in a narrow range, with Pb-206/Pb-204 from 18.02 to 19.95, (207) Pb/(204) Pb from 15.56 to 15.67 and Pb-208/Pb-204 from 37.76 to 39.04 (n = 28). In the regional context of previous studies, these results confirm early Jurassic ages and an arc affinity for the widespread magmatism exposed in the eastern and northeastern Colombian Andes. We also note patterns in the distribution and composition of magmas. The magmatic activity in the Central Cordillera tends to be younger than in the Eastern Cordillera and is spatially more restricted to the vicinity of regional fault systems. In terms of composition, Jurassic igneous rocks in the Eastern Cordillera have systematically lower epsilon Nd-(180) values than those from the Central Cordillera, whereas the Pb isotope ratios overlap. We ascribe the Nd isotope variations to heterogeneity in the mantle source and/or degree of crustal contamination, whereas the Pb isotope ratios are crust-dominated and similar throughout the region. The spatio-temporal and compositional evolution of Jurassic magmatic rocks in the Northern Andes reflect the major plate kinematic readjustment between the Triassic and the Early Jurassic in the proto-Andean margin.
The article provides a contribution to the glacial geomorphology of tropical high-altitude mountains. It focuses on the diversity of glacial landforms in the Cordillera Blanca (Peru). The landforms are classified as individual landform types as well as in their hypsometric sequence. Their spatial arrangement is depicted in vertical series and summarised in a glacial-morphological altitudinal belt model. The glacial landforms are analysed in a genetic succession from the moraine-delimited glacier forefields to the excessive pedestal moraines and taking account of the topographical conditions in the verticality. The Pleistocene glacial relief formation is considered as one of the dominating factors in the development of the moraine types. The sediments are on the one hand in situ deposition forms and on the other depositions dislocated by mass movements and represent paraglacial landforms. Using the example of the Paron Valley in the Artesonraju-Huandoy Massif and neighbouring valleys, the article provides a typology of the range of glacial and paraglacial landforms based on a holistic glacial landsystem approach with special reference to pedestal moraines. (C) 2018 Elsevier B.V. All rights reserved.
Risk-based insurance is a commonly proposed and discussed flood risk adaptation mechanism in policy debates across the world such as in the United Kingdom and the United States of America. However, both risk-based premiums and growing risk pose increasing difficulties for insurance to remain affordable. An empirical concept of affordability is required as the affordability of adaption strategies is an important concern for policymakers, yet such a concept is not often examined. Therefore, a robust metric with a commonly acceptable affordability threshold is required. A robust metric allows for a previously normative concept to be quantified in monetary terms, and in this way, the metric is rendered more suitable for integration into public policy debates. This paper investigates the degree to which risk-based flood insurance premiums are unaffordable in Europe. In addition, this paper compares the outcomes generated by three different definitions of unaffordability in order to investigate the most robust definition. In doing so, the residual income definition was found to be the least sensitive to changes in the threshold. While this paper focuses on Europe, the selected definition can be employed elsewhere in the world and across adaption measures in order to develop a common metric for indicating the potential unaffordability problem.
Tomato (Solanum lycopersicum L.) being an important vegetable is cultivated and used throughout the world. It not only contributes in fulfilling the basic nutritional requirements of the human body but also has many health benefits due to its rich biochemical composition. However, its production at large scale is hampered by many limiting factors such as biotic and abiotic stresses. Among the different abiotic stresses, drought poses drastic impact on tomato yield. Drought stress is genetically regulated by many transcription factors that not only regulate the stress responsive mechanism but also facilitate the growth and development of tomato plants. NAC is an important stress related transcription factor genes family, and the ATAF1 gene, a member of this family, is involved in ABA signaling and stress response. In this study, tomato variety Rio Drande was transformed with drought resistant ATAF1 gene via Agrobacterium mediated gene transformation method. The ATAF1 gene was first cloned in the pK7WFG2 vector having kanamycin selectable marker and then it was introduced in the Agrobacterium tumefaciens strain GV3101 through heat shock method. The tomato cotyledon and hypocotyl ex-plants of variety "Rio Ggrande" were cultured on callus induction medium (MS + 2.5 mg/L IAA + 2 mg/L BAP). The calli were then infected with Agrobacterium tumefaciens strain GV3101 containing ATAF1 gene and selection was carried out on the kanamycin selectable medium (MS + 100 mg/L Kan), and were regenerated on MS medium with 1 mg/L IAA + 1 mg/L BAP. Out of 216 putative transformed calli, 13 calli were able to regenerate on the selection medium. Of the 13 calli, three transgenic tomato plantlets were recovered, and these were confirmed through PCR analysis for the presence of 432 bp fragment of ATAF1 gene. The transformation protocol reported here can be used to generate drought resistant tomato plants in future.
It is of major interest to estimate the feedback of arctic ecosystems to the global warming we expect in upcoming decades. The speed of this response is driven by the potential of species to migrate, tracking their climate optimum. For this, sessile plants have to produce and disperse seeds to newly available habitats, and pollination of ovules is needed for the seeds to be viable. These two processes are also the vectors that pass genetic information through a population. A restricted exchange among subpopulations might lead to a maladapted population due to diversity losses. Hence, a realistic implementation of these dispersal processes into a simulation model would allow an assessment of the importance of diversity for the migration of plant species in various environments worldwide. To date, dynamic global vegetation models have been optimized for a global application and overestimate the migration of biome shifts in currently warming temperatures. We hypothesize that this is caused by neglecting important fine-scale processes, which are necessary to estimate realistic vegetation trajectories. Recently, we built and parameterized a simulation model LAVESI for larches that dominate the latitudinal treelines in the northernmost areas of Siberia. In this study, we updated the vegetation model by including seed and pollen dispersal driven by wind speed and direction. The seed dispersal is modelled as a ballistic flight, and for the pollination of ovules of seeds produced, we implemented a wind-determined and distance-dependent probability distribution function using a von Mises distribution to select the pollen donor. A local sensitivity analysis of both processes supported the robustness of the model's results to the parameterization, although it highlighted the importance of recruitment and seed dispersal traits for migration rates. This individual-based and spatially explicit implementation of both dispersal processes makes it easily feasible to inherit plant traits and genetic information to assess the impact of migration processes on the genetics. Finally, we suggest how the final model can be applied to substantially help in unveiling the important drivers of migration dynamics and, with this, guide the improvement of recent global vegetation models.
The equatorial electrojet occasionally reverses during morning and afternoon hours, leading to periods of westward current in the ionospheric E region that are known as counter electrojet (CEJ) events. We present the first analysis of CEJ climatology and CEJ dependence on solar flux and lunar phase for the Brazilian sector, based on an extensive ground-based data set for the years 2008 to 2017 from the geomagnetic observatory Tatuoca (1.2 degrees S, 48.5 degrees W), and we compare it to the results found for Huancayo (12.0 degrees S, 75.3 degrees W) observatory in the Peruvian sector. We found a predominance of morning CEJ events for both sectors. The afternoon CEJ occurrence rate in the Brazilian sector is twice as high as in the Peruvian sector. The afternoon CEJ occurrence rate strongly depends on season, with maximum rates occurring during the northern-hemisphere summer for the Brazilian sector and during the northern-hemisphere winter for the Peruvian sector. Significant discrepancies between the two sectors are also found for morning CEJ rates during the northern-hemisphere summer. These longitudinal differences are in agreement with a CEJ climatology derived from contemporary Swarm satellite data and can be attributed in part to the well-known longitudinal wave-4 structure in the background equatorial electrojet strength that results from nonmigrating solar tides and stationary planetary waves. Simulations with the Thermosphere-Ionosphere-Electrodynamics General Circulation Model show that the remaining longitudinal variability in CEJ during northern summer can be explained by the effect of migrating tides in the presence of the varying geomagnetic field in the South Atlantic Anomaly.
TDS-1 GNSS Reflectometry
(2018)
This study presents the development and a systematic evaluation study of GNSS reflectometry wind speeds. After establishing a wind speed retrieval algorithm, UK TechDemoSat-1 (TDS-1) derived winds, from May 2015 to July 2017, are compared to the Advanced Scatterometer (ASCAT). ERA-Interim wind fields of the European Centre for Medium-range Weather Forecasts (ECMWF) and in situ observation from Tropical Atmosphere Ocean buoy array in the Pacific are taken as reference. One-year averaged TDS-1 global winds demonstrate small differences with ECMWF in a majority of areas as well as discuss under- and overestimations. The pioneering TDS-1 winds demonstrate a root-mean-squared error (RMSE) and bias of 2.77 and -0.33 m/s, which are comparable to the RMSE and bias derived by ASCAT winds, as large as 2.31 and 0.25 m/s, respectively. Using buoys measurements as reference, RMSE and bias of 2.23 and -0.03 m/s for TDS-1 as well as 1.40 and -0.68 m/s for ASCAT are obtained. Utilizing rain microwave-infrared estimates of the Tropical Rainfall Measuring Mission, rain-affected observation of both ASCAT and TDS-1 are collected and evaluated. Although ASCAT winds show a significant performance degradation resulting in an RMSE and bias of 3.16 and 1.03 m/s, respectively, during rain condition, TDS-1 shows a more reliable performance with an RMSE and bias of 2.94 and -0.21 m/s, respectively, which indicates the promising capability of GNSS forward scattering for wind retrievals during rain. A decrease in TDS-1-derived bistatic radar cross sections during rain events, at weak winds, is also demonstrated.
Recent studies indicate that the tropical belt has been expanding during recent decades, which can significantly influence precipitation in subtropical climates. Often the location of the tropical border is identified using the Hadley cell edge (HCE) or the subtropical jet stream (STJ), but most studies concentrated on the zonal-mean state, thereby missing regional impacts. Here we detect longitudinal-resolved trends in STJ cores and HCEs over 1979-2016 in both hemispheres at a higher spatial and temporal resolution than previous studies. Besides pronounced regional trend differences in both sign and magnitude, we show that winter HCE and STJ variability increased in the Mediterranean region and decreased over the American and Asian continents. Rainfall variability in these regions changed likewise, and most of those changes can be explained by changes in HCE/STJ variability. This highlights the importance of understanding future tropical belt changes both regionally and in terms of variability. Plain Language Summary We applied a new network-based method to detect motion of the tropical climate border with longitudinal resolution. Depending on the longitudinal position, there are differences in both direction and magnitude of the border motion. In addition, we demonstrate that the rainfall variability is increasing in the Mediterranean region and decreasing over the American and Asian continents, which can be explained by the variability of the tropical belt location.
Coupling land-use change and hydrologic models for quantification of catchment ecosystem services
(2018)
Representation of land-use and hydrologic interactions in respective models has traditionally been problematic. The use of static land-use in most hydrologic models or that of the use of simple hydrologic proxies in land-use change models call for more integrated approaches. The objective of this study is to assess whether dynamic feedback between land-use change and hydrology can (1) improve model performances, and/or (2) produce a more realistic quantification of ecosystem services. To test this, we coupled a land-use change model and a hydrologic mode. First, the land-use change and the hydrologic models were separately developed and calibrated. Then, the two models were dynamically coupled to exchange data at yearly time-steps. The approach is applied to a catchment in South Africa. Performance of coupled models when compared to the uncoupled models were marginal, but the coupled models excelled at the quantification of catchment ecosystem services more robustly.
Mechanical and/or chemical removal of material from the subsurface may generate large subsurface cavities, the destabilisation of which can lead to ground collapse and the formation of sinkholes. Numerical simulation of the interaction of cavity growth, host material deformation and overburden collapse is desirable to better understand the sinkhole hazard but is a challenging task due to the involved high strains and material discontinuities. Here, we present 2-D distinct element method numerical simulations of cavity growth and sinkhole development. Firstly, we simulate cavity formation by quasi-static, stepwise removal of material in a single growing zone of an arbitrary geometry and depth. We benchmark this approach against analytical and boundary element method models of a deep void space in a linear elastic material. Secondly, we explore the effects of properties of different uniform materials on cavity stability and sinkhole development. We perform simulated biaxial tests to calibrate macroscopic geotechnical parameters of three model materials representative of those in which sinkholes develop at the Dead Sea shoreline: mud, alluvium and salt. We show that weak materials do not support large cavities, leading to gradual sagging or suffusion-style subsidence. Strong materials support quasi-stable to stable cavities, the overburdens of which may fail suddenly in a caprock or bedrock collapse style. Thirdly, we examine the consequences of layered arrangements of weak and strong materials. We find that these are more susceptible to sinkhole collapse than uniform materials not only due to a lower integrated strength of the overburden but also due to an inhibition of stabilising stress arching. Finally, we compare our model sinkhole geometries to observations at the Ghor Al-Haditha sinkhole site in Jordan. Sinkhole depth / diameter ratios of 0.15 in mud, 0.37 in alluvium and 0.33 in salt are reproduced successfully in the calibrated model materials. The model results suggest that the observed distribution of sinkhole depth / diameter values in each material type may partly reflect sinkhole growth trends.
Linking deep seismic profiles with regional-scale gravity inversion is a powerful tool to deduce the architecture of rifted margins and their structural evolution. Here we map upper and lower crustal thicknesses of the northern South China Sea (SCS) margin in order to investigate the occurrence of depth-dependent crustal extension from the proximal to the distal margin. By comparing upper and lower crustal stretching factors, we find that the northern margin of the SCS is segmented in three parts: (1) sedimentary basins where upper crust is stretched more than lower crust, (2) distal margin where lower crust is stretched more than upper crust, (3) mostly proximal margin regions where the two layers have similar stretching factors. Our results suggest that sedimentary basins and distal margin prominently feature depth-dependent extension, however accommodated by different processes. While differential thinning within sedimentary basins appears to be governed by lateral pressure variations inducing lower crustal flow, we suggest the distal margin to be affected by a combination of mantle flow-induced lower crustal shearing and sequential fault activity during crustal hyper-extension.
Himalayan rivers are frequently hit by catastrophic floods that are caused by the failure of glacial lake and landslide dams; however, the dynamics and long-term impacts of such floods remain poorly understood. We present a comprehensive set of observations that capture the July 2016 glacial lake outburst flood (GLOF) in the Bhotekoshi/Sunkoshi River of Nepal. Seismic records of the flood provide new insights into GLOF mechanics and their ability to mobilize large boulders that otherwise prevent channel erosion. Because of this boulder mobilization, GLOF impacts far exceed those of the annual summer monsoon, and GLOFs may dominate fluvial erosion and channel-hillslope coupling many tens of kilometers downstream of glaciated areas. Long-term valley evolution in these regions may therefore be driven by GLOF frequency and magnitude, rather than by precipitation.
The ability of water to transport and transform soil materials is one of the main drivers of soil and landscape development. In turn, soil and landscape properties determine how water is distributed in soil landscapes. Understanding the complex dynamics of this co-evolution of soils, landscapes and the hydrological system is fundamental in adapting land management to changes in climate. Soil-Landscape Evolution Models (SLEMs) are used to simulate the development and evolution of soils and landscapes. However, many hydrologic processes, such as preferential flow and subsurface lateral flow, are currently absent in these models. This limits the applicability of SLEMs to improve our understanding of feedbacks in the hydro-pedo-geomorphological system. Implementation of these hydrologic processes in SLEMs faces several complications related to calculation demands, limited methods for linking pedogenic and hydrologic processes, and limited data on quantification of changes in the hydrological system over time. In this contribution, we first briefly review processes and feedbacks in soil-landscape-hydrological systems. Next, we elaborate on the development required to include these processes in SLEMs. We discuss the state-of-the-art knowledge, identify complications, give partial solutions and suggest important future development. The main requirements for incorporating hydrologic processes in SLEMs are: (1) designing a model framework that can deal with varying timescales for different sets of processes, (2) developing and implementing methods for simulating pedogenesis as a function of water flow, (3) improving and implementing knowledge on the evolution and dynamics of soil hydraulic properties over different timescales, and (4) improving the database on temporal changes and dynamics of flow paths.
Thermal conductivity (lambda) is an essential physical property of minerals and rocks and fundamental in constraining the thermal field of the lithosphere. In case that adequate samples to measure lambda are not available, it could be indirectly inferred from calculation. One of the most widely applied indirect methods for rocks involve modal mineralogy and porosity as parameters that are incorporated into mathematical mean or mixing models. Robust inferences from these approaches for crystalline rocks were impeded by a small number of studied samples or restriction to certain rock types. We employ this method and examine its applicability to low-porosity plutonic rocks by calculating bulk thermal conductivity lambda(b) for 45 samples covering the entire range from gabbro/diorite to granite. We show that the use of the harmonic-mean model for both rock matrix and porosity provided a good match between lambda(b.meas) and lambda(b.calc) of <10% deviation (2 sigma), with relative and absolute errors amounting to 1.49.7% and 4.44.9%, respectively. The results of our study constitute a big step forward to a robust conclusion on the overall applicability of the harmonic-mean model for inferring lambda(b) of isotropic, low-porosity, mafic to silicic plutonic and metamorphic rocks with an acceptable magnitude of error. Drill cuttings and enclaves form particularly interesting objects for application of this method, as they are poorly suited for direct measurement. Well-derived lambda values for those rocks would permit to calculate heat flow and to model more profoundly the thermal state of the deeper lithosphere.
Two ground motion prediction equation models for subduction zones have been tested using a public ground motion database of the KiK-net records obtained by automated processing protocols (Dawood et al., 2016, https://doi.org/10.1193/071214EQS106). The database contains records of more than 700 interface earthquakes that occurred on the Japan subduction between 1998 and 2012. The Zhao et al. (2006, https://doi.org/10.1785/0120050122) ground motion prediction equation was shown to be the best suited model for the region. It was then used as backbone to analyze the variability of ground motion records. The residuals between observed and predicted ground motions have been analyzed to study the spatial variation of the earthquakes' ground motion frequency content on the Japan megathrust. This analysis revealed a depth dependency of generated ground motions consistent with the downdip segmentation proposed for subduction interfaces (Lay et al., 2012, https://doi.org/10.1029/2011JB009133), a regional ground motion dependency that may be related with lateral variations of the mechanical properties of the subduction interface and a high-frequency radiations drop in the earthquake sequence that preceded the Tohoku-Oki earthquake Mw 9.0. The regional ground motion dependency suggests the existence of different domains along trench of the Japan subduction megathrust that control the ground motions and the wave radiation patterns of interface earthquakes. The location of their boundaries is consistent with the extension of the rupture of the 2011 Tohoku-Oki earthquake, with pre-Tohoku interseismic coupling, and with the free air gravity anomalies.
Evidence for Quaternary climate change in East Africa has been derived from outcrops on land and lake cores and from marine dust, leaf wax, and pollen records. These data have previously been used to evaluate the impact of climate change on hominin evolution, but correlations have proved to be difficult, given poor data continuity and the great distances between marine cores and terrestrial basins where fossil evidence is located. Here, we present continental coring evidence for progressive aridification since about 575 thousand years before present (ka), based on Lake Magadi (Kenya) sediments. This long-term drying trend was interrupted by many wet-dry cycles, with the greatest variability developing during times of high eccentricity-modulated precession. Intense aridification apparent in the Magadi record took place between 525 and 400 ka, with relatively persistent arid conditions after 350 ka and through to the present. Arid conditions in the Magadi Basin coincide with the Mid-Brunhes Event and overlap with mammalian extinctions in the South Kenya Rift between 500 and 400 ka. The 525 to 400 ka arid phase developed in the South Kenya Rift between the period when the last Acheulean tools are reported (at about 500 ka) and before the appearance of Middle Stone Age artifacts (by about 320 ka). Our data suggest that increasing Middle- to Late-Pleistocene aridification and environmental variability may have been drivers in the physical and cultural evolution of Homo sapiens in East Africa.
Runoff predictions in ungauged arctic basins using conceptual models forced by reanalysis data
(2018)
Due to global warming, the problem of assessing water resources and their vulnerability to climate drivers in the Arctic region has become a focus in the recent years. This study is aimed at investigating three lumped hydrological models to predict daily runoff of large-scale Arctic basins in the case of substantial data scarcity. All models were driven only by meteorological forcing reanalysis dataset without any additional information about landscape, soil, or vegetation cover properties of the studied basins. Model parameter regionalization based on transferring the whole parameter set showed good efficiency for predictions in ungauged basins. We run a blind test of the proposed methodology for ensemble runoff predictions on five sub-basins, for which only monthly observations were available, and obtained promising results for current water resources assessment for a broad domain of ungauged basins in the Russian Arctic.
The lateral and vertical temperature distribution in Oman is so far only poorly understood, particularly in the area between Muscat and the Batinah coast, which is the area of this study and which is composed of Cenozoic sediments developed as part of a foreland basin of the Makran Thrust Zone. Temperature logs (T-logs) were run and physical rock properties of the sediments were analyzed to understand the temperature distribution, thermal and hydraulic properties, and heat-transport processes within the sedimentary cover of northern Oman. An advective component is evident in the otherwise conduction-dominated geothermal play system, and is caused by both topography and density driven flow. Calculated temperature gradients (T-gradients) in two wells that represent conductive conditions are 18.7 and 19.5 °C km−1, corresponding to about 70–90 °C at 2000–3000 m depth. This indicates a geothermal potential that can be used for energy intensive applications like cooling or water desalinization. Sedimentation in the foreland basin was initiated after the obduction of the Semail Ophiolite in the late Campanian, and reflects the complex history of alternating periods of transgressive and regressive sequences with erosion of the Oman Mountains. Thermal and hydraulic parameters were analyzed of the basin’s heterogeneous clastic and carbonate sedimentary sequence. Surface heat-flow values of 46.4 and 47.9 mW m−2 were calculated from the T-logs and calculated thermal conductivity values in two wells. The results of this study serve as a starting point for assessing different geothermal applications that may be suitable for northern Oman.
The phi(ev) is calculated from high-resolution discharge and precipitation data for several rain events with a cumulative precipitation P-cum ranging from less than 5mm to more than 80 mm. Because of the high uncertainty of phi(ev) associated with the hydrograph separation method, phi(ev) is calculated with several methods, including graphical methods, digital filters and a tracer-based method. The results indicate that the hydrological response depends on (theta) over bar (ini): during dry conditions phi(ev) is consistently below 0.1, even for events with high and intense precipitation. Above a threshold of (theta) over bar (ini) = 34 vol % phi(ev) can reach values up to 0.99 but there is a high scatter. Some variability can be explained with a weak correlation of phi(ev) with P-cum and rain intensity, but a considerable part of the variability remains unexplained. It is concluded that threshold-based methods can be helpful to prevent overestimation of the hydrological response during dry catchment conditions. The impact of soil moisture on the hydrological response during wet catchment conditions, however, is still insufficiently understood and cannot be generalized based on the present results.
For several decades the Sant Llorenc de Montgai reservoir has experienced different problems that could affect the safety of the engineering structure. For this reason, several corrective actions have been taken over the years. Here, we present a study involving complementary geophysical methods including electrical resistivity tomography, seismic refraction tomography and frequency-domain electromagnetic surveys. The analysis of the inverted electrical resistivity tomography cross-sections combined with the seismic refraction results and land subsidence monitoring data show the likely mechanism of abnormal seepage. The areas where mortar injections were applied as a corrective measure are also clearly delineated. In addition, the evolution of the state of the embankment has been established from two successive electrical resistivity tomography surveys in the last two decades. The results show areas where corrective mortar injections have been effective, while in other areas new abnormal seepage has been detected. The lithological heterogeneity of the bedrock, especially the dissolution of gypsum-rich rocks, induced subsidence effects and caused abnormal seepage in different areas along the embankment. Our results indicate how corrective solutions can be optimized to reduce the cost of corrective engineering interventions.
For the first time, a rain signature in Global Navigation Satellite System Reflectometry (GNSS-R) observations is demonstrated. Based on the argument that the forward quasi-specular scattering relies upon surface gravity waves with lengths larger than several wavelengths of the reflected signal, a commonly made conclusion is that the scatterometric GNSS-R measurements are not sensitive to the surface small-scale roughness generated by raindrops impinging on the ocean surface. On the contrary, this study presents an evidence that the bistatic radar cross section sigma(0) derived from TechDemoSat-1 data is reduced due to rain at weak winds, lower than approximate to 6 m/s. The decrease is as large as approximate to 0.7 dB at the wind speed of 3 m/s due to a precipitation of 0-2 mm/hr. The simulations based on the recently published scattering theory provide a plausible explanation for this phenomenon which potentially enables the GNSS-R technique to detect precipitation over oceans at low winds.
Meteorological conditions during the ACLOUD/PASCAL field campaign near Svalbard in early summer 2017
(2018)
The two concerted field campaigns, Arctic CLoud Observations Using airborne measurements during polar Day (ACLOUD) and the Physical feedbacks of Arctic planetary boundary level Sea ice, Cloud and AerosoL (PASCAL), took place near Svalbard from 23 May to 26 June 2017. They were focused on studying Arctic mixed-phase clouds and involved observations from two airplanes (ACLOUD), an icebreaker (PASCAL) and a tethered balloon, as well as ground-based stations. Here, we present the synoptic development during the 35-day period of the campaigns, using near-surface and upper-air meteorological observations, as well as operational satellite, analysis, and reanalysis data. Over the campaign period, short-term synoptic variability was substantial, dominating over the seasonal cycle. During the first campaign week, cold and dry Arctic air from the north persisted, with a distinct but seasonally unusual cold air outbreak. Cloudy conditions with mostly low-level clouds prevailed. The subsequent 2 weeks were characterized by warm and moist maritime air from the south and east, which included two events of warm air advection. These synoptical disturbances caused lower cloud cover fractions and higher-reaching cloud systems. In the final 2 weeks, adiabatically warmed air from the west dominated, with cloud properties strongly varying within the range of the two other periods. Results presented here provide synoptic information needed to analyze and interpret data of upcoming studies from ACLOUD/PASCAL, while also offering unprecedented measurements in a sparsely observed region.
Great megathrust earthquakes arise from the sudden release of energy accumulated during centuries of interseismic plate convergence. The moment deficit (energy available for future earthquakes) is commonly inferred by integrating the rate of interseismic plate locking over the time since the previous great earthquake. But accurate integration requires knowledge of how interseismic plate locking changes decades after earthquakes, measurements not available for most great earthquakes. Here we reconstruct the post-earthquake history of plate locking at Guafo Island, above the seismogenic zone of the giant 1960 (M-w = 9.5) Chile earthquake, through forward modeling of land-level changes inferred from aerial imagery (since 1974) and measured by GPS (since 1994). We find that interseismic locking increased to similar to 70% in the decade following the 1960 earthquake and then gradually to 100% by 2005. Our findings illustrate the transient evolution of plate locking in Chile, and suggest a similarly complex evolution elsewhere, with implications for the time- and magnitude-dependent probability of future events.
Subduction zone magmas are more oxidised on eruption than those at mid-ocean ridges. This is attributed either to oxidising components, derived from subducted lithosphere (slab) and added to the mantle wedge, or to oxidation processes occurring during magma ascent via differentiation. Here we provide direct evidence for contributions of oxidising slab agents to melts trapped in the sub-arc mantle. Measurements of sulfur (S) valence state in sub-arc mantle peridotites identify sulfate, both as crystalline anhydrite (CaSO4) and dissolved SO42- in spinel-hosted glass (formerly melt) inclusions. Copper-rich sulfide precipitates in the inclusions and increased Fe3+/Sigma Fe in spinel record a S6+-Fe2+ redox coupling during melt percolation through the sub-arc mantle. Sulfate-rich glass inclusions exhibit high U/Th, Pb/Ce, Sr/Nd and delta S-34 (+ 7 to + 11%), indicating the involvement of dehydration products of serpentinised slab rocks in their parental melt sources. These observations provide a link between liberated slab components and oxidised arc magmas.
Micropetrology
(2018)
Inclusions in minerals, whether fluids, melts or crystalline phases, are small pieces of the large-scale puzzle of Nature, time-consuming to investigate and often of difficult interpretation. Yet they are windows into the past of their host mineral. Mineral inclusions provide the opportunity to unravel the genesis of their host, and the increasingly refined understanding of their elastic behaviour provides the basis for alternative, equilibrium-independent geobarometry. Fluid and melt inclusions reveal information about material transfer in the Earth system, from shallow mineralization to mantle re-fertilization via subduction. The study of inclusions is thus one of the most intriguing and fertile branches of micropetrology. In this contribution, we focus on two recent developments: the use of elasticity models to extract the formation conditions of the host crystal, and the discovery and investigation of melt inclusions in metamorphic rocks. We also discuss how to evaluate the information provided by inclusions, given that they are no longer at the pressure and temperature conditions of entrapment. We discuss how to understand and quantify the changes undergone during cooling and depressurization, and how metastability-related phenomena in inclusions, such as crystallization of rare polymorphs and preservation of the original content of volatiles in fluid and melt inclusions, provide direct evidence that inclusions represent closed systems. The field of study of inclusions in minerals still has a largely untapped potential. The most fruitful avenues for future research will emerge from continuous technological innovation in analytical and imaging techniques, the application of experimental petrology, and the development and application of new theoretical models for coupled mineral behaviour under changing P-T conditions.
Originally developed for use in the petroleum industry, Rock-Eval pyrolysis is a technique commonly applied to lake sediments to infer paleoenvironmental reconstructions. The standard Rock-Eval parameters provide information on the amount of total organic and inorganic carbon (TOC and MinC, respectively), and are usually interpreted as proxies for the source (aquatic or terrestrial) of the primary production of organic matter (Hydrogen Index vs Oxygen Index). Although this method usually provides valuable evidence, the common presence of siderite in tropical lake sediments can alter the primary signal of the sedimentary organic matter (SOM). Indeed, the CO2 and CO released by the pyrolysis of siderite are integral to the calculation of the SOM-related standard Rock-Eval parameters. In this study, we analyze sediments from a core collected in the Lake Barombi (southwest Cameroon) and describe the impact of siderite on standard Rock-Eval parameters. We propose a workflow that allows standard Rock-Eval parameters to be corrected, based on the analysis of thermograms. The proposed corrections provide siderite-effect-free parameters, accurately reflecting the changes in sedimentary organic matter composition. (C) 2018 Elsevier Ltd. All rights reserved.
Abrupt or gradual?
(2018)
We used a change point analysis on a late Pleistocene-Holocene lake-sediment record from the Chew Bahir basin in the southern Ethiopian Rift to determine the amplitude and duration of past climate transitions. The most dramatic changes occurred over 240 yr (from similar to 15,700 to 15,460 yr) during the onset of the African Humid Period (AHP), and over 990 yr (from similar to 4875 to 3885 yr) during its protracted termination. The AHP was interrupted by a distinct dry period coinciding with the high-latitude Younger Dryas stadial, which had an abrupt onset (less than similar to 100 yr) at similar to 13,260 yr and lasted until similar to 11,730 yr. Wet-dry-wet transitions prior to the AHP may reflect the high-latitude Dansgaard-Oeschger cycles, as indicated by cross-correlation of the potassium record with the NorthGRIP ice core record between similar to 45-20 ka. These findings may contribute to the debates regarding the amplitude, and duration and mechanisms of past climate transitions, and their possible influence on the development of early modern human cultures.
By using 3-year global positioning system (GPS) measurements from December 2013 to November 2016, we provide in this study a detailed survey on the climatology of the GPS signal loss of Swarm onboard receivers. Our results show that the GPS signal losses prefer to occur at both low latitudes between +/- 5 and +/- 20 degrees magnetic latitude (MLAT) and high latitudes above 60 degrees MLAT in both hemispheres. These events at all latitudes are observed mainly during equinoxes and December solstice months, while totally absent during June solstice months. At low latitudes the GPS signal losses are caused by the equatorial plasma irregularities shortly after sunset, and at high latitude they are also highly related to the large density gradients associated with ionospheric irregularities. Additionally, the high-latitude events are more often observed in the Southern Hemisphere, occurring mainly at the cusp region and along nightside auroral latitudes. The signal losses mainly happen for those GPS rays with elevation angles less than 20 degrees, and more commonly occur when the line of sight between GPS and Swarm satellites is aligned with the shell structure of plasma irregularities. Our results also confirm that the capability of the Swarm receiver has been improved after the bandwidth of the phase-locked loop (PLL) widened, but the updates cannot radically avoid the interruption in tracking GPS satellites caused by the ionospheric plasma irregularities. Additionally, after the PLL bandwidth increased larger than 0.5 Hz, some unexpected signal losses are observed even at middle latitudes, which are not related to the ionospheric plasma irregularities. Our results suggest that rather than 1.0 Hz, a PLL bandwidth of 0.5 Hz is a more suitable value for the Swarm receiver.
Numerous charcoal accumulations discovered in the submerged Chan Hol cave near Tulum, Quintana Roo, Mexico, have been C-14-dated revealing ages between 8110 +/- 28 C-14 a BP (9122-8999 cal a BP) and 7177 +/- 27 C-14 a BP (8027-7951 cal a BP). These charcoal concentrations, interpreted here as ancient illumination sites, provide strong evidence that the Chan Hol cave was dry and accessible during that time interval. Humans used the cave for at least 1200 years during the early and middle Holocene, before access was successively interrupted by global sea level rise and flooding of the cave system. Our data thus narrow the gap between an early settlement in the Tulum area reaching from the late Pleistocene (similar to 13 000 a) to middle Holocene (e.g. 7177 C-14 a BP), and the Maya Formative period at approximately 3000 a bp. Yet, no evidence has been presented to date for human settlement during the similar to 4000-year interval between 7000 and 3000 a. This is remarkable as settlement in other areas of south-eastern Mexico (e.g. Chiapas, Tabasco) and in Guatemala was apparently continuous.
Modern mobile devices (i.e. smartphones and tablet computers) are widespread, everyday tools, which are equipped with a variety of sensors including three-axis magnetometers. Here, we investigate the feasibility and the potential of using such mobile devices to mimic geophysical experiments in the classroom in a table-top setup. We focus on magnetic surveying and present a basic setup of a table-top experiment for collecting three-component magnetic data across well-defined source bodies and structures. Our results demonstrate that the quality of the recorded data is sufficient to address a number of important basic concepts in the magnetic method. The shown examples cover the analysis of magnetic data recorded across different kinds of dipole sources, thus illustrating the complexity of magnetic anomalies. In addition, we analyze the horizontal resolution capabilities using a pair of dipole sources placed at different horizontal distances to each other. Furthermore, we demonstrate that magnetic data recorded with a mobile device can even be used to introduce filtering, transformation, and inversion approaches as they are typically used when processing magnetic data sets recorded for real-world field applications. Thus, we conclude that such table-top experiments represent an easy-to-implement experimental procedure (as student exercise or classroom demonstration) and can provide first hands-on experience in the basic principles of magnetic surveying including the fundamentals of data acquisition, analysis and processing, as well as data evaluation and interpretation.
Recent improvements in the development of VLBI (very long baseline interferometry) and other space geodetic techniques such as the global navigation satellite systems (GNSS) require very precise a-priori information of short-period (daily and sub-daily) Earth rotation variations. One significant contribution to Earth rotation is caused by the diurnal and semi-diurnal ocean tides. Within this work, we developed a new model for the short-period ocean tidal variations in Earth rotation, where the ocean tidal angular momentum model and the Earth rotation variation have been setup jointly. Besides the model of the short-period variation of the Earth’s rotation parameters (ERP), based on the empirical ocean tide model EOT11a, we developed also ERP models, that are based on the hydrodynamic ocean tide models FES2012 and HAMTIDE. Furthermore, we have assessed the effect of uncertainties in the elastic Earth model on the resulting ERP models. Our proposed alternative ERP model to the IERS 2010 conventional model considers the elastic model PREM and 260 partial tides. The choice of the ocean tide model and the determination of the tidal velocities have been identified as the main uncertainties. However, in the VLBI analysis all models perform on the same level of accuracy. From these findings, we conclude that the models presented here, which are based on a re-examined theoretical description and long-term satellite altimetry observation only, are an alternative for the IERS conventional model but do not improve the geodetic results.
Detections of pP and sP phase arrivals (the so-called depth phases) at teleseismic distance provide one of the best ways to estimate earthquake focal depth, as the P-pP and the P-sP delays are strongly dependent on the depth. Based on a new processing workflow and using a single seismic array at teleseismic distance, we can estimate the depth of clusters of small events down to magnitude M-b 3.5. Our method provides a direct view of the relative variations of the seismicity depth from an active area. This study focuses on the application of this new methodology to study the lateral variations of the Guerrero subduction zone (Mexico) using the Eielson seismic array in Alaska (USA). After denoising the signals, 1232 M-b 3.5 + events were detected, with clear P, pP, sP and PcP arrivals. A high-resolution view of the lateral variations of the depth of the seismicity of the Guerero-Oaxaca area is thus obtained. The seismicity is shown to be mainly clustered along the interface, coherently following the geometry of the plate as constrained by the receiver-function analysis along the Meso America Subduction Experiment profile. From this study, the hypothesis of tears on the western part of Guerrero and the eastern part of Oaxaca are strongly confirmed by dramatic lateral changes in the depth of the earthquake clusters. The presence of these two tears might explain the observed lateral variations in seismicity, which is correlated with the boundaries of the slow slip events.
A transparent and data-driven global tectonic regionalization model for seismic hazard assessment
(2018)
A key concept that is common to many assumptions inherent within seismic hazard assessment is that of tectonic similarity. This recognizes that certain regions of the globe may display similar geophysical characteristics, such as in the attenuation of seismic waves, the magnitude scaling properties of seismogenic sources or the seismic coupling of the lithosphere. Previous attempts at tectonic regionalization, particularly within a seismic hazard assessment context, have often been based on expert judgements; in most of these cases, the process for delineating tectonic regions is neither reproducible nor consistent from location to location. In this work, the regionalization process is implemented in a scheme that is reproducible, comprehensible from a geophysical rationale, and revisable when new relevant data are published. A spatial classification-scheme is developed based on fuzzy logic, enabling the quantification of concepts that are approximate rather than precise. Using the proposed methodology, we obtain a transparent and data-driven global tectonic regionalization model for seismic hazard applications as well as the subjective probabilities (e.g. degree of being active/degree of being cratonic) that indicate the degree to which a site belongs in a tectonic category.
Dumortierite was synthesized in piston-cylinder experiments at 2.5-4.0 GPa, 650-700 degrees C in the Al2O3 -B2O3-SiO2-H2O (ABSH) system. Electron-microprobe (EMP) analyses reveal significant boron-excess (up to 0.26 B-[4] per formula unit, pfu) and silicon-deficiency relative to the ideal anhydrous dumortierite stoichiometry Al7BSi3O18 . The EMP data in conjunction with results from single-crystal Raman spectroscopy and powder X-ray diffraction provide evidence that silicon at the tetrahedral site is replaced by excess boron via the substitution Si-[4] <--> B-[4] + H. The Raman spectrum of synthetic dumortierite in the frequency region 2000 4000 cm(-1) comprises eight bands, of which six are located at frequencies below 3400 cm(-1). This points to strong hydrogen bonding, most likely O2-H center dot center dot center dot O7 and O7-H center dot center dot center dot O2, arising from a high number of octahedral vacancies at the All site and substitution of trivalent Al3+ and B3+ for Si4+ at Si1 and Si2 sites, causing decreasing acceptor-donor distances and lower incident valence at the acceptor oxygen. Contrary to the synthetic high-pressure ABSH-dumortierite, magnesiodumortierite from the Dora-Maira Massif, which is assumed to have formed at similar conditions (2.5-3.0 GPa, 700 degrees C), does not show any B-excess. Tourmaline shows an analogous behaviour in that magnesium-rich (e.g., dravitic) tourmaline formed at high pressure shows no or only minor amounts of tetrahedral boron, whereas natural aluminum-rich tourmaline and synthetic olenitic tourmaline formed at high pressures can incorporate significant amounts of tetrahedral boron. Two mechanisms might account for this discrepancy: (i) Structural avoidance of Mg-[6]-(OR3+)-R-[4] configurations in magnesiodumortierite due to charge deficieny at the oxygens O2 and O7 and strong local distortion of M1 due to decreased O2-O7 bond length, and/or (ii) decreasing fluid mobility of boron in Al-rich systems at high pressures.
Python is at the forefront of scientific computation for seismologists and therefore should be introduced to students interested in becoming seismologists. On its own, Python is open source and well designed with extensive libraries. However, Python code can also be executed, visualized, and communicated to others with "Jupyter Notebooks". Thus, Jupyter Notebooks are ideal for teaching students Python and scientific computation. In this article, we designed an openly available Python library and collection of Jupyter Notebooks based on defined scientific computation learning goals for seismology students. The Notebooks cover topics from an introduction to Python to organizing data, earthquake catalog statistics, linear regression, and making maps. Our Python library and collection of Jupyter Notebooks are meant to be used as course materials for an upper-division data analysis course in an Earth Science Department, and the materials were tested in a Probabilistic Seismic Hazard course. However, seismologists or anyone else who is interested in Python for data analysis and map making can use these materials.
The equatorial electrojet is an enhanced eastward current in the dayside E region ionosphere flowing along the magnetic equator. The equatorial electrojet is highly variable as it is subject to various forcing mechanisms including atmospheric waves from the lower layers of the atmosphere. There are occasionally times when the intensity of the equatorial electrojet at a fixed longitude shows an oscillatory variation with a period of approximately 6days. We present case studies of such events based on the equatorial electrojet measurements from the CHAMP and Swarm satellites. The spatial and temporal variability of the equatorial electrojet intensity during these events reveals characteristics of a westward propagating wave with zonal wavenumber 1, consistent with the effect of the quasi-6-day planetary wave. Analyses of the geopotential height data from the Aura satellite confirm the presence of the quasi-6-day planetary wave in the lower thermosphere during the events. The amplitude of the quasi-6-day variation in the equatorial electrojet intensity depends on longitude, but no systematic longitudinal dependence is found for different events. During the event of August 2010, quasi-6-day variations are also observed by ground-based magnetometers and a radar in the Peruvian sector. The effect of the quasi-6-day wave accounts for up to +/- 5.9m/s in the equatorial vertical plasma velocity at noon, which is much larger than previously predicted by a numerical model. These results suggest that the quasi-6-day planetary wave is an important source of short-term variability in the equatorial ionosphere.
The Gongjue basin from the eastern Qiangtang terrane is located in the transition region where the regional structural lineation curves from east-west-oriented in Tibet to north-south-oriented in Yunnan. In this study, we sampled the red beds in the basin from the lower Gongjue to upper Ranmugou formations for the first time covering the entire stratigraphic profile. The stratigraphic ages are bracketed within 53-43Ma by new detrital zircon U-Pb ages constraining the maximum deposition age to 52.51.5Ma. Rock magnetic and petrographic studies indicate that detrital magnetite and hematite are the magnetic carriers. Positive reversals and fold tests demonstrate that the characteristic remanent magnetization has a primary origin. The Gongjue and Ranmugou formations yield mean characteristic remanent magnetization directions of D-s/I-s=31.0 degrees/21.3 degrees and D-s/I-s=15.9 degrees/22.0 degrees, respectively. The magnetic inclination of these characteristic remanent magnetizations is significantly shallowed compared to the expected inclination for the locality. However, the elongation/inclination correction method does not provide a meaningful correction, likely because of syn-depositional rotation. Rotations relative to the Eurasian apparent polar wander path occurred in three stages: Stage I, 33.33.4 degrees clockwise rotation during the deposition of the Gongjue and lower Ranmugou formations; Stage II, 26.93.7 degrees counterclockwise rotation during deposition of the lower and middle Ranmugou formation; and Stage III, 17.73.3 degrees clockwise rotation after 43Ma. The complex rotation history recorded in the basin is possibly linked to sinistral shear along the Qiangtang block during India indentation into Asia and the early stage of the extrusion of the northwestern Indochina blocks away from eastern Tibet.
Fibroblast growth factor 23 (FGF23) is a proteohormone regulating renal phosphate transport and vitamin D metabolism as well as inducing left heart hypertrophy. FGF23-deficient mice suffer from severe tissue calcification, accelerated aging and a myriad of aging-associated diseases. Bone cells produce FGF23 upon store-operated calcium ion entry (SOCE) through the calcium selective ion channel Orai1. AMP-activated kinase (AMPK) is a powerful energy sensor helping cells survive states of energy deficiency, and AMPK down-regulates Orai1. Here we investigated the role of AMPK in FGF23 production. Fgf23 gene transcription was analyzed by qRT-PCR and SOCE by fluorescence optics in UMR106 osteoblast-like cells while the serum FGF23 concentration and phosphate metabolism were assessed in AMPKa1-knockout and wild-type mice. The AMPK activator, 5-aminoimidazole-4-carboxamide ribonucleotide (AICAR) down-regulated, whereas the AMPK inhibitor, dorsomorphin dihydrochloride (compound C) and AMPK gene silencing induced Fgf23 transcription. AICAR decreased membrane abundance of Orai1 and SOCE. SOCE inhibitors lowered Fgf23 gene expression induced by AMPK inhibition. AMPKa1-knockout mice had a higher serum FGF23 concentration compared to wild-type mice. Thus, AMPK participates in the regulation of FGF23 production in vitro and in vivo. The inhibitory effect of AMPK on FGF23 production is at least in part mediated by Orai1-involving SOCE.
A submerged pine forest from the early Holocene in the Mecklenburg Lake District, northern Germany
(2018)
For the first time, evidence of a submerged pine forest from the early Holocene can be documented in a central European lake. Subaquatic tree stumps were discovered in Lake Giesenschlagsee at a depth of between 2 and 5m using scuba divers, side-scan sonar and a remotely operated vehicle. Several erect stumps, anchored to the ground by roots, represent an insitu record of this former forest. Botanical determination revealed the stumps to be Scots pine (Pinus sylvestris) with an individual tree age of about 80years. The trees could not be dated by means of dendrochronology, as they are older than the regional reference chronology for pine. Radiocarbon ages from the wood range from 10880 +/- 210 to 10370 +/- 130cal. a BP, which is equivalent to the mid-Preboreal to early Boreal biozones. The trees are rooted in sedge peat, which can be dated to this period as well, using pollen stratigraphical analysis. Tilting of the peat bed by 4m indicates subsidence of the ground due to local dead ice melting, causing the trees to become submerged and preserved for millennia. Together with recently detected Lateglacial insitu tree occurrences in nearby lakes, the submerged pine forest at Giesenschlagsee represents a new and highly promising type of geo-bio-archive for the wider region. Comparable insitu pine remnants occur at some terrestrial (buried setting) and marine (submerged setting) sites in northern central Europe and beyond, but they partly differ in age. In general, the insitu pine finds document shifts of the zonal boreal forest ecosystem during the late Quaternary.
The potential link between erosion rates at the Earth’s surface and changes in global climate has intrigued geoscientists for decades1,2 because such a coupling has implications for the influence of silicate weathering3,4 and organic-carbon burial5 on climate and for the role of Quaternary glaciations in landscape evolution1,6. A global increase in late-Cenozoic erosion rates in response to a cooling, more variable climate has been proposed on the basis of worldwide sedimentation rates7. Other studies have indicated, however, that global erosion rates may have remained steady, suggesting that the reported increases in sediment-accumulation rates are due to preservation biases, depositional hiatuses and varying measurement intervals8,9,10. More recently, a global compilation of thermochronology data has been used to infer a nearly twofold increase in the erosion rate in mountainous landscapes over late-Cenozoic times6. It has been contended that this result is free of the biases that affect sedimentary records11, although others have argued that it contains biases related to how thermochronological data are averaged12 and to erosion hiatuses in glaciated landscapes13. Here we investigate the 30 locations with reported accelerated erosion during the late Cenozoic6. Our analysis shows that in 23 of these locations, the reported increases are a result of a spatial correlation bias—that is, combining data with disparate exhumation histories, thereby converting spatial erosion-rate variations into temporal increases. In four locations, the increases can be explained by changes in tectonic boundary conditions. In three cases, climatically induced accelerations are recorded, driven by localized glacial valley incision. Our findings suggest that thermochronology data currently have insufficient resolution to assess whether late-Cenozoic climate change affected erosion rates on a global scale. We suggest that a synthesis of local findings that include location-specific information may help to further investigate drivers of global erosion rates.
The effects of varying LiPF6 salt concentration and the presence of lithium bis(oxalate)borate additive on the electronic structure of commonly used lithium-ion battery electrolyte solvents (ethylene carbonate-dimethyl carbonate and propylene carbonate) have been investigated. X-ray Raman scattering spectroscopy (a non-resonant inelastic X-ray scattering method) was utilized together with a closed-circle flow cell. Carbon and oxygen K-edges provide characteristic information on the electronic structure of the electrolyte solutions, which are sensitive to local chemistry. Higher Li+ ion concentration in the solvent manifests itself as a blue-shift of both the pi* feature in the carbon edge and the carbonyl pi* feature in the oxygen edge. While these oxygen K-edge results agree with previous soft X-ray absorption studies on LiBF4 salt concentration in propylene carbonate, carbon K-edge spectra reveal a shift in energy, which can be explained with differing ionic conductivities of the electrolyte solutions.
In catastrophe risk modeling, a defensible estimation of impact severity and its likelihood of occurrence to a portfolio of assets can only be made through a rigorous treatment of uncertainty and the consideration of multiple alternative models. This approach, however, requires repeating lengthy calculations multiple times. To limit the demand on computational time and resources, a frequent practice in the industry is to estimate the distribution of earthquake-induced portfolio losses using a simulated catalog of events from a single representative mean ground motion hazard model for the region. This simplified approach is faster but may provide biased estimates of the likelihood of occurrence of the large and infrequent losses that drive many risk mitigation decisions. Investigation through case studies of different portfolios of assets located in the San Francisco Bay Region shows the potential for both a bias in the mean loss estimates and an underestimation of their central 70% interpercentile. We propose a simplified and computationally practical approach that reduces the bias in the mean portfolio loss estimates. This approach does not improve the estimate of the interpercentile range, however, a quantity of no direct practical use.
To evaluate the spatiotemporal variations of ground motions in northern Chile, we built a high-quality rock seismic acceleration database and an interface earthquakes catalog. Two ground-motion prediction equation (GMPE) models for subduction zones have been tested and validated for the area. They were then used as backbone models to describe the time-space variations of earthquake frequency content (Fourier and response spectra). Consistent with previous studies of large subduction earthquakes, moderate interface earthquakes in northern Chile show an increase of the high-frequency energy released with depth. A regional variability of earthquake frequency content is also observed, which may be related to a lateral segmentation of the mechanical properties of the subduction interface. Finally, interface earthquakes show a temporal evolution of their frequency content in the earthquake sequence associated with the 2014 Iquique M-w 8.1 megathrust earthquake. Surprisingly, the change does not occur with the mainshock but is associated with an 8 month slow slip preceding the megathrust. Electronic Supplement: Strong-motion database.
The Timpa delle Murge ophiolite in the North Calabrian Unit is part of the Liguride Complex (southern Apennines). The study is concentrated on the gabbroic part of the ophiolite of the Pollino area. They preserve the high-grade ocean floor metamorphic and locally developed flaser textures under ocean floor conditions. The primary magmatic assemblages are clinopyroxene, plagioclase, and opaques. Brown hornblende is a late magmatic phase. Green hornblende, actinolite, albite, chlorite and epidote display metamorphic recrystallization under lower amphibolite facies conditions, followed by greenschist facies. The gabbros show subalkaline near to alkaline character with a tendency to a more calkalkaline trend. The normalization to primitive mantle and mid-ocean ridge basalt (N-MORB) compositions indicates a considerable depletion in Nb, P, Zr and Ti and an enrichment in Ba, Rb, K, Sr and Eu. This points to a mantle source, which is not compatible with a "normal" mid-ocean ridge situation. Rather, the gabbros are generated from a N-MORB-like melt with a strong crustal component, which was influenced by subduction related fluids and episodic melting during mid-ocean-ridge processes. Plausible geodynamic settings of the Timpa delle Murge gabbros are oceanic back-arc positions with embryonic MORB-activities. Similar slab contaminated magmatism is also known from the early stage of island arc formation in supra-subduction zone environments like the Izu-Bonin-Mariana island arc.
In this study, we developed an autodetection technique for the equatorial plasma depletions (EPDs) and their occurrence and depletion amplitudes based on in situ electron density measurements gathered by Swarm A satellite. For the first time, comparisons are made among the detected EPDs and their amplitudes with the loss of Global Positioning System (GPS) signal of receivers onboard Swarm A, and the Swarm Level-2 product, Ionospheric Bubble Index (IBI). It has been found that the highest rate of EPD occurrence takes place generally between 2200 and 0000 magnetic local time (MLT), in agreement with the IBI. However, the largest amplitudes of EPD are detected earlier at about 1900-2100 MLT. This coincides with the moment of higher background electron density and the largest occurrence of GPS signal loss. From a longitudinal perspective, the higher depletion amplitude is always witnessed in spatial bins with higher background electron density. At most longitudes, the occurrence rate of postmidnight EPDs is reduced compared to premidnight ones; while more postmidnight EPDs are observed at African longitudes. CHAMP observations confirm this point regardless of high or low solar activity condition. Further by comparing with previous studies and the plasma vertical drift velocity from ROCSAT-1, we suggest that while the F region vertical plasma drift plays a key role in dominating the occurrence of EPDs during premidnight hours, the postmidnight EPDs are the combined results from the continuing of former EPDs and newborn EPDs, especially during June solstice. And these newborn EPDs during postmidnight hours seem to be less related to the plasma vertical drift.
One of the main purposes of detrital thermochronology is to provide constraints on the regional-scale exhumation rate and its spatial variability in actively eroding mountain ranges. Procedures that use cooling age distributions coupled with hypsometry and thermal models have been developed in order to extract quantitative estimates of erosion rate and its spatial distribution, assuming steady state between tectonic uplift and erosion. This hypothesis precludes the use of these procedures to assess the likely transient response of mountain belts to changes in tectonic or climatic forcing. Other methods are based on an a priori knowledge of the in situ distribution of ages to interpret the detrital age distributions. In this paper, we describe a simple method that, using the observed detrital mineral age distributions collected along a river, allows us to extract information about the relative distribution of erosion rates in an eroding catchment without relying on a steady-state assumption, the value of thermal parameters or an a priori knowledge of in situ age distributions. The model is based on a relatively low number of parameters describing lithological variability among the various sub-catchments and their sizes and only uses the raw ages. The method we propose is tested against synthetic age distributions to demonstrate its accuracy and the optimum conditions for it use. In order to illustrate the method, we invert age distributions collected along the main trunk of the Tsangpo-Siang-Brahmaputra river system in the eastern Himalaya. From the inversion of the cooling age distributions we predict present-day erosion rates of the catchments along the Tsangpo-Siang-Brahmaputra river system, as well as some of its tributaries. We show that detrital age distributions contain dual information about present-day erosion rate, i. e., from the predicted distribution of surface ages within each catchment and from the relative contribution of any given catchment to the river distribution. The method additionally allows comparing modern erosion rates to long-term exhumation rates. We provide a simple implementation of the method in Python code within a Jupyter Notebook that includes the data used in this paper for illustration purposes.
Stable isotope ratios delta O-18 and delta D in polar ice provide a wealth of information about past climate evolution. Snow-pit studies allow us to relate observed weather and climate conditions to the measured isotope variations in the snow. They therefore offer the possibility to test our understanding of how isotope signals are formed and stored in firn and ice. As delta O-18 and delta D in the snowfall are strongly correlated to air temperature, isotopes in the near-surface snow are thought to record the seasonal cycle at a given site. Accordingly, the number of seasonal cycles observed over a given depth should depend on the accumulation rate of snow. However, snow-pit studies from different accumulation conditions in East Antarctica reported similar isotopic variability and comparable apparent cycles in the delta O-18 and delta D profiles with typical wavelengths of similar to 20 cm. These observations are unexpected as the accumulation rates strongly differ between the sites, ranging from 20 to 80mmw.e.yr(-1) (similar to 6-21 cm of snow per year). Various mechanisms have been proposed to explain the isotopic variations individually at each site; however, none of these are consistent with the similarity of the different profiles independent of the local accumulation conditions. Here, we systematically analyse the properties and origins of delta O-18 and delta D variations in high-resolution firn profiles from eight East Antarctic sites. First, we confirm the suggested cycle length (mean distance between peaks) of similar to 20 cm by counting the isotopic maxima. Spectral analysis further shows a strong similarity between the sites but indicates no dominant periodic features. Furthermore, the appar-ent cycle length increases with depth for most East Antarctic sites, which is inconsistent with burial and compression of a regular seasonal cycle. We show that these results can be explained by isotopic diffusion acting on a noise-dominated isotope signal. The firn diffusion length is rather stable across the Antarctic Plateau and thus leads to similar power spectral densities of the isotopic variations. This in turn implies a similar distance between isotopic maxima in the firn profiles. Our results explain a large set of observations discussed in the literature, providing a simple explanation for the interpretation of apparent cycles in shallow isotope records, without invoking complex mechanisms. Finally, the results underline previous suggestions that isotope signals in single ice cores from low-accumulation regions have a small signal-to-noise ratio and thus likely do not allow the reconstruction of interannual to decadal climate variations.
There has been recent progress in the understanding of the evolution of Quaternary climate. Simultaneously, there have been improvements in the understanding of glacial erosion processes, with better parameter constraints. Despite this, there remains much debate about whether or not the observed cooling over the Quaternary has driven an increase in glacial erosion rates. Most studies agree that the erosional response to climate change must be transient; therefore, the time scale of the climatic change and the response time of glacial erosion must be accounted for. Here we analyze the equations governing glacial erosion in a steadily uplifting landscape with variable climatic forcing and derive expressions for two fundamental response time scales. The first time scale describes the response of the glacier and the second one the glacial erosion response. We find that glaciers have characteristic time scales of the order of 10 to 10,000 years, while the characteristic time scale for glacial erosion is of the order of a few tens of thousands to a few million years. We then use a numerical model to validate the approximations made to derive the analytical solutions. The solutions show that short period forcing is dampened by the glacier response time, and long period forcing (>1 Myr) may be dampened by erosional response of glaciers when the rock uplift rates are high. In most tectonic and climatic conditions, we expect to see the strongest response of glacial erosion to periodic climatic forcing corresponding to Plio-Pleistocene climatic cycles. Finally, we use the numerical model to predict the response of glacial systems to the observed climatic forcing of the Quaternary, including, but not limited to, the Milankovich periods and the long-term secular cooling trend. We conclude that an increase of glacial erosion in response to Quaternary cooling is physically plausible, and we show that the magnitude of the increase depends on rock uplift and ice accumulation rates.
In the Pan-African belt of the Dronning Maud Land, Antarctica, crystallized melt inclusions (nanogranitoids) occur in garnet from ultramafic granulites. The granulites contain the peak assemblage pargasite+garnet+clinopyroxene with rare relict orthopyroxene and biotite, and retrograde symplectites at contacts between garnet and amphibole. Garnet contains two generations of melt inclusions. Type 1 inclusions, interpreted as primary, are isolated, < 10 mu m in size, and generally have negative crystal shapes. They contain kokchetavite, kumdykolite, and phlogopite, with quartz and zoisite as minor phases, and undevitrified glass was identified in one inclusion. Type 2 inclusions are < 30 mu m in size, secondary, and contain amphibole, feldspars, and zoisite. Type 2 inclusions appear to be the crystallization products of a melt that coexisted with an immiscible CO2-rich fluid. The nanogranitoids were re-homogenized after heating in a piston-cylinder in a series of four experiments to investigate their composition. The conditions ranged between 900 and 950 degrees C at 1.5-2.4 GPa. Type 1 inclusions are trachytic and ultrapotassic, whereas type 2 melts are dacitic to rhyolitic. Thermodynamic modeling of the ultramafic composition in the MnNCKFMASHTO system shows that anatexis occurred at the end of the prograde P-T path, between the solidus (at ca. 860 degrees C-1.4 GPa) and the peak conditions (at ca. 960 degrees C-1.7 GPa). The model melt composition is felsic and similar to that of type 1 inclusions, particularly when the melting degree is low (< 1 mol%), close to the solidus. However the modeling fails to reproduce the highly potassic signature of the melt and its low H2O content. The combination of petrology, melt inclusion study, and thermodynamic modeling supports the interpretation that melt was produced by anatexis of the ultramafic boudins near peak P-T conditions, and that type 1 inclusions contain the anatectic melt that was present during garnet growth. The felsic, ultrapotassic composition of the primary anatectic melts is compatible with low melting degrees in the presence of biotite and amphibole as reactants.
Recent global warming is pronounced in high-latitude regions (e.g. northern Asia), and will cause the vegetation to change. Future vegetation trends (e.g. the "arctic greening") will feed back into atmospheric circulation and the global climate system. Understanding the nature and causes of past vegetation changes is important for predicting the composition and distribution of future vegetation communities. Fossil pollen records from 468 sites in northern and eastern Asia were biomised at selected times between 40 cal ka bp and today. Biomes were also simulated using a climate-driven biome model and results from the two approaches compared in order to help understand the mechanisms behind the observed vegetation changes. The consistent biome results inferred by both approaches reveal that long-term and broad-scale vegetation patterns reflect global- to hemispheric-scale climate changes. Forest biomes increase around the beginning of the late deglaciation, become more widespread during the early and middle Holocene, and decrease in the late Holocene in fringe areas of the Asian Summer Monsoon. At the southern and southwestern margins of the taiga, forest increases in the early Holocene and shows notable species succession, which may have been caused by winter warming at ca. 7 cal ka bp. At the northeastern taiga margin (central Yakutia and northeastern Siberia), shrub expansion during the last deglaciation appears to prevent the permafrost from thawing and hinders the northward expansion of evergreen needle-leaved species until ca. 7 cal ka bp. The vegetation-climate disequilibrium during the early Holocene in the taiga-tundra transition zone suggests that projected climate warming will not cause a northward expansion of evergreen needle-leaved species.
Tracking changes in biodiversity through time requires an understanding of the relationship between modern diversity and how this diversity is preserved in the fossil record. Fossil pollen is one way in which past vegetation diversity can be reconstructed. However, there is limited understanding of modern pollen-vegetation diversity relationships from biodiverse tropical ecosystems. Here, pollen (palynological) richness and diversity (Hill N (1)) are compared with vegetation richness and diversity from forest and savannah ecosystems in the New World and Old World tropics (Neotropics and Palaeotropics). Modern pollen data were obtained from artificial pollen traps deployed in 1-ha vegetation study plots from which vegetation inventories had been completed in Bolivia and Ghana. Pollen counts were obtained from 15 to 22 traps per plot, and aggregated pollen sums for each plot were > 2,500. The palynological richness/diversity values from the Neotropics were moist evergreen forest = 86/6.8, semi-deciduous dry forest = 111/21.9, wooded savannah = 138/31.5, and from the Palaeotropics wet evergreen forest = 144/28.3, semi-deciduous moist forest = 104/4.4, forest-savannah transition = 121/14.1; the corresponding vegetation richness/diversity was 100/36.7, 80/38.7 and 71/39.4 (Neotropics), and 101/54.8, 87/45.5 and 71/34.5 (Palaeotropics). No consistent relationship was found between palynological richness/diversity, and plot vegetation richness/diversity, due to the differential influence of other factors such as landscape diversity, pollination strategy, and pollen source area. Palynological richness exceeded vegetation richness, while pollen diversity was lower than vegetation diversity. The relatively high global diversity of tropical vegetation was found to be reflected in the pollen rain.
Global patterns of declining temperature variability from the Last Glacial Maximum to the Holocene
(2018)
Changes in climate variability are as important for society to address as are changes in mean climate(1). Contrasting temperature variability during the Last Glacial Maximum and the Holocene can provide insights into the relationship between the mean state of the climate and its variability(2,3). However, although glacial-interglacial changes in variability have been quantified for Greenland(2), a global view remains elusive. Here we use a network of marine and terrestrial temperature proxies to show that temperature variability decreased globally by a factor of four as the climate warmed by 3-8 degrees Celsius from the Last Glacial Maximum (around 21,000 years ago) to the Holocene epoch (the past 11,500 years). This decrease had a clear zonal pattern, with little change in the tropics (by a factor of only 1.6-2.8) and greater change in the mid-latitudes of both hemispheres (by a factor of 3.3-14). By contrast, Greenland ice-core records show a reduction in temperature variability by a factor of 73, suggesting influences beyond local temperature or a decoupling of atmospheric and global surface temperature variability for Greenland. The overall pattern of reduced variability can be explained by changes in the meridional temperature gradient, a mechanism that points to further decreases in temperature variability in a warmer future.
Sub-seasonal thaw slump mass wasting is not consistently energy limited at the landscape scale
(2018)
Predicting future thaw slump activity requires a sound understanding of the atmospheric drivers and geomorphic controls on mass wasting across a range of timescales. On sub-seasonal timescales, sparse measurements indicate that mass wasting at active slumps is often limited by the energy available for melting ground ice, but other factors such as rainfall or the formation of an insulating veneer may also be relevant. To study the sub-seasonal drivers, we derive topographic changes from single-pass radar interferometric data acquired by the TanDEM-X satellites. The estimated elevation changes at 12m resolution complement the commonly observed planimetric retreat rates by providing information on volume losses. Their high vertical precision (around 30 cm), frequent observations (11 days) and large coverage (5000 km(2)) allow us to track mass wasting as drivers such as the available energy change during the summer of 2015 in two study regions. We find that thaw slumps in the Tuktoyaktuk coastlands, Canada, are not energy limited in June, as they undergo limited mass wasting (height loss of around 0 cm day 1) despite the ample available energy, suggesting the widespread presence of early season insulating snow or debris veneer. Later in summer, height losses generally increase (around 3 cm day 1), but they do so in distinct ways. For many slumps, mass wasting tracks the available energy, a temporal pattern that is also observed at coastal yedoma cliffs on the Bykovsky Peninsula, Russia. However, the other two common temporal trajectories are asynchronous with the available energy, as they track strong precipitation events or show a sudden speed-up in late August respectively. The observed temporal patterns are poorly related to slump characteristics like the headwall height. The contrasting temporal behaviour of nearby thaw slumps highlights the importance of complex local and temporally varying controls on mass wasting.
The Northern Central Iranian Micro-continent (CIM) represents Neotethys-related oceanic crust remnants, emplaced due to convergence between CIM and Eurasia plates during Eocene. Mafic and ultramafic units are exposed along the northern part of the CIM in the Sabzevar area. The geology and field relation of Sabzevar ophiolite indicate northward subduction of the Sabzevar basin. The average whole rock chemistry of mafic (gabbros) and ultramafic samples (lherzolite, harzburgite and dunite) is characterized by a range of MgO of 11.16-31.88, CaO 5.22-11.53 and Al2O3 2.77-14.57, respectively. Low LREE/HREE ratio of ultramafic samples is accompanied by enrichment of large ion lithophile elements (LILE) such as Sr, Pb and K. Mafic samples show two distinct groups with low and high LREE/HREE ratios. The spider diagram of mafic samples indicates enrichment in Sr, Pb and K and depletion in REE. Petrological and geochemical evidence and field relations show that the mafic rocks formed in a supra-subduction zone setting. Petrological studies reveal the role of fractional crystallization and assimilation effect by released fluids during subduction related generation of the Sabzevar mafic rocks. We suggest that the studied mafic rocks likely represent the basement of an initial island arc, which was generated in a supra-subduction zone setting within the Neotethys branch of the Sabzevar Ocean at the north of CIM. Copper, gold and chromite mineralizations are studied in relation to island arc setting and supra-subduction environment. Similarities in lithology, ophiolite age and mineralization between Sabzevar ophiolite and Bardaskan-Torbat Heydariyeh ophiolites testify for their separation due to rotation (or faulting) of the Central Iranian Micro-continent.
Arctic tundra ecosystems have experienced unprecedented change associated with climate warming over recent decades. Across the Pan-Arctic, vegetation productivity and surface greenness have trended positively over the period of satellite observation. However, since 2011 these trends have slowed considerably, showing signs of browning in many regions. It is unclear what factors are driving this change and which regions/landforms will be most sensitive to future browning. Here we provide evidence linking decadal patterns in arctic greening and browning with regional climate change and local permafrost-driven landscape heterogeneity. We analyzed the spatial variability of decadal-scale trends in surface greenness across the Arctic Coastal Plain of northern Alaska (similar to 60,000 km(2)) using the Landsat archive (1999-2014), in combination with novel 30 m classifications of polygonal tundra and regional watersheds, finding landscape heterogeneity and regional climate change to be the most important factors controlling historical greenness trends. Browning was linked to increased temperature and precipitation, with the exception of young landforms (developed following lake drainage), which will likely continue to green. Spatiotemporal model forecasting suggests carbon uptake potential to be reduced in response to warmer and/or wetter climatic conditions, potentially increasing the net loss of carbon to the atmosphere, at a greater degree than previously expected.
Microbial oxidation of lithospheric organic carbon in rapidly eroding tropical mountain soils
(2018)
Lithospheric organic carbon ("petrogenic"; OCpetro) is oxidized during exhumation and subsequent erosion of mountain ranges. This process is a considerable source of carbon dioxide (CO2) to the atmosphere over geologic time scales, but the mechanisms that govern oxidation rates in mountain landscapes are poorly constrained. We demonstrate that, on average, 67 +/- 11% of the OCpetro initially present in bedrock exhumed from the tropical, rapidly eroding Central Range of Taiwan is oxidized in soils, leading to CO2 emissions of 6.1 to 18.6 metric tons of carbon per square kilometer per year. The molecular and isotopic evolution of bulk OC and lipid biomarkers during soil formation reveals that OCpetro remineralization is microbially mediated. Rapid oxidation in mountain soils drives CO2 emission fluxes that increase with erosion rate, thereby counteracting CO2 drawdown by silicate weathering and biospheric OC burial.
Soil pollutants such as hydrocarbons can induce toxic effects in plants and associated arbuscular mycorrhizal fungi (AMF). This study was conducted to evaluate if the legume Lotus corniculatus and the grass Elymus trachycaulus and arbuscular mycorrhizal fungi could grow in two oil sands processing by-products after bitumen extraction from the oil sands in northern Alberta, Canada. Substrate treatments were coarse tailings sand (CTS), a mix of dry mature fine tailings (MFT) with CTS (1: 1) and Pleistocene sandy soil (hydrocarbon free); microbial treatments were without AMF, with AMF and AMF plus soil bacteria isolated from oil sands reclamation sites. Plant biomass, root morphology, leaf water content, shoot tissue phosphorus content and mycorrhizal colonization were evaluated. Both plant species had reduced growth in CTS and tailings mix relative to sandy soil. AMF frequency and intensity in roots of E. trachycaulus was not influenced by soil hydrocarbons; however, it decreased significantly over time in roots of L. corniculatus without bacteria in CTS. Mycorrhizal inoculation alone did not significantly improve plant growth in CTS and tailings mix; however, inoculation with mycorrhizae plus bacteria led to a significantly positive response of both plant species in CTS. Thus, combined inoculation with selected mycorrhizae and bacteria led to synergistic effects. Such combinations may be used in future to improve plant growth in reclamation of CTS and tailings mix.
The Atlantic meridional overturning circulation (AMOC)—a system of ocean currents in the North Atlantic—has a major impact on climate, yet its evolution during the industrial era is poorly known owing to a lack of direct current measurements. Here we provide evidence for a weakening of the AMOC by about 3 ± 1 sverdrups (around 15 per cent) since the mid-twentieth century. This weakening is revealed by a characteristic spatial and seasonal sea-surface temperature ‘fingerprint’—consisting of a pattern of cooling in the subpolar Atlantic Ocean and warming in the Gulf Stream region—and is calibrated through an ensemble of model simulations from the CMIP5 project. We find this fingerprint both in a high-resolution climate model in response to increasing atmospheric carbon dioxide concentrations, and in the temperature trends observed since the late nineteenth century. The pattern can be explained by a slowdown in the AMOC and reduced northward heat transport, as well as an associated northward shift of the Gulf Stream. Comparisons with recent direct measurements from the RAPID project and several other studies provide a consistent depiction of record-low AMOC values in recent years.
Salt marshes filter pollutants, protect coastlines against storm surges, and sequester carbon, yet are under threat from sea level rise and anthropogenic modification. The sustained existence of the salt marsh ecosystem depends on the topographic evolution of marsh platforms. Quantifying marsh platform topography is vital for improving the management of these valuable landscapes. The determination of platform boundaries currently relies on supervised classification methods requiring near-infrared data to detect vegetation, or demands labour-intensive field surveys and digitisation. We propose a novel, unsupervised method to reproducibly isolate salt marsh scarps and platforms from a digital elevation model (DEM), referred to as Topographic Identification of Platforms (TIP). Field observations and numerical models show that salt marshes mature into subhorizontal platforms delineated by subvertical scarps. Based on this premise, we identify scarps as lines of local maxima on a slope raster, then fill landmasses from the scarps upward, thus isolating mature marsh platforms. We test the TIP method using lidar-derived DEMs from six salt marshes in England with varying tidal ranges and geometries, for which topographic platforms were manually isolated from tidal flats. Agreement between manual and unsupervised classification exceeds 94% for DEM resolutions of 1 m, with all but one site maintaining an accuracy superior to 90% for resolutions up to 3 m. For resolutions of 1 m, platforms detected with the TIP method are comparable in surface area to digitised platforms and have similar elevation distributions. We also find that our method allows for the accurate detection of local block failures as small as 3 times the DEM resolution. Detailed inspection reveals that although tidal creeks were digitised as part of the marsh platform, unsupervised classification categorises them as part of the tidal flat, causing an increase in false negatives and overall platform perimeter. This suggests our method may benefit from combination with existing creek detection algorithms. Fallen blocks and high tidal flat portions, associated with potential pioneer zones, can also lead to differences between our method and supervised mapping. Although pioneer zones prove difficult to classify using a topographic method, we suggest that these transition areas should be considered when analysing erosion and accretion processes, particularly in the case of incipient marsh platforms. Ultimately, we have shown that unsupervised classification of marsh platforms from high-resolution topography is possible and sufficient to monitor and analyse topographic evolution.
Agricultural soil landscapes of hummocky ground moraines are characterized by 3D spatial patterns of soil types that result from profile modifications due to the combined effect of water and tillage erosion. We hypothesize that crops reflect such soil landscape patterns by increased or reduced plant and root growth. Root development may depend on the thickness and vertical sequence of soil horizons as well as on the structural development state of these horizons at different landscape positions. The hypotheses were tested using field data of the root density (RD) and the root lengths (RL) of winter wheat using the minirhizotron technique. We compared data from plots at the CarboZALF-D site (NE Germany) that are representing a non-eroded reference soil profile (Albic Luvisol) at a plateau position, a strongly eroded profile at steep slope (Calcaric Regosol), and a depositional profile at the footslope (Anocolluvic Regosol). At each of these plots, three Plexiglas access tubes were installed down to approx. 1.5 m soil depth. Root measurements were carried out during the growing season of winter wheat (September 2014-August 2015) on six dates. The root length density (RLD) and the root biomass density were derived from RD values assuming a mean specific root length of 100 m g(-1). Values of RD and RLD were highest for the Anocolluvic Regosol and lowest for the Calcaric Regosol. The maximum root penetration depth was lower in the Anocolluvic Regosol because of a relatively high and fluctuating water table at this landscape position. Results revealed positive relations between below-ground (root) and above-ground crop parameters (i.e., leaf area index, plant height, biomass, and yield) for the three soil types. Observed root densities and root lengths in soils at the three landscape positions corroborated the hypothesis that the root system was reflecting erosion-induced soil profile modifications. Soil landscape position dependent root growth should be considered when attempting to quantify landscape scale water and element balances as well as agricultural productivity.
Management intensity modifies soil properties, e.g., organic carbon (C-org) concentrations and soil pH with potential feedbacks on plant diversity. These changes might influence microbial P concentrations (P-mic) in soil representing an important component of the Pcycle. Our objectives were to elucidate whether abiotic and biotic variables controlling P-mic concentrations in soil are the same for forests and grasslands, and to assess the effect of region and management on P-mic concentrations in forest and grassland soils as mediated by the controlling variables. In three regions of Germany, Schwabische Alb, Hanich-Dun, and Schorfheide-Chorin, we studied forest and grassland plots (each n=150) differing in plant diversity and land-use intensity. In contrast to controls of microbial biomass carbon (C-mic), P-mic was strongly influenced by soil pH, which in turn affected phosphorus (P) availability and thus microbial Puptake in forest and grassland soils. Furthermore, P-mic concentrations in forest and grassland soils increased with increasing plant diversity. Using structural equation models, we could show that soil C-org is the profound driver of plant diversity effects on P-mic in grasslands. For both forest and grassland, we found regional differences in P-mic attributable to differing environmental conditions (pH, soil moisture). Forest management and tree species showed no effect on P-mic due to a lack of effects on controlling variables (e.g., C-org). We also did not find management effects in grassland soils which might be caused by either compensation of differently directed effects across sites or by legacy effects of former fertilization constraining the relevance of actual practices. We conclude that variables controlling P-mic or C-mic in soil differ in part and that regional differences in controlling variables are more important for P-mic in soil than those induced by management.
Many hydrological models have been calibrated and validated using hydrographs alone. Because streamflow integrates water fluxes in space, many distributed hydrological models tend to have multiple feasible descriptions of hydrological processes. This equifinality usually leads to substantial prediction uncertainty. In this study, additional constraintsnamely, the spatial patterns of long-term average evapotranspiration (ET), shallow groundwater level, and land cover changewere used to investigate the reduction of equifinality and prediction uncertainty in the Soil and Water Assessment Tool (SWAT) in the Wami River basin in Tanzania. The additional constraints were used in the set-up, parameter emulation and calibration of the SWAT model termed an improved hydrological model (IHM). The IHM was then compared with a classical hydrological model (CHM) that was also developed using the SWAT model but without additional constraints. In the calibration, the CHM used only the hydrograph, but the IHM used the hydrograph and the spatial pattern of long-term average ET as an additional constraint. The IHM produced a single, unique behavioural simulation, whereas the CHM produced many behavioural simulations that resulted in prediction uncertainty. The performance of the IHM with respect to the hydrograph was more consistent than that of the CHM, and the former clearly captured the mean behaviour of ET in the river basin. Therefore, we conclude that additional constraints substantially reduce equifinality and prediction uncertainty in a distributed hydrological model.
West Africa has been afflicted by droughts since the declining rains of the 1970s. Therefore, this study examines the characteristics of drought over the Niger River Basin (NRB), investigates the influence of the drought on the river flow, and projects the impacts of future climate change on drought. A combination of observation data and regional climate simulations of past (1986-2005) and future climates (2046-2065 and 2081-2100) were analyzed. The standardized precipitation index (SPI) and standardized precipitation and evapotranspiration index (SPEI) were used to characterize drought while the standardized runoff index (SRI) was used to quantify river flow. Results of the study show that the historical pattern of drought is consistent with previous studies over the Basin and most part of West Africa. RCA4 ensemble gives realistic simulations of the climatology of the Basin in the past climate. Generally, an increase in drought intensity and frequency are projected over NRB. The coupling between SRI and drought indices was very strong (P < 0.05). The dominant peaks can be classified into three distinct drought cycles with periods 1-2, 2-4, 4-8 years. These cycles may be associated with Quasi-Biennial Oscillation (QBO) and El-Nino Southern Oscillation (ENSO). River flow was highly sensitive to precipitation in the NRB and a 1-3 month lead time was found between drought indices and SRI. Under RCP4.5, changes in the SPEI drought frequency range from 1.8 (2046-2065) to 2.4 (2081-2100) month year(-1) while under RCP8.5, the change ranges from 2.2 (2046-2065) to 3.0 month year(-1) (2081-2100). Niger Middle sub-basin is likely to be mostly impacted in the future while the Upper Niger was projected to be least impacted. Results of this study may guide policymakers to evolve strategies to facilitate vulnerability assessment and adaptive capacity of the basin in order to minimize the negative impacts of climate change.
Tillage erosion on arable land is a very important process leading to a net downslope movement of soil and soil constitutes. Tillage erosion rates are commonly in the same order of magnitude as water erosion rates and can be even higher, especially under highly mechanized agricultural soil management. Despite its prevalence and magnitude, tillage erosion is still understudied compared to water erosion. The goal of this study was to bring together experts using different techniques to determine tillage erosion and use the different results to discuss and quantify uncertainties associated with tillage erosion measurements. The study was performed in northeastern Germany on a 10 m by 50 m plot with a mean slope of 8%. Tillage erosion was determined after two sequences of seven tillage operations. Two different micro-tracers (magnetic iron oxide mixed with soil and fluorescent sand) and one macro-tracer (passive radio-frequency identification transponders (RFIDs), size: 4 x 22 mm) were used to directly determine soil fluxes. Moreover, tillage induced changes in topography were measured for the entire plot with two different terrestrial laser scanners and an unmanned aerial system for structure from motion topography analysis. Based on these elevation differences, corresponding soil fluxes were calculated. The mean translocation distance of all techniques was 0.57 m per tillage pass, with a relatively wide range of mean soil translocation distances ranging from 039 to 0.72 m per pass. A benchmark technique could not be identified as all used techniques have individual error sources, which could not be quantified. However, the translocation distances of the macro-tracers used were consistently smaller than the translocation distances of the micro-tracers (mean difference = -26 +/- 12%), which questions the widely used assumption of non-selective soil transport via tillage operations. This study points out that tillage erosion measurements, carried out under almost optimal conditions, are subject to major uncertainties that are far from negligible. (C) 2018 Elsevier B.V. All rights reserved.
Sulfur is an important component in volcanic gases at the Earth surface but also present in the deep Earth in hydrothermal or magmatic fluids. Little is known about the evolution of such fluids during ascent in the crust. A new optical cell was developed for in situ Raman spectroscopic investigations on fluids allowing abrupt or continuous changes of pressure up to 200 MPa at temperatures up to 750 degrees C. The concept is based on a flexible gold bellow, which separates the sample fluid from the pressure medium water. To avoid reactions between aggressive fluids and the pressure cell, steel components in contact with the fluid are shielded by gold foil. The cell was tested to study redox reactions in fluids using aqueous ammonium sulfate solutions as a model system. During heating at constant pressure of 130 MPa, sulfate ions transform first to HSO4- ions and then to molecular units such as H2SO4. Variation of pressure shows that the stability of sulfate species relies on fluid density, i.e., highly charged species are stable only in high-density fluids. Partial decomposition of ammonium was evident above 550 degrees C by the occurrence of a nitrogen peak in the Raman spectra. Reduced sulfur species were observed above 700 degrees C by Raman signals near 2590 cm(-1) assigned to HS- and H2S. No clear evidence for the formation of sulfur dioxide was found in contrary to previous studies on aqueous H2SO4, suggesting very reducing conditions in our experiments. Fluid-mineral interaction was studied by inserting into the cell a small, semi-open capsule filled with a mixture of pyrite and pyrrhotite. Oxidation of the sample assembly was evident by transformation of pyrite to pyrrhotite. As a consequence, sulfide species were observed in the fluid already at temperatures of similar to 600 degrees C.
An overview is given on the current state of X-ray absorption measurements on silicate melts and glasses. The challenges, limitations, and achievements of analyzing X-ray absorption spectra measured in liquids to determine structural properties of major and minor elements in magmas are described, with particular focus on describing non-Gaussian pair distribution functions in highly disordered glasses and melts, measured at in situ conditions. This includes a discussion on the progress of combining experiments with data from molecular dynamics simulations. For the measurements at conditions of the deep Earth, various experimental approaches and necessities are discussed and two examples are described in more detail. Finally, the achievements and prospects are presented for measuring X-ray absorption spectra indirectly by X-ray Raman scattering.
Crop yield variations are strongly influenced by the spatial and temporal availabilities of water and nitrogen in the soil during the crop growth season. To estimate the quantities and distributions of water and nitrogen within a given soil, process-oriented soil models have often been used. These models require detailed information about the soil characteristics and profile architecture (e.g., soil depth, clay content, bulk density, field capacity and wilting point), but high resolution information about these soil properties, both vertically and laterally, is difficult to obtain through conventional approaches. However, on-the-go electrical resistivity tomography (ERT) measurements of the soil and data inversion tools have recently improved the lateral resolutions of the vertically distributed measurable information. Using these techniques, nearly 19,000 virtual soil profiles with defined layer depths were successfully created for a 30 ha silty cropped soil over loamy and sandy substrates in Central Germany, which were used to initialise the CArbon and Nitrogen DYnamics (CANDY) model. The soil clay content was derived from the electrical resistivity (ER) and the collected soil samples using a simple linear regression approach (the mean R-2 of clay = 0.39). The additional required structural and hydrological properties were derived from pedotransfer functions. The modelling results, derived soil texture distributions and original ER data were compared with the spatial winter wheat yield distribution in a relatively dry year using regression and boundary line analysis. The yield variation was best explained by the simulated soil water content (R-2 = 0.18) during the grain filling and was additionally validated by the measured soil water content with a root mean square error (RMSE) of 7.5 Vol%.
Ice-rich yedoma-dominated landscapes store considerable amounts of organic carbon (C) and nitrogen (N) and are vulnerable to degradation under climate warming. We investigate the C and N pools in two thermokarst-affected yedoma landscapes - on Sobo-Sise Island and on Bykovsky Peninsula in the north of eastern Siberia. Soil cores up to 3m depth were collected along geomorphic gradients and analysed for organic C and N contents. A high vertical sampling density in the profiles allowed the calculation of C and N stocks for short soil column intervals and enhanced understanding of within-core parameter variability. Profile-level C and N stocks were scaled to the landscape level based on landform classifications from 5 m resolution, multispectral RapidEye satellite imagery. Mean landscape C and N storage in the first metre of soil for Sobo-Sise Island is estimated to be 20.2 kg C m(-2) and 1.8 kg N m(-2) and for Bykovsky Peninsula 25.9 kg C m(-2) and 2.2 kg N m(-2). Radiocarbon dating demonstrates the Holocene age of thermokarst basin deposits but also suggests the presence of thick Holoceneage cover layers which can reach up to 2 m on top of intact yedoma landforms. Reconstructed sedimentation rates of 0.10-0.57 mm yr(-1) suggest sustained mineral soil accumulation across all investigated landforms. Both yedoma and thermokarst landforms are characterized by limited accumulation of organic soil layers (peat). We further estimate that an active layer deepening of about 100 cm will increase organic C availability in a seasonally thawed state in the two study areas by similar to 5.8 Tg (13.2 kg C m(-2)). Our study demonstrates the importance of increasing the number of C and N storage inventories in ice-rich yedoma and thermokarst environments in order to account for high variability of permafrost and thermokarst environments in pan-permafrost soil C and N pool estimates.
Although hydrologic models provide hypothesis testing of complex dynamics occurring at catchments, fresh-water quality modeling is still incipient at many subtropical headwaters. In Brazil, a few modeling studies assess freshwater nutrients, limiting policies on hydrologic ecosystem services. This paper aims to compare freshwater quality scenarios under different land-use and land-cover (LULC) change, one of them related to ecosystem-based adaptation (EbA), in Brazilian headwaters. Using the spatially semi-distributed Soil and Water Assessment Tool (SWAT) model, nitrate, total phosphorous (TP) and sediment were modeled in catchments ranging from 7.2 to 1037 km(2). These head-waters were eligible areas of the Brazilian payment for ecosystem services (PES) projects in the Cantareira water supply system, which had supplied water to 9 million people in the Sao Paulo metropolitan region (SPMR). We considered SWAT modeling of three LULC scenarios: (i) recent past scenario (S1), with historical LULC in 1990; (ii) current land-use scenario (S2), with LULC for the period 2010-2015 with field validation; and (iii) future land-use scenario with PES (S2 + EbA). This latter scenario proposed forest cover restoration through EbA following the river basin plan by 2035. These three LULC scenarios were tested with a selected record of rainfall and evapotranspiration observed in 2006-2014, with the occurrence of extreme droughts. To assess hydrologic services, we proposed the hydrologic service index (HSI), as a new composite metric comparing water pollution levels (WPL) for reference catchments, related to the grey water footprint (greyWF) and water yield. On the one hand, water quality simulations allowed for the regionalization of greyWF at spatial scales under LULC scenarios. According to the critical threshold, HSI identified areas as less or more sustainable catchments. On the other hand, conservation practices simulated through the S2 + EbA scenario envisaged not only additional and viable best management practices (BMP), but also preventive decision-making at the headwaters of water supply systems.
Understanding the variability of the ionosphere is important for the prediction of space weather and climate. Recent studies have shown that forcing from the lower atmosphere plays a significant role for the short-term (day-to-day) variability of the low-latitude ionosphere. The present study aims to assess the importance of atmospheric forcing for the variability of the daytime equatorial ionospheric electric field on the interannual (year-to-year) time scale. Magnetic field measurements from Huancayo (12.05 degrees S, 75.33 degrees W) are used to augment the equatorial vertical plasma drift velocity (V-Z) measurements from the Jicamarca Unattended Long-term Investigations of the Ionosphere and Atmosphere radar during 2001-2016. V-Z can be regarded as a measure of the zonal electric field. After removing the seasonal variation of similar to 10m/s, midday values of V-Z show an interannual variation of similar to 2m/s with an oscillation period of 2-3years. No evidence of solar cycle influence is found. The Ground-to-topside Atmosphere-Ionosphere model for Aeronomy, which takes into account realistic atmospheric variability below 30km, reproduces the pattern of the observed interannual variation without having to include variable forcing from the magnetosphere. The results indicate that lower atmospheric forcing plays a dominant role for the observed interannual variability of V-Z at 1200 local time.
Germination marks a critical transition in plant life that is prone to high mortality. Strong selection pressure is therefore expected to finely tune it to environmental conditions. Our study on the common Mediterranean grass Brachypodium hybridum assessed whether germination behavior changes systematically along a steep natural rainfall gradient ranging from harsh desert to rather mild mesic-Mediterranean conditions, We specifically tested hypotheses that germination behavior confers greater risk-spreading in populations from drier, unpredictable environments, and that seeds from wetter populations are better competitors. In 14 populations (spanning 114-954 mm annual rainfall) we assessed three alternative key parameters of germination in a greenhouse experiment: between-year dormancy, days to emergence within a season, and temporal spread. Addition of neighbor seeds accounted for competition as another crucial environmental factor. In six of the 14 populations, we also compared seeds originating from corresponding north (more mesic) and south (more arid) exposed hill slopes to test whether germination patterns along the large-scale rainfall gradient are paralleled at this smaller scale. B. hybridum exhibited generally high germination fractions and rapid emergence with very little temporal spread, indicating overall little risk-spreading germination. Surprisingly, none of the three parameters changed systematically with increasing aridity, neither at large scale along the rainfall gradient nor at small scale between north and south exposures. Neighbor seeds, however, mildly suppressed germination. Germination of neighbor seeds, in turn, was more strongly suppressed by B. hybridum seeds from drier populations, and this effect was stronger for forb than for grass neighbor species. Our results provide strong evidence that increased risk-spreading germination is not a universal, essential strategy to persist in increasingly dry, unpredictable environments. They also highlight that competition with neighbors occurs even at the earliest plant life stage. Since neighbor effects were species-specific, competition among seeds can affect community composition at later plant stages.
Effects of upper mantle heterogeneities on the lithospheric stress field and dynamic topography
(2018)
The orientation and tectonic regime of the observed crustal/lithospheric stress field contribute to our knowledge of different deformation processes occurring within the Earth's crust and lithosphere. In this study, we analyze the influence of the thermal and density structure of the upper mantle on the lithospheric stress field and topography. We use a 3-D lithosphere–asthenosphere numerical model with power-law rheology, coupled to a spectral mantle flow code at 300 km depth. Our results are validated against the World Stress Map 2016 (WSM2016) and the observation-based residual topography. We derive the upper mantle thermal structure from either a heat flow model combined with a seafloor age model (TM1) or a global S-wave velocity model (TM2). We show that lateral density heterogeneities in the upper 300 km have a limited influence on the modeled horizontal stress field as opposed to the resulting dynamic topography that appears more sensitive to such heterogeneities. The modeled stress field directions, using only the mantle heterogeneities below 300 km, are not perturbed much when the effects of lithosphere and crust above 300 km are added. In contrast, modeled stress magnitudes and dynamic topography are to a greater extent controlled by the upper mantle density structure. After correction for the chemical depletion of continents, the TM2 model leads to a much better fit with the observed residual topography giving a good correlation of 0.51 in continents, but this correction leads to no significant improvement of the fit between the WSM2016 and the resulting lithosphere stresses. In continental regions with abundant heat flow data, TM1 results in relatively small angular misfits. For example, in western Europe the misfit between the modeled and observation-based stress is 18.3°. Our findings emphasize that the relative contributions coming from shallow and deep mantle dynamic forces are quite different for the lithospheric stress field and dynamic topography.
Lithospheric plates move over the low-viscosity asthenosphere balancing several forces, which generate plate motions. We use a global 3-D lithosphere-asthenosphere model (SLIM3D) with visco-elasto-plastic rheology coupled to a spectral model of mantle flow at 300 km depth to quantify the influence of intra-plate friction and asthenospheric viscosity on plate velocities. We account for the brittle-ductile deformation at plate boundaries (yield stress) using a plate boundary friction coefficient to predict the present-day plate motion and net rotation of the lithospheric plates. Previous modeling studies have suggested that small friction coefficients (mu < 0.1, yield stress similar to 100 MPa) can lead to plate tectonics in models of mantle convection. Here we show that in order to match the observed present-day plate motion and net rotation, the frictional parameter must be less than 0.05. We obtain a good fit with the magnitude and orientation of the observed plate velocities (NUVEL-1A) in a no-net-rotation (NNR) reference frame with mu < 0.05 and a minimum asthenosphere viscosity of similar to 5 . 10(19) Pas to 10(20) Pas. Our estimates of net rotation (NR) of the lith-osphere suggest that amplitudes similar to 0.1-0.2 (degrees/Ma), similar to most observation-based estimates, can be obtained with asthenosphere viscosity cutoff values of similar to 10(19) Pas to 5 . 10(19) Pas and friction coefficients mu < 0.05.
Landslide hazard motivates the need for a deeper understanding of the events that occur before, during, and after catastrophic slope failures. Due to the destructive nature of such events, in situ observation is often difficult or impossible. Here, we use data from a network of 58 seismic stations to characterise a large landslide at the Askja caldera, Iceland, on 21 July 2014. High data quality and extensive network coverage allow us to analyse both long- and short-period signals associated with the landslide, and thereby obtain information about its triggering, initiation, timing, and propagation. At long periods, a landslide force history inversion shows that the Askja landslide was a single, large event starting at the SE corner of the caldera lake at 23:24:05 UTC and propagating to the NW in the following 2 min The bulk sliding mass was 7-16 x 10(10) kg, equivalent to a collapsed volume of 35-80 x 10(6) m(3). The sliding mass was displaced downslope by 1260 +/- 250 m. At short periods, a seismic tremor was observed for 30 min before the landslide. The tremor is approximately harmonic with a fundamental frequency of 2.3 Hz and shows time-dependent changes of its frequency content. We attribute the seismic tremor to stick-slip motion along the landslide failure plane. Accelerating motion leading up to the catastrophic slope failure culminated in an aseismic quiescent period for 2 min before the landslide. We propose that precursory seismic signals may be useful in landslide early-warning systems. The 8 h after the main landslide failure are characterised by smaller slope failures originating from the destabilised caldera wall decaying in frequency and magnitude. We introduce the term "afterslides" for this subsequent, declining slope activity after a large landslide.
The impact of reservoir heterogeneities on High-Temperature aquifer thermal energy storage systems
(2018)
We conducted a geoscientific feasibility study for the development of a high-temperature thermal aquifer energy storage system (HT-ATES) outside the capital of Muscat, northern Oman. The aquifer storage is part of a solar geothermal cooling project for the sustainable and continuous cooling of office buildings. The main concept is that excess solar energy will be stored in the subsurface through hot water injection and subsequently utilised as auxiliary energy source during peak demand times. The characterisation of aquifer heterogeneities is thus essential to predict subsurface thermal heat plume development and recovery efficiency of the storage system. We considered two aquifer systems as potential storage horizons, (i) a clastic-dominated alluvial fan system where individual channel systems in combination with diagenetic alterations constitute the main heterogeneities and (ii) a carbonate-dominated system represented by a homogenous layer-cake architecture. The feasibility study included a multidisciplinary approach from initial field work, geocellular reservoir modelling to finite element fluid flow and thermal modelling. Our results show that for the HT-ATES system, with a high frequency of injection and production cycles, heat loss mainly occurs due to heterogeneities in the permeability field of the aquifer in combination with buoyancy driven vertical fluid flow. An impermeable cap-rock is needed to keep the heat plume in place. Conductive heat loss is a minor issue. Highly complex heat plume geometries are apparent in the clastic channel system and ATES well planning is challenging due to the complex and interconnected high permeable channels. The carbonate sequence shows uniform plume geometries due to the layer cake architecture of the system and is tentatively more suitable for ATES development. Based on our findings we propose the general concept of HT-ATES traps, incorporating and building on expertise and knowledge from petroleum and reservoir geology regarding reservoir rocks and suitable trap&seal geometries. The concept can be used as guideline for future high-temperature aquifer storage exploration and development.
With increasing amount of strong motion data, Ground Motion Prediction Equation (GMPE) developers are able to quantify empirical site amplification functions (delta S2S(s)) from GMPE residuals, for use in site-specific Probabilistic Seismic Hazard Assessment. In this study, we first derive a GMPE for 5% damped Pseudo Spectral Acceleration (g) of Active Shallow Crustal earthquakes in Japan with 3.4 <= M-w <= 7.3 and 0 <= R-JB <= 600km. Using k-mean spectral clustering technique, we then classify our estimated delta S2S(s)(T = 0.01 - 2s) of 588 wellcharacterized sites, into 8 site clusters with distinct mean site amplification functions, and within-cluster site-tosite variability similar to 50% smaller than the overall dataset variability (phi(S2S)). Following an evaluation of existing schemes, we propose a revised data-driven site classification characterized by kernel density distributions of V-s30, V-s10, H-800, and predominant period (T-G) of the site clusters.
For over a century, geomorphologists have attempted to unravel information about landscape evolution, and processes that drive it, using river profiles. Many studies have combined new topographic datasets with theoretical models of channel incision to infer erosion rates, identify rock types with different resistance to erosion, and detect potential regions of tectonic activity. The most common metric used to analyse river profile geometry is channel steepness, or k(s). However, the calculation of channel steepness requires the normalisation of channel gradient by drainage area. This normalisation requires a power law exponent that is referred to as the channel concavity index. Despite the concavity index being crucial in determining channel steepness, it is challenging to constrain. In this contribution, we compare both slope-area methods for calculating the concavity index and methods based on integrating drainage area along the length of the channel, using so-called "chi" (chi) analysis. We present a new chi-based method which directly compares chi values of tributary nodes to those on the main stem; this method allows us to constrain the concavity index in transient landscapes without assuming a linear relationship between chi and elevation. Patterns of the concavity index have been linked to the ratio of the area and slope exponents of the stream power incision model (m/n); we therefore construct simple numerical models obeying detachment-limited stream power and test the different methods against simulations with imposed m and n. We find that chi-based methods are better than slope-area methods at reproducing imposed m/n ratios when our numerical landscapes are subject to either transient uplift or spatially varying uplift and fluvial erodibility. We also test our methods on several real landscapes, including sites with both lithological and structural heterogeneity, to provide examples of the methods' performance and limitations. These methods are made available in a new software package so that other workers can explore how the concavity index varies across diverse landscapes, with the aim to improve our understanding of the physics behind bedrock channel incision.
The use of chironomids (Diptera: Chironomidae) as indicators of changes in natural and climatic settings allows one to reconstruct paleoclimatic and paleoenvironmental conditions. Based on the results of lithological and chironomid analyses, the loss on ignition (LOI) values, as well as the results of radiocarbon dating (C-14 AMS), new data on the paleoclimate in the Karelian Isthmus, in particular the quantitative reconstruction of the mean July temperature (T (July), A degrees C) in the Late Pleistocene and Holocene and the evolution of the ecosystem of Medvedevskoe Lake have been obtained.
Quantifying the pace of ice-sheet growth is critical to understanding ice-age climate and dynamics. Here, we show that the diversion of the Hudson River (northeastern North America) late in the last glaciation phase (ca. 30 ka), which some previous studies have speculated was due to glacial isostatic adjustment (GIA), can be used to infer the timing of the Laurentide Ice Sheet’s growth to its maximum extent. Landscapes in the vicinity of glaciated regions have likely responded to crustal deformation produced by ice-sheet growth and decay through river drainage reorganization, given that rates of uplift and subsidence are on the order of tens of meters per thousand years. We perform global, gravitationally self-consistent simulations of GIA and input the predicted crustal deformation field into a landscape evolution model. Our calculations indicate that the eastward diversion of the Hudson River at 30 ka is consistent with exceptionally rapid growth of the Laurentide Ice Sheet late in the glaciation phase, beginning at 50–35 ka.
Glaciers and rivers are the main agents of mountain erosion. While in the fluvial realm empirical relationships and their mathematical description, such as the stream power law, improved the understanding of fundamental controls on landscape evolution, simple constraints on glacial topography and governing scaling relations are widely lacking. We present a steady state solution for longitudinal profiles along eroding glaciers in a coupled system that includes tectonics and climate. We combined the shallow ice approximation and a glacial erosion rule to calculate ice surface and bed topography from prescribed glacier mass balance gradient and rock uplift rate. Our approach is inspired by the classic application of the stream power law for describing a fluvial steady state but with the striking difference that, in the glacial realm, glacier mass balance is added as an altitude-dependent variable. From our analyses we find that ice surface slope and glacial relief scale with uplift rate with scaling exponents indicating that glacial relief is less sensitive to uplift rate than relief in most fluvial landscapes. Basic scaling relations controlled by either basal sliding or internal deformation follow a power law with the exponent depending on the exponents for the glacial erosion rule and Glen's flow law. In a mixed scenario of sliding and deformation, complicated scaling relations with variable exponents emerge. Furthermore, a cutoff in glacier mass balance or cold ice in high elevations can lead to substantially larger scaling exponents which may provide an explanation for high relief in high latitudes.
In areas characterized by many volcanoes, identifying the source of a deposit may not be trivial. This becomes much more complicated when looking for the source of a debris avalanche deposit (DAD), the common products of catastrophic volcanic edifice collapses. To overcome this problem, in this work a methodology is proposed based on the integration of texture features and areal distribution of the deposit, comparison between the petrography of the coarser clasts within the DAD and of the proximal products, grain-size analysis, and the volumetric estimations of the deposit and the volume missing from the volcanic edifice. This methodology has been tested to a DAD occurred near the city of Ibarra (Imbabura Province; Northern Ecuador), having a controversial source. Two main volcanic edifice are located in proximity of the DAD, the Cubilche volcano (3826 m.a.s.l.), located immediately south of and east of the colossal dormant Imbabura volcano. The former displays a sharp horseshoe shaped scar towards the north and inside this post-collapse edifice, that we name old Cubilche volcano (OCV), is located the young Cubilche volcano (YCV) that refilled a portion of the collapse scar and partially covered the southern flank of the OCV. Detailed knowledge of Cubilche volcano is critical because of its close proximity and interspersed activity with Imbabura volcano. In fact, Imbabura most recent edifice was built over the northwestern slope of the OCV and partially covered it. Recent studies linked the studied DAD to both Imbabura volcano as a product of its northern sector collapse, as well as neighboring Cubilche volcano. Our data points to Cubilche as the most likely source for this DAD. A perspective view of the shaded relief image of the present day OCV shows that the morphology of the volcano is well-preserved on its southern, eastern, and western flanks. This allows us to reconstruct the morphology of the OCV previous to the collapse through interpolation of elevation and altitude data of preserved flanks. A DEM of the present day topography was used for extrapolating the morphology. Using similar methodology, the post collapse base of the amphitheater was reconstructed by removing the relief of the present day YCV. The reconstructed topography of the OCV shows that it could have been a symmetric cone, reaching a maximum elevation of similar to 4100 m.a.s.l. with a lack volume of similar to 3.5 km(3). Based on this scenario, the deposit originated from the OCV main collapse should have a volume >3-3.5 km(3) in accordance to the volume calculated for the studied DAD. (C) 2018 Elsevier B.V. All rights reserved.
Climate and environmental changes since the Last Glacial Maximum in the tropical zone of West Africa are usually inferred from marine and continental records. In this study, the potential of carbonate pedo-sedimentary geosystems, i.e. Vertisol relics, to record paleoenvironmental changes in the southwestern part of Chad Basin are investigated. A multi-dating approach was applied on different pedogenic organo-mineral constituents. Optically stimulated luminescence (OSL) dating was performed on the soil K-rich feldspars and was combined with radiocarbon dating on both the inorganic (C-14(inorg)) and organic carbon (C-14(org)) soil fractions. Three main pedo-sedimentary processes were assessed over the last 20 ka BP: 1) the soil parent material deposition, from 18 ka to 12 ka BP (OSL), 2) the soil organic matter integration, from 11 cal ka to 8 cal ka BP (C-14(org)), and 3) the pedogenic carbonate nodule precipitation, from 7 cal ka to 5 cal ka BP (C-14(inorg)). These processes correlate well with the Chad Basin stratigraphy and West African records and are shown to be related to significant changes in the soil water balance responding to the evolution of continental hydrology during the Late Quaternary. The last phase affecting the Vertisol relics is the increase of erosion, which is hypothesized to be due to a decrease of the vegetation cover triggered by (i) the onset of drier conditions, possibly strengthened by (ii) anthropogenic pressure. Archaeological data from Far North Cameroon and northern Nigeria, as well as sedimentation times in Lake Tilla (northeastern Nigeria), were used to test these relationships. The increase of erosion is suggested to possibly occur between c. 3 cal ka and 1 cal ka BP. Finally, satellite images revealed similar geosystems all along the Sudano-Sahelian belt, and initial C-14(inorg) ages of the samples collected in four sites gave similar ages to those reported in this study. Consequently, the carbonate pedo-sedimentary geosystems are valuable continental paleoenvironmental archives and soil water balance proxies of the semiarid tropics of West Africa. (C) 2018 Elsevier Ltd. All rights reserved.
The literature contains a sizable number of publications where weather types are used to decompose climate shifts or trends into contributions of frequency and mean of those types. They are all based on the product rule, that is, a transformation of a product of sums into a sum of products, the latter providing the decomposition. While there is nothing to argue about the transformation itself, its interpretation as a climate shift or trend decomposition is bound to fail. While the case of a climate shift may be viewed as an incomplete description of a more complex behaviour, trend decomposition indeed produces bogus trends, as demonstrated by a synthetic counterexample with well-defined trends in type frequency and mean. Consequently, decompositions based on that transformation, be it for climate shifts or trends, must not be used.
It is well-known that prey species often face trade-offs between defense against predation and competitiveness, enabling predator-mediated coexistence. However, we lack an understanding of how the large variety of different defense traits with different competition costs affects coexistence and population dynamics. Our study focusses on two general defense mechanisms, that is, pre-attack (e.g., camouflage) and post-attack defenses (e.g., weaponry) that act at different phases of the predator—prey interaction. We consider a food web model with one predator, two prey types and one resource. One prey type is undefended, while the other one is pre- or post-attack defended paying costs either by a higher half-saturation constant for resource uptake or a lower maximum growth rate. We show that post-attack defenses promote prey coexistence and stabilize the population dynamics more strongly than pre-attack defenses by interfering with the predator's functional response: Because the predator spends time handling “noncrackable” prey, the undefended prey is indirectly facilitated. A high half-saturation constant as defense costs promotes coexistence more and stabilizes the dynamics less than a low maximum growth rate. The former imposes high costs at low resource concentrations but allows for temporally high growth rates at predator-induced resource peaks preventing the extinction of the defended prey. We evaluate the effects of the different defense mechanisms and costs on coexistence under different enrichment levels in order to vary the importance of bottom-up and top-down control of the prey community.
Ground‐motion prediction equations (GMPEs) are calibrated to predict the intensity of ground shaking at any given location, based on earthquake magnitude, source‐to‐site distance, local soil amplifications, and other parameters. GMPEs are generally assumed to be independent of time; however, evidence is increasing that large earthquakes modify the shallow soil conditions and those of the fault zone for months or years. These changes may affect the intensity of shaking and result in time‐dependent effects that can potentially be resolved by analyzing between‐event residuals (residuals between observed and predicted ground motion for individual earthquakes averaged over all stations). Here, we analyze a data set of about 65,000 recordings for about 1400 earthquakes in the moment magnitude range 2.5–6.5 that occurred in central Italy from 2008 to 2017 to capture the temporal variability of the ground shaking at high frequency. We first compute between‐event residuals for each earthquake in the Fourier domain with respect to a GMPE developed ad hoc for the analyzed data set. The between‐events show large changes after the occurrence of mainshocks such as the 2009 Mw 6.3 L'Aquila, the 2016 Mw 6.2 Amatrice, and Mw 6.5 Norcia earthquakes. Within the time span of a few months after the mainshocks, the between‐event contribution to the ground shaking varies by a factor 7. In particular, we find a large drop in the between‐events in the aftermath of the L'Aquila earthquake, followed by a slow positive trend that leads to a recovery interrupted by a new drop at the beginning of 2014. We also quantify the frequency‐dependent correlation between the Brune stress drop Δσ and the between‐events. We find that the temporal changes of Δσ resemble those of the between‐event residuals; in particular, during the period when the between‐events show the positive trend, the average logarithm of Δσ increases with an annual rate of 0.19 (i.e., the amplification factor for Δσ is 1.56 per year). Breakpoint analysis located a change in the linear trend coefficients of Δσ versus time in February 2014, although no large earthquakes occurred at that time. Finally, the temporal variability of Δσ mirrors the relative seismic‐velocity variations observed in previous studies for the same area and period, suggesting that both crack healing along the main fault system and healing of microcracks distributed at shallow depths throughout the surrounding region might be necessary to explain the wider observations of postearthquake recovery.
Nitraria is a halophytic taxon (i.e., adapted to saline environments) that belongs to the plant family Nitrariaceae and is distributed from the Mediterranean, across Asia into the south-eastern tip of Australia. This taxon is thought to have originated in Asia during the Paleogene (66-23 Ma), alongside the proto-Paratethys epicontinental sea. The evolutionary history of Nitraria might hold important clues on the links between climatic and biotic evolution but limited taxonomic documentation of this taxon has thus far hindered this line of research. Here we investigate if the pollen morphology and the chemical composition of the pollen wall are informative of the evolutionary history of Nitraria and could explain if origination along the proto-Paratethys and dispersal to the Tibetan Plateau was simultaneous or a secondary process. To answer these questions, we applied a novel approach consisting of a combination of Fourier Transform Infrared spectroscopy (FTIR), to determine the chemical composition of the pollen wall, and pollen morphological analyses using Light Microscopy (LM) and Scanning Electron Microscopy (SEM). We analysed our data using ordinations (principal components analysis and non-metric multidimensional scaling), and directly mapped it on the Nitrariaceae phylogeny to produce a phylomorphospace and a phylochemospace. Our LM, SEM and FTIR analyses show clear morphological and chemical differences between the sister groups Peganum and Nitraria. Differences in the morphological and chemical characteristics of highland species (Nitraria schoberi, N. sphaerocarpa, N. sibirica and N. tangutorum) and lowland species (Nitraria billardierei and N. retusa) are very subtle, with phylogenetic history appearing to be a more important control on Nitraria pollen than local environmental conditions. Our approach shows a compelling consistency between the chemical and morphological characteristics of the eight studied Nitrariaceae species, and these traits are in agreement with the phylogenetic tree. Taken together, this demonstrates how novel methods for studying fossil pollen can facilitate the evolutionary investigation of living and extinct taxa, and the environments they represent.
The impressive number of stream gauges in Chile, combined with a suite of past and recent large earthquakes, makes Chile a unique natural laboratory to study several streams that recorded responses to multiple seismic events. We document changes in discharge in eight streams in Chile following two or more large earthquakes. In all cases, discharge increases. Changes in discharge occur for peak ground velocities greater than about 7-11cm/s. Above that threshold, the magnitude of both the increase in discharge and the total excess water do not increase with increasing peak ground velocities. While these observations are consistent with previous work in California, they conflict with lab experiments that show that the magnitude of permeability changes increases with increasing amplitude of ground motion. Instead, our study suggests that streamflow responses are binary. Plain Language Summary Earthquakes deform and shake the surface and the ground below. These changes may affect groundwater flows by increasing the permeability along newly formed cracks and/or clearing clogged pores. As a result, groundwater flow may substantially increase after earthquakes and remain elevated for several months. Here we document streamflow anomalies following multiple high magnitude earthquakes in multiple streams in one of the most earthquake prone regions worldwide, Chile. We take advantage of the dense monitoring network in Chile that recorded streamflow since the 1940s. We show that once a critical ground motion is exceeded, streamflow responses to earthquakes can be expected.
A comprehensive understanding of the regional vegetation responses to long-term climate change will help to forecast Earth system dynamics. Based on a new well-dated pollen data set from Kanas Lake and a review on the published pollen records in and around the Altai Mountains, the regional vegetation dynamics and forcing mechanisms are discussed. In the Altai Mountains, the forest optimum occurred during 10-7ka for the upper forest zone and the tree line decline and/or ecological shifts were caused by climatic cooling from around 7ka. In the lower forest zone, the forest reached an optimum in the middle Holocene, and then increased openness of the forest, possibly caused by both climate cooling and human activities, took place in the late Holocene. In the lower basins or plains around the Altai Mountains, the development of protograssland or forest benefited from increasing humidity in the middle to late Holocene. Plain Language Summary In the Altai Mountains and surrounding area of central Asia, the previous studies of the Holocene paleovegetation and paleoclimate studies did not discuss the different ecological limiting factors for the vegetation in high mountains and low-elevation areas due to limited data. With accumulating fossil pollen data and surface pollen data, it is possible to understand better the geomorphological effect on the vegetation and discrepancies of vegetation/forest responses to large-scale climate forcing, and it is also possible to get reliable quantitative reconstructions of climate. Here our new pollen data and review on the published fossil pollen data will help us to look into the past climate change and vertical evolution of vegetation in this important area of the Northern Hemisphere. Based on our study, it can be concluded that the growth of taiga forest in the wetter areas may be promoted under a future warmer climate, while the forest in the relatively dry areas is liable to decline, and the different vegetation dynamics will contribute to future high-resolution coupled vegetation-climate model for Earth system modelling.