Gold Open-Access
Refine
Document Type
- Article (239)
- Review (4)
- Other (3)
- Doctoral Thesis (1)
Is part of the Bibliography
- yes (247)
Keywords
- climate change (9)
- permafrost (9)
- reactive transport (5)
- PHREEQC (4)
- digital rock physics (3)
- micro-CT (3)
- model (3)
- precipitation (3)
- remote sensing (3)
- time series analysis (3)
Institute
- Institut für Geowissenschaften (247) (remove)
Retrogressive thaw slumps (RTSs) are among the most active thermokarst landforms in the Arctic and deliver a large amount of material to the Arctic Ocean. However, their contribution to the organic carbon (OC) budget is unknown. We provide the first estimate of the contribution of RTSs to the nearshore OC budget of the Yukon Coast, Canada, and describe the evolution of coastal RTSs between 1952 and 2011 in this area. We (1) describe the evolution of RTSs between 1952 and 2011; (2) calculate the volume of eroded material and stocks of OC mobilized through slumping, including soil organic carbon (SOC) and dissolved organic carbon (DOC); and (3) estimate the OC fluxes mobilized through slumping between 1972 and 2011. We identified RTSs using high- resolution satellite imagery from 2011 and geocoded aerial photographs from 1952 and 1972. To estimate the volume of eroded material, we applied spline interpolation on an airborne lidar dataset acquired in July 2013. We inferred the stocks of mobilized SOC and DOC from existing related literature. Our results show a 73% increase in the number of RTSs and 14% areal expansion between 1952 and 2011. In the study area, RTSs displaced at least 16.6 x 10(6) m(3) of material, 53% of which was ice, and mobilized 145.9 x 10(6) kg of OC. Between 1972 and 2011, 49 RTSs displaced 8.6 x 10(3) m(3) yr(-1) of material, adding 0.6% to the OC flux released by coastal retreat along the Yukon Coast. Our results show that the contribution of RTSs to the nearshore OC budget is non-negligible and should be included when estimating the quantity of OC released from the Arctic coast to the ocean.
Salt marshes filter pollutants, protect coastlines against storm surges, and sequester carbon, yet are under threat from sea level rise and anthropogenic modification. The sustained existence of the salt marsh ecosystem depends on the topographic evolution of marsh platforms. Quantifying marsh platform topography is vital for improving the management of these valuable landscapes. The determination of platform boundaries currently relies on supervised classification methods requiring near-infrared data to detect vegetation, or demands labour-intensive field surveys and digitisation. We propose a novel, unsupervised method to reproducibly isolate salt marsh scarps and platforms from a digital elevation model (DEM), referred to as Topographic Identification of Platforms (TIP). Field observations and numerical models show that salt marshes mature into subhorizontal platforms delineated by subvertical scarps. Based on this premise, we identify scarps as lines of local maxima on a slope raster, then fill landmasses from the scarps upward, thus isolating mature marsh platforms. We test the TIP method using lidar-derived DEMs from six salt marshes in England with varying tidal ranges and geometries, for which topographic platforms were manually isolated from tidal flats. Agreement between manual and unsupervised classification exceeds 94% for DEM resolutions of 1 m, with all but one site maintaining an accuracy superior to 90% for resolutions up to 3 m. For resolutions of 1 m, platforms detected with the TIP method are comparable in surface area to digitised platforms and have similar elevation distributions. We also find that our method allows for the accurate detection of local block failures as small as 3 times the DEM resolution. Detailed inspection reveals that although tidal creeks were digitised as part of the marsh platform, unsupervised classification categorises them as part of the tidal flat, causing an increase in false negatives and overall platform perimeter. This suggests our method may benefit from combination with existing creek detection algorithms. Fallen blocks and high tidal flat portions, associated with potential pioneer zones, can also lead to differences between our method and supervised mapping. Although pioneer zones prove difficult to classify using a topographic method, we suggest that these transition areas should be considered when analysing erosion and accretion processes, particularly in the case of incipient marsh platforms. Ultimately, we have shown that unsupervised classification of marsh platforms from high-resolution topography is possible and sufficient to monitor and analyse topographic evolution.
One of the main purposes of detrital thermochronology is to provide constraints on the regional-scale exhumation rate and its spatial variability in actively eroding mountain ranges. Procedures that use cooling age distributions coupled with hypsometry and thermal models have been developed in order to extract quantitative estimates of erosion rate and its spatial distribution, assuming steady state between tectonic uplift and erosion. This hypothesis precludes the use of these procedures to assess the likely transient response of mountain belts to changes in tectonic or climatic forcing. Other methods are based on an a priori knowledge of the in situ distribution of ages to interpret the detrital age distributions. In this paper, we describe a simple method that, using the observed detrital mineral age distributions collected along a river, allows us to extract information about the relative distribution of erosion rates in an eroding catchment without relying on a steady-state assumption, the value of thermal parameters or an a priori knowledge of in situ age distributions. The model is based on a relatively low number of parameters describing lithological variability among the various sub-catchments and their sizes and only uses the raw ages. The method we propose is tested against synthetic age distributions to demonstrate its accuracy and the optimum conditions for it use. In order to illustrate the method, we invert age distributions collected along the main trunk of the Tsangpo-Siang-Brahmaputra river system in the eastern Himalaya. From the inversion of the cooling age distributions we predict present-day erosion rates of the catchments along the Tsangpo-Siang-Brahmaputra river system, as well as some of its tributaries. We show that detrital age distributions contain dual information about present-day erosion rate, i. e., from the predicted distribution of surface ages within each catchment and from the relative contribution of any given catchment to the river distribution. The method additionally allows comparing modern erosion rates to long-term exhumation rates. We provide a simple implementation of the method in Python code within a Jupyter Notebook that includes the data used in this paper for illustration purposes.
Quantifying rock weakening due to decreasing calcite mineral content by numerical simulations
(2018)
The quantification of changes in geomechanical properties due to chemical reactions is of paramount importance for geological subsurface utilisation, since mineral dissolution generally reduces rock stiffness. In the present study, the effective elastic moduli of two digital rock samples, the Fontainebleau and Bentheim sandstones, are numerically determined based on micro-CT images. Reduction in rock stiffness due to the dissolution of 10% calcite cement by volume out of the pore network is quantified for three synthetic spatial calcite distributions (coating, partial filling and random) using representative sub-cubes derived from the digital rock samples. Due to the reduced calcite content, bulk and shear moduli decrease by 34% and 38% in maximum, respectively. Total porosity is clearly the dominant parameter, while spatial calcite distribution has a minor impact, except for a randomly chosen cement distribution within the pore network. Moreover, applying an initial stiffness reduced by 47% for the calcite cement results only in a slightly weaker mechanical behaviour. Using the quantitative approach introduced here substantially improves the accuracy of predictions in elastic rock properties compared to general analytical methods, and further enables quantification of uncertainties related to spatial variations in porosity and mineral distribution.
By using 3-year global positioning system (GPS) measurements from December 2013 to November 2016, we provide in this study a detailed survey on the climatology of the GPS signal loss of Swarm onboard receivers. Our results show that the GPS signal losses prefer to occur at both low latitudes between +/- 5 and +/- 20 degrees magnetic latitude (MLAT) and high latitudes above 60 degrees MLAT in both hemispheres. These events at all latitudes are observed mainly during equinoxes and December solstice months, while totally absent during June solstice months. At low latitudes the GPS signal losses are caused by the equatorial plasma irregularities shortly after sunset, and at high latitude they are also highly related to the large density gradients associated with ionospheric irregularities. Additionally, the high-latitude events are more often observed in the Southern Hemisphere, occurring mainly at the cusp region and along nightside auroral latitudes. The signal losses mainly happen for those GPS rays with elevation angles less than 20 degrees, and more commonly occur when the line of sight between GPS and Swarm satellites is aligned with the shell structure of plasma irregularities. Our results also confirm that the capability of the Swarm receiver has been improved after the bandwidth of the phase-locked loop (PLL) widened, but the updates cannot radically avoid the interruption in tracking GPS satellites caused by the ionospheric plasma irregularities. Additionally, after the PLL bandwidth increased larger than 0.5 Hz, some unexpected signal losses are observed even at middle latitudes, which are not related to the ionospheric plasma irregularities. Our results suggest that rather than 1.0 Hz, a PLL bandwidth of 0.5 Hz is a more suitable value for the Swarm receiver.
Effects of upper mantle heterogeneities on the lithospheric stress field and dynamic topography
(2018)
The orientation and tectonic regime of the observed crustal/lithospheric stress field contribute to our knowledge of different deformation processes occurring within the Earth's crust and lithosphere. In this study, we analyze the influence of the thermal and density structure of the upper mantle on the lithospheric stress field and topography. We use a 3-D lithosphere–asthenosphere numerical model with power-law rheology, coupled to a spectral mantle flow code at 300 km depth. Our results are validated against the World Stress Map 2016 (WSM2016) and the observation-based residual topography. We derive the upper mantle thermal structure from either a heat flow model combined with a seafloor age model (TM1) or a global S-wave velocity model (TM2). We show that lateral density heterogeneities in the upper 300 km have a limited influence on the modeled horizontal stress field as opposed to the resulting dynamic topography that appears more sensitive to such heterogeneities. The modeled stress field directions, using only the mantle heterogeneities below 300 km, are not perturbed much when the effects of lithosphere and crust above 300 km are added. In contrast, modeled stress magnitudes and dynamic topography are to a greater extent controlled by the upper mantle density structure. After correction for the chemical depletion of continents, the TM2 model leads to a much better fit with the observed residual topography giving a good correlation of 0.51 in continents, but this correction leads to no significant improvement of the fit between the WSM2016 and the resulting lithosphere stresses. In continental regions with abundant heat flow data, TM1 results in relatively small angular misfits. For example, in western Europe the misfit between the modeled and observation-based stress is 18.3°. Our findings emphasize that the relative contributions coming from shallow and deep mantle dynamic forces are quite different for the lithospheric stress field and dynamic topography.
Subsurface lateral flow in hillslope soils depends on lower permeability or texture-contrasting soil horizons. In the arable hummocky soil landscape, erosion processes caused glacial till appearance closer to the soil surface at upslope positions. The objective of this work was to quantify the potential for subsurface lateral flow depending on the erosion-affected spatial hydropedological complexity. The eroded Haplic Luvisol profile was studied due to the presence of a relatively dense C horizon that varied in depth, thickness, and sloping angle. A two-dimensional numerical modeling and sensitivity analysis for the saturated hydraulic conductivity (K-s) of the C horizon and the depth to C horizon (i.e., soil solum thickness) was performed for rainstorms in 2011 and 2012 using HYDRUS-2D. A K-s value of <2.5 cm d(-1) for the C horizon was required for lateral flow initiation. Lateral flow was (i) increasing with decreasing solum thickness, indicating an erosion-induced feedback on subsurface lateral flow, and (ii) dependent on the soil moisture regime prior to rainstorms. The effect of lateral flow on the movement of a conservative tracer was simulated in the form of a "virtual experiment". Simulation scenarios revealed only a relatively small lateral shift of the tracer plume caused by a local decoupling of water (already lateral) from subsequent tracer movement (still vertical). Longer term simulations suggested that, for the present conditions, lateral flow was limited mostly to occasional summer storm events. Even without considering preferential flow contribution to lateral flow, highly complex hydropedologic interactions are present in erosion-affected heterogeneous soil profiles.
The oldest ice core records are obtained from the East Antarctic Plateau. Water isotopes are key proxies to reconstructing past climatic conditions over the ice sheet and at the evaporation source. The accuracy of climate reconstructions depends on knowledge of all processes affecting water vapour, precipitation and snow isotopic compositions. Fractionation processes are well understood and can be integrated in trajectory-based Rayleigh distillation and isotope-enabled climate models. However, a quantitative understanding of processes potentially altering snow isotopic composition after deposition is still missing. In low-accumulation sites, such as those found in East Antarctica, these poorly constrained processes are likely to play a significant role and limit the interpretability of an ice core's isotopic composition.
By combining observations of isotopic composition in vapour, precipitation, surface snow and buried snow from Dome C, a deep ice core site on the East Antarctic Plateau, we found indications of a seasonal impact of metamorphism on the surface snow isotopic signal when compared to the initial precipitation. Particularly in summer, exchanges of water molecules between vapour and snow are driven by the diurnal sublimation-condensation cycles. Overall, we observe in between precipitation events modification of the surface snow isotopic composition. Using high-resolution water isotopic composition profiles from snow pits at five Antarctic sites with different accumulation rates, we identified common patterns which cannot be attributed to the seasonal variability of precipitation. These differences in the precipitation, surface snow and buried snow isotopic composition provide evidence of post-deposition processes affecting ice core records in low-accumulation areas.
Most large-scale hydrologic models fall short in reproducing groundwater head dynamics and simulating transport process due to their oversimplified representation of groundwater flow. In this study, we aim to extend the applicability of the mesoscale Hydrologic Model (mHM v5.7) to subsurface hydrology by coupling it with the porous media simulator OpenGeoSys (OGS). The two models are one-way coupled through model interfaces GIS2FEM and RIV2FEM, by which the grid-based fluxes of groundwater recharge and the river-groundwater exchange generated by mHM are converted to fixed-flux boundary conditions of the groundwater model OGS. Specifically, the grid-based vertical reservoirs in mHM are completely preserved for the estimation of land-surface fluxes, while OGS acts as a plug-in to the original mHM modeling framework for groundwater flow and transport modeling. The applicability of the coupled model (mHM-OGS v1.0) is evaluated by a case study in the central European mesoscale river basin - Nagelstedt. Different time steps, i.e., daily in mHM and monthly in OGS, are used to account for fast surface flow and slow groundwater flow. Model calibration is conducted following a two-step procedure using discharge for mHM and long-term mean of groundwater head measurements for OGS. Based on the model summary statistics, namely the Nash-Sutcliffe model efficiency (NSE), the mean absolute error (MAE), and the interquartile range error (QRE), the coupled model is able to satisfactorily represent the dynamics of discharge and groundwater heads at several locations across the study basin. Our exemplary calculations show that the one-way coupled model can take advantage of the spatially explicit modeling capabilities of surface and groundwater hydrologic models and provide an adequate representation of the spatiotemporal behaviors of groundwater storage and heads, thus making it a valuable tool for addressing water resources and management problems.
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.
Landslide hazard motivates the need for a deeper understanding of the events that occur before, during, and after catastrophic slope failures. Due to the destructive nature of such events, in situ observation is often difficult or impossible. Here, we use data from a network of 58 seismic stations to characterise a large landslide at the Askja caldera, Iceland, on 21 July 2014. High data quality and extensive network coverage allow us to analyse both long- and short-period signals associated with the landslide, and thereby obtain information about its triggering, initiation, timing, and propagation. At long periods, a landslide force history inversion shows that the Askja landslide was a single, large event starting at the SE corner of the caldera lake at 23:24:05 UTC and propagating to the NW in the following 2 min The bulk sliding mass was 7-16 x 10(10) kg, equivalent to a collapsed volume of 35-80 x 10(6) m(3). The sliding mass was displaced downslope by 1260 +/- 250 m. At short periods, a seismic tremor was observed for 30 min before the landslide. The tremor is approximately harmonic with a fundamental frequency of 2.3 Hz and shows time-dependent changes of its frequency content. We attribute the seismic tremor to stick-slip motion along the landslide failure plane. Accelerating motion leading up to the catastrophic slope failure culminated in an aseismic quiescent period for 2 min before the landslide. We propose that precursory seismic signals may be useful in landslide early-warning systems. The 8 h after the main landslide failure are characterised by smaller slope failures originating from the destabilised caldera wall decaying in frequency and magnitude. We introduce the term "afterslides" for this subsequent, declining slope activity after a large landslide.
A comprehensive hydro-sedimentological dataset for the Isabena catchment, northeastern (NE) Spain, for the period 2010-2018 is presented to analyse water and sediment fluxes in a Mediterranean mesoscale catchment. The dataset includes rainfall data from 12 rain gauges distributed within the study area complemented by meteorological data of 12 official meteo-stations. It comprises discharge data derived from water stage measurements as well as suspended sediment concentrations (SSCs) at six gauging stations of the River Isabena and its sub-catchments. Soil spectroscopic data from 351 suspended sediment samples and 152 soil samples were collected to characterize sediment source regions and sediment properties via fingerprinting analyses. The Isabena catchment (445 km(2)) is located in the southern central Pyrenees ranging from 450 m to 2720 m a.s.l.; together with a pronounced topography, this leads to distinct temperature and precipitation gradients. The River Isabena shows marked discharge variations and high sediment yields causing severe siltation problems in the downstream Barasona Reservoir. The main sediment source is badland areas located on Eocene marls that are well connected to the river network. The dataset features a comprehensive set of variables in a high spatial and temporal resolution suitable for the advanced process understanding of water and sediment fluxes, their origin and connectivity and sediment budgeting and for the evaluation and further development of hydro-sedimentological models in Mediterranean mesoscale mountainous catchments.
It is well-known that prey species often face trade-offs between defense against predation and competitiveness, enabling predator-mediated coexistence. However, we lack an understanding of how the large variety of different defense traits with different competition costs affects coexistence and population dynamics. Our study focusses on two general defense mechanisms, that is, pre-attack (e.g., camouflage) and post-attack defenses (e.g., weaponry) that act at different phases of the predator—prey interaction. We consider a food web model with one predator, two prey types and one resource. One prey type is undefended, while the other one is pre- or post-attack defended paying costs either by a higher half-saturation constant for resource uptake or a lower maximum growth rate. We show that post-attack defenses promote prey coexistence and stabilize the population dynamics more strongly than pre-attack defenses by interfering with the predator's functional response: Because the predator spends time handling “noncrackable” prey, the undefended prey is indirectly facilitated. A high half-saturation constant as defense costs promotes coexistence more and stabilizes the dynamics less than a low maximum growth rate. The former imposes high costs at low resource concentrations but allows for temporally high growth rates at predator-induced resource peaks preventing the extinction of the defended prey. We evaluate the effects of the different defense mechanisms and costs on coexistence under different enrichment levels in order to vary the importance of bottom-up and top-down control of the prey community.
The intensification of Northern Hemisphere glaciation (iNHG) and uplift of the Tibetan Plateau have been argued to be among the main drivers of climate change in midlatitude Central Asia during the Pliocene/Pleistocene. While most proxy records that support this hypothesis are from regions outside the Tibetan Plateau (such as from the Chinese Loess Plateau), detailed paleoclimatic information for the plateau itself during that time has yet remained elusive. Here we present a temporally highly resolved (similar to 500 years) sedimentological record from the Qaidam Basin situated on the northeastern Tibetan Plateau that shows pronounced glacial-interglacial climate variability during the interval from 2.7 to 2.1 Ma. Glacial (interglacial) intervals are generally characterized by coarser (finer) grain size, minima (maxima) in organic matter content, and maxima (minima) in carbonate content. Comparison of our results with Earth's orbital parameters and proxy records from the Chinese Loess Plateau suggests that the observed climate fluctuations were mainly driven by changes in the Siberian High/East Asian winter monsoon system as a response to the iNHG. They are further proposed to be enhanced by the topography of the Tibetan Plateau and its impact on the position and intensity of the westerlies.
The Cluster mission has produced a large data set of electron flux measurements in the Earth's magnetosphere since its launch in late 2000. Electron fluxes are measured using Research with Adaptive Particle Imaging Detector (RAPID)/Imaging Electron Spectrometer (IES) detector as a function of energy, pitch angle, spacecraft position, and time. However, no adiabatic invariants have been calculated for Cluster so far. In this paper we present a step-by-step guide to calculations of adiabatic invariants and conversion of the electron flux to phase space density (PSD) in these coordinates. The electron flux is measured in two RAPID/IES energy channels providing pitch angle distribution at energies 39.2-50.5 and 68.1-94.5 keV in nominal mode since 2004. A fitting method allows to expand the conversion of the differential fluxes to the range from 40 to 150 keV. Best data coverage for phase space density in adiabatic invariant coordinates can be obtained for values of second adiabatic invariant, K, similar to 10(2), and values of the first adiabatic invariant mu in the range approximate to 5-20 MeV/G. Furthermore, we describe the production of a new data product "LSTAR," equivalent to the third adiabatic invariant, available through the Cluster Science Archive for years 2001-2018 with 1-min resolution. The produced data set adds to the availability of observations in Earth's radiation belts region and can be used for long-term statistical purposes.
Millennial-scale palaeoclimate variability has been documented in various terrestrial and marine palaeoclimate proxy records throughout the Northern Hemisphere for the last glacial cycle. Its clear expression and rapid shifts between different states of climate (Greenland Interstadials and Stadials) represents a correlation tool beyond the resolution of e.g. luminescence dating, especially relevant for terrestrial deposits. Usually, comparison of terrestrial proxy datasets and the Greenland ice cores indicates a complex expression of millennial-scale climate variability as recorded in terrestrial geoarchives including loess. Loess is the most widespread terrestrial geoarchive of the Quaternary and especially widespread over Eurasia. However, loess often records a smoothed representation of millennial-scale variability without all fidelity when compared to the Greenland data, this being a relevant limiting feature in integrating loess with other palaeoclimate records. To better understand the loess proxy-response to millennial-scale climate variability, we simulate a proxy signal smoothing by natural processes through application of low-pass filters of delta O-18 data from Greenland, a high-resolution palaeoclimate reference record, alongside speleothem isotope records from the Black Sea-Mediterranean region. We show that low-pass filters represent rather simple models for better constraining the expression of millennial-scale climate variability in low sedimentation environments, and in sediments where proxy-response signals are most likely affected by natural smoothing (by e.g. bioturbation). Interestingly, smoothed datasets from Greenland and the Black Sea-Mediterranean region are most similar in the last similar to 15 ka and between similar to 50-30 ka. Between similar to 30-15 ka, roughly corresponding to the Last Glacial Maximum and the deglaciation, the records show dissimilarities, challenging the construction of robust correlative time-scales in this age range. From our analysis it becomes apparent that patterns of palaeoclimate signals in loess-palaeosol sequences often might be better explained by smoothed Greenland reference data than the original high-resolution Greenland dataset, or other reference data. This opens the possibility to better assess the temporal resolution and palaeoclimate potential of loess-palaeosol sequences in recording supra-regional climate patterns, as well as to securely integrate loess with other chronologically better-resolved palaeoclimate records.
Observations of rift and rifted margin architecture suggest that significant spatial and temporal structural heterogeneity develops during the multiphase evolution of continental rifting. Inheritance is often invoked to explain this heterogeneity, such as preexisting anisotropies in rock composition, rheology, and deformation. Here, we use high-resolution 3-D thermal-mechanical numerical models of continental extension to demonstrate that rift-parallel heterogeneity may develop solely through fault network evolution during the transition from distributed to localized deformation. In our models, the initial phase of distributed normal faulting is seeded through randomized initial strength perturbations in an otherwise laterally homogeneous lithosphere extending at a constant rate. Continued extension localizes deformation onto lithosphere-scale faults, which are laterally offset by tens of km and discontinuous along-strike. These results demonstrate that rift- and margin-parallel heterogeneity of large-scale fault patterns may in-part be a natural byproduct of fault network coalescence.
Temporal variation of natural light sources such as airglow limits the ability of night light sensors to detect changes in small sources of artificial light (such as villages). This study presents a method for correcting for this effect globally, using the satellite radiance detected from regions without artificial light emissions. We developed a routine to define an approximate grid of locations worldwide that do not have regular light emission. We apply this method with a 5 degree equally spaced global grid (total of 2016 individual locations), using data from the Visible Infrared Imaging Radiometer Suite (VIIRS) Day-Night Band (DNB). This code could easily be adapted for other future global sensors. The correction reduces the standard deviation of data in the Earth Observation Group monthly DNB composites by almost a factor of two. The code and datasets presented here are available under an open license by GFZ Data Services, and are implemented in the Radiance Light Trends web application.
Nitraria is a halophytic taxon (i.e., adapted to saline environments) that belongs to the plant family Nitrariaceae and is distributed from the Mediterranean, across Asia into the south-eastern tip of Australia. This taxon is thought to have originated in Asia during the Paleogene (66-23 Ma), alongside the proto-Paratethys epicontinental sea. The evolutionary history of Nitraria might hold important clues on the links between climatic and biotic evolution but limited taxonomic documentation of this taxon has thus far hindered this line of research. Here we investigate if the pollen morphology and the chemical composition of the pollen wall are informative of the evolutionary history of Nitraria and could explain if origination along the proto-Paratethys and dispersal to the Tibetan Plateau was simultaneous or a secondary process. To answer these questions, we applied a novel approach consisting of a combination of Fourier Transform Infrared spectroscopy (FTIR), to determine the chemical composition of the pollen wall, and pollen morphological analyses using Light Microscopy (LM) and Scanning Electron Microscopy (SEM). We analysed our data using ordinations (principal components analysis and non-metric multidimensional scaling), and directly mapped it on the Nitrariaceae phylogeny to produce a phylomorphospace and a phylochemospace. Our LM, SEM and FTIR analyses show clear morphological and chemical differences between the sister groups Peganum and Nitraria. Differences in the morphological and chemical characteristics of highland species (Nitraria schoberi, N. sphaerocarpa, N. sibirica and N. tangutorum) and lowland species (Nitraria billardierei and N. retusa) are very subtle, with phylogenetic history appearing to be a more important control on Nitraria pollen than local environmental conditions. Our approach shows a compelling consistency between the chemical and morphological characteristics of the eight studied Nitrariaceae species, and these traits are in agreement with the phylogenetic tree. Taken together, this demonstrates how novel methods for studying fossil pollen can facilitate the evolutionary investigation of living and extinct taxa, and the environments they represent.
Relative abundances of 157 diatom taxa from Yakutian lake surface-sediments were investigated for their potential to indicate certain environmental conditions. Data from 206 sites from Arctic, sub-Arctic and boreal environments were included. Redundancy analyses were performed to assess the explanatory power of mean July temperature (T-July), conductivity, pH, dissolved silica concentration, phosphate concentration, lake depth and vegetation type on diatom species composition. Boosted regression tree analyses were performed to infer the most relevant environmental variables for abundances of individual taxa and weighted average regression was applied to infer their respective optimum and tolerance. Electrical conductivity was best indicated by diatom taxa. In contrast, only few taxa were indicative of Si and water depth. Few taxa were related to specific pH values. Although T-July, explained the highest proportion of variance in the diatom spectra and was, after conductivity, the second-most selected splitting variable, we a priori decided not to present indicator taxa because of the poorly understood relationship between diatom occurrences and T-July. In total, 92 diatom taxa were reliable indicators of a certain vegetation type or a combination of several types. The high numbers of indicative species for open vegetation sites and for forested sites suggest that the principal turnover is the transition from forest-tundra to northern taiga. Overall, our results reveal that preference ranges of diatom taxa for environmental variables are mostly broad, and the use of indicator taxa for the purposes of environmental reconstruction or environmental monitoring is therefore restricted to marked rather than subtle environmental transitions.
The Sesia zone in the Italian Western Alps is a piece of continental crust that has been subducted to eclogite-facies conditions and records a complex metamorphic history. The exact timing of events and the significance of geochronological information are debated due to the interplay of tectonic, metamorphic, and metasomatic processes. Here we present new geochronological data using Rb-Sr internal mineral isochrons and in situ Ar-40/Ar-39 laser ablation data to provide constraints on the relative importance of fluid-mediated mineral replacement reactions and diffusion for the interpretation of radiogenic isotope signatures, and on the use of these isotopic systems for dating metamorphic and variably deformed rocks. Our study focuses on the shear zone at the contact between two major lithological units of the Sesia zone, the eclogitic micaschists and the gneiss minuti. Metasedimentary rocks of the eclogitic micaschists unit contain phengite with step-like zoning in major element chemistry as evidence for petrologic disequilibrium. Distinct Ar-40/Ar-39 spot ages of relict phengite cores and over-printed rims demonstrate the preservation of individual age domains in the crystals. The eclogitic micaschists also show systematic Sr isotope disequilibria among different phengite populations, so that minimum ages of relict assemblage crystallization can be differentiated from the timing of late increments of deformation. The preservation of these disequilibrium features shows the lack of diffusive re-equilibration and underpins that fluid-assisted dissolution and recrystallization reactions are the main factors controlling the isotope record in these subduction-related metamorphic rocks. Blueschist-facies mylonites record deformation along the major shear zone that separates the eclogitic micaschists from the gneiss minuti. Two Rb-Sr isochrones that comprise several white mica fractions and glaucophane constrain the timing of this deformation and accompanying near-complete blueschist-facies re-equilibration of the Rb-Sr system to 60.1 +/- 0.9 Ma and 60.9 +/- 2.1 Ma, respectively. Overlapping ages in eclogitic micaschists of 60.1 +/- 1.1 (Rb-Sr isochron of sheared matrix assemblage), 58.6 +/- 0.8, and 60.9 +/- 0.4 Ma (white mica Ar-40/Ar-39 inverse isochron ages) support the significance of this age and show that fluid-rock interaction and partial re-equilibration occurred as much as several kilometers away from the shear zone. An earlier equilibration during high-pressure conditions in the eclogitic mica schists is recorded in minimum Rb-Sr ages for relict assemblages (77.2 +/- 0.8 and 72.4 +/- 1.1 Ma) and an Ar-40/Ar-39 inverse isochron age of 75.4 +/- 0.8 Ma for white mica cores, again demonstrating that the two isotope systems provide mutually supporting geochronological information. Local reactivation and recrystallization along the shear zone lasted >15 m.y., as late increments of deformation are recorded in a greenschist-facies mylonite by a Rb-Sr isochron age of 46.5 +/- 0.7 Ma.
Atmospheric water vapour content is a key variable that controls the development of deep convective storms and rainfall extremes over the central Andes. Direct measurements of water vapour are challenging; however, recent developments in microwave processing allow the use of phase delays from L-band radar to measure the water vapour content throughout the atmosphere: Global Navigation Satellite System (GNSS)-based integrated water vapour (IWV) monitoring shows promising results to measure vertically integrated water vapour at high temporal resolutions. Previous works also identified convective available potential energy (CAPE) as a key climatic variable for the formation of deep convective storms and rainfall in the central Andes. Our analysis relies on GNSS data from the Argentine Continuous Satellite Monitoring Network, Red Argentina de Monitoreo Satelital Continuo (RAMSAC) network from 1999 to 2013. CAPE is derived from version 2.0 of the ECMWF’s (European Centre for Medium-Range Weather Forecasts) Re-Analysis (ERA-interim) and rainfall from the TRMM (Tropical Rainfall Measuring Mission) product. In this study, we first analyse the rainfall characteristics of two GNSS-IWV stations by comparing their complementary cumulative distribution function (CCDF). Second, we separately derive the relation between rainfall vs. CAPE and GNSS-IWV. Based on our distribution fitting analysis, we observe an exponential relation of rainfall to GNSS-IWV. In contrast, we report a power-law relationship between the daily mean value of rainfall and CAPE at the GNSS-IWV station locations in the eastern central Andes that is close to the theoretical relationship based on parcel theory. Third, we generate a joint regression model through a multivariable regression analysis using CAPE and GNSS-IWV to explain the contribution of both variables in the presence of each other to extreme rainfall during the austral summer season. We found that rainfall can be characterised with a higher statistical significance for higher rainfall quantiles, e.g., the 0.9 quantile based on goodness-of-fit criterion for quantile regression. We observed different contributions of CAPE and GNSS-IWV to rainfall for each station for the 0.9 quantile. Fourth, we identify the temporal relation between extreme rainfall (the 90th, 95th, and 99th percentiles) and both GNSS-IWV and CAPE at 6 h time steps. We observed an increase before the rainfall event and at the time of peak rainfall—both for GNSS-integrated water vapour and CAPE. We show higher values of CAPE and GNSS-IWV for higher rainfall percentiles (99th and 95th percentiles) compared to the 90th percentile at a 6-h temporal scale. Based on our correlation analyses and the dynamics of the time series, we show that both GNSS-IWV and CAPE had comparable magnitudes, and we argue to consider both climatic variables when investigating their effect on rainfall extremes.
Time series of groundwater and stream water quality often exhibit substantial temporal and spatial variability, whereas typical existing monitoring data sets, e.g. from environmental agencies, are usually characterized by relatively low sampling frequency and irregular sampling in space and/or time. This complicates the differentiation between anthropogenic influence and natural variability as well as the detection of changes in water quality which indicate changes in single drivers. We suggest the new term "dominant changes" for changes in multivariate water quality data which concern (1) multiple variables, (2) multiple sites and (3) long-term patterns and present an exploratory framework for the detection of such dominant changes in data sets with irregular sampling in space and time. Firstly, a non-linear dimension-reduction technique was used to summarize the dominant spatiotemporal dynamics in the multivariate water quality data set in a few components. Those were used to derive hypotheses on the dominant drivers influencing water quality. Secondly, different sampling sites were compared with respect to median component values. Thirdly, time series of the components at single sites were analysed for long-term patterns. We tested the approach with a joint stream water and groundwater data set quality consisting of 1572 samples, each comprising sixteen variables, sampled with a spatially and temporally irregular sampling scheme at 29 sites in northeast Germany from 1998 to 2009. The first four components were interpreted as (1) an agriculturally induced enhancement of the natural background level of solute concentration, (2) a redox sequence from reducing conditions in deep groundwater to post-oxic conditions in shallow groundwater and oxic conditions in stream water, (3) a mixing ratio of deep and shallow groundwater to the streamflow and (4) sporadic events of slurry application in the agricultural practice. Dominant changes were observed for the first two components. The changing intensity of the first component was interpreted as response to the temporal variability of the thickness of the unsaturated zone. A steady increase in the second component at most stream water sites pointed towards progressing depletion of the denitrification capacity of the deep aquifer.
Improving society's ability to prepare for, respond to and recover from flooding requires integrated, anticipatory flood risk management (FRM). However, most countries still focus their efforts on responding to flooding events if and when they occur rather than addressing their current and future vulnerability to flooding. Flood insurance is one mechanism that could promote a more ex ante approach to risk by supporting risk reduction activities. This paper uses an adapted version of Easton's System Theory to investigate the role of insurance for FRM in Germany and England. We introduce an anticipatory FRM framework, which allows flood insurance to be considered as part of a broader policy field. We analyze if and how flood insurance can catalyze a change toward a more anticipatory approach to FRM. In particular we consider insurance's role in influencing five key components of anticipatory FRM: risk knowledge, prevention through better planning, property‐level protection measures, structural protection and preparedness (for response). We find that in both countries FRM is still a reactive, event‐driven process, while anticipatory FRM remains underdeveloped. Collaboration between insurers and FRM decision‐makers has already been successful, for example in improving risk knowledge and awareness, while in other areas insurance acts as a disincentive for more risk reduction action. In both countries there is evidence that insurance can play a significant role in encouraging anticipatory FRM, but this remains underutilized. Effective collaboration between insurers and government should not be seen as a cost, but as an investment to secure future insurability through flood resilience.
Adaptation to flood risk
(2017)
As flood impacts are increasing in large parts of the world, understanding the primary drivers of changes in risk is essential for effective adaptation. To gain more knowledge on the basis of empirical case studies, we analyze eight paired floods, that is, consecutive flood events that occurred in the same region, with the second flood causing significantly lower damage. These success stories of risk reduction were selected across different socioeconomic and hydro-climatic contexts. The potential of societies to adapt is uncovered by describing triggered societal changes, as well as formal measures and spontaneous processes that reduced flood risk. This novel approach has the potential to build the basis for an international data collection and analysis effort to better understand and attribute changes in risk due to hydrological extremes in the framework of the IAHSs Panta Rhei initiative. Across all case studies, we find that lower damage caused by the second event was mainly due to significant reductions in vulnerability, for example, via raised risk awareness, preparedness, and improvements of organizational emergency management. Thus, vulnerability reduction plays an essential role for successful adaptation. Our work shows that there is a high potential to adapt, but there remains the challenge to stimulate measures that reduce vulnerability and risk in periods in which extreme events do not occur.
Although hydrologic models provide hypothesis testing of complex dynamics occurring at catchments, fresh-water quality modeling is still incipient at many subtropical headwaters. In Brazil, a few modeling studies assess freshwater nutrients, limiting policies on hydrologic ecosystem services. This paper aims to compare freshwater quality scenarios under different land-use and land-cover (LULC) change, one of them related to ecosystem-based adaptation (EbA), in Brazilian headwaters. Using the spatially semi-distributed Soil and Water Assessment Tool (SWAT) model, nitrate, total phosphorous (TP) and sediment were modeled in catchments ranging from 7.2 to 1037 km(2). These head-waters were eligible areas of the Brazilian payment for ecosystem services (PES) projects in the Cantareira water supply system, which had supplied water to 9 million people in the Sao Paulo metropolitan region (SPMR). We considered SWAT modeling of three LULC scenarios: (i) recent past scenario (S1), with historical LULC in 1990; (ii) current land-use scenario (S2), with LULC for the period 2010-2015 with field validation; and (iii) future land-use scenario with PES (S2 + EbA). This latter scenario proposed forest cover restoration through EbA following the river basin plan by 2035. These three LULC scenarios were tested with a selected record of rainfall and evapotranspiration observed in 2006-2014, with the occurrence of extreme droughts. To assess hydrologic services, we proposed the hydrologic service index (HSI), as a new composite metric comparing water pollution levels (WPL) for reference catchments, related to the grey water footprint (greyWF) and water yield. On the one hand, water quality simulations allowed for the regionalization of greyWF at spatial scales under LULC scenarios. According to the critical threshold, HSI identified areas as less or more sustainable catchments. On the other hand, conservation practices simulated through the S2 + EbA scenario envisaged not only additional and viable best management practices (BMP), but also preventive decision-making at the headwaters of water supply systems.
Evaluating climate geoengineering proposals in the context of the Paris Agreement temperature goals
(2018)
Current mitigation efforts and existing future commitments are inadequate to accomplish the Paris Agreement temperature goals. In light of this, research and debate are intensifying on the possibilities of additionally employing proposed climate geoengineering technologies, either through atmospheric carbon dioxide removal or farther-reaching interventions altering the Earth’s radiative energy budget. Although research indicates that several techniques may eventually have the physical potential to contribute to limiting climate change, all are in early stages of development, involve substantial uncertainties and risks, and raise ethical and governance dilemmas. Based on present knowledge, climate geoengineering techniques cannot be relied on to significantly contribute to meeting the Paris Agreement temperature goals.
The SISAL database
(2018)
Stable isotope records from speleothems provide information on past climate changes, most particularly information that can be used to reconstruct past changes in precipitation and atmospheric circulation. These records are increasingly being used to provide "out-of-sample" evaluations of isotope-enabled climate models. SISAL (Speleothem Isotope Synthesis and Analysis) is an international working group of the Past Global Changes (PAGES) project. The working group aims to provide a comprehensive compilation of speleothem isotope records for climate reconstruction and model evaluation. The SISAL database contains data for individual speleothems, grouped by cave system. Stable isotopes of oxygen and carbon (delta O-18, delta C-13) measurements are referenced by distance from the top or bottom of the speleothem. Additional tables provide information on dating, including information on the dates used to construct the original age model and sufficient information to assess the quality of each data set and to erect a standardized chronology across different speleothems. The metadata table provides location information, information on the full range of measurements carried out on each speleothem and information on the cave system that is relevant to the interpretation of the records, as well as citations for both publications and archived data.
Nearly 13,000 years ago, the warming trend into the Holocene was sharply interrupted by a reversal to near glacial conditions. Climatic causes and ecological consequences of the Younger Dryas (YD) have been extensively studied, however proxy archives from the Mediterranean basin capturing this period are scarce and do not provide annual resolution. Here, we report a hydroclimatic reconstruction from stable isotopes (delta O-18, delta C-13) in subfossil pines from southern France. Growing before and during the transition period into the YD (12 900-12 600 cal BP), the trees provide an annually resolved, continuous sequence of atmospheric change. Isotopic signature of tree sourcewater (delta O-18(sw)) and estimates of relative air humidity were reconstructed as a proxy for variations in air mass origin and precipitation regime. We find a distinct increase in inter-annual variability of sourcewater isotopes (delta O-18(sw)), with three major downturn phases of increasing magnitude beginning at 12 740 cal BP. The observed variation most likely results from an amplified intensity of North Atlantic (low delta O-18(sw)) versus Mediterranean (high delta O-18(sw)) precipitation. This marked pattern of climate variability is not seen in records from higher latitudes and is likely a consequence of atmospheric circulation oscillations at the margin of the southward moving polar front.
More than 41% of the Earth’s land area is covered by permanent or seasonally arid dryland ecosystems. Global development and human activity have led to an increase in aridity, resulting in ecosystem degradation and desertification around the world. The objective of the present work was to investigate and compare the microbial community structure and geochemical characteristics of two geographically distinct saline pan sediments in the Kalahari Desert of southern Africa. Our data suggest that these microbial communities have been shaped by geochemical drivers, including water content, salinity, and the supply of organic matter. Using Illumina 16S rRNA gene sequencing, this study provides new insights into the diversity of bacteria and archaea in semi-arid, saline, and low-carbon environments. Many of the observed taxa are halophilic and adapted to water-limiting conditions. The analysis reveals a high relative abundance of halophilic archaea (primarily Halobacteria), and the bacterial diversity is marked by an abundance of Gemmatimonadetes and spore-forming Firmicutes. In the deeper, anoxic layers, candidate division MSBL1, and acetogenic bacteria (Acetothermia) are abundant. Together, the taxonomic information and geochemical data suggest that acetogenesis could be a prevalent form of metabolism in the deep layers of a saline pan.
Near-surface geophysical imaging of alluvial fan settings is a challenging task but crucial for understating geological processes in such settings. The alluvial fan of Ghor Al-Haditha at the southeast shore of the Dead Sea is strongly affected by localized subsidence and destructive sinkhole collapses, with a significantly increasing sinkhole formation rate since ca. 1983. A similar increase is observed also on the western shore of the Dead Sea, in correlation with an ongoing decline in the Dead Sea level. Since different structural models of the upper 50 m of the alluvial fan and varying hypothetical sinkhole processes have been suggested for the Ghor Al-Haditha area in the past, this study aimed to clarify the subsurface characteristics responsible for sinkhole development. For this purpose, high-frequency shear wave reflection vibratory seismic surveys were carried out in the Ghor Al-Haditha area along several crossing and parallel profiles with a total length of 1.8 and 2.1 km in 2013 and 2014, respectively. The sedimentary architecture of the alluvial fan at Ghor Al-Haditha is resolved down to a depth of nearly 200 m at a high resolution and is calibrated with the stratigraphic profiles of two boreholes located inside the survey area. The most surprising result of the survey is the absence of evidence of a thick (>2-10 m) compacted salt layer formerly suggested to lie at ca. 35-40 m depth. Instead, seismic reflection amplitudes and velocities image with good continuity a complex interlocking of alluvial fan deposits and lacustrine sediments of the Dead Sea between 0 and 200 m depth. Furthermore, the underground section of areas affected by sinkholes is characterized by highly scattering wave fields and reduced seismic interval velocities. We propose that the Dead Sea mud layers, which comprise distributed inclusions or lenses of evaporitic chloride, sulfate, and carbonate minerals as well as clay silicates, become increasingly exposed to unsaturated water as the sea level declines and are consequently destabilized and mobilized by both dissolution and physical erosion in the subsurface. This new interpretation of the underlying cause of sinkhole development is supported by surface observations in nearby channel systems. Overall, this study shows that shear wave seismic reflection technique is a promising method for enhanced near-surface imaging in such challenging alluvial fan settings.
The location and magnitude of Himalayan tectonic activity has been debated for decades, and several aspects remain unknown. For instance, the spatial distribution of crustal shortening that ultimately sustains Himalayan topography and the activity of major fault zones remain unknown at Ma timescales. In this study, we address the spatial deformation pattern in the data-scarce western Himalaya. We calculated catchment averaged, normalized river-steepness indices of non-glaciated drainage basins with tributary catchment areas between 5 and 200 km(2) (n = 2138). We analyzed the spatial distribution of the relative change of river steepness both along and across strike to gain information about the regional distribution of differential uplift pattern and relate this to the activity of distinctive fault segments. For our study area, we observe a positive correlation of averaged k(sn) values with long-term exhumation rates derived from previously published thermochronologic datasets combined with thermal modeling as well as with millennial timescale denudation rates based on cosmogenic nuclide dating. Our results indicate three tectono-geomorphic segments with distinctive landscape morphology, structural architecture, and fault geometry along the western Himalaya: Garhwal-Sutlej, Chamba, and Kashmir Himalaya (from east to west). Moreover, our data recognize distinctive fault segments showing varying thrust activity along strike of the Main Frontal Thrust, the Main Boundary Thrust, and in the vicinity of the steep topographic transition between the Lesser and Greater Himalaya. In this region, we relate out-of-sequence deformation along major basement thrust ramps, such as the Munsiari Thrust with deformation along a mid-crustal ramp along the basal decollement. We suggest that during the Quaternary, all major fault zones in the Western Himalaya experienced out-of-sequence faulting and have accommodated some portion of crustal shortening.
Rotational motions play a key role in measuring seismic wavefield properties. Using newly developed portable rotational instruments, it is now possible to directly measure rotational motions in a broad frequency range. Here, we investigated the instrumental self-noise and data quality in a huddle test in Fürstenfeldbruck, Germany, in August 2019. We compare the data from six rotational and three translational sensors. We studied the recorded signals using correlation, coherence analysis, and probabilistic power spectral densities. We sorted the coherent noise into five groups with respect to the similarities in frequency content and shape of the signals. These coherent noises were most likely caused by electrical devices, the dehumidifier system in the building, humans, and natural sources such as wind. We calculated self-noise levels through probabilistic power spectral densities and by applying the Sleeman method, a three-sensor method. Our results from both methods indicate that self-noise levels are stable between 0.5 and 40 Hz. Furthermore, we recorded the 29 August 2019 ML 3.4 Dettingen earthquake. The calculated source directions are found to be realistic for all sensors in comparison to the real back azimuth. We conclude that the five tested blueSeis-3A rotational sensors, when compared with respect to coherent noise, self-noise, and source direction, provide reliable and consistent results. Hence, field experiments with single rotational sensors can be undertaken.
Models for the predictions of monetary losses from floods mainly blend data deemed to represent a single flood type and region. Moreover, these approaches largely ignore indicators of preparedness and how predictors may vary between regions and events, challenging the transferability of flood loss models. We use a flood loss database of 1812 German flood-affected households to explore how Bayesian multilevel models can estimate normalised flood damage stratified by event, region, or flood process type. Multilevel models acknowledge natural groups in the data and allow each group to learn from others. We obtain posterior estimates that differ between flood types, with credibly varying influences of water depth, contamination, duration, implementation of property-level precautionary measures, insurance, and previous flood experience; these influences overlap across most events or regions, however. We infer that the underlying damaging processes of distinct flood types deserve further attention. Each reported flood loss and affected region involved mixed flood types, likely explaining the uncertainty in the coefficients. Our results emphasise the need to consider flood types as an important step towards applying flood loss models elsewhere. We argue that failing to do so may unduly generalise the model and systematically bias loss estimations from empirical data.
It was lower during the monsoon due to rainfall, which reduces open burning activities within the valley and in the surrounding regions and thus reduces sources of CO. The meteorology of the valley also played a key role in determining the CO mixing ratios. The wind is calm and easterly in the shallow mixing layer, with a mixing layer height (MLH) of about 250 m, during the night and early morning. The MLH slowly increases after sunrise and decreases in the afternoon. As a result, the westerly wind becomes active and reduces the mixing ratio during the daytime. Furthermore, there was evidence of an increase in the O-3 mixing ratios in the Kathmandu Valley as a result of emissions in the Indo-Gangetic Plain (IGP) region, particularly from biomass burning including agroresidue burning. A top-down estimate of the CO emission flux was made by using the CO mixing ratio and mixing layer height measured at Bode. The estimated annual CO flux at Bode was 4.9 mu g M-2 s(-1), which is 2-14 times higher than that in widely used emission inventory databases (EDGAR HTAP, REAS and INTEX-B). This difference in CO flux between Bode and other emission databases likely arises from large uncertainties in both the top-down and bottom-up approaches to estimating the emission flux. The O-3 mixing ratio was found to be highest during the premonsoon season at all sites, while the timing of the seasonal minimum varied across the sites. The daily maximum 8 h average O-3 exceeded the WHO recommended guideline of 50 ppb on more days at the hilltop station of Nagarkot (159 out of 357 days) than at the urban valley bottom sites of Paknajol (132 out of 354 days) and Bode (102 out of 353 days), presumably due to the influence of free-tropospheric air at the high-altitude site (as also indicated by Putero et al., 2015, for the Paknajol site in the Kathmandu Valley) as well as to titration of O-3 by fresh NOx emissions near the urban sites. More than 78 % of the exceedance days were during the premonsoon period at all sites. The high O-3 mixing ratio observed during the premonsoon period is of a concern for human health and ecosystems, including agroecosystems in the Kathmandu Valley and surrounding regions.
Organic matter characteristics in yedoma and thermokarst deposits on Baldwin Peninsula, west Alaska
(2018)
As Arctic warming continues and permafrost thaws, more soil and sedimentary organic matter (OM) will be decomposed in northern high latitudes. Still, uncertainties remain in the quality of the OM and the size of the organic carbon (OC) pools stored in different deposit types of permafrost landscapes. This study presents OM data from deep permafrost and lake deposits on the Baldwin Peninsula which is located in the southern portion of the continuous permafrost zone in west Alaska. Sediment samples from yedoma and drained thermokarst lake basin (DTLB) deposits as well as thermokarst lake sediments were analyzed for cryostratigraphical and biogeochemical parameters and their lipid biomarker composition to identify the below-ground OC pool size and OM quality of ice-rich permafrost on the Baldwin Peninsula. We provide the first detailed characterization of yedoma deposits on Baldwin Peninsula. We show that three-quarters of soil OC in the frozen deposits of the study region (total of 68 Mt) is stored in DTLB deposits (52 Mt) and one-quarter in the frozen yedoma deposits (16 Mt). The lake sediments contain a relatively small OC pool (4 Mt), but have the highest volumetric OC content (93 kgm(-3)) compared to the DTLB (35 kgm(-3)) and yedoma deposits (8 kgm(-3)), largely due to differences in the ground ice content. The biomarker analysis indicates that the OM in both yedoma and DTLB deposits is mainly of terrestrial origin. Nevertheless, the relatively high carbon preference index of plant leaf waxes in combination with a lack of a degradation trend with depth in the yedoma deposits indi-cates that OM stored in yedoma is less degraded than that stored in DTLB deposits. This suggests that OM in yedoma has a higher potential for decomposition upon thaw, despite the relatively small size of this pool. These findings show that the use of lipid biomarker analysis is valuable in the assessment of the potential future greenhouse gas emissions from thawing permafrost, especially because this area, close to the discontinuous permafrost boundary, is projected to thaw substantially within the 21st century.
A tale of shifting relations
(2021)
Understanding the dynamics between the East Asian summer (EASM) and winter monsoon (EAWM) is needed to predict their variability under future global warming scenarios. Here, we investigate the relationship between EASM and EAWM as well as the mechanisms driving their variability during the last 10,000 years by stacking marine and terrestrial (non-speleothem) proxy records from the East Asian realm. This provides a regional and proxy independent signal for both monsoonal systems. The respective signal was subsequently analysed using a linear regression model. We find that the phase relationship between EASM and EAWM is not time-constant and significantly depends on orbital configuration changes. In addition, changes in the Atlantic Meridional Overturning circulation, Arctic sea-ice coverage, El Niño-Southern Oscillation and Sun Spot numbers contributed to millennial scale changes in the EASM and EAWM during the Holocene. We also argue that the bulk signal of monsoonal activity captured by the stacked non-speleothem proxy records supports the previously argued bias of speleothem climatic archives to moisture source changes and/or seasonality.
Movements of tectonic plates often induce oblique deformation at divergent plate boundaries. This is in striking contrast with traditional conceptual models of rifting and rifted margin formation, which often assume 2-D deformation where the rift velocity is oriented perpendicular to the plate boundary. Here we quantify the validity of this assumption by analysing the kinematics of major continent-scale rift systems in a global plate tectonic reconstruction from the onset of Pangea breakup until the present day. We evaluate rift obliquity by joint examination of relative extension velocity and local rift trend using the script-based plate reconstruction software pyGPlates. Our results show that the global mean rift obliquity since 230 Ma amounts to 34 degrees with a standard deviation of 24 degrees, using the convention that the angle of obliquity is spanned by extension direction and rift trend normal. We find that more than similar to 70 % of all rift segments exceeded an obliquity of 20 degrees demonstrating that oblique rifting should be considered the rule, not the exception. In many cases, rift obliquity and extension velocity increase during rift evolution (e.g. Australia-Antarctica, Gulf of California, South Atlantic, India-Antarctica), which suggests an underlying geodynamic correlation via obliquity-dependent rift strength. Oblique rifting produces 3-D stress and strain fields that cannot be accounted for in simplified 2-D plane strain analysis. We therefore highlight the importance of 3-D approaches in modelling, surveying, and interpretation of most rift segments on Earth where oblique rifting is the dominant mode of deformation.
High-latitude treeless ecosystems represent spatially highly heterogeneous landscapes with small net carbon fluxes and a short growing season. Reliable observations and process understanding are critical for projections of the carbon balance of the climate-sensitive tundra. Space-borne remote sensing is the only tool to obtain spatially continuous and temporally resolved information on vegetation greenness and activity in remote circumpolar areas. However, confounding effects from persistent clouds, low sun elevation angles, numerous lakes, widespread surface inundation, and the sparseness of the vegetation render it highly challenging. Here, we conduct an extensive analysis of the timing of peak vegetation productivity as shown by satellite observations of complementary indicators of plant greenness and photosynthesis. We choose to focus on productivity during the peak of the growing season, as it importantly affects the total annual carbon uptake. The suite of indicators are as follows: (1) MODIS-based vegetation indices (VIs) as proxies for the fraction of incident photosynthetically active radiation (PAR) that is absorbed (fPAR), (2) VIs combined with estimates of PAR as a proxy of the total absorbed radiation (APAR), (3) sun-induced chlorophyll fluorescence (SIF) serving as a proxy for photosynthesis, (4) vegetation optical depth (VOD), indicative of total water content and (5) empirically upscaled modelled gross primary productivity (GPP). Averaged over the pan-Arctic we find a clear order of the annual peak as APAR <= GPP < SIF < VIs/VOD. SIF as an indicator of photosynthesis is maximised around the time of highest annual temperatures. The modelled GPP peaks at a similar time to APAR. The time lag of the annual peak between APAR and instantaneous SIF fluxes indicates that the SIF data do contain information on light-use efficiency of tundra vegetation, but further detailed studies are necessary to verify this. Delayed peak greenness compared to peak photosynthesis is consistently found across years and land-cover classes. A particularly late peak of the normalised difference vegetation index (NDVI) in regions with very small seasonality in greenness and a high amount of lakes probably originates from artefacts. Given the very short growing season in circumpolar areas, the average time difference in maximum annual photosynthetic activity and greenness or growth of 3 to 25 days (depending on the data sets chosen) is important and needs to be considered when using satellite observations as drivers in vegetation models.
Geochemical processes such as mineral dissolution and precipitation alter the microstructure of rocks, and thereby affect their hydraulic and mechanical behaviour. Quantifying these property changes and considering them in reservoir simulations is essential for a sustainable utilisation of the geological subsurface. Due to the lack of alternatives, analytical methods and empirical relations are currently applied to estimate evolving hydraulic and mechanical rock properties associated with chemical reactions. However, the predictive capabilities of analytical approaches remain limited, since they assume idealised microstructures, and thus are not able to reflect property evolution for dynamic processes. Hence, aim of the present thesis is to improve the prediction of permeability and stiffness changes resulting from pore space alterations of reservoir sandstones.
A detailed representation of rock microstructure, including the morphology and connectivity of pores, is essential to accurately determine physical rock properties. For that purpose, three-dimensional pore-scale models of typical reservoir sandstones, obtained from highly resolved micro-computed tomography (micro-CT), are used to numerically calculate permeability and stiffness. In order to adequately depict characteristic distributions of secondary minerals, the virtual samples are systematically altered and resulting trends among the geometric, hydraulic, and mechanical rock properties are quantified. It is demonstrated that the geochemical reaction regime controls the location of mineral precipitation within the pore space, and thereby crucially affects the permeability evolution. This emphasises the requirement of determining distinctive porosity-permeability relationships
by means of digital pore-scale models. By contrast, a substantial impact of spatial alterations patterns on the stiffness evolution of reservoir sandstones are only observed in case of certain microstructures, such as highly porous granular rocks or sandstones comprising framework-supporting cementations. In order to construct synthetic granular samples a process-based approach is proposed including grain deposition and diagenetic cementation. It is demonstrated that the generated samples reliably represent the microstructural complexity of natural sandstones. Thereby, general limitations of imaging techniques can be overcome and various realisations of granular rocks can be flexibly produced. These can be further altered by virtual experiments, offering a fast and cost-effective way to examine the impact of precipitation, dissolution or fracturing on various petrophysical correlations.
The presented research work provides methodological principles to quantify trends in permeability and stiffness resulting from geochemical processes. The calculated physical property relations are directly linked to pore-scale alterations, and thus have a higher accuracy than commonly applied analytical approaches. This will considerably improve the predictive capabilities of reservoir models, and is further relevant to assess and reduce potential risks, such as productivity or injectivity losses as well as reservoir compaction or fault reactivation. Hence, the proposed method is of paramount importance for a wide range of natural and engineered subsurface applications, including geothermal energy systems, hydrocarbon reservoirs, CO2 and energy storage as well as hydrothermal deposit exploration.
One common approach to cope with floods is the implementation of structural flood protection measures, such as levees or flood-control reservoirs, which substantially reduce the probability of flooding at the time of implementation. Numerous scholars have problematized this approach. They have shown that increasing the levels of flood protection can attract more settlements and high-value assets in the areas protected by the new measures. Other studies have explored how structural measures can generate a sense of complacency, which can act to reduce preparedness. These paradoxical risk changes have been described as "levee effect", "safe development paradox" or "safety dilemma". In this commentary, we briefly review this phenomenon by critically analysing the intended benefits and unintended effects of structural flood protection, and then we propose an interdisciplinary research agenda to uncover these paradoxical dynamics of risk.
Mechanical and/or chemical removal of material from the subsurface may generate large subsurface cavities, the destabilisation of which can lead to ground collapse and the formation of sinkholes. Numerical simulation of the interaction of cavity growth, host material deformation and overburden collapse is desirable to better understand the sinkhole hazard but is a challenging task due to the involved high strains and material discontinuities. Here, we present 2-D distinct element method numerical simulations of cavity growth and sinkhole development. Firstly, we simulate cavity formation by quasi-static, stepwise removal of material in a single growing zone of an arbitrary geometry and depth. We benchmark this approach against analytical and boundary element method models of a deep void space in a linear elastic material. Secondly, we explore the effects of properties of different uniform materials on cavity stability and sinkhole development. We perform simulated biaxial tests to calibrate macroscopic geotechnical parameters of three model materials representative of those in which sinkholes develop at the Dead Sea shoreline: mud, alluvium and salt. We show that weak materials do not support large cavities, leading to gradual sagging or suffusion-style subsidence. Strong materials support quasi-stable to stable cavities, the overburdens of which may fail suddenly in a caprock or bedrock collapse style. Thirdly, we examine the consequences of layered arrangements of weak and strong materials. We find that these are more susceptible to sinkhole collapse than uniform materials not only due to a lower integrated strength of the overburden but also due to an inhibition of stabilising stress arching. Finally, we compare our model sinkhole geometries to observations at the Ghor Al-Haditha sinkhole site in Jordan. Sinkhole depth / diameter ratios of 0.15 in mud, 0.37 in alluvium and 0.33 in salt are reproduced successfully in the calibrated model materials. The model results suggest that the observed distribution of sinkhole depth / diameter values in each material type may partly reflect sinkhole growth trends.
The phi(ev) is calculated from high-resolution discharge and precipitation data for several rain events with a cumulative precipitation P-cum ranging from less than 5mm to more than 80 mm. Because of the high uncertainty of phi(ev) associated with the hydrograph separation method, phi(ev) is calculated with several methods, including graphical methods, digital filters and a tracer-based method. The results indicate that the hydrological response depends on (theta) over bar (ini): during dry conditions phi(ev) is consistently below 0.1, even for events with high and intense precipitation. Above a threshold of (theta) over bar (ini) = 34 vol % phi(ev) can reach values up to 0.99 but there is a high scatter. Some variability can be explained with a weak correlation of phi(ev) with P-cum and rain intensity, but a considerable part of the variability remains unexplained. It is concluded that threshold-based methods can be helpful to prevent overestimation of the hydrological response during dry catchment conditions. The impact of soil moisture on the hydrological response during wet catchment conditions, however, is still insufficiently understood and cannot be generalized based on the present results.
Coastal ecosystems in the Arctic are affected by climate change. As summer rainfall frequency and intensity are projected to increase in the future, more organic matter, nutrients and sediment could bemobilized and transported into the coastal nearshore zones. However, knowledge of current processes and future changes is limited. We investigated streamflow dynamics and the impacts of summer rainfall on lateral fluxes in a small coastal catchment on Herschel Island in the western Canadian Arctic. For the summer monitoring periods of 2014-2016, mean dissolved organic matter flux over 17 days amounted to 82.7 +/- 30.7 kg km(-2) and mean total dissolved solids flux to 5252 +/- 1224 kg km(-2). Flux of suspended sediment was 7245 kg km(-2) in 2015, and 369 kg km(-2) in 2016. We found that 2.0% of suspended sediment was composed of particulate organic carbon. Data and hysteresis analysis suggest a limited supply of sediments; their interannual variability is most likely caused by short-lived localized disturbances. In contrast, our results imply that dissolved organic carbon is widely available throughout the catchment and exhibits positive linear relationship with runoff. We hypothesize that increased projected rainfall in the future will result in a similar increase of dissolved organic carbon fluxes.
Alluvial and transport-limited bedrock rivers constitute the majority of fluvial systems on Earth. Their long profiles hold clues to their present state and past evolution. We currently possess first-principles-based governing equations for flow, sediment transport, and channel morphodynamics in these systems, which we lack for detachment-limited bedrock rivers. Here we formally couple these equations for transport-limited gravel-bed river long-profile evolution. The result is a new predictive relationship whose functional form and parameters are grounded in theory and defined through experimental data. From this, we produce a power-law analytical solution and a finite-difference numerical solution to long-profile evolution. Steady-state channel concavity and steepness are diagnostic of external drivers: concavity decreases with increasing uplift rate, and steepness increases with an increasing sediment-to-water supply ratio. Constraining free parameters explains common observations of river form: to match observed channel concavities, gravel-sized sediments must weather and fine - typically rapidly - and valleys typically should widen gradually. To match the empirical square-root width-discharge scaling in equilibrium-width gravel-bed rivers, downstream fining must occur. The ability to assign a cause to such observations is the direct result of a deductive approach to developing equations for landscape evolution.
The thermal structure of subduction zones exerts a major influence on deep-seated mechanical and chemical processes controlling arc magmatism, seismicity, and global element cycles. Accretionary complexes exposed inland may comprise tectonic blocks with contrasting pressure-temperature (P-T) histories, making it possible to investigate the dynamics and thermal evolution of former subduction interfaces. With this aim, we present new Lu-Hf geochronological results for mafic rocks of the Halilbagi Complex (Anatolia) that evolved along different thermal gradients. Samples include a lawsonite-epidote blueschist, a lawsonite-epidote eclogite, and an epidote eclogite (all with counter-clockwise P-T paths), a prograde lawsonite blueschist with a "hairpin"-type P-T path, and a garnet amphibolite from the overlying sub-ophiolitic metamorphic sole. Equilibrium phase diagrams suggest that the garnet amphibolite formed at similar to 0.6-0.7 GPa and 800-850 degrees C, whereas the prograde lawsonite blueschist records burial from 2.1 GPa and 420 degrees C to 2.6 GPa and 520 degrees C. Well-defined Lu-Hf isochrons were obtained for the epidote eclogite (92.38 +/- 0.22 Ma) and the lawsonite-epidote blueschist (90.19 +/- 0.54 Ma), suggesting rapid garnet growth. The lawsonite-epidote eclogite (87.30 +/- 0.39 Ma) and the prograde lawsonite blueschist (ca. 86 Ma) are younger, whereas the garnet amphibolite (104.5 +/- 3.5 Ma) is older. Our data reveal a consistent trend of progressively decreasing geothermal gradient from granulite-facies conditions at similar to 104 Ma to the epidote-eclogite facies around 92 Ma, and the lawsonite blueschist-facies between 90 Ma and 86 Ma. Three Lu-Hf garnet dates (between 92 Ma and 87 Ma) weighted toward the growth of post-peak rims (as indicated by Lu distribution in garnet) suggest that the HP/LT rocks were exhumed continuously and not episodically. We infer that HP/LT metamorphic rocks within the Halilbagi Complex were subjected to continuous return flow, with "warm" rocks being exhumed during the tectonic burial of "cold" ones. Our results, combined with regional geological constraints, allow us to speculate that subduction started at a transform fault near a mid-oceanic spreading centre. Following its formation, this ancient subduction interface evolved thermally over more than 15 Myr, most likely as a result of heat dissipation rather than crustal underplating. (C) 2018, China University of Geosciences (Beijing) and Peking University. Production and hosting by Elsevier B.V.
SiO(2 )is the main component of silicate melts and thus controls their network structure and physical properties. The compressibility and viscosities of melts at depth are governed by their short range atomic and electronic structure. We measured the O K-edge and the Si L-2,L-3-edge in silica up to 110 GPa using X-ray Raman scattering spectroscopy, and found a striking match to calculated spectra based on structures from molecular dynamic simulations. Between 20 and 27 GPa, Si-[4] species are converted into a mixture of Si-[5] and Si-[6] species and between 60 and 70 GPa, Si-[6] becomes dominant at the expense of Si-[5] with no further increase up to at least 110 GPa. Coordination higher than 6 is only reached beyond 140 GPa, corroborating results from Brillouin scattering. Network modifying elements in silicate melts may shift this change in coordination to lower pressures and thus magmas could be denser than residual solids at the depth of the core-mantle boundary.
In active mountain belts with steep terrain, bedrock landsliding is a major erosional agent. In the Himalayas, landsliding is driven by annual hydro-meteorological forcing due to the summer monsoon and by rarer, exceptional events, such as earthquakes. Independent methods yield erosion rate estimates that appear to increase with sampling time, suggesting that rare, high-magnitude erosion events dominate the erosional budget. Nevertheless, until now, neither the contribution of monsoon and earthquakes to landslide erosion nor the proportion of erosion due to rare, giant landslides have been quantified in the Himalayas. We address these challenges by combining and analysing earthquake- and monsoon-induced landslide inventories across different timescales. With time series of 5 m satellite images over four main valleys in central Nepal, we comprehensively mapped landslides caused by the monsoon from 2010 to 2018. We found no clear correlation between monsoon properties and landsliding and a similar mean landsliding rate for all valleys, except in 2015, where the valleys affected by the earthquake featured similar to 5-8 times more landsliding than the pre-earthquake mean rate. The longterm size-frequency distribution of monsoon-induced landsliding (MIL) was derived from these inventories and from an inventory of landslides larger than similar to 0.1 km(2) that occurred between 1972 and 2014. Using a published landslide inventory for the Gorkha 2015 earthquake, we derive the size-frequency distribution for earthquakeinduced landsliding (EQIL). These two distributions are dominated by infrequent, large and giant landslides but under-predict an estimated Holocene frequency of giant landslides (> 1 km(3)) which we derived from a literature compilation. This discrepancy can be resolved when modelling the effect of a full distribution of earthquakes of variable magnitude and when considering that a shallower earthquake may cause larger landslides. In this case, EQIL and MIL contribute about equally to a total long-term erosion of similar to 2 +/- 0.75 mm yr(-1) in agreement with most thermo-chronological data. Independently of the specific total and relative erosion rates, the heavy-tailed size-frequency distribution from MIL and EQIL and the very large maximal landslide size in the Himalayas indicate that mean landslide erosion rates increase with sampling time, as has been observed for independent erosion estimates. Further, we find that the sampling timescale required to adequately capture the frequency of the largest landslides, which is necessary for deriving long-term mean erosion rates, is often much longer than the averaging time of cosmogenic Be-10 methods. This observation presents a strong caveat when interpreting spatial or temporal variability in erosion rates from this method. Thus, in areas where a very large, rare landslide contributes heavily to long-term erosion (as the Himalayas), we recommend Be-10 sample in catchments with source areas > 10 000 km(2) to reduce the method mean bias to below similar to 20 % of the long-term erosion.
Shear-waves are the most energetic body-waves radiated from an earthquake, and are responsible for the destruction of engineered structures. In both short-term emergency response and long-term risk forecasting of disaster-resilient built environment, it is critical to predict spatially accurate distribution of shear-wave amplitudes. Although decades’ old theory proposes a deterministic, highly anisotropic, four-lobed shear-wave radiation pattern, from lack of convincing evidence, most empirical ground-shaking prediction models settled for an oversimplified stochastic radiation pattern that is isotropic on average. Today, using the large datasets of uniformly processed seismograms from several strike, normal, reverse, and oblique-slip earthquakes across the globe, compiled specifically for engineering applications, we could reveal, quantify, and calibrate the frequency-, distance-, and style-of-faulting dependent transition of shear-wave radiation between a stochastic-isotropic and a deterministic-anisotropic phenomenon. Consequent recalibration of empirical ground-shaking models dramatically improved their predictions: with isodistant anisotropic variations of ±40%, and 8% reduction in uncertainty. The outcomes presented here can potentially trigger a reappraisal of several practical issues in engineering seismology, particularly in seismic ground-shaking studies and seismic hazard and risk assessment.
A comprehensive description of electromagnetic processes related to equatorial plasma depletions (EPDs) is essential for understanding their evolution and day-to-day variability. Recently, field-aligned currents (FACs) flowing at both western and eastern edges of EPDs were observed to be interhemispheric rather than anti-parallel about the dip equator, as suggested by previous theoretical studies. In this paper, we investigate the spatial and temporal variability of the FACs orientation using simultaneous measurements of electron density and magnetic field gathered by ESA’s Swarm constellation mission. By using empirical models, we assess the role of the Pedersen conductance in the preference of the FACs to close either in the northern or southern magnetic hemisphere. Here we show that the closure of the FACs agrees with an electrostatic regime determined by a hemispherical asymmetry of the Pedersen conductance. That is, the EPD-related FACs close at lower altitudes in the hemisphere with the highest conductivity. The evidence of this conclusion stands on the general agreement between the longitudinal and seasonal variability of both the conductivity and the FACs orientation.