Refine
Has Fulltext
- no (151) (remove)
Year of publication
- 2020 (151) (remove)
Document Type
- Article (142)
- Doctoral Thesis (4)
- Monograph/Edited Volume (2)
- Other (2)
- Review (1)
Is part of the Bibliography
- yes (151)
Keywords
- climate change (5)
- Himalaya (3)
- model (3)
- models (3)
- precipitation (3)
- Arctic (2)
- CO2 (2)
- China (2)
- Chinese loess (2)
- Europe (2)
- LGM (2)
- Siberia (2)
- biomarker (2)
- body waves (2)
- carbonate (2)
- climate-change (2)
- diversity (2)
- dynamics (2)
- early warning (2)
- earthquake source observations (2)
- erosion (2)
- extension (2)
- glacier (2)
- hydrological drought (2)
- induced seismicity (2)
- isotopes (2)
- modelling (2)
- patterns (2)
- permafrost (2)
- permeability (2)
- plateau (2)
- plume (2)
- reactive transport (2)
- remote sensing (2)
- rivers (2)
- salt diffusion (2)
- site effects (2)
- snow (2)
- subduction (2)
- subduction zone (2)
- tectonics (2)
- thermokarst (2)
- time series analysis (2)
- 239+240 Plutonium (1)
- 3D modeling (1)
- Affordability (1)
- Aftershocks (1)
- Algorithm (1)
- Andean Plateau (1)
- Andes (1)
- Antarctic ice (1)
- Apatite (U-TH)/HE (1)
- Applied geophysics (1)
- Ar-40/Ar-39 age (1)
- Arctic lakes (1)
- Band (1)
- Basanite (1)
- Bayesian networks (1)
- Bedded chert (1)
- Bentheim sandstone (1)
- Biomedical engineering (1)
- Bohemian Massif (1)
- Brain injuries (1)
- Brazilian Earth System Model (1)
- Brittle precursors (1)
- Bulk compressibility (1)
- C-Q relationship (1)
- CCSEM (1)
- CMIP5 (1)
- CMIP5 models (1)
- COH fluid (1)
- Cantera (1)
- Carnian Pluvial Event (1)
- Central andes (1)
- Central-Asia (1)
- Clay Minerals (1)
- Compound flood event (1)
- Computational models (1)
- Computational neuroscience (1)
- Contact model (1)
- Copernicus (1)
- Coulomb failure stress (1)
- Crustal density (1)
- DNA preservation (1)
- Data acquisition (1)
- Data processing (1)
- Deep ocean redox condition (1)
- ECHSE (1)
- Eastern Cordillera (1)
- Eger Rift (1)
- Eisenbuhl (1)
- Ejina Basin (1)
- Ejina Becken (1)
- Electrical and electronic engineering (1)
- Eocene deformation (1)
- FT-ICR-MS (1)
- Fission-track thermochronology (1)
- Fluid Flow (1)
- Forest (1)
- Fracture and flow (1)
- Gabbroic rocks (1)
- Gaxun Nur (1)
- Germanic Basin (1)
- Granular materials (1)
- Gravity (1)
- Greenland ice cores (1)
- Ground-penetrating radar (1)
- Groundwater (1)
- HVSR (1)
- Heavy Minerals (1)
- Heihe (1)
- Heinrich event (1)
- Herkunftsanalyse (1)
- Herschel Island Qikiqtaruk (1)
- High frequency back-projection (1)
- High-temperature shear zones (1)
- Himalayas (1)
- Holocene (1)
- IGP (1)
- IPCC AR5 scenarios (1)
- IRD (1)
- InSAR (1)
- Induced seismicity (1)
- Ion-Beam (FIB) (1)
- Jafnayn formation (1)
- Kashmir earthquake (1)
- KiK-net (1)
- Kwajalein (1)
- L-Asterisk (1)
- Lake Malawi basin (1)
- Larix (1)
- Late Devonian (1)
- Late Triassic (1)
- Levee (1)
- Limestone assimilation in basanite (1)
- Lithospheric structure (1)
- Loss modeling (1)
- Lu-Hf in garnet (1)
- Major and trace elements (1)
- Meerfelder Maar (1)
- Methodology (1)
- Micromechanical model (1)
- Middle Eocene deformation (1)
- Mont Terri (1)
- Moraine landscape (1)
- NW Morocco (1)
- Natural hazard (1)
- Neo-Tethys (1)
- Neodymium isotopes (1)
- Neural circuits (1)
- Northern Sierras Pampeanas (1)
- P-T-d-tpaths (1)
- P-waves (1)
- PHREEQC (1)
- Palaeoceanography (1)
- Paleoceanography (1)
- Paris Agreement (1)
- Patagonia (1)
- Permafrost carbon feedback (1)
- Provenance Analysis (1)
- Puna Plateau (1)
- Qilian mountains (1)
- Quaternary (1)
- RCP scenarios (1)
- RCP4.5 (1)
- RCP8.5 (1)
- Reflectivity (1)
- Riverine floods (1)
- Roughness (1)
- SAR (1)
- SPEI (1)
- SPI (1)
- SWAT (1)
- Salt Range (1)
- Santa-Barbara system (1)
- Schwerminerale (1)
- Seismo-tectonics (1)
- Sentinel-1 (1)
- Shire River basin (1)
- Siberian arctic (1)
- Soil catena (1)
- Soil erosion (1)
- Source (1)
- South American monsoon system (1)
- South Asia (1)
- Southern Norway (1)
- Spatial analysis (1)
- Statistical methods (1)
- Strain localization Carrara marble (1)
- Subduction (1)
- Surface floods (1)
- Surface waves and free oscillations (1)
- Svalbard (1)
- Theoretical seismology (1)
- Time-series analysis (1)
- Tonminerale (1)
- Triaxial deformation (1)
- Uk'37 (1)
- Uncertainties (1)
- VIIRS DNB (1)
- Variscan (1)
- Vinoren (1)
- Visual system (1)
- Volcanic activity (1)
- WASA-SED (1)
- Wave (1)
- Weather (1)
- Western Europe (1)
- Western Interior Basin (1)
- Western Qaidam Basin (1)
- Yedoma Ice Complex (1)
- Zelezna hurka (1)
- Zircon U-Pb dating (1)
- absorption (1)
- accountability (1)
- acoustic emission (1)
- active layer (1)
- adaptation (1)
- aillikite (1)
- airglow (1)
- amphibole (1)
- archaeology (1)
- architecture (1)
- argon dating (1)
- artificial light (1)
- astronomical calibration (1)
- atlantic (1)
- atmospheric aerosols (1)
- automatic processing (1)
- bandgap (1)
- bedload transport (1)
- bedrock-alluvial channels (1)
- biodiversity (1)
- black carbon (1)
- borehole leakage (1)
- breaches (1)
- brown (1)
- calibration (1)
- carbon (1)
- carbon cycle (1)
- carbon density (1)
- carbon-cycle (1)
- case history (1)
- cement (1)
- central south island (1)
- change (1)
- channel width (1)
- characteristics (1)
- chronology (AICC2012) (1)
- classification (1)
- climate (1)
- climate action (1)
- climate networks (1)
- climate warming (1)
- clonal growth (1)
- coastal floods (1)
- cold surges (1)
- companies (1)
- composition (1)
- compound flood (1)
- conservative (1)
- continental collision (1)
- continental lithosphere (1)
- critical infrastructure resilience (1)
- crust (1)
- crustal magma chamber (1)
- daily-discharge time (1)
- damage model (1)
- damage surveys (1)
- data fusion (1)
- data-mining (1)
- decision making (1)
- decomposition (1)
- deep biosphere (1)
- deep carbon cycle (1)
- deep low-frequency earthquakes (1)
- deformation monitoring (1)
- design (1)
- diagnostics (1)
- diamond (1)
- diatoms (1)
- digital rock physics (1)
- discharge variability (1)
- distributed volcanic fields (1)
- distribution (1)
- disturbance tracking (1)
- drained lake basins (1)
- dust (1)
- dust sources (1)
- early eocene (1)
- earth system model (1)
- earthflow (1)
- earthquake (1)
- earthquake hazard analysis (1)
- earthquake modeling (1)
- eco-hydrology (1)
- ecology (1)
- ecosystem-based adaptation (1)
- ecosystems (1)
- edge-driven convection (1)
- elastic properties (1)
- electrical resistivity (1)
- electrochemistry (1)
- electrons (1)
- end-member mixing (1)
- energy security (1)
- eocene thermal maximum (1)
- epistemic (1)
- epsilon Nd (1)
- equifinality (1)
- evaporites (1)
- events (1)
- exhumation (1)
- extinction event (1)
- extreme discharge (1)
- extremes (1)
- facies (1)
- fatty acids (1)
- fault bend (1)
- fault bend fold (1)
- fault creep (1)
- flat-slab subduction (1)
- flexible model (1)
- flood (1)
- flood loss estimation (1)
- flood risk (1)
- floodplains (1)
- floods (1)
- flow-through experiment (1)
- fluid injection (1)
- fluvial floods (1)
- flux tracking (1)
- food security (1)
- forearc (1)
- forecast (1)
- formation damage (1)
- fracture aperture (1)
- freshwater ecosystems (1)
- frictional properties (1)
- frozen soil (1)
- functional traits (1)
- gas hydrate (1)
- gas-solid (1)
- generic sampling (1)
- geo-bio interaction (1)
- geochemical cycling (1)
- geochronology (1)
- geomorphology (1)
- geothermal energy (1)
- giant earthquakes (1)
- glacial and interglacial permafrost (1)
- glacial lakes (1)
- global (1)
- global monsoon (1)
- global sensitivity analysis (1)
- grain-size distribution (1)
- granulite (1)
- gravitational collapse (1)
- ground motion (1)
- ground motion models (1)
- ground response analysis (1)
- ground-motion modelling (1)
- groundwater vulnerability (1)
- growth rate (1)
- growth structures (1)
- gulf coastal plain (1)
- hanging wall (1)
- heat-stress (1)
- heavy-tailed distributions (1)
- hemolysins (1)
- heterogeneity (1)
- high-Pmetamorphism (1)
- high-frequency data (1)
- high-resolution biomarkers (1)
- high-resolution record (1)
- hybrid modelling (1)
- hydrological signatures (1)
- hysteresis (1)
- hysteresis pattern (1)
- ice (1)
- ice sheet (1)
- ice-grain mixture (1)
- ice-rafted debris (1)
- identifiability analysis (1)
- identifying influential nodes (1)
- ilmenite (1)
- imaging (1)
- impact forecasting (1)
- inference (1)
- insects (1)
- insights (1)
- internal structure (1)
- ionic liquids (1)
- lagoon (1)
- lagoons (1)
- lake Donggi Cona (1)
- lake drainage (1)
- lakes (1)
- landscape effect (1)
- landscape evolution (1)
- landscape transience (1)
- landslide (1)
- landslides (1)
- last glacial (1)
- last glacial period (1)
- lateglacial (1)
- lateral spread (1)
- lattice preferred orientations (1)
- leading edge (1)
- length measurements (1)
- liquids (1)
- lithosphere-mantle wind interactions (1)
- location of scatterers (1)
- low velocity zone (1)
- low-relief (1)
- machine learning (1)
- mafic migmatite (1)
- magma reservoirs (1)
- magma-poor (1)
- magnesium silicate skarn (1)
- magnetics (1)
- magnetosphere (1)
- mantle discontinuities (1)
- mantle flow modeling (1)
- margins (1)
- mass (1)
- mass spectrometry (1)
- maximum (1)
- mechanism (1)
- mega-monsoon (1)
- melilitite (1)
- meridional overturning circulation (1)
- mesoscale convective systems (1)
- metagenomic analysis (1)
- metal-containing ionic (1)
- metaophiolite (1)
- meteoric 10Be (1)
- meteorological drought (1)
- methane (1)
- methanogenic archaea (1)
- mica (1)
- micro-CT (1)
- microbial diversity (1)
- microbial lipid (1)
- migration (1)
- model structure (1)
- modeling (1)
- mofette systems (1)
- moment tensor (1)
- morphological traits (1)
- motivation theory (1)
- multi-slab (1)
- multi-year insurance (1)
- multivariable (1)
- natural hazards (1)
- near surface (1)
- network analysis (1)
- nightlights (1)
- nitrate export dynamic (1)
- nitrogen (1)
- non-state actors (1)
- north-western Argentina (1)
- northern Puna (1)
- numerical model (1)
- numerical modeling (1)
- numerics (1)
- ocean acidification (1)
- ocean circulation (1)
- oligocene climate (1)
- open data (1)
- optimization (1)
- orbital forcing (1)
- organic carbon cycle (1)
- origination (1)
- orogeny (1)
- orographic barrier (1)
- oxidizer (1)
- palaeoecology (1)
- paleoclimatology (1)
- paleosol sequence (1)
- pargasite (1)
- part 1 (1)
- partial melting (1)
- particulate organic carbon (1)
- passive microwave (1)
- pathogens (1)
- peatlands (1)
- pelitic migmatite (1)
- periglacial landscapes (1)
- permafrost degradation (1)
- permafrost regions (1)
- permafrost thaw (1)
- permeability test (1)
- petrostructural analysis (1)
- phlogopite (1)
- phreeqc (1)
- plasmidome (1)
- playa (1)
- playa-lake (1)
- plume-induced (1)
- pluvial floods (1)
- pollution (1)
- porphyry Cu-Au deposit (1)
- pressure solution (1)
- private governance (1)
- probabilistic modeling (1)
- propagation (1)
- protection (1)
- protists (1)
- quartz sandstone (1)
- radar interferometry (1)
- radial flow (1)
- radiation-damage (1)
- rainfall (1)
- ramp (1)
- random vibration (1)
- reaction (1)
- reactivation (1)
- record (1)
- reflection seismic (1)
- regional mapping (1)
- regionalisation (1)
- remote (1)
- residential buildings (1)
- response time (1)
- reversible scavenging (1)
- risk perception (1)
- risk reduction (1)
- river (1)
- river floods (1)
- river sediment (1)
- saline groundwater (1)
- saline pan cycle (1)
- saline permafrost (1)
- salt pan (1)
- san andreas fault (1)
- scale (1)
- scaling (1)
- scenario (1)
- sea (1)
- sea level rise (1)
- sea level trends (1)
- sea surface height (1)
- sea-level projections (1)
- seasonality (1)
- seasonality effect (1)
- security and economy (1)
- sediment (1)
- sediment routing (1)
- sediment transit time (1)
- sedimentary DNA (1)
- sedimentary ancient DNA (1)
- seismic anisotropy (1)
- seismic moment release (1)
- seismic risk (1)
- seismicity and tectonics (1)
- sensing (1)
- sequestration (1)
- series ERA-interim (1)
- serpentinite carbonation (1)
- shallow earthquakes (1)
- shotgun sequencing (1)
- signal propagation (1)
- silica concentration (1)
- singleslab (1)
- site amplification (1)
- skill (1)
- slip partitioning (1)
- social justice (1)
- soft law (1)
- soils (1)
- source array design (1)
- source time function (1)
- source-to-sink (1)
- southern Central Andes (1)
- southern hemisphere oceans (1)
- spatial correlation (1)
- special sensor microwave imager (1)
- special sensor microwave imager/sounder (1)
- speciation (1)
- spectral analysis (1)
- spinel (1)
- splitting parameters (1)
- stable (1)
- stable craton (1)
- stable isotopes (1)
- statistical seismology (1)
- stick slip (1)
- storm surge (1)
- strain localization (1)
- stratigraphy (1)
- streamflow (1)
- stress-tolerance genes (1)
- subduction initiation (PISI) (1)
- submarine (1)
- subsea (1)
- subsea permafrost (1)
- subsidence (1)
- sulphate metasomatism (1)
- summer monsoon (1)
- supercontinent Nuna (1)
- surprise (1)
- sustainable development (1)
- sustained casing pressure (1)
- synchronization (1)
- synthesis gas (1)
- systems (1)
- talik (1)
- tectonic inversion (1)
- temperature (1)
- tensile strength (1)
- tetrahalido metallates (1)
- theory (1)
- thermal (1)
- thermodynamic equilibrium (1)
- thermokarst lake (1)
- thermokarst lakes (1)
- thickness-displacement relationships (1)
- tibetan plateau (1)
- tide gauge observations (1)
- time-scale (1)
- time-series mapping (1)
- titanite (1)
- trait-based approaches (1)
- transfer (1)
- transit-time (1)
- transport processes (1)
- travel time distributions (1)
- treeline (1)
- trend analysis (1)
- tropical biodiversity (1)
- tropical swamp (1)
- uncertainty (1)
- underground coal gasification (1)
- upper tail behaviour (1)
- valley distribution (1)
- variance (1)
- virulence (1)
- water (1)
- water age (1)
- water security (1)
- water-energy-food nexus (1)
- waveform inversion (1)
- weakening mechanisms (1)
- weather (1)
Institute
- Institut für Geowissenschaften (151) (remove)
Width control on event-scale deposition and evacuation of sediment in bedrock-confined channels
(2020)
In mixed bedrock-alluvial rivers, the response of the system to a flood event can be affected by a number of factors, including coarse sediment availability in the channel, sediment supply from the hillslopes and upstream, flood sequencing and coarse sediment grain size distribution. However, the impact of along-stream changes in channel width on bedload transport dynamics remains largely unexplored. We combine field data, theory and numerical modelling to address this gap. First, we present observations from the Daan River gorge in western Taiwan, where the river flows through a 1 km long 20-50 m wide bedrock gorge bounded upstream and downstream by wide braidplains. We documented two flood events during which coarse sediment evacuation and redeposition appear to cause changes of up to several metres in channel bed elevation. Motivated by this case study, we examined the relationships between discharge, channel width and bedload transport capacity, and show that for a given slope narrow channels transport bedload more efficiently than wide ones at low discharges, whereas wider channels are more efficient at high discharges. We used the model sedFlow to explore this effect, running a random sequence of floods through a channel with a narrow gorge section bounded upstream and downstream by wider reaches. Channel response to imposed floods is complex, as high and low discharges drive different spatial patterns of erosion and deposition, and the channel may experience both of these regimes during the peak and recession periods of each flood. Our modelling suggests that width differences alone can drive substantial variations in sediment flux and bed response, without the need for variations in sediment supply or mobility. The fluctuations in sediment transport rates that result from width variations can lead to intermittent bed exposure, driving incision in different segments of the channel during different portions of the hydrograph.
Giant earthquakes with magnitudes above 8.5 occur only in subduction zones. Despite the developments made in observing large subduction zone earthquakes with geophysical instruments, the factors controlling the maximum size of these earthquakes are still poorly understood. Previous studies have suggested the importance of slab shape, roughness of the plate interface contact, state of the strain in the upper plate, thickness of sediments filling the trenches, and subduction rate. Here, we present 2-D cross-scale numerical models of seismic cycles for subduction zones with various geometries, subduction channel friction configurations, and subduction rates. We found that low-angle subduction and thick sediments in the subduction channel are the necessary conditions for generating giant earthquakes, while the subduction rate has a negligible effect. We suggest that these key parameters determine the maximum magnitude of a subduction earthquake by controlling the seismogenic zone width and smoothness of the subduction interface. This interpretation supports previous studies that are based upon observations and scaling laws. Our modeling results also suggest that low static friction in the sediment-filled subduction channel results in neutral or moderate compressive deformation in the overriding plate for low-angle subduction zones hosting giant earthquakes. These modeling results agree well with observations for the largest earthquakes. Based on our models we predict maximum magnitudes of subduction earthquakes worldwide, demonstrating the fit to magnitudes of all giant earthquakes of the 20th and 21st centuries and good agreement with the predictions based on statistical analyses of observations.
The Victoria microplate between the Eastern and Western Branches of the East African Rift System is one of the largest continental microplates on Earth. In striking contrast to its neighboring plates, Victoria rotates counterclockwise with respect to Nubia. The underlying cause of this distinctive rotation has remained elusive so far. Using 3D numerical models, we investigate the role of pre-existing lithospheric heterogeneities in continental microplate rotation. We find that Victoria's rotation is primarily controlled by the distribution of rheologically stronger zones that transmit the drag of the major plates to the microplate and of the mechanically weaker mobile belts surrounding Victoria that facilitate rotation. Our models reproduce Victoria's GPS-derived counterclockwise rotation as well as key complexities of the regional tectonic stress field. These results reconcile competing ideas on the opening of the rift system by highlighting differences in orientation of the far-field divergence, local extension, and the minimum horizontal stress. One of the largest continental microplates on Earth is situated in the center of the East African Rift System, and oddly, the Victoria microplate rotates counterclockwise with respect to the neighboring African tectonic plate. Here, the authors' modelling results suggest that Victoria microplate rotation is caused by edge-driven lithospheric processes related to the specific geometry of rheologically weak and strong regions.
Rivers transfer particulate organic carbon (POC) from eroding mountains into geological sinks. Organic carbon source composition and selective mobilization have been shown to affect the type and quantity of POC export, but their combined effects across complex mountain ranges remain underexplored. Here, we examine the variation in organic carbon sourcing and transport in the trans-Himalayan Kali Gandaki River catchment, along strong gradients in precipitation, rock type and vegetation. Combining bulk stable nitrogen, and stable and radioactive organic carbon isotopic composition of bedrock, litter, soil and river sediment samples with a Bayesian end-member mixing approach, we differentiate POC sources along the river and quantify their export. Our analysis shows that POC export from the Tibetan segment of the catchment, where carbon bearing shales are partially covered by aged and modern soils, is dominated by petrogenic POC. Based on our data we re-assess the presence of aged biospheric OC in this part of the catchment, and its contribution to the river load. In the High Himalayan segment, we observed low inputs of petrogenic and biospheric POC, likely due to very low organic carbon concentrations in the metamorphic bedrock, combined with erosion dominated by deep-seated landslides. Our findings show that along the Kali Gandaki River, the sourcing of sediment and organic carbon are decoupled, due to differences in rock organic carbon content, soil and above ground carbon stocks, and geomorphic process activity. While the fast eroding High Himalayas are the principal source of river sediment, the Tibetan headwaters, where erosion rates are lower, are the principal source of organic carbon. To robustly estimate organic carbon export from the Himalayas, the mountain range should be divided into tectono-physiographic zones with distinct organic carbon yields due to differences in substrate and erosion processes and rates.
Variation of deuterium excess in surface waters across a 5000-m elevation gradient in eastern Nepal
(2020)
The strong elevation gradient of the Himalaya allows for investigation of altitude and orographic impacts on surface water delta O-18 and delta D stable isotope values. This study differentiates the time- and altitude-variable contributions of source waters to the Arun River in eastern Nepal. It provides isotope data along a 5000-m gradient collected from tributaries as well as groundwater, snow, and glacial-sourced surface waters and time-series data from April to October 2016. We find nonlinear trends in delta O-18 and delta D lapse rates with high-elevation lapse rates (4000-6000 masl) 5-7 times more negative than low-elevation lapse rates (1000-3000 masl). A distinct seasonal signal in delta O-18 and delta D lapse rates indicates time-variable source-water contributions from glacial and snow meltwater as well as precipitation transitions between the Indian Summer Monsoon and Winter Westerly Disturbances. Deuterium excess correlates with the extent of snowpack and tracks melt events during the Indian Summer Monsoon season. Our analysis identifies the influence of snow and glacial melt waters on river composition during low-flow conditions before the monsoon (April/May 2016) followed by a 5-week transition to the Indian Summer Monsoon-sourced rainfall around mid-June 2016. In the post-monsoon season, we find continued influence from glacial melt waters as well as ISM-sourced groundwater.
The southern margin of the Central Anatolian Plateau (CAP) records a strong uplift phase after the early Middle Pleistocene, which has been related to the slab break-off of the subducting Arabian plate beneath the Anatolian microplate. During the last 450 kyr the area underwent an uplift phase at a mean rate of similar to 3.2 m/kyr, as suggested by Middle Pleistocene marine sediments exposed at similar to 1,500 m above sea level. These values are significantly higher than the 1.0-1.5 m/kyr estimated since the Late Pleistocene, suggesting temporal variations in uplift rate. To estimate changes in uplift rate during the Pleistocene we studied the marine terraces along the CAP southern margin, mapping the remnants of the platforms and their associated deposits in the field, and used the TerraceM software to identify the position and elevation of associated shoreline angles. We used shoreline angles and the timing of Quaternary marine sedimentation as constrains for a Landscape Evolution Model that simulates wave erosion of an uplifting coast. We applied random optimization algorithms and minimization statistics to find the input parameters that better reproduce the morphology of CAP marine terraces. The best-fitting uplift rate history suggests a significative increase from 1.9 to 3.5 m/kyr between 500 and 200 kyr, followed by an abrupt decrease to 1.4 m/kyr until the present. Our results agree with slab break-off models, which suggest a strong uplift pulse during slab rupture followed by a smoother decrease.
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.
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.
Cities will play a key role in the grand challenge of nourishing a growing global population, because, due to their population density, they set the demand. To ensure that food systems are sustainable, as well as nourishing, one solution often suggested is to shorten their supply chains toward a regional rather than a global basis. While such regional systems may have a range of costs and benefits, we investigate the mitigation potential of regionalized urban food systems by examining the greenhouse gas emissions associated with food transport. Using data on food consumption for 7108 urban administrative units (UAUs), we simulate total transport emissions for both regionalized and globalized supply chains. In regionalized systems, the UAUs' demands are fulfilled by peripheral food production, whereas to simulate global supply chains, food demand is met from an international pool (where the origin can be any location globally). We estimate that regionalized systems could reduce current emissions from food transport. However, because longer supply chains benefit from maximizing comparative advantage, this emission reduction would require closing yield gaps, reducing food waste, shifting toward diversified farming, and consuming seasonal produce. Regionalization of food systems will be an essential component to limit global warming to well below 2 degrees C in the future.
The chapter illustrates how cities can be understood from a system–theory perspective as complex social systems. It argues that the classical and linear intervention methods are often no (longer) suitable for the complex structure, temporal dynamics and multifaceted processuality of urban development. It offers a systemic and systems theory-inspired method as an alternative approach, which allows for extended possibilities that are more appropriate for dealing with urban development processes. The method was developed on the basis of practical experience and theoretical insights. The approach should demonstrate for local decision-makers potential areas of activity for organising urban changes through co-design.
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 Indus Molasse records orogenic sedimentation associated with uplift and erosion of the southern margin of Asia in the course of ongoing India-Eurasia collision. Detailed field investigation clarifies the nature and extent of the depositional contact between this molasse and the underlying basement units. We report the first dataset on detrital zircon U-Pb ages, Hf isotopes and apatite U-Pb ages for the autochthonous molasse in the Indus Suture Zone. A latest Oligocene depositional age is proposed on the basis of the youngest detrital zircon U-Pb age peak and is consistent with published biostratigraphic data. Multiple provenance indicators suggest exclusively northerly derivation with no input from India in the lowermost parts of the section. The results provide constraints on the uplift and erosion history of the Ladakh Range following the initial India-Asia collision.
In recent years, nature-based solutions are receiving increasing attention in the field of disaster risk reduction and climate change adaptation as inclusive, no regret approaches. Ecosystem-based adaptation (EbA) can mitigate the impacts of climate change, build resilience and tackle environmental degradation thereby supporting the targets set by the 2030 Agenda, the Paris Agreement and the Sendai Framework. Despite these benefits, EbA is still rarely implemented in practice. To better understand the barriers to implementation, this research examines policy-makers' perceptions of EbA, using an extended version of Protection Motivation Theory as an analytical framework. Through semi-structured interviews with policy-makers at regional and provincial level in Central Vietnam, it was found that EbA is generally considered a promising response option, mainly due to its multiple ecosystem-service benefits. The demand for EbA measures was largely driven by the perceived consequences of natural hazards and climate change. Insufficient perceived response efficacy and time-lags in effectiveness for disaster risk reduction were identified as key impediments for implementation. Pilot projects and capacity building on EbA are important means to overcome these perceptual barriers. This paper contributes to bridging the knowledge-gap on political decision-making regarding EbA and can, thereby, promote its mainstreaming into policy plans.
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.
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.
As the Arctic coast erodes, it drains thermokarst lakes, transforming them into lagoons, and, eventually, integrates them into subsea permafrost. Lagoons represent the first stage of a thermokarst lake transition to a marine setting and possibly more saline and colder upper boundary conditions. In this research, borehole data, electrical resistivity surveying, and modeling of heat and salt diffusion were carried out at Polar Fox Lagoon on the Bykovsky Peninsula, Siberia. Polar Fox Lagoon is a seasonally isolated water body connected to Tiksi Bay through a channel, leading to hypersaline waters under the ice cover. The boreholes in the center of the lagoon revealed floating ice and a saline cryotic bed underlain by a saline cryotic talik, a thin ice-bearing permafrost layer, and unfrozen ground. The bathymetry showed that most of the lagoon had bedfast ice in spring. In bedfast ice areas, the electrical resistivity profiles suggested that an unfrozen saline layer was underlain by a thick layer of refrozen talik. The modeling showed that thermokarst lake taliks can refreeze when submerged in saltwater with mean annual bottom water temperatures below or slightly above 0 degrees C. This occurs, because the top-down chemical degradation of newly formed ice-bearing permafrost is slower than the refreezing of the talik. Hence, lagoons may precondition taliks with a layer of ice-bearing permafrost before encroachment by the sea, and this frozen layer may act as a cap on gas migration out of the underlying talik.
Thermal erosion is a major mechanism of permafrost degradation, resulting in characteristic landforms. We inventory thermo-erosional valleys in ice-rich coastal lowlands adjacent to the Siberian Laptev Sea based on remote sensing, Geographic Information System (GIS), and field investigations for a first regional assessment of their spatial distribution and characteristics. Three study areas with similar geological (Yedoma Ice Complex) but diverse geomorphological conditions vary in valley areal extent, incision depth, and branching geometry. The most extensive valley networks are incised deeply (up to 35 m) into the broad inclined lowland around Mamontov Klyk. The flat, low-lying plain forming the Buor Khaya Peninsula is more degraded by thermokarst and characterized by long valleys of lower depth with short tributaries. Small, isolated Yedoma Ice Complex remnants in the Lena River Delta predominantly exhibit shorter but deep valleys. Based on these hydrographical network and topography assessments, we discuss geomorphological and hydrological connections to erosion processes. Relative catchment size along with regional slope interact with other Holocene relief-forming processes such as thermokarst and neotectonics. Our findings suggest that thermo-erosional valleys are prominent, hitherto overlooked permafrost degradation landforms that add to impacts on biogeochemical cycling, sediment transport, and hydrology in the degrading Siberian Yedoma Ice Complex.
Theories and concepts
(2020)
Urban Change and Citizenship in Times of Crisis addresses the fact that in the beginning of the twenty-first century the majority of the world’s population is urbanised, a social fact that has turned cities more than ever into focal sites of social change. Multiple economic and political strategies, employed by a variety of individual and collective actors, on a number of scales, constitute cities as contested spaces that hold opportunities as well as restrictions for their inhabitants. While cities and urban spaces have long been of central concern for the social sciences, today, classical sociological questions about the city acquire new meaning: Can cities be spaces of emancipation, or does life in the modern city entail a corrosion of citizenship rights? Is the city the focus of societal transformation processes, or do urban environments lose importance in shaping social reality and economic relationships? Furthermore, new questions urgently need to be asked: What is the impact of different historical phenomena such as neo-liberal restructuring, financial and economic crises, or migration flows, as well as their respective counter-movements, on the structure of contemporary cities and on the citizenship rights of city inhabitants? The three volumes address such crucial questions thereby opening up new spaces of debate on both the city and new developments of urbanism.
The contributions to Theories and Concepts offer new theoretical reflections on the city in a philosophical and historical perspective as well as fresh empirical analyses of social life in urban contexts. Chapters not only critically revisit classical and modern philosophical considerations about the nature of cities but no less discuss normative philosophical reflections of urban life and the role of religion in historical processes of the emergence of cities. Composed around the question whether there can be such a thing as a ‘successful city’, this volume addresses issues of urban political subjectivities by considering the city’s role in historical processes of emancipation, the fight for citizenship rights, and today’s challenges and opportunities with regard to promoting social justice, integration, and diversity. Consequentially, theory-driven empirical analyses offer new insight into ways of solving problems in urban contexts and a genuine approach to analyse the Social Quality in cities.
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.
During the period 750-600 Ma ago, prior to the final break-up of the supercontinent Rodinia, the crust of both the North American Craton and Baltica was intruded by significant amounts of rift-related magmas originating from the mantle. In the Proterozoic crust of Southern Norway, the 580 Ma old Fen carbonatite-ultramafic complex is a representative of this type of rocks. In this paper, we report the occurrence of an ultramafic lamprophyre dyke which possibly is linked to the Fen complex, although Ar-40/Ar-39 data from phenocrystic phlogopite from the dyke gave an age of 686 +/- 9 Ma. The lamprophyre dyke was recently discovered in one of the Kongsberg silver mines at Vinoren, Norway. Whole rock geochemistry, geochronological and mineralogical data from the ultramafic lamprophyre dyke are presented aiming to elucidate its origin and possible geodynamic setting. From the whole-rock composition of the Vinoren dyke, the rock could be recognized as transitional between carbonatite and kimberlite-II (orangeite). From its diagnostic mineralogy, the rock is classified as aillikite. The compositions and xenocrystic nature of several of the major and accessory minerals from the Vinoren aillikite are characteristic for diamondiferous rocks (kimberlites/lamproites/UML): Phlogopite with kinoshitalite-rich rims, chromite-spinel-ulvospinel series, Mg- and Mn-rich ilmenites, rutile and lucasite-(Ce). We suggest that the aillikite melt formed during partial melting of a MARID (mica-amphibole-rutile-ilmenite-diopside)-like source under CO2 fluxing. The pre-rifting geodynamic setting of the Vinoren aillikite before the Rodinia supercontinent breakup suggests a relatively thick SCLM (Subcontinental Lithospheric Mantle) during this stage and might indicate a diamond-bearing source for the parental melt. This is in contrast to the about 100 Ma younger Fen complex, which were derived from a thin SCLM.
Tensile strength is an important parameter when it comes to predictions of potential fracturing of sediments by natural processes such as the emplacement of ice or gas hydrate lenses, as well as anthropogenic fracturing or else the stability of engineering constructions such as boreholes. Yet, tensile strength (sigma(tau)) measurements of unconsolidated ice-bearing or gas hydrate-bearing sands are scarce and affected by a large variability. <br /> In the course of the SUGAR project we successfully used ice as a model for pore-filling and "load-bearing" gas hydrate in sand to determine compressional wave velocity. We were thus able to verify comparable formation characteristics and morphologies of ice and gas hydrate within the pore space. As these are important values for the tensile strength of ice/hydrate-bearing sands, ice was also used as a model for hydrate-bearing sands, despite differences in the mechanical behavior and strength of pure ice and gas hydrate. Water-saturated sand cores with ice saturations (S-ice) between 0 and 100% were tested at -6.8 degrees C. The varying S-ice were a result of the freezing point depression caused by saline solutions of different concentrations. The sigma(tau) was directly determined using a sleeve-fracturing test with an internal pressure that was created within the frozen samples. The setup was also adapted to fit a pressure vessel for tests using confining pressure. <br /> The correlation of S-ice - sigma(tau) shows an exponential increase of sigma(tau) with S-ice. Whereas at S-ice < 60% the increase is small, it is large at S-ice > 80%. In conjunction with the change in strength, the viscoelastic behavior changes. A clear peak strength occurs at S-ice > 80%. We conclude that given 60% < S-ice < 80% the pore-filling morphology of the ice converts into a frame-building habitus and at S-ice > 80% the frame gains strength while the amount of residual water decreases. Tensile failure and cracking now exceed grain boundary sliding as the prevailing failure mode. The ice morphology in the sand is non-cementing and comparable to a gas hydrate-sand mixture.
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.
The affordability of property-level adaptation measures against flooding is crucial due to the movement toward integrated flood risk management, which requires the individuals threatened by flooding to actively manage flooding. It is surprising to find that affordability is not often discussed, given the important roles that affordability and social justice play regarding flood risk management. This article provides a starting point for investigating the potential rate of unaffordability of flood risk property-level adaptation measures across Europe using two definitions of affordability, which are combined with two different affordability thresholds from within flood risk research. It uses concepts of investment and payment affordability, with affordability thresholds based on residual income and expenditure definitions of unaffordability. These concepts, in turn, are linked with social justice through fairness concerns, in that, all should have equal capability to act, of which affordability is one avenue. In doing so, it was found that, for a large proportion of Europe, property owners generally cannot afford to make one-time payment of the cost of protective measures. These can be made affordable with installment payment mechanisms or similar mechanisms that spread costs over time. Therefore, the movement toward greater obligations for flood-prone residents to actively adapt to flooding should be accompanied by socially accessible financing mechanisms.
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.
This review provides a synthesis of current knowledge on the morphological and functional traits of testate amoebae, a polyphyletic group of protists commonly used as proxies of past hydrological changes in paleoecological investigations from peatland, lake sediment and soil archives. A trait-based approach to understanding testate amoebae ecology and paleoecology has gained in popularity in recent years, with research showing that morphological characteristics provide complementary information to the commonly used environmental inferences based on testate amoeba (morpho-)species data. We provide a broad overview of testate amoeba morphological and functional traits and trait-environment relationships in the context of ecology, evolution, genetics, biogeography, and paleoecology. As examples we report upon previous ecological and paleoecological studies that used trait-based approaches, and describe key testate amoebae traits that can be used to improve the interpretation of environmental studies. We also highlight knowledge gaps and speculate on potential future directions for the application of trait-based approaches in testate amoeba research.
Anthropogenic activities have led to a global decline in biodiversity, and monitoring studies indicate that both insect communities and wetland ecosystems are particularly affected. However, there is a need for long-term data (over centennial or millennial timescales) to better understand natural community dynamics and the processes that govern the observed trends. Chironomids (Insecta: Diptera: Chironomidae) are often the most abundant insects in lake ecosystems, sensitive to environmental change, and, because their larval exoskeleton head capsules preserve well in lake sediments, they provide a unique record of insect community dynamics through time. Here, we provide the results of a metadata analysis of chironomid diversity across a range of spatial and temporal scales. First, we analyse spatial trends in chironomid diversity using Northern Hemispheric data sets overall consisting of 837 lakes. Our results indicate that in most of our data sets, summer temperature (T-jul) is strongly associated with spatial trends in modern-day chironomid diversity. We observe a strong increase in chironomid alpha diversity with increasing T-jul in regions with present-day T-jul between 2.5 and 14 degrees C. In some areas with T-jul > 14 degrees C, chironomid diversity stabilizes or declines. Second, we demonstrate that the direction and amplitude of change in alpha diversity in a compilation of subfossil chironomid records spanning the last glacial-interglacial transition (similar to 15,000-11,000 years ago) are similar to those observed in our modern data. A compilation of Holocene records shows that during phases when the amplitude of temperature change was small, site-specific factors had a greater influence on the chironomid fauna obscuring the chironomid diversity-temperature relationship. Our results imply expected overall chironomid diversity increases in colder regions such as the Arctic under sustained global warming, but with complex and not necessarily predictable responses for individual sites.
Underground coal gasification (UCG) is an in situ conversion technique that enables the production of high-calorific synthesis gas from resources that are economically not minable by conventional methods. A broad range of end-use options is available for the synthesis gas, including fuels and chemical feedstock production. Furthermore, UCG also offers a high potential for integration with Carbon Capture and Storage (CCS) to mitigate greenhouse gas emissions. In the present study, a stoichiometric equilibrium model, based on minimization of the Gibbs function has been used to estimate the equilibrium composition of the synthesis gas. Thereto, we further developed and applied a proven thermodynamic equilibrium model to simulate the relevant thermochemical coal conversion processes (pyrolysis and gasification). Our modeling approach has been validated against thermodynamic models, laboratory gasification experiments and UCG field trial data reported in the literature. The synthesis gas compositions have been found to be in good agreement under a wide range of different operating conditions. Consequently, the presented modeling approach enables an efficient quantification of synthesis gas quality resulting from UCG, considering varying coal and oxidizer compositions at deposit-specific pressures and temperatures.
It has recently been demonstrated that the interaction of a mantle plume with sufficiently old oceanic lithosphere can initiate subduction. However, the existence of large lithospheric heterogeneities, such as a buoyant plateau, in proximity to a rising plume head may potentially hinder the formation of a new subduction zone. Here, we investigate this scenario by means of 3-D numerical thermomechanical modeling. We explore how plume-lithosphere interaction is affected by lithospheric age, relative location of plume head and plateau border, and the strength of the oceanic crust. Our numerical experiments suggest four different geodynamic regimes: (a) oceanic trench formation, (b) circular oceanic-plateau trench formation, (c) plateau trench formation, and (d) no trench formation. We show that regardless of the age and crustal strength of the oceanic lithosphere, subduction can initiate when the plume head is either below the plateau border or at a distance less than the plume radius from the plateau edge. Crustal heterogeneity facilitates subduction initiation of old oceanic lithosphere. High crustal strength hampers the formation of a new subduction zone when the plume head is located below a young lithosphere containing a thick and strong plateau. We suggest that plume-plateau interaction in the western margin of the Caribbean could have resulted in subduction initiation when the plume head impinged onto the oceanic lithosphere close to the border between plateau and oceanic crust.
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.
Groningen is the largest onshore gas field under production in Europe. The pressure depletion of the gas field started in 1963. In 1991, the first induced micro-earthquakes have been located at reservoir level with increasing rates in the following decades. Most of these events are of magnitude less than 2.0 and cannot be felt. However, maximum observed magnitudes continuously increased over the years until the largest, significant event with ML=3.6 was recorded in 2014, which finally led to the decision to reduce the production. This causal sequence displays the crucial role of understanding and modeling the relation between production and induced seismicity for economic planing and hazard assessment. Here we test whether the induced seismicity related to gas exploration can be modeled by the statistical response of fault networks with rate-and-state-dependent frictional behavior. We use the long and complete local seismic catalog and additionally detailed information on production-induced changes at the reservoir level to test different seismicity models. Both the changes of the fluid pressure and of the reservoir compaction are tested as input to approximate the Coulomb stress changes. We find that the rate-and-state model with a constant tectonic background seismicity rate can reproduce the observed long delay of the seismicity onset. In contrast, so-called Coulomb failure models with instantaneous earthquake nucleation need to assume that all faults are initially far from a critical state of stress to explain the delay. Our rate-and-state model based on the fluid pore pressure fits the spatiotemporal pattern of the seismicity best, where the fit further improves by taking the fault density and orientation into account. Despite its simplicity with only three free parameters, the rate-and-state model can reproduce the main statistical features of the observed activity.
Forested areas are assumed not to be influenced by erosion processes. However, forest soils of Northern Germany in a hummocky ground moraine landscape can sometimes exhibit a very shallow thickness on crest positions and buried soils on slope positions. The question consequently is: Are these on-going or ancient erosional and depositional processes? Plutonium isotopes act as soil erosion/deposition tracers for recent (last few decades) processes. Here, we quantified the 239+240PU inventories in a small, forested catchment (ancient forest "Melzower Forst", deciduous trees), which is characterised by a hummocky terrain including a kettle hole. Soil development depths (depth to C horizon) and 239+240PU inventories along a catena of sixteen different profiles were determined and correlated to relief parameters. Moreover, we compared different modelling approaches to derive erosion rates from Pu data. <br /> We find a strong relationship between soil development depths, distance-to-sink and topography along the catena. Fully developed Retisols (thicknesses > 1 m) in the colluvium overlay old land surfaces as documented by fossil Ah horizons. However, we found no relationship of Pu-based erosion rates to any relief parameter. Instead, 239+240PU inventories showed a very high local, spatial variability (36-70 Bq m(-2)). Low annual rainfall, spatially distributed interception and stem flow might explain the high variability of the 239+240PU inventories, giving rise to a patchy input pattern. Different models resulted in quite similar erosion and deposition rates (max: -5 t ha(-1) yr(-1) to +7.3 t ha(-1) yr(-1)). Although some rates are rather high, the magnitude of soil erosion and deposition - in terms of soil thickness change - is negligible during the last 55 years. The partially high values are an effect of the patchy Pu deposition on the forest floor. This forest has been protected for at least 240 years. Therefore rather natural events and anthropogenic activities during medieval times or even earlier must have caused the observed soil pattern, which documents strong erosion and deposition processes.
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.
AM(W)= 5.1 earthquake on January 21st, 2016 marked the beginning of a significant seismic sequence in the southern Alboran Sea, culminating in aM(W)= 6.3 earthquake on January 25th, and continuing with further moderate magnitude earthquakes until March. We use data from 35 seismic broadband stations in Spain, Morocco and Portugal to relocate the seismicity, estimate seismic moment tensors, and isolate regional apparent source time functions for the main earthquake. Relocation and regional moment tensor inversion consistently yield very shallow depths for the majority of events. We obtain 50 moment tensors for the sequence, showing a mixture of strike-slip faulting for the foreshock and the main event and reverse faulting for the major aftershocks. The leading role of reverse focal mechanisms among the aftershocks may be explained by the geometry of the fault network. The mainshock nucleates at a bend along the left-lateral Al-Idrisi fault, introducing local transpression within the transtensional Alboran Basin. The shallow depths of the 2016 Alboran Sea earthquakes may favor slip-partitioning on the involved faults. Apparent source durations for the main event suggest a similar to 21 km long, asymmetric rupture that propagates primarily toward NE into the restraining fault segment, with fast rupture speed of similar to 3.0 km/s. Consistently, the inversion for laterally variable fault displacement situates the main slip in the restraining segment. The partitioning into strike-slip rupture and dip-slip aftershocks confirms a non-optimal orientation of this segment, and suggests that the 2016 event settled a slip deficit from previous ruptures that could not propagate into the stronger restraining segment.
Diffusive transport and sorption processes of uranium in the Swiss Opalinus Clay were investigated as a function of partial pressure of carbon dioxide pCO(2), varying mineralogy in the facies and associated changes in porewater composition. Simulations were conducted in one-dimensional diffusion models on the 100 m-scale for a time of one million years using a bottom-up approach based on mechanistic surface complexation models as well as cation exchange to quantify sorption. Speciation calculations have shown, uranium is mainly present as U(VI) and must therefore be considered as mobile for in-situ conditions. Uranium migrated up to 26 m in both, the sandy and the carbonate-rich facies, whereas in the shaly facies 16 m was the maximum. The main species was the anionic complex CaUO2(CO3)(3)(2-) . Hence, anion exclusion was taken into account and further reduced the migration distances by 30 %. The concentrations of calcium and carbonates reflected by the set pCO(2) determine speciation and activity of uranium and consequently the sorption behaviour. Our simulation results allow for the first time to prioritize on the far-field scale the governing parameters for diffusion and sorption of uranium and hence outline the sensitivity of the system. Sorption processes are controlled in descending priority by the carbonate and calcium concentrations, pH, pe and the clay mineral content. Therefore, the variation in porewater composition resulting from the heterogeneity of the facies in the Opalinus Clay formation needs to be considered in the assessment of uranium migration in the far field of a potential repository.
Earthquakes often rupture across more than one fault segment. If such rupture segmentation occurs on a significant scale, a simple point-source or one-fault model may not represent the rupture process well. As a consequence earthquake characteristics inferred, based on one-source assumptions, may become systematically wrong. This might have effects on follow-up analyses, for example regional stress field inversions and seismic hazard assessments. While rupture segmentation is evident for most M-w > 7 earthquakes, also smaller ones with 5.5 < M-w < 7 can be segmented. We investigate the sensitivity of globally available data sets to rupture segmentation and their resolution to reliably estimate the mechanisms in presence of segmentation. We focus on the sensitivity of InSAR (Interferometric Synthetic Aperture Radar) data in the static near-field and seismic waveforms in the far-field of the rupture and carry out non-linear and Bayesian optimizations of single-source and two-sources kinematic models (double-couple point sources and finite, rectangular sources) using InSAR and teleseismic waveforms separately. Our case studies comprises of four M-w 6-7 earthquakes: the 2009 L'Aquila and 2016 Amatrice (Italy) and the 2005 and 2008 Zhongba (Tibet) earthquakes. We contrast the data misfits of different source complexity by using the Akaike informational criterion (AIC). We find that the AIC method is well suited for data-driven inferences on significant rupture segmentation for the given data sets. This is based on our observation that an AIC-stated significant improvement of data fit for two-segment models over one-segment models correlates with significantly different mechanisms of the two source segments and their average compared to the single-segment mechanism. We attribute these modelled differences to a sufficient sensitivity of the data to resolve rupture segmentation. Our results show that near-field data are generally more sensitive to rupture segmentation of shallow earthquakes than far-field data but that also teleseismic data can resolve rupture segmentation in the studied magnitude range. We further conclude that a significant difference in the modelled source mechanisms for different segmentations shows that an appropriate choice of model segmentation matters for a robust estimation of source mechanisms. It reduces systematic biases and trade-off and thereby improves the knowledge on the rupture. Our study presents a strategy and method to detect significant rupture segmentation such that an appropriate model complexity can be used in the source mechanism inference. A similar, systematic investigation of earthquakes in the range of M-w 5.5-7 could provide important hazard-relevant statistics on rupture segmentation. In these cases single-source models introduce a systematic bias. Consideration of rupture segmentation therefore matters for a robust estimation of source mechanisms of the studied earthquakes.
The Quaternary volcanic fields of the Eifel (Rhineland-Palatinate, Germany) had their last eruptions less than 13,000 years ago. Recently, deep low-frequency (DLF) earthquakes were detected beneath one of the volcanic fields showing evidence of ongoing magmatic activity in the lower crust and upper mantle. In this work, seismic wide- and steep-angle experiments from 1978/1979 and 1987/1988 are compiled, partially reprocessed and interpreted, together with other data to better determine the location, size, shape, and state of magmatic reservoirs in the Eifel region near the crust-mantle boundary. We discuss seismic evidence for a low-velocity gradient layer from 30-36 km depth, which has developed over a large region under all Quaternary volcanic fields of the Rhenish Massif and can be explained by the presence of partial melts. We show that the DLF earthquakes connect the postulated upper mantle reservoir with the upper crust at a depth of about 8 km, directly below one of the youngest phonolitic volcanic centers in the Eifel, where CO(2)originating from the mantle is massively outgassing. A bright spot in the West Eifel between 6 and 10 km depth represents a Tertiary magma reservoir and is seen as a model for a differentiated reservoir beneath the young phonolitic center today. We find that the distribution of volcanic fields is controlled by the Variscan lithospheric structures and terrane boundaries as a whole, which is reflected by an offset of the Moho depth, a wedge-shaped transparent zone in the lower crust and the system of thrusts over about 120 km length.
Analysis of past and present stimulation projects reveals that the temporal evolution and growth of maximum observed moment magnitudes may be linked directly to the injected fluid volume and hydraulic energy. Overall evolution of seismic moment seems independent of the tectonic stress regime and is most likely governed by reservoir specific parameters, such as the preexisting structural inventory. Data suggest that magnitudes can grow either in a stable way, indicating the constant propagation of self-arrested ruptures, or unbound, for which the maximum magnitude is only limited by the size of tectonic faults and fault connectivity. Transition between the two states may occur at any time during injection or not at all. Monitoring and traffic light systems used during stimulations need to account for the possibility of unstable rupture propagation from the very beginning of injection by observing the entire seismicity evolution in near-real time and at high resolution for an immediate reaction in injection strategy.
Plain Language Summary Predicting and controlling the size of earthquakes caused by fluid injection is currently the major concern of many projects associated with geothermal energy production. Here, we analyze the magnitude and seismic moment evolution with injection parameters for prominent geothermal and scientific projects to date. Evolution of seismicity seems to be largely independent of the tectonic stress background and seemingly depends on reservoir specific characteristics. We find that the maximum observed magnitudes relate linearly to the injected volume or hydraulic energy. A linear relation suggests stable growth of induced ruptures, as predicted by current models, or rupture growth may no longer depend on the stimulated volume but on tectonics. A system may change between the two states during the course of fluid injection. Close-by and high-resolution monitoring of seismic and hydraulic parameters in near-real time may help identify these fundamental changes in ample time to change injection strategy and manage maximum magnitudes.
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.
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.
The Brazilian Earth System Model (BESM-OA2.5), while simulating the historical period proposed by the fifth phase of the Coupled Model Intercomparison Project (CMIP5), detects an increasing trend in the sea surface height (SSH) on the southern hemisphere oceans relative to that of the pre-industrial era. The increasing trend is accentuated in the CMIP5 RCP4.5 and RCP8.5 future scenarios with higher concentrations of greenhouse gases in the atmosphere. This study sheds light on the sources of such trends in these regions. The results suggest an association with the thermal expansion of the oceans in the upper 700 m due to a gradual warming inflicted by those future scenarios. BESM-OA2.5 presents a surface height increase of 0.11 m in the historical period of 1850-2005. Concerning future projections, BESM-OA2.5 projects SSH increases of 0.14 and 0.23 m (relative to the historical 2005 value) for RCP4.5 and RCP8.5, respectively, by the end of 2100. These increases are predominantly in a band of latitude within 35-60 degrees S in the Atlantic and Indian oceans. The reproducibility of the trend signal detected in the BESM-OA2.5 simulations is confirmed by the results of three other CMIP5 models.
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.
In the comment on "Varves of the Dead Sea sedimentary record." Quaternary Science Reviews 215 (Ben Dor et al., 2019): 173-184. by R. Bookman, two recently published papers are suggested to prove that the interpretation of the laminated sedimentary sequence of the Dead Sea, deposited mostly during MIS2 and Holocene pluvials, as annual deposits (i.e., varves) is wrong. In the following response, we delineate several lines of evidence which coalesce to demonstrate that based on the vast majority of evidence, including some of the evidence provided in the comment itself, the interpretation of these sediments as varves is the more likely scientific conclusion. We further discuss the evidence brought up in the comment and its irrelevance and lack of robustness for addressing the question under discussion.
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.
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.
Probabilistic assessment of seismic hazard and risk over a geographical region presents the modeler with challenges in the characterization of the site amplification that are not present in site-specific assessment. Using site-to-site residuals from a ground motion model fit to observations from the Japanese KiK-net database, correlations between measured local amplifications and mappable proxies such as topographic slope and geology are explored. These are used subsequently to develop empirical models describing amplification as a direct function of slope, conditional upon geological period. These correlations also demonstrate the limitations of inferring 30-m shearwave velocity from slope and applying them directly into ground motion models. Instead, they illustrate the feasibility of deriving spectral acceleration amplification factors directly from sets of observed records, which are calibrated to parameters that can be mapped uniformly on a regional scale. The result is a geologically calibrated amplification model that can be incorporated into national and regional seismic hazard and risk assessment, ensuring that the corresponding total aleatory variability reflects the predictive capability of the mapped site proxy.
In near- surface geophysics, ground-based mapping surveys are routinely used in a variety of applications including those from archaeology, civil engineering, hydrology, and soil science. The resulting geophysical anomaly maps of, for example, magnetic or electrical parameters are usually interpreted to laterally delineate subsurface structures such as those related to the remains of past human activities, subsurface utilities and other installations, hydrological properties, or different soil types. To ease the interpretation of such data sets, we have developed a multiscale processing, analysis, and visualization strategy. Our approach relies on a discrete redundant wavelet transform (RWT) implemented using cubic-spline filters and the a trous algorithm, which allows to efficiently compute a multiscale decomposition of 2D data using a series of 1D convolutions. The basic idea of the approach is presented using a synthetic test image, whereas our archaeogeophysical case study from northeast Germany demonstrates its potential to analyze and process rather typical geophysical anomaly maps including magnetic and topographic data. Our vertical-gradient magnetic data show amplitude variations over several orders of magnitude, complex anomaly patterns at various spatial scales, and typical noise patterns, whereas our topographic data show a distinct hill structure superimposed by a microtopographic stripe pattern and random noise. Our results demonstrate that the RWT approach is capable to successfully separate these components and that selected wavelet planes can be scaled and combined so that the reconstructed images allow for a detailed, multiscale structural interpretation also using integrated visualizations of magnetic and topographic data. Because our analysis approach is straightforward to implement without laborious parameter testing and tuning, computationally efficient, and easily adaptable to other geophysical data sets, we believe that it can help to rapidly analyze and interpret different geophysical mapping data collected to address a variety of near-surface applications from engineering practice and research.
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.
Projected changes in compound flood hazard from riverine and coastal floods in northwestern Europe
(2020)
Compound flooding in coastal regions, that is, the simultaneous or successive occurrence of high sea levels and high river flows, is expected to increase in a warmer world. To date, however, there is no robust evidence on projected changes in compound flooding for northwestern Europe. We combine projected storm surges and river floods with probabilistic, localized relative sea-level rise (SLR) scenarios to assess the future compound flood hazard over northwestern coastal Europe in the high (RCP8.5) emission scenario. We use high-resolution, dynamically downscaled regional climate models (RCM) to drive a storm surge model and a hydrological model, and analyze the joint occurrence of high coastal water levels and associated river peaks in a multivariate copula-based approach. The RCM-forced multimodel mean reasonably represents the observed spatial pattern of the dependence strength between annual maxima surge and peak river discharge, although substantial discrepancies exist between observed and simulated dependence strength. All models overestimate the dependence strength, possibly due to limitations in model parameterizations. This bias affects compound flood hazard estimates and requires further investigation. While our results suggest decreasing compound flood hazard over the majority of sites by 2050s (2040-2069) compared to the reference period (1985-2005), an increase in projected compound flood hazard is limited to around 34% of the sites. Further, we show the substantial role of SLR, a driver of compound floods, which has frequently been neglected. Our findings highlight the need to be aware of the limitations of the current generation of Earth system models in simulating coastal compound floods.
The 2.7-2.9 Ma Ertsberg East Skarn System (EESS) is a world-class Cu-Au skarn that formed within and adjacent to an intrusion within a paleodepth of 0.5 km and > 2.5 km. Its economic mineralisation developed by sustained reaction of magmatic fluid with contact metamorphosed siliciclastic and carbonate rocks at the margin of the adjacent Ertsberg quartz monzodiorite intrusion. Based on high-resolution mineral mapping, chemical analysis and thermodynamic calculations, the multistage formation processes of the exoskarn components of the EESS are examined in the context of changing pressure, temperature, fluid composition and fluid phase. We show that contact metamorphism of dolomitic sediments occurred at 51 +/- 5 MPa, between 700 degrees C and 800 degrees C and in the presence of a H2O-CO2-fluid containing similar to 10 to similar to 70 mol% CO2. This prograde metamorphism formed a forsterite + diopside + calcite + phlogopite + spinel assemblage. Such forsterite-dominated skarns account for similar to 55 vol% of the EESS exoskarns. Rare pargasite (previously unrecognized in this deposit) formed locally in the metamorphosed carbonate sequence where the protolith was composed of supratidal evaporites with dolomitic carbonate and interlayered calc-silicate rocks. The subsequent flux of a lower pressure magmatic gas containing SO2(g) caused sulphate metasomatism. This high temperature gas alteration of the metamorphic assemblage also caused skarn Cu-Fe-sulphide mineralisation. The influx of a SO2 gas through fracture permeability occurred at a temperature between similar to 600 and 700 degrees C and caused calcite to be replaced by anhydrite, with the coupled release of H2S(g). This in-situ release of H2S(g) scavenged trace Cu from the gas phase to deposit Cu-Fe-sulphides, which make the economic value of the distinct. We demonstrate that the formation of metal sulphides within forsterite skarns of the Ertsberg East Skarn System required a minimum flux of similar to 1,050 Mt SO2(g) and show that volcanic degassing may have occurred over a time span of similar to 3,900 years. As the system waned, the ambient fluid resulted in partial retrograde serpentinization of olivine and diopside without carbonation, and at temperatures sufficiently high to preserve anhydrite.
Flood loss modeling is a central component of flood risk analysis. Conventionally, this involves univariable and deterministic stage-damage functions. Recent advancements in the field promote the use of multivariable and probabilistic loss models, which consider variables beyond inundation depth and account for prediction uncertainty. Although companies contribute significantly to total loss figures, novel modeling approaches for companies are lacking. Scarce data and the heterogeneity among companies impede the development of company flood loss models. We present three multivariable flood loss models for companies from the manufacturing, commercial, financial, and service sector that intrinsically quantify prediction uncertainty. Based on object-level loss data (n = 1,306), we comparatively evaluate the predictive capacity of Bayesian networks, Bayesian regression, and random forest in relation to deterministic and probabilistic stage-damage functions, serving as benchmarks. The company loss data stem from four postevent surveys in Germany between 2002 and 2013 and include information on flood intensity, company characteristics, emergency response, private precaution, and resulting loss to building, equipment, and goods and stock. We find that the multivariable probabilistic models successfully identify and reproduce essential relationships of flood damage processes in the data. The assessment of model skill focuses on the precision of the probabilistic predictions and reveals that the candidate models outperform the stage-damage functions, while differences among the proposed models are negligible. Although the combination of multivariable and probabilistic loss estimation improves predictive accuracy over the entire data set, wide predictive distributions stress the necessity for the quantification of uncertainty.
The change of the mechanical properties of granular materials with pressure is an important topic associated with many industrial applications. In this paper we investigate the influence of hydrostatic pressure (P-e) on the effective bulk compressibility (C-eff) of a granular material by applying two modified theoretical approaches that are based on contact mechanics and micromechanics, respectively. For a granular material composed of rough grains, an extended contact model is developed to elucidate the effect of roughness of grain surfaces on bulk compressibility. At relatively low pressures, the model predicts that the decrease of bulk compressibility with pressure may be described by a power law with an exponent of -1/2 (i.e., C-eff proportional to P-e(1/2) ), but deviates at intermediate pressures. At elevated pressures beyond full contact, bulk compressibility remains almost unchanged, which may be roughly evaluated by continuum contact mechanics. As an alternative explanation of pressure-dependent bulk compressibility, we suggest a micromechanical model that accounts for effects of different types of pore space present in granular materials. Narrow and compliant inter-granular cracks are approximated by three-dimensional oblate spheroidal cracks with rough surfaces, whereas the equant and stiff pores surrounded by three and four neighboring grains are modeled as tubular pores with cross sections of three and four cusp-like corners, respectively. In this model, bulk compressibility is strongly reduced with increasing pressure by progressive closure of rough-walled cracks. At pressures exceeding crack closure pressure, deformation of the remaining equant pores is largely insensitive to pressure, with almost no further change in bulk compressibility. To validate these models, we performed hydrostatic compression tests on Bentheim sandstone (a granular rock consisting of quartz with high porosity) under a wide range of pressure. The relation between observed microstructures and measured pressure-dependent bulk compressibility is well explained by both suggested models.
The Arctic ecosystem, a reservoir of genetic microbial diversity, represents a virtually unlimited source of microorganisms that could interact with human beings. Despite continuous exploration of Arctic habitats and description of their microbial communities, bacterial phenotypes commonly associated with pathogenicity, such as hemolytic activity, have rarely been reported. In this study, samples of snow, fresh and marine water, soil, and sediment from several habitats in the Arctic archipelago of Svalbard were collected during Summer, 2017. Bacterial isolates were obtained after incubation on oligotrophic media at different temperatures and their hemolytic potential was assessed on sheep blood agar plates. Partial (alpha) or true (beta) hemolysis was observed in 32 out of 78 bacterial species. Genes expressing cytolytic compounds, such as hemolysins, likely increase the general fitness of the producing microorganisms and confer a competitive advantage over the availability of nutrients in natural habitats. In environmental species, the nutrient-acquisition function of these compounds presumably precedes their function as toxins for mammalian erythrocytes. However, in the light of global warming, the presence of hemolytic bacteria in Arctic environments highlights the possible risks associated with these microorganisms in the event of habitat melting/destruction, ecosystem transition, and re-colonization.
Many of the volcanic plateau margins of the eastern, formerly glaciated, foreland of the Patagonian Andes are undermined by giant landslides (>= 10(8) m(3)). One cluster of such landslides extends along the margin of the Meseta del Lago Buenos Aires (MLBA) plateau that is formed mainly by Neogene-Quaternary basalts. The dry climate is at odds with numerous >2-km long earthflows nested within older and larger compound landslides. We present a hydrological analysis, a detailed geomorphic map, interpretations of exposed landslide interiors, and radiocarbon dating of the El Mirador landslide, which is one of the largest and morphologically most representative landslide. We find that the presence of lakes on top of the plateau, causing low infiltration rates, correlates negatively with the abundance of earthflows on compound landslides along the plateau margins. Field outcrops show that the pattern of compound landslides and earthflows is likely controlled by groundwater seepage at the contact between the basalts and underlying soft Miocene molasse. Numerous peat bogs store water and sediment and are more abundant in earthflow-affected areas than in their contributing catchment areas. <br /> Radiocarbon dates indicate that these earthflows displaced metre-thick layers of peat in the late Holocene (<2.5 ka). We conclude that earthflows of the MLBA plateau might be promising proxies of past hydroclimatic conditions in the Patagonian foreland, if strong earthquakes or gradual crustal stress changes due to glacioisostatic rebound can be ruled out.
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.
The Frasnian-Famennian (F-F) transition of Late Devonian was a critical episode in geological history, recording a major mass extinction event. In this study, we focus on an F-F succession from a deep marine context in Bancheng, southern Guangxi, South China, to investigate coeval changes in pelagic environments of the Paleo-Tethys Ocean. The studied succession is exclusively composed of bedded cherts intercalated with multiple siliceous volcanic ash beds. A SIMS zircon U-Pb Concordia age of 367.8 +/- 2.5 Ma is reported for a tuffaceous layer slightly above the F-F boundary. Geochemical ratios of Al/(Al + Fe + Mn), Ce/Ce*, Y/Ho, and Al, Fe contents in bedded cherts indicate that they are of predominantly biogenic/chemical origin with some terrigenous inputs. Negligible enrichment of redox sensitive elements (Mo, U, V) and low V/Cr ratios (<2) suggest persistently oxic conditions existed in the deep pelagic basin at Bancheng, South China during the F-F transition. These findings call into question the widely held hypothesis that marine anoxia was the primary killing mechanism for the F-F crisis. In contrast, multiple tuffaceous layers throughout the F-F boundary succession indicate frequent volcanic activity, which could have released massive amounts of greenhouse gases into the atmosphere, inducing climate warming. This scenario may have increased continental weathering and riverine fluxes into the ocean, reconciling the increases in Al2O3 content and Al/(Al + Fe + Mn) ratio across the F-F boundary. Documentation of persistently oxic conditions and frequent volcanic activitiy provides new perspectives on the inter-relationship between volcanism, climate, and oceanic redox fluctuation during the F-F biotic crisis.
Temperature changes and variations in pore fluid salinity may negatively affect the permeability of clay-bearing sandstones with implications for natural fluid flow and geotechnical applications alike. In this study these factors are investigated for a sandstone dominated by illite as the clay phase. Systematic long-term flow-through experiments were conducted and complemented with comprehensive microstructural investigations and the application of Derjaguin-Landau-Verwey-Overbeek (DLVO) theory to explain mechanistically the observed permeability changes. Initially, sample permeability was not affected by low pore fluid salinity indicating strong attraction of the illite particles to the pore walls as supported by electron microprobe analysis (EMPA). Increasing temperature up to 145 degrees C resulted in an irreversible permeability decrease by 1.5 orders of magnitude regardless of the pore fluid composition (i.e., deionized water and 2 M NaCl solution). Subsequently diluting the high salinity pore fluid to below 0.5 M yielded an additional permeability decline by 1.5 orders of magnitude, both at 145 degrees C and after cooling to room temperature. By applying scanning electron microscopy (SEM) and mercury intrusion porosimetry (MIP) thermo-mechanical pore throat closure and illite particle migration were identified as independently operating mechanisms responsible for observed permeability changes during heating and dilution, respectively. These observations indicate that permeability of illite-bearing sandstones will be impaired by heating and exposure to low salinity pore fluids. In addition, chemically induced permeability variations proved to be path dependent with respect to the applied succession of fluid salinity changes.
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 transfer of particulate organic carbon from continents to the ocean is an important component of the global carbon cycle. Transfer to and burial of photosynthetically fixed biospheric organic carbon in marine sediments can effectively sequester atmospheric carbon dioxide over geological timescales. The exhumation and erosion of fossil organic carbon contained in sedimentary rocks, i.e. petrogenic carbon, can result in remineralization, releasing carbon to the atmosphere. In contrast, eroded petrogenic organic carbon that gets transferred back to the ocean and reburied does not affect atmospheric carbon content.
Mountain ranges play a key role in this transfer since they can source vast amounts of sediment including particulate organic carbon. Globally, the export of both, biospheric and petrogenic organic carbon has been linked to sediment export. Additionally, short transfer times from mountains to the ocean and high sediment concentrations have been shown to increase the likelihood of organic carbon burial. While the importance of mountain ranges in the organic carbon cycle is now widely recognized, the processes acting within mountain ranges to influence the storage, cycling and mobilization of organic carbon, as well as carbon fluxes from mountain ranges remain poorly constrained.
In this thesis, I employ different methods to assess the nature and fate of particulate organic carbon in mountain belts, ranging from the molecular to regional landscape scale. These studies are located along the Trans-Himalayan Kali Gandaki River in Central Nepal. This river traverses all major geological and climatic zones of the Himalaya, from the dry northern Tibetan plateau to the high-relief, monsoon dominated steep High Himalaya and the lower relief and abundant vegetation of the Lesser Himalayan region.
First, I document how biospheric organic matter has accumulated during the Holocene in the headwaters of the Kali Gandaki River valley, by combining compound specific isotope measurements with different dating methods and grain size data, and investigate the stability of this organic carbon reservoir on millennial timescales. I show, that around 1.6 ka an eco-geomorphic tipping point occurred leading to a destabilization of the landscape resulting in today’s high erosion rates and the excavation of the aged organic carbon reservoir. This study highlights the climatic and geomorphic controls on biospheric organic carbon storage and release from mountain ranges.
Second, I systematically investigate the spatial variation of particulate organic carbon fluxes across the Himalaya along the Kali Gandaki River, using bulk stable and radioactive isotopes combined with a new Bayesian modeling approach. The detailed dataset allows the distinction of aged and modern biospheric organic carbon as well as petrogenic organic carbon across the Himalayan mountain range and the investigation of the role of climatic and geomorphic factors in their riverine export. The data suggest a decoupling of the particulate organic carbon from the sediment yield along the Kali Gandaki River, partially driven by climatic and geomorphic processes. In contrast to the suspended sediment, a large part of the particulate organic carbon exported by the river originates from the Tibetan part of the catchment and is dominated by petrogenic organic carbon derived from Jurassic shales with only minor contributions of modern and aged biospheric organic carbon. These findings emphasize the importance of organic carbon source distribution and erosion mechanisms in determining the organic carbon export from mountain ranges.
In a third step, I explore the potential of ultra-high resolution mass spectrometry for particulate organic carbon transport studies. I have generated a novel and unprecedented high-resolution molecular dataset, which contains up to 103 molecular formulas of the lipid fraction of particulate organic matter for modern and aged biospheric carbon, petrogenic organic carbon and river sediments. First, I test if this dataset can be used to better resolve different organic carbon sources and to identify new geochemical tracers. Using multivariate statistics, I identify up to 10² characteristic molecular formulas for the major organic carbon sources in the upper part of the Kali Gandaki catchment, and trace their transfer from the surrounding landscape into the river sediment. Second, I test the potential of the molecular dataset to trace molecular transformations along source-to-sink pathways. I identify changes in molecular metrics derived from the dataset, which are characteristic of transformation processes during incorporation of litter into soil, the aging of soil material, and the mobilization of the organic carbon into the river. These two studies demonstrate that high-resolution molecular datasets open a promising analytical window on particulate organic carbon and can provide novel insights into the composition, sourcing and transformation of riverine particulate organic carbon.
Collectively, these studies advance our understanding of the processes contributing to the storage and mobilization of organic carbon in the Central Himalaya, the mountain belt that dominates global erosional fluxes. They do so by identifying the major sources of particulate organic carbon to the Trans-Himalayan Kali Gandaki River, by elucidating their sensitivity to climate and geomorphic processes, and by identifying some of the transformations of this material on the molecular scale. As a result, the thesis demonstrates that the amount and composition of organic carbon routed from mountain belts is a function of the dynamic interactions of geologic, biologic, geomorphic and climatic processes within the mountain belt. This understanding will ultimately help in answering whether the build-up and erosion of mountain ranges over geological time represents a net carbon source or sink to the atmosphere. Beyond this, the thesis contributes to our technical ability to characterize organic matter and attribute it to sources by scoping the potential of high-end molecular analysis.
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.
Una extraordinaria edición actualizada de Observaciones sobre el clima de Lima y sus influencias en los seres organizados, en especial el hombre de José Hipólito Unanue y Pavón ha sido publicado en el 2018 y merece ser anunciado a la comunidad académica peruana e internacional. Este libro incluye un exhaustivo estudio in troductorio del destacado historiador Lizardo Seiner Lizárraga. Esta introducción está conformada de tres partes encabezadas por los siguientes títulos: a) Unanue: una historiografía de amplio espectro b) Textos, registros y referencias: una indagación en la biblioteca de Unanue y c) Observaciones sobre el clima de Lima: una obra en dos tiempos.
Over the past few decades, azimuthal seismic anisotropy measurements have been widely used proxy to study past and present-day deformation of the lithosphere and to characterize convection in the mantle. Beneath continental regions, distinguishing between shallow and deep sources of anisotropy remains difficult due to poor depth constraints of measurements and a lack of regional-scale geodynamic modeling. Here, we constrain the sources of seismic anisotropy beneath Madagascar where a complex pattern cannot be explained by a single process such as absolute plate motion, global mantle flow, or geology. We test the hypotheses that either Edge-Driven Convection (EDC) or mantle flow derived from mantle wind interactions with lithospheric topography is the dominant source of anisotropy beneath Madagascar. We, therefore, simulate two sets of mantle convection models using regional-scale 3-D computational modeling. We then calculate Lattice Preferred Orientation that develops along pathlines of the mantle flow models and use them to calculate synthetic splitting parameters. Comparison of predicted with observed seismic anisotropy shows a good fit in northern and southern Madagascar for the EDC model, but the mantle wind case only fits well in northern Madagascar. This result suggests the dominant control of the measured anisotropy may be from EDC, but the role of localized fossil anisotropy in narrow shear zones cannot be ruled out in southern Madagascar. Our results suggest that the asthenosphere beneath northern and southern Madagascar is dominated by dislocation creep. Dislocation creep rheology may be dominant in the upper asthenosphere beneath other regions of continental lithosphere.
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.
Multiple P-T-d-t paths reveal the evolution of the final Nuna assembly in northeast Australia
(2020)
The final assembly of the Mesoproterozoic supercontinent Nuna was marked by the collision of Laurentia and Australia at 1.60 Ga, which is recorded in the Georgetown Inlier of NE Australia. Here, we decipher the metamorphic evolution of this final Nuna collisional event using petrostructural analysis, major and trace element compositions of key minerals, thermodynamic modelling, and multi-method geochronology. The Georgetown Inlier is characterised by deformed and metamorphosed 1.70-1.62 Ga sedimentary and mafic rocks, which were intruded byc. 1.56 Ga old S-type granites. Garnet Lu-Hf and monazite U-Pb isotopic analyses distinguish two major metamorphic events (M1 atc. 1.60 Ga and M2 atc. 1.55 Ga), which allows at least two composite fabrics to be identified at the regional scale-c. 1.60 Ga S1 (consisting in fabrics S1a and S1b) andc. 1.55 Ga S2 (including fabrics S2a and S2b). Also, three tectono-metamorphic domains are distinguished: (a) the western domain, with S1 defined by low-P(LP) greenschist facies assemblages; (b) the central domain, where S1 fabric is preserved as medium-P(MP) amphibolite facies relicts, and locally as inclusion trails in garnet wrapped by the regionally dominant low-Pamphibolite facies S2 fabric; and (c) the eastern domain dominated by upper amphibolite to granulite facies S2 foliation. In the central domain, 1.60 GaMP-medium-T(MT) metamorphism (M1) developed within the staurolite-garnet stability field, with conditions ranging from 530-550 degrees C at 6-7 kbar (garnet cores) to 620-650 degrees C at 8-9 kbar (garnet rims), and it is associated with S1 fabric. The onset of 1.55 GaLP-high-T(HT) metamorphism (M2) is marked by replacement of staurolite by andalusite (M2a/D2a), which was subsequently pseudomorphed by sillimanite (M2b/D2b) where granite and migmatite are abundant.P-Tconditions ranged from 600 to 680 degrees C and 4-6 kbar for the M2b sillimanite stage. 1.60 Ga garnet relicts within the S2 foliation highlight the progressive obliteration of the S1 fabric by regional S2 in the central zone during peak M2 metamorphism. In the eastern migmatitic complex, partial melting of paragneiss and amphibolite occurred syn- to post-S2, at 730-770 degrees C and 6-8 kbar, and at 750-790 degrees C and 6 kbar, respectively. The pressure-temperature-deformation-time paths reconstructed for the Georgetown Inlier suggest ac. 1.60 Ga M1/D1 event recorded under greenschist facies conditions in the western domain and under medium-Pand medium-Tconditions in the central domain. This event was followed by the regional 1.56-1.54 Ga low-Pand high-Tphase (M2/D2), extensively recorded in the central and eastern domains. Decompression between these two metamorphic events is ascribed to an episode of exhumation. The two-stage evolution supports the previous hypothesis that the Georgetown Inlier preserves continental collisional and subsequent thermal perturbation associated with granite emplacement.
Multi-year index-based insurance for adapting Water Utility Companies to hydrological drought
(2020)
The sustainability of water utility companies is threatened by non-stationary drivers, such as climate and anthropogenic changes. To cope with potential economic losses, instruments such as insurance are useful for planning scenarios and mitigating impacts, but data limitations and risk uncertainties affect premium estimation and, consequently, business sustainability. This research estimated the possible economic impacts of business interruption to the Sao Paulo Water Utility Company derived from hydrological drought and how this could be mitigated with an insurance scheme. Multi-year insurance (MYI) was proposed through a set of "change" drivers: the climate driver, through forcing the water evaluation and planning system (WEAP) hydrological tool; the anthropogenic driver, through water demand projections; and the economic driver, associated with recent water price policies adopted by the utility company during water scarcity periods. In our study case, the evaluated indices showed that MYI contracts that cover only longer droughts, regardless of the magnitude, offer better financial performance than contracts that cover all events (in terms of drought duration). Moreover, through MYI contracts, we demonstrate solvency for the insurance fund in the long term and an annual average actuarially fair premium close to the total expected revenue reduction.
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
The Hartousov mofette system is a natural CO2 degassing site in the central Cheb Basin (Eger Rift, Central Europe). In early 2016 a 108 m deep core was obtained from this system to investigate the impact of ascending mantle-derived CO2 on indigenous deep microbial communities and their surrounding life habitat. During drilling, a CO2 blow out occurred at a depth of 78.5 meter below surface (mbs) suggesting a CO2 reservoir associated with a deep low-permeable CO2-saturated saline aquifer at the transition from Early Miocene terrestrial to lacustrine sediments. Past microbial communities were investigated by hopanoids and glycerol dialkyl glycerol tetraethers (GDGTs) reflecting the environmental conditions during the time of deposition rather than showing a signal of the current deep biosphere. The composition and distribution of the deep microbial community potentially stimulated by the upward migration of CO2 starting during Mid Pleistocene time was investigated by intact polar lipids (IPLs), quantitative polymerase chain reaction (qPCR), and deoxyribonucleic acid (DNA) analysis. The deep biosphere is characterized by microorganisms that are linked to the distribution and migration of the ascending CO2-saturated groundwater and the availability of organic matter instead of being linked to single lithological units of the investigated rock profile. Our findings revealed high relative abundances of common soil and water bacteria, in particular the facultative, anaerobic and potential iron-oxidizing Acidovorax and other members of the family Comamonadaceae across the whole recovered core. The results also highlighted the frequent detection of the putative sulfate-oxidizing and CO2-fixating genus Sulfuricurvum at certain depths. A set of new IPLs are suggested to be indicative for microorganisms associated to CO2 accumulation in the mofette system.
This thesis investigates how the permafrost microbiota responds to global warming. In detail, the constraints behind methane production in thawing permafrost were linked to methanogenic activity, abundance and composition. Furthermore, this thesis offers new insights into microbial adaptions to the changing environmental conditions during global warming. This was assesed by investigating the potential ecological relevant functions encoded by plasmid DNA within the permafrost microbiota. Permafrost of both interglacial and glacial origin spanning the Holocene to the late Pleistocene, including Eemian, were studied during long-term thaw incubations. Furthermore, several permafrost cores of different stratigraphy, soil type and vegetation cover were used to target the main constraints behind methane production during short-term thaw simulations. Short- and long-term incubations simulating thaw with and without the addition of substrate were combined with activity measurements, amplicon and metagenomic sequencing of permanently frozen and seasonally thawed active layer. Combined, it allowed to address the following questions. i) What constraints methane production when permafrost thaws and how is this linked to methanogenic activity, abundance and composition? ii) How does the methanogenic community composition change during long-term thawing conditions? iii) Which potential ecological relevant functions are encoded by plasmid DNA in active layer soils?
The major outcomes of this thesis are as follows. i) Methane production from permafrost after long-term thaw simulation was found to be constrained mainly by the abundance of methanogens and the archaeal community composition. Deposits formed during periods of warmer temperatures and increased precipitation, (here represented by deposits from the Late Pleistocene of both interstadial and interglacial periods) were found to respond strongest to thawing conditions and to contain an archaeal community dominated by methanogenic archaea (40% and 100% of all detected archaea). Methanogenic population size and carbon density were identified as main predictors for potential methane production in thawing permafrost in short-term incubations when substrate was sufficiently available.
ii) Besides determining the methanogenic activity after long-term thaw, the paleoenvironmental conditions were also found to influence the response of the methanogenic community composition. Substantial shifts within methanogenic community structure and a drop in diversity were observed in deposits formed during warmer periods, but not in deposits from stadials, when colder and drier conditions occurred. Overall, a shift towards a dominance of hydrogenotrophic methanogens was observed in all samples, except for the oldest interglacial deposits from the Eemian, which displayed a potential dominance of acetoclastic methanogens. The Eemian, which is discussed to serve as an analogue to current climate conditions, contained highly active methanogenic communities. However, all potential limitation of methane production after permafrost thaw, it means methanogenic community structure, methanogenic population size, and substrate pool might be overcome after permafrost had thawed on the long-term. iii) Enrichments with soil from the seasonally thawed active layer revealed that its plasmid DNA (‘metaplasmidome’) carries stress-response genes. In particular it encoded antibiotic resistance genes, heavy metal resistance genes, cold shock proteins and genes encoding UV-protection. Those are functions that are directly involved in the adaptation of microbial communities to stresses in polar environments. It was further found that metaplasmidomes from the Siberian active layer originate mainly from Gammaproteobacteria. By applying enrichment cultures followed by plasmid DNA extraction it was possible to obtain a higher average contigs length and significantly higher recovery of plasmid sequences than from extracting plasmid sequences from metagenomes. The approach of analyzing ‘metaplasmidomes’ established in this thesis is therefore suitable for studying the ecological role of plasmids in polar environments in general.
This thesis emphasizes that including microbial community dynamics have the potential to improve permafrost-carbon projections. Microbially mediated methane release from permafrost environments may significantly impact future climate change. This thesis identified drivers of methanogenic composition, abundance and activity in thawing permafrost landscapes. Finally, this thesis underlines the importance to study how the current warming Arctic affects microbial communities in order to gain more insight into microbial response and adaptation strategies.
Meteorological and hydrological drought assessment in Lake Malawi and Shire River basins (1970-2013)
(2020)
The study assesses the variability and trends of both meteorological and hydrological droughts from 1970 to 2013 in Lake Malawi and Shire River basins using the standardized precipitation index (SPI) and standardized precipitation and evaporation index (SPEI) for meteorological droughts and the lake level change index (LLCI) for hydrological droughts. Trends and slopes in droughts and drought drivers are estimated using Mann-Kendall test and Sen's slope, respectively. Results suggest that meteorological droughts are increasing due to a decrease in precipitation which is exacerbated by an increase in temperature (potential evapotranspiration). The hydrological system of Lake Malawi seems to have a >24-month memory towards meteorological conditions, since the 36-month SPEI can predict hydrological droughts 10 months in advance. The study has found the critical lake level that would trigger hydrological drought to be 474.1 m a.s.l. The increase in drought is a concern as this will have serious impacts on water resources and hydropower supply in Malawi.
The formation of the supercontinent Pangaea during the Permo-Triassic gave rise to an extreme monsoonal climate (often termed 'mega-monsoon') that has been documented by numerous palaeo-records. However, considerable debate exists about the role of orbital forcing in causing humid intervals in an otherwise arid climate. To shed new light on the forcing of monsoonal variability in subtropical Pangaea, this study focuses on sediment facies and colour variability of playa and alluvial fan deposits in an outcrop from the late Carnian (ca 225 Ma) in the southern Germanic Basin, south-western Germany. The sediments were deposited against a background of increasingly arid conditions following the humid Carnian Pluvial Event (ca 234 to 232 Ma). The ca 2 center dot 4 Myr long sedimentary succession studied shows a tripartite long-term evolution, starting with a distal mud-flat facies deposited under arid conditions. This phase was followed by a highly variable playa-lake environment that documents more humid conditions and finally a regression of the playa-lake due to a return of arid conditions. The red-green (a*) and lightness (L*) records show that this long-term variability was overprinted by alternating wet/dry cycles driven by orbital precession and ca 405 kyr eccentricity, without significant influence of obliquity. The absence of obliquity in this record indicates that high-latitude forcing played only a minor role in the southern Germanic Basin during the late Carnian. This is different from the subsequent Norian when high-latitude signals became more pronounced, potentially related to the northward drift of the Germanic Basin. The recurring pattern of pluvial events during the late Triassic demonstrates that orbital forcing, in particular eccentricity, stimulated the occurrence and intensity of wet phases. It also highlights the possibility that the Carnian Pluvial Event, although most likely triggered by enhanced volcanic activity, may also have been modified by an orbital stimulus.
Subsea permafrost is perennially cryotic earth material that lies offshore. Most submarine permafrost is relict terrestrial permafrost beneath the Arctic shelf seas, was inundated after the last glaciation, and has been warming and thawing ever since. It is a reservoir and confining layer for gas hydrates and has the potential to release greenhouse gases and affect global climate change. Furthermore, subsea permafrost thaw destabilizes coastal infrastructure. While numerous studies focus on its distribution and rate of thaw over glacial timescales, these studies have not been brought together and examined in their entirety to assess rates of thaw beneath the Arctic Ocean. In addition, there is still a large gap in our understanding of sub-aquatic permafrost processes on finer spatial and temporal scales. The degradation rate of subsea permafrost is influenced by the initial conditions upon submergence. Terrestrial permafrost that has already undergone warming, partial thawing or loss of ground ice may react differently to inundation by seawater compared to previously undisturbed ice-rich permafrost. Heat conduction models are sufficient to model the thaw of thick subsea permafrost from the bottom, but few studies have included salt diffusion for top-down chemical degradation in shallow waters characterized by mean annual cryotic conditions on the seabed. Simulating salt transport is critical for assessing degradation rates for recently inundated permafrost, which may accelerate in response to warming shelf waters, a lengthening open water season, and faster coastal erosion rates. In the nearshore zone, degradation rates are also controlled by seasonal processes like bedfast ice, brine injection, seasonal freezing under floating ice conditions and warm freshwater discharge from large rivers. The interplay of all these variables is complex and needs further research. To fill this knowledge gap, this thesis investigates sub-aquatic permafrost along the southern coast of the Bykovsky Peninsula in eastern Siberia. Sediment cores and ground temperature profiles were collected at a freshwater thermokarst lake and two thermokarst lagoons in 2017. At this site, the coastline is retreating, and seawater is inundating various types of permafrost: sections of ice-rich Pleistocene permafrost (Yedoma) cliffs at the coastline alternate with lagoons and lower elevation previously thawed and refrozen permafrost basins (Alases). Electrical resistivity surveys with floating electrodes were carried out to map ice-bearing permafrost and taliks (unfrozen zones in the permafrost, usually formed beneath lakes) along the diverse coastline and in the lagoons. Combined with the borehole data, the electrical resistivity results permit estimation of contemporary ice-bearing permafrost characteristics, distribution, and occasionally, thickness. To conceptualize possible geomorphological and marine evolutionary pathways to the formation of the observed layering, numerical models were applied. The developed model incorporates salt diffusion and seasonal dynamics at the seabed, including bedfast ice. Even along coastlines with mean annual non-cryotic boundary conditions like the Bykovsky Peninsula, the modelling results show that salt diffusion minimizes seasonal freezing of the seabed, leading to faster degradation rates compared to models without salt diffusion. Seasonal processes are also important for thermokarst lake to lagoon transitions because lagoons can generate cold hypersaline conditions underneath the ice cover. My research suggests that ice-bearing permafrost can form in a coastal lagoon environment, even under floating ice. Alas basins, however, may degrade more than twice as fast as Yedoma permafrost in the first several decades of inundation. In addition to a lower ice content compared to Yedoma permafrost, Alas basins may be pre-conditioned with salt from adjacent lagoons. Considering the widespread distribution of thermokarst in the Arctic, its integration into geophysical models and offshore surveys is important to quantify and understand subsea permafrost degradation and aggradation. Through numerical modelling, fieldwork, and a circum-Arctic review of subsea permafrost literature, this thesis provides new insights into sub-aquatic permafrost evolution in saline coastal environments.
A mid-aperture broad-band test array (OBS array DOCTAR) was deployed from June 2011 to April 2012 about 100 km north of the Gloria fault in the Eastern North Atlantic in about 5000 m water depth. In addition arrays were installed on Madeira Island and in western Portugal mainland. For the first time in the Eastern North Atlantic, we recorded a large number of high frequency Po and So waves from local and regional small and moderate earthquakes (M-L < 4). An incoherent beamforming method was adapted to scan continuous data for such Po and So arrivals applying a sliding window waveform migration and frequency-wavenumber technique. We identify about 320 Po and 1550 So arrivals and compare the phase onsets with the ISC catalogue (ISC 2015) for the same time span. Up to a distance of 6 degrees to the DOCTAR stations all events listed in the ISC catalogue could be associated to Po and So phases. Arrivals from events in more than 10 degrees distance could be identified only in some cases. Only few Po and/or So arrivals were detected for earthquakes from the European and African continental area, the continental shelf regions and for earthquakes within or northwest of the Azores plateau. Unexpectedly, earthquake clusters are detected within the oceanic plates north and south of the Gloria fault and far from plate boundaries, indicating active intraplate structures. We also observe and locate numerous small magnitude earthquakes on the segment of the Gloria fault directly south of DOCTAR, which likely coincides with the rupture of the 25 November 1941 event. Local small magnitude earthquakes located beneath DOCTAR show hypocentres up to 30 km depth and strike-slip focal mechanisms. A comparison with detections at temporary mid-aperture arrays on Madeira and in western Portugal shows that the deep ocean array performs much better than the island and the continental array regarding the detection threshold for events in the oceanic plates.
We conclude that sparsely distributed mid-aperture seismic arrays in the deep ocean could decrease the detection and location threshold for seismicity with M-L < 4 in the oceanic plate and might constitute a valuable tool to monitor oceanic plate seismicity.
Reported magnitudes ranging between 7.8 and 8.7 highlight a confusion about the true size of the 1920 Haiyuan earthquake, the largest earthquake recorded in the northeast Tibetan Plateau. We compiled a global data set of previously unlooked-at historical seismograms and performed modern computational analyses on the digitized seismic records. We found the epicenter to be near Haiyuan town and obtained a moment magnitude of M-W=7.90.2. Following traditional approaches, we obtained m(B)=7.90.3 with data from 21 stations and M-S(20)=8.10.2 with data from three stations. Geomorphologically, we mapped the surface rupture and horizontal offsets on high-resolution Pleiades satellite and drone imagery that covered the entire western and middle sections of the 1920 Haiyuan earthquake rupture and compiled offsets reported on the eastern section from field measurements in the 1980s. Careful discrimination between single-event and cumulative offsets suggests average horizontal slips of 3.01.0m on the western section, 4.51.5m on the middle section, and 3.5 +/- 0.5m on the eastern section, indicating a total moment magnitude of M-W=7.8 +/- 0.1. Thus, the seismological and geomorphological results agree within the uncertainties, a weighted average giving a moment magnitude of M-W=7.9 +/- 0.2 for the 1920 Haiyuan earthquake. It is likely that earthquake magnitudes based on the historical M were systematically overestimated. <br /> Plain Language Summary Earthquakes are the main mechanism by which elastic energy accumulating due to tectonic motion is released. As the earthquake magnitude scale is logarithmic, major earthquakes control the bulk of this energy budget and are often the most destructive, like the 1920 Haiyuan earthquake with similar to 230,000 casualties. However, major earthquakes tend to have recurrence periods of several hundred years, longer than our instrumental records. To obtain knowledge of historic major earthquakes, paleoseismologists measure geomorphic offsets and map surface ruptures left by past events and estimate the shaking intensity from historical writings. However, in the case of the Haiyuan earthquake, which happened in the late historic, early instrumental period, the magnitudes reported from these two communities differed significantly. In order to constrain the magnitude of this earthquake for seismic hazard assessment and to reconcile the differences between published magnitudes, we reestimated its magnitude from both newly compiled and digitized seismological records and modern satellite and drone imagery. The results show that the early seismological magnitudes were overestimated, which may affect historical magnitudes systematically. The 1920 Haiyuan earthquake was of a similar magnitude to the 2001 Kokoxili and 2008 Wenchuan earthquakes that also occurred in and around the Tibetan Plateau, instead of more than half a magnitude larger.
The occurrence of refugia beyond the arctic treeline and genetic adaptation therein play a crucial role of largely unknown effect size. While refugia have potential for rapidly colonizing the tundra under global warming, the taxa may be maladapted to the new environmental conditions. Understanding the genetic composition and age of refugia is thus crucial for predicting any migration response.
Here, we genotype 194 larch individuals from an similar to 1.8 km(2)area in northcentral Siberia on the southern Taimyr Peninsula by applying an assay of 16 nuclear microsatellite markers. For estimating the age of clonal individuals, we counted tree rings at sections along branches to establish a lateral growth rate that was then combined with geographic distance.
Findings reveal that the predominant reproduction type is clonal (58.76%) by short distance spreading of ramets. One outlier of clones 1 km apart could have been dispersed by reindeer. In clonal groups and within individuals, we find that somatic mutations accumulate with geographic distance. Clonal groups of two or more individuals are observed. Clonal age estimates regularly suggest individuals as old as 2,200 years, which coincides with a major environmental change that forced a treeline retreat in the region.
We conclude that individuals with clonal growth mode were naturally selected as it lowers the likely risk of extinction under a harsh environment. We discuss this legacy from the past that might now be a maladaptation and hinder expansion under currently strongly increasing temperatures.
Receiver function approaches have proven to be valuable for the investigation of crustal and upper mantle discontinuities whose sharp changes in seismic velocities cause wave conversions. While the crustal and mantle transition zone discontinuities are largely understood, the X-discontinuity at 250-350 km depth is still an object of controversial debate. The origin and global distribution of this structure with a velocity jump of 1.5-4.8% for compressional and shear waves is still unexplained. Although the crustal and mantle transition zone discontinuities beneath SW Morocco and surroundings have been investigated, only a few studies observed the X-discontinuity and place the depth at 260-370 km beneath the region of western Morocco. In order to better locate and characterize the X-discontinuity beneath southwest Morocco, we create P-wave receiver functions using data recorded by the Morocco-Munster array and detect the X-discontinuity at apparent depths of 285-350 km. In the western part of our study region we find apparent depths of similar to 310-340 km. The eastern part of the study area appears more complex: we locate two velocity jumps at apparent depths of around 285-295 km and 330-350 km in the northeast, and in the southeast we find a discontinuity at apparent depths of 340-350 km. Due to the large depth range and the twofold appearance of the X-discontinuity, we suggest that two different phase transitions cause the X-discontinuity beneath SW Morocco. The velocity contrasts at larger depths likely point to the coesite-stishovite phase transition occurring in deep eclogitic pools. The shallower depths can be explained by the transition from orthoenstatite to high-pressure clinoenstatite which requires the reaction between eclogite and peridotite to form orthopyroxene-rich peridotite. This reaction is likely related to previously proposed small-scale mantle upwellings beneath SW Morocco. Since both phase transitions require eclogite occurrence, the location of the X-discontinuity in this region can be used to indicate the location of recycled oceanic crust.
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.
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 (similar to km) relative displacements with an accuracy of <1 cm/epoch and similar to 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 thick-skinned fold-and-thrust belt on the eastern flank of the Andean Plateau in northwestern Argentina (NWA) is a zone of active contractional deformation characterized by fault-bounded mountain ranges with no systematic spatiotemporal pattern of tectonic activity. In contrast, the thin-skinned Subandean fold-and-thrust belt of northern Argentina and southern Bolivia is characterized primarily by in-sequence (i.e., west to east) fault progression, with a narrow zone of Quaternary deformation focused at the front of the orogenic wedge. To better understand how recent deformation is accommodated across these mountain ranges and the Argentinian portion of the orogen in particular, estimating and comparing deformation rates and patterns across different timescales is essential. We present Late Pleistocene shortening rates for the central Calchaqui intermontane valley in NWA associated with at least three episodes of deformation. Global Positioning System data for the same region reveal a gradual decrease in horizontal surface velocities from the Eastern Cordillera toward the foreland, which contrasts with the rapid velocity gradient associated with a locked decollement in the Subandean Ranges of southern Bolivia. Our new results represent a small view of regional deformation that, when considered in combination with the shallow crustal seismicity and decadal-scale surface velocities, support the notion that strain release in NWA is associated with numerous slowly deforming structures that are distributed throughout the orogen.
A silica undersaturated alkali-olivine basanitic magma intruded the late Paleocene/early Eocene Jafnayn Formation near Muscat. Geochemical analyses indicate that a significant amount of host rock (limestone) was assimilated into the magma. We dated the basanite as 42.7 +/- 1.0 Ma (2 sigma error; late Lutetian), using the whole rock Ar-40/Ar-39 step-wise heating technique. Intrusion occurred in the hanging wall of a major regional extensional shear zone (Frontal Range Fault, FRF) bounding the northern margin of two domes within the Oman Mountains (Jabal Akhdar and Saih Hatat domes). Two shear intervals along the FRF have been documented. The first interval lasted immediately after emplacement of the Semail Ophiolite (latest Cretaceous-early Eocene) while the second and poorly constrained interval was assumed to have occurred during the Oligocene.
The proximity of the basanite to the FRF suggests that magma used extensional faults for the upper part of its ascent path. Reactivated Permian rift faults of the Pangaea rift or other preexisting faults may have been used for the lower ascent part.
We conclude that the basanite intrusion coincided with the onset of the second deformation interval along the FRF, because (1) the position of the basanite is near a dextral releasing bend, associated with the second shear interval, (2) the overlap of our Ar-40/Ar-39 age with the cooling curves for rocks from the nearby Jabal Akhdar Dome, and (3) the basanite postdates the first FRF deformation episode by > 10 Ma. Thus, the second interval along the FRF had started already during the late Lutetian and probably lasted into the Miocene.
The Coulomb failure stress (CFS) criterion is the most commonly used method for predicting spatial distributions of aftershocks following large earthquakes. However, large uncertainties are always associated with the calculation of Coulomb stress change. The uncertainties mainly arise due to nonunique slip inversions and unknown receiver faults; especially for the latter, results are highly dependent on the choice of the assumed receiver mechanism. Based on binary tests (aftershocks yes/no), recent studies suggest that alternative stress quantities, a distance-slip probabilistic model as well as deep neural network (DNN) approaches, all are superior to CFS with predefined receiver mechanism. To challenge this conclusion, which might have large implications, we use 289 slip inversions from SRCMOD database to calculate more realistic CFS values for a layered half-space and variable receiver mechanisms. We also analyze the effect of the magnitude cutoff, grid size variation, and aftershock duration to verify the use of receiver operating characteristic (ROC) analysis for the ranking of stress metrics. The observations suggest that introducing a layered half-space does not improve the stress maps and ROC curves. However, results significantly improve for larger aftershocks and shorter time periods but without changing the ranking. We also go beyond binary testing and apply alternative statistics to test the ability to estimate aftershock numbers, which confirm that simple stress metrics perform better than the classic Coulomb failure stress calculations and are also better than the distance-slip probabilistic model.
Thirteen N-butylpyridinium salts, including three monometallic [C4Py](2)[MCl4], nine bimetallic [C4Py](2)[(M1-xMxCl4)-M-a-Cl-b] and one trimetallic compound [C4Py](2)[(M1-y-zMyMz (c) Cl4)-M-a-M-b] (M=Co, Cu, Mn; x=0.25, 0.50 or 0.75 and y=z=0.33), were synthesized and their structure and thermal and electrochemical properties were studied. All compounds are ionic liquids (ILs) with melting points between 69 and 93 degrees C. X-ray diffraction proves that all ILs are isostructural. The conductivity at room temperature is between 10(-4) and 10(-8) S cm(-1). Some Cu-based ILs reach conductivities of 10(-2) S cm(-1), which is, however, probably due to IL dec. This correlates with the optical bandgap measurements indicating the formation of large bandgap semiconductors. At elevated temperatures approaching the melting points, the conductivities reach up to 1.47x10(-1) S cm(-1) at 70 degrees C. The electrochemical stability windows of the ILs are between 2.5 and 3.0 V.
High-pressure experiments were performed to investigate the effectiveness, rate and mechanism of carbonation of serpentinites by a carbon-saturated COH fluid at 1.5-2.5 GPa and 375-700 degrees C. This allows a better understanding of the fate and redistribution of slab-derived carbonic fluids when they react with the partially hydrated mantle within and above the subducting slab under pressure and temperature conditions corresponding to the forearc mantle. Interactions between carbon-saturated CO2-H2O-CH4 fluids and serpentinite were investigated using natural serpentinite cylinders with natural grain sizes and shapes in piston-cylinder experiments. The volatile composition of post-run fluids was quantified by gas chromatography. Solid phases were examined by Raman spectroscopy, electron microscopy and laser ablation inductively coupled plasma mass spectrometry. Textures, porosity and phase abundances of recovered rock cores were visualized and quantified by three-dimensional, high-resolution computed tomography. We find that carbonation of serpentinites is efficient at sequestering CO2 from the interacting fluid into newly formed magnesite. Time-series experiments demonstrate that carbonation is completed within similar to 96 h at 2 GPa and 600 degrees C. With decreasing CO2, aq antigorite is replaced first by magnesite + quartz followed by magnesite + talc + chlorite in distinct, metasomatic fronts. Above antigorite stability magnesite + enstatite + talc + chlorite occur additionally. The formation of fluid-permeable reaction zones enhances the reaction rate and efficiency of carbonation. Carbonation probably occurs via an interface-coupled replacement process, whereby interconnected porosity is present within reaction zones after the experiment. Consequently, carbonation of serpentinites is self-promoting and efficient even if fluid flow is channelized into veins. We conclude that significant amounts of carbonates may accumulate, over time, in the hydrated forearc mantle.
Understanding the relation between injection-induced seismic moment release and operational parameters is crucial for early identification of possible seismic hazards associated with fluid-injection projects. We conducted laboratory fluid-injection experiments on permeable sandstone samples containing a critically stressed fault at different fluid pressurization rates. The observed fluid-induced fault deformation is dominantly aseismic. Fluid-induced stick-slip and fault creep reveal that total seismic moment release of acoustic emission (AE) events is related to total injected volume, independent of respective fault slip behavior. Seismic moment release rate of AE scales with measured fault slip velocity. For injection-induced fault slip in a homogeneous pressurized region, released moment shows a linear scaling with injected volume for stable slip (steady slip and fault creep), while we find a cubic relation for dynamic slip. Our results highlight that monitoring evolution of seismic moment release with injected volume in some cases may assist in discriminating between stable slip and unstable runaway ruptures.
Forecasting and early warning systems are important investments to protect lives, properties, and livelihood. While early warning systems are frequently used to predict the magnitude, location, and timing of potentially damaging events, these systems rarely provide impact estimates, such as the expected amount and distribution of physical damage, human consequences, disruption of services, or financial loss. Complementing early warning systems with impact forecasts has a twofold advantage: It would provide decision makers with richer information to take informed decisions about emergency measures and focus the attention of different disciplines on a common target. This would allow capitalizing on synergies between different disciplines and boosting the development of multihazard early warning systems. This review discusses the state of the art in impact forecasting for a wide range of natural hazards. We outline the added value of impact-based warnings compared to hazard forecasting for the emergency phase, indicate challenges and pitfalls, and synthesize the review results across hazard types most relevant for Europe.
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.
Neodymium isotopic composition (epsilon Nd) has enjoyed widespread use as a palaeotracer, principally because it behaves quasi-conservatively in the modern ocean. However, recent bottom water epsilon Nd reconstructions from the eastern North Atlantic are difficult to interpret under assumptions of conservative behaviour. The observation that this apparent departure from conservative behaviour increases with enhanced ice-rafted debris (IRD) fluxes has resulted in the suggestion that IRD leads to the overprinting of bottom water epsilon Nd through reversible scavenging. In this study, a simple water column model successfully reproduces epsilon Nd reconstructions from the eastern North Atlantic at the Last Glacial Maximum and Heinrich Stadial 1, and demonstrates that the changes in scavenging intensity required for good model-data fit is in good agreement with changes in the observed IRD flux. Although uncertainties in model parameters preclude a more definitive conclusion, the results indicate that the suggestion of IRD as a source of non-conservative behaviour in the epsilon Nd tracer is reasonable and that further research into the fundamental chemistry underlying the marine neodymium cycle is necessary to increase confidence in assumptions of conservative epsilon Nd behaviour in the past.
Accelerograms are the primary source for characterizing strong ground motion. It is therefore of paramount interest to have high-quality recordings free from any nonphysical contamination. Frequently, accelerograms are affected by baseline jumps and drifts, either related to the instrument and/or a major earthquake. In this work, I propose a correction method for these undesired baseline drifts based on segmented linear least squares. The algorithm operates on the integrated waveforms and combines all three instrument components to estimate a model that modifies the baseline to be at zero continuously. The procedure consists of two steps: first a suite of models with variable numbers of discontinuities is derived for all three instrument components. During this process, the number of discontinuities is reduced in a parsimonious way, for example, two very close discontinuities are merged into a single one. In the second step, the optimal model is selected on the basis of the Bayesian information criterion. I exemplify the application on synthetic waveforms with known discontinuities and on observed waveforms from a unified strong-motion database of the Japan Meteorological Agency (JMA) and the National Research Institute for Earth Science and Disaster Prevention (NIED, Japan) networks for the major events of the 2016 Kumamoto earthquakes. After the baseline jump correction, the waveforms are furthermore corrected for displacement according to Wang et al.(2011). The resulting displacements are comparable to the Interferometric Synthetic Aperture Radar-derived displacement estimates for the Kumamoto earthquake sequence.
Geochemical processes change the microstructure of rocks and thereby affect their physical behaviour at the macro scale. A micro-computer tomography (micro-CT) scan of a typical reservoir sandstone is used to numerically examine the impact of three spatial alteration patterns on pore morphology, permeability and elastic moduli by correlating precipitation with the local flow velocity magnitude. The results demonstrate that the location of mineral growth strongly affects the permeability decrease with variations by up to four orders in magnitude. Precipitation in regions of high flow velocities is characterised by a predominant clogging of pore throats and a drastic permeability reduction, which can be roughly described by the power law relation with an exponent of 20. A continuous alteration of the pore structure by uniform mineral growth reduces the permeability comparable to the power law with an exponent of four or the Kozeny-Carman relation. Preferential precipitation in regions of low flow velocities predominantly affects smaller throats and pores with a minor impact on the flow regime, where the permeability decrease is considerably below that calculated by the power law with an exponent of two. Despite their complete distinctive impact on hydraulics, the spatial precipitation patterns only slightly affect the increase in elastic rock properties with differences by up to 6.3% between the investigated scenarios. Hence, an adequate characterisation of the spatial precipitation pattern is crucial to quantify changes in hydraulic rock properties, whereas the present study shows that its impact on elastic rock parameters is limited. The calculated relations between porosity and permeability, as well as elastic moduli can be applied for upscaling micro-scale findings to reservoir-scale models to improve their predictive capabilities, what is of paramount importance for a sustainable utilisation of the geological subsurface.
In the field of hydrological modeling, many alternative representations of natural processes exist. Choosing specific process formulations when building a hydrological model is therefore associated with a high degree of ambiguity and subjectivity. In addition, the numerical integration of the underlying differential equations and parametrization of model structures influence model performance. Identifiability analysis may provide guidance by constraining the a priori range of alternatives based on observations. In this work, a flexible simulation environment is used to build an ensemble of semidistributed, process-based hydrological model configurations with alternative process representations, numerical integration schemes, and model parametrizations in an integrated manner. The flexible simulation environment is coupled with an approach for dynamic identifiability analysis. The objective is to investigate the applicability of the framework to identify the most adequate model. While an optimal model configuration could not be clearly distinguished, interesting results were obtained when relating model identifiability with hydro-meteorological boundary conditions. For instance, we tested the Penman-Monteith and Shuttleworth & Wallace evapotranspiration models and found that the former performs better under wet and the latter under dry conditions. Parametrization of model structures plays a dominant role as it can compensate for inadequate process representations and poor numerical solvers. Therefore, it was found that numerical solvers of high order of accuracy do often, though not necessarily, lead to better model performance. The proposed coupled framework proved to be a straightforward diagnostic tool for model building and hypotheses testing and shows potential for more in-depth analysis of process implementations and catchment functioning.
High-resolution surface velocities and strain for Anatolia from Sentinel-1 InSAR and GNSS data
(2020)
Measurements of present-day surface deformation are essential for the assessment of long-term seismic hazard. The European Space Agency's Sentinel-1 satellites enable global, high-resolution observation of crustal motion from Interferometric Synthetic Aperture Radar (InSAR). We have developed automated InSAR processing systems that exploit the first similar to 5 years of Sentinel-1 data to measure surface motions for the similar to 800,000-km(2) Anatolian region. Our new 3-D velocity and strain rate fields illuminate deformation patterns dominated by westward motion of Anatolia relative to Eurasia, localized strain accumulation along the North and East Anatolian Faults, and rapid vertical signals associated with anthropogenic activities and to a lesser extent extension across the grabens of western Anatolia. We show that automatically processed Sentinel-1 InSAR data can characterize details of the velocity and strain rate fields with high resolution and accuracy over large regions. These results are important for assessing the relationship between strain accumulation and release in earthquakes. <br /> Plain Language Summary Satellite-based measurements of small rates of motion of the Earth's surface made at high spatial resolutions and over large areas are important for many geophysical applications including improving earthquake hazard models. We take advantage of recent advances in geodetic techniques in order to measure surface velocities and tectonic strain accumulation across the Anatolia region, including the highly seismogenic and often deadly North Anatolian Fault. We show that by combining Sentinel-1 Interferometric Synthetic Aperture Radar (InSAR) data with Global Navigation Satellite System (GNSS) measurements we can enhance our view of surface deformation associated with active tectonics, the earthquake cycle, and anthropogenic processes.
Subduction of Neo-Tethys oceanic lithosphere beneath the Iranian plate during the Mesozoic formed several igneous bodies of ultramafic to intermediate and felsic composition. Intrusion of these magmas into a regional metamorphic sequence (the Sanandaj-Sirjan Zone) caused partial melting and formation of migmatites with meta-pelitic protoliths. The Alvand complex (west Iran) is a unique area comprising migmatites of both mafic and pelitic protoliths. In this area, the gabbroic rocks contain veins of leucosome at their contact with pyroxenite and olivine gabbro. These leucosomes are geochemically and mineralogically different from leucosomes of the meta-pelitic migmatites and clearly show properties of I-type granites. Microscopic observations and whole rock compositions of the mafic migmatite leucosomes show that migmatization occurred through partial melting of biotite, hornblende and plagioclase. Thermobarometric calculations indicate 800 degrees C and 3.7 kbar for partial melting, although phase diagram modeling demonstrates that the presence of water could decrease the solidus temperature by about 40 degrees C. Our results suggest an asthenospheric magma upwelling as the source of heat for partial melting of the gabbroic rock during subduction of Neo-Tethys oceanic crust under the western edge of the Iranian plate. The present study also reveals relationships between migmatization and formation of S- and I -type granites in the area.
Der zentralasiatische Naturraum, wie er sich uns heute präsentiert, ist das Ergebnis eines Zusammenwirkens vieler verschiedener Faktoren über Jahrmillionen hinweg. Im aktuellen Kontext des Klimawandels zeigt sich jedoch, wie stark sich Stoffflüsse auch kurzfristig ändern und dabei das Gesicht der Landschaft verwandeln können. Die Gobi-Wüste in der Inneren Mongolei (China), als Teil der gleichnamigen Trockenregionen Nordwestchinas, ist aufgrund der Ausgestaltung ihrer landschaftsprägenden Elemente sowie ihrer Landschaftsdynamik, im Zusammenhang mit der Lage zum Tibet-Plateau, in den Fokus der klimageschichtlichen Grundlagenforschung gerückt. Als großes Langzeitarchiv unterschiedlichster fluvialer, lakustriner und äolischer Sedimente stellt sie eine bedeutende Lokalität zur Rekonstruktion von lokalen und regionalen Stoffflüssen dar.. Andererseits ist die Gobi-Wüste zugleich auch eine bedeutende Quelle für den überregionalen Staubtransport, da sie aufgrund der klimatischen Bedingungen insbesondere der Erosion durch Ausblasung preisgegeben wird. Vor diesem Hintergrund erfolgten zwischen 2011 und 2014, im Rahmen des BMBF-Verbundprogramms WTZ Zentralasien – Monsundynamik & Geoökosysteme (Förderkennzeichen 03G0814), mehrere deutsch-chinesische Expeditionen in das Ejina-Becken (Innere Mongolei) und das Qilian Shan-Vorland. Im Zuge dieser Expeditionen wurden für eine Bestimmung potenzieller Sedimentquellen erstmals zahlreiche Oberflächenproben aus dem gesamten Einzugsgebiet des Heihe (schwarzer Fluss) gesammelt. Zudem wurden mit zwei Bohrungen im inneren des Ejina-Beckens, ergänzende Sedimentbohrkerne zum bestehenden Bohrkern D100 (siehe Wünnemann (2005)) abgeteuft, um weit reichende, ergänzende Informationen zur Landschaftsgeschichte und zum überregionalen Sedimenttransfer zu erhalten. Gegenstand und Ziel der vorliegenden Doktorarbeit ist die sedimentologisch-mineralogische Charakterisierung des Untersuchungsgebietes in Bezug auf potenzielle Sedimentquellen und Stoffflüsse des Ejina-Beckens sowie die Rekonstruktion der Ablagerungsgeschichte eines dort erbohrten, 19m langen Sedimentbohrkerns (GN100). Schwerpunkt ist hierbei die Klärung der Sedimentherkunft innerhalb des Bohrkerns sowie die Ausweisung von Herkunftssignalen und möglichen Sedimentquellen bzw. Sedimenttransportpfaden. Die methodische Herangehensweise basiert auf einem Multi-Proxy-Ansatz zur Charakterisierung der klastischen Sedimentfazies anhand von Geländebeobachtungen, lithologisch-granulometrischen und mineralogisch-geochemischen Analysen sowie statistischen Verfahren. Für die mineralogischen Untersuchungen der Sedimente wurde eine neue, rasterelektronenmikroskopische Methode zur automatisierten Partikelanalyse genutzt und den traditionellen Methoden gegenübergestellt. Die synoptische Betrachtung der granulometrischen, geochemischen und mineralogischen Befunde der Oberflächensedimente ergibt für das Untersuchungsgebiet ein logisches Kaskadenmodell mit immer wiederkehrenden Prozessbereichen und ähnlichen Prozesssignalen. Die umfangreichen granulometrischen Analysen deuten dabei auf abnehmende Korngrößen mit zunehmender Entfernung vom Qilian Shan hin und ermöglichen die Identifizierung von vier texturellen Signalen: den fluvialen Sanden, den Dünensanden, den Stillwassersedimenten und Stäuben. Diese Ergebnisse können als Interpretationsgrundlage für die Korngrößenanalysen des Bohrkerns genutzt werden. Somit ist es möglich, die Ablagerungsgeschichte der Bohrkernsedimente zu rekonstruieren und in Verbindung mit eigenen und literaturbasierten Datierungen in einen Gesamtkontext einzuhängen. Für das Untersuchungsgebiet werden somit vier Ablagerungsphasen ausgewiesen, die bis in die Zeit des letzten glazialen Maximums (LGM) zurückreichen. Während dieser Ablagerungsphasen kam es im Zuge unterschiedlicher Aktivitäts- und Stabilitätsphasen zu einer kontinuierlichen Progradation und Überprägung des Schwemmfächers. Eine besonders aktive Phase kann zwischen 8 ka und 4 ka BP festgestellt werden, während der es aufgrund zunehmender fluvialer Aktivitäten zu einer deutlich verstärkten Schwemmfächerdynamik gekommen zu sein scheint. In den Abschnitten davor und danach waren es vor allem äolische Prozesse, die zu einer Überprägung des Schwemmfächers geführt haben. Hinsichtlich der mineralogischen Herkunftssignale gibt es eine große Variabilität. Dies spiegelt die enorme Heterogenität der Geologie des Untersuchungsgebietes wider, wodurch die räumlichen Signale nicht sehr stark ausgeprägt sind. Dennoch, können für das Einzugsgebiet drei größere Bereiche deklariert werden, die als Herkunftsgebiet in Frage kommen. Das östliche Qilian Shan Vorland zeichnet sich dabei durch deutlich höhere Chloritgehalte als primäre Quelle für die Sedimente im Ejina-Becken aus. Sie unterscheiden sich insbesondere durch stark divergierende Chloritgehalte in der Tonmineral- und Gesamtmineralfraktion, was das östliche Qilian Shan Vorland als primäre Quelle für die Sedimente im Ejina-Becken auszeichnet. Dies steht in Zusammenhang mit den Grünschiefern, Ophioliten und Serpentiniten in diesem Bereich. Geochemisch deutet vor allem das Cr/Rb-Verhältnis eine große Variabilität innerhalb des Einzugsgebietes an. Auch hier ist es das östliche Vorland, welches aufgrund seines hohen Anteils an mafischen Gesteinen reich an Chromiten und Spinellen ist und sich somit vom restlichen Untersuchungsgebiet abhebt. Die zeitliche aber auch die generelle Variabilität der Sedimentherkunft lässt sich in den Bohrkernsedimenten nicht so deutlich nachzeichnen. Die mineralogisch-sedimentologischen Eigenschaften der erbohrten klastischen Sedimente zeugen zwar von zwischenzeitlichen Änderungen bei der Sedimentherkunft, diese sind jedoch nicht so deutlich ausgeprägt, wie es die Quellsignale in den Oberflächensedimenten vermuten lassen. Ein Grund dafür scheint die starke Vermischung unterschiedlichster Sedimente während des Transportes zu sein. Die Kombination der Korngrößenergebnisse mit den Befunden der Gesamt- und Schwermineralogie deuten darauf hin, dass es zwischenzeitlich eine Phase mit überwiegend äolischen Prozessen gegeben hat, die mit einem Sedimenteintrag aus dem westlichen Bei Shan in Verbindung stehen. Neben der Zunahme ultrastabiler Schwerminerale wie Zirkon und Granat und der Abnahme opaker Schwerminerale, weisen vor allem die heutigen Verhältnisse darauf hin. Der Vergleich der traditionellen Schwermineralanalyse mit der Computer-Controlled-Scanning-Electron-Microscopy (kurz: CCSEM), die eine automatisierte Partikelauswertung der Proben ermöglicht, zeigt den deutlichen Vorteil der modernen Analysemethode. Neben einem zeitlichen Vorteil, den man durch die automatisierte Abarbeitung der vorbereiteten Proben erlangen kann, steht vor allem die deutlich größere statistische Signifikanz des Ergebnisses im Vordergrund. Zudem können mit dieser Methode auch chemische Varietäten einiger Schwerminerale bestimmt werden, die eine noch feinere Klassifizierung und sicherere Aussagen zu einer möglichen Sedimentherkunft ermöglichen. Damit ergeben sich außerdem verbesserte Aussagen zu Zusammensetzungen und Entstehungsprozessen der abgelagerten Sedimente. Die Studie verdeutlicht, dass die Sedimentherkunft innerhalb des Untersuchungsgebietes sowie die ablaufenden Prozesse zum Teil stark von lokalen Gegebenheiten abhängen. Die Heterogenität der Geologie und die Größe des Einzugsgebietes sowie die daraus resultierende Komplexität der Sedimentgenese, machen exakte Zuordnungen zu klar definierten Sedimentquellen sehr schwer. Dennoch zeigen die Ergebnisse, dass die Sedimentzufuhr in das Ejina-Becken in erster Linie durch fluviale klastische Sedimente des Heihe aus dem Qilian Shan erfolgt sein muss. Die Untersuchungsergebnisse zeigen jedoch ebenso die Notwendigkeit einer ergänzenden Bearbeitung angrenzender Untersuchungsgebiete, wie beispielsweise den Gobi-Altai im Norden oder den Beishan im Westen, sowie die Verdichtung der Oberflächenbeprobung zur feineren Auflösung von lokalen Sedimentquellen.
Earth and environmental sciences rely on detailed information about subsurface processes. Whereas geophysical techniques typically provide highly resolved spatial images, monitoring subsurface processes is often associated with enormous effort and, therefore, is usually limited to point information in time or space. Thus, the development of spatial and temporal continuous field monitoring methods is a major challenge for the understanding of subsurface processes. We have developed a novel method for ground-penetrating-radar (GPR) reflection monitoring of subsurface flow processes under unsaturated conditions and applied it to a hydrological infiltration experiment performed across a periglacial slope deposit in northwest Luxembourg. Our approach relies on a spatial and temporal quasicontinuous data recording and processing, followed by an attribute analysis based on analyzing differences between individual time steps. The results demonstrate the ability of time-lapse GPR monitoring to visualize the spatial and temporal dynamics of preferential flow processes with a spatial resolution in the order of a few decimeters and temporal resolution in the order of a few minutes. We observe excellent agreement with water table information originating from different boreholes. This demonstrates the potential of surface-based GPR reflection monitoring to observe the spatiotemporal dynamics of water movements in the subsurface. It provides valuable, and so far not accessible, information for example in the field of hydrology and pedology that allows studying the actual subsurface processes rather than deducing them from point information.
Glacial-cycle simulations of the Antarctic Ice Sheet with the Parallel Ice Sheet Model (PISM)
(2020)
Simulations of the glacial-interglacial history of the Antarctic Ice Sheet provide insights into dynamic threshold behavior and estimates of the ice sheet's contributions to global sea-level changes for the past, present and future. However, boundary conditions are weakly constrained, in particular at the interface of the ice sheet and the bedrock. Also climatic forcing covering the last glacial cycles is uncertain, as it is based on sparse proxy data. <br /> We use the Parallel Ice Sheet Model (PISM) to investigate the dynamic effects of different choices of input data, e.g., for modern basal heat flux or reconstructions of past changes of sea level and surface temperature. As computational resources are limited, glacial-cycle simulations are performed using a comparably coarse model grid of 16 km and various parameterizations, e.g., for basal sliding, iceberg calving, or for past variations in precipitation and ocean temperatures. In this study we evaluate the model's transient sensitivity to corresponding parameter choices and to different boundary conditions over the last two glacial cycles and provide estimates of involved uncertainties. We also discuss isolated and combined effects of climate and sea-level forcing. Hence, this study serves as a "cookbook" for the growing community of PISM users and paleo-ice sheet modelers in general. <br /> For each of the different model uncertainties with regard to climatic forcing, ice and Earth dynamics, and basal processes, we select one representative model parameter that captures relevant uncertainties and motivates corresponding parameter ranges that bound the observed ice volume at present. The four selected parameters are systematically varied in a parameter ensemble analysis, which is described in a companion paper.
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.
Lakes cover large parts of the climatically sensitive Arctic landscape and respond rapidly to environmental change. Arctic lakes have different origins and include the predominant thermokarst lakes, which are small, young and highly dynamic, as well as large, old and stable glacial lakes. Freshwater diatoms dominate the primary producer community in these lakes and can be used to detect biotic responses to climate and environmental change. We used specific diatom metabarcoding on sedimentary DNA, combined with next-generation sequencing and diatom morphology, to assess diatom diversity in five glacial and 15 thermokarst lakes within the easternmost expanse of the Siberian treeline ecotone in Chukotka, Russia. We obtained 163 verified diatom sequence types and identified 176 diatom species morphologically. Although there were large differences in taxonomic assignment using the two approaches, they showed similar high abundances and diversity of Fragilariceae and Aulacoseiraceae. In particular, the genetic approach detected hidden within-lake variations of fragilarioids in glacial lakes and dominance of centric Aulacoseira species, whereas Lindavia ocellata was predominant using morphology. In thermokarst lakes, sequence types and valve counts also detected high diversity of Fragilariaceae, which followed the vegetation gradient along the treeline. Ordination analyses of the genetic data from glacial and thermokarst lakes suggest that concentrations of sulfate (SO42-), an indicator of the activity of sulfate-reducing microbes under anoxic conditions, and bicarbonate (HCO3-), which relates to surrounding vegetation, have a significant influence on diatom community composition. For thermokarst lakes, we also identified lake depth as an important variable, but SO42- best explains diatom diversity derived from genetic data, whereas HCO3- best explains the data from valve counts. Higher diatom diversity was detected in glacial lakes, most likely related to greater lake age and different edaphic settings, which gave rise to diversification and endemism. In contrast, small, dynamic thermokarst lakes are inhabited by stress-tolerant fragilarioids and are related to different vegetation types along the treeline ecotone. Our study demonstrated that genetic investigations of lake sediments can be used to interpret climate and environmental responses of diatoms. It also showed how lake type affects diatom diversity, and that such genetic analyses can be used to track diatom community changes under ongoing warming in the Arctic.
Full-waveform-based characterization of acoustic emission activity in a mine-scale experiment
(2020)
Understanding fracturing processes and the hydromechanical relation to induced seismicity is a key question for enhanced geothermal systems (EGS). Commonly massive fluid injection, predominately causing hydroshearing, are used in large-scale EGS but also hydraulic fracturing approaches were discussed. To evaluate the applicability of hydraulic fracturing techniques in EGS, six in situ, multistage hydraulic fracturing experiments with three different injection schemes were performed under controlled conditions in crystalline rock at the Aspo Hard Rock Laboratory (Sweden). During the experiments the near-field ground motion was continuously recorded by 11 piezoelectric borehole sensors with a sampling rate of 1 MHz. The sensor network covered a volume of 30x30x30 m around a horizontal, 28-m-long injection borehole at a depth of 410 m. To extract and characterize massive, induced, high-frequency acoustic emission (AE) activity from continuous recordings, a semi-automated workflow was developed relying on full waveform based detection, classification and location procedures. The approach extended the AE catalogue from 196 triggered events in previous studies to more than 19600 located AEs. The enhanced catalogue, for the first time, allows a detailed analysis of induced seismicity during single hydraulic fracturing experiments, including the individual fracturing stages and the comparison between injection schemes. Beside the detailed study of the spatio-temporal patterns, event clusters and the growth of seismic clouds, we estimate relative magnitudes and b-values of AEs for conventional, cyclic progressive and dynamic pulse injection schemes, the latter two being fatigue hydraulic fracturing techniques. While the conventional fracturing leads to AE patterns clustered in planar regions, indicating the generation of a single main fracture plane, the cyclic progressive injection scheme results in a more diffuse, cloud-like AE distribution, indicating the activation of a more complex fracture network. For a given amount of hydraulic energy (pressure multiplied by injected volume) pumped into the system, the cyclic progressive scheme is characterized by a lower rate of seismicity, lower maximum magnitudes and significantly larger b-values, implying an increased number of small events relative to the large ones. To our knowledge, this is the first direct comparison of high resolution seismicity in a mine-scale experiment induced by different hydraulic fracturing schemes.
In this paper, we present the new frequency spectrum recurrence analysis technique by means of electro-encephalon signals (EES) analyses. The technique is suitable for time series analysis with noise and disturbances. EES were collected, and alpha waves of the occipital region were analysed by comparing the signals from participants in two states, eyes open and eyes closed. Firstly, EES were characterized and analysed by means of techniques already known to compare with the results of the innovative technique that we present here. We verified that, standard recurrence quantification analysis by means of EES time series cannot statistically distinguish the two states. However, the new frequency spectrum recurrence quantification exhibit quantitatively whether the participants have their eyes open or closed. In sequence, new quantifiers are created for analysing the recurrence concentration on frequency bands. These analyses show that EES with similar frequency spectrum have different recurrence levels revealing different behaviours of the nervous system. The technique can be used to deepen the study on depression, stress, concentration level and other neurological issues and also can be used in any complex system.
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.
Filling in for Governments?
(2020)
The 2015 Paris Agreement on climate change abandons the Kyoto Protocol's paradigm of binding emissions targets and relies instead on countries' voluntary contributions. However, the Paris Agreement encourages not only governments but also sub-national governments, corporations and civil society to contribute to reaching ambitious climate goals. In a transition from the regulated architecture of the Kyoto Protocol to the open system of the Paris Agreement, the Agreement seeks to integrate non-state actors into the treaty-based climate regime. In 2014 the secretariat of the United Nations Framework Convention on Climate Change Peru and France created the Non-State Actor Zone for Climate Action (and launched the Global Climate Action portal). In December 2019, this portal recorded more than twenty thousand climate-commitments of private and public non-state entities, making the non-state venues of international climate meetings decisively more exciting than the formal negotiation space. This level engagement and governments' response to it raises a flurry of questions in relation to the evolving nature of the climate regime and climate change governance, including the role of private actors as standard setters and the lack of accountability mechanisms for non-state actions. This paper takes these developments as occasion to discuss the changing role of private actors in the climate regime.