Refine
Has Fulltext
- yes (139)
Year of publication
Document Type
- Postprint (139) (remove)
Keywords
- climate change (6)
- remote sensing (5)
- permafrost (4)
- temperature (4)
- Holocene (3)
- climate-change (3)
- deep biosphere (3)
- dynamics (3)
- events (3)
- model (3)
- record (3)
- subsidence (3)
- time-series (3)
- Chinese loess (2)
- Lake Van (2)
- Sentinel-1 (2)
- UAV (2)
- Uncertainties (2)
- climate (2)
- climate extremes (2)
- correlation (2)
- damage (2)
- ecosystems (2)
- extreme rainfall (2)
- governance (2)
- hyperspectral (2)
- imaging spectroscopy (2)
- impacts (2)
- karst (2)
- lake sediments (2)
- modelling (2)
- models (2)
- playa (2)
- preparedness (2)
- runoff (2)
- salt pan (2)
- snow (2)
- streamflow (2)
- time series analysis (2)
- uncertainty (2)
- water (2)
- "Little Ice Age' (LIA) (1)
- "Medieval Warm Period' (MWP) (1)
- 3D printing (1)
- ALOS World 3D (1)
- ALOS/PALSAR (1)
- ASTER GDEM (1)
- Accuracy Asseessment (1)
- Accuracy Assessment (1)
- Acidithiobacillus (1)
- Acidobactetiaceae (1)
- Acidothermus (1)
- AgI (1)
- Air-pollution (1)
- Algorithm (1)
- Andes (1)
- Antarctic ice (1)
- Anthropogenic sources (1)
- Apatite (U-TH)/HE (1)
- Aral Sea (1)
- Arctic Ocean (1)
- Arctic Siberia (1)
- Arctic lakes (1)
- Arctic ocean (1)
- Arctic tundra ecosystems (1)
- Argentine margine (1)
- August 2002 flood (1)
- Baldeggersee (1)
- Band (1)
- Barents Sea (1)
- Bayesian classification (1)
- Bayesianism (1)
- Big Naryn complex (1)
- CDOM (1)
- CMIP5 models (1)
- CU (1)
- Carbo-Iron (1)
- Central Asia (1)
- Central Europe (1)
- Central andes (1)
- Central-Asia (1)
- China (1)
- Climate (1)
- Code_Aster (1)
- Connectivity (1)
- Copernicus DEM (1)
- Crustal density (1)
- DEM noise (1)
- DNA preservation (1)
- Digital Elevation Model (1)
- Digital Elevation Models (1)
- ERA5 (1)
- EROEI (1)
- EnGeoMAP 2.0 (1)
- EnMAP (1)
- Environmental sciences (1)
- Equatorial Pacific (1)
- Equatorial plasma irregularities (1)
- Error (1)
- European basin system (1)
- Expression (1)
- Extraction (1)
- Fen complex (1)
- Firmicutes (1)
- Fission-track thermochronology (1)
- Floods Directive (1)
- Fluid Flow (1)
- Forecasting Framework (1)
- Forschungsmethodik (1)
- Fourier analysis (1)
- GEDI (1)
- GLDAS (1)
- GNSS-integrated water vapour (1)
- GRACE (1)
- Gemmatimonadetes (1)
- Geomechanical Model (1)
- Germany (1)
- Glacial refugia (1)
- Glacier Mass Balances (1)
- Global Navigation Satellite System (GNSS) (1)
- Gravel-bed (1)
- Gravity (1)
- Greenland ice cores (1)
- Grundwassersanierung (1)
- Gulf of Mexico (1)
- HVSR (1)
- Halobacteria (1)
- Himalayas (1)
- HySpex (1)
- Hydrus-2D (1)
- Hyperion (1)
- ICESat-2 (1)
- IMPRESSIONS (1)
- InSAR (1)
- Ion-Beam (FIB) (1)
- Ionic liquids (1)
- Jaguaribe Basin (1)
- June 2013 flood (1)
- Kalahari (1)
- Kuilyu complex (1)
- Kwajalein (1)
- L-Asterisk (1)
- LGM (1)
- Lake Mead (1)
- Lake Towuti (1)
- Landsat (1)
- Landscape Response (1)
- Landslide inventory (1)
- Landslide susceptibility (1)
- Laptev Sea region (1)
- Larix larch (1)
- Lena delta (1)
- Liguride Complex (1)
- Lithospheric structure (1)
- Long-Term (1)
- MO (1)
- Methodology (1)
- Middle Eocene deformation (1)
- Middle Strand (1)
- NAO (1)
- NE Pamir (1)
- NEG (1)
- NO2 (1)
- NW Argentina (1)
- NW Turkey (1)
- Nanoeisen (1)
- Nearshore zone (1)
- Nevado Coropuna (1)
- Niger River (1)
- Nordeste (1)
- North Anatolian Fault (1)
- North Calabrian Unit (1)
- Northern Sierras Pampeanas (1)
- Northwestern Anatolia (1)
- Ocean color remote sensing (1)
- Ocean sciences (1)
- OpenForecast (1)
- PLA (1)
- Palaeoecology (1)
- Patterns (1)
- Pinus pinaster (1)
- Precipitation (1)
- Pull-Apart (1)
- Puna Plateau (1)
- Qilian mountains (1)
- Quercus ilex (1)
- Quercus pubescens (1)
- RCP scenarios (1)
- Raman spectroscopy (1)
- Reflectivity (1)
- Reservoir Networks (1)
- River (1)
- River Incision Model (1)
- Russia (1)
- SIO₂ (1)
- SRTM (1)
- SWAT (1)
- SWIM (1)
- Sahel (1)
- Salt Range (1)
- Samara (1)
- Santa-Barbara system (1)
- Sediment Flux (1)
- Sediment Transport (1)
- Sentinel 2 (1)
- Sentinel-2 MSI (1)
- Shear Zone (1)
- Siberia (1)
- Siberian larch (1)
- Slip Distribution (1)
- Slope (1)
- Solid Earth sciences (1)
- Southern Apennines (1)
- Southern Kyrgyzstan (1)
- Soziale Medien (1)
- Structure-from-motion (1)
- Subduction (1)
- Swarm constellation (1)
- TIME-GCM (1)
- TRMM (1)
- TanDEM-X (1)
- Tarim Basin (1)
- Taylor’s law (1)
- Tectonic Evolution (1)
- Teteriv (1)
- Tianshan orogenic belt (1)
- Tonian–Cryogenian (1)
- Tropospheric nitrogen-dioxide (1)
- Turbidity retrieval (1)
- Twitter (1)
- Ukraine (1)
- Uk’37 (1)
- Uncertainty Processor (1)
- Unterrichtsmethoden (1)
- VI (1)
- VIIRS DNB (1)
- Validation (1)
- Variscan orogeny (1)
- WGHM (1)
- Water Availability (1)
- Weather (1)
- Western Bug (1)
- Western Interior Basin (1)
- Western Qaidam Basin (1)
- WorldDEM (1)
- WorldView-2 (1)
- Yamal (1)
- Younger Dryas (1)
- Zipf’s law (1)
- accumulation (1)
- acidophilic microorganisms (1)
- adaption measure (1)
- agriculture (1)
- airglow (1)
- alaskan tundra (1)
- alpine (1)
- alpine catchments (1)
- alps (1)
- annually laminated sediments (1)
- architecture (1)
- area (1)
- artificial light (1)
- assimilation (1)
- astronomical calibration (1)
- automatic processing (1)
- barrow (1)
- bat fatalities (1)
- biodiversity decline (1)
- bioenergy (1)
- biofuel (1)
- biomaterials (1)
- blueSeis-3A sensors (1)
- boreal forests (1)
- borehole leakage (1)
- bottom waters (1)
- calc-alkaline (1)
- calcium phosphate (1)
- calcium phosphate hybrid material (1)
- calibration (1)
- canopy height (1)
- carbon dioxide exchange (1)
- carbon export (1)
- carbon-cycle (1)
- catchments (1)
- cell counts (1)
- cellulose (1)
- cement (1)
- central Mediterranean (1)
- central Poland (1)
- central south island (1)
- change detection (1)
- chemical-mechanical interaction (1)
- chitosan (1)
- chronology (AICC2012) (1)
- circulation (1)
- cirrus detection (1)
- city clusters (1)
- classification (1)
- climate change adaptation (1)
- climate change impact (1)
- climate networks (1)
- climate variability (1)
- clinopyroxenite (1)
- cloud detection (1)
- coast (1)
- coherence (1)
- coherency (1)
- coherent noise (1)
- collapse (1)
- complexes (1)
- composite properties (1)
- conflict (1)
- continental lithosphere (1)
- convective available potential energy (1)
- convective available potential energy (CAPE) (1)
- cosmogenic AL-26 (1)
- critical meteorological condition (1)
- crust (1)
- crystalline (1)
- data fusion (1)
- data science (1)
- decision tree (1)
- decision trees (1)
- deep convection (1)
- deep crust (1)
- deep permafrost carbon (1)
- deep seated landslide (1)
- deformation monitoring (1)
- degradation (1)
- delta-c-13 (1)
- dendrometer measurements (1)
- detrital carbonate (1)
- dew-point temperature (1)
- digital rock physics (1)
- dip-coating (1)
- dissolution (1)
- ditch irrigation (1)
- diversity (1)
- drained lake basins (1)
- drought (1)
- dust sources (1)
- early eocene (1)
- early warning (1)
- earth system model (1)
- earth system models (1)
- eastern south–central Andes (1)
- eclogite (1)
- eco-hydrology (1)
- ecological genetics (1)
- ecological modelling (1)
- elastic properties (1)
- electrons (1)
- elevated CO2 concentration (1)
- emergency response (1)
- energy efficiency (1)
- environmental seismology (1)
- eocene thermal maximum (1)
- equifinality (1)
- erosion (1)
- erosion rates (1)
- evaporite minerals (1)
- evaporites (1)
- evapotranspiration (1)
- event coincidence analysis (1)
- evolution (1)
- exhumation (1)
- extension (1)
- extinction event (1)
- extracellular DNA (1)
- fatigue (1)
- fault (1)
- fault bend fold (1)
- fire emissions (1)
- fire history (1)
- flat-slab subduction (1)
- flood adaptation (1)
- flood mitigation (1)
- floodplain (1)
- floodplains (1)
- fluctuations (1)
- flux tracking (1)
- food web (1)
- forecasting (1)
- foreland basin (1)
- forest (1)
- forest structure (1)
- frequency analysis (1)
- frictional properties (1)
- fundamental frequency (1)
- gabbros (1)
- gas production (1)
- gelatin (1)
- geo-bio interaction (1)
- geochronology (1)
- geomagnetic observatory data (1)
- geomagnetic storm drivers (1)
- glacial hazards (1)
- glacial lake outburst floods (1)
- glacier (1)
- glass (1)
- global (1)
- global monsoon (1)
- global vegetation model (1)
- grain-size distribution (1)
- granulite (1)
- graphite (1)
- green-green dilemma (1)
- gross primary production (1)
- ground-motion modelling (1)
- groundwater level (1)
- groundwater surface water interaction (1)
- groundwater vulnerability (1)
- growing season (1)
- growth response (1)
- gulf coastal plain (1)
- half-life (1)
- hanging wall (1)
- hazard (1)
- heat-stress (1)
- heuristic analysis (1)
- high latitudes (1)
- high-P metamorphism (1)
- high-intensity industrialized agricultural production systems (1)
- high-resolution record (1)
- historical geomagnetic storms (1)
- huddle test (1)
- hybrid modelling (1)
- hydro-meterological hazards (1)
- hydroclimatology (1)
- hydrogel (1)
- hydrogenase (1)
- hydrogeology (1)
- hydrological modeling (1)
- hydrological modelling (1)
- hydrological signatures (1)
- hydropower (1)
- ice (1)
- ice sheet (1)
- identifying influential nodes (1)
- induced chlorophyll Fluorescence (1)
- innovation (1)
- insights (1)
- insurance (1)
- interferometric SAR (InSAR) (1)
- internal structure (1)
- introgression (1)
- inundation (1)
- ionic conductivity (1)
- ionospheric scale lengths (1)
- iron reduction (1)
- iron-rich sediment (1)
- isotopic composition (1)
- isturbance tracking (1)
- lacustrine sediment (1)
- lake (1)
- lake Donggi Cona (1)
- lake Meerfelder Maar (1)
- lake basins (1)
- lake catchments (1)
- lake drainage (1)
- lakes (1)
- land cover change (1)
- landscape transience (1)
- larch species (1)
- last glacial period (1)
- latitude ionosphere (1)
- length measurements (1)
- lidar (1)
- localized flooding (1)
- logistic regression (1)
- low-intensity eco-agricultural production systems (1)
- lower crustal bodies (1)
- lower thermosphere (1)
- machine learning (1)
- magma-poor (1)
- magmatic arc (1)
- magnetic field (1)
- magnetosphere (1)
- margins (1)
- mass (1)
- melt inclusions (1)
- metagenomic analysis (1)
- metal recycling plants (1)
- metaophiolite (1)
- metasomatism (1)
- meteoric 10Be (1)
- micro-CT (1)
- microbial activity (1)
- microtopography (1)
- mid-latitude Westerlies (1)
- middle atmosphere (1)
- mineral mapping (1)
- mitigation (1)
- mitochondrial haplotypes (1)
- mountain geomorphology (1)
- multi-slab (1)
- multispectral (1)
- multitemporal (1)
- nZVI (1)
- nano zero-valent iron (1)
- natural hazard (1)
- natural hazard management (1)
- natural hazards (1)
- nearshore zone (1)
- neodymium (1)
- net ecosystem exchange (1)
- network (1)
- nightlights (1)
- nitrogen (1)
- noise cross (1)
- nonstationarity (1)
- north-atlantic (1)
- northern high latitudes (1)
- northern seward peninsula (1)
- numerical model (1)
- numerical simulation (1)
- object-based damage modeling (1)
- object-based image analysis (1)
- ocean acidification (1)
- ocean color remote sensing (1)
- oligocene climate (1)
- open (1)
- operational service (1)
- ophiolite (1)
- opinion (1)
- optical data (1)
- organic carbon (1)
- origination (1)
- orogenic peridotite (1)
- orogeny (1)
- oxygen (1)
- oxygen isotope (1)
- palaeoecology (1)
- paleo-sediment (1)
- paleosol sequence (1)
- part 1 (1)
- passive microwave (1)
- patterns (1)
- percolation (1)
- permafrost disturbances (1)
- permafrost ecosystems (1)
- permafrost regions (1)
- permafrost thaw (1)
- permeability test (1)
- photogrammetry (1)
- plant ecology (1)
- plateau (1)
- plume (1)
- pluvial floods (1)
- point cloud (1)
- pollution (1)
- pollution indices (1)
- poly(lactic acid) (1)
- potentially affected areas (PAA) (1)
- potentially toxic elements (1)
- precipitation (1)
- pressures (1)
- principal components (1)
- prior derivation (1)
- probabilistic approaches (1)
- production rates (1)
- productivity (1)
- radiation-damage (1)
- radiocarbon (1)
- radiogenic heat-production (1)
- railway transportation (1)
- rainfall (1)
- ramp (1)
- random vibration (1)
- range shifts (1)
- rapid damage mapping (1)
- rare earth elements (1)
- reactivation (1)
- recovery (1)
- recycling (1)
- reduction (1)
- reflection (1)
- regimes (1)
- release (1)
- remediation (1)
- remote sensing data (1)
- renewable energy (1)
- required minimum runoff (1)
- reservoir (1)
- response time (1)
- risk assessment (1)
- risk management cycle (1)
- river discharge (1)
- river sediment (1)
- rock mass (1)
- rotation (1)
- rotational seismology (1)
- rutile-type (1)
- saline pan (1)
- saline pan cycle (1)
- salt dissolution (1)
- san andreas fault (1)
- savanna (1)
- scale (1)
- scenario (1)
- sea (1)
- sea-level projections (1)
- seasonality (1)
- sediment (1)
- sediment routing (1)
- sediment transit time (1)
- sedimentary ancient DNA (1)
- sedimentary basins (1)
- sedimentary microbes (1)
- sedimentation pattern (1)
- seismic monitoring (1)
- shadow detection (1)
- shoreline (1)
- shotgun sequencing (1)
- signal propagation (1)
- silicon (1)
- simulation (1)
- simulation-based attribution (1)
- singleslab (1)
- site effects (1)
- skill (1)
- slags (1)
- slope failure (1)
- small baseline subset (SBAS) (1)
- snow detection (1)
- social (1)
- soil landscape (1)
- soils (1)
- source direction (1)
- source-to-sink (1)
- spatial correlation (1)
- spatial scales (1)
- spatial-distribution (1)
- special sensor microwave imager (1)
- special sensor microwave imager/sounder (1)
- spectral adjustment (1)
- spectral analysis (1)
- strain localization (1)
- stratigraphy (1)
- streamflow variability (1)
- structure from motion (1)
- subduction zone (1)
- subrosion (1)
- subsidence analysis (1)
- subsurface biosphere (1)
- sulfate reduction (1)
- summer monsoon (1)
- supervised machine learning (1)
- surface reflectance (1)
- surface water flooding (1)
- suspended sediment (1)
- sustained casing pressure (1)
- switzerland (1)
- synchronization (1)
- systems (1)
- tectonic evolution (1)
- tectonic inversion (1)
- tectonics (1)
- temperature changes (1)
- theory (1)
- thermal field (1)
- thermal properties (1)
- thermokarst lakes (1)
- thickness-displacement relationships (1)
- thrust belts (1)
- tibetan plateau (1)
- tide gauge observations (1)
- time series (1)
- time-scale (1)
- time-series mapping (1)
- transferability (1)
- transit-time (1)
- transport processes (1)
- travel time distributions (1)
- treeline (1)
- trend attribution (1)
- trend detection (1)
- trend drivers (1)
- tritium assay (1)
- tropical biodiversity (1)
- tropical swamp (1)
- tundra (1)
- tundra–taiga ecotone (1)
- uncertainty analysis (1)
- united states (1)
- urban infill development (1)
- variability (1)
- varved sediments (1)
- vegetation (1)
- vegetation history (1)
- vegetation states (1)
- vegetation trajectories (1)
- vegetation-climate feedbacks (1)
- vegetation‐climate‐fire‐soil feedbacks (1)
- very slow moving landslide (1)
- voring basin (1)
- wasting (1)
- water age (1)
- water balance (1)
- water balance model (1)
- water budget (1)
- water fluxes (1)
- water resources (1)
- weakening mechanisms (1)
- weather (1)
- wetland (1)
- wetlands (1)
- wind energy production (1)
- wind energy-biodiversity (1)
- winds (1)
- winter (1)
Institute
- Institut für Geowissenschaften (139) (remove)
Strong hydroclimatic controls on vulnerability to subsurface nitrate contamination across Europe
(2020)
Subsurface contamination due to excessive nutrient surpluses is a persistent and widespread problem in agricultural areas across Europe. The vulnerability of a particular location to pollution from reactive solutes, such as nitrate, is determined by the interplay between hydrologic transport and biogeochemical transformations. Current studies on the controls of subsurface vulnerability do not consider the transient behaviour of transport dynamics in the root zone. Here, using state-of-the-art hydrologic simulations driven by observed hydroclimatic forcing, we demonstrate the strong spatiotemporal heterogeneity of hydrologic transport dynamics and reveal that these dynamics are primarily controlled by the hydroclimatic gradient of the aridity index across Europe. Contrasting the space-time dynamics of transport times with reactive timescales of denitrification in soil indicate that similar to 75% of the cultivated areas across Europe are potentially vulnerable to nitrate leaching for at least onethird of the year. We find that neglecting the transient nature of transport and reaction timescale results in a great underestimation of the extent of vulnerable regions by almost 50%. Therefore, future vulnerability and risk assessment studies must account for the transient behaviour of transport and biogeochemical transformation processes.
We present a new set of global and local sea‐level projections at example tide gauge locations under the RCP2.6, RCP4.5, and RCP8.5 emissions scenarios. Compared to the CMIP5‐based sea‐level projections presented in IPCC AR5, we introduce a number of methodological innovations, including (i) more comprehensive treatment of uncertainties, (ii) direct traceability between global and local projections, and (iii) exploratory extended projections to 2300 based on emulation of individual CMIP5 models. Combining the projections with observed tide gauge records, we explore the contribution to total variance that arises from sea‐level variability, different emissions scenarios, and model uncertainty. For the period out to 2300 we further breakdown the model uncertainty by sea‐level component and consider the dependence on geographic location, time horizon, and emissions scenario. Our analysis highlights the importance of local variability for sea‐level change in the coming decades and the potential value of annual‐to‐decadal predictions of local sea‐level change. Projections to 2300 show a substantial degree of committed sea‐level rise under all emissions scenarios considered and highlight the reduced future risk associated with RCP2.6 and RCP4.5 compared to RCP8.5. Tide gauge locations can show large ( > 50%) departures from the global average, in some cases even reversing the sign of the change. While uncertainty in projections of the future Antarctic ice dynamic response tends to dominate post‐2100, we see substantial differences in the breakdown of model variance as a function of location, time scale, and emissions scenario.
Lake sediments are increasingly explored as reliable paleoflood archives. In addition to established flood proxies including detrital layer thickness, chemical composition, and grain size, we explore stable oxygen and carbon isotope data as paleoflood proxies for lakes in catchments with carbonate bedrock geology. In a case study from Lake Mondsee (Austria), we integrate high-resolution sediment trapping at a proximal and a distal location and stable isotope analyses of varved lake sediments to investigate flood-triggered detrital sediment flux. First, we demonstrate a relation between runoff, detrital sediment flux, and isotope values in the sediment trap record covering the period 2011-2013 CE including 22 events with daily (hourly) peak runoff ranging from 10 (24) m(3) s(-1) to 79 (110) m(3) s(-1). The three- to ten-fold lower flood-triggered detrital sediment deposition in the distal trap is well reflected by attenuated peaks in the stable isotope values of trapped sediments. Next, we show that all nine flood-triggered detrital layers deposited in a sediment record from 1988 to 2013 have elevated isotope values compared with endogenic calcite. In addition, even two runoff events that did not cause the deposition of visible detrital layers are distinguished by higher isotope values. Empirical thresholds in the isotope data allow estimation of magnitudes of the majority of floods, although in some cases flood magnitudes are overestimated because local effects can result in too-high isotope values. Hence we present a proof of concept for stable isotopes as reliable tool for reconstructing flood frequency and, although with some limitations, even for flood magnitudes.
Diet analysis of bats killed at wind turbines suggests large-scale losses of trophic interactions
(2022)
Agricultural practice has led to landscape simplification and biodiversity decline, yet recently, energy-producing infrastructures, such as wind turbines, have been added to these simplified agroecosystems, turning them into multi-functional energy-agroecosystems. Here, we studied the trophic interactions of bats killed at wind turbines using a DNA metabarcoding approach to shed light on how turbine-related bat fatalities may possibly affect local habitats. Specifically, we identified insect DNA in the stomachs of common noctule bats (Nyctalus noctula) killed by wind turbines in Germany to infer in which habitats these bats hunted. Common noctule bats consumed a wide variety of insects from different habitats, ranging from aquatic to terrestrial ecosystems (e.g., wetlands, farmland, forests, and grasslands). Agricultural and silvicultural pest insects made up about 20% of insect species consumed by the studied bats. Our study suggests that the potential damage of wind energy production goes beyond the loss of bats and the decline of bat populations. Bat fatalities at wind turbines may lead to the loss of trophic interactions and ecosystem services provided by bats, which may add to the functional simplification and impaired crop production, respectively, in multi-functional ecosystems.
The sediment profile from Lake Goscia(z) over dot in central Poland comprises a continuous, seasonally resolved and exceptionally well-preserved archive of the Younger Dryas (YD) climate variation. This provides a unique opportunity for detailed investigation of lake system responses during periods of rapid climate cooling (YD onset) and warming (YD termination). The new varve record of Lake Goscia(z) over dot presented here spans 1662 years from the late Allerod (AL) to the early Preboreal (PB). Microscopic varve counting provides an independent chronology with a YD duration of 1149+14/-22 years, which confirms previous results of 1140 +/- 40 years. We link stable oxygen isotopes and chironomid-based air temperature reconstructions with the response of various geochemical and varve microfacies proxies especially focusing on the onset and termination of the YD. Cooling at the YD onset lasted similar to 180 years, which is about a century longer than the terminal warming that was completed in similar to 70 years. During the AL/YD transition, environmental proxy data lagged the onset of cooling by similar to 90 years and revealed an increase of lake productivity and internal lake re-suspension as well as slightly higher detrital sediment input. In contrast, rapid warming and environmental changes during the YD/PB transition occurred simultaneously. However, initial changes such as declining diatom deposition and detrital input occurred already a few centuries before the rapid warming at the YD/PB transition. These environmental changes likely reflect a gradual increase in summer air temperatures already during the YD. Our data indicate complex and differing environmental responses to the major climate changes related to the YD, which involve different proxy sensitivities and threshold processes.
Large rock slope failures play a pivotal role in long-term landscape evolution and are a major concern in land use planning and hazard aspects. While the failure phase and the time immediately prior to failure are increasingly well studied, the nature of the preparation phase remains enigmatic. This knowledge gap is due, to a large degree, to difficulties associated with instrumenting high mountain terrain and the local nature of classic monitoring methods, which does not allow integral observation of large rock volumes. Here, we analyse data from a small network of up to seven seismic sensors installed during July-October 2018 (with 43 days of data loss) at the summit of the Hochvogel, a 2592 m high Alpine peak. We develop proxy time series indicative of cyclic and progressive changes of the summit. Modal analysis, horizontal-to-vertical spectral ratio data and end-member modelling analysis reveal diurnal cycles of increasing and decreasing coupling stiffness of a 260,000 m(3) large, instable rock volume, due to thermal forcing. Relative seismic wave velocity changes also indicate diurnal accumulation and release of stress within the rock mass. At longer time scales, there is a systematic superimposed pattern of stress increased over multiple days and episodic stress release within a few days, expressed in an increased emission of short seismic pulses indicative of rock cracking. Our data provide essential first order information on the development of large-scale slope instabilities towards catastrophic failure. (c) 2020 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd
The Big Naryn Complex (BNC) in the East Djetim-Too Range of the Kyrgyz Middle Tianshan block is a tectonized, at least 2 km thick sequence of predominantly felsic to intermediate volcanic rocks intruded by porphyric rhyolite sills. It overlies a basement of metamorphic rocks and is overlain by late Neoproterozoic Djetim-Too Formation sediments; these also occur as tectonic intercalations in the BNC. The up to ca. 1100 m thick Lower Member is composed of predominantly rhyolites-to-dacites and minor basalts, while the at least 900 m thick pyroclastic Upper Member is dominated by rhyolitic-to-dacitic ignimbrites. Porphyric rhyolite sills are concentrated at the top of the Lower Member. A Lower Member rhyolite and a sill sample have LA-ICP-MS U-Pb zircon crystallization ages of 726.1 +/- 2.2 Ma and 720.3 +/- 6.5 Ma, respectively, showing that most of the magmatism occurred within a short time span in the late Tonian-early Cryogenian. Inherited zircons in the sill sample have Neoarchean (2.63, 2.64 Ga), Paleo- (2.33-1.81 Ga), Meso- (1.55 Ga), and Neoproterozoic (ca. 815 Ma) ages, and were derived from a heterogeneous Kuilyu Complex basement. A 1751 +/- 7 Ma Ar-40/Ar-39 age for amphibole from metagabbro is the age of cooling subsequent to Paleoproterozoic metamorphism of the Kuilyu Complex. The large amount of pyroclastic rocks, and their major and trace element compositions, the presence of Neoarchean to Neoproterozoic inherited zircons and a depositional basement of metamorphic rocks point to formation of the BNC in a continental magmatic arc setting.
The Salt Range in Pakistan exposes Precambrian to Pleistocene strata outcropping along the Salt Range Thrust (SRT). To better understand the in-situ Cambrian and Pliocene tectonic evolution of the Pakistan Subhimalaya, we have conducted low-temperature thermochronological analysis using apatite (U-Th-Sm)/He and fission track dating. We combine cooling ages from different samples located along the thrust front of the SRT into a thermal model that shows two major cooling events associated with rifting and regional erosion in the Late Palaeozoic and SRT activity since the Pliocene. Our results suggest that the SRT maintained a long-term average shortening rate of similar to 5-6 mm/yr and a high exhumation rate above the SRT ramp since similar to 4 Ma.
The simulation of broad-band (0.1 to 10 + Hz) ground-shaking over deep and spatially extended sedimentary basins at regional scales is challenging. We evaluate the ground-shaking of a potential M 6.5 earthquake in the southern Lower Rhine Embayment, one of the most important areas of earthquake recurrence north of the Alps, close to the city of Cologne in Germany. In a first step, information from geological investigations, seismic experiments and boreholes is combined for deriving a harmonized 3D velocity and attenuation model of the sedimentary layers. Three alternative approaches are then applied and compared to evaluate the impact of the sedimentary cover on ground-motion amplification. The first approach builds on existing response spectra ground-motion models whose amplification factors empirically take into account the influence of the sedimentary layers through a standard parameterization. In the second approach, site-specific 1D amplification functions are computed from the 3D basin model. Using a random vibration theory approach, we adjust the empirical response spectra predicted for soft rock conditions by local site amplification factors: amplifications and associated ground-motions are predicted both in the Fourier and in the response spectra domain. In the third approach, hybrid physics-based ground-motion simulations are used to predict time histories for soft rock conditions which are subsequently modified using the 1D site-specific amplification functions computed in method 2. For large distances and at short periods, the differences between the three approaches become less notable due to the significant attenuation of the sedimentary layers. At intermediate and long periods, generic empirical ground-motion models provide lower levels of amplification from sedimentary soils compared to methods taking into account site-specific 1D amplification functions. In the near-source region, hybrid physics-based ground-motions models illustrate the potentially large variability of ground-motion due to finite source effects.
Borehole leakage is a common and complex issue. Understanding the fluid flow characteristics of a cemented area inside a borehole is crucial to monitor and quantify the wellbore integrity as well as to find solutions to minimise existing leakages. In order to improve our understanding of the flow behaviour of cemented boreholes, we investigated experimental data of a large-scale borehole leakage tests by means of numerical modelling using three different conceptual models. The experiment was performed with an autoclave system consisting of two vessels bridged by a cement-filled casing. After a partial bleed-off at the well-head, a sustained casing pressure was observed due to fluid flow through the cementsteel composite. The aim of our simulations is to investigate and quantify the permeability of the cement-steel composite. From our model results, we conclude that the flow occurred along a preferential flow path at the cement-steel interface. Thus, the inner part of the cement core was impermeable during the duration of the experiment. The preferential flow path can be described as a highly permeable and highly porous area with an aperture of about 5 mu m and a permeability of 3 . 10(-12) m(2) (3 Darcy). It follows that the fluid flow characteristics of a cemented area inside a borehole cannot be described using one permeability value for the entire cement-steel composite. Furthermore, it can be concluded that the quality of the cement and the filling process regarding the cement-steel interface is crucial to minimize possible well leakages.
The southern Central Andes (SCA) (between 27 degrees S and 40 degrees S) is bordered to the west by the convergent margin between the continental South American Plate and the oceanic Nazca Plate. The subduction angle along this margin is variable, as is the deformation of the upper plate. Between 33 degrees S and 35 degrees S, the subduction angle of the Nazca plate increases from sub-horizontal (< 5 degrees) in the north to relatively steep (similar to 30 degrees) in the south. The SCA contain inherited lithological and structural heterogeneities within the crust that have been reactivated and overprinted since the onset of subduction and associated Cenozoic deformation within the Andean orogen. The distribution of the deformation within the SCA has often been attributed to the variations in the subduction angle and the reactivation of these inherited heterogeneities. However, the possible influence that the thickness and composition of the continental crust have had on both short-term and long-term deformation of the SCA is yet to be thoroughly investigated. For our investigations, we have derived density distributions and thicknesses for various layers that make up the lithosphere and evaluated their relationships with tectonic events that occurred over the history of the Andean orogeny and, in particular, investigated the short- and long-term nature of the present-day deformation processes. We established a 3D model of lithosphere beneath the orogen and its foreland (29 degrees S-39 degrees S) that is consistent with currently available geological and geophysical data, including the gravity data. The modelled crustal configuration and density distribution reveal spatial relationships with different tectonic domains: the crystalline crust in the orogen (the magmatic arc and the main orogenic wedge) is thicker (similar to 55 km) and less dense (similar to 2900 kg/m(3)) than in the forearc (similar to 35 km, similar to 2975 kg/m(3)) and foreland (similar to 30 km, similar to 3000 kg/m(3)). Crustal thickening in the orogen probably occurred as a result of stacking of low-density domains, while density and thickness variations beneath the forearc and foreland most likely reflect differences in the tectonic evolution of each area following crustal accretion. No clear spatial relationship exists between the density distribution within the lithosphere and previously proposed boundaries of crustal terranes accreted during the early Paleozoic. Areas with ongoing deformation show a spatial correlation with those areas that have the highest topographic gradients and where there are abrupt changes in the average crustal-density contrast. This suggests that the short-term deformation within the interior of the Andean orogen and its foreland is fundamentally influenced by the crustal composition and the relative thickness of different crustal layers. A thicker, denser, and potentially stronger lithosphere beneath the northern part of the SCA foreland is interpreted to have favoured a strong coupling between the Nazca and South American plates, facilitating the development of a sub-horizontal slab.
Arctic lakes located in permafrost regions are susceptible to catastrophic drainage. In this study, we reconstructed historical lake drainage events on the western Arctic Coastal Plain of Alaska between 1955 and 2017 using USGS topographic maps, historical aerial photography (1955), and Landsat Imagery (ca. 1975, ca. 2000, and annually since 2000). We identified 98 lakes larger than 10 ha that partially (>25% of area) or completely drained during the 62-year period. Decadal-scale lake drainage rates progressively declined from 2.0 lakes/yr (1955-1975), to 1.6 lakes/yr (1975-2000), and to 1.2 lakes/yr (2000-2017) in the ~30,000-km(2) study area. Detailed Landsat trend analysis between 2000 and 2017 identified two years, 2004 and 2006, with a cluster (five or more) of lake drainages probably associated with bank overtopping or headward erosion. To identify future potential lake drainages, we combined the historical lake drainage observations with a geospatial dataset describing lake elevation, hydrologic connectivity, and adjacent lake margin topographic gradients developed with a 5-m-resolution digital surface model. We identified ~1900 lakes likely to be prone to drainage in the future. Of the 20 lakes that drained in the most recent study period, 85% were identified in this future lake drainage potential dataset. Our assessment of historical lake drainage magnitude, mechanisms and pathways, and identification of potential future lake drainages provides insights into how arctic lowland landscapes may change and evolve in the coming decades to centuries.
Most hydrological studies rely on a model calibrated using discharge alone. However, judging the model reliability based on such calibration is problematic, as it does not guarantee the correct representation of internal hydrological processes. This study aims (a) to develop a comprehensive multi-objective calibration framework using remote sensing vegetation data and hydrological signatures (flow duration curve - FDC, and baseflow index) in addition to discharge, and (b) to apply this framework for calibration of the Soil and Water Assessment Tool (SWAT) in a typical Andean catchment. Overall, our calibration approach outperformed traditional discharge-based and FDC signature-based calibration strategies in terms of vegetation, streamflow, and flow partitioning simulation. New hydrological insights for the region are the following: baseflow is the main component of the streamflow sustaining the long dry-season flow, and pasture areas offer higher water yield and baseflow than other land-cover types. The proposed approach could be used in other data-scarce regions with complex topography.
Due to the high concentration of people and infrastructures in European cities, the possible impacts of climate change are particularly high (cities' social, economic and technical vulnerabilities). Adaptation measures to reduce the sensitivity of a city to climate risks are therefore of particular importance. Nevertheless, it is also common to develop compact and dense urban areas to reduce urban sprawl. Urban infill development and sustainable spatial climate policies are thus in apparent conflict with each other. This article examines how German cities deal with the tensions between these two policy fields. Using six case studies, a new heuristic analysis method is applied. This study identifies three key governance aspects that are essential for promoting the joint implementation: instruments, organisation and interaction. Based on our case studies, we conclude that successful implementation can only be achieved through integrative governance including all three domains.
The hyperthermal events of the Cenozoic, including the Paleocene-Eocene Thermal Maximum, provide an opportunity to investigate the potential effects of climate warming on marine ecosystems. Here, we examine the shallow benthic marine communities preserved in the late Cretaceous to Eocene strata on the Gulf Coastal Plain (United States). In stark contrast to the ecological shifts following the end-Cretaceous mass extinction, our data show that the early Cenozoic hyperthermals did not have a long-term impact on the generic diversity nor composition of the Gulf Coastal Plain molluscan communities. We propose that these communities were resilient to climate change because molluscs are better adapted to high temperatures than other taxa, as demonstrated by their physiology and evolutionary history. In terms of resilience, these communities differ from other shallow-water carbonate ecosystems, such as reef communities, which record significant changes during the early Cenozoic hyperthermals. These data highlight the strikingly different responses of community types, i.e., the almost imperceptible response of molluscs versus the marked turnover of foraminifera and reef faunas. The impact on molluscan communities may have been low because detrimental conditions did not devastate the entire Gulf Coastal Plain, allowing molluscs to rapidly recolonise vacated areas once harsh environmental conditions ameliorated.
The origin of Asian monsoons
(2020)
The Cenozoic inception and development of the Asian monsoons remain unclear and have generated much debate, as several hypotheses regarding circulation patterns at work in Asia during the Eocene have been proposed in the few last decades. These include (a) the existence of modern-like monsoons since the early Eocene; (b) that of a weak South Asian monsoon (SAM) and little to no East Asian monsoon (EAM); or (c) a prevalence of the Intertropical Convergence Zone (ITCZ) migrations, also referred to as Indonesian-Australian monsoon (I-AM). As SAM and EAM are supposed to have been triggered or enhanced primarily by Asian palaeogeographic changes, their possible inception in the very dynamic Eocene palaeogeographic context remains an open question, both in the modelling and field-based communities. We investigate here Eocene Asian climate conditions using the IPSL-CM5A2 (Sepulchre et al., 2019) earth system model and revised palaeogeographies. Our Eocene climate simulation yields atmospheric circulation patterns in Asia substantially different from modern conditions. A large high-pressure area is simulated over the Tethys ocean, which generates intense low tropospheric winds blowing southward along the western flank of the proto-Himalayan-Tibetan plateau (HTP) system. This low-level wind system blocks, to latitudes lower than 10 degrees N, the migration of humid and warm air masses coming from the Indian Ocean. This strongly contrasts with the modern SAM, during which equatorial air masses reach a latitude of 20-25 degrees N over India and southeastern China. Another specific feature of our Eocene simulation is the widespread subsidence taking place over northern India in the midtroposphere (around 5000 m), preventing deep convective updraught that would transport water vapour up to the condensation level. Both processes lead to the onset of a broad arid region located over northern India and over the HTP. More humid regions of high seasonality in precipitation encircle this arid area, due to the prevalence of the Intertropical Convergence Zone (ITCZ) migrations (or Indonesian-Australian monsoon, I-AM) rather than monsoons. Although the existence of this central arid region may partly result from the specifics of our simulation (model dependence and palaeogeographic uncertainties) and has yet to be confirmed by proxy records, most of the observational evidence for Eocene monsoons are located in the highly seasonal transition zone between the arid area and the more humid surroundings. We thus suggest that a zonal arid climate prevailed over Asia before the initiation of monsoons that most likely occurred following Eocene palaeogeographic changes. Our results also show that precipitation seasonality should be used with caution to infer the presence of a monsoonal circulation and that the collection of new data in this arid area is of paramount importance to allow the debate to move forward.
Groundwater levels are monitored by environmental agencies to support the sustainable use of groundwater resources. For this purpose continuous and spatially comprehensive monitoring in high spatial and temporal resolution is desired. This leads to large datasets that have to be checked for quality and analysed to distinguish local anthropogenic influences from natural variability of the groundwater level dynamics at each well. Both technical problems with the measurements as well as local anthropogenic influences can lead to local anomalies in the hydrographs. We suggest a fast and efficient screening method for the identification of well-specific peculiarities in hydrographs of groundwater head monitoring networks. The only information required is a set of time series of groundwater heads all measured at the same instants of time. For each well of the monitoring network a reference hydrograph is calculated, describing expected “normal” behaviour at the respective well as is typical for the monitored region. The reference hydrograph is calculated by multiple linear regression of the observed hydrograph with the “stable” principal components (PCs) of a principal component analysis of all groundwater head series of the network as predictor variables. The stable PCs are those PCs which were found in a random subsampling procedure to be rather insensitive to the specific selection of the analysed observation wells, i.e. complete series, and to the specific selection of measurement dates. Hence they can be considered to be representative for the monitored region in the respective period. The residuals of the reference hydrograph describe local deviations from the normal behaviour. Peculiarities in the residuals allow the data to be checked for measurement errors and the wells with a possible anthropogenic influence to be identified. The approach was tested with 141 groundwater head time series from the state authority groundwater monitoring network in northeastern Germany covering the period from 1993 to 2013 at an approximately weekly frequency of measurement.
New Zealand's Alpine Fault is a large, platebounding strike-slip fault, which ruptures in large (M-w > 8) earthquakes. We conducted field and laboratory analyses of fault rocks to assess its fault zone architecture. Results reveal that the Alpine Fault Zone has a complex geometry, comprising an anastomosing network of multiple slip planes that have accommodated different amounts of displacement. This contrasts with the previous perception of the Alpine Fault Zone, which assumes a single principal slip zone accommodated all displacement. This interpretation is supported by results of drilling projects and geophysical investigations. Furthermore, observations presented here show that the young, largely unconsolidated sediments that constitute the footwall at shallow depths have a significant influence on fault gouge rheological properties and structure.
Much of contemporary landslide research is concerned with predicting and mapping susceptibility to slope failure. Many studies rely on generalised linear models with environmental predictors that are trained with data collected from within and outside of the margins of mapped landslides. Whether and how the performance of these models depends on sample size, location, or time remains largely untested. We address this question by exploring the sensitivity of a multivariate logistic regression-one of the most widely used susceptibility models-to data sampled from different portions of landslides in two independent inventories (i.e. a historic and a multi-temporal) covering parts of the eastern rim of the Fergana Basin, Kyrgyzstan. We find that considering only areas on lower parts of landslides, and hence most likely their deposits, can improve the model performance by >10% over the reference case that uses the entire landslide areas, especially for landslides of intermediate size. Hence, using landslide toe areas may suffice for this particular model and come in useful where landslide scars are vague or hidden in this part of Central Asia. The model performance marginally varied after progressively updating and adding more landslides data through time. We conclude that landslide susceptibility estimates for the study area remain largely insensitive to changes in data over about a decade. Spatial or temporal stratified sampling contributes only minor variations to model performance. Our findings call for more extensive testing of the concept of dynamic susceptibility and its interpretation in data-driven models, especially within the broader framework of landslide risk assessment under environmental and land-use change.
A new solid-state material, N-butyl pyridinium diiodido argentate(I), is synthesized using a simple and effective one-pot approach. In the solid state, the compound exhibits 1D ([AgI2](-))(n) chains that are stabilized by the N-butyl pyridinium cation. The 1D structure is further manifested by the formation of long, needle-like crystals, as revealed from electron microscopy. As the general composition is derived from metal halide-based ionic liquids, the compound has a low melting point of 100-101 degrees C, as confirmed by differential scanning calorimetry. Most importantly, the compound has a conductivity of 10(-6) S cm(-1) at room temperature. At higher temperatures the conductivity increases and reaches to 10(-4 )S cm(-1) at 70 degrees C. In contrast to AgI, however, the current material has a highly anisotropic 1D arrangement of the ionic domains. This provides direct and tuneable access to fast and anisotropic ionic conduction. The material is thus a significant step forward beyond current ion conductors and a highly promising prototype for the rational design of highly conductive ionic solid-state conductors for battery or solar cell applications.
OpenForecast
(2019)
The development and deployment of new operational runoff forecasting systems are a strong focus of the scientific community due to the crucial importance of reliable and timely runoff predictions for early warnings of floods and flashfloods for local businesses and communities. OpenForecast, the first operational runoff forecasting system in Russia, open for public use, is presented in this study. We developed OpenForecast based only on open-source software and data-GR4J hydrological model, ERA-Interim meteorological reanalysis, and ICON deterministic short-range meteorological forecasts. Daily forecasts were generated for two basins in the European part of Russia. Simulation results showed a limited efficiency in reproducing the spring flood of 2019. Although the simulations managed to capture the timing of flood peaks, they failed in estimating flood volume. However, further implementation of the parsimonious data assimilation technique significantly alleviates simulation errors. The revealed limitations of the proposed operational runoff forecasting system provided a foundation to outline its further development and improvement.
Forest structure is a crucial component in the assessment of whether a forest is likely to act as a carbon sink under changing climate. Detailed 3D structural information about the tundra–taiga ecotone of Siberia is mostly missing and still underrepresented in current research due to the remoteness and restricted accessibility. Field based, high-resolution remote sensing can provide important knowledge for the understanding of vegetation properties and dynamics. In this study, we test the applicability of consumer-grade Unmanned Aerial Vehicles (UAVs) for rapid calculation of stand metrics in treeline forests. We reconstructed high-resolution photogrammetric point clouds and derived canopy height models for 10 study sites from NE Chukotka and SW Yakutia. Subsequently, we detected individual tree tops using a variable-window size local maximum filter and applied a marker-controlled watershed segmentation for the delineation of tree crowns. With this, we successfully detected 67.1% of the validation individuals. Simple linear regressions of observed and detected metrics show a better correlation (R2) and lower relative root mean square percentage error (RMSE%) for tree heights (mean R2 = 0.77, mean RMSE% = 18.46%) than for crown diameters (mean R2 = 0.46, mean RMSE% = 24.9%). The comparison between detected and observed tree height distributions revealed that our tree detection method was unable to representatively identify trees <2 m. Our results show that plot sizes for vegetation surveys in the tundra–taiga ecotone should be adapted to the forest structure and have a radius of >15–20 m to capture homogeneous and representative forest stands. Additionally, we identify sources of omission and commission errors and give recommendations for their mitigation. In summary, the efficiency of the used method depends on the complexity of the forest’s stand structure.
In this study, we analyze interactions in lake and lake catchment systems of a continuous permafrost area. We assessed colored dissolved organic matter (CDOM) absorption at 440 nm (a(440)(CDOM)) and absorption slope (S300-500) in lakes using field sampling and optical remote sensing data for an area of 350 km(2) in Central Yamal, Siberia. Applying a CDOM algorithm (ratio of green and red band reflectance) for two high spatial resolution multispectral GeoEye-1 and Worldview-2 satellite images, we were able to extrapolate the a()(CDOM) data from 18 lakes sampled in the field to 356 lakes in the study area (model R-2 = 0.79). Values of a(440)(CDOM) in 356 lakes varied from 0.48 to 8.35 m(-1) with a median of 1.43 m(-1). This a()(CDOM) dataset was used to relate lake CDOM to 17 lake and lake catchment parameters derived from optical and radar remote sensing data and from digital elevation model analysis in order to establish the parameters controlling CDOM in lakes on the Yamal Peninsula. Regression tree model and boosted regression tree analysis showed that the activity of cryogenic processes (thermocirques) in the lake shores and lake water level were the two most important controls, explaining 48.4% and 28.4% of lake CDOM, respectively (R-2 = 0.61). Activation of thermocirques led to a large input of terrestrial organic matter and sediments from catchments and thawed permafrost to lakes (n = 15, mean a(440)(CDOM) = 5.3 m(-1)). Large lakes on the floodplain with a connection to Mordy-Yakha River received more CDOM (n = 7, mean a(440)(CDOM) = 3.8 m(-1)) compared to lakes located on higher terraces.
Many widely used observational data sets are comprised of several overlapping instrument records. While data inter-calibration techniques often yield continuous and reliable data for trend analysis, less attention is generally paid to maintaining higher-order statistics such as variance and autocorrelation. A growing body of work uses these metrics to quantify the stability or resilience of a system under study and potentially to anticipate an approaching critical transition in the system. Exploring the degree to which changes in resilience indicators such as the variance or autocorrelation can be attributed to non-stationary characteristics of the measurement process – rather than actual changes in the dynamical properties of the system – is important in this context. In this work we use both synthetic and empirical data to explore how changes in the noise structure of a data set are propagated into the commonly used resilience metrics lag-one autocorrelation and variance. We focus on examples from remotely sensed vegetation indicators such as vegetation optical depth and the normalized difference vegetation index from different satellite sources. We find that time series resulting from mixing signals from sensors with varied uncertainties and covering overlapping time spans can lead to biases in inferred resilience changes. These biases are typically more pronounced when resilience metrics are aggregated (for example, by land-cover type or region), whereas estimates for individual time series remain reliable at reasonable sensor signal-to-noise ratios. Our work provides guidelines for the treatment and aggregation of multi-instrument data in studies of critical transitions and resilience.
Records from ocean bottom seismometers (OBSs) are highly contaminated by noise, which is much stronger
compared to data from most land stations, especially on the horizontal components. As a consequence, the high energy of the oceanic noise at frequencies below 1 Hz considerably complicates the analysis of the teleseismic earthquake signals recorded by OBSs.
Previous studies suggested different approaches to remove low-frequency noises from OBS recordings but mainly focused on the vertical component. The records of horizontal components, which are crucial for the application of many methods in passive seismological analysis of body and surface waves, could not be much improved in the teleseismic frequency band. Here we introduce a noise reduction method, which is derived from the harmonic–percussive separation algorithms used in Zali et al. (2021), in order to separate long-lasting narrowband signals from broadband transients in the OBS signal. This leads to significant noise reduction of OBS records on both the vertical and horizontal components and increases the earthquake signal-to-noise ratio (SNR) without distortion of the broadband earthquake waveforms. This is demonstrated through tests with synthetic data. Both SNR and cross-correlation coefficients showed significant improvements for different realistic noise realizations. The application of denoised signals in surface wave analysis and receiver functions is discussed through tests with synthetic and real data.
The Permo-Triassic period marks the time interval between Hercynian (Variscan) orogenic events in the Tien Shan and the North Pamir, and the Cimmerian accretion of the Gondwana-derived Central and South Pamir to the southern margin of the Paleo-Asian continent. A well-preserved Permo-Triassic volcano-sedimentary sequence from the Chinese North Pamir yields important information on the geodynamic evolution of Asia’s pre-Cimmerian southern margin. The oldest volcanic rocks from that section are dated to the late Guadalupian epoch by a rhyolite and a dacitic dike that gave zircon U-Pb ages of ~260 Ma. Permian volcanism was largely pyroclastic and mafic to intermediate. Upsection, a massive ignimbritic crystal tuff in the Chinese Qimgan valley was dated to 244.1 +/- 1.1 Ma, a similar unit in the nearby Gez valley to 245 +/- 11 Ma, and an associated rhyolite to 233.4 +/- 1.1 Ma. Deposition of the locally ~200 m thick crystal tuff unit follows an unconformity and marks the onset of intense, mainly mafic to intermediate, calc-alkaline magmatic activity. Triassic volcanic activity in the North Pamir was coeval with the major phase of Cimmerian intrusive activity in the Karakul-Mazar arc-accretionary complex to the south, caused by northward subduction of the Paleo-Tethys. It also coincided with the emplacement of basanitic and carbonatitic dikes and a thermal event in the South Tien Shan, to the north of our study area. Evidence for arc-related magmatic activity in a back-arc position provides strong arguments for back-arc extension or transtension and basin formation. This puts the Qimgan succession in line with a more than 1000 km long realm of extensional Triassic back-arc basins known from the North Pamir in the Kyrgyz Altyn Darya valley (Myntekin formation), the North Pamir of Tajikistan and Afghanistan, and the Afghan Hindukush (Doab formation) and further west from the Paropamisus and Kopet Dag (Aghdarband, NE Iran).
Glaciated high-alpine areas are fundamentally altered by climate change, with well-known implications for hydrology, e.g., due to glacier retreat, longer snow-free periods, and more frequent and intense summer rainstorms. While knowledge on how these hydrological changes will propagate to suspended sediment dynamics is still scarce, it is needed to inform mitigation and adaptation strategies. To understand the processes and source areas most relevant to sediment dynamics, we analyzed discharge and sediment dynamics in high temporal resolution as well as their patterns on several spatial scales, which to date few studies have done.
We used a nested catchment setup in the Upper Ötztal in Tyrol, Austria, where high-resolution (15 min) time series of discharge and suspended sediment concentrations are available for up to 15 years (2006–2020). The catchments of the gauges in Vent, Sölden and Tumpen range from 100 to almost 800 km2 with 10 % to 30 % glacier cover and span an elevation range of 930 to 3772 m a.s.l. We analyzed discharge and suspended sediment yields (SSY), their distribution in space, their seasonality and spatial differences therein, and the relative importance of short-term events. We complemented our analysis by linking the observations to satellite-based snow cover maps, glacier inventories, mass balances and precipitation data.
Our results indicate that the areas above 2500 m a.s.l., characterized by glacier tongues and the most recently deglaciated areas, are crucial for sediment generation in all sub-catchments. This notion is supported by the synchronous spring onset of sediment export at the three gauges, which coincides with snowmelt above 2500 m but lags behind spring discharge onsets. This points at a limitation of suspended sediment supply as long as the areas above 2500 m are snow-covered. The positive correlation of annual SSY with glacier cover (among catchments) and glacier mass balances (within a catchment) further supports the importance of the glacier-dominated areas. The analysis of short-term events showed that summer precipitation events were associated with peak sediment concentrations and yields but on average accounted for only 21 % of the annual SSY in the headwaters. These results indicate that under current conditions, thermally induced sediment export (through snow and glacier melt) is dominant in the study area.
Our results extend the scientific knowledge on current hydro-sedimentological conditions in glaciated high-alpine areas and provide a baseline for studies on projected future changes in hydro-sedimentological system dynamics.
The first step towards assessing hazards in seismically active regions involves mapping capable faults and estimating their recurrence times. While the mapping of active faults is commonly based on distinct geologic and geomorphic features evident at the surface, mapping blind seismogenic faults is complicated by the absence of on-fault diagnostic features. Here we investigated the Pichilemu Fault in coastal Chile, unknown until it generated a Mw 7.0 earthquake in 2010. The lack of evident surface faulting suggests activity along a partly-hidden blind fault. We used off-fault deformed marine terraces to estimate a fault-slip rate of 0.52 ± 0.04 m/ka, which, when integrated with satellite geodesy suggests a 2.12 ± 0.2 ka recurrence time for Mw~7.0 normal-faulting earthquakes. We propose that extension in the Pichilemu region is associated with stress changes during megathrust earthquakes and accommodated by sporadic slip during upper-plate earthquakes, which has implications for assessing the seismic potential of cryptic faults along convergent margins and elsewhere.
River-valley morphology preserves information on tectonic and climatic conditions that shape landscapes. Observations suggest that river discharge and valley-wall lithology are the main controls on valley width. Yet, current models based on these observations fail to explain the full range of cross-sectional valley shapes in nature, suggesting hitherto unquantified controls on valley width. In particular, current models cannot explain the existence of paired terrace sequences that form under cyclic climate forcing. Paired river terraces are staircases of abandoned floodplains on both valley sides, and hence preserve past valley widths. Their formation requires alternating phases of predominantly river incision and predominantly lateral planation, plus progressive valley narrowing. While cyclic Quaternary climate changes can explain shifts between incision and lateral erosion, the driving mechanism of valley narrowing is unknown. Here, we extract valley geometries from climatically formed, alluvial river-terrace sequences and show that across our dataset, the total cumulative terrace height (here: total valley height) explains 90%–99% of the variance in valley width at the terrace sites. This finding suggests that valley height, or a parameter that scales linearly with valley height, controls valley width in addition to river discharge and lithology. To explain this valley-width-height relationship, we reformulate existing valley-width models and suggest that, when adjusting to new boundary conditions, alluvial valleys evolve to a width at which sediment removal from valley walls matches lateral sediment supply from hillslope erosion. Such a hillslope-channel coupling is not captured in current valley-evolution models. Our model can explain the existence of paired terrace sequences under cyclic climate forcing and relates valley width to measurable field parameters. Therefore, it facilitates the reconstruction of past climatic and tectonic conditions from valley topography.
Sea level rise and coastal erosion have inundated large areas of Arctic permafrost. Submergence by warm and saline waters increases the rate of inundated permafrost thaw compared to sub-aerial thawing on land. Studying the contact between the unfrozen and frozen sediments below the seabed, also known as the ice-bearing permafrost table (IBPT), provides valuable information to understand the evolution of sub-aquatic permafrost, which is key to improving and understanding coastal erosion prediction models and potential greenhouse gas emissions. In this study, we use data from 2D electrical resistivity tomography (ERT) collected in the nearshore coastal zone of two Arctic regions that differ in their environmental conditions (e.g., seawater depth and resistivity) to image and study the subsea permafrost. The inversion of 2D ERT data sets is commonly performed using deterministic approaches that favor smoothed solutions, which are typically interpreted using a user-specified resistivity threshold to identify the IBPT position. In contrast, to target the IBPT position directly during inversion, we use a layer-based model parameterization and a global optimization approach to invert our ERT data. This approach results in ensembles of layered 2D model solutions, which we use to identify the IBPT and estimate the resistivity of the unfrozen and frozen sediments, including estimates of uncertainties. Additionally, we globally invert 1D synthetic resistivity data and perform sensitivity analyses to study, in a simpler way, the correlations and influences of our model parameters. The set of methods provided in this study may help to further exploit ERT data collected in such permafrost environments as well as for the design of future field experiments.
The Ice, Cloud, and Land Elevation Satellite-2 (ICESat-2) with its land and vegetation height data product (ATL08), and Global Ecosystem Dynamics Investigation (GEDI) with its terrain elevation and height metrics data product (GEDI Level 2A) missions have great potential to globally map ground and canopy heights. Canopy height is a key factor in estimating above-ground biomass and its seasonal changes; these satellite missions can also improve estimated above-ground carbon stocks. This study presents a novel Sparse Vegetation Detection Algorithm (SVDA) which uses ICESat-2 (ATL03, geolocated photons) data to map tree and vegetation heights in a sparsely vegetated savanna ecosystem. The SVDA consists of three main steps: First, noise photons are filtered using the signal confidence flag from ATL03 data and local point statistics. Second, we classify ground photons based on photon height percentiles. Third, tree and grass photons are classified based on the number of neighbors. We validated tree heights with field measurements (n = 55), finding a root-mean-square error (RMSE) of 1.82 m using SVDA, GEDI Level 2A (Geolocated Elevation and Height Metrics product): 1.33 m, and ATL08: 5.59 m. Our results indicate that the SVDA is effective in identifying canopy photons in savanna ecosystems, where ATL08 performs poorly. We further identify seasonal vegetation height changes with an emphasis on vegetation below 3 m; widespread height changes in this class from two wet-dry cycles show maximum seasonal changes of 1 m, possibly related to seasonal grass-height differences. Our study shows the difficulties of vegetation measurements in savanna ecosystems but provides the first estimates of seasonal biomass changes.
Cosmic-ray neutron sensing (CRNS) has become an effective method to measure soil moisture at a horizontal scale of hundreds of metres and a depth of decimetres. Recent studies proposed operating CRNS in a network with overlapping footprints in order to cover root-zone water dynamics at the small catchment scale and, at the same time, to represent spatial heterogeneity. In a joint field campaign from September to November 2020 (JFC-2020), five German research institutions deployed 15 CRNS sensors in the 0.4 km2 Wüstebach catchment (Eifel mountains, Germany). The catchment is dominantly forested (but includes a substantial fraction of open vegetation) and features a topographically distinct catchment boundary. In addition to the dense CRNS coverage, the campaign featured a unique combination of additional instruments and techniques: hydro-gravimetry (to detect water storage dynamics also below the root zone); ground-based and, for the first time, airborne CRNS roving; an extensive wireless soil sensor network, supplemented by manual measurements; and six weighable lysimeters. Together with comprehensive data from the long-term local research infrastructure, the published data set (available at https://doi.org/10.23728/b2share.756ca0485800474e9dc7f5949c63b872; Heistermann et al., 2022) will be a valuable asset in various research contexts: to advance the retrieval of landscape water storage from CRNS, wireless soil sensor networks, or hydrogravimetry; to identify scale-specific combinations of sensors and methods to represent soil moisture variability; to improve the understanding and simulation of land–atmosphere exchange as well as hydrological and hydrogeological processes at the hillslope and the catchment scale; and to support the retrieval of soil water content from airborne and spaceborne remote sensing platforms.
The Arctic is greatly affected by climate change. Increasing air temperatures drive permafrost thaw and an increase in coastal erosion and river discharge. This results in a greater input of sediment and organic matter into nearshore waters, impacting ecosystems by reducing light transmission through the water column and altering biogeochemistry. This potentially results in impacts on the subsistence economy of local people as well as the climate due to the transformation of suspended organic matter into greenhouse gases. Even though the impacts of increased suspended sediment concentrations and turbidity in the Arctic nearshore zone are well-studied, the mechanisms underpinning this increase are largely unknown. Wave energy and tides drive the level of turbidity in the temperate and tropical parts of the world, and this is generally assumed to also be the case in the Arctic. However, the tidal range is considerably lower in the Arctic, and processes related to the occurrence of permafrost have the potential to greatly contribute to nearshore turbidity. In this study, we use high-resolution satellite imagery alongside in situ and ERA5 reanalysis data of ocean and climate variables in order to identify the drivers of nearshore turbidity, along with its seasonality in the nearshore waters of Herschel Island Qikiqtaruk, in the western Canadian Arctic. Nearshore turbidity correlates well to wind direction, wind speed, significant wave height, and wave period. Nearshore turbidity is superiorly correlated to wind speed at the Beaufort Shelf compared to in situ measurements at Herschel Island Qikiqtaruk, showing that nearshore turbidity, albeit being of limited spatial extent, is influenced by large-scale weather and ocean phenomenons. We show that, in contrast to the temperate and tropical ocean, freshly eroded material is the predominant driver of nearshore turbidity in the Arctic, rather than resuspension, which is caused by the vulnerability of permafrost coasts to thermo-erosion.
The Arctic is greatly impacted by climate change. The increase in air temperature drives the thawing of permafrost and an increase in coastal erosion and river discharge. This leads to a greater input of sediment and organic matter into coastal waters, which substantially impacts the ecosystems by reducing light transmission through the water column and altering the biogeochemistry, but also the subsistence economy of local people, and changes in climate because of the transformation of organic matter into greenhouse gases. Yet, the quantification of suspended sediment in Arctic coastal and nearshore waters remains unsatisfactory due to the absence of dedicated algorithms to resolve the high loads occurring in the close vicinity of the shoreline. In this study we present the Arctic Nearshore Turbidity Algorithm (ANTA), the first reflectance-turbidity relationship specifically targeted towards Arctic nearshore waters that is tuned with in-situ measurements from the nearshore waters of Herschel Island Qikiqtaruk in the western Canadian Arctic. A semi-empirical model was calibrated for several relevant sensors in ocean color remote sensing, including MODIS, Sentinel 3 (OLCI), Landsat 8 (OLI), and Sentinel 2 (MSI), as well as the older Landsat sensors TM and ETM+. The ANTA performed better with Landsat 8 than with Sentinel 2 and Sentinel 3. The application of the ANTA to Sentinel 2 imagery that matches in-situ turbidity samples taken in Adventfjorden, Svalbard, shows transferability to nearshore areas beyond Herschel Island Qikiqtaruk.
Soziale Medien sind ein wesentlicher Bestandteil des Alltags von Schüler*innen und gleichzeitig zunehmend wichtig in Wirtschaft, Politik und Wissenschaft. Am Beispiel von Twitter zeigt dieser Beitrag, dass soziale Medien im Unterricht auch für die Beantwortung geographischer Fragestellungen verwendet werden können. Hierfür eignen sich Twitter-Daten aufgrund ihrer Georeferenzierung und weiterer interessanter Inhalte besonders. Der Beitrag gibt einen Überblick über die Verwendung von Twitter für sozialwissenschaftliche und humangeographische Fragestellungen und reflektiert die Nutzung von Twitter im Unterricht. Für die Unterrichtspraxis werden Beispiele zu den Themen Braunkohle, Flutereignisse und Raumwahrnehmungen sowie Anleitungen zur Auswertung, Anwendung und Reflexion von Twitter-Analysen vorgestellt.
Quantitative geomorphic research depends on accurate topographic data often collected via remote sensing. Lidar, and photogrammetric methods like structure-from-motion, provide the highest quality data for generating digital elevation models (DEMs). Unfortunately, these data are restricted to relatively small areas, and may be expensive or time-consuming to collect. Global and near-global DEMs with 1 arcsec (∼30 m) ground sampling from spaceborne radar and optical sensors offer an alternative gridded, continuous surface at the cost of resolution and accuracy. Accuracy is typically defined with respect to external datasets, often, but not always, in the form of point or profile measurements from sources like differential Global Navigation Satellite System (GNSS), spaceborne lidar (e.g., ICESat), and other geodetic measurements. Vertical point or profile accuracy metrics can miss the pixel-to-pixel variability (sometimes called DEM noise) that is unrelated to true topographic signal, but rather sensor-, orbital-, and/or processing-related artifacts. This is most concerning in selecting a DEM for geomorphic analysis, as this variability can affect derivatives of elevation (e.g., slope and curvature) and impact flow routing. We use (near) global DEMs at 1 arcsec resolution (SRTM, ASTER, ALOS, TanDEM-X, and the recently released Copernicus) and develop new internal accuracy metrics to assess inter-pixel variability without reference data. Our study area is in the arid, steep Central Andes, and is nearly vegetation-free, creating ideal conditions for remote sensing of the bare-earth surface. We use a novel hillshade-filtering approach to detrend long-wavelength topographic signals and accentuate short-wavelength variability. Fourier transformations of the spatial signal to the frequency domain allows us to quantify: 1) artifacts in the un-projected 1 arcsec DEMs at wavelengths greater than the Nyquist (twice the nominal resolution, so > 2 arcsec); and 2) the relative variance of adjacent pixels in DEMs resampled to 30-m resolution (UTM projected). We translate results into their impact on hillslope and channel slope calculations, and we highlight the quality of the five DEMs. We find that the Copernicus DEM, which is based on a carefully edited commercial version of the TanDEM-X, provides the highest quality landscape representation, and should become the preferred DEM for topographic analysis in areas without sufficient coverage of higher-quality local DEMs.
The intensification of Northern Hemisphere glaciations at the end of the Pliocene epoch marks one of the most substantial climatic shifts of the Cenozoic. Despite global cooling, sea surface temperatures in the high latitude North Atlantic Ocean rose between 2.9–2.7 million years ago. Here we present sedimentary geochemical proxy data from the Gulf of Cadiz to reconstruct the variability of Mediterranean Outflow Water, an important heat source to the North Atlantic. We find evidence for enhanced production of Mediterranean Outflow from the mid-Pliocene to the late Pliocene which we infer could have driven a sub-surface heat channel into the high-latitude North Atlantic. We then use Earth System Models to constrain the impact of enhanced Mediterranean Outflow production on the northward heat transport in the North Atlantic. In accord with the proxy data, the numerical model results support the formation of a sub-surface channel that pumped heat from the subtropics into the high latitude North Atlantic. We further suggest that this mechanism could have delayed ice sheet growth at the end of the Pliocene.
Sediment archives in the terrestrial and marine realm are regularly analyzed to infer changes in climate, tectonic, or anthropogenic boundary conditions of the past. However, contradictory observations have been made regarding whether short period events are faithfully preserved in stratigraphic archives; for instance, in marine sediments offshore large river systems. On the one hand, short period events are hypothesized to be non-detectable in the signature of terrestrially derived sediments due to buffering during sediment transport along large river systems. On the other hand, several studies have detected signals of short period events in marine records offshore large river systems. We propose that this apparent discrepancy is related to the lack of a differentiation between different types of signals and the lack of distinction between river response times and signal propagation times. In this review, we (1) expand the definition of the term ‘signal’ and group signals in sub-categories related to hydraulic grain size characteristics, (2) clarify the different types of ‘times’ and suggest a precise and consistent terminology for future use, and (3) compile and discuss factors influencing the times of signal transfer along sediment routing systems and how those times vary with hydraulic grain size characteristics. Unraveling different types of signals and distinctive time periods related to signal propagation addresses the discrepancies mentioned above and allows a more comprehensive exploration of event preservation in stratigraphy – a prerequisite for reliable environmental reconstructions from terrestrially derived sedimentary records.
Controversy over the plate tectonic affinity and evolution of the Saxon granulites in a two- or multi-plate setting during inter- or intracontinental collision makes the Saxon Granulite Massif a key area for the understanding of the Palaeozoic Variscan orogeny. The massif is a large dome structure in which tectonic slivers of metapelite and metaophiolite units occur along a shear zone separating a diapir-like body of high-Pgranulite below from low-Pmetasedimentary rocks above. Each of the upper structural units records a different metamorphic evolution until its assembly with the exhuming granulite body. New age and petrologic data suggest that the metaophiolites developed from early Cambrian protoliths during high-Pamphibolite facies metamorphism in the mid- to late-Devonian and thermal overprinting by the exhuming hot granulite body in the early Carboniferous. A correlation of new Ar-Ar biotite ages with publishedP-T-tdata for the granulites implies that exhumation and cooling of the granulite body occurred at average rates of similar to 8 mm/year and similar to 80 degrees C/Ma, with a drop in exhumation rate from similar to 20 to similar to 2.5 mm/year and a slight rise in cooling rate between early and late stages of exhumation. A time lag ofc. 2 Ma between cooling through the closure temperatures for argon diffusion in hornblende and biotite indicates a cooling rate of 90 degrees C/Ma when all units had assembled into the massif. A two-plate model of the Variscan orogeny in which the above evolution is related to a short-lived intra-Gondwana subduction zone conflicts with the oceanic affinity of the metaophiolites and the timescale ofc. 50 Ma for the metamorphism. Alternative models focusing on the internal Variscan belt assume distinctly different material paths through the lower or upper crust for strikingly similar granulite massifs. An earlier proposed model of bilateral subduction below the internal Variscan belt may solve this problem.
Graphite forms the endpoint for organic carbon metamorphism; it is extremely resilient to physical, biological and chemical degradation. Carbonaceous materials (CM) contained within sediments, collected across Taiwan and from the Gaoping submarine canyon, were analyzed using Raman spectroscopy to determine the crystallinity. This allowed the erosional and orogenic movements of petrogenic organic carbon (OCpetro) during the Taiwanese orogeny to be deduced. After automatically fitting and classifying spectra, the distribution of four groups of CM within the sediments provides evidence that many forms of OCpetro have survived at least one previous cycle of erosion, transport and burial before forming rocks in the Western Foothills of the island. There is extensive detrital graphite present in rocks that have not experienced high-grade metamorphism, and graphite flakes are also found in recently deposited marine sediments off Taiwan. The tectonic and geological history of the island shows that these graphite flakes must have survived at least three episodes of recycling. Therefore, transformation to graphite during burial and orogeny is a mechanism for stabilizing organic carbon over geological time, removing biospheric carbon from the active carbon cycle and protecting it from oxidation during future erosion events.
Holocene temperature proxy records are commonly used in quantitative synthesis and model-data comparisons. However, comparing correlations between time series from records collected in proximity to one another with the expected correlations based on climate model simulations indicates either regional or noisy climate signals in Holocene temperature proxy records. In this study, we evaluate the consistency of spatial correlations present in Holocene proxy records with those found in data from the Last Glacial Maximum (LGM). Specifically, we predict correlations expected in LGM proxy records if the only difference to Holocene correlations would be due to more time uncertainty and more climate variability in the LGM. We compare this simple prediction to the actual correlation structure in the LGM proxy records. We found that time series data of ice-core stable isotope records and planktonic foraminifera Mg/Ca ratios were consistent between the Holocene and LGM periods, while time series of Uk'37 proxy records were not as we found no correlation between nearby LGM records. Our results support the finding of highly regional or noisy marine proxy records in the compilation analysed here and suggest the need for further studies on the role of climate proxies and the processes of climate signal recording and preservation.
Marked along-strike changes in stratigraphy, mountain belt morphology, basement exhumation, and deformation styles characterize the Andean retroarc; these changes have previously been related to spatiotemporal variations in the subduction angle. We modeled new apatite fission track and apatite (U-Th-Sm)/He data from nine ranges located between 26 degrees S and 28 degrees S. Using new and previously published data, we constructed a Cretaceous to Pliocene paleogeographic model that delineates a four-stage tectonic evolution: extensional tectonics during the Cretaceous (120-75 Ma), the formation of a broken foreland basin between 55 and 30 Ma, reheating due to burial beneath sedimentary rocks (18-13 Ma), and deformation, exhumation, and surface uplift during the Late Miocene and the Pliocene (13-3 Ma). Our model highlights how preexisting upper plate structures control the deformation patterns of broken foreland basins. Because retroarc deformation predates flat-slab subduction, we propose that slab anchoring may have been the precursor of Eocene-Oligocene compression in the Andean retroarc. Our model challenges models which consider broken foreland basins and retroarc deformation in the NW Argentinian Andes to be directly related to Miocene flat subduction.
The intensification of Northern Hemisphere glaciation (iNHG) and uplift of the Tibetan Plateau have been argued to be among the main drivers of climate change in midlatitude Central Asia during the Pliocene/Pleistocene. While most proxy records that support this hypothesis are from regions outside the Tibetan Plateau (such as from the Chinese Loess Plateau), detailed paleoclimatic information for the plateau itself during that time has yet remained elusive. Here we present a temporally highly resolved (similar to 500 years) sedimentological record from the Qaidam Basin situated on the northeastern Tibetan Plateau that shows pronounced glacial-interglacial climate variability during the interval from 2.7 to 2.1 Ma. Glacial (interglacial) intervals are generally characterized by coarser (finer) grain size, minima (maxima) in organic matter content, and maxima (minima) in carbonate content. Comparison of our results with Earth's orbital parameters and proxy records from the Chinese Loess Plateau suggests that the observed climate fluctuations were mainly driven by changes in the Siberian High/East Asian winter monsoon system as a response to the iNHG. They are further proposed to be enhanced by the topography of the Tibetan Plateau and its impact on the position and intensity of the westerlies.
The Cluster mission has produced a large data set of electron flux measurements in the Earth's magnetosphere since its launch in late 2000. Electron fluxes are measured using Research with Adaptive Particle Imaging Detector (RAPID)/Imaging Electron Spectrometer (IES) detector as a function of energy, pitch angle, spacecraft position, and time. However, no adiabatic invariants have been calculated for Cluster so far. In this paper we present a step-by-step guide to calculations of adiabatic invariants and conversion of the electron flux to phase space density (PSD) in these coordinates. The electron flux is measured in two RAPID/IES energy channels providing pitch angle distribution at energies 39.2-50.5 and 68.1-94.5 keV in nominal mode since 2004. A fitting method allows to expand the conversion of the differential fluxes to the range from 40 to 150 keV. Best data coverage for phase space density in adiabatic invariant coordinates can be obtained for values of second adiabatic invariant, K, similar to 10(2), and values of the first adiabatic invariant mu in the range approximate to 5-20 MeV/G. Furthermore, we describe the production of a new data product "LSTAR," equivalent to the third adiabatic invariant, available through the Cluster Science Archive for years 2001-2018 with 1-min resolution. The produced data set adds to the availability of observations in Earth's radiation belts region and can be used for long-term statistical purposes.
Issue Despite their rather similar climatic conditions, eastern Eurasia and northern North America are largely covered by different plant functional types (deciduous or evergreen boreal forest) composed of larch or pine, spruce and fir, respectively. I propose that these deciduous and evergreen boreal forests represent alternative quasi-stable states, triggered by their different northern tree refugia that reflect the different environmental conditions experienced during the Last Glacial. Evidence This view is supported by palaeoecological and environmental evidence. Once established, Asian larch forests are likely to have stabilized through a complex vegetation-fire-permafrost soil-climate feedback system. Conclusion With respect to future forest developments, this implies that Asian larch forests are likely to be governed by long-term trajectories and are therefore largely resistant to natural climate variability on time-scales shorter than millennia. The effects of regional human impact and anthropogenic global warming might, however, cause certain stability thresholds to be crossed, meaning that irreversible transitions occur and resulting in marked consequences for ecosystem services on these human-relevant time-scales.
Millennial-scale palaeoclimate variability has been documented in various terrestrial and marine palaeoclimate proxy records throughout the Northern Hemisphere for the last glacial cycle. Its clear expression and rapid shifts between different states of climate (Greenland Interstadials and Stadials) represents a correlation tool beyond the resolution of e.g. luminescence dating, especially relevant for terrestrial deposits. Usually, comparison of terrestrial proxy datasets and the Greenland ice cores indicates a complex expression of millennial-scale climate variability as recorded in terrestrial geoarchives including loess. Loess is the most widespread terrestrial geoarchive of the Quaternary and especially widespread over Eurasia. However, loess often records a smoothed representation of millennial-scale variability without all fidelity when compared to the Greenland data, this being a relevant limiting feature in integrating loess with other palaeoclimate records. To better understand the loess proxy-response to millennial-scale climate variability, we simulate a proxy signal smoothing by natural processes through application of low-pass filters of delta O-18 data from Greenland, a high-resolution palaeoclimate reference record, alongside speleothem isotope records from the Black Sea-Mediterranean region. We show that low-pass filters represent rather simple models for better constraining the expression of millennial-scale climate variability in low sedimentation environments, and in sediments where proxy-response signals are most likely affected by natural smoothing (by e.g. bioturbation). Interestingly, smoothed datasets from Greenland and the Black Sea-Mediterranean region are most similar in the last similar to 15 ka and between similar to 50-30 ka. Between similar to 30-15 ka, roughly corresponding to the Last Glacial Maximum and the deglaciation, the records show dissimilarities, challenging the construction of robust correlative time-scales in this age range. From our analysis it becomes apparent that patterns of palaeoclimate signals in loess-palaeosol sequences often might be better explained by smoothed Greenland reference data than the original high-resolution Greenland dataset, or other reference data. This opens the possibility to better assess the temporal resolution and palaeoclimate potential of loess-palaeosol sequences in recording supra-regional climate patterns, as well as to securely integrate loess with other chronologically better-resolved palaeoclimate records.
Observations of rift and rifted margin architecture suggest that significant spatial and temporal structural heterogeneity develops during the multiphase evolution of continental rifting. Inheritance is often invoked to explain this heterogeneity, such as preexisting anisotropies in rock composition, rheology, and deformation. Here, we use high-resolution 3-D thermal-mechanical numerical models of continental extension to demonstrate that rift-parallel heterogeneity may develop solely through fault network evolution during the transition from distributed to localized deformation. In our models, the initial phase of distributed normal faulting is seeded through randomized initial strength perturbations in an otherwise laterally homogeneous lithosphere extending at a constant rate. Continued extension localizes deformation onto lithosphere-scale faults, which are laterally offset by tens of km and discontinuous along-strike. These results demonstrate that rift- and margin-parallel heterogeneity of large-scale fault patterns may in-part be a natural byproduct of fault network coalescence.
Temporal variation of natural light sources such as airglow limits the ability of night light sensors to detect changes in small sources of artificial light (such as villages). This study presents a method for correcting for this effect globally, using the satellite radiance detected from regions without artificial light emissions. We developed a routine to define an approximate grid of locations worldwide that do not have regular light emission. We apply this method with a 5 degree equally spaced global grid (total of 2016 individual locations), using data from the Visible Infrared Imaging Radiometer Suite (VIIRS) Day-Night Band (DNB). This code could easily be adapted for other future global sensors. The correction reduces the standard deviation of data in the Earth Observation Group monthly DNB composites by almost a factor of two. The code and datasets presented here are available under an open license by GFZ Data Services, and are implemented in the Radiance Light Trends web application.
Atmospheric water vapour content is a key variable that controls the development of deep convective storms and rainfall extremes over the central Andes. Direct measurements of water vapour are challenging; however, recent developments in microwave processing allow the use of phase delays from L-band radar to measure the water vapour content throughout the atmosphere: Global Navigation Satellite System (GNSS)-based integrated water vapour (IWV) monitoring shows promising results to measure vertically integrated water vapour at high temporal resolutions. Previous works also identified convective available potential energy (CAPE) as a key climatic variable for the formation of deep convective storms and rainfall in the central Andes. Our analysis relies on GNSS data from the Argentine Continuous Satellite Monitoring Network, Red Argentina de Monitoreo Satelital Continuo (RAMSAC) network from 1999 to 2013. CAPE is derived from version 2.0 of the ECMWF’s (European Centre for Medium-Range Weather Forecasts) Re-Analysis (ERA-interim) and rainfall from the TRMM (Tropical Rainfall Measuring Mission) product. In this study, we first analyse the rainfall characteristics of two GNSS-IWV stations by comparing their complementary cumulative distribution function (CCDF). Second, we separately derive the relation between rainfall vs. CAPE and GNSS-IWV. Based on our distribution fitting analysis, we observe an exponential relation of rainfall to GNSS-IWV. In contrast, we report a power-law relationship between the daily mean value of rainfall and CAPE at the GNSS-IWV station locations in the eastern central Andes that is close to the theoretical relationship based on parcel theory. Third, we generate a joint regression model through a multivariable regression analysis using CAPE and GNSS-IWV to explain the contribution of both variables in the presence of each other to extreme rainfall during the austral summer season. We found that rainfall can be characterised with a higher statistical significance for higher rainfall quantiles, e.g., the 0.9 quantile based on goodness-of-fit criterion for quantile regression. We observed different contributions of CAPE and GNSS-IWV to rainfall for each station for the 0.9 quantile. Fourth, we identify the temporal relation between extreme rainfall (the 90th, 95th, and 99th percentiles) and both GNSS-IWV and CAPE at 6 h time steps. We observed an increase before the rainfall event and at the time of peak rainfall—both for GNSS-integrated water vapour and CAPE. We show higher values of CAPE and GNSS-IWV for higher rainfall percentiles (99th and 95th percentiles) compared to the 90th percentile at a 6-h temporal scale. Based on our correlation analyses and the dynamics of the time series, we show that both GNSS-IWV and CAPE had comparable magnitudes, and we argue to consider both climatic variables when investigating their effect on rainfall extremes.
Initiation of subduction following the impingement of a hot buoyant mantle plume is one of the few scenarios that allow breaking the lithosphere and recycling a stagnant lid without requiring any preexisting weak zones. Here, we investigate factors controlling the number and shape of retreating subducting slabs formed by plume-lithosphere interaction. Using 3-D thermomechanical models we show that the deformation regime, which defines formation of single-slab or multi-slab subduction, depends on several parameters such as age of oceanic lithosphere, thickness of the crust and large-scale lithospheric extension rate. Our model results indicate that on present-day Earth multi-slab plume-induced subduction is initiated only if the oceanic lithosphere is relatively young (<30-40 Myr, but >10 Myr), and the crust has a typical thickness of 8 km. In turn, development of single-slab subduction is facilitated by older lithosphere and pre-imposed extensional stresses. In early Earth, plume-lithosphere interaction could have led to formation of either episodic short-lived circular subduction when the oceanic lithosphere was young or to multi-slab subduction when the lithosphere was old.
Climatic change alters the frequency and intensity of natural hazards. In order to assess potential future changes in flood seasonality in the Rhine River Basin, we analyse changes in streamflow, snowmelt, precipitation, and evapotranspiration at 1.5, 2.0 and 3.0 ◦C global warming levels. The mesoscale Hydrological Model (mHM) forced with an ensemble of climate projection scenarios (five general circulation models under three representative concentration pathways) is used to simulate the present and future climate conditions of both, pluvial and nival hydrological regimes. Our results indicate that the interplay between changes in snowmelt- and rainfall-driven runoff is crucial to understand changes in streamflow maxima in the Rhine River. Climate projections suggest that future changes in flood characteristics in the entire Rhine River are controlled by both, more intense precipitation events and diminishing snow packs. The nature of this interplay defines the type of change in runoff peaks. On the sub-basin level (the Moselle River), more intense rainfall during winter is mostly counterbalanced by reduced snowmelt contribution to the streamflow. In the High Rhine (gauge at Basel), the strongest increases in streamflow maxima show up during winter, when strong increases in liquid precipitation intensity encounter almost unchanged snowmelt-driven runoff. The analysis of snowmelt events suggests that at no point in time during the snowmelt season, a warming climate results in an increase in the risk of snowmelt-driven flooding. We do not find indications of a transient merging of pluvial and nival floods due to climate warming.
Glacial lakes in the Hindu Kush–Karakoram–Himalayas–Nyainqentanglha (HKKHN) region have grown rapidly in number and area in past decades, and some dozens have drained in catastrophic glacial lake outburst floods (GLOFs). Estimating regional susceptibility of glacial lakes has largely relied on qualitative assessments by experts, thus motivating a more systematic and quantitative appraisal. Before the backdrop of current climate-change projections and the potential of elevation-dependent warming, an objective and regionally consistent assessment is urgently needed. We use an inventory of 3390 moraine-dammed lakes and their documented outburst history in the past four decades to test whether elevation, lake area and its rate of change, glacier-mass balance, and monsoonality are useful inputs to a probabilistic classification model. We implement these candidate predictors in four Bayesian multi-level logistic regression models to estimate the posterior susceptibility to GLOFs. We find that mostly larger lakes have been more prone to GLOFs in the past four decades regardless of the elevation band in which they occurred. We also find that including the regional average glacier-mass balance improves the model classification. In contrast, changes in lake area and monsoonality play ambiguous roles. Our study provides first quantitative evidence that GLOF susceptibility in the HKKHN scales with lake area, though less so with its dynamics. Our probabilistic prognoses offer improvement compared to a random classification based on average GLOF frequency. Yet they also reveal some major uncertainties that have remained largely unquantified previously and that challenge the applicability of single models. Ensembles of multiple models could be a viable alternative for more accurately classifying the susceptibility of moraine-dammed lakes to GLOFs.
Rotational motions play a key role in measuring seismic wavefield properties. Using newly developed portable rotational instruments, it is now possible to directly measure rotational motions in a broad frequency range. Here, we investigated the instrumental self-noise and data quality in a huddle test in Fürstenfeldbruck, Germany, in August 2019. We compare the data from six rotational and three translational sensors. We studied the recorded signals using correlation, coherence analysis, and probabilistic power spectral densities. We sorted the coherent noise into five groups with respect to the similarities in frequency content and shape of the signals. These coherent noises were most likely caused by electrical devices, the dehumidifier system in the building, humans, and natural sources such as wind. We calculated self-noise levels through probabilistic power spectral densities and by applying the Sleeman method, a three-sensor method. Our results from both methods indicate that self-noise levels are stable between 0.5 and 40 Hz. Furthermore, we recorded the 29 August 2019 ML 3.4 Dettingen earthquake. The calculated source directions are found to be realistic for all sensors in comparison to the real back azimuth. We conclude that the five tested blueSeis-3A rotational sensors, when compared with respect to coherent noise, self-noise, and source direction, provide reliable and consistent results. Hence, field experiments with single rotational sensors can be undertaken.
Models for the predictions of monetary losses from floods mainly blend data deemed to represent a single flood type and region. Moreover, these approaches largely ignore indicators of preparedness and how predictors may vary between regions and events, challenging the transferability of flood loss models. We use a flood loss database of 1812 German flood-affected households to explore how Bayesian multilevel models can estimate normalised flood damage stratified by event, region, or flood process type. Multilevel models acknowledge natural groups in the data and allow each group to learn from others. We obtain posterior estimates that differ between flood types, with credibly varying influences of water depth, contamination, duration, implementation of property-level precautionary measures, insurance, and previous flood experience; these influences overlap across most events or regions, however. We infer that the underlying damaging processes of distinct flood types deserve further attention. Each reported flood loss and affected region involved mixed flood types, likely explaining the uncertainty in the coefficients. Our results emphasise the need to consider flood types as an important step towards applying flood loss models elsewhere. We argue that failing to do so may unduly generalise the model and systematically bias loss estimations from empirical data.
A tale of shifting relations
(2021)
Understanding the dynamics between the East Asian summer (EASM) and winter monsoon (EAWM) is needed to predict their variability under future global warming scenarios. Here, we investigate the relationship between EASM and EAWM as well as the mechanisms driving their variability during the last 10,000 years by stacking marine and terrestrial (non-speleothem) proxy records from the East Asian realm. This provides a regional and proxy independent signal for both monsoonal systems. The respective signal was subsequently analysed using a linear regression model. We find that the phase relationship between EASM and EAWM is not time-constant and significantly depends on orbital configuration changes. In addition, changes in the Atlantic Meridional Overturning circulation, Arctic sea-ice coverage, El Niño-Southern Oscillation and Sun Spot numbers contributed to millennial scale changes in the EASM and EAWM during the Holocene. We also argue that the bulk signal of monsoonal activity captured by the stacked non-speleothem proxy records supports the previously argued bias of speleothem climatic archives to moisture source changes and/or seasonality.
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.
Sediment Transit Time and Floodplain Storage Dynamics in Alluvial Rivers Revealed by Meteoric 10Be
(2020)
Quantifying the time scales of sediment transport and storage through river systems is fundamental for understanding weathering processes, biogeochemical cycling, and improving watershed management, but measuring sediment transit time is challenging. Here we provide the first systematic test of measuring cosmogenic meteoric Beryllium-10 (10Bem) in the sediment load of a large alluvial river to quantify sediment transit times. We take advantage of a natural experiment in the Rio Bermejo, a lowland alluvial river traversing the east Andean foreland basin in northern Argentina. This river has no tributaries along its trunk channel for nearly 1,300 km downstream from the mountain front. We sampled suspended sediment depth profiles along the channel and measured the concentrations of 10Bem in the chemically extracted grain coatings. We calculated depth-integrated 10Bem concentrations using sediment flux data and found that 10Bem concentrations increase 230% from upstream to downstream, indicating a mean total sediment transit time of 8.4 ± 2.2 kyr. Bulk sediment budget-based estimates of channel belt and fan storage times suggest that the 10Bem tracer records mixing of old and young sediment reservoirs. On a reach scale, 10Bem transit times are shorter where the channel is braided and superelevated above the floodplain, and longer where the channel is incised and meandering, suggesting that transit time is controlled by channel morphodynamics. This is the first systematic application of 10Bem as a sediment transit time tracer and highlights the method's potential for inferring sediment routing and storage dynamics in large river systems.
In this study, we investigated the scale sizes of equatorial plasma irregularities (EPIs) using measurements from the Swarm satellites during its early mission and final constellation phases. We found that with longitudinal separation between Swarm satellites larger than 0.4°, no significant correlation was found any more. This result suggests that EPI structures include plasma density scale sizes less than 44 km in the zonal direction. During the Swarm earlier mission phase, clearly better EPI correlations are obtained in the northern hemisphere, implying more fragmented irregularities in the southern hemisphere where the ambient magnetic field is low. The previously reported inverted-C shell structure of EPIs is generally confirmed by the Swarm observations in the northern hemisphere, but with various tilt angles. From the Swarm spacecrafts with zonal separations of about 150 km, we conclude that larger zonal scale sizes of irregularities exist during the early evening hours (around 1900 LT).
The Cheb Basin (CZ) is a shallow Neogene intracontinental basin filled with fluvial and lacustrine sediments that is located in the western part of the Eger Rift. The basin is situated in a seismically active area and is characterized by diffuse degassing of mantle-derived CO2 in mofette fields. The Hartousov mofette field shows a daily CO2 flux of 23-97 tons of CO2 released over an area of 0.35 km(2) and a soil gas concentration of up to 100% CO2. The present study aims to explore the geo-bio interactions provoked by the influence of elevated CO2 concentrations on the geochemistry and microbial community of soils and sediments. To sample the strata, two 3-m cores were recovered. One core stems from the center of the degassing structure, whereas the other core was taken 8 m from the ENE and served as an undisturbed reference site. The sites were compared regarding their geochemical features, microbial abundances, and microbial community structures. The mofette site is characterized by a low pH and high TOC/sulfate contents. Striking differences in the microbial community highlight the substantial impact of elevated CO2 concentrations and their associated side effects on microbial processes. The abundance of microbes did not show a typical decrease with depth, indicating that the uprising CO2-rich fluid provides sufficient substrate for chemolithoautotrophic anaerobic microorganisms. Illumine MiSeq sequencing of the 16S rRNA genes and multivariate statistics reveals that the pH strongly influences microbial composition and explains around 38.7% of the variance at the mofette site and 22.4% of the variance between the mofette site and the undisturbed reference site. Accordingly, acidophilic microorganisms (e.g., OTUs assigned to Acidobacteriaceae and Acidithiobacillus) displayed a much higher relative abundance at the mofette site than at the reference site. The microbial community at the mofette site is characterized by a high relative abundance of methanogens and taxa involved in sulfur cycling. The present study provides intriguing insights into microbial life and geo-bio interactions in an active seismic region dominated by emanating mantle-derived CO2-rich fluids, and thereby builds the basis for further studies, e.g., focusing on the functional repertoire of the communities. However, it remains open if the observed patterns can be generalized for different time-points or sites.
Quantifying rock weakening due to decreasing calcite mineral content by numerical simulations
(2018)
The quantification of changes in geomechanical properties due to chemical reactions is of paramount importance for geological subsurface utilisation, since mineral dissolution generally reduces rock stiffness. In the present study, the effective elastic moduli of two digital rock samples, the Fontainebleau and Bentheim sandstones, are numerically determined based on micro-CT images. Reduction in rock stiffness due to the dissolution of 10% calcite cement by volume out of the pore network is quantified for three synthetic spatial calcite distributions (coating, partial filling and random) using representative sub-cubes derived from the digital rock samples. Due to the reduced calcite content, bulk and shear moduli decrease by 34% and 38% in maximum, respectively. Total porosity is clearly the dominant parameter, while spatial calcite distribution has a minor impact, except for a randomly chosen cement distribution within the pore network. Moreover, applying an initial stiffness reduced by 47% for the calcite cement results only in a slightly weaker mechanical behaviour. Using the quantitative approach introduced here substantially improves the accuracy of predictions in elastic rock properties compared to general analytical methods, and further enables quantification of uncertainties related to spatial variations in porosity and mineral distribution.
LiCSBAS
(2020)
For the past five years, the 2-satellite Sentinel-1 constellation has provided abundant and useful Synthetic Aperture Radar (SAR) data, which have the potential to reveal global ground surface deformation at high spatial and temporal resolutions. However, for most users, fully exploiting the large amount of associated data is challenging, especially over wide areas. To help address this challenge, we have developed LiCSBAS, an open-source SAR interferometry (InSAR) time series analysis package that integrates with the automated Sentinel-1 InSAR processor (LiCSAR). LiCSBAS utilizes freely available LiCSAR products, and users can save processing time and disk space while obtaining the results of InSAR time series analysis. In the LiCSBAS processing scheme, interferograms with many unwrapping errors are automatically identified by loop closure and removed. Reliable time series and velocities are derived with the aid of masking using several noise indices. The easy implementation of atmospheric corrections to reduce noise is achieved with the Generic Atmospheric Correction Online Service for InSAR (GACOS). Using case studies in southern Tohoku and the Echigo Plain, Japan, we demonstrate that LiCSBAS applied to LiCSAR products can detect both large-scale (>100 km) and localized (~km) relative displacements with an accuracy of <1 cm/epoch and ~2 mm/yr. We detect displacements with different temporal characteristics, including linear, periodic, and episodic, in Niigata, Ojiya, and Sanjo City, respectively. LiCSBAS and LiCSAR products facilitate greater exploitation of globally available and abundant SAR datasets and enhance their applications for scientific research and societal benefit.
The variabilities of the semidiurnal solar and lunar tides of the equatorial electrojet (EEJ) are investigated during the 2003, 2006, 2009 and 2013 major sudden stratospheric warming (SSW) events in this study. For this purpose, ground-magnetometer recordings at the equatorial observatories in Huancayo and Fúquene are utilized. Results show a major enhancement in the amplitude of the EEJ semidiurnal lunar tide in each of the four warming events. The EEJ semidiurnal solar tidal amplitude shows an amplification prior to the onset of warmings, a reduction during the deceleration of the zonal mean zonal wind at 60∘ N and 10 hPa, and a second enhancement a few days after the peak reversal of the zonal mean zonal wind during all four SSWs. Results also reveal that the amplitude of the EEJ semidiurnal lunar tide becomes comparable or even greater than the amplitude of the EEJ semidiurnal solar tide during all these warming events. The present study also compares the EEJ semidiurnal solar and lunar tidal changes with the variability of the migrating semidiurnal solar (SW2) and lunar (M2) tides in neutral temperature and zonal wind obtained from numerical simulations at E-region heights. A better agreement between the enhancements of the EEJ semidiurnal lunar tide and the M2 tide is found in comparison with the enhancements of the EEJ semidiurnal solar tide and the SW2 tide in both the neutral temperature and zonal wind at the E-region altitudes.
Background. Metal recycling factories (MRFs) have developed rapidly in Nigeria as recycling policies have been increasingly embraced. These MRFs are point sources for introducing potentially toxic elements (PTEs) into environmental media. Objectives. The aim of this study was to determine the constituents (elemental and mineralogy) of the wastes (slag and particulate matter, (PM)) and soils around the MRFs and to determine the level of pollution within the area. Methods. Sixty samples (30 slag samples, 15 soil samples and 15 PM samples) were collected for this study. The soils, slag and PM samples were analyzed for elemental constituents using inductively coupled plasma optical emission spectrometry. Mineralogy of the PM was determined using scanning electron microscope-energy dispersive spectroscopy (SEM-EDS), and soil mineralogy was determined by an X-ray diffractometer (XRD). Results. The results of the soil analyses revealed the following concentrations for the selected metals in mg/kg include lead (Pb) (21.0-2399.0), zinc (Zn) (56.0-4188.0), copper (Cu) (10.0-1470.0), nickel (Ni) (6.0-215.0), chromium (Cr) (921.0-1737.0) and cadmium (Cd) (below detectable limit (Bdl)-18.1). For the slags the results were Pb (68.0-.333.0), Zn (1364.0-3062), Cu (119.0-1470.0), Ni (12.0-675.0), Cr (297-1737) and Cd (Bdl-15.8). The results in µg/g for the metal analysis in PM were Pb (4.6-160.0), Zn (18.0-471.0), Cu (2.5-11.0), Ni (0.8-4.2), and Cr (2.5-11.0), while Cd was undetected. The slags are currently utilized for filling the foundations of buildings and roads, providing additional pathways for the introduction of PTEs into the environment from the suspended materials generated from mechanical breakdown of the slags. Conclusions. The MRFs were found to have impacted the quality of environmental media through the introduction of PTEs, impairing soil quality, in addition to PM, which can have detrimental health consequences. Further studies on the health implications of these pollutants and their impacts on human health are needed. Competing Interests. The authors declare no competing financial interests
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.
The article describes the surface modification of 3D printed poly(lactic acid) (PLA) scaffolds with calcium phosphate (CP)/gelatin and CP/chitosan hybrid coating layers. The presence of gelatin or chitosan significantly enhances CP co-deposition and adhesion of the mineral layer on the PLA scaffolds. The hydrogel/CP coating layers are fairly thick and the mineral is a mixture of brushite, octacalcium phosphate, and hydroxyapatite. Mineral formation is uniform throughout the printed architectures and all steps (printing, hydrogel deposition, and mineralization) are in principle amenable to automatization. Overall, the process reported here therefore has a high application potential for the controlled synthesis of biomimetic coatings on polymeric biomaterials.
Changes in species' distributions are classically projected based on their climate envelopes. For Siberian forests, which have a tremendous significance for vegetation-climate feedbacks, this implies future shifts of each of the forest-forming larch (Larix) species to the north-east. However, in addition to abiotic factors, reliable projections must assess the role of historical biogeography and biotic interactions. Here, we use sedimentary ancient DNA and individual-based modelling to investigate the distribution of larch species and mitochondrial haplotypes through space and time across the treeline ecotone on the southern Taymyr peninsula, which at the same time presents a boundary area of two larch species. We find spatial and temporal patterns, which suggest that forest density is the most influential driver determining the precise distribution of species and mitochondrial haplotypes. This suggests a strong influence of competition on the species' range shifts. These findings imply possible climate change outcomes that are directly opposed to projections based purely on climate envelopes. Investigations of such fine-scale processes of biodiversity change through time are possible using paleoenvironmental DNA, which is available much more readily than visible fossils and can provide information at a level of resolution that is not reached in classical palaeoecology.
Dating growth strata and basin fill by combining 26Al/10Be burial dating and magnetostratigraphy
(2018)
Cosmogenic burial dating enables dating of coarse-grained, Pliocene-Pleistocene sedimentary units that are typically difficult to date with traditional methods, such as magnetostratigraphy. In the actively deforming western Tarim Basin in NW China, Pliocene-Pleistocene conglomerates were dated at eight sites, integrating Al-26/Be-10 burial dating with previously published magnetostratigraphic sections. These samples were collected from growth strata on the flanks of growing folds and from sedimentary units beneath active faults to place timing constraints on the initiation of deformation of structures within the basin and on shortening rates on active faults. These new basin-fill and growthstrata ages document the late Neogene and Quaternary growth of the Pamir and Tian Shan orogens between >5 and 1 Ma and delineate the eastward propagation of deformation at rates up to 115 km/m.y. and basinward growth of both mountain belts at rates up to 12 km/m.y.
Arctic tundra landscapes are composed of a complex mosaic of patterned ground features, varying in soil moisture, vegetation composition, and surface hydrology over small spatial scales (10-100 m). The importance of microtopography and associated geomorphic landforms in influencing ecosystem structure and function is well founded, however, spatial data products describing local to regional scale distribution of patterned ground or polygonal tundra geomorphology are largely unavailable. Thus, our understanding of local impacts on regional scale processes (e.g., carbon dynamics) may be limited. We produced two key spatiotemporal datasets spanning the Arctic Coastal Plain of northern Alaska (similar to 60,000 km(2)) to evaluate climate-geomorphological controls on arctic tundra productivity change, using (1) a novel 30m classification of polygonal tundra geomorphology and (2) decadal-trends in surface greenness using the Landsat archive (1999-2014). These datasets can be easily integrated and adapted in an array of local to regional applications such as (1) upscaling plot-level measurements (e.g., carbon/energy fluxes), (2) mapping of soils, vegetation, or permafrost, and/or (3) initializing ecosystem biogeochemistry, hydrology, and/or habitat modeling.
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.
High Mountain Asia (HMA) is dependent upon both the amount and timing of snow and glacier meltwater. Previous model studies and coarse resolution (0.25° × 0.25°, ∼25 km × 25 km) passive microwave assessments of trends in the volume and timing of snowfall, snowmelt, and glacier melt in HMA have identified key spatial and seasonal heterogeneities in the response of snow to changes in regional climate. Here we use recently developed, continuous, internally consistent, and high-resolution passive microwave data (3.125 km × 3.125 km, 1987–2016) from the special sensor microwave imager instrument family to refine and extend previous estimates of changes in the snow regime of HMA. We find an overall decline in snow volume across HMA; however, there exist spatially contiguous regions of increasing snow volume—particularly during the winter season in the Pamir, Karakoram, Hindu Kush, and Kunlun Shan. Detailed analysis of changes in snow-volume trends through time reveal a large step change from negative trends during the period 1987–1997, to much more positive trends across large regions of HMA during the periods 1997–2007 and 2007–2016. We also find that changes in high percentile monthly snow-water volume exhibit steeper trends than changes in low percentile snow-water volume, which suggests a reduction in the frequency of high snow-water volumes in much of HMA. Regions with positive snow-water storage trends generally correspond to regions of positive glacier mass balances.
Permafrost is warming in the northern high latitudes, inducing highly dynamic thaw-related permafrost disturbances across the terrestrial Arctic. Monitoring and tracking of permafrost disturbances is important as they impact surrounding landscapes, ecosystems and infrastructure. Remote sensing provides the means to detect, map, and quantify these changes homogeneously across large regions and time scales. Existing Landsat-based algorithms assess different types of disturbances with similar spatiotemporal requirements. However, Landsat-based analyses are restricted in northern high latitudes due to the long repeat interval and frequent clouds, in particular at Arctic coastal sites. We therefore propose to combine Landsat and Sentinel-2 data for enhanced data coverage and present a combined annual mosaic workflow, expanding currently available algorithms, such as LandTrendr, to achieve more reliable time series analysis. We exemplary test the workflow for twelve sites across the northern high latitudes in Siberia. We assessed the number of images and cloud-free pixels, the spatial mosaic coverage and the mosaic quality with spectral comparisons. The number of available images increased steadily from 1999 to 2019 but especially from 2016 onward with the addition of Sentinel-2 images. Consequently, we have an increased number of cloud-free pixels even under challenging environmental conditions, which then serve as the input to the mosaicking process. In a comparison of annual mosaics, the Landsat+Sentinel-2 mosaics always fully covered the study areas (99.9–100 %), while Landsat-only mosaics contained data-gaps in the same years, only reaching coverage percentages of 27.2 %, 58.1 %, and 69.7 % for Sobo Sise, East Taymyr, and Kurungnakh in 2017, respectively. The spectral comparison of Landsat image, Sentinel-2 image, and Landsat+Sentinel-2 mosaic showed high correlation between the input images and mosaic bands (e.g., for Kurungnakh 0.91–0.97 between Landsat and Landsat+Sentinel-2 mosaic and 0.92–0.98 between Sentinel-2 and Landsat+Sentinel-2 mosaic) across all twelve study sites, testifying good quality mosaic results. Our results show that especially the results for northern, coastal areas was substantially improved with the Landsat+Sentinel-2 mosaics. By combining Landsat and Sentinel-2 data we accomplished to create reliably high spatial resolution input mosaics for time series analyses. Our approach allows to apply a high temporal continuous time series analysis to northern high latitude permafrost regions for the first time, overcoming substantial data gaps, and assess permafrost disturbance dynamics on an annual scale across large regions with algorithms such as LandTrendr by deriving the location, timing and progression of permafrost thaw disturbances
Uplift in the broken Andean foreland of the Argentine Santa Bárbara System (SBS) is associated with the contractional reactivation of basement anisotropies, similar to those reported from the thick-skinned Cretaceous-Eocene Laramide province of North America. Fault scarps, deformed Quaternary deposits and landforms, disrupted drainage patterns, and medium-sized earthquakes within the SBS suggest that movement along these structures may be a recurring phenomenon, with yet to be defined repeat intervals and rupture lengths. In contrast to the Subandes thrust belt farther north, where eastward-migrating deformation has generated a well-defined thrust front, the SBS records spatiotemporally disparate deformation along structures that are only known to the first order. We present herein the results of geomorphic desktop analyses, structural field observations, and 2D electrical resistivity tomography and seismic-refraction tomography surveys and an interpretation of seismic reflection profiles across suspected fault scarps in the sedimentary basins adjacent to the Candelaria Range (CR) basement uplift, in the south-central part of the SBS. Our analysis in the CR piedmont areas reveals consistency between the results of near-surface electrical resistivity and seismic-refraction tomography surveys, the locations of prominent fault scarps, and structural geometries at greater depth imaged by seismic reflection data. We suggest that this deformation is driven by deep-seated blind thrusting beneath the CR and associated regional warping, while shortening involving Mesozoic and Cenozoic sedimentary strata in the adjacent basins was accommodated by layer-parallel folding and flexural-slip faults that cut through Quaternary landforms and deposits at the surface.
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.
More than 41% of the Earth’s land area is covered by permanent or seasonally arid dryland ecosystems. Global development and human activity have led to an increase in aridity, resulting in ecosystem degradation and desertification around the world. The objective of the present work was to investigate and compare the microbial community structure and geochemical characteristics of two geographically distinct saline pan sediments in the Kalahari Desert of southern Africa. Our data suggest that these microbial communities have been shaped by geochemical drivers, including water content, salinity, and the supply of organic matter. Using Illumina 16S rRNA gene sequencing, this study provides new insights into the diversity of bacteria and archaea in semi-arid, saline, and low-carbon environments. Many of the observed taxa are halophilic and adapted to water-limiting conditions. The analysis reveals a high relative abundance of halophilic archaea (primarily Halobacteria), and the bacterial diversity is marked by an abundance of Gemmatimonadetes and spore-forming Firmicutes. In the deeper, anoxic layers, candidate division MSBL1, and acetogenic bacteria (Acetothermia) are abundant. Together, the taxonomic information and geochemical data suggest that acetogenesis could be a prevalent form of metabolism in the deep layers of a saline pan.
Salt pans are highly dynamic environments that are difficult to study by in situ methods because of their harsh climatic conditions and large spatial areas. Remote sensing can help to elucidate their environmental dynamics and provide important constraints regarding their sedimentological, mineralogical, and hydrological evolution. This study utilizes spaceborne multitemporal multispectral optical data combined with spectral endmembers to document spatial distribution of surface crust types over time on the Omongwa pan located in the Namibian Kalahari. For this purpose, 49 surface samples were collected for spectral and mineralogical characterization during three field campaigns (2014–2016) reflecting different seasons and surface conditions of the salt pan. An approach was developed to allow the spatiotemporal analysis of the salt pan crust dynamics in a dense time-series consisting of 77 Landsat 8 cloud-free scenes between 2014 and 2017, covering at least three major wet–dry cycles. The established spectral analysis technique Sequential Maximum Angle Convex Cone (SMACC) extraction method was used to derive image endmembers from the Landsat time-series stack. Evaluation of the extracted endmember set revealed that the multispectral data allowed the differentiation of four endmembers associated with mineralogical mixtures of the crust’s composition in dry conditions and three endmembers associated with flooded or muddy pan conditions. The dry crust endmember spectra have been identified in relation to visible, near infrared, and short-wave infrared (VNIR–SWIR) spectroscopy and X-ray diffraction (XRD) analyses of the collected surface samples. According these results, the spectral endmembers are interpreted as efflorescent halite crust, mixed halite–gypsum crust, mixed calcite quartz sepiolite crust, and gypsum crust. For each Landsat scene the spatial distribution of these crust types was mapped with the Spectral Angle Mapper (SAM) method and significant spatiotemporal dynamics of the major surface crust types were observed. Further, the surface crust dynamics were analyzed in comparison with the pan’s moisture regime and other climatic parameters. The results show that the crust dynamics are mainly driven by flooding events in the wet season, but are also influenced by temperature and aeolian activity in the dry season. The approach utilized in this study combines the advantages of multitemporal satellite data for temporal event characterization with advantages from hyperspectral methods for the image and ground data analyses that allow improved mineralogical differentiation and characterization.
This study combines spaceborne multitemporal and hyperspectral data to analyze the spatial distribution of surface evaporite minerals and changes in a semi-arid depositional environment associated with episodic flooding events, the Omongwa salt pan (Kalahari, Namibia). The dynamic of the surface crust is evaluated by a change-detection approach using the Iterative-reweighted Multivariate Alteration Detection (IR-MAD) based on the Landsat archive imagery from 1984 to 2015. The results show that the salt pan is a highly dynamic and heterogeneous landform. A change gradient is observed from very stable pan border to a highly dynamic central pan. On the basis of hyperspectral EO-1 Hyperion images, the current distribution of surface evaporite minerals is characterized using Spectral Mixture Analysis (SMA). Assessment of field and image endmembers revealed that the pan surface can be categorized into three major crust types based on diagnostic absorption features and mineralogical ground truth data. The mineralogical crust types are related to different zones of surface change as well as pan morphology that influences brine flow during the pan inundation and desiccation cycles. These combined information are used to spatially map depositional environments where the more dynamic halite crust concentrates in lower areas although stable gypsum and calcite/sepiolite crusts appear in higher elevated areas.
Organic matter characteristics in yedoma and thermokarst deposits on Baldwin Peninsula, west Alaska
(2018)
As Arctic warming continues and permafrost thaws, more soil and sedimentary organic matter (OM) will be decomposed in northern high latitudes. Still, uncertainties remain in the quality of the OM and the size of the organic carbon (OC) pools stored in different deposit types of permafrost landscapes. This study presents OM data from deep permafrost and lake deposits on the Baldwin Peninsula which is located in the southern portion of the continuous permafrost zone in west Alaska. Sediment samples from yedoma and drained thermokarst lake basin (DTLB) deposits as well as thermokarst lake sediments were analyzed for cryostratigraphical and biogeochemical parameters and their lipid biomarker composition to identify the below-ground OC pool size and OM quality of ice-rich permafrost on the Baldwin Peninsula. We provide the first detailed characterization of yedoma deposits on Baldwin Peninsula. We show that three-quarters of soil OC in the frozen deposits of the study region (total of 68 Mt) is stored in DTLB deposits (52 Mt) and one-quarter in the frozen yedoma deposits (16 Mt). The lake sediments contain a relatively small OC pool (4 Mt), but have the highest volumetric OC content (93 kgm(-3)) compared to the DTLB (35 kgm(-3)) and yedoma deposits (8 kgm(-3)), largely due to differences in the ground ice content. The biomarker analysis indicates that the OM in both yedoma and DTLB deposits is mainly of terrestrial origin. Nevertheless, the relatively high carbon preference index of plant leaf waxes in combination with a lack of a degradation trend with depth in the yedoma deposits indi-cates that OM stored in yedoma is less degraded than that stored in DTLB deposits. This suggests that OM in yedoma has a higher potential for decomposition upon thaw, despite the relatively small size of this pool. These findings show that the use of lipid biomarker analysis is valuable in the assessment of the potential future greenhouse gas emissions from thawing permafrost, especially because this area, close to the discontinuous permafrost boundary, is projected to thaw substantially within the 21st century.
Solar wind observations show that geomagnetic storms are mainly driven by interplanetary coronal mass ejections (ICMEs) and corotating or stream interaction regions (C/SIRs). We present a binary classifier that assigns one of these drivers to 7,546 storms between 1930 and 2015 using ground‐based geomagnetic field observations only. The input data consists of the long‐term stable Hourly Magnetospheric Currents index alongside the corresponding midlatitude geomagnetic observatory time series. This data set provides comprehensive information on the global storm time magnetic disturbance field, particularly its spatial variability, over eight solar cycles. For the first time, we use this information statistically with regard to an automated storm driver identification. Our supervised classification model significantly outperforms unskilled baseline models (78% accuracy with 26[19]% misidentified interplanetary coronal mass ejections [corotating or stream interaction regions]) and delivers plausible driver occurrences with regard to storm intensity and solar cycle phase. Our results can readily be used to advance related studies fundamental to space weather research, for example, studies connecting galactic cosmic ray modulation and geomagnetic disturbances. They are fully reproducible by means of the underlying open‐source software (Pick, 2019, http://doi.org/10.5880/GFZ.2.3.2019.003)
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.
Garnet of eclogite (formerly termed garnet clinopyroxenite) hosted in lenses of orogenic garnet peridotite from the Granulitgebirge, NW Bohemian Massif, contains unique inclusions of granitic melt, now either glassy or crystallized. Analysed glasses and re‐homogenized inclusions are hydrous, peraluminous, and enriched in highly incompatible elements characteristic of the continental crust such as Cs, Li, B, Pb, Rb, Th, and U. The original melt thus represents a pristine, chemically evolved metasomatic agent, which infiltrated the mantle via deep continental subduction during the Variscan orogeny. The bulk chemical composition of the studied eclogites is similar to that of Fe‐rich basalt and the enrichment in LILE and U suggest a subduction‐related component. All these geochemical features confirm metasomatism. In comparison with many other garnet+clinopyroxene‐bearing lenses in peridotites of the Bohemian Massif, the studied samples from Rubinberg and Klatschmühle are more akin to eclogite than pyroxenites, as reflected in high jadeite content in clinopyroxene, relatively low Mg, Cr, and Ni but relatively high Ti. However, trace elements of both bulk rock and individual mineral phases show also important differences making these samples rather unique. Metasomatism involving a melt requiring a trace element pattern very similar to the composition reported here has been suggested for the source region of rocks of the so‐called durbachite suite, that is, ultrapotassic melanosyenites, which are found throughout the high‐grade Variscan basement. Moreover, the Th, U, Pb, Nb, Ta, and Ti patterns of these newly studied melt inclusions (MI) strongly resemble those observed for peridotite and its enclosed pyroxenite from the T‐7 borehole (Staré, České Středhoři Mountains) in N Bohemia. This suggests that a similar kind of crustal‐derived melt also occurred here. This study of granitic MI in eclogites from peridotites has provided the first direct characterization of a preserved metasomatic melt, possibly responsible for the metasomatism of several parts of the mantle in the Variscides.
Sedimentary ancient DNA has been proposed as a key methodology for reconstructing biodiversity over time. Yet, despite the concentration of Earth’s biodiversity in the tropics, this method has rarely been applied in this region. Moreover, the taphonomy of sedimentary DNA, especially in tropical environments, is poorly understood. This study elucidates challenges and opportunities of sedimentary ancient DNA approaches for reconstructing tropical biodiversity. We present shotgun-sequenced metagenomic profiles and DNA degradation patterns from multiple sediment cores from Mubwindi Swamp, located in Bwindi Impenetrable Forest (Uganda), one of the most diverse forests in Africa. We describe the taxonomic composition of the sediments covering the past 2200 years and compare the sedimentary DNA data with a comprehensive set of environmental and sedimentological parameters to unravel the conditions of DNA degradation. Consistent with the preservation of authentic ancient DNA in tropical swamp sediments, DNA concentration and mean fragment length declined exponentially with age and depth, while terminal deamination increased with age. DNA preservation patterns cannot be explained by any environmental parameter alone, but age seems to be the primary driver of DNA degradation in the swamp. Besides degradation, the presence of living microbial communities in the sediment also affects DNA quantity. Critically, 92.3% of our metagenomic data of a total 81.8 million unique, merged reads cannot be taxonomically identified due to the absence of genomic references in public databases. Of the remaining 7.7%, most of the data (93.0%) derive from Bacteria and Archaea, whereas only 0–5.8% are from Metazoa and 0–6.9% from Viridiplantae, in part due to unbalanced taxa representation in the reference data. The plant DNA record at ordinal level agrees well with local pollen data but resolves less diversity. Our animal DNA record reveals the presence of 41 native taxa (16 orders) including Afrotheria, Carnivora, and Ruminantia at Bwindi during the past 2200 years. Overall, we observe no decline in taxonomic richness with increasing age suggesting that several-thousand-year-old information on past biodiversity can be retrieved from tropical sediments. However, comprehensive genomic surveys of tropical biota need prioritization for sedimentary DNA to be a viable methodology for future tropical biodiversity studies.
RainNet v1.0
(2020)
In this study, we present RainNet, a deep convolutional neural network for radar-based precipitation nowcasting. Its design was inspired by the U-Net and SegNet families of deep learning models, which were originally designed for binary segmentation tasks. RainNet was trained to predict continuous precipitation intensities at a lead time of 5min, using several years of quality-controlled weather radar composites provided by the German Weather Service (DWD). That data set covers Germany with a spatial domain of 900km × 900km and has a resolution of 1km in space and 5min in time. Independent verification experiments were carried out on 11 summer precipitation events from 2016 to 2017. In order to achieve a lead time of 1h, a recursive approach was implemented by using RainNet predictions at 5min lead times as model inputs for longer lead times. In the verification experiments, trivial Eulerian persistence and a conventional model based on optical flow served as benchmarks. The latter is available in the rainymotion library and had previously been shown to outperform DWD's operational nowcasting model for the same set of verification events.
RainNet significantly outperforms the benchmark models at all lead times up to 60min for the routine verification metrics mean absolute error (MAE) and the critical success index (CSI) at intensity thresholds of 0.125, 1, and 5mm h⁻¹. However, rainymotion turned out to be superior in predicting the exceedance of higher intensity thresholds (here 10 and 15mm h⁻¹). The limited ability of RainNet to predict heavy rainfall intensities is an undesirable property which we attribute to a high level of spatial smoothing introduced by the model. At a lead time of 5min, an analysis of power spectral density confirmed a significant loss of spectral power at length scales of 16km and below. Obviously, RainNet had learned an optimal level of smoothing to produce a nowcast at 5min lead time. In that sense, the loss of spectral power at small scales is informative, too, as it reflects the limits of predictability as a function of spatial scale. Beyond the lead time of 5min, however, the increasing level of smoothing is a mere artifact – an analogue to numerical diffusion – that is not a property of RainNet itself but of its recursive application. In the context of early warning, the smoothing is particularly unfavorable since pronounced features of intense precipitation tend to get lost over longer lead times. Hence, we propose several options to address this issue in prospective research, including an adjustment of the loss function for model training, model training for longer lead times, and the prediction of threshold exceedance in terms of a binary segmentation task. Furthermore, we suggest additional input data that could help to better identify situations with imminent precipitation dynamics. The model code, pretrained weights, and training data are provided in open repositories as an input for such future studies.
Hydrometric networks play a vital role in providing information for decision-making in water resource management. They should be set up optimally to provide as much information as possible that is as accurate as possible and, at the same time, be cost-effective. Although the design of hydrometric networks is a well-identified problem in hydrometeorology and has received considerable attention, there is still scope for further advancement. In this study, we use complex network analysis, defined as a collection of nodes interconnected by links, to propose a new measure that identifies critical nodes of station networks. The approach can support the design and redesign of hydrometric station networks. The science of complex networks is a relatively young field and has gained significant momentum over the last few years in different areas such as brain networks, social networks, technological networks, or climate networks. The identification of influential nodes in complex networks is an important field of research. We propose a new node-ranking measure – the weighted degree–betweenness (WDB) measure – to evaluate the importance of nodes in a network. It is compared to previously proposed measures used on synthetic sample networks and then applied to a real-world rain gauge network comprising 1229 stations across Germany to demonstrate its applicability. The proposed measure is evaluated using the decline rate of the network efficiency and the kriging error. The results suggest that WDB effectively quantifies the importance of rain gauges, although the benefits of the method need to be investigated in more detail.
Many institutions struggle to tap into the potential of their large archives of radar reflectivity: these data are often affected by miscalibration, yet the bias is typically unknown and temporally volatile. Still, relative calibration techniques can be used to correct the measurements a posteriori. For that purpose, the usage of spaceborne reflectivity observations from the Tropical Rainfall Measuring Mission (TRMM) and Global Precipitation Measurement (GPM) platforms has become increasingly popular: the calibration bias of a ground radar (GR) is estimated from its average reflectivity difference to the spaceborne radar (SR). Recently, Crisologo et al. (2018) introduced a formal procedure to enhance the reliability of such estimates: each match between SR and GR observations is assigned a quality index, and the calibration bias is inferred as a quality-weighted average of the differences between SR and GR. The relevance of quality was exemplified for the Subic S-band radar in the Philippines, which is greatly affected by partial beam blockage.
The present study extends the concept of quality-weighted averaging by accounting for path-integrated attenuation (PIA) in addition to beam blockage. This extension becomes vital for radars that operate at the C or X band. Correspondingly, the study setup includes a C-band radar that substantially overlaps with the S-band radar. Based on the extended quality-weighting approach, we retrieve, for each of the two ground radars, a time series of calibration bias estimates from suitable SR overpasses. As a result of applying these estimates to correct the ground radar observations, the consistency between the ground radars in the region of overlap increased substantially. Furthermore, we investigated if the bias estimates can be interpolated in time, so that ground radar observations can be corrected even in the absence of prompt SR overpasses. We found that a moving average approach was most suitable for that purpose, although limited by the absence of explicit records of radar maintenance operations.
Introducing PebbleCounts
(2019)
Grain-size distributions are a key geomorphic metric of gravel-bed rivers. Traditional measurement methods include manual counting or photo sieving, but these are achievable only at the 1–10 ㎡ scale. With the advent of drones and increasingly high-resolution cameras, we can now generate orthoimagery over hectares at millimeter to centimeter resolution. These scales, along with the complexity of high-mountain rivers, necessitate different approaches for photo sieving. As opposed to other image segmentation methods that use a watershed approach, our open-source algorithm, PebbleCounts, relies on k-means clustering in the spatial and spectral domain and rapid manual selection of well-delineated grains. This improves grain-size estimates for complex riverbed imagery, without post-processing. We also develop a fully automated method, PebbleCountsAuto, that relies on edge detection and filtering suspect grains, without the k-means clustering or manual selection steps. The algorithms are tested in controlled indoor conditions on three arrays of pebbles and then applied to 12 × 1 ㎡ orthomosaic clips of high-energy mountain rivers collected with a camera-on-mast setup (akin to a low-flying drone). A 20-pixel b-axis length lower truncation is necessary for attaining accurate grain-size distributions. For the k-means PebbleCounts approach, average percentile bias and precision are 0.03 and 0.09 ψ, respectively, for ∼1.16 mm pixel⁻¹ images, and 0.07 and 0.05 ψ for one 0.32 mm pixel⁻¹ image. The automatic approach has higher bias and precision of 0.13 and 0.15 ψ, respectively, for ∼1.16 mm pixel⁻¹ images, but similar values of −0.06 and 0.05 ψ for one 0.32 mm pixel⁻¹ image. For the automatic approach, only at best 70 % of the grains are correct identifications, and typically around 50 %. PebbleCounts operates most effectively at the 1 ㎡ patch scale, where it can be applied in ∼5–10 min on many patches to acquire accurate grain-size data over 10–100 ㎡ areas. These data can be used to validate PebbleCountsAuto, which may be applied at the scale of entire survey sites (102–104 ㎡ ). We synthesize results and recommend best practices for image collection, orthomosaic generation, and grain-size measurement using both algorithms.
The interactions between atmosphere and steep topography in the eastern south–central Andes result in complex relations with inhomogenous rainfall distributions. The atmospheric conditions leading to deep convection and extreme rainfall and their spatial patterns—both at the valley and mountain-belt scales—are not well understood. In this study, we aim to identify the dominant atmospheric conditions and their spatial variability by analyzing the convective available potential energy (CAPE) and dew-point temperature (Td). We explain the crucial effect of temperature on extreme rainfall generation along the steep climatic and topographic gradients in the NW Argentine Andes stretching from the low-elevation eastern foreland to the high-elevation central Andean Plateau in the west. Our analysis relies on version 2.0 of the ECMWF’s (European Centre for Medium-RangeWeather Forecasts) Re-Analysis (ERA-interim) data and TRMM (Tropical Rainfall Measuring Mission) data. We make the following key observations: First, we observe distinctive gradients along and across strike of the Andes in dew-point temperature and CAPE that both control rainfall distributions. Second, we identify a nonlinear correlation between rainfall and a combination of dew-point temperature and CAPE through a multivariable regression analysis. The correlation changes in space along the climatic and topographic gradients and helps to explain controlling factors for extreme-rainfall generation. Third, we observe more contribution (or higher importance) of Td in the tropical low-elevation foreland and intermediate-elevation areas as compared to the high-elevation central Andean Plateau for 90th percentile rainfall. In contrast, we observe a higher contribution of CAPE in the intermediate-elevation area between low and high elevation, especially in the transition zone between the tropical and subtropical areas for the 90th percentile rainfall. Fourth, we find that the parameters of the multivariable regression using CAPE and Td can explain rainfall with higher statistical significance for the 90th percentile compared to lower rainfall percentiles. Based on our results, the spatial pattern of rainfall-extreme events during the past ∼16 years can be described by a combination of dew-point temperature and CAPE in the south–central Andes.
The Arctic-Boreal regions experience strong changes of air temperature and precipitation regimes, which affect the thermal state of the permafrost. This results in widespread permafrost-thaw disturbances, some unfolding slowly and over long periods, others occurring rapidly and abruptly. Despite optical remote sensing offering a variety of techniques to assess and monitor landscape changes, a persistent cloud cover decreases the amount of usable images considerably. However, combining data from multiple platforms promises to increase the number of images drastically. We therefore assess the comparability of Landsat-8 and Sentinel-2 imagery and the possibility to use both Landsat and Sentinel-2 images together in time series analyses, achieving a temporally-dense data coverage in Arctic-Boreal regions. We determined overlapping same-day acquisitions of Landsat-8 and Sentinel-2 images for three representative study sites in Eastern Siberia. We then compared the Landsat-8 and Sentinel-2 pixel-pairs, downscaled to 60 m, of corresponding bands and derived the ordinary least squares regression for every band combination. The acquired coefficients were used for spectral bandpass adjustment between the two sensors. The spectral band comparisons showed an overall good fit between Landsat-8 and Sentinel-2 images already. The ordinary least squares regression analyses underline the generally good spectral fit with intercept values between 0.0031 and 0.056 and slope values between 0.531 and 0.877. A spectral comparison after spectral bandpass adjustment of Sentinel-2 values to Landsat-8 shows a nearly perfect alignment between the same-day images. The spectral band adjustment succeeds in adjusting Sentinel-2 spectral values to Landsat-8 very well in Eastern Siberian Arctic-Boreal landscapes. After spectral adjustment, Landsat and Sentinel-2 data can be used to create temporally-dense time series and be applied to assess permafrost landscape changes in Eastern Siberia. Remaining differences between the sensors can be attributed to several factors including heterogeneous terrain, poor cloud and cloud shadow masking, and mixed pixels.
Hydrometeorological hazards caused losses of approximately 110 billion U.S. Dollars in 2016 worldwide. Current damage estimations do not consider the uncertainties in a comprehensive way, and they are not consistent between spatial scales. Aggregated land use data are used at larger spatial scales, although detailed exposure data at the object level, such as openstreetmap.org, is becoming increasingly available across the globe.We present a probabilistic approach for object-based damage estimation which represents uncertainties and is fully scalable in space. The approach is applied and validated to company damage from the flood of 2013 in Germany. Damage estimates are more accurate compared to damage models using land use data, and the estimation works reliably at all spatial scales. Therefore, it can as well be used for pre-event analysis and risk assessments. This method takes hydrometeorological damage estimation and risk assessments to the next level, making damage estimates and their uncertainties fully scalable in space, from object to country level, and enabling the exploitation of new exposure data.
Abstract. The Sea of Marmara, in northwestern Turkey, is a transition zone where the dextral North Anatolian Fault zone (NAFZ) propagates westward from the Anatolian Plate to the Aegean Sea Plate. The area is of interest in the context of seismic hazard of Istanbul, a metropolitan area with about 15 million inhabitants. Geophysical observations indicate that the crust is heterogeneous beneath the Marmara basin, but a detailed characterization of the crustal heterogeneities is still missing. To assess if and how crustal heterogeneities are related to the NAFZ segmentation below the Sea of Marmara, we develop new crustal-scale 3-D density models which integrate geological and seismological data and that are additionally constrained by 3-D gravity modeling. For the latter, we use two different gravity datasets including global satellite data and local marine gravity observation. Considering the two different datasets and the general non-uniqueness in potential field modeling, we suggest three possible “end-member” solutions that are all consistent with the observed gravity field and illustrate the spectrum of possible solutions. These models indicate that the observed gravitational anomalies originate from significant density heterogeneities within the crust. Two layers of sediments, one syn-kinematic and one pre-kinematic with respect to the Sea of Marmara formation are underlain by a heterogeneous crystalline crust. A felsic upper crystalline crust (average density of 2720 kgm⁻³) and an intermediate to mafic lower crystalline crust (average density of 2890 kgm⁻³) appear to be cross-cut by two large, dome-shaped mafic highdensity bodies (density of 2890 to 3150 kgm⁻³) of considerable thickness above a rather uniform lithospheric mantle (3300 kgm⁻³). The spatial correlation between two major bends of the main Marmara fault and the location of the highdensity bodies suggests that the distribution of lithological heterogeneities within the crust controls the rheological behavior along the NAFZ and, consequently, maybe influences fault segmentation and thus the seismic hazard assessment in the region.
Sea surface temperature (SST) patterns can – as surface climate forcing – affect weather and climate at large distances. One example is El Niño-Southern Oscillation (ENSO) that causes climate anomalies around the globe via teleconnections. Although several studies identified and characterized these teleconnections, our understanding of climate processes remains incomplete, since interactions and feedbacks are typically exhibited at unique or multiple temporal and spatial scales. This study characterizes the interactions between the cells of a global SST data set at different temporal and spatial scales using climate networks. These networks are constructed using wavelet multi-scale correlation that investigate the correlation between the SST time series at a range of scales allowing instantaneously deeper insights into the correlation patterns compared to traditional methods like empirical orthogonal functions or classical correlation analysis. This allows us to identify and visualise regions of – at a certain timescale – similarly evolving SSTs and distinguish them from those with long-range teleconnections to other ocean regions. Our findings re-confirm accepted knowledge about known highly linked SST patterns like ENSO and the Pacific Decadal Oscillation, but also suggest new insights into the characteristics and origins of long-range teleconnections like the connection between ENSO and Indian Ocean Dipole.
Determining the optimal grid resolution for topographic analysis on an airborne lidar dataset
(2019)
Digital elevation models (DEMs) are a gridded representation of the surface of the Earth and typically contain uncertainties due to data collection and processing. Slope and aspect estimates on a DEM contain errors and uncertainties inherited from the representation of a continuous surface as a grid (referred to as truncation error; TE) and from any DEM uncertainty. We analyze in detail the impacts of TE and propagated elevation uncertainty (PEU) on slope and aspect.
Using synthetic data as a control, we define functions to quantify both TE and PEU for arbitrary grids. We then develop a quality metric which captures the combined impact of both TE and PEU on the calculation of topographic metrics. Our quality metric allows us to examine the spatial patterns of error and uncertainty in topographic metrics and to compare calculations on DEMs of different sizes and accuracies.
Using lidar data with point density of ∼10 pts m−2 covering Santa Cruz Island in southern California, we are able to generate DEMs and uncertainty estimates at several grid resolutions. Slope (aspect) errors on the 1 m dataset are on average 0.3∘ (0.9∘) from TE and 5.5∘ (14.5∘) from PEU. We calculate an optimal DEM resolution for our SCI lidar dataset of 4 m that minimizes the error bounds on topographic metric calculations due to the combined influence of TE and PEU for both slope and aspect calculations over the entire SCI. Average slope (aspect) errors from the 4 m DEM are 0.25∘ (0.75∘) from TE and 5∘ (12.5∘) from PEU. While the smallest grid resolution possible from the high-density SCI lidar is not necessarily optimal for calculating topographic metrics, high point-density data are essential for measuring DEM uncertainty across a range of resolutions.
The ICDP "PaleoVan" drilling campaign at Lake Van, Turkey, provided a long (> 100 m) record of lacustrine subsurface sedimentary microbial cell abundance. After the ICDP campaign at Potrok Aike, Argentina, this is only the second time deep lacustrine cell counts have been documented. Two sites were cored and revealed a strikingly similar cell distribution despite differences in organic matter content and microbial activity. Although shifted towards higher values, cell counts from Lake Potrok Aike, Argentina, reveal very similar distribution patterns with depth. The lacustrine cell count data are significantly different from published marine records; the most probable cause is differences in sedimentary organic matter composition with marine sediments containing a higher fraction of labile organic matter. Previous studies showed that microbial activity and abundance increase centimetres to metres around geologic interfaces. The finely laminated Lake Van sediment allowed studying this phenomenon on the microscale. We sampled at the scale of individual laminae, and in some depth intervals, we found large differences in microbial abundance between the different laminae. This small-scale heterogeneity is normally overlooked due to much larger sampling intervals that integrate over several centimetres. However, not all laminated intervals exhibit such large differences in microbial abundance, and some non-laminated horizons show large variability on the millimetre scale as well. The reasons for such contrasting observations remain elusive, but indicate that heterogeneity of microbial abundance in subsurface sediments has not been taken into account sufficiently. These findings have implications not just for microbiological studies but for geochemistry as well, as the large differences in microbial abundance clearly show that there are distinct microhabitats that deviate considerably from the surrounding layers.
The semiarid northeast of Brazil is one of the most densely populated dryland regions in the world and recurrently affected by severe droughts. Thus, reliable seasonal forecasts of streamflow and reservoir storage are of high value for water managers. Such forecasts can be generated by applying either hydrological models representing underlying processes or statistical relationships exploiting correlations among meteorological and hydrological variables. This work evaluates and compares the performances of seasonal reservoir storage forecasts derived by a process-based hydrological model and a statistical approach.
Driven by observations, both models achieve similar simulation accuracies. In a hindcast experiment, however, the accuracy of estimating regional reservoir storages was considerably lower using the process-based hydrological model, whereas the resolution and reliability of drought event predictions were similar by both approaches. Further investigations regarding the deficiencies of the process-based model revealed a significant influence of antecedent wetness conditions and a higher sensitivity of model prediction performance to rainfall forecast quality.
Within the scope of this study, the statistical model proved to be the more straightforward approach for predictions of reservoir level and drought events at regionally and monthly aggregated scales. However, for forecasts at finer scales of space and time or for the investigation of underlying processes, the costly initialisation and application of a process-based model can be worthwhile. Furthermore, the application of innovative data products, such as remote sensing data, and operational model correction methods, like data assimilation, may allow for an enhanced exploitation of the advanced capabilities of process-based hydrological models.
We measure valence-to-core x-ray emission spectra of compressed crystalline GeO₂ up to 56 GPa and of amorphous GeO₂ up to 100 GPa. In a novel approach, we extract the Ge coordination number and mean Ge-O distances from the emission energy and the intensity of the Kβ'' emission line. The spectra of high-pressure polymorphs are calculated using the Bethe-Salpeter equation. Trends observed in the experimental and calculated spectra are found to match only when utilizing an octahedral model. The results reveal persistent octahedral Ge coordination with increasing distortion, similar to the compaction mechanism in the sequence of octahedrally coordinated crystalline GeO₂ high-pressure polymorphs.
Extreme weather events are likely to occur more often under climate change and the resulting effects on ecosystems could lead to a further acceleration of climate change. But not all extreme weather events lead to extreme ecosystem response. Here, we focus on hazardous ecosystem behaviour and identify coinciding weather conditions. We use a simple probabilistic risk assessment based on time series of ecosystem behaviour and climate conditions. Given the risk assessment terminology, vulnerability and risk for the previously defined hazard are estimated on the basis of observed hazardous ecosystem behaviour.
We apply this approach to extreme responses of terrestrial ecosystems to drought, defining the hazard as a negative net biome productivity over a 12-month period. We show an application for two selected sites using data for 1981-2010 and then apply the method to the pan-European scale for the same period, based on numerical modelling results (LPJmL for ecosystem behaviour; ERA-Interim data for climate).
Our site-specific results demonstrate the applicability of the proposed method, using the SPEI to describe the climate condition. The site in Spain provides an example of vulnerability to drought because the expected value of the SPEI is 0.4 lower for hazardous than for non-hazardous ecosystem behaviour. In northern Germany, on the contrary, the site is not vulnerable to drought because the SPEI expectation values imply wetter conditions in the hazard case than in the non-hazard case.
At the pan-European scale, ecosystem vulnerability to drought is calculated in the Mediterranean and temperate region, whereas Scandinavian ecosystems are vulnerable under conditions without water shortages. These first model- based applications indicate the conceptual advantages of the proposed method by focusing on the identification of critical weather conditions for which we observe hazardous ecosystem behaviour in the analysed data set. Application of the method to empirical time series and to future climate would be important next steps to test the approach.
In the arctic and high mountains it is common to measure vertical changes of ice sheets and glaciers via digital elevation model (DEM) differencing. This requires the signal of change to outweigh the noise associated with the datasets. Excluding large landslides, on the ice-free earth the land-level change is smaller in vertical magnitude and thus requires more accurate DEMs for differencing and identification of change. Previously, this has required meter to submeter data at small spatial scales. Following careful corrections, we are able to measure land-level changes in gravel-bed channels and steep hillslopes in the south-central Andes using the SRTM-C (collected in 2000) and the TanDEM-X (collected from 2010 to 2015) near-global 12–30m DEMs. Long-standing errors in the SRTM-C are corrected using the TanDEM-X as a control surface and applying cosine-fit co-registration to remove ∼ 1∕10 pixel (∼ 3m) shifts, fast Fourier transform (FFT) and filtering to remove SRTM-C short- and long-wavelength stripes, and blocked shifting to remove remaining complex biases. The datasets are then differenced and outlier pixels are identified as a potential signal for the case of gravel-bed channels and hillslopes. We are able to identify signals of incision and aggradation (with magnitudes down to ∼ 3m in the best case) in two > 100km river reaches, with increased geomorphic activity downstream of knickpoints. Anthropogenic gravel excavation and piling is prominently measured, with magnitudes exceeding ±5m (up to > 10m for large piles). These values correspond to conservative average rates of 0.2 to > 0.5myr−1 for vertical changes in gravel-bed rivers. For hillslopes, since we require stricter cutoffs for noise, we are only able to identify one major landslide in the study area with a deposit volume of 16±0.15×106m3. Additional signals of change can be garnered from TanDEM-X auxiliary layers; however, these are more difficult to quantify. The methods presented can be extended to any region of the world with SRTM-C and TanDEM-X coverage where vertical land-level changes are of interest, with the caveat that remaining vertical uncertainties in primarily the SRTM-C limit detection in steep and complex topography.
We explore the potential of spaceborne radar (SR) observations from the Ku-band precipitation radars onboard the Tropical Rainfall Measuring Mission (TRMM) and Global Precipitation Measurement (GPM) satellites as a reference to quantify the ground radar (GR) reflectivity bias. To this end, the 3-D volume-matching algorithm proposed by Schwaller and Morris (2011) is implemented and applied to 5 years (2012–2016) of observations. We further extend the procedure by a framework to take into account the data quality of each ground radar bin. Through these methods, we are able to assign a quality index to each matching SR–GR volume, and thus compute the GR calibration bias as a quality-weighted average of reflectivity differences in any sample of matching GR–SR volumes. We exemplify the idea of quality-weighted averaging by using the beam blockage fraction as the basis of a quality index. As a result, we can increase the consistency of SR and GR observations, and thus the precision of calibration bias estimates. The remaining scatter between GR and SR reflectivity as well as the variability of bias estimates between overpass events indicate, however, that other error sources are not yet fully addressed. Still, our study provides a framework to introduce any other quality variables that are considered relevant in a specific context. The code that implements our analysis is based on the wradlib open-source software library, and is, together with the data, publicly available to monitor radar calibration or to scrutinize long series of archived radar data back to December 1997, when TRMM became operational.
Mapping Damage-Affected Areas after Natural Hazard Events Using Sentinel-1 Coherence Time Series
(2018)
The emergence of the Sentinel-1A and 1B satellites now offers freely available and widely accessible Synthetic Aperture Radar (SAR) data. Near-global coverage and rapid repeat time (6–12 days) gives Sentinel-1 data the potential to be widely used for monitoring the Earth’s surface. Subtle land-cover and land surface changes can affect the phase and amplitude of the C-band SAR signal, and thus the coherence between two images collected before and after such changes. Analysis of SAR coherence therefore serves as a rapidly deployable and powerful tool to track both seasonal changes and rapid surface disturbances following natural disasters. An advantage of using Sentinel-1 C-band radar data is the ability to easily construct time series of coherence for a region of interest at low cost. In this paper, we propose a new method for Potentially Affected Area (PAA) detection following a natural hazard event. Based on the coherence time series, the proposed method (1) determines the natural variability of coherence within each pixel in the region of interest, accounting for factors such as seasonality and the inherent noise of variable surfaces; and (2) compares pixel-by-pixel syn-event coherence to temporal coherence distributions to determine where statistically significant coherence loss has occurred. The user can determine to what degree the syn-event coherence value (e.g., 1st, 5th percentile of pre-event distribution) constitutes a PAA, and integrate pertinent regional data, such as population density, to rank and prioritise PAAs. We apply the method to two case studies, Sarpol-e, Iran following the 2017 Iran-Iraq earthquake, and a landslide-prone region of NW Argentina, to demonstrate how rapid identification and interpretation of potentially affected areas can be performed shortly following a natural hazard event.
Observed recent and expected future increases in frequency and intensity of climatic extremes in central Europe may pose critical challenges for domestic tree species. Continuous dendrometer recordings provide a valuable source of information on tree stem radius variations, offering the possibility to study a tree's response to environmental influences at a high temporal resolution. In this study, we analyze stem radius variations (SRV) of three domestic tree species (beech, oak, and pine) from 2012 to 2014. We use the novel statistical approach of event coincidence analysis (ECA) to investigate the simultaneous occurrence of extreme daily weather conditions and extreme SRVs, where extremes are defined with respect to the common values at a given phase of the annual growth period. Besides defining extreme events based on individual meteorological variables, we additionally introduce conditional and joint ECA as new multivariate extensions of the original methodology and apply them for testing 105 different combinations of variables regarding their impact on SRV extremes. Our results reveal a strong susceptibility of all three species to the extremes of several meteorological variables. Yet, the inter-species differences regarding their response to the meteorological extremes are comparatively low. The obtained results provide a thorough extension of previous correlation-based studies by emphasizing on the timings of climatic extremes only. We suggest that the employed methodological approach should be further promoted in forest research regarding the investigation of tree responses to changing environmental conditions.