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 (10)
- 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)
How biased are our models?
(2021)
Geophysical process simulations play a crucial role in the understanding of the subsurface. This understanding is required to provide, for instance, clean energy sources such as geothermal energy. However, the calibration and validation of the physical models heavily rely on state measurements such as temperature. In this work, we demonstrate that focusing analyses purely on measurements introduces a high bias. This is illustrated through global sensitivity studies. The extensive exploration of the parameter space becomes feasible through the construction of suitable surrogate models via the reduced basis method, where the bias is found to result from very unequal data distribution. We propose schemes to compensate for parts of this bias. However, the bias cannot be entirely compensated. Therefore, we demonstrate the consequences of this bias with the example of a model calibration.
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.
How insoluble inclusions and intersecting layers affect the leaching process within potash seams
(2021)
Potash seams are a valuable resource containing several economically interesting, but also highly soluble minerals. In the presence of water, uncontrolled leaching can occur, endangering subsurface mining operations. In the present study, the influence of insoluble inclusions and intersecting layers on leaching zone evolution was examined by means of a reactive transport model. For that purpose, a scenario analysis was carried out, considering different rock distributions within a carnallite-bearing potash seam. The results show that reaction-dominated systems are not affected by heterogeneities at all, whereas transport-dominated systems exhibit a faster advance in homogeneous rock compositions. In return, the ratio of permeated rock in vertical direction is higher in heterogeneous systems. Literature data indicate that most natural potash systems are transport-dominated. Accordingly, insoluble inclusions and intersecting layers can usually be seen as beneficial with regard to reducing hazard potential as long as the mechanical stability of leaching zones is maintained. Thereby, the distribution of insoluble areas is of minor impact unless an inclined, intersecting layer occurs that accelerates leaching zone growth in one direction. Moreover, it is found that the saturation dependency of dissolution rates increases the growth rate in the long term, and therefore must be considered in risk assessments.
Hungry cities: how local food self-sufficiency relates to climate change, diets, and urbanisation
(2019)
Using a newly developed model approach and combining it with remote sensing, population, and climate data, first insights are provided into how local diets, urbanisation, and climate change relates to local urban food self-sufficiency. In plain terms, by utilizing the global peri-urban (PU) food production potential approximately lbn urban residents (30% of global urban population) can be locally nourished, whereby further urbanisation is by far the largest pressure factor on PU agriculture, followed by a change of diets, and climate change. A simple global food transport model which optimizes transport and neglects differences in local emission intensities indicates that CO2 emissions related to food transport can be reduced by a factor of 10.
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.
Ore precipitation in porphyry copper systems is generally characterized by metal zoning (Cu-Mo to Zn-Pb-Ag), which is suggested to be variably related to solubility decreases during fluid cooling, fluid-rock interactions, partitioning during fluid phase separation and mixing with external fluids. Here, we present new advances of a numerical process model by considering published constraints on the temperature- and salinity-dependent solubility of Cu, Pb and Zn in the ore fluid. We quantitatively investigate the roles of vapor-brine separation, halite saturation, initial metal contents, fluid mixing and remobilization as first-order controls of the physical hydrology on ore formation. The results show that the magmatic vapor and brine phases ascend with different residence times but as miscible fluid mixtures, with salinity increases generating metal-undersaturated bulk fluids. The release rates of magmatic fluids affect the location of the thermohaline fronts, leading to contrasting mechanisms for ore precipitation: higher rates result in halite saturation without significant metal zoning, lower rates produce zoned ore shells due to mixing with meteoric water. Varying metal contents can affect the order of the final metal precipitation sequence. Redissolution of precipitated metals results in zoned ore shell patterns in more peripheral locations and also decouples halite saturation from ore precipitation.
The Argentine-German Geodetic Observatory (AGGO) is one of the very few sites in the Southern Hemisphere equipped with comprehensive cutting-edge geodetic instrumentation. The employed observation techniques are used for a wide range of geophysical applications. The data set provides gravity time series and selected gravity models together with the hydrometeorological monitoring data of the observatory. These parameters are of great interest to the scientific community, e.g. for achieving accurate realization of terrestrial and celestial reference frames. Moreover, the availability of the hydrometeorological products is beneficial to inhabitants of the region as they allow for monitoring of environmental changes and natural hazards including extreme events. The hydrological data set is composed of time series of groundwater level, modelled and observed soil moisture content, soil temperature, and physical soil properties and aquifer properties. The meteorological time series include air temperature, humidity, pressure, wind speed, solar radiation, precipitation, and derived reference evapotranspiration. These data products are extended by gravity models of hydrological, oceanic, La Plata estuary, and atmospheric effects. The quality of the provided meteorological time series is tested via comparison to the two closest WMO (World Meteorological Organization) sites where data are available only in an inferior temporal resolution. The hydrological series are validated by comparing the respective forward-modelled gravity effects to independent gravity observations reduced up to a signal corresponding to local water storage variation. Most of the time series cover the time span between April 2016 and November 2018 with either no or only few missing data points. The data set is available at https://doi.org/10.588/GFZ.5.4.2018.001 (Mikolaj et al., 2018).
The new in situ geodynamic laboratory established in the framework of the ICDP Eger project aims to develop the most modern, comprehensive, multiparameter laboratory at depth for studying earthquake swarms, crustal fluid flow, mantle-derived CO2 and helium degassing, and processes of the deep biosphere. In order to reach a new level of high-frequency, near-source and multiparameter observation of earthquake swarms and related phenomena, such a laboratory comprises a set of shallow boreholes with high-frequency 3-D seismic arrays as well as modern continuous real-time fluid monitoring at depth and the study of the deep biosphere.
This laboratory is located in the western part of the Eger Rift at the border of the Czech Republic and Germany (in the West Bohemia–Vogtland geodynamic region) and comprises a set of five boreholes around the seismoactive zone. To date, all monitoring boreholes have been drilled. This includes the seismic monitoring boreholes S1, S2 and S3 in the crystalline units north and east of the major Nový Kostel seismogenic zone, borehole F3 in the Hartoušov mofette field and borehole S4 in the newly discovered Bažina maar near Libá. Supplementary borehole P1 is being prepared in the Neualbenreuth maar for paleoclimate and biological research. At each of these sites, a borehole broadband seismometer will be installed, and sites S1, S2 and S3 will also host a 3-D seismic array composed of a vertical geophone chain and surface seismic array. Seismic instrumenting has been completed in the S1 borehole and is in preparation in the remaining four monitoring boreholes. The continuous fluid monitoring site of Hartoušov includes three boreholes, F1, F2 and F3, and a pilot monitoring phase is underway. The laboratory also enables one to analyze microbial activity at CO2 mofettes and maar structures in the context of changes in habitats. The drillings into the maar volcanoes contribute to a better understanding of the Quaternary paleoclimate and volcanic activity.
The new in situ geodynamic laboratory established in the framework of the ICDP Eger project aims to develop the most modern, comprehensive, multiparameter laboratory at depth for studying earthquake swarms, crustal fluid flow, mantle-derived CO2 and helium degassing, and processes of the deep biosphere. In order to reach a new level of high-frequency, near-source and multiparameter observation of earthquake swarms and related phenomena, such a laboratory comprises a set of shallow boreholes with high-frequency 3-D seismic arrays as well as modern continuous real-time fluid monitoring at depth and the study of the deep biosphere.
This laboratory is located in the western part of the Eger Rift at the border of the Czech Republic and Germany (in the West Bohemia-Vogtland geodynamic region) and comprises a set of five boreholes around the seismoactive zone. To date, all monitoring boreholes have been drilled. This includes the seismic monitoring boreholes S1, S2 and S3 in the crystalline units north and east of the major Novy Kostel seismogenic zone, borehole F3 in the Hartousov mofette field and borehole S4 in the newly discovered Bazina maar near Liba. Supplementary borehole P1 is being prepared in the Neualbenreuth maar for paleoclimate and biological research. At each of these sites, a borehole broadband seismometer will be installed, and sites S1, S2 and S3 will also host a 3-D seismic array composed of a vertical geophone chain and surface seismic array. Seismic instrumenting has been completed in the S1 borehole and is in preparation in the remaining four monitoring boreholes. The continuous fluid monitoring site of Hartousov includes three boreholes, F1, F2 and F3, and a pilot monitoring phase is underway. The laboratory also enables one to analyze microbial activity at CO2 mofettes and maar structures in the context of changes in habitats. The drillings into the maar volcanoes contribute to a better understanding of the Quaternary paleoclimate and volcanic activity.
Through its magnetic activity, the Sun governs the conditions in Earth's vicinity, creating space weather events, which have drastic effects on our space- and ground-based technology.
One of the most important solar magnetic features creating the space weather is the solar wind that originates from the coronal holes (CHs).
The identification of the CHs on the Sun as one of the source regions of the solar wind is therefore crucial to achieve predictive capabilities.
In this study, we used an unsupervised machine-learning method, k-means, to pixel-wise cluster the passband images of the Sun taken by the Atmospheric Imaging Assembly on the Solar Dynamics Observatory in 171, 193, and 211 angstrom in different combinations.
Our results show that the pixel-wise k-means clustering together with systematic pre- and postprocessing steps provides compatible results with those from complex methods, such as convolutional neural networks.
More importantly, our study shows that there is a need for a CH database where a consensus about the CH boundaries is reached by observers independently.
This database then can be used as the "ground truth," when using a supervised method or just to evaluate the goodness of the models.
Deformation associated with plate convergence at subduction zones is accommodated by a complex system involving fault slip and viscoelastic flow. These processes have proven difficult to disentangle. The 2010 M-w 8.8 Maule earthquake occurred close to the Chilean coast within a dense network of continuously recording Global Positioning System stations, which provide a comprehensive history of surface strain. We use these data to assemble a detailed picture of a structurally controlled megathrust fault frictional patchwork and the three-dimensional rheological and time-dependent viscosity structure of the lower crust and upper mantle, all of which control the relative importance of afterslip and viscoelastic relaxation during postseismic deformation. These results enhance our understanding of subduction dynamics including the interplay of localized and distributed deformation during the subduction zone earthquake cycle.
Soil bacteria play a fundamental role in pedogenesis. However, knowledge about both the impact of climate and slope aspects on microbial communities and the consequences of these items in pedogenesis is lacking. Therefore, soil-bacterial communities from four sites and two different aspects along the climate gradient of the Chilean Coastal Cordillera were investigated. Using a combination of microbiological and physicochemical methods, soils that developed in arid, semi-arid, mediterranean, and humid climates were analyzed. Proteobacteria, Acidobacteria, Chloroflexi, Verrucomicrobia, and Planctomycetes were found to increase in abundance from arid to humid climates, while Actinobacteria and Gemmatimonadetes decreased along the transect. Bacterial-community structure varied with climate and aspect and was influenced by pH, bulk density, plant-available phosphorus, clay, and total organic-matter content. Higher bacterial specialization was found in arid and humid climates and on the south-facing slope and was likely promoted by stable microclimatic conditions. The presence of specialists was associated with ecosystem-functional traits, which shifted from pioneers that accumulated organic matter in arid climates to organic decomposers in humid climates. These findings provide new perspectives on how climate and slope aspects influence the composition and functional capabilities of bacteria, with most of these capabilities being involved in pedogenetic processes.
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.
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.
Natural gas hydrate occurrences contain predominantly methane; however, there are increasing reports of complex mixed gas hydrates and coexisting hydrate phases. Changes in the feed gas composition due to the preferred incorporation of certain components into the hydrate phase and an inadequate gas supply is often assumed to be the cause of coexisting hydrate phases. This could also be the case for the gas hydrate system in Qilian Mountain permafrost (QMP), which is mainly controlled by pores and fractures with complex gas compositions. This study is dedicated to the experimental investigations on the formation process of mixed gas hydrates based on the reservoir conditions in QMP. Hydrates were synthesized from water and a gas mixture under different gas supply conditions to study the effects on the hydrate formation process. In situ Raman spectroscopic measurements and microscopic observations were applied to record changes in both gas and hydrate phase over the whole formation process. The results demonstrated the effects of gas flow on the composition of the resulting hydrate phase, indicating a competitive enclathration of guest molecules into the hydrate lattice depending on their properties. Another observation was that despite significant changes in the gas composition, no coexisting hydrate phases were formed.
Understanding the key factors influencing the water quality of large river systems forms an important basis for the assessment and protection of cross-regional ecosystems and the implementation of adapted water management concepts. However, identifying these factors requires in-depth comprehension of the unique environmental systems, which can only be achieved by detailed water quality monitoring.
Within the scope of the joint science and sports event "Elbschwimmstaffel" (swimming relay on the river Elbe) in June/July 2017 organized by the German Ministry of Education and Research, water quality data were acquired along a 550 km long stretch of the Elbe River in Germany. During the survey, eight physiochemical water quality parameters were recorded in high spatial and temporal resolution with the BIOFISH multisensor system. Multivariate statistical methods were applied to identify and delineate processes influencing the water quality.
The BIOFISH dataset revealed that phytoplankton activity has a major impact on the water quality of the Elbe River in the summer months. The results suggest that phytoplankton biomass constitutes a substantial proportion of the suspended particles and that photosynthetic activity of phytoplankton is closely related to significant temporal changes in pH and oxygen saturation.
An evaluation of the BIOFISH data based on the combination of statistical analysis with weather and discharge data shows that the hydrological and meteorological history of the sampled water body was the main driver of phytoplankton dynamics. This study demonstrates the capacity of longitudinal river surveys with the BIOFISH or similar systems for water quality assessment, the identification of pollution sources and their utilization for online in situ monitoring of rivers.
Intercomparison of cosmic-ray neutron sensors and water balance monitoring in an urban environment
(2018)
Sensor-to-sensor variability is a source of error common to all geoscientific instruments that needs to be assessed before comparative and applied research can be performed with multiple sensors. Consistency among sensor systems is especially critical when subtle features of the surrounding terrain are to be identified. Cosmic-ray neutron sensors (CRNSs) are a recent technology used to monitor hectometre-scale environmental water storages, for which a rigorous comparison study of numerous co-located sensors has not yet been performed. In this work, nine stationary CRNS probes of type "CRS1000" were installed in relative proximity on a grass patch surrounded by trees, buildings, and sealed areas. While the dynamics of the neutron count rates were found to be similar, offsets of a few percent from the absolute average neutron count rates were found. Technical adjustments of the individual detection parameters brought all instruments into good agreement. Furthermore, we found a critical integration time of 6 h above which all sensors showed consistent dynamics in the data and their RMSE fell below 1% of gravimetric water content. The residual differences between the nine signals indicated local effects of the complex urban terrain on the scale of several metres. Mobile CRNS measurements and spatial simulations with the URANOS neutron transport code in the surrounding area (25 ha) have revealed substantial sub-footprint heterogeneity to which CRNS detectors are sensitive despite their large averaging volume. The sealed and constantly dry structures in the footprint furthermore damped the dynamics of the CRNS-derived soil moisture. We developed strategies to correct for the sealed-area effect based on theoretical insights about the spatial sensitivity of the sensor. This procedure not only led to reliable soil moisture estimation during dry-out periods, it further revealed a strong signal of intercepted water that emerged over the sealed surfaces during rain events. The presented arrangement offered a unique opportunity to demonstrate the CRNS performance in complex terrain, and the results indicated great potential for further applications in urban climate research.
The Northeast German Lowland Observatory (TERENO-NE) was established to investigate the regional impact of climate and land use change. TERENO-NE focuses on the Northeast German lowlands, for which a high vulnerability has been determined due to increasing temperatures and decreasing amounts of precipitation projected for the coming decades. To facilitate in-depth evaluations of the effects of climate and land use changes and to separate the effects of natural and anthropogenic drivers in the region, six sites were chosen for comprehensive monitoring. In addition, at selected sites, geoarchives were used to substantially extend the instrumental records back in time. It is this combination of diverse disciplines working across different time scales that makes the observatory TERENO-NE a unique observation platform. We provide information about the general characteristics of the observatory and its six monitoring sites and present examples of interdisciplinary research activities at some of these sites. We also illustrate how monitoring improves process understanding, how remote sensing techniques are fine-tuned by the most comprehensive ground-truthing site DEMMIN, how soil erosion dynamics have evolved, how greenhouse gas monitoring of rewetted peatlands can reveal unexpected mechanisms, and how proxy data provides a long-term perspective of current ongoing changes.
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.
Ice-rich permafrost has been subject to abrupt thaw and thermokarst formation in the past and is vulnerable to current global warming. The ice-rich permafrost domain includes Yedoma sediments that have never thawed since deposition during the late Pleistocene and Alas sediments that were formed by previous thermokarst processes during the Lateglacial and Holocene warming. Permafrost thaw unlocks organic carbon (OC) and minerals from these deposits and exposes OC to mineralization. A portion of the OC can be associated with iron (Fe), a redox-sensitive element acting as a trap for OC. Post-depositional thaw processes may have induced changes in redox conditions in these deposits and thereby affected Fe distribution and interactions between OC and Fe, with knock-on effects on the role that Fe plays in mediating present day OC mineralization. To test this hypothesis, we measured Fe concentrations and proportion of Fe oxides and Fe complexed with OC in unthawed Yedoma and previously thawed Alas deposits. Total Fe concentrations were determined on 1,292 sediment samples from the Yedoma domain using portable X-ray fluorescence; these concentrations were corrected for trueness using a calibration based on a subset of 144 samples measured by inductively coupled plasma optical emission spectrometry after alkaline fusion (R (2) = 0.95). The total Fe concentration is stable with depth in Yedoma deposits, but we observe a depletion or accumulation of total Fe in Alas deposits, which experienced previous thaw and/or flooding events. Selective Fe extractions targeting reactive forms of Fe on unthawed and previously thawed deposits highlight that about 25% of the total Fe is present as reactive species, either as crystalline or amorphous oxides, or complexed with OC, with no significant difference in proportions of reactive Fe between Yedoma and Alas deposits. These results suggest that redox driven processes during past thermokarst formation impact the present-day distribution of total Fe, and thereby the total amount of reactive Fe in Alas versus Yedoma deposits. This study highlights that ongoing thermokarst lake formation and drainage dynamics in the Arctic influences reactive Fe distribution and thereby interactions between Fe and OC, OC mineralization rates, and greenhouse gas emissions.
LegacyPollen 1.0
(2022)
Here we describe the LegacyPollen 1.0, a dataset of 2831 fossil pollen records with metadata, a harmonized taxonomy, and standardized chronologies.
A total of 1032 records originate from North America, 1075 from Europe, 488 from Asia, 150 from Latin America, 54 from Africa, and 32 from the Indo-Pacific.
The pollen data cover the late Quaternary (mostly the Holocene). The original 10 110 pollen taxa names (including variations in the notations) were harmonized to 1002 terrestrial taxa (including Cyperaceae), with woody taxa and major herbaceous taxa harmonized to genus level and other herbaceous taxa to family level.
The dataset is valuable for synthesis studies of, for example, taxa areal changes, vegetation dynamics, human impacts (e.g., deforestation), and climate change at global or continental scales.
The harmonized pollen and metadata as well as the harmonization table are available from PANGAEA (https://doi.org/10.1594/PANGAEA.929773; Herzschuh et al., 2021). R code for the harmonization is provided at Zenodo (https://doi.org/10.5281/zenodo.5910972; Herzschuh et al., 2022) so that datasets at a customized harmonization level can be easily established.
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.
Proglacial environments are ideal for studying the development of soils through the changes of rocks exposed by glacier retreat to weathering and microbial processes. Carbon (C) and nitrogen (N) contents as well as soil pH and soil elemental compositions are thought to be dominant factors structuring the bacterial, archaeal and fungal communities in the early stages of soil ecosystem formation. However, the functional linkages between C and N contents, soil composition and microbial community structures remain poorly understood. Here, we describe a multivariate analysis of geochemical properties and associated microbial community structures between a moraine and a glaciofluvial outwash in the proglacial area of a High Arctic glacier (Longyearbreen, Svalbard). Our results reveal distinct differences in developmental stages and heterogeneity between the moraine and the glaciofluvial outwash. We observed significant relationships between C and N contents, delta C-13(org) and delta N-15 isotopic ratios, weathering and microbial abundance and community structures. We suggest that the observed differences in microbial and geochemical parameters between the moraine and the glaciofluvial outwash are primarily a result of geomorphological variations of the proglacial terrain.
Geomorphic strain markers accumulating the effects of many earthquake cycles help to constrain the mechanical behaviour of continental rift systems as well as the related seismic hazards. In the Corinth Rift (Greece), the unique record of onshore and offshore markers of Pleistocene similar to 100-ka climate cycles provides an outstanding possibility to constrain rift mechanics over a range of timescales. Here we use high-resolution topography to analyse the 3D geometry of a sequence of Pleistocene emerged marine terraces associated with flexural rift-flank uplift. We integrate this onshore dataset with offshore seismic data to provide a synoptic view of the flexural deformation across the rift. This allows us to derive an average slip rate of 4.5-9.0 mm.yr(-1) on the master fault over the past similar to 610 ka and an uplift/ subsidence ratio of 1:1.1-2.4. We reproduce the observed flexure patterns, using 3 and 5-layered lithospheric scale finite element models. Modelling results imply that the observed elastic flexure is produced by coseismic slip along 40-60 degrees planar normal faults in the elastic upper crust, followed by postseismic viscous relaxation occurring within the basal lower crust or upper mantle. We suggest that such a mechanism may typify rapid localised extension of continental lithosphere.
Dryland xeric conditions exert a deterministic effect on microbial communities, forcing life into refuge niches. Deposited rocks can form a lithic niche for microorganisms in desert regions. Mineral weathering is a key process in soil formation and the importance of microbial-driven mineral weathering for nutrient extraction is increasingly accepted. Advances in geobiology provide insight into the interactions between microorganisms and minerals that play an important role in weathering processes. In this study, we present the examination of the microbial diversity in dryland rocks from the Tsauchab River banks in Namibia. We paired culture-independent 16S rRNA gene amplicon sequencing with culture-dependent (isolation of bacteria) techniques to assess the community structure and diversity patterns. Bacteria isolated from dryland rocks are typical of xeric environments and are described as being involved in rock weathering processes. For the first time, we extracted extra- and intracellular DNA from rocks to enhance our understanding of potentially rock-weathering microorganisms. We compared the microbial community structure in different rock types (limestone, quartz-rich sandstone and quartz-rich shale) with adjacent soils below the rocks. Our results indicate differences in the living lithic and sublithic microbial communities.
During eruptive activity of andesitic stratovolcanoes, the extrusion of lava domes, their collapse and intermittent explosions are common volcanic hazards. Many lava domes grow in a preferred direction, in turn affecting the direction of lava flows and pyroclastic density currents. Access to active lava domes is difficult and hazardous, so detailed data characterizing lava dome growth are typically limited, keeping the processes controlling the directionality of extrusions unclear. Here we combine TerraSAR-X satellite radar observations with high-resolution airborne photogrammetry to assess morphological changes, and perform finite element modeling to investigate the impact of loading stress on shallow magma ascent directions associated with lava dome extrusion and crater formation at Volcan de Colima, Mexico. The TerraSAR-X data, acquired in similar to 1-m resolution spotlight mode, enable us to derive a chronology of the eruptive processes from intensity-based time-lapse observations of the general crater and dome evolution. The satellite images are complemented by close-range airborne photos, processed by the Structure-from-Motion workflow. This allows the derivation of high-resolution digital elevation models, providing insight into detailed loading and unloading features. During the observation period from Jan-2013 to Feb-2016, we identify a dominantly W-directed dome growth and lava flow production until Jan-2015. In Feb-2015, following the removal of the active summit dome, the surface crater widened and elongated along a NE-SW axis. Later in May-2015, a new dome grew toward the SW of the crater while a separate vent developed in the NE of the crater, reflecting a change in the direction of magma ascent and possible conduit bifurcation. Finite element models show a significant stress change in agreement with the observed magma ascent direction changes in response to the changing surface loads, both for loading (dome growth) and unloading (crater forming excavation) cases. These results allow insight into shallow dome growth dynamics and the migration of magma ascent in response to changing volcano summit morphology. They further highlight the importance of detailed volcano summit morphology surveillance, as changes in direction or location of dome extrusion may have major implications regarding the directions of potential volcanic hazards, such as pyroclastic density currents generated by dome collapse.
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.
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.
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.
We examined the relationship between the mechanical strength of the lithosphere and the distribution of seismicity within the overriding continental plate of the southern Central Andes (SCA, 29 degrees-39 degrees S), where the oceanic Nazca Plate changes its subduction angle between 33 degrees S and 35 degrees S, from subhorizontal in the north (<5 degrees) to steep in the south (similar to 30 degrees). We computed the long-term lithospheric strength based on an existing 3D model describing variations in thickness, density, and temperature of the main geological units forming the lithosphere of the SCA and adjacent forearc and foreland regions. The comparison between our results and seismicity within the overriding plate (upper-plate seismicity) shows that most of the events occur within the modeled brittle domain of the lithosphere. The depth where the deformation mode switches from brittle frictional to thermally activated ductile creep provides a conservative lower bound to the seismogenic zone in the overriding plate of the study area. We also found that the majority of upper-plate earthquakes occurs within the realm of first-order contrasts in integrated strength (12.7-13.3 log Pam in the Andean orogen vs. 13.5-13.9 log Pam in the forearc and the foreland). Specific conditions characterize the mechanically strong northern foreland of the Andes, where seismicity is likely explained by the effects of slab steepening.
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.
The current awareness of the high importance of urban green leads to a stronger need for tools to comprehensively represent urban green and its benefits. A common scientific approach is the development of urban ecosystem services (UES) based on remote sensing methods at the city or district level. Urban planning, however, requires fine-grained data that match local management practices. Hence, this study linked local biotope and tree mapping methods to the concept of ecosystem services. The methodology was tested in an inner-city district in SW Germany, comparing publicly accessible areas and non-accessible courtyards. The results provide area-specific [m(2)] information on the green inventory at the microscale, whereas derived stock and UES indicators form the basis for comparative analyses regarding climate adaptation and biodiversity. In the case study, there are ten times more micro-scale green spaces in private courtyards than in the public space, as well as twice as many trees. The approach transfers a scientific concept into municipal planning practice, enables the quantitative assessment of urban green at the microscale and illustrates the importance for green stock data in private areas to enhance decision support in urban development. Different aspects concerning data collection and data availability are critically discussed.
Tectonically active coasts are dynamic environments characterized by the presence of multiple marine terraces formed by the combined effects of wave erosion, tectonic uplift, and sea-level oscillations at glacialcycle timescales. Well-preserved erosional terraces from the last interglacial sea-level highstand are ideal marker horizons for reconstructing past sea-level positions and calculating vertical displacement rates. We carried out an almost continuous mapping of the last interglacial marine terrace along similar to 5000 km of the western coast of South America between 1 degrees N and 40 degrees S. We used quantitatively replicable approaches constrained by published terrace-age estimates to ultimately compare elevations and patterns of uplifted terraces with tectonic and climatic parameters in order to evaluate the controlling mechanisms for the formation and preservation of marine terraces and crustal deformation. Uncertainties were estimated on the basis of measurement errors and the distance from referencing points. Overall, our results indicate a median elevation of 30.1 m, which would imply a median uplift rate of 0.22 m kyr(-1) averaged over the past similar to 125 kyr. The patterns of terrace elevation and uplift rate display high-amplitude (similar to 100-200 m) and long-wavelength (similar to 10(2) km) structures at the Manta Peninsula (Ecuador), the San Juan de Marcona area (central Peru), and the Arauco Peninsula (south-central Chile). Medium-wavelength structures occur at the Mejillones Peninsula and Topocalma in Chile, while short-wavelength (< 10 km) features are for instance located near Los Vilos, Valparaiso, and Carranza, Chile. We interpret the long-wavelength deformation to be controlled by deep-seated processes at the plate interface such as the subduction of major bathymetric anomalies like the Nazca and Carnegie ridges. In contrast, short-wavelength deformation may be primarily controlled by sources in the upper plate such as crustal faulting, which, however, may also be associated with the subduction of topographically less pronounced bathymetric anomalies. Latitudinal differences in climate additionally control the formation and preservation of marine terraces. Based on our synopsis we propose that increasing wave height and tidal range result in enhanced erosion and morphologically well-defined marine terraces in south-central Chile. Our study emphasizes the importance of using systematic measurements and uniform, quantitative methodologies to characterize and correctly interpret marine terraces at regional scales, especially if they are used to unravel the tectonic and climatic forcing mechanisms of their formation. This database is an integral part of the World Atlas of Last Interglacial Shorelines (WALIS), published online at https://doi.org/10.5281/zenodo.4309748 (Freisleben et al., 2020).
An earthquake swarm affected the Bransfield Strait, Antarctica, a unique rift basin in transition from intra-arc rifting to ocean spreading. The swarm, counting similar to 85,000 volcano-tectonic earthquakes since August 2020, is located close to the Orca submarine volcano, previously considered inactive. Simultaneously, geodetic data reported up to similar to 11 cm north-westward displacement over King George Island. We use a broad variety of geophysical data and methods to reveal the complex migration of seismicity, accompanying the intrusion of 0.26-0.56 km(3) of magma. Strike-slip earthquakes mark the intrusion at depth, while shallower normal faulting the similar to 20 km long lateral growth of a dike. Seismicity abruptly decreased after a Mw 6.0 earthquake, suggesting the magmatic dike lost pressure with the slipping of a large fault. A seafloor eruption is likely, but not confirmed by sea surface temperature anomalies. The unrest documents episodic magmatic intrusion in the Bransfield Strait, providing unique insights into active continental rifting.
The Ice, Cloud, and Land Elevation Satellite-2 (ICESat-2) with its land and vegetation height data product (ATL08), and Global Ecosystem Dynamics Investigation (GEDI) with its terrain elevation and height metrics data product (GEDI Level 2A) missions have great potential to globally map ground and canopy heights. Canopy height is a key factor in estimating above-ground biomass and its seasonal changes; these satellite missions can also improve estimated above-ground carbon stocks. This study presents a novel Sparse Vegetation Detection Algorithm (SVDA) which uses ICESat-2 (ATL03, geolocated photons) data to map tree and vegetation heights in a sparsely vegetated savanna ecosystem. The SVDA consists of three main steps: First, noise photons are filtered using the signal confidence flag from ATL03 data and local point statistics. Second, we classify ground photons based on photon height percentiles. Third, tree and grass photons are classified based on the number of neighbors. We validated tree heights with field measurements (n = 55), finding a root-mean-square error (RMSE) of 1.82 m using SVDA, GEDI Level 2A (Geolocated Elevation and Height Metrics product): 1.33 m, and ATL08: 5.59 m. Our results indicate that the SVDA is effective in identifying canopy photons in savanna ecosystems, where ATL08 performs poorly. We further identify seasonal vegetation height changes with an emphasis on vegetation below 3 m; widespread height changes in this class from two wet-dry cycles show maximum seasonal changes of 1 m, possibly related to seasonal grass-height differences. Our study shows the difficulties of vegetation measurements in savanna ecosystems but provides the first estimates of seasonal biomass changes.
The intensification of Northern Hemisphere glaciations at the end of the Pliocene epoch marks one of the most substantial climatic shifts of the Cenozoic. Despite global cooling, sea surface temperatures in the high latitude North Atlantic Ocean rose between 2.9–2.7 million years ago. Here we present sedimentary geochemical proxy data from the Gulf of Cadiz to reconstruct the variability of Mediterranean Outflow Water, an important heat source to the North Atlantic. We find evidence for enhanced production of Mediterranean Outflow from the mid-Pliocene to the late Pliocene which we infer could have driven a sub-surface heat channel into the high-latitude North Atlantic. We then use Earth System Models to constrain the impact of enhanced Mediterranean Outflow production on the northward heat transport in the North Atlantic. In accord with the proxy data, the numerical model results support the formation of a sub-surface channel that pumped heat from the subtropics into the high latitude North Atlantic. We further suggest that this mechanism could have delayed ice sheet growth at the end of the Pliocene.
Meteorological conditions during the ACLOUD/PASCAL field campaign near Svalbard in early summer 2017
(2018)
The two concerted field campaigns, Arctic CLoud Observations Using airborne measurements during polar Day (ACLOUD) and the Physical feedbacks of Arctic planetary boundary level Sea ice, Cloud and AerosoL (PASCAL), took place near Svalbard from 23 May to 26 June 2017. They were focused on studying Arctic mixed-phase clouds and involved observations from two airplanes (ACLOUD), an icebreaker (PASCAL) and a tethered balloon, as well as ground-based stations. Here, we present the synoptic development during the 35-day period of the campaigns, using near-surface and upper-air meteorological observations, as well as operational satellite, analysis, and reanalysis data. Over the campaign period, short-term synoptic variability was substantial, dominating over the seasonal cycle. During the first campaign week, cold and dry Arctic air from the north persisted, with a distinct but seasonally unusual cold air outbreak. Cloudy conditions with mostly low-level clouds prevailed. The subsequent 2 weeks were characterized by warm and moist maritime air from the south and east, which included two events of warm air advection. These synoptical disturbances caused lower cloud cover fractions and higher-reaching cloud systems. In the final 2 weeks, adiabatically warmed air from the west dominated, with cloud properties strongly varying within the range of the two other periods. Results presented here provide synoptic information needed to analyze and interpret data of upcoming studies from ACLOUD/PASCAL, while also offering unprecedented measurements in a sparsely observed region.
Microbial community composition and abundance after millennia of submarine permafrost warming
(2019)
Warming of the Arctic led to an increase in permafrost temperatures by about 0.3 degrees C during the last decade. Permafrost warming is associated with increasing sediment water content, permeability, and diffusivity and could in the long term alter microbial community composition and abundance even before permafrost thaws. We studied the long-term effect (up to 2500 years) of submarine permafrost warming on microbial communities along an onshore-offshore transect on the Siberian Arctic Shelf displaying a natural temperature gradient of more than 10 degrees C. We analysed the in situ development of bacterial abundance and community composition through total cell counts (TCCs), quantitative PCR of bacterial gene abundance, and amplicon sequencing and correlated the microbial community data with temperature, pore water chemistry, and sediment physicochemical parameters. On timescales of centuries, permafrost warming coincided with an overall decreasing microbial abundance, whereas millennia after warming microbial abundance was similar to cold onshore permafrost. In addition, the dissolved organic carbon content of all cores was lowest in submarine permafrost after millennial-scale warming. Based on correlation analysis, TCC, unlike bacterial gene abundance, showed a significant rank-based negative correlation with increasing temperature, while bacterial gene copy numbers showed a strong negative correlation with salinity. Bacterial community composition correlated only weakly with temperature but strongly with the pore water stable isotopes delta O-18 and delta D, as well as with depth. The bacterial community showed substantial spatial variation and an overall dominance of Actinobacteria, Chloroflexi, Firmicutes, Gemmatimonadetes, and Proteobacteria, which are amongst the microbial taxa that were also found to be active in other frozen permafrost environments. We suggest that, millennia after permafrost warming by over 10 degrees C, microbial community composition and abundance show some indications for proliferation but mainly reflect the sedimentation history and paleoenvironment and not a direct effect through warming.
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.
Magmatic continental rifts often constitute nascent plate boundaries, yet long-term extension rates and transient rate changes associated with these early stages of continental breakup remain difficult to determine. Here, we derive a time-averaged minimum extension rate for the inner graben of the Northern Kenya Rift (NKR) of the East African Rift System for the last 0.5 m.y. We use the TanDEM-X science digital elevation model to evaluate fault-scarp geometries and determine fault throws across the volcano-tectonic axis of the inner graben of the NKR. Along rift-perpendicular profiles, amounts of cumulative extension are determined, and by integrating four new Ar-40/Ar-39 radiometric dates for the Silali volcano into the existing geochronology of the faulted volcanic units, time-averaged extension rates are calculated. This study reveals that in the inner graben of the NKR, the long-term extension rate based on mid-Pleistocene to recent brittle deformation has minimum values of 1.0-1.6 mm yr(-1), locally with values up to 2.0 mm yr(-1). A comparison with the decadal, geodetically determined extension rate reveals that at least 65% of the extension must be accommodated within a narrow, 20-km-wide zone of the inner rift. In light of virtually inactive border faults of the NKR, we show that extension is focused in the region of the active volcano-tectonic axis in the inner graben, thus highlighting the maturing of continental rifting in the NKR.
The Issyk-Kul Basin (Kyrgyzstan), situated in the central Tian Shan Mountains, hosts the largest and deepest mountain lake in Central Asia. Erosion of the surrounding Terskey and Kungey ranges led to the accumulation of up to 4 km of sediment in the adjacent depression. Creation of the basin from regional shortening and uplift likely initiated around the Oligocene-Miocene, yet precise age control is sparse. To better understand the timing of these processes, we obtained magnetostratigraphic age constraints on fossil-poor, fluvio-lacustrine sediments exposed south of Lake Issyk-Kul, that agree well with previous age constraints of the equivalent strata outside the Issyk-Kul Basin. Two 500-650 m thick sections comprised mainly of Chu Group sediments were dated at 6.3-2.8 Ma and 7.0-2.4 Ma (late Miocene to early Pleistocene). Together with reinterpreted magnetostratigraphic constraints from underlying strata, we find that syn-tectonic deposition commenced at similar to 22 Ma with average sedimentation rates <10 cm/ka. Sedimentation rates increased to 10-30 cm/ka at 7 Ma, concurrent with accelerated uplift in the Terskey Range to the south. A deformation event in one section (Kaji-Say) between 5 and 3 Ma together with concurrent shifts of depositional centers throughout the basin signal the onset of substantial uplift of the Kungey Range to the north at similar to 5 Ma. This uplift and deformation transformed the Issyk-Kul area into a closed basin that facilitated the formation of a deep lake. Lacustrine facies deposited around 3 Ma mark the existence of Lake Issyk-Kul by that time.
Eight d-metal-containing N-butylpyridinium ionic liquids (ILs) with the nominal composition (C4Py)2[Ni0.5M0.5Cl4] or (C4Py)2[Zn0.5M0.5Cl4] (M = Cu, Co, Mn, Ni, Zn; C4Py = N-butylpyridinium) were synthesized, characterized, and investigated for their optical properties. Single crystal and powder X-ray analysis shows that the compounds are isostructural to existing examples based on other d-metal ions. Inductively coupled plasma optical emission spectroscopy measurements confirm that the metal/metal ratio is around 50 : 50. UV-Vis spectroscopy shows that the optical absorption can be tuned by selection of the constituent metals. Moreover, the compounds can act as an optical sensor for the detection of gases such as ammonia as demonstrated via a simple prototype setup.
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.
Environmental factors shape the spatial distribution and dynamics of populations. Understanding how these factors interact with movement behavior is critical for efficient conservation, in particular for migratory species. Adult female green sea turtles, Chelonia mydas, migrate between foraging and nesting sites that are generally separated by thousands of kilometers. As an emblematic endangered species, green turtles have been intensively studied, with a focus on nesting, migration, and foraging. Nevertheless, few attempts integrated these behaviors and their trade‐offs by considering the spatial configurations of foraging and nesting grounds as well as environmental heterogeneity like oceanic currents and food distribution. We developed an individual‐based model to investigate the impact of local environmental conditions on emerging migratory corridors and reproductive output and to thereby identify conservation priority sites. The model integrates movement, nesting, and foraging behavior. Despite being largely conceptual, the model captured realistic movement patterns which confirm field studies. The spatial distribution of migratory corridors and foraging hot spots was mostly constrained by features of the regional landscape, such as nesting site locations, distribution of feeding patches, and oceanic currents. These constraints also explained the mixing patterns in regional forager communities. By implementing alternative decision strategies of the turtles, we found that foraging site fidelity and nesting investment, two characteristics of green turtles' biology, are favorable strategies under unpredictable environmental conditions affecting their habitats. Based on our results, we propose specific guidelines for the regional conservation of green turtles as well as future research suggestions advancing spatial ecology of sea turtles. Being implemented in an easy to learn open‐source software, our model can coevolve with the collection and analysis of new data on energy budget and movement into a generic tool for sea turtle research and conservation. Our modeling approach could also be useful for supporting the conservation of other migratory marine animals.
Germany and the United Kingdom have domestic shale gas reserves which they may exploit in the future to complement their national energy strategies. However gas production releases volatile organic compounds (VOC) and nitrogen oxides (NOx), which through photochemical reaction form ground-level ozone, an air pollutant that can trigger adverse health effects e.g. on the respiratory system. This study explores the range of impacts of a potential shale gas industry in these two countries on local and regional ambient ozone. To this end, comprehensive emission scenarios are used as the basis for input to an online-coupled regional chemistry transport model (WRF-Chem). Here we simulate shale gas scenarios over summer (June, July, August) 2011, exploring the effects of varying VOC emissions, gas speciation, and concentration of NOx emissions over space and time, on ozone formation. An evaluation of the model setup is performed, which exhibited the model’s ability to predict surface meteorological and chemical variables well compared with observations, and consistent with other studies. When different shale gas scenarios were employed, the results show a peak increase in maximum daily 8-hour average ozone from 3.7 to 28.3 μg m–3. In addition, we find that shale gas emissions can force ozone exceedances at a considerable percentage of regulatory measurement stations locally (up to 21% in Germany and 35% in the United Kingdom) and in distant countries through long-range transport, and increase the cumulative health-related metric SOMO35 (maximum percent increase of ~28%) throughout the region. Findings indicate that VOC emissions are important for ozone enhancement, and to a lesser extent NOx, meaning that VOC regulation for a future European shale gas industry will be of especial importance to mitigate unfavorable health outcomes. Overall our findings demonstrate that shale gas production in Europe can worsen ozone air quality on both the local and regional scales.
Modern pollen-vegetation-climate relationships underpin palaeovegetation and palaeoclimate reconstructions from fossil pollen records. East Siberia is an ideal area for investigating the relationships between modern pollen assemblages and near natural vegetation under cold continental climate conditions. Reliable pollen-based quantitative vegetation and climate reconstructions are still scarce due to the limited number of modern pollen datasets. Furthermore, differences in pollen representation of samples from lake sediments and soils are not well understood. Here, we present a new pollen dataset of 48 moss/soil and 24 lake surface-sediment samples collected in Chukotka and central Yakutia in East Siberia. The pollen-vegetation-climate relationships were investigated by ordination analyses. Generally, tundra and taiga vegetation types can be well distinguished in the surface pollen assemblages. Moss/soil and lake samples contain generally similar pollen assemblages as revealed by a Procrustes comparison with some exceptions. Overall, modern pollen assemblages reflect the temperature and precipitation gradients in the study areas as revealed by constrained ordination analysis. We estimate the relative pollen productivity (RPP) of major taxa and the relevant source area of pollen (RSAP) for moss/soil samples from Chukotka and central Yakutia using Extended R-Value (ERV) analysis. The RSAP of the tundra-forest transition area in Chukotka and taiga area in central Yakutia are ca. 1300 and 360 m, respectively. For Chukotka, RPPs relative to both Poaceae and Ericaceae were estimated while RPPs for central Yakutia were relative only to Ericaceae. Relative to Ericaceae (reference taxon, RPP = 1), Larix, Betula, Picea, and Pinus are overrepresented while Alnus, Cyperaceae, Poaceae, and Salix are underrepresented in the pollen spectra. Our estimates are in general agreement with previously published values and provide the basis for reliable quantitative reconstructions of East Siberian vegetation.
Monitoring the response of volcanic CO2 emissions to changes in the Los Humeros hydrothermal system
(2021)
Carbon dioxide is the most abundant, non-condensable gas in volcanic systems, released into the atmosphere through either diffuse or advective fluid flow. The emission of substantial amounts of CO2 at Earth's surface is not only controlled by volcanic plumes during periods of eruptive activity or fumaroles, but also by soil degassing along permeable structures in the subsurface. Monitoring of these processes is of utmost importance for volcanic hazard analyses, and is also relevant for managing geothermal resources. Fluid-bearing faults are key elements of economic value for geothermal power generation. Here, we describe for the first time how sensitively and quickly natural gas emissions react to changes within a deep hydrothermal system due to geothermal fluid reinjection. For this purpose, we deployed an automated, multi-chamber CO2 flux monitoring system within the damage zone of a deep-rooted major normal fault in the Los Humeros Volcanic Complex (LHVC) in Mexico and recorded data over a period of five months. After removing the atmospheric effects on variations in CO2 flux, we calculated correlation coefficients between residual CO2 emissions and reinjection rates, identifying an inverse correlation of rho = - 0.51 to - 0.66. Our results indicate that gas emissions respond to changes in reinjection rates within 24 h, proving an active hydraulic communication between the hydrothermal system and Earth's surface. This finding is a promising indication not only for geothermal reservoir monitoring but also for advanced long-term volcanic risk analysis. Response times allow for estimation of fluid migration velocities, which is a key constraint for conceptual and numerical modelling of fluid flow in fracture-dominated systems.
Water bodies are a highly abundant feature of Arctic permafrost ecosystems and strongly influence their hydrology, ecology and biogeochemical cycling. While very high resolution satellite images enable detailed mapping of these water bodies, the increasing availability and abundance of this imagery calls for fast, reliable and automatized monitoring. This technical work presents a largely automated and scalable workflow that removes image noise, detects water bodies, removes potential misclassifications from infrastructural features, derives lake shoreline geometries and retrieves their movement rate and direction on the basis of ortho-ready very high resolution satellite imagery from Arctic permafrost lowlands. We applied this workflow to typical Arctic lake areas on the Alaska North Slope and achieved a successful and fast detection of water bodies. We derived representative values for shoreline movement rates ranging from 0.40-0.56 m yr(-1) for lake sizes of 0.10 ha-23.04 ha. The approach also gives an insight into seasonal water level changes. Based on an extensive quantification of error sources, we discuss how the results of the automated workflow can be further enhanced by incorporating additional information on weather conditions and image metadata and by improving the input database. The workflow is suitable for the seasonal to annual monitoring of lake changes on a sub-meter scale in the study areas in northern Alaska and can readily be scaled for application across larger regions within certain accuracy limitations.
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
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.
A new solid-state material, N-butyl pyridinium diiodido argentate(I), is synthesized using a simple and effective one-pot approach. In the solid state, the compound exhibits 1D ([AgI2](-))(n) chains that are stabilized by the N-butyl pyridinium cation. The 1D structure is further manifested by the formation of long, needle-like crystals, as revealed from electron microscopy. As the general composition is derived from metal halide-based ionic liquids, the compound has a low melting point of 100-101 degrees C, as confirmed by differential scanning calorimetry. Most importantly, the compound has a conductivity of 10(-6) S cm(-1) at room temperature. At higher temperatures the conductivity increases and reaches to 10(-4 )S cm(-1) at 70 degrees C. In contrast to AgI, however, the current material has a highly anisotropic 1D arrangement of the ionic domains. This provides direct and tuneable access to fast and anisotropic ionic conduction. The material is thus a significant step forward beyond current ion conductors and a highly promising prototype for the rational design of highly conductive ionic solid-state conductors for battery or solar cell applications.
The global geological volatile cycle (H, C, N) plays an important role in the long term self-regulation of the Earth system. However, the complex interaction between its deep, solid Earth components (i.e. crust and mantle), Earth's fluid envelopes (i.e. atmosphere and hydrosphere) and plate tectonic processes is a subject of ongoing debate. In this study we want to draw attention to how the presence of primary melt (MI) and fluid (FI) inclusions in high-grade metamorphic minerals could help constrain the crustal component of the volatile cycle. To that end, we review the distribution of MI and FI throughout Earth's history, from ca. 3.0 Ga ago up to the present day. We argue that the lower crust might constitute an important, long-term, volatile storage unit, capable to influence the composition of the surface envelopes through the mean of weathering, crustal thickening, partial melting and crustal assimilation during volcanic activity. Combined with thermodynamic modelling, our compilation indicates that periods of well-established plate tectonic regimes at <0.85 Ga and 1.7-2.1 Ga, might be more prone to the reworking of supracrustal lithologies and the storage of volatiles in the lower crust. Such hypothesis has implication beyond the scope of metamorphic petrology as it potentially links geodynamic mechanisms to habitable surface conditions. MI and FI in metamorphic crustal rocks then represent an invaluable archive to assess and quantify the co-joint evolution of plate tectonics and Earth's external processes. (C) 2021 China University of Geosciences (Beijing) and Peking University. Production and hosting by Elsevier B.V. This is an open access article under the CC BY-NC-ND license (http://creativecommons.org/licenses/by-nc-nd/4.0/).
We present two new empirical models of radiation belt electron flux at geostationary orbit. GOES-15 measurements of 0.8 MeV electrons were used to train a Nonlinear Autoregressive with Exogenous input (NARX) neural network for both modeling GOES-15 flux values and an upper boundary condition scaling factor (BF). The GOES-15 flux model utilizes an input and feedback delay of 2 and 2 time steps (i.e., 5 min time steps) with the most efficient number of hidden layers set to 10. Magnetic local time, Dst, Kp, solar wind dynamic pressure, AE, and solar wind velocity were found to perform as predicative indicators of GOES-15 flux and therefore were used as the exogenous inputs. The NARX-derived upper boundary condition scaling factor was used in conjunction with the Versatile Electron Radiation Belt (VERB) code to produce reconstructions of the radiation belts during the period of July-November 1990, independent of in-situ observations. Here, Kp was chosen as the sole exogenous input to be more compatible with the VERB code. This Combined Release and Radiation Effects Satellite-era reconstruction showcases the potential to use these neural network-derived boundary conditions as a method of hindcasting the historical radiation belts. This study serves as a companion paper to another recently published study on reconstructing the radiation belts during Solar Cycles 17-24 (Saikin et al., 2021, ), for which the results featured in this paper were used.
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.
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 Mackenzie Delta (MD) is a permafrost-bearing region along the coasts of the Canadian Arctic which exhibits high sub-permafrost gas hydrate (GH) reserves. The GH occurring at the Mallik site in the MD is dominated by thermogenic methane (CH4), which migrated from deep conventional hydrocarbon reservoirs, very likely through the present fault systems. Therefore, it is assumed that fluid flow transports dissolved CH4 upward and out of the deeper overpressurized reservoirs via the existing polygonal fault system and then forms the GH accumulations in the Kugmallit-Mackenzie Bay Sequences. We investigate the feasibility of this mechanism with a thermo-hydraulic-chemical numerical model, representing a cross section of the Mallik site. We present the first simulations that consider permafrost formation and thawing, as well as the formation of GH accumulations sourced from the upward migrating CH4-rich formation fluid. The simulation results show that temperature distribution, as well as the thickness and base of the ice-bearing permafrost are consistent with corresponding field observations. The primary driver for the spatial GH distribution is the permeability of the host sediments. Thus, the hypothesis on GH formation by dissolved CH4 originating from deeper geological reservoirs is successfully validated. Furthermore, our results demonstrate that the permafrost has been substantially heated to 0.8-1.3 degrees C, triggered by the global temperature increase of about 0.44 degrees C and further enhanced by the Arctic Amplification effect at the Mallik site from the early 1970s to the mid-2000s.
The LArge-scale Reservoir Simulator (LARS) has been previously developed to study hydrate dissociation in hydrate-bearing systems under in-situ conditions. In the present study, a numerical framework of equations of state describing hydrate formation at equilibrium conditions has been elaborated and integrated with a numerical flow and transport simulator to investigate a multi-stage hydrate formation experiment undertaken in LARS. A verification of the implemented modeling framework has been carried out by benchmarking it against another established numerical code. Three-dimensional (3D) model calibration has been performed based on laboratory data available from temperature sensors, fluid sampling, and electrical resistivity tomography. The simulation results demonstrate that temperature profiles, spatial hydrate distribution, and bulk hydrate saturation are consistent with the observations. Furthermore, our numerical framework can be applied to calibrate geophysical measurements, optimize post-processing workflows for monitoring data, improve the design of hydrate formation experiments, and investigate the temporal evolution of sub-permafrost methane hydrate reservoirs.
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.
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.
Records from ocean bottom seismometers (OBSs) are highly contaminated by noise, which is much stronger compared to data from most land stations, especially on the horizontal components. As a consequence, the high energy of the oceanic noise at frequencies below 1 Hz considerably complicates the analysis of the teleseismic earthquake signals recorded by OBSs.
Previous studies suggested different approaches to remove low-frequency noises from OBS recordings but mainly focused on the vertical component. The records of horizontal components, which are crucial for the application of many methods in passive seismological analysis of body and surface waves, could not be much improved in the teleseismic frequency band. Here we introduce a noise reduction method, which is derived from the harmonic–percussive separation algorithms used in Zali et al. (2021), in order to separate long-lasting narrowband signals from broadband transients in the OBS signal. This leads to significant noise reduction of OBS records on both the vertical and horizontal components and increases the earthquake signal-to-noise ratio (SNR) without distortion of the broadband earthquake waveforms. This is demonstrated through tests with synthetic data. Both SNR and cross-correlation coefficients showed significant improvements for different realistic noise realizations. The application of denoised signals in surface wave analysis and receiver functions is discussed through tests with synthetic and real data.
OpenForecast
(2019)
The development and deployment of new operational runoff forecasting systems are a strong focus of the scientific community due to the crucial importance of reliable and timely runoff predictions for early warnings of floods and flashfloods for local businesses and communities. OpenForecast, the first operational runoff forecasting system in Russia, open for public use, is presented in this study. We developed OpenForecast based only on open-source software and data-GR4J hydrological model, ERA-Interim meteorological reanalysis, and ICON deterministic short-range meteorological forecasts. Daily forecasts were generated for two basins in the European part of Russia. Simulation results showed a limited efficiency in reproducing the spring flood of 2019. Although the simulations managed to capture the timing of flood peaks, they failed in estimating flood volume. However, further implementation of the parsimonious data assimilation technique significantly alleviates simulation errors. The revealed limitations of the proposed operational runoff forecasting system provided a foundation to outline its further development and improvement.
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.
Organic carbon (OC) stored in Arctic permafrost represents one of Earth's largest and most vulnerable terrestrial carbon pools. Amplified climate warming across the Arctic results in widespread permafrost thaw. Permafrost deposits exposed at river cliffs and coasts are particularly susceptible to thawing processes. Accelerating erosion of terrestrial permafrost along shorelines leads to increased transfer of organic matter (OM) to nearshore waters. However, the amount of terrestrial permafrost carbon and nitrogen as well as the OM quality in these deposits is still poorly quantified. We define the OM quality as the intrinsic potential for further transformation, decomposition and mineralisation. Here, we characterise the sources and the quality of OM supplied to the Lena River at a rapidly eroding permafrost river shoreline cliff in the eastern part of the delta (Sobo-Sise Island). Our multi-proxy approach captures bulk elemental, molecu- lar geochemical and carbon isotopic analyses of Late Pleistocene Yedoma permafrost and Holocene cover deposits, discontinuously spanning the last similar to 52 kyr. We showed that the ancient permafrost exposed in the Sobo-Sise cliff has a high organic carbon content (mean of about 5 wt %). The oldest sediments stem from Marine Isotope Stage (MIS) 3 interstadial deposits (dated to 52 to 28 cal ka BP) and are overlaid by last glacial MIS 2 (dated to 28 to 15 cal ka BP) and Holocene MIS 1 (dated to 7-0 cal ka BP) deposits. The relatively high average chain length (ACL) index of n-alkanes along the cliff profile indicates a predominant contribution of vascular plants to the OM composition. The elevated ratio of isoand anteiso-branched fatty acids (FAs) relative to mid- and long-chain (C >= 20) n-FAs in the interstadial MIS 3 and the interglacial MIS 1 deposits suggests stronger microbial activity and consequently higher input of bacterial biomass during these climatically warmer periods. The overall high carbon preference index (CPI) and higher plant fatty acid (HPFA) values as well as high C/N ratios point to a good quality of the preserved OM and thus to a high potential of the OM for decomposition upon thaw. A decrease in HPFA values downwards along the profile probably indicates stronger OM decomposition in the oldest (MIS 3) deposits of the cliff. The characterisation of OM from eroding permafrost leads to a better assessment of the greenhouse gas potential of the OC released into river and nearshore waters in the future.
The Arctic is rich in aquatic systems and experiences rapid warming due to climate change. The accelerated warming causes permafrost thaw and the mobilization of organic carbon.
When dissolved organic carbon is mobilized, this DOC can be transported to aquatic systems and degraded in the water bodies and further downstream. Here, we analyze the influence of different landscape components on DOC concentrations and export in a small (6.45 km(2)) stream catchment in the Lena River Delta.
The catchment includes lakes and ponds, with the flow path from Pleistocene yedoma deposits across Holocene non-yedoma deposits to the river outlet. In addition to DOC concentrations, we use radiocarbon dating of DOC as well as stable oxygen and hydrogen isotopes (delta O-18 and delta D) to assess the origin of DOC.
We find significantly higher DOC concentrations in the Pleistocene yedoma area of the catchment compared to the Holocene non-yedoma area with medians of 5 and 4.5 mg L-1 (p < 0.05), respectively. When yedoma thaw streams with high DOC concentration reach a large yedoma thermokarst lake, we observe an abrupt decrease in DOC concentration, which we attribute to dilution and lake processes such as mineralization. The DOC ages in the large thermokarst lake (between 3,428 and 3,637 C-14 y BP) can be attributed to a mixing of mobilized old yedoma and Holocene carbon.
Further downstream after the large thermokarst lake, we find progressively younger DOC ages in the stream water to its mouth, paired with decreasing DOC concentrations. This process could result from dilution with leaching water from Holocene deposits and/or emission of ancient yedoma carbon to the atmosphere. Our study shows that thermokarst lakes and ponds may act as DOC filters, predominantly by diluting incoming waters of higher DOC concentrations or by re-mineralizing DOC to CO2 and CH4.
Nevertheless, our results also confirm that the small catchment still contributes DOC on the order of 1.2 kg km(-2) per day from a permafrost landscape with ice-rich yedoma deposits to the Lena River.
With the advent of the two Sentinel-1 (S1) satellites, Synthetic Aperture Radar (SAR) data with high temporal and spatial resolution are freely available. This provides a promising framework to facilitate detailed investigations of surface instabilities and movements on large scales with high temporal resolution, but also poses substantial processing challenges because of storage and computation requirements. Methods are needed to efficiently detect short term changes in dynamic environments. Approaches considering pair-wise processing of a series of consecutive scenes to retain maximum temporal resolution in conjunction with time series analyses are required. Here we present OSARIS, the “Open Source SAR Investigation System,” as a framework to process large stacks of S1 data on high-performance computing clusters. Based on Generic Mapping Tools SAR, shell scripts, and the workload manager Slurm, OSARIS provides an open and modular framework combining parallelization of high-performance C programs, flexible processing schemes, convenient configuration, and generation of geocoded stacks of analysis-ready base data, including amplitude, phase, coherence, and unwrapped interferograms. Time series analyses can be conducted by applying automated modules to the data stacks. The capabilities of OSARIS are demonstrated in a case study from the northwestern Tien Shan, Central Asia. After merging of slices, a total of 80 scene pairs were processed from 174 total input scenes. The coherence time series exhibits pronounced seasonal variability, with relatively high coherence values prevailing during the summer months in the nival zone. As an example of a time series analysis module, we present OSARIS' “Unstable Coherence Metric” which identifies pixels affected by significant drops from high to low coherence values. Measurements of motion provided by LOSD measurements require careful evaluation because interferometric phase unwrapping is prone to errors. Here, OSARIS provides a series of modules to detect and mask unwrapping errors, correct for atmospheric disturbances, and remove large-scale trends. Wall clock processing time for the case study (area ~9,000 km2) was ~12 h 4 min on a machine with 400 cores and 2 TB RAM. In total, ~12 d 10 h 44 min (~96%) were saved through parallelization. A comparison of selected OSARIS datasets to results from two state-of-the-art SAR processing suites, ISCE and SNAP, shows that OSARIS provides products of competitive quality despite its high level of automatization. OSARIS thus facilitates efficient S1-based region-wide investigations of surface movement events over multiple years.
Rapid local adaptation frequently occurs during the spread of invading species. It remains unclear, however, how consistent, and therefore potentially predictable, such patterns of local adaptation are. One approach to this question is to measure patterns of local differentiation in functional traits and plasticity levels in invasive species in multiple regions. Finding consistent patterns of local differentiation in replicate regions suggests that these patterns are adaptive. Further, this outcome indicates that the invading species likely responds predictably to selection along environmental gradients, even though standing genetic variation is likely to have been reduced during introduction. We studied local differentiation in the invasive annual plant Erodium cicutarium in two invaded regions, California and Chile. We collected seeds from across strong gradients in precipitation and temperature in Mediterranean-climate parts of the two regions (10 populations per region). We grew seeds from maternal families from these populations through two generations and exposed the second generation to contrasting levels of water and nutrient availability. We measured growth, flowering time and leaf functional traits across these treatments to obtain trait means and plasticity measures. We found strong differentiation among populations in all traits. Plants from drier environments flowered earlier, were less plastic in flowering time and reached greater size in all treatments. Correlations among traits within regions suggested a coordinated evolutionary response along environmental gradients associated with growing season length. There was little divergence in traits and trait intercorrelations between regions, but strongly parallel divergence in traits within regions. Similar, statistically consistent patterns of local trait differentiation across two regions suggest that local adaptation to environmental gradients has aided the spread of this invasive species, and that the formation of ecotypes in newly invaded environments has been relatively consistent and predictable.
Past climate and continentality inferred from ice wedges at Batagay Highlands, interior Yakutia
(2019)
Ice wedges in the Yana Highlands of interior Yakutia - the most continental region of the Northern Hemisphere - were investigated to elucidate changes in winter climate and continentality that have taken place since the Middle Pleistocene. The Batagay megaslump exposes ice wedges and composite wedges that were sampled from three cryostratigraphic units: the lower ice complex of likely pre-Marine Isotope Stage (MIS) 6 age, the upper ice complex (Yedoma) and the upper sand unit (both MIS 3 to 2). A terrace of the nearby Adycha River provides a Late Holocene (MIS 1) ice wedge that serves as a modern reference for interpretation. The stable-isotope composition of ice wedges in the MIS 3 upper ice complex at Batagay is more depleted (mean delta O-18 about -35 parts per thousand) than those from 17 other ice-wedge study sites across coastal and central Yakutia. This observation points to lower winter temperatures and therefore higher continentality in the Yana Highlands during MIS 3. Likewise, more depleted isotope values are found in Holocene wedge ice (mean delta O-18 about -29 parts per thousand) compared to other sites in Yakutia. Ice-wedge isotopic signatures of the lower ice complex mean delta O-18 about -33 parts per thousand) and of the MIS 3-2 upper sand unit (mean delta O-18 from about -33 parts per thousand to -30 parts per thousand) are less distinctive regionally. The latter unit preserves traces of fast formation in rapidly accumulating sand sheets and of post-depositional isotopic fractionation.
James Ross Island (JRI) offers the exceptional opportunity to study microbial-driven pedogenesis without the influence of vascular plants or faunal activities (e.g., penguin rookeries). In this study, two soil profiles from JRI (one at Santa Martha Cove - SMC, and another at Brandy Bay BB) were investigated, in order to gain information about the initial state of soil formation and its interplay with prokaryotic activity, by combining pedological, geochemical and microbiological methods. The soil profiles are similar with respect to topographic position and parent material but are spatially separated by an orographic barrier and therefore represent windward and leeward locations towards the mainly southwesterly winds. These different positions result in differences in electric conductivity of the soils caused by additional input of bases by sea spray at the windward site and opposing trends in the depth functions of soil pH and electric conductivity. Both soils are classified as Cryosols, dominated by bacterial taxa such as Actinobacteria, Proteobacteria, Acidobacteria, Gemmatimonadetes and Chloroflexi. A shift in the dominant taxa was observed below 20 cm in both soils as well as an increased abundance of multiple operational taxonomic units (OTUs) related to potential chemolithoautotrophic Acidiferrobacteraceae. This shift is coupled by a change in microstructure. While single/pellicular grain microstructure (SMC) and platy microstructure (BB) are dominant above 20 cm, lenticular microstructure is dominant below 20 cm in both soils. The change in microstructure is caused by frequent freeze-thaw cycles and a relative high water content, and it goes along with a development of the pore spacing and is accompanied by a change in nutrient content. Multivariate statistics revealed the influence of soil parameters such as chloride, sulfate, calcium and organic carbon contents, grain size distribution and pedogenic oxide ratios on the overall microbial community structure and explained 49.9% of its variation. The correlation of the pedogenic oxide ratios with the compositional distribution of microorganisms as well as the relative abundance certain microorganisms such as potentially chemolithotrophic Acidiferrobacteraceae-related OTUs could hint at an interplay between soil-forming processes and microorganisms.
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.
In-depth understanding of the reorganization of the hydrological cycle in response to global climate change is crucial in highly sensitive regions like the eastern Mediterranean, where water availability is a major factor for socioeconomic and political development.
The sediments of Lake Lisan provide a unique record of hydroclimatic change during the last glacial to Holocene transition (ca. 24-11 ka) with its tremendous water level drop of similar to 240 m that finally led to its transition into the present hypersaline water body-the Dead Sea.
Here we utilize high-resolution sedimentological analyses from the marginal terraces and deep lake to reconstruct an unprecedented seasonal record of the last millennia of Lake Lisan. Aragonite varve formation in intercalated intervals of our record demonstrates that a stepwise long-term lake level decline was interrupted by almost one millennium of rising or stable water level.
Even periods of pronounced water level drops indicated by gypsum deposition were interrupted by decades of positive water budgets.
Our results thus highlight that even during major climate change at the end of the last glacial, decadal to millennial periods of relatively stable or positive moisture supply occurred which could have been an important premise for human sedentism.
Woody plants are expanding into the Arctic in response to the warming climate. The impact on arctic plant communities is not well understood due to the limited knowledge about plant assembly rules.
Records of past plant diversity over long time series are rare. Here, we applied sedimentary ancient DNA metabarcoding targeting the P6 loop of the chloroplast trnL gene to a sediment record from Lake Ilirney (central Chukotka, Far Eastern Russia) covering the last 28 thousand years.
Our results show that forb-rich steppe-tundra and dwarf-shrub tundra dominated during the cold climate before 14 ka, while deciduous erect-shrub tundra was abundant during the warm period since 14 ka. Larix invasion during the late Holocene substantially lagged behind the likely warmest period between 10 and 6 ka, where the vegetation biomass could be highest.
We reveal highest richness during 28-23 ka and a second richness peak during 13-9 ka, with both periods being accompanied by low relative abundance of shrubs. During the cold period before 14 ka, rich plant assemblages were phylogenetically clustered, suggesting low genetic divergence in the assemblages despite the great number of species. This probably originates from environmental filtering along with niche differentiation due to limited resources under harsh environmental conditions. In contrast, during the warmer period after 14 ka, rich plant assemblages were phylogenetically overdispersed.
This results from a high number of species which were found to harbor high genetic divergence, likely originating from an erratic recruitment process in the course of warming. Some of our evidence may be of relevance for inferring future arctic plant assembly rules and diversity changes. By analogy to the past, we expect a lagged response of tree invasion. Plant richness might overshoot in the short term; in the long-term, however, the ongoing expansion of deciduous shrubs will eventually result in a phylogenetically more diverse community.
POET (v0.1): speedup of many-core parallel reactive transport simulations with fast DHT lookups
(2021)
Coupled reactive transport simulations are extremely demanding in terms of required computational power, which hampers their application and leads to coarsened and oversimplified domains. The chemical sub-process represents the major bottleneck: its acceleration is an urgent challenge which gathers increasing interdisciplinary interest along with pressing requirements for subsurface utilization such as spent nuclear fuel storage, geothermal energy and CO2 storage. In this context we developed POET (POtsdam rEactive Transport), a research parallel reactive transport simulator integrating algorithmic improvements which decisively speed up coupled simulations. In particular, POET is designed with a master/worker architecture, which ensures computational efficiency in both multicore and cluster compute environments. POET does not rely on contiguous grid partitions for the parallelization of chemistry but forms work packages composed of grid cells distant from each other. Such scattering prevents particularly expensive geochemical simulations, usually concentrated in the vicinity of a reactive front, from generating load imbalance between the available CPUs (central processing units), as is often the case with classical partitions. Furthermore, POET leverages an original implementation of the distributed hash table (DHT) mechanism to cache the results of geochemical simulations for further reuse in subsequent time steps during the coupled simulation. The caching is hence particularly advantageous for initially chemically homogeneous simulations and for smooth reaction fronts. We tune the rounding employed in the DHT on a 2D benchmark to validate the caching approach, and we evaluate the performance gain of POET's master/worker architecture and the DHT speedup on a 3D benchmark comprising around 650 000 grid elements. The runtime for 200 coupling iterations, corresponding to 960 simulation days, reduced from about 24 h on 11 workers to 29 min on 719 workers. Activating the DHT reduces the runtime further to 2 h and 8 min respectively. Only with these kinds of reduced hardware requirements and computational costs is it possible to realistically perform the longterm complex reactive transport simulations, as well as perform the uncertainty analyses required by pressing societal challenges connected with subsurface utilization.
As the recent permafrost thawing of northern Asia proceeds due to anthropogenic climate change, precise and detailed palaeoecological records from past warm periods are essential to anticipate the extent of future permafrost variations. Here, based on the modern relationship between permafrost and vegetation (represented by pollen assemblages), we trained a Random Forest model using pollen and permafrost data and verified its reliability to reconstruct the history of permafrost in northern Asia during the Holocene. An early Holocene (12-8 cal ka BP) strong thawing trend, a middle-to-late Holocene (8-2 cal ka BP) relatively slow thawing trend, and a late Holocene freezing trend of permafrost in northern Asia are consistent with climatic proxies such as summer solar radiation and Northern Hemisphere temperature. The extensive distribution of permafrost in northern Asia inhibited the spread of evergreen coniferous trees during the early Holocene warming and might have decelerated the enhancement of the East Asian summer monsoon (EASM) by altering hydrological processes and albedo. Based on these findings, we suggest that studies of the EASM should consider more the state of permafrost and vegetation in northern Asia, which are often overlooked and may have a profound impact on climate change in this region.
Pollen-based quantitative land-cover reconstruction for northern Asia covering the last 40 ka cal BP
(2019)
We collected the available relative pollen productivity estimates (PPEs) for 27 major pollen taxa from Eurasia and applied them to estimate plant abundances during the last 40 ka cal BP (calibrated thousand years before present) using pollen counts from 203 fossil pollen records in northern Asia (north of 40 degrees N). These pollen records were organized into 42 site groups and regional mean plant abundances calculated using the REVEALS (Regional Estimates of Vegetation Abundance from Large Sites) model. Time-series clustering, constrained hierarchical clustering, and detrended canonical correspondence analysis were performed to investigate the regional pattern, time, and strength of vegetation changes, respectively. Reconstructed regional plant functional type (PFT) components for each site group are generally consistent with modern vegetation in that vegetation changes within the regions are characterized by minor changes in the abundance of PFTs rather than by an increase in new PFTs, particularly during the Holocene. We argue that pollen-based REVEALS estimates of plant abundances should be a more reliable reflection of the vegetation as pollen may overestimate the turnover, particularly when a high pollen producer invades areas dominated by low pollen producers. Comparisons with vegetation-independent climate records show that climate change is the primary factor driving land-cover changes at broad spatial and temporal scales. Vegetation changes in certain regions or periods, however, could not be explained by direct climate change, e.g. inland Siberia, where a sharp increase in evergreen conifer tree abundance occurred at ca. 7-8 ka cal BP despite an unchanging climate, potentially reflecting their response to complex climate-permafrost-fire-vegetation interactions and thus a possible long-term lagged climate response.
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.
Changes in snowpack associated with climatic warming has drastic impacts on surface energy balance in the cryosphere. Yet, traditional monitoring techniques, such as punctual measurements in the field, do not cover the full snowpack spatial and temporal variability, which hampers efforts to upscale measurements to the global scale. This variability is one of the primary constraints in model development. In terms of spatial resolution, active microwaves (synthetic aperture radar - SAR) can address the issue and outperform methods based on passive microwaves. Thus, high-spatial-resolution monitoring of snow depth (SD) would allow for better parameterization of local processes that drive the spatial variability of snow. The overall objective of this study is to evaluate the potential of the TerraSAR-X (TSX) SAR sensor and the wave co-polar phase difference (CPD) method for characterizing snow cover at high spatial resolution. Consequently, we first (1) investigate SD and depth hoar fraction (DHF) variability between different vegetation classes in the Ice Creek catchment (Qikiqtaruk/Herschel Island, Yukon, Canada) using in situ measurements collected over the course of a field campaign in 2019; (2) evaluate linkages between snow characteristics and CPD distribution over the 2019 dataset; and (3) determine CPD seasonality considering meteorological data over the 2015-2019 period. SD could be extracted using the CPD when certain conditions are met. A high incidence angle (>30 circle) with a high topographic wetness index (TWI) (>7.0) showed correlation between SD and CPD (R2 up to 0.72). Further, future work should address a threshold of sensitivity to TWI and incidence angle to map snow depth in such environments and assess the potential of using interpolation tools to fill in gaps in SD information on drier vegetation types.
Transport properties of potential host rocks for nuclear waste disposal are typically determined in laboratory or in-situ experiments under geochemically controlled and constant conditions. Such a homogeneous assumption is no longer applicable on the host rock scale as can be seen from the pore water profiles of the potential host rock Opalinus Clay at Mont Terri (Switzerland). The embedding aquifers are the hydro-geological boundaries, that established gradients in the 210 m thick low permeable section through diffusive exchange over millions of years. Present-day pore water profiles were confirmed by a data-driven as well as by a conceptual scenario. Based on the modelled profiles, the influence of the geochemical gradient on uranium migration was quantified by comparing the distances after one million years with results of common homogeneous models. Considering the heterogeneous system, uranium migrated up to 24 m farther through the formation depending on the source term position within the gradient and on the partial pressure of carbon dioxide pCO2 of the system. Migration lengths were almost equal for single- and multicomponent diffusion. Differences can predominantly be attributed to changes in the sorption capacity, whereby pCO2 governs how strong uranium migration is affected by the geochemical gradient. Thus, the governing parameters for uranium migration in the Opalinus Clay can be ordered in descending priority: pCO2, geochemical gradients, mineralogical heterogeneity.</p>
The subsurface is a temporally dynamic and spatially heterogeneous compartment of the Earth's critical zone, and biogeochemical transformations taking place in this compartment are crucial for the cycling of nutrients.
The impact of spatial heterogeneity on such microbially mediated nutrient cycling is not well known, which imposes a severe challenge in the prediction of in situ biogeochemical transformation rates and further of nutrient loading contributed by the groundwater to the surface water bodies.
Therefore, we used a numerical modelling approach to evaluate the sensitivity of groundwater microbial biomass distribution and nutrient cycling to spatial heterogeneity in different scenarios accounting for various residence times.
The model results gave us an insight into domain characteristics with respect to the presence of oxic niches in predominantly anoxic zones and vice versa depending on the extent of spatial heterogeneity and the flow regime.
The obtained results show that microbial abundance, distribution, and activity are sensitive to the applied flow regime and that the mobile (i.e. observable by groundwater sampling) fraction of microbial biomass is a varying, yet only a small, fraction of the total biomass in a domain. Furthermore, spatial heterogeneity resulted in anaerobic niches in the domain and shifts in microbial biomass between active and inactive states. The lack of consideration of spatial heterogeneity, thus, can result in inaccurate estimation of microbial activity. In most cases this leads to an overestimation of nutrient removal (up to twice the actual amount) along a flow path.
We conclude that the governing factors for evaluating this are the residence time of solutes and the Damkohler number (Da) of the biogeochemical reactions in the domain. We propose a relationship to scale the impact of spatial heterogeneity on nutrient removal governed by the logioDa.
This relationship may be applied in upscaled descriptions of microbially mediated nutrient cycling dynamics in the subsurface thereby resulting in more accurate predictions of, for example, carbon and nitrogen cycling in groundwater over long periods at the catchment scale.
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.
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.
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.
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 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.
Knowledge on the response of sediment export to recent climate change in glacierized areas in the European Alps is limited, primarily because long-term records of suspended sediment concentrations (SSCs) are scarce. Here we tested the estimation of sediment export of the past five decades using quantile regression forest (QRF), a nonparametric, multivariate regression based on random forest. The regression builds on short-term records of SSCs and long records of the most important hydroclimatic drivers (discharge, precipitation and air temperature - QPT). We trained independent models for two nested and partially glacier-covered catchments, Vent (98 km(2)) and Vernagt (11.4 km(2)), in the upper otztal in Tyrol, Austria (1891 to 3772 m a.s.l.), where available QPT records start in 1967 and 1975. To assess temporal extrapolation ability, we used two 2-year SSC datasets at gauge Vernagt, which are almost 20 years apart, for a validation. For Vent, we performed a five-fold cross-validation on the 15 years of SSC measurements. Further, we quantified the number of days where predictors exceeded the range represented in the training dataset, as the inability to extrapolate beyond this range is a known limitation of QRF. Finally, we compared QRF performance to sediment rating curves (SRCs). We analyzed the modeled sediment export time series, the predictors and glacier mass balance data for trends (Mann-Kendall test and Sen's slope estimator) and step-like changes (using the widely applied Pettitt test and a complementary Bayesian approach).Our validation at gauge Vernagt demonstrated that QRF performs well in estimating past daily sediment export (Nash-Sutcliffe efficiency (NSE) of 0.73) and satisfactorily for SSCs (NSE of 0.51), despite the small training dataset. The temporal extrapolation ability of QRF was superior to SRCs, especially in periods with high-SSC events, which demonstrated the ability of QRF to model threshold effects. Days with high SSCs tended to be underestimated, but the effect on annual yields was small. Days with predictor exceedances were rare, indicating a good representativity of the training dataset. Finally, the QRF reconstruction models outperformed SRCs by about 20 percent points of the explained variance.Significant positive trends in the reconstructed annual suspended sediment yields were found at both gauges, with distinct step-like increases around 1981. This was linked to increased glacier melt, which became apparent through step-like increases in discharge at both gauges as well as change points in mass balances of the two largest glaciers in the Vent catchment. We identified exceptionally high July temperatures in 1982 and 1983 as a likely cause. In contrast, we did not find coinciding change points in precipitation. Opposing trends at the two gauges after 1981 suggest different timings of "peak sediment". We conclude that, given large-enough training datasets, the presented QRF approach is a promising tool with the ability to deepen our understanding of the response of high-alpine areas to decadal climate change.
Arctic tundra ecosystems have experienced unprecedented change associated with climate warming over recent decades. Across the Pan-Arctic, vegetation productivity and surface greenness have trended positively over the period of satellite observation. However, since 2011 these trends have slowed considerably, showing signs of browning in many regions. It is unclear what factors are driving this change and which regions/landforms will be most sensitive to future browning. Here we provide evidence linking decadal patterns in arctic greening and browning with regional climate change and local permafrost-driven landscape heterogeneity. We analyzed the spatial variability of decadal-scale trends in surface greenness across the Arctic Coastal Plain of northern Alaska (similar to 60,000 km(2)) using the Landsat archive (1999-2014), in combination with novel 30 m classifications of polygonal tundra and regional watersheds, finding landscape heterogeneity and regional climate change to be the most important factors controlling historical greenness trends. Browning was linked to increased temperature and precipitation, with the exception of young landforms (developed following lake drainage), which will likely continue to green. Spatiotemporal model forecasting suggests carbon uptake potential to be reduced in response to warmer and/or wetter climatic conditions, potentially increasing the net loss of carbon to the atmosphere, at a greater degree than previously expected.
Colombia's agriculture, forestry and other land use sector accounts for nearly half of its total greenhouse gas (GHG) emissions. The importance of smallholder deforestation is comparatively high in relation to its regional counterparts, and livestock agriculture represents the largest driver of primary forest depletion. Silvopastoral systems (SPSs) are presented as agroecological solutions that synergistically enhance livestock productivity, improve local farmers' livelihoods and hold the potential to reduce pressure on forest conversion. The department of Caquetá represents Colombia's most important deforestation hotspot. Targeting smallholder livestock farms through survey data, in this work we investigate the GHG mitigation potential of implementing SPSs for smallholder farms in this region. Specifically, we assess whether the carbon sequestration taking place in the soil and biomass of SPSs is sufficient to offset the per-hectare increase in livestock GHG emissions resulting from higher stocking rates. To address these questions we use data on livestock population characteristics and historic land cover changes reported from a survey covering 158 farms and model the carbon sequestration occurring in three different scenarios of progressively-increased SPS complexity using the CO2 fix model. We find that, even with moderate tree planting densities, the implementation of SPSs can reduce GHG emissions by 2.6 Mg CO2e ha−1 yr−1 in relation to current practices, while increasing agriculture productivity and contributing to the restoration of severely degraded landscapes.
Nature-based solutions (NBS) have recently received attention due to their potential ability to sustainably reduce hydro-meteorological risks, providing co-benefits for both ecosystems and affected people. Therefore, pioneering research has dedicated efforts to optimize the design of NBS, to evaluate their wider co-benefits and to understand promoting and/or hampering governance conditions for the uptake of NBS. In this article, we aim to complement this research by conducting a comprehensive literature review of factors shaping people’s perceptions of NBS as a means to reduce hydro-meteorological risks. Based on 102 studies, we identified six topics shaping the current discussion in this field of research: (1) valuation of the co-benefits (including those related to ecosystems and society); (2) evaluation of risk reduction efficacy; (3) stakeholder participation; (4) socio-economic and location-specific conditions; (5) environmental attitude, and (6) uncertainty. Our analysis reveals that concerned empirical insights are diverse and even contradictory, they vary in the depth of the insights generated and are often not comparable for a lack of a sound theoretical-methodological grounding. We, therefore, propose a conceptual model outlining avenues for future research by indicating potential inter-linkages between constructs underlying perceptions of NBS to hydro-meteorological risks.
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.
Regional flood frequency analysis (RFFA) is a powerful method for interrogating hydrological series since it combines observational time series from several sites within a region to estimate risk-relevant statistical parameters with higher accuracy than from single-site series. Since RFFA extreme value estimates depend on the shape of the selected distribution of the data-generating stochastic process, there is need for a suitable goodness-of-distributional-fit measure in order to optimally utilize given data. Here we present a novel, least-squares-based measure to select the optimal fit from a set of five distributions, namely Generalized Extreme Value (GEV), Generalized Logistic, Gumbel, Log-Normal Type III and Log-Pearson Type III. The fit metric is applied to annual maximum discharge series from six hydrological stations along the Sava River in South-eastern Europe, spanning the years 1961 to 2020. Results reveal that (1) the Sava River basin can be assessed as hydrologically homogeneous and (2) the GEV distribution provides typically the best fit. We offer hydrological-meteorological insights into the differences among the six stations. For the period studied, almost all stations exhibit statistically insignificant trends, which renders the conclusions about flood risk as relevant for hydrological sciences and the design of regional flood protection infrastructure.
Many widely used observational data sets are comprised of several overlapping instrument records. While data inter-calibration techniques often yield continuous and reliable data for trend analysis, less attention is generally paid to maintaining higher-order statistics such as variance and autocorrelation. A growing body of work uses these metrics to quantify the stability or resilience of a system under study and potentially to anticipate an approaching critical transition in the system. Exploring the degree to which changes in resilience indicators such as the variance or autocorrelation can be attributed to non-stationary characteristics of the measurement process – rather than actual changes in the dynamical properties of the system – is important in this context. In this work we use both synthetic and empirical data to explore how changes in the noise structure of a data set are propagated into the commonly used resilience metrics lag-one autocorrelation and variance. We focus on examples from remotely sensed vegetation indicators such as vegetation optical depth and the normalized difference vegetation index from different satellite sources. We find that time series resulting from mixing signals from sensors with varied uncertainties and covering overlapping time spans can lead to biases in inferred resilience changes. These biases are typically more pronounced when resilience metrics are aggregated (for example, by land-cover type or region), whereas estimates for individual time series remain reliable at reasonable sensor signal-to-noise ratios. Our work provides guidelines for the treatment and aggregation of multi-instrument data in studies of critical transitions and resilience.
In this study, we analyzed a large seismological dataset from temporary and permanent networks in the southern and eastern Alps to establish high-precision hypocenters and 1-D V-P and V-P/V-S models. The waveform data of a subset of local earthquakes with magnitudes in the range of 1-4.2 M-L were recorded by the dense, temporary SWATH-D network and selected stations of the AlpArray network between September 2017 and the end of 2018. The first arrival times of P and S waves of earthquakes are determined by a semi-automatic procedure. We applied a Markov chain Monte Carlo inversion method to simultaneously calculate robust hypocenters, a 1-D velocity model, and station corrections without prior assumptions, such as initial velocity models or earthquake locations. A further advantage of this method is the derivation of the model parameter uncertainties and noise levels of the data. The precision estimates of the localization procedure is checked by inverting a synthetic travel time dataset from a complex 3-D velocity model and by using the real stations and earthquakes geometry. The location accuracy is further investigated by a quarry blast test. The average uncertainties of the locations of the earthquakes are below 500m in their epicenter and similar to 1.7 km in depth. The earthquake distribution reveals seismicity in the upper crust (0-20 km), which is characterized by pronounced clusters along the Alpine frontal thrust, e.g., the Friuli-Venetia (FV) region, the Giudicarie-Lessini (GL) and Schio-Vicenza domains, the Austroalpine nappes, and the Inntal area. Some seismicity also occurs along the Periadriatic Fault. The general pattern of seismicity reflects head-on convergence of the Adriatic indenter with the Alpine orogenic crust. The seismicity in the FV and GL regions is deeper than the modeled frontal thrusts, which we interpret as indication for southward propagation of the southern Alpine deformation front (blind thrusts).
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.
The hyperthermal events of the Cenozoic, including the Paleocene-Eocene Thermal Maximum, provide an opportunity to investigate the potential effects of climate warming on marine ecosystems. Here, we examine the shallow benthic marine communities preserved in the late Cretaceous to Eocene strata on the Gulf Coastal Plain (United States). In stark contrast to the ecological shifts following the end-Cretaceous mass extinction, our data show that the early Cenozoic hyperthermals did not have a long-term impact on the generic diversity nor composition of the Gulf Coastal Plain molluscan communities. We propose that these communities were resilient to climate change because molluscs are better adapted to high temperatures than other taxa, as demonstrated by their physiology and evolutionary history. In terms of resilience, these communities differ from other shallow-water carbonate ecosystems, such as reef communities, which record significant changes during the early Cenozoic hyperthermals. These data highlight the strikingly different responses of community types, i.e., the almost imperceptible response of molluscs versus the marked turnover of foraminifera and reef faunas. The impact on molluscan communities may have been low because detrimental conditions did not devastate the entire Gulf Coastal Plain, allowing molluscs to rapidly recolonise vacated areas once harsh environmental conditions ameliorated.
The imagination of clearly separated core-shell structures is already outdated by the fact, that the nanoparticle core-shell structures remain in terms of efficiency behind their respective bulk material due to intermixing between core and shell dopant ions. In order to optimize the photoluminescence of core-shell UCNP the intermixing should be as small as possible and therefore, key parameters of this process need to be identified. In the present work the Ln(III) ion migration in the host lattices NaYF4 and NaGdF4 was monitored. These investigations have been performed by laser spectroscopy with help of lanthanide resonance energy transfer (LRET) between Eu(III) as donor and Pr(III) or Nd(III) as acceptor. The LRET is evaluated based on the Forster theory. The findings corroborate the literature and point out the migration of ions in the host lattices. Based on the introduced LRET model, the acceptor concentration in the surrounding of one donor depends clearly on the design of the applied core-shell-shell nanoparticles. In general, thinner intermediate insulating shells lead to higher acceptor concentration, stronger quenching of the Eu(III) donor and subsequently stronger sensitization of the Pr(III) or the Nd(III) acceptors. The choice of the host lattice as well as of the synthesis temperature are parameters to be considered for the intermixing process.
The Sentinel Application Platform (SNAP) architecture facilitates Earth Observation data processing. In this work, we present results from a new Snow Processor for SNAP. We also describe physical principles behind the developed snow property retrieval technique based on the analysis of Ocean and Land Colour Instrument (OLCI) onboard Sentinel-3A/B measurements over clean and polluted snow fields. Using OLCI spectral reflectance measurements in the range 400-1020 nm, we derived important snow properties such as spectral and broadband albedo, snow specific surface area, snow extent and grain size on a spatial grid of 300 m. The algorithm also incorporated cloud screening and atmospheric correction procedures over snow surfaces. We present validation results using ground measurements from Antarctica, the Greenland ice sheet and the French Alps. We find the spectral albedo retrieved with accuracy of better than 3% on average, making our retrievals sufficient for a variety of applications. Broadband albedo is retrieved with the average accuracy of about 5% over snow. Therefore, the uncertainties of satellite retrievals are close to experimental errors of ground measurements. The retrieved surface grain size shows good agreement with ground observations. Snow specific surface area observations are also consistent with our OLCI retrievals. We present snow albedo and grain size mapping over the inland ice sheet of Greenland for areas including dry snow, melted/melting snow and impurity rich bare ice. The algorithm can be applied to OLCI Sentinel-3 measurements providing an opportunity for creation of long-term snow property records essential for climate monitoring and data assimilation studies-especially in the Arctic region, where we face rapid environmental changes including reduction of snow/ice extent and, therefore, planetary albedo.
Rupture directivity, implying a predominant earthquake rupture propagation direction, is typically inferred upon the identification of 2D azimuthal patterns of seismic observations for weak to large earthquakes using surface-monitoring networks. However, the recent increase of 3D monitoring networks deployed in the shallow subsurface and underground laboratories toward the monitoring of microseismicity allows to extend the directivity analysis to 3D modeling, beyond the usual range of magnitudes. The high-quality full waveforms recorded for the largest, decimeter-scale acoustic emission (AE) events during a meter-scale hydraulic fracturing experiment in granites at similar to 410 m depth allow us to resolve the apparent durations observed at each AE sensor to analyze 3D-directivity effects. Unilateral and (asymmetric) bilateral ruptures are then characterized by the introduction of a parameter kappa, representing the angle between the directivity vector and the station vector. While the cloud of AE activity indicates the planes of the hydrofractures, the resolved directivity vectors show off-plane orientations, indicating that rupture planes of microfractures on a scale of centimeters have different geometries. Our results reveal a general alignment of the rupture directivity with the orientation of the minimum horizontal stress, implying that not only the slip direction but also the fracture growth produced by the fluid injections is controlled by the local stress conditions.
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.
At the interface between the lithosphere and the atmosphere, the critical zone records the complex interactions between erosion, climate, geologic substrate, and life and can be directly monitored. Long data records (30 consecutive years for sediment yields) collected in the sparsely vegetated, steep, and small marly badland catchments of the Draix-Bleone Critical Zone Observatory (CZO), SE France, allow analyzing potential climatic controls on regolith dynamics and sediment export. Although widely accepted as a first-order control, rainfall variability does not fully explain the observed interannual variability in sediment export. Previous studies in this area have suggested that frost-weathering processes could drive regolith production and potentially modulate the observed pattern of sediment export. Here, we define sediment export anomalies as the residuals from a predictive model with annual rainfall intensity above a threshold as the control. We then use continuous soil temperature data recorded at different locations over multiple years to highlight the role of different frost-weathering processes (i.e., ice segregation versus volumetric expansion) in regolith production. Several proxies for different frost-weathering processes have been calculated from these data and compared to the sediment export anomalies, with careful consideration of field data quality. Our results suggest that frost-cracking intensity (linked to ice segregation) can explain about half (47 %-64 %) of the sediment export anomalies. In contrast, the number of freeze-thaw cycles (linked to volumetric expansion) has only a minor impact on catchment sediment response. The time spent below 0 degrees C also correlates well with the sediment export anomalies and requires fewer field data to be calculated than the frost-cracking intensity. Thus, frost-weathering processes modulate sediment export by controlling regolith production in these catchments and should be taken into account when building predictive models of sediment export from these badlands under a changing climate.