Gold Open-Access
Refine
Document Type
- Article (200)
- Review (4)
- Other (2)
- Doctoral Thesis (1)
Is part of the Bibliography
- yes (207)
Keywords
- climate change (8)
- permafrost (7)
- 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 (207) (remove)
White mica and tourmaline are the dominant hydrothermal alteration minerals at the world-class Panasqueira W-Sn-Cu deposit in Portugal. Thus, understanding the controls on their chemical composition helps to constrain ore formation processes at this deposit and determine their usefulness as pathfinder minerals for mineralization in general. We combine whole-rock geochemistry of altered and unaltered metasedimentary host rocks with in situ LA-ICP-MS measurements of tourmaline and white mica from the alteration halo. Principal component analysis (PCA) is used to better identify geochemical patterns and trends of hydrothermal alteration in the datasets. The hydrothermally altered metasediments are enriched in As, Sn, Cs, Li, W, F, Cu, Rb, Zn, Tl, and Pb relative to unaltered samples. In situ mineral analyses show that most of these elements preferentially partition into white mica over tourmaline (Li, Rb, Cs, Tl, W, and Sn), whereas Zn is enriched in tourmaline. White mica has distinct compositions in different settings within the deposit (greisen, vein selvages, wall rock alteration zone, late fault zone), indicating a compositional evolution with time. In contrast, tourmaline from different settings overlaps in composition, which is ascribed to a stronger dependence on host rock composition and also to the effects of chemical zoning and microinclusions affecting the LA-ICP-MS analyses. Hence, in this deposit, white mica is the better recorder of the fluid composition. The calculated trace-element contents of the Panasqueira mineralizing fluid based on the mica data and estimates of mica-fluid partition coefficients are in good agreement with previous fluid-inclusion analyses. A compilation of mica and tourmaline trace-element compositions from Panasqueira and other W-Sn deposits shows that white mica has good potential as a pathfinder mineral, with characteristically high Li, Cs, Rb, Sn, and W contents. The trace-element contents of hydrothermal tourmaline are more variable. Nevertheless, the compiled data suggest that high Sn and Li contents are distinctive for tourmaline from W-Sn deposits.
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.
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.
Cyanobacteria are important primary producers in temperate freshwater ecosystems. However, studies on the seasonal and spatial distribution of cyanobacteria in deep lakes based on high-throughput DNA sequencing are still rare. In this study, we combined monthly water sampling and monitoring in 2019, amplicon sequence variants analysis (ASVs; a proxy for different species) and quantitative PCR targeting overall cyanobacteria abundance to describe the seasonal and spatial dynamics of cyanobacteria in the deep hard-water oligo-mesotrophic Lake Tiefer See, NE Germany. We observed significant seasonal variation in the cyanobacterial community composition (p < 0.05) in the epi- and metalimnion layers, but not in the hypolimnion. In winter-when the water column is mixed-picocyanobacteria (Synechococcus and Cyanobium) were dominant. With the onset of stratification in late spring, we observed potential niche specialization and coexistence among the cyanobacteria taxa driven mainly by light and nutrient dynamics. Specifically, ASVs assigned to picocyanobacteria and the genus Planktothrix were the main contributors to the formation of deep chlorophyll maxima along a light gradient. While Synechococcus and different Cyanobium ASVs were abundant in the epilimnion up to the base of the euphotic zone from spring to fall, Planktothrix mainly occurred in the metalimnetic layer below the euphotic zone where also overall cyanobacteria abundance was highest in summer. Our data revealed two potentially psychrotolerant (cold-adapted) Cyanobium species that appear to cope well under conditions of lower hypolimnetic water temperature and light as well as increasing sediment-released phosphate in the deeper waters in summer. The potential cold-adapted Cyanobium species were also dominant throughout the water column in fall and winter. Furthermore, Snowella and Microcystis-related ASVs were abundant in the water column during the onset of fall turnover. Altogether, these findings suggest previously unascertained and considerable spatiotemporal changes in the community of cyanobacteria on the species level especially within the genus Cyanobium in deep hard-water temperate lakes.
Strong hydroclimatic controls on vulnerability to subsurface nitrate contamination across Europe
(2020)
Subsurface contamination due to excessive nutrient surpluses is a persistent and widespread problem in agricultural areas across Europe. The vulnerability of a particular location to pollution from reactive solutes, such as nitrate, is determined by the interplay between hydrologic transport and biogeochemical transformations. Current studies on the controls of subsurface vulnerability do not consider the transient behaviour of transport dynamics in the root zone. Here, using state-of-the-art hydrologic simulations driven by observed hydroclimatic forcing, we demonstrate the strong spatiotemporal heterogeneity of hydrologic transport dynamics and reveal that these dynamics are primarily controlled by the hydroclimatic gradient of the aridity index across Europe. Contrasting the space-time dynamics of transport times with reactive timescales of denitrification in soil indicate that similar to 75% of the cultivated areas across Europe are potentially vulnerable to nitrate leaching for at least onethird of the year. We find that neglecting the transient nature of transport and reaction timescale results in a great underestimation of the extent of vulnerable regions by almost 50%. Therefore, future vulnerability and risk assessment studies must account for the transient behaviour of transport and biogeochemical transformation processes.
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.
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.
Nocardioides alcanivorans sp. nov., a novel hexadecane-degrading species isolated from plastic waste
(2022)
Strain NGK65(T), a novel hexadecane degrading, non-motile, Gram-positive, rod-to-coccus shaped, aerobic bacterium, was isolated from plastic polluted soil sampled at a landfill.
Strain NGK65(T) hydrolysed casein, gelatin, urea and was catalase-positive. It optimally grew at 28 degrees C. in 0-1% NaCl and at pH 7.5-8.0. Glycerol, D-glucose, arbutin, aesculin, salicin, potassium 5-ketogluconate. sucrose, acetate, pyruvate and hexadecane were used as sole carbon sources.
The predominant membrane fatty acids were iso-C-16:0 followed by iso-C(17:)0 and C-18:1 omega 9c. The major polar lipids were phosphatidylglycerol, phosphatidylethanolamine, phosphatidylinositol and hydroxyphosphatidylinositol.
The cell-wall peptidoglycan type was A3 gamma, with LL-diaminopimelic acid and glycine as the diagnostic amino acids. MK 8 (H-4) was the predominant menaquinone. Phylogenetic analysis based on 16S rRNA gene sequences indicated that strain NGK65(T) belongs to the genus Nocardioides (phylum Actinobacteria). appearing most closely related to Nocardioides daejeonensis MJ31(T) (98.6%) and Nocardioides dubius KSL-104(T) (98.3%).
The genomic DNA G+C content of strain NGK65(T) was 68.2%.
Strain NGK65(T) and the type strains of species involved in the analysis had average nucleotide identity values of 78.3-71.9% as well as digital DNA-DNA hybridization values between 22.5 and 19.7%, which clearly indicated that the isolate represents a novel species within the genus Nocardioides.
Based on phenotypic and molecular characterization, strain NGK65(T) can clearly be differentiated from its phylogenetic neighbours to establish a novel species, for which the name Nocardioides alcanivorans sp. nov. is proposed.
The type strain is NGK65(T) (=DSM 113112(T)=NCCB 100846(T)).
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.
We present a new set of global and local sea‐level projections at example tide gauge locations under the RCP2.6, RCP4.5, and RCP8.5 emissions scenarios. Compared to the CMIP5‐based sea‐level projections presented in IPCC AR5, we introduce a number of methodological innovations, including (i) more comprehensive treatment of uncertainties, (ii) direct traceability between global and local projections, and (iii) exploratory extended projections to 2300 based on emulation of individual CMIP5 models. Combining the projections with observed tide gauge records, we explore the contribution to total variance that arises from sea‐level variability, different emissions scenarios, and model uncertainty. For the period out to 2300 we further breakdown the model uncertainty by sea‐level component and consider the dependence on geographic location, time horizon, and emissions scenario. Our analysis highlights the importance of local variability for sea‐level change in the coming decades and the potential value of annual‐to‐decadal predictions of local sea‐level change. Projections to 2300 show a substantial degree of committed sea‐level rise under all emissions scenarios considered and highlight the reduced future risk associated with RCP2.6 and RCP4.5 compared to RCP8.5. Tide gauge locations can show large ( > 50%) departures from the global average, in some cases even reversing the sign of the change. While uncertainty in projections of the future Antarctic ice dynamic response tends to dominate post‐2100, we see substantial differences in the breakdown of model variance as a function of location, time scale, and emissions scenario.
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.
The hyperthermal events of the Cenozoic, including the Paleocene-Eocene Thermal Maximum, provide an opportunity to investigate the potential effects of climate warming on marine ecosystems. Here, we examine the shallow benthic marine communities preserved in the late Cretaceous to Eocene strata on the Gulf Coastal Plain (United States). In stark contrast to the ecological shifts following the end-Cretaceous mass extinction, our data show that the early Cenozoic hyperthermals did not have a long-term impact on the generic diversity nor composition of the Gulf Coastal Plain molluscan communities. We propose that these communities were resilient to climate change because molluscs are better adapted to high temperatures than other taxa, as demonstrated by their physiology and evolutionary history. In terms of resilience, these communities differ from other shallow-water carbonate ecosystems, such as reef communities, which record significant changes during the early Cenozoic hyperthermals. These data highlight the strikingly different responses of community types, i.e., the almost imperceptible response of molluscs versus the marked turnover of foraminifera and reef faunas. The impact on molluscan communities may have been low because detrimental conditions did not devastate the entire Gulf Coastal Plain, allowing molluscs to rapidly recolonise vacated areas once harsh environmental conditions ameliorated.
The origin of Asian monsoons
(2020)
The Cenozoic inception and development of the Asian monsoons remain unclear and have generated much debate, as several hypotheses regarding circulation patterns at work in Asia during the Eocene have been proposed in the few last decades. These include (a) the existence of modern-like monsoons since the early Eocene; (b) that of a weak South Asian monsoon (SAM) and little to no East Asian monsoon (EAM); or (c) a prevalence of the Intertropical Convergence Zone (ITCZ) migrations, also referred to as Indonesian-Australian monsoon (I-AM). As SAM and EAM are supposed to have been triggered or enhanced primarily by Asian palaeogeographic changes, their possible inception in the very dynamic Eocene palaeogeographic context remains an open question, both in the modelling and field-based communities. We investigate here Eocene Asian climate conditions using the IPSL-CM5A2 (Sepulchre et al., 2019) earth system model and revised palaeogeographies. Our Eocene climate simulation yields atmospheric circulation patterns in Asia substantially different from modern conditions. A large high-pressure area is simulated over the Tethys ocean, which generates intense low tropospheric winds blowing southward along the western flank of the proto-Himalayan-Tibetan plateau (HTP) system. This low-level wind system blocks, to latitudes lower than 10 degrees N, the migration of humid and warm air masses coming from the Indian Ocean. This strongly contrasts with the modern SAM, during which equatorial air masses reach a latitude of 20-25 degrees N over India and southeastern China. Another specific feature of our Eocene simulation is the widespread subsidence taking place over northern India in the midtroposphere (around 5000 m), preventing deep convective updraught that would transport water vapour up to the condensation level. Both processes lead to the onset of a broad arid region located over northern India and over the HTP. More humid regions of high seasonality in precipitation encircle this arid area, due to the prevalence of the Intertropical Convergence Zone (ITCZ) migrations (or Indonesian-Australian monsoon, I-AM) rather than monsoons. Although the existence of this central arid region may partly result from the specifics of our simulation (model dependence and palaeogeographic uncertainties) and has yet to be confirmed by proxy records, most of the observational evidence for Eocene monsoons are located in the highly seasonal transition zone between the arid area and the more humid surroundings. We thus suggest that a zonal arid climate prevailed over Asia before the initiation of monsoons that most likely occurred following Eocene palaeogeographic changes. Our results also show that precipitation seasonality should be used with caution to infer the presence of a monsoonal circulation and that the collection of new data in this arid area is of paramount importance to allow the debate to move forward.
In an ocean-continent subduction zone, the assessment of the lithospheric thermal state is essential to determine the controls of the deformation within the upper plate and the dip angle of the subducting lithosphere. In this study, we evaluate the degree of influence of both the configuration of the upper plate (i.e., thickness and composition of the rock units) and variations of the subduction angle on the lithospheric thermal field of the southern Central Andes (29 degrees-39 degrees S). Here, the subduction angle increases from subhorizontal (5 degrees) north of 33 degrees S to steep (similar to 30 degrees) in the south. We derived the 3D temperature and heat flow distribution of the lithosphere in the southern Central Andes considering conversion of S wave tomography to temperatures together with steady-state conductive thermal modeling. We found that the orogen is overall warmer than the forearc and the foreland and that the lithosphere of the northern part of the foreland appears colder than its southern counterpart. Sedimentary blanketing and the thickness of the radiogenic crust exert the main control on the shallow thermal field (<50km depth). Specific conditions are present where the oceanic slab is relatively shallow (<85 km depth) and the radiogenic crust is thin. This configuration results in relatively colder temperatures compared to regions where the radiogenic crust is thick and the slab is steep. At depths >50km, the temperatures of the overriding plate are mainly controlled by the mantle heat input and the subduction angle. The thermal field of the upper plate likely preserves the flat subduction angle and influences the spatial distribution of shortening.
New Zealand's Alpine Fault is a large, platebounding strike-slip fault, which ruptures in large (M-w > 8) earthquakes. We conducted field and laboratory analyses of fault rocks to assess its fault zone architecture. Results reveal that the Alpine Fault Zone has a complex geometry, comprising an anastomosing network of multiple slip planes that have accommodated different amounts of displacement. This contrasts with the previous perception of the Alpine Fault Zone, which assumes a single principal slip zone accommodated all displacement. This interpretation is supported by results of drilling projects and geophysical investigations. Furthermore, observations presented here show that the young, largely unconsolidated sediments that constitute the footwall at shallow depths have a significant influence on fault gouge rheological properties and structure.
We produce climate projections through the 21st century using the fractional energy balance equation (FEBE): a generalization of the standard energy balance equation (EBE). The FEBE can be derived from Budyko-Sellers models or phenomenologically through the application of the scaling symmetry to energy storage processes, easily implemented by changing the integer order of the storage (derivative) term in the EBE to a fractional value.
The FEBE is defined by three parameters: a fundamental shape parameter, a timescale and an amplitude, corresponding to, respectively, the scaling exponent h, the relaxation time tau and the equilibrium climate sensitivity (ECS). Two additional parameters were needed for the forcing: an aerosol recalibration factor alpha to account for the large aerosol uncertainty and a volcanic intermittency correction exponent upsilon. A Bayesian framework based on historical temperatures and natural and anthropogenic forcing series was used for parameter estimation. Significantly, the error model was not ad hoc but rather predicted by the model itself: the internal variability response to white noise internal forcing.
The 90 % credible interval (CI) of the exponent and relaxation time were h = [0.33, 0.44] (median = 0.38) and tau = [2.4, 7.0] (median = 4.7) years compared to the usual EBE h = 1, and literature values of tau typically in the range 2-8 years. Aerosol forcings were too strong, requiring a decrease by an average factor alpha = [0.2, 1.0] (median = 0.6); the volcanic intermittency correction exponent was upsilon = [0.15, 0.41] (median = 0.28) compared to standard values alpha = upsilon = 1. The overpowered aerosols support a revision of the global modern (2005) aerosol forcing 90 % CI to a narrower range [ -1.0, -0.2] W m(-2). The key parameter ECS in comparison to IPCC AR5 (and to the CMIP6 MME), the 90 % CI range is reduced from [1.5, 4.5] K ([2.0, 5.5] K) to [1.6, 2.4] K ([1.5, 2.2] K), with median value lowered from 3.0 K (3.7 K) to 2.0 K (1.8 K) Similarly we found for the transient climate response (TCR), the 90 % CI range shrinks from [1.0, 2.5] K ([1.2, 2.8] K) to [1.2, 1.8] K ([1.1, 1.6] K) and the median estimate decreases from 1.8 K (2.0 K) to 1.5 K (1.4 K). As often seen in other observational-based studies, the FEBE values for climate sensitivities are therefore somewhat lower but still consistent with those in IPCC AR5 and the CMIP6 MME. <br /> Using these parameters, we made projections to 2100 using both the Representative Concentration Pathway (RCP) and Shared Socioeconomic Pathway (SSP) scenarios, and compared them to the corresponding CMIP5 and CMIP6 multi-model ensembles (MMEs). The FEBE historical reconstructions (1880-2020) closely follow observations, notably during the 1998-2014 slowdown ("hiatus"). We also reproduce the internal variability with the FEBE and statistically validate this against centennial-scale temperature observations. Overall, the FEBE projections were 10 %-15 % lower but due to their smaller uncertainties, their 90 % CIs lie completely within the GCM 90 % CIs. This agreement means that the FEBE validates the MME, and vice versa.
The effect of lithology on the relationship between denudation rate and chemical weathering pathways
(2022)
The denudation of rocks in mountain belts exposes a range of fresh minerals to the surface of the Earth that are chemically weathered by acidic and oxygenated fluids. The impact of the resulting coupling between denudation and weathering rates fundamentally depends on the types of minerals that are weathering. Whereas silicate weathering sequesters CO2, the combination of sulfide oxidation and carbonate dissolution emits CO2 to the atmosphere. Here, we combine the concentrations of dissolved major elements in stream waters with Be-10 basin-wide denudation rates from 35 small catchments in eastern Tibet to elucidate the importance of lithology in modulating the relationships between denudation rate, chemical weathering pathways, and CO2 consumption or release. Our catchments span 3 orders of magnitude in denudation rate in low-grade flysch, high-grade metapelites, and granitoid rocks. For each stream, we estimate the concentrations of solutes sourced from silicate weathering, carbonate dissolution, and sulfide oxidation using a mixing model. We find that for all lithologies, cation concentrations from silicate weathering are largely independent of denudation rate, but solute concentrations from carbonates and, where present, sulfides increase with increasing denudation rate. With increasing denudation rates, weathering may therefore shift from consuming to releasing CO2 in both (meta)sedimentary and granitoid lithologies. For a given denudation rate, we report dissolved solid concentrations and inferred weathering fluxes in catchments underlain by (meta)sedimentary rock that are 2-10 times higher compared to catchments containing granitoid lithologies, even though climatic and topographic parameters do not vary systematically between these catchments. Thus, varying proportions of exposed (meta)sedimentary and igneous rocks during orogenesis could lead to changes in the sequestration and release of CO2 that are independent of denudation rate.
A review of source models to further the understanding of the seismicity of the Groningen field
(2022)
The occurrence of felt earthquakes due to gas production in Groningen has initiated numerous studies and model attempts to understand and quantify induced seismicity in this region. The whole bandwidth of available models spans the range from fully deterministic models to purely empirical and stochastic models. In this article, we summarise the most important model approaches, describing their main achievements and limitations. In addition, we discuss remaining open questions and potential future directions of development.
40Ar/39Ar dating of a hydrothermal pegmatitic buddingtonite–muscovite assemblage from Volyn, Ukraine
(2022)
We determined Ar-40/Ar-39 ages of buddingtonite, occurring together with muscovite, with the laser-ablation method. This is the first attempt to date the NH4-feldspar buddingtonite, which is typical for sedimentary-diagenetic environments of sediments, rich in organic matter, or in hydrothermal environments, associated with volcanic geyser systems. The sample is a hydrothermal breccia, coming from the Paleoproterozoic pegmatite field of the Korosten Plutonic Complex, Volyn, Ukraine. A detailed characterization by optical methods, electron microprobe analyses, backscattered electron imaging, and IR analyses showed that the buddingtonite consists of euhedral-appearing platy crystals of tens of micrometers wide, 100 or more micrometers in length, which consist of fine-grained fibers of <= 1 mu m thickness. The crystals are sector and growth zoned in terms of K-NH4-H3O content. The content of K allows for an age determination with the Ar-40/Ar-39 method, as well as in the accompanying muscovite, intimately intergrown with the buddingtonite. The determinations on muscovite yielded an age of 1491 +/- 9 Ma, interpreted as the hydrothermal event forming the breccia. However, buddingtonite apparent ages yielded a range of 563 +/- 14 Ma down to 383 +/- 12 Ma, which are interpreted as reset ages due to Ar loss of the fibrous buddingtonite crystals during later heating. We conclude that buddingtonite is suited for Ar-40/Ar-39 age determinations as a supplementary method, together with other methods and minerals; however, it requires a detailed mineralogical characterization, and the ages will likely represent minimum ages.
Leaching zones within potash seams generally represent a significant risk to subsurface mining operations and the construction of technical caverns in salt rocks, but their temporal and spatial formation has been investigated only rudimentarily to date. To the knowledge of the authors, current reactive transport simulation implementations are not capable to address hydraulic-chemical interactions within potash salt. For this reason, a reactive transport model has been developed and complemented by an innovative approach to calculate the interchange of minerals and solution at the water-rock interface. Using this model, a scenario analysis was carried out based on a carnallite-bearing potash seam. The results show that the evolution of leaching zones depends on the mineral composition and dissolution rate of the original salt rock, and that the formation can be classified by the dimensionless parameters of Peclet (Pe) and Damkohler (Da). For Pe > 2 and Da > 1, a funnel-shaped leaching zone is formed, otherwise the dissolution front is planar. Additionally, Da > 1 results in the formation of a sylvinitic zone and a flow barrier. Most scenarios represent hybrid forms of these cases. The simulated shapes and mineralogies are confirmed by literature data and can be used to assess the hazard potential.
Barite scalings are a common cause of permanent formation damage to deep geothermal reservoirs. Well injectivity can be impaired because the ooling of saline fluids reduces the solubility of barite, and the continuous re-injection of supersaturated fluids forces barite to precipitate in the host rock. Stimulated reservoirs in the Upper Rhine Graben often have multiple relevant flow paths in the porous matrix and fracture zones, sometimes spanning multiple stratigraphical units to achieve the economically necessary injectivity. While the influence of barite scaling on injectivity has been investigated for purely porous media, the role of fractures within reservoirs consisting of both fractured and porous sections is still not well understood. Here, we present hydro-chemical simulations of a dual-layer geothermal reservoir to study the long-term impact of barite scale formation on well injectivity. Our results show that, compared to purely porous reservoirs, fractured porous reservoirs have a significantly reduced scaling risk by up to 50%, depending on the flow rate ratio of fractures. Injectivity loss is doubled, however, if the amount of active fractures is increased by one order of magnitude, while the mean fracture aperture is decreased, provided the fractured aquifer dictates the injection rate. We conclude that fractured, and especially hydraulically stimulated, reservoirs are generally less affected by barite scaling and that large, but few, fractures are favourable. We present a scaling score for fractured-porous reservoirs, which is composed of easily derivable quantities such as the radial equilibrium length and precipitation potential. This score is suggested for use approximating the scaling potential and its impact on injectivity of a fractured-porous reservoir for geothermal exploitation.
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 spatial pattern of extreme precipitation from 40 years of gauge data in the central Himalaya
(2022)
The topography of the Himalaya exerts a substantial control on the spatial distribution of monsoonal rainfall, which is a vital water source for the regional economy and population. But the occurrence of short-lived and high-intensity precipitation results in socio-economic losses. This study relies on 40 years of daily data from 204 ground stations in Nepal to derive extreme precipitation thresholds, amounts, and days at the 95th percentile. We additionally determine the precipitation magnitude-frequency relation. We observe that extreme precipitation amounts follow an almost uniform band parallel to topographic contour lines in the southern Himalaya mountains in central and eastern Nepal but not in western Nepal. The relationship of extreme precipitation indices with topographic relief shows that extreme precipitation thresholds decrease with increasing elevation, but extreme precipitation days increase in higher elevation areas. Furthermore, stations above 1 km elevation exhibit a power-law relation in the rainfall magnitude-frequency framework. Stations at higher elevations generally have lower values of power-law exponents than low elevation areas. This suggests a fundamentally different behaviour of the rainfall distribution and an increased occurrence of extreme rainfall storms in the high elevation areas of Nepal.
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.
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.
Sulfate reduction is the quantitatively most important process to degrade organic matter in anoxic marine sediment and has been studied intensively in a variety of settings. Guaymas Basin, a young marginal ocean basin, offers the unique opportunity to study sulfate reduction in an environment characterized by organic-rich sediment, high sedimentation rates, and high geothermal gradients (100-958 degrees C km(-1)). We measured sulfate reduction rates (SRR) in samples taken during the International Ocean Discovery Program (IODP) Expedition 385 using incubation experiments with radiolabeled (SO42-)-S-35 carried out at in situ pressure and temperature. The highest SRR (387 nmol cm(-3) d(-1)) was recorded in near-surface sediments from Site U1548C, which had the steepest geothermal gradient (958 degrees C km(-1)). At this site, SRR were generally over an order of magnitude higher than at similar depths at other sites (e.g., 387-157 nmol cm(-3) d(-1) at 1.9 mbsf from Site U1548C vs. 46-1.0 nmol cm(-3) d(-1) at 2.1 mbsf from Site U1552B). Site U1546D is characterized by a sill intrusion, but it had already reached thermal equilibrium and SRR were in the same range as nearby Site U1545C, which is minimally affected by sills. The wide temperature range observed at each drill site suggests major shifts in microbial community composition with very different temperature optima but awaits confirmation by molecular biological analyses. At the transition between the mesophilic and thermophilic range around 40 degrees C-60 degrees C, sulfate-reducing activity appears to be decreased, particularly in more oligotrophic settings, but shows a slight recovery at higher temperatures.
Barite scales in geothermal installations are a highly unwanted effect of circulating deep saline fluids. They build up in the reservoir if supersaturated fluids are re-injected, leading to irreversible loss of injectivity. A model is presented for calculating the total expected barite precipitation. To determine the related injectivity decline over time, the spatial precipitation distribution in the subsurface near the injection well is assessed by modelling barite growth kinetics in a radially diverging Darcy flow domain. Flow and reservoir properties as well as fluid chemistry are chosen to represent reservoirs subject to geothermal exploration located in the North German Basin (NGB) and the Upper Rhine Graben (URG) in Germany. Fluids encountered at similar depths are hotter in the URG, while they are more saline in the NGB. The associated scaling amount normalised to flow rate is similar for both regions. The predicted injectivity decline after 10 years, on the other hand, is far greater for the NGB (64%) compared to the URG (24%), due to the temperature- and salinity-dependent precipitation rate. The systems in the NGB are at higher risk. Finally, a lightweight score is developed for approximating the injectivity loss using the Damkohler number, flow rate and total barite scaling potential. This formula can be easily applied to geothermal installations without running complex reactive transport simulations.
This study, based on new and high quality in situ observations, quantifies for the first time, the individual contributions of light-absorbing aerosols (black carbon (BC), brown carbon (BrC) and dust) to aerosol absorption over the Indo-Gangetic Plain (IGP) and the Himalayan foothill region, a relatively poorly studied region with several sensitive ecosystems of global importance, as well as highly vulnerable populations. The annual and seasonal average single scattering albedo (SSA) over Kathmandu is the lowest of all the locations. The SSA over Kathmandu is < 0.89 during all seasons, which confirms the dominance of light-absorbing carbonaceous aerosols from local and regional sources over Kathmandu. It is observed here that the SSA decreases with increasing elevation, confirming the dominance of light absorbing carbonaceous aerosols at higher elevations. In contrast, the SSA over the IGP does not exhibit a pronounced spatial variation. BC dominates (>= 75%) the aerosol absorption over the IGP and the Himalayan foothills throughout the year. Higher BC concentration at elevated locations in the Himalayas leads to lower SSA at elevated locations in the Himalayas. The contribution of dust to aerosol absorption is higher throughout the year over the IGP than over the Himalayan foothills. The aerosol absorption over South Asia is very high, exceeding available observations over East Asia, and also exceeds previous model estimates. This quantification will be valuable as observational constraints to help improve regional simulations of climate change, impacts on the glaciers and the hydrological cycle, and will help to direct the focus towards BC as the main contributor to aerosol-induced warming in the region.
The SiDroForest (Siberian drone-mapped forest inventory) data collection is an attempt to remedy the scarcity of forest structure data in the circumboreal region by providing adjusted and labeled tree-level and vegetation plot-level data for machine learning and upscaling purposes. We present datasets of vegetation composition and tree and plot level forest structure for two important vegetation transition zones in Siberia, Russia; the summergreen-evergreen transition zone in Central Yakutia and the tundra-taiga transition zone in Chukotka (NE Siberia). The SiDroForest data collection consists of four datasets that contain different complementary data types that together support in-depth analyses from different perspectives of Siberian Forest plot data for multi-purpose applications. i. Dataset 1 provides unmanned aerial vehicle (UAV)-borne data products covering the vegetation plots surveyed during fieldwork (Kruse et al., 2021, ). The dataset includes structure-from-motion (SfM) point clouds and red-green-blue (RGB) and red-green-near-infrared (RGN) orthomosaics. From the orthomosaics, point-cloud products were created such as the digital elevation model (DEM), canopy height model (CHM), digital surface model (DSM) and the digital terrain model (DTM). The point-cloud products provide information on the three-dimensional (3D) structure of the forest at each plot. Dataset 2 contains spatial data in the form of point and polygon shapefiles of 872 individually labeled trees and shrubs that were recorded during fieldwork at the same vegetation plots (van Geffen et al., 2021c, ). The dataset contains information on tree height, crown diameter, and species type. These tree and shrub individually labeled point and polygon shapefiles were generated on top of the RGB UVA orthoimages. The individual tree information collected during the expedition such as tree height, crown diameter, and vitality are provided in table format. This dataset can be used to link individual information on trees to the location of the specific tree in the SfM point clouds, providing for example, opportunity to validate the extracted tree height from the first dataset. The dataset provides unique insights into the current state of individual trees and shrubs and allows for monitoring the effects of climate change on these individuals in the future. Dataset 3 contains a synthesis of 10 000 generated images and masks that have the tree crowns of two species of larch ( and ) automatically extracted from the RGB UAV images in the common objects in context (COCO) format (van Geffen et al., 2021a, ). As machine-learning algorithms need a large dataset to train on, the synthetic dataset was specifically created to be used for machine-learning algorithms to detect Siberian larch species. Larix gmeliniiLarix cajanderiDataset 4 contains Sentinel-2 (S-2) Level-2 bottom-of-atmosphere processed labeled image patches with seasonal information and annotated vegetation categories covering the vegetation plots (van Geffen et al., 2021b, ). The dataset is created with the aim of providing a small ready-to-use validation and training dataset to be used in various vegetation-related machine-learning tasks. It enhances the data collection as it allows classification of a larger area with the provided vegetation classes. The SiDroForest data collection serves a variety of user communities. <br /> The detailed vegetation cover and structure information in the first two datasets are of use for ecological applications, on one hand for summergreen and evergreen needle-leaf forests and also for tundra-taiga ecotones. Datasets 1 and 2 further support the generation and validation of land cover remote-sensing products in radar and optical remote sensing. In addition to providing information on forest structure and vegetation composition of the vegetation plots, the third and fourth datasets are prepared as training and validation data for machine-learning purposes. For example, the synthetic tree-crown dataset is generated from the raw UAV images and optimized to be used in neural networks. Furthermore, the fourth SiDroForest dataset contains S-2 labeled image patches processed to a high standard that provide training data on vegetation class categories for machine-learning classification with JavaScript Object Notation (JSON) labels provided. The SiDroForest data collection adds unique insights into remote hard-to-reach circumboreal forest regions.
Mineral resource exploration and mining is an essential part of today's high-tech industry. Elements such as rare-earth elements (REEs) and copper are, therefore, in high demand. Modern exploration techniques from multiple platforms (e.g., spaceborne and airborne), to detect and map the spectral characteristics of the materials of interest, require spectral libraries as an essential reference. They include field and laboratory spectral information in combination with geochemical analyses for validation. Here, we present a collection of REE- and copper-related hyperspectral spectra with associated geochemical information. The libraries contain reflectance spectra from rare-earth element oxides, REE-bearing minerals, copper-bearing minerals and mine surface samples from the Apliki copper-gold-pyrite mine in the Republic of Cyprus. The samples were measured with the HySpex imaging spectrometers in the visible and near infrared (VNIR) and shortwave infrared (SWIR) range (400-2500 nm). The geochemical validation of each sample is provided with the reflectance spectra. The spectral libraries are openly available to assist future mineral mapping campaigns and laboratory spectroscopic analyses. The spectral libraries and corresponding geochemistry are published via GFZ Data Services with the following DOIs: https://doi.org/10.5880/GFZ.1.4.2019.004 (13 REE-bearing minerals and 16 oxide powders, Koerting et al., 2019a), https://doi.org/10.5880/GFZ.1.4.2019.003 (20 copper-bearing minerals, Koellner et al., 2019), and https://doi.org/10.5880/GFZ.1.4.2019.005 (37 copper-bearing surface material samples from the Apliki coppergold-pyrite mine in Cyprus, Koerting et al., 2019b). All spectral libraries are united and comparable by the internally consistent method of hyperspectral data acquisition in the laboratory.
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).
Nature-based solutions (NBS) are seen as a promising adaptation measure that sustainably deals with diverse societal challenges, while simultaneously delivering multiple benefits. Nature-based solutions have been highlighted as a resilient and sustainable means of mitigating floods and other hazards globally. This study examined diverging conceptualizations of NBS, as well as the attitudinal (for example, emotions and beliefs) and contextual (for example, legal and political aspects) barriers and drivers of NBS for flood risks in South Korea. Semistructured interviews were conducted with 11 experts and focused on the topic of flood risk measures and NBS case studies. The analysis found 11 barriers and five drivers in the attitudinal domain, and 13 barriers and two drivers in the contextual domain. Most experts see direct monetary benefits as an important attitudinal factor for the public. Meanwhile, the cost-effectiveness of NBS and their capacity to cope with flood risks were deemed influential factors that could lead decision makers to opt for NBS. Among the contextual factors, insufficient systems to integrate NBS in practice and the ideologicalization of NBS policy were found to be peculiar barriers, which hinder consistent realization of initiatives and a long-term national plan for NBS. Understanding the barriers and drivers related to the mainstreaming of NBS is critical if we are to make the most of such solutions for society and nature. It is also essential that we have a shared definition, expectation, and vision of NBS.
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.
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>
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.
Here I present a comparison between two of the most widely used reduced-complexity models for the representation of sediment transport and deposition processes, namely the transport-limited (or TL) model and the under-capacity (or xi-q) model more recently developed by Davy and Lague (2009). Using both models, I investigate the behavior of a sedimentary continental system of length L fed by a fixed sedimentary flux from a catchment of size A(0) in a nearby active orogen through which sediments transit to a fixed base level representing a large river, a lake or an ocean. This comparison shows that the two models share the same steady-state solution, for which I derive a simple 1D analytical expression that reproduces the major features of such sedimentary systems: a steep fan that connects to a shallower alluvial plain. The resulting fan geometry obeys basic observational constraints on fan size and slope with respect to the upstream drainage area, A(0). The solution is strongly dependent on the size of the system, L, in comparison to a distance L-0, which is determined by the size of A(0), and gives rise to two fundamentally different types of sedimentary systems: a constrained system where L < L-0 and open systems where L > L-0. I derive simple expressions that show the dependence of the system response time on the system characteristics, such as its length, the size of the upstream catchment area, the amplitude of the incoming sedimentary flux and the respective rate parameters (diffusivity or erodibility) for each of the two models. I show that the xi-q model predicts longer response times. I demonstrate that although the manner in which signals propagates through the sedimentary system differs greatly between the two models, they both predict that perturbations that last longer than the response time of the system can be recorded in the stratigraphy of the sedimentary system and in particular of the fan. Interestingly, the xi-q model predicts that all perturbations in the incoming sedimentary flux will be transmitted through the system, whereas the TL model predicts that rapid perturbations cannot. I finally discuss why and under which conditions these differences are important and propose observational ways to determine which of the two models is most appropriate to represent natural systems.
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.
Draft Genome Sequence of Nocardioides alcanivorans NGK65(T), a Hexadecane-Degrading Bacterium
(2022)
The Gram-positive bacterium Nocardioides alcanivorans NGK65(T) was isolated from plastic-polluted soil and cultivated on medium with polyethylene as the single carbon source. Nanopore sequencing revealed the presence of candidate enzymes for the biodegradation of polyethylene. Here, we report the draft genome of this newly described member of the terrestrial plastisphere.
During the period 750-600 Ma ago, prior to the final break-up of the supercontinent Rodinia, the crust of both the North American Craton and Baltica was intruded by significant amounts of rift-related magmas originating from the mantle. In the Proterozoic crust of Southern Norway, the 580 Ma old Fen carbonatite-ultramafic complex is a representative of this type of rocks. In this paper, we report the occurrence of an ultramafic lamprophyre dyke which possibly is linked to the Fen complex, although Ar-40/Ar-39 data from phenocrystic phlogopite from the dyke gave an age of 686 +/- 9 Ma. The lamprophyre dyke was recently discovered in one of the Kongsberg silver mines at Vinoren, Norway. Whole rock geochemistry, geochronological and mineralogical data from the ultramafic lamprophyre dyke are presented aiming to elucidate its origin and possible geodynamic setting. From the whole-rock composition of the Vinoren dyke, the rock could be recognized as transitional between carbonatite and kimberlite-II (orangeite). From its diagnostic mineralogy, the rock is classified as aillikite. The compositions and xenocrystic nature of several of the major and accessory minerals from the Vinoren aillikite are characteristic for diamondiferous rocks (kimberlites/lamproites/UML): Phlogopite with kinoshitalite-rich rims, chromite-spinel-ulvospinel series, Mg- and Mn-rich ilmenites, rutile and lucasite-(Ce). We suggest that the aillikite melt formed during partial melting of a MARID (mica-amphibole-rutile-ilmenite-diopside)-like source under CO2 fluxing. The pre-rifting geodynamic setting of the Vinoren aillikite before the Rodinia supercontinent breakup suggests a relatively thick SCLM (Subcontinental Lithospheric Mantle) during this stage and might indicate a diamond-bearing source for the parental melt. This is in contrast to the about 100 Ma younger Fen complex, which were derived from a thin SCLM.
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.
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.
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.
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.
Compound natural hazards likeEl Ninoevents cause high damage to society, which to manage requires reliable risk assessments. Damage modelling is a prerequisite for quantitative risk estimations, yet many procedures still rely on expert knowledge, and empirical studies investigating damage from compound natural hazards hardly exist. A nationwide building survey in Peru after theEl Ninoevent 2017 - which caused intense rainfall, ponding water, flash floods and landslides - enables us to apply data-mining methods for statistical groundwork, using explanatory features generated from remote sensing products and open data. We separate regions of different dominant characteristics through unsupervised clustering, and investigate feature importance rankings for classifying damage via supervised machine learning. Besides the expected effect of precipitation, the classification algorithms select the topographic wetness index as most important feature, especially in low elevation areas. The slope length and steepness factor ranks high for mountains and canyons. Partial dependence plots further hint at amplified vulnerability in rural areas. An example of an empirical damage probability map, developed with a random forest model, is provided to demonstrate the technical feasibility.
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/).
Multi-component (MC) diffusion simulations enable a process based and more precise approach to calculate transport and sorption compared to the commonly used single-component (SC) models following Fick's law. The MC approach takes into account the interaction of chemical species in the porewater with the diffuse double layer (DDL) adhering clay mineral surfaces. We studied the shaly, sandy and carbonate-rich facies of the Opalinus Clay. High clay contents dominate diffusion and sorption of uranium. The MC simulations show shorter diffusion lengths than the SC models due to anion exclusion from the DDL. This hampers diffusion of the predominant species CaUO2(CO3)32-. On the one side, species concentrations and ionic strengths of the porewater and on the other side surface charge of the clay minerals control the composition and behaviour of the DDL. For some instances, it amplifies the diffusion of uranium. We developed a workflow to transfer computationally intensive MC simulations to SC models via calibrated effective diffusion and distribution coefficients. Simulations for one million years depict maximum uranium diffusion lengths between 10 m and 35 m. With respect to the minimum requirement of a thickness of 100 m, the Opalinus Clay seems to be a suitable host rock for nuclear waste repositories.
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.
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.
Increasing arctic coastal erosion rates imply a greater release of sediments and organic matter into the coastal zone. With 213 sediment samples taken around Herschel Island-Qikiqtaruk, Canadian Beaufort Sea, we aimed to gain new insights on sediment dynamics and geochemical properties of a shallow arctic nearshore zone. Spatial characteristics of nearshore sediment texture (moderately to poorly sorted silt) are dictated by hydrodynamic processes, but ice-related processes also play a role. We determined organic matter (OM) distribution and inferred the origin and quality of organic carbon by C/N ratios and stable carbon isotopes delta C-13. The carbon content was higher offshore and in sheltered areas (mean: 1.0 wt.%., S.D.: 0.9) and the C/N ratios also showed a similar spatial pattern (mean: 11.1, S.D.: 3.1), while the delta C-13 (mean: -26.4 parts per thousand VPDB, S.D.: 0.4) distribution was more complex. We compared the geochemical parameters of our study with terrestrial and marine samples from other studies using a bootstrap approach. Sediments of the current study contained 6.5 times and 1.8 times less total organic carbon than undisturbed and disturbed terrestrial sediments, respectively. Therefore, degradation of OM and separation of carbon pools take place on land and continue in the nearshore zone, where OM is leached, mineralized, or transported beyond the study area.
Diagenetic trends of synthetic reservoir sandstone properties assessed by digital rock physics
(2021)
Quantifying interactions and dependencies among geometric, hydraulic and mechanical properties of reservoir sandstones is of particular importance for the exploration and utilisation of the geological subsurface and can be assessed by synthetic sandstones comprising the microstructural complexity of natural rocks. In the present study, three highly resolved samples of the Fontainebleau, Berea and Bentheim sandstones are generated by means of a process-based approach, which combines the gravity-driven deposition of irregularly shaped grains and their diagenetic cementation by three different schemes. The resulting evolution in porosity, permeability and rock stiffness is examined and compared to the respective micro-computer tomographic (micro-CT) scans. The grain contact-preferential scheme implies a progressive clogging of small throats and consequently produces considerably less connected and stiffer samples than the two other schemes. By contrast, uniform quartz overgrowth continuously alters the pore space and leads to the lowest elastic properties. The proposed stress-dependent cementation scheme combines both approaches of contact-cement and quartz overgrowth, resulting in granulometric, hydraulic and elastic properties equivalent to those of the respective micro-CT scans, where bulk moduli slightly deviate by 0.8%, 4.9% and 2.5% for the Fontainebleau, Berea and Bentheim sandstone, respectively. The synthetic samples can be further altered to examine the impact of mineral dissolution or precipitation as well as fracturing on various petrophysical correlations, which is of particular relevance for numerous aspects of a sustainable subsurface utilisation.
ArcticBeach v1.0
(2022)
In the Arctic, air temperatures are increasing and sea ice is declining, resulting in larger waves and a longer open water season, all of which intensify the thaw and erosion of ice-rich coasts. Climate change has been shown to increase the rate of Arctic coastal erosion, causing problems for Arctic cultural heritage, existing industrial, military, and civil infrastructure, as well as changes in nearshore biogeochemistry. Numerical models that reproduce historical and project future Arctic erosion rates are necessary to understand how further climate change will affect these problems, and no such model yet exists to simulate the physics of erosion on a pan-Arctic scale. We have coupled a bathystrophic storm surge model to a simplified physical erosion model of a permafrost coastline. This Arctic erosion model, called ArcticBeach v1.0, is a first step toward a physical parameterization of Arctic shoreline erosion for larger-scale models. It is forced by wind speed and direction, wave period and height, sea surface temperature, all of which are masked during times of sea ice cover near the coastline. Model tuning requires observed historical retreat rates (at least one value), as well as rough nearshore bathymetry. These parameters are already available on a pan-Arctic scale. The model is validated at three study sites at 1) Drew Point (DP), Alaska, 2) Mamontovy Khayata (MK), Siberia, and 3) Veslebogen Cliffs, Svalbard. Simulated cumulative retreat rates for DP and MK respectively (169 and 170 m) over the time periods studied at each site (2007-2016, and 1995-2018) are found to the same order of magnitude as observed cumulative retreat (172 and 120 m). The rocky Veslebogen cliffs have small observed cumulative retreat rates (0.05 m over 2014-2016), and our model was also able to reproduce this same order of magnitude of retreat (0.08 m). Given the large differences in geomorphology between the study sites, this study provides a proof-of-concept that ArcticBeach v1.0 can be applied on very different permafrost coastlines. ArcticBeach v1.0 provides a promising starting point to project retreat of Arctic shorelines, or to evaluate historical retreat in places that have had few observations.
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 computational costs associated with coupled reactive transport simulations are mostly due to the chemical subsystem: replacing it with a pre-trained statistical surrogate is a promising strategy to achieve decisive speedups at the price of small accuracy losses and thus to extend the scale of problems which can be handled. We introduce a hierarchical coupling scheme in which "full-physics" equation-based geochemical simulations are partially replaced by surrogates. Errors in mass balance resulting from multivariate surrogate predictions effectively assess the accuracy of multivariate regressions at runtime: inaccurate surrogate predictions are rejected and the more expensive equation-based simulations are run instead. Gradient boosting regressors such as XGBoost, not requiring data standardization and being able to handle Tweedie distributions, proved to be a suitable emulator. Finally, we devise a surrogate approach based on geochemical knowledge, which overcomes the issue of robustness when encountering previously unseen data and which can serve as a basis for further development of hybrid physics-AI modelling.
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.
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.
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.
Flood risk management in Germany follows an integrative approach in which both private households and businesses can make an important contribution to reducing flood damage by implementing property-level adaptation measures. While the flood adaptation behavior of private households has already been widely researched, comparatively less attention has been paid to the adaptation strategies of businesses. However, their ability to cope with flood risk plays an important role in the social and economic development of a flood-prone region. Therefore, using quantitative survey data, this study aims to identify different strategies and adaptation drivers of 557 businesses damaged by a riverine flood in 2013 and 104 businesses damaged by pluvial or flash floods between 2014 and 2017. Our results indicate that a low perceived self-efficacy may be an important factor that can reduce the motivation of businesses to adapt to flood risk. Furthermore, property-owners tended to act more proactively than tenants. In addition, high experience with previous flood events and low perceived response costs could strengthen proactive adaptation behavior. These findings should be considered in business-tailored risk communication.
The Fram Strait is an area with a relatively low and irregular distribution of diatom microfossils in surface sediments, and thus microfossil records are scarce, rarely exceed the Holocene, and contain sparse information about past richness and taxonomic composition. These attributes make the Fram Strait an ideal study site to test the utility of sedimentary ancient DNA (sedaDNA) metabarcoding. Amplifying a short, partial rbcL marker from samples of sediment core MSM05/5-712-2 resulted in 95.7% of our sequences being assigned to diatoms across 18 different families, with 38.6% of them being resolved to species and 25.8% to genus level. Independent replicates show a high similarity of PCR products, especially in the oldest samples. Diatom sedaDNA richness is highest in the Late Weichselian and lowest in Mid- and Late Holocene samples. Taxonomic composition is dominated by cold-water and sea-ice-associated diatoms and suggests several reorganisations - after the Last Glacial Maximum, after the Younger Dryas, and after the Early and after the Mid-Holocene. Different sequences assigned to, amongst others, Chaetoceros socialis indicate the detectability of intra-specific diversity using sedaDNA. We detect no clear pattern between our diatom sedaDNA record and the previously published IP25 record of this core, although proportions of pennate diatoms increase with higher IP25 concentrations and proportions of Nitzschia cf. frigida exceeding 2% of the assemblage point towards past sea-ice presence.
The Victoria microplate between the Eastern and Western Branches of the East African Rift System is one of the largest continental microplates on Earth. In striking contrast to its neighboring plates, Victoria rotates counterclockwise with respect to Nubia. The underlying cause of this distinctive rotation has remained elusive so far. Using 3D numerical models, we investigate the role of pre-existing lithospheric heterogeneities in continental microplate rotation. We find that Victoria's rotation is primarily controlled by the distribution of rheologically stronger zones that transmit the drag of the major plates to the microplate and of the mechanically weaker mobile belts surrounding Victoria that facilitate rotation. Our models reproduce Victoria's GPS-derived counterclockwise rotation as well as key complexities of the regional tectonic stress field. These results reconcile competing ideas on the opening of the rift system by highlighting differences in orientation of the far-field divergence, local extension, and the minimum horizontal stress. One of the largest continental microplates on Earth is situated in the center of the East African Rift System, and oddly, the Victoria microplate rotates counterclockwise with respect to the neighboring African tectonic plate. Here, the authors' modelling results suggest that Victoria microplate rotation is caused by edge-driven lithospheric processes related to the specific geometry of rheologically weak and strong regions.
(40)A/Ar-39 step-heating of mica and amphibole megacrysts from hauyne-bearing olivine melilitite scoria/tephra from the Zelezna hurka yielded a 435 +/- 108 ka isotope correlation age for phlogopite and a more imprecise 1.55 Ma total gas age of the kaersutite megacryst. The amphibole megacrysts may constitute the first, and the younger phlogopite megacrysts the later phase of mafic, hydrous melilitic magma crystallization. It cannot be ruled out that the amphibole megacrysts are petrogenetically unrelated to tephra and phlogopite megacrysts and were derived from mantle xenoliths or disaggregated older, deep crustal pegmatites. This is in line both with the rarity of amphibole at Zelezna hurka and with the observed signs of magmatic resorption at the edges of amphibole crystals.
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.
Drainage divide networks
(2020)
Drainage divides are organized into tree-like networks that may record information about drainage divide mobility. However, views diverge about how to best assess divide mobility. Here, we apply a new approach of automatically extracting and ordering drainage divide networks from digital elevation models to results from landscape evolution model experiments. We compared landscapes perturbed by strike-slip faulting and spatiotemporal variations in erodibility to a reference model to assess which topographic metrics (hillslope relief, flow distance, and chi) are diagnostic of divide mobility. Results show that divide segments that are a minimum distance of similar to 5 km from river confluences strive to attain constant values of hillslope relief and flow distance to the nearest stream. Disruptions of such patterns can be related to mobile divides that are lower than stable divides, closer to streams, and often asymmetric in shape. In general, we observe that drainage divides high up in the network, i.e., at great distances from river confluences, are more susceptible to disruptions than divides closer to these confluences and are thus more likely to record disturbance for a longer time period. We found that across-divide differences in hillslope relief proved more useful for assessing divide migration than other tested metrics. However, even stable drainage divide networks exhibit across-divide differences in any of the studied topographic metrics. Finally, we propose a new metric to quantify the connectivity of divide junctions.
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.
In this study, we analyze interactions in lake and lake catchment systems of a continuous permafrost area. We assessed colored dissolved organic matter (CDOM) absorption at 440 nm (a(440)(CDOM)) and absorption slope (S300-500) in lakes using field sampling and optical remote sensing data for an area of 350 km(2) in Central Yamal, Siberia. Applying a CDOM algorithm (ratio of green and red band reflectance) for two high spatial resolution multispectral GeoEye-1 and Worldview-2 satellite images, we were able to extrapolate the a()(CDOM) data from 18 lakes sampled in the field to 356 lakes in the study area (model R-2 = 0.79). Values of a(440)(CDOM) in 356 lakes varied from 0.48 to 8.35 m(-1) with a median of 1.43 m(-1). This a()(CDOM) dataset was used to relate lake CDOM to 17 lake and lake catchment parameters derived from optical and radar remote sensing data and from digital elevation model analysis in order to establish the parameters controlling CDOM in lakes on the Yamal Peninsula. Regression tree model and boosted regression tree analysis showed that the activity of cryogenic processes (thermocirques) in the lake shores and lake water level were the two most important controls, explaining 48.4% and 28.4% of lake CDOM, respectively (R-2 = 0.61). Activation of thermocirques led to a large input of terrestrial organic matter and sediments from catchments and thawed permafrost to lakes (n = 15, mean a(440)(CDOM) = 5.3 m(-1)). Large lakes on the floodplain with a connection to Mordy-Yakha River received more CDOM (n = 7, mean a(440)(CDOM) = 3.8 m(-1)) compared to lakes located on higher terraces.
Glacial-cycle simulations of the Antarctic Ice Sheet with the Parallel Ice Sheet Model (PISM)
(2020)
Simulations of the glacial-interglacial history of the Antarctic Ice Sheet provide insights into dynamic threshold behavior and estimates of the ice sheet's contributions to global sea-level changes for the past, present and future. However, boundary conditions are weakly constrained, in particular at the interface of the ice sheet and the bedrock. Also climatic forcing covering the last glacial cycles is uncertain, as it is based on sparse proxy data. <br /> We use the Parallel Ice Sheet Model (PISM) to investigate the dynamic effects of different choices of input data, e.g., for modern basal heat flux or reconstructions of past changes of sea level and surface temperature. As computational resources are limited, glacial-cycle simulations are performed using a comparably coarse model grid of 16 km and various parameterizations, e.g., for basal sliding, iceberg calving, or for past variations in precipitation and ocean temperatures. In this study we evaluate the model's transient sensitivity to corresponding parameter choices and to different boundary conditions over the last two glacial cycles and provide estimates of involved uncertainties. We also discuss isolated and combined effects of climate and sea-level forcing. Hence, this study serves as a "cookbook" for the growing community of PISM users and paleo-ice sheet modelers in general. <br /> For each of the different model uncertainties with regard to climatic forcing, ice and Earth dynamics, and basal processes, we select one representative model parameter that captures relevant uncertainties and motivates corresponding parameter ranges that bound the observed ice volume at present. The four selected parameters are systematically varied in a parameter ensemble analysis, which is described in a companion paper.
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.
Fractures efficiently affect fluid flow in geological formations, and thereby determine mass and energy transport in reservoirs, which are not least exploited for economic resources. In this context, their response to mechanical and thermal changes, as well as fluid-rock interactions, is of paramount importance. In this study, a two-stage flow-through experiment was conducted on a pure quartz sandstone core of low matrix permeability, containing one single macroscopic tensile fracture. In the first short-term stage, the effects of mechanical and hydraulic aperture on pressure and temperature cycles were investigated. The purpose of the subsequent intermittent-flow long-term (140 days) stage was to constrain the evolution of the geometrical and hydraulic fracture properties resulting from pressure solution. Deionized water was used as the pore fluid, and permeability, as well as the effluent Si concentrations, were systematically measured. Overall, hydraulic aperture was shown to be significantly less affected by pressure, temperature and time, in comparison to mechanical aperture. During the long-term part of the experiment at 140 degrees C, the effluent Si concentrations likely reached a chemical equilibrium state within less than 8 days of stagnant flow, and exceeded the corresponding hydrostatic quartz solubility at this temperature. This implies that the pressure solution was active at the contacting fracture asperities, both at 140 degrees C and after cooling to 33 degrees C. The higher temperature yielded a higher dissolution rate and, consequently, a faster attainment of chemical equilibrium within the contact fluid. X-ray mu CT observations evidenced a noticeable increase in fracture contact area ratio, which, in combination with theoretical considerations, implies a significant decrease in mechanical aperture. In contrast, the sample permeability, and thus the hydraulic fracture aperture, virtually did not vary. In conclusion, pressure solution-induced fracture aperture changes are affected by the degree of time-dependent variations in pore fluid composition. In contrast to the present case of a quasi-closed system with mostly stagnant flow, in an open system with continuous once-through fluid flow, the activity of the pressure solution may be amplified due to the persistent fluid-chemical nonequilibrium state, thus possibly enhancing aperture and fracture permeability changes.
Relative pollen productivity (RPP) estimates are fractionate values, often in relation to Poaceae, that allow vegetation cover to be estimated from pollen counts with the help of models. RPP estimates are especially used in the scientific community in Europe and China, with a few studies in North America. Here we present a comprehensive compilation of available northern hemispheric RPP studies and their results arising from 51 publications with 60 sites and 131 taxa. This compilation allows scientists to identify data gaps in need of further RPP analyses but can also aid them in finding an RPP set for their study region. We also present a taxonomically harmonised, unified RPP dataset for the Northern Hemisphere and subsets for North America (including Greenland), Europe (including arctic Russia), and China, which we generated from the available studies. The unified dataset gives the mean RPP for 55 harmonised taxa as well as fall speeds, which are necessary to reconstruct vegetation cover from pollen counts and RPP values. Data are openly available at https://doi.org/10.1594/PANGAEA.922661 (Wieczorek and Herzschuh, 2020).
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.
Diet analysis of bats killed at wind turbines suggests large-scale losses of trophic interactions
(2022)
Agricultural practice has led to landscape simplification and biodiversity decline, yet recently, energy-producing infrastructures, such as wind turbines, have been added to these simplified agroecosystems, turning them into multi-functional energy-agroecosystems. Here, we studied the trophic interactions of bats killed at wind turbines using a DNA metabarcoding approach to shed light on how turbine-related bat fatalities may possibly affect local habitats. Specifically, we identified insect DNA in the stomachs of common noctule bats (Nyctalus noctula) killed by wind turbines in Germany to infer in which habitats these bats hunted. Common noctule bats consumed a wide variety of insects from different habitats, ranging from aquatic to terrestrial ecosystems (e.g., wetlands, farmland, forests, and grasslands). Agricultural and silvicultural pest insects made up about 20% of insect species consumed by the studied bats. Our study suggests that the potential damage of wind energy production goes beyond the loss of bats and the decline of bat populations. Bat fatalities at wind turbines may lead to the loss of trophic interactions and ecosystem services provided by bats, which may add to the functional simplification and impaired crop production, respectively, in multi-functional ecosystems.
Underground coal gasification (UCG) is an in situ conversion technique that enables the production of high-calorific synthesis gas from resources that are economically not minable by conventional methods. A broad range of end-use options is available for the synthesis gas, including fuels and chemical feedstock production. Furthermore, UCG also offers a high potential for integration with Carbon Capture and Storage (CCS) to mitigate greenhouse gas emissions. In the present study, a stoichiometric equilibrium model, based on minimization of the Gibbs function has been used to estimate the equilibrium composition of the synthesis gas. Thereto, we further developed and applied a proven thermodynamic equilibrium model to simulate the relevant thermochemical coal conversion processes (pyrolysis and gasification). Our modeling approach has been validated against thermodynamic models, laboratory gasification experiments and UCG field trial data reported in the literature. The synthesis gas compositions have been found to be in good agreement under a wide range of different operating conditions. Consequently, the presented modeling approach enables an efficient quantification of synthesis gas quality resulting from UCG, considering varying coal and oxidizer compositions at deposit-specific pressures and temperatures.
Many widely used observational data sets are comprised of several overlapping instrument records. While data inter-calibration techniques often yield continuous and reliable data for trend analysis, less attention is generally paid to maintaining higher-order statistics such as variance and autocorrelation. A growing body of work uses these metrics to quantify the stability or resilience of a system under study and potentially to anticipate an approaching critical transition in the system. Exploring the degree to which changes in resilience indicators such as the variance or autocorrelation can be attributed to non-stationary characteristics of the measurement process – rather than actual changes in the dynamical properties of the system – is important in this context. In this work we use both synthetic and empirical data to explore how changes in the noise structure of a data set are propagated into the commonly used resilience metrics lag-one autocorrelation and variance. We focus on examples from remotely sensed vegetation indicators such as vegetation optical depth and the normalized difference vegetation index from different satellite sources. We find that time series resulting from mixing signals from sensors with varied uncertainties and covering overlapping time spans can lead to biases in inferred resilience changes. These biases are typically more pronounced when resilience metrics are aggregated (for example, by land-cover type or region), whereas estimates for individual time series remain reliable at reasonable sensor signal-to-noise ratios. Our work provides guidelines for the treatment and aggregation of multi-instrument data in studies of critical transitions and resilience.
Records from ocean bottom seismometers (OBSs) are highly contaminated by noise, which is much stronger compared to data from most land stations, especially on the horizontal components. As a consequence, the high energy of the oceanic noise at frequencies below 1 Hz considerably complicates the analysis of the teleseismic earthquake signals recorded by OBSs.
Previous studies suggested different approaches to remove low-frequency noises from OBS recordings but mainly focused on the vertical component. The records of horizontal components, which are crucial for the application of many methods in passive seismological analysis of body and surface waves, could not be much improved in the teleseismic frequency band. Here we introduce a noise reduction method, which is derived from the harmonic–percussive separation algorithms used in Zali et al. (2021), in order to separate long-lasting narrowband signals from broadband transients in the OBS signal. This leads to significant noise reduction of OBS records on both the vertical and horizontal components and increases the earthquake signal-to-noise ratio (SNR) without distortion of the broadband earthquake waveforms. This is demonstrated through tests with synthetic data. Both SNR and cross-correlation coefficients showed significant improvements for different realistic noise realizations. The application of denoised signals in surface wave analysis and receiver functions is discussed through tests with synthetic and real data.
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.
The Permo-Triassic period marks the time interval between Hercynian (Variscan) orogenic events in the Tien Shan and the North Pamir, and the Cimmerian accretion of the Gondwana-derived Central and South Pamir to the southern margin of the Paleo-Asian continent. A well-preserved Permo-Triassic volcano-sedimentary sequence from the Chinese North Pamir yields important information on the geodynamic evolution of Asia’s pre-Cimmerian southern margin. The oldest volcanic rocks from that section are dated to the late Guadalupian epoch by a rhyolite and a dacitic dike that gave zircon U-Pb ages of ~260 Ma. Permian volcanism was largely pyroclastic and mafic to intermediate. Upsection, a massive ignimbritic crystal tuff in the Chinese Qimgan valley was dated to 244.1 +/- 1.1 Ma, a similar unit in the nearby Gez valley to 245 +/- 11 Ma, and an associated rhyolite to 233.4 +/- 1.1 Ma. Deposition of the locally ~200 m thick crystal tuff unit follows an unconformity and marks the onset of intense, mainly mafic to intermediate, calc-alkaline magmatic activity. Triassic volcanic activity in the North Pamir was coeval with the major phase of Cimmerian intrusive activity in the Karakul-Mazar arc-accretionary complex to the south, caused by northward subduction of the Paleo-Tethys. It also coincided with the emplacement of basanitic and carbonatitic dikes and a thermal event in the South Tien Shan, to the north of our study area. Evidence for arc-related magmatic activity in a back-arc position provides strong arguments for back-arc extension or transtension and basin formation. This puts the Qimgan succession in line with a more than 1000 km long realm of extensional Triassic back-arc basins known from the North Pamir in the Kyrgyz Altyn Darya valley (Myntekin formation), the North Pamir of Tajikistan and Afghanistan, and the Afghan Hindukush (Doab formation) and further west from the Paropamisus and Kopet Dag (Aghdarband, NE Iran).
The first step towards assessing hazards in seismically active regions involves mapping capable faults and estimating their recurrence times. While the mapping of active faults is commonly based on distinct geologic and geomorphic features evident at the surface, mapping blind seismogenic faults is complicated by the absence of on-fault diagnostic features. Here we investigated the Pichilemu Fault in coastal Chile, unknown until it generated a Mw 7.0 earthquake in 2010. The lack of evident surface faulting suggests activity along a partly-hidden blind fault. We used off-fault deformed marine terraces to estimate a fault-slip rate of 0.52 ± 0.04 m/ka, which, when integrated with satellite geodesy suggests a 2.12 ± 0.2 ka recurrence time for Mw~7.0 normal-faulting earthquakes. We propose that extension in the Pichilemu region is associated with stress changes during megathrust earthquakes and accommodated by sporadic slip during upper-plate earthquakes, which has implications for assessing the seismic potential of cryptic faults along convergent margins and elsewhere.
Many institutions struggle to tap into the potential of their large archives of radar reflectivity: these data are often affected by miscalibration, yet the bias is typically unknown and temporally volatile. Still, relative calibration techniques can be used to correct the measurements a posteriori. For that purpose, the usage of spaceborne reflectivity observations from the Tropical Rainfall Measuring Mission (TRMM) and Global Precipitation Measurement (GPM) platforms has become increasingly popular: the calibration bias of a ground radar (GR) is estimated from its average reflectivity difference to the spaceborne radar (SR). Recently, Crisologo et al. (2018) introduced a formal procedure to enhance the reliability of such estimates: each match between SR and GR observations is assigned a quality index, and the calibration bias is inferred as a quality-weighted average of the differences between SR and GR. The relevance of quality was exemplified for the Subic S-band radar in the Philippines, which is greatly affected by partial beam blockage. The present study extends the concept of quality-weighted averaging by accounting for path-integrated attenuation (PIA) in addition to beam blockage. This extension becomes vital for radars that operate at the C or X band. Correspondingly, the study setup includes a C-band radar that substantially overlaps with the S-band radar. Based on the extended quality-weighting approach, we retrieve, for each of the two ground radars, a time series of calibration bias estimates from suitable SR overpasses. As a result of applying these estimates to correct the ground radar observations, the consistency between the ground radars in the region of overlap increased substantially. Furthermore, we investigated if the bias estimates can be interpolated in time, so that ground radar observations can be corrected even in the absence of prompt SR overpasses. We found that a moving average approach was most suitable for that purpose, although limited by the absence of explicit records of radar maintenance operations.
River-valley morphology preserves information on tectonic and climatic conditions that shape landscapes. Observations suggest that river discharge and valley-wall lithology are the main controls on valley width. Yet, current models based on these observations fail to explain the full range of cross-sectional valley shapes in nature, suggesting hitherto unquantified controls on valley width. In particular, current models cannot explain the existence of paired terrace sequences that form under cyclic climate forcing. Paired river terraces are staircases of abandoned floodplains on both valley sides, and hence preserve past valley widths. Their formation requires alternating phases of predominantly river incision and predominantly lateral planation, plus progressive valley narrowing. While cyclic Quaternary climate changes can explain shifts between incision and lateral erosion, the driving mechanism of valley narrowing is unknown. Here, we extract valley geometries from climatically formed, alluvial river-terrace sequences and show that across our dataset, the total cumulative terrace height (here: total valley height) explains 90%–99% of the variance in valley width at the terrace sites. This finding suggests that valley height, or a parameter that scales linearly with valley height, controls valley width in addition to river discharge and lithology. To explain this valley-width-height relationship, we reformulate existing valley-width models and suggest that, when adjusting to new boundary conditions, alluvial valleys evolve to a width at which sediment removal from valley walls matches lateral sediment supply from hillslope erosion. Such a hillslope-channel coupling is not captured in current valley-evolution models. Our model can explain the existence of paired terrace sequences under cyclic climate forcing and relates valley width to measurable field parameters. Therefore, it facilitates the reconstruction of past climatic and tectonic conditions from valley topography.
Sea level rise and coastal erosion have inundated large areas of Arctic permafrost. Submergence by warm and saline waters increases the rate of inundated permafrost thaw compared to sub-aerial thawing on land. Studying the contact between the unfrozen and frozen sediments below the seabed, also known as the ice-bearing permafrost table (IBPT), provides valuable information to understand the evolution of sub-aquatic permafrost, which is key to improving and understanding coastal erosion prediction models and potential greenhouse gas emissions. In this study, we use data from 2D electrical resistivity tomography (ERT) collected in the nearshore coastal zone of two Arctic regions that differ in their environmental conditions (e.g., seawater depth and resistivity) to image and study the subsea permafrost. The inversion of 2D ERT data sets is commonly performed using deterministic approaches that favor smoothed solutions, which are typically interpreted using a user-specified resistivity threshold to identify the IBPT position. In contrast, to target the IBPT position directly during inversion, we use a layer-based model parameterization and a global optimization approach to invert our ERT data. This approach results in ensembles of layered 2D model solutions, which we use to identify the IBPT and estimate the resistivity of the unfrozen and frozen sediments, including estimates of uncertainties. Additionally, we globally invert 1D synthetic resistivity data and perform sensitivity analyses to study, in a simpler way, the correlations and influences of our model parameters. The set of methods provided in this study may help to further exploit ERT data collected in such permafrost environments as well as for the design of future field experiments.
The Ice, Cloud, and Land Elevation Satellite-2 (ICESat-2) with its land and vegetation height data product (ATL08), and Global Ecosystem Dynamics Investigation (GEDI) with its terrain elevation and height metrics data product (GEDI Level 2A) missions have great potential to globally map ground and canopy heights. Canopy height is a key factor in estimating above-ground biomass and its seasonal changes; these satellite missions can also improve estimated above-ground carbon stocks. This study presents a novel Sparse Vegetation Detection Algorithm (SVDA) which uses ICESat-2 (ATL03, geolocated photons) data to map tree and vegetation heights in a sparsely vegetated savanna ecosystem. The SVDA consists of three main steps: First, noise photons are filtered using the signal confidence flag from ATL03 data and local point statistics. Second, we classify ground photons based on photon height percentiles. Third, tree and grass photons are classified based on the number of neighbors. We validated tree heights with field measurements (n = 55), finding a root-mean-square error (RMSE) of 1.82 m using SVDA, GEDI Level 2A (Geolocated Elevation and Height Metrics product): 1.33 m, and ATL08: 5.59 m. Our results indicate that the SVDA is effective in identifying canopy photons in savanna ecosystems, where ATL08 performs poorly. We further identify seasonal vegetation height changes with an emphasis on vegetation below 3 m; widespread height changes in this class from two wet-dry cycles show maximum seasonal changes of 1 m, possibly related to seasonal grass-height differences. Our study shows the difficulties of vegetation measurements in savanna ecosystems but provides the first estimates of seasonal biomass changes.
Cosmic-ray neutron sensing (CRNS) has become an effective method to measure soil moisture at a horizontal scale of hundreds of metres and a depth of decimetres. Recent studies proposed operating CRNS in a network with overlapping footprints in order to cover root-zone water dynamics at the small catchment scale and, at the same time, to represent spatial heterogeneity. In a joint field campaign from September to November 2020 (JFC-2020), five German research institutions deployed 15 CRNS sensors in the 0.4 km2 Wüstebach catchment (Eifel mountains, Germany). The catchment is dominantly forested (but includes a substantial fraction of open vegetation) and features a topographically distinct catchment boundary. In addition to the dense CRNS coverage, the campaign featured a unique combination of additional instruments and techniques: hydro-gravimetry (to detect water storage dynamics also below the root zone); ground-based and, for the first time, airborne CRNS roving; an extensive wireless soil sensor network, supplemented by manual measurements; and six weighable lysimeters. Together with comprehensive data from the long-term local research infrastructure, the published data set (available at https://doi.org/10.23728/b2share.756ca0485800474e9dc7f5949c63b872; Heistermann et al., 2022) will be a valuable asset in various research contexts: to advance the retrieval of landscape water storage from CRNS, wireless soil sensor networks, or hydrogravimetry; to identify scale-specific combinations of sensors and methods to represent soil moisture variability; to improve the understanding and simulation of land–atmosphere exchange as well as hydrological and hydrogeological processes at the hillslope and the catchment scale; and to support the retrieval of soil water content from airborne and spaceborne remote sensing platforms.
In this paper, we present the new frequency spectrum recurrence analysis technique by means of electro-encephalon signals (EES) analyses. The technique is suitable for time series analysis with noise and disturbances. EES were collected, and alpha waves of the occipital region were analysed by comparing the signals from participants in two states, eyes open and eyes closed. Firstly, EES were characterized and analysed by means of techniques already known to compare with the results of the innovative technique that we present here. We verified that, standard recurrence quantification analysis by means of EES time series cannot statistically distinguish the two states. However, the new frequency spectrum recurrence quantification exhibit quantitatively whether the participants have their eyes open or closed. In sequence, new quantifiers are created for analysing the recurrence concentration on frequency bands. These analyses show that EES with similar frequency spectrum have different recurrence levels revealing different behaviours of the nervous system. The technique can be used to deepen the study on depression, stress, concentration level and other neurological issues and also can be used in any complex system.
AM(W)= 5.1 earthquake on January 21st, 2016 marked the beginning of a significant seismic sequence in the southern Alboran Sea, culminating in aM(W)= 6.3 earthquake on January 25th, and continuing with further moderate magnitude earthquakes until March. We use data from 35 seismic broadband stations in Spain, Morocco and Portugal to relocate the seismicity, estimate seismic moment tensors, and isolate regional apparent source time functions for the main earthquake. Relocation and regional moment tensor inversion consistently yield very shallow depths for the majority of events. We obtain 50 moment tensors for the sequence, showing a mixture of strike-slip faulting for the foreshock and the main event and reverse faulting for the major aftershocks. The leading role of reverse focal mechanisms among the aftershocks may be explained by the geometry of the fault network. The mainshock nucleates at a bend along the left-lateral Al-Idrisi fault, introducing local transpression within the transtensional Alboran Basin. The shallow depths of the 2016 Alboran Sea earthquakes may favor slip-partitioning on the involved faults. Apparent source durations for the main event suggest a similar to 21 km long, asymmetric rupture that propagates primarily toward NE into the restraining fault segment, with fast rupture speed of similar to 3.0 km/s. Consistently, the inversion for laterally variable fault displacement situates the main slip in the restraining segment. The partitioning into strike-slip rupture and dip-slip aftershocks confirms a non-optimal orientation of this segment, and suggests that the 2016 event settled a slip deficit from previous ruptures that could not propagate into the stronger restraining segment.
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.
This review provides a synthesis of current knowledge on the morphological and functional traits of testate amoebae, a polyphyletic group of protists commonly used as proxies of past hydrological changes in paleoecological investigations from peatland, lake sediment and soil archives. A trait-based approach to understanding testate amoebae ecology and paleoecology has gained in popularity in recent years, with research showing that morphological characteristics provide complementary information to the commonly used environmental inferences based on testate amoeba (morpho-)species data. We provide a broad overview of testate amoeba morphological and functional traits and trait-environment relationships in the context of ecology, evolution, genetics, biogeography, and paleoecology. As examples we report upon previous ecological and paleoecological studies that used trait-based approaches, and describe key testate amoebae traits that can be used to improve the interpretation of environmental studies. We also highlight knowledge gaps and speculate on potential future directions for the application of trait-based approaches in testate amoeba research.
The origins and development of the arid and highly seasonal steppe-desert biome in Central Asia, the largest of its kind in the world, remain largely unconstrained by existing records. It is unclear how Cenozoic climatic, geological, and biological forces, acting at diverse spatial and temporal scales, shaped Central Asian ecosystems through time. Our synthesis shows that the Central Asian steppe-desert has existed since at least Eocene times but experienced no less than two regime shifts, one at the Eocene-Oligocene Transition and one in the mid-Miocene. These shifts separated three successive "stable states," each characterized by unique floral and faunal structures. Past responses to disturbance in the Asian steppe-desert imply that modern ecosystems are unlikely to recover their present structures and diversity if forced into a new regime. This is of concern for Asian steppes today, which are being modified for human use and lost to desertification at unprecedented rates.
The Arctic is greatly affected by climate change. Increasing air temperatures drive permafrost thaw and an increase in coastal erosion and river discharge. This results in a greater input of sediment and organic matter into nearshore waters, impacting ecosystems by reducing light transmission through the water column and altering biogeochemistry. This potentially results in impacts on the subsistence economy of local people as well as the climate due to the transformation of suspended organic matter into greenhouse gases. Even though the impacts of increased suspended sediment concentrations and turbidity in the Arctic nearshore zone are well-studied, the mechanisms underpinning this increase are largely unknown. Wave energy and tides drive the level of turbidity in the temperate and tropical parts of the world, and this is generally assumed to also be the case in the Arctic. However, the tidal range is considerably lower in the Arctic, and processes related to the occurrence of permafrost have the potential to greatly contribute to nearshore turbidity. In this study, we use high-resolution satellite imagery alongside in situ and ERA5 reanalysis data of ocean and climate variables in order to identify the drivers of nearshore turbidity, along with its seasonality in the nearshore waters of Herschel Island Qikiqtaruk, in the western Canadian Arctic. Nearshore turbidity correlates well to wind direction, wind speed, significant wave height, and wave period. Nearshore turbidity is superiorly correlated to wind speed at the Beaufort Shelf compared to in situ measurements at Herschel Island Qikiqtaruk, showing that nearshore turbidity, albeit being of limited spatial extent, is influenced by large-scale weather and ocean phenomenons. We show that, in contrast to the temperate and tropical ocean, freshly eroded material is the predominant driver of nearshore turbidity in the Arctic, rather than resuspension, which is caused by the vulnerability of permafrost coasts to thermo-erosion.
The Arctic is greatly impacted by climate change. The increase in air temperature drives the thawing of permafrost and an increase in coastal erosion and river discharge. This leads to a greater input of sediment and organic matter into coastal waters, which substantially impacts the ecosystems by reducing light transmission through the water column and altering the biogeochemistry, but also the subsistence economy of local people, and changes in climate because of the transformation of organic matter into greenhouse gases. Yet, the quantification of suspended sediment in Arctic coastal and nearshore waters remains unsatisfactory due to the absence of dedicated algorithms to resolve the high loads occurring in the close vicinity of the shoreline. In this study we present the Arctic Nearshore Turbidity Algorithm (ANTA), the first reflectance-turbidity relationship specifically targeted towards Arctic nearshore waters that is tuned with in-situ measurements from the nearshore waters of Herschel Island Qikiqtaruk in the western Canadian Arctic. A semi-empirical model was calibrated for several relevant sensors in ocean color remote sensing, including MODIS, Sentinel 3 (OLCI), Landsat 8 (OLI), and Sentinel 2 (MSI), as well as the older Landsat sensors TM and ETM+. The ANTA performed better with Landsat 8 than with Sentinel 2 and Sentinel 3. The application of the ANTA to Sentinel 2 imagery that matches in-situ turbidity samples taken in Adventfjorden, Svalbard, shows transferability to nearshore areas beyond Herschel Island Qikiqtaruk.
Soziale Medien sind ein wesentlicher Bestandteil des Alltags von Schüler*innen und gleichzeitig zunehmend wichtig in Wirtschaft, Politik und Wissenschaft. Am Beispiel von Twitter zeigt dieser Beitrag, dass soziale Medien im Unterricht auch für die Beantwortung geographischer Fragestellungen verwendet werden können. Hierfür eignen sich Twitter-Daten aufgrund ihrer Georeferenzierung und weiterer interessanter Inhalte besonders. Der Beitrag gibt einen Überblick über die Verwendung von Twitter für sozialwissenschaftliche und humangeographische Fragestellungen und reflektiert die Nutzung von Twitter im Unterricht. Für die Unterrichtspraxis werden Beispiele zu den Themen Braunkohle, Flutereignisse und Raumwahrnehmungen sowie Anleitungen zur Auswertung, Anwendung und Reflexion von Twitter-Analysen vorgestellt.
Quantitative geomorphic research depends on accurate topographic data often collected via remote sensing. Lidar, and photogrammetric methods like structure-from-motion, provide the highest quality data for generating digital elevation models (DEMs). Unfortunately, these data are restricted to relatively small areas, and may be expensive or time-consuming to collect. Global and near-global DEMs with 1 arcsec (∼30 m) ground sampling from spaceborne radar and optical sensors offer an alternative gridded, continuous surface at the cost of resolution and accuracy. Accuracy is typically defined with respect to external datasets, often, but not always, in the form of point or profile measurements from sources like differential Global Navigation Satellite System (GNSS), spaceborne lidar (e.g., ICESat), and other geodetic measurements. Vertical point or profile accuracy metrics can miss the pixel-to-pixel variability (sometimes called DEM noise) that is unrelated to true topographic signal, but rather sensor-, orbital-, and/or processing-related artifacts. This is most concerning in selecting a DEM for geomorphic analysis, as this variability can affect derivatives of elevation (e.g., slope and curvature) and impact flow routing. We use (near) global DEMs at 1 arcsec resolution (SRTM, ASTER, ALOS, TanDEM-X, and the recently released Copernicus) and develop new internal accuracy metrics to assess inter-pixel variability without reference data. Our study area is in the arid, steep Central Andes, and is nearly vegetation-free, creating ideal conditions for remote sensing of the bare-earth surface. We use a novel hillshade-filtering approach to detrend long-wavelength topographic signals and accentuate short-wavelength variability. Fourier transformations of the spatial signal to the frequency domain allows us to quantify: 1) artifacts in the un-projected 1 arcsec DEMs at wavelengths greater than the Nyquist (twice the nominal resolution, so > 2 arcsec); and 2) the relative variance of adjacent pixels in DEMs resampled to 30-m resolution (UTM projected). We translate results into their impact on hillslope and channel slope calculations, and we highlight the quality of the five DEMs. We find that the Copernicus DEM, which is based on a carefully edited commercial version of the TanDEM-X, provides the highest quality landscape representation, and should become the preferred DEM for topographic analysis in areas without sufficient coverage of higher-quality local DEMs.
The intensification of Northern Hemisphere glaciations at the end of the Pliocene epoch marks one of the most substantial climatic shifts of the Cenozoic. Despite global cooling, sea surface temperatures in the high latitude North Atlantic Ocean rose between 2.9–2.7 million years ago. Here we present sedimentary geochemical proxy data from the Gulf of Cadiz to reconstruct the variability of Mediterranean Outflow Water, an important heat source to the North Atlantic. We find evidence for enhanced production of Mediterranean Outflow from the mid-Pliocene to the late Pliocene which we infer could have driven a sub-surface heat channel into the high-latitude North Atlantic. We then use Earth System Models to constrain the impact of enhanced Mediterranean Outflow production on the northward heat transport in the North Atlantic. In accord with the proxy data, the numerical model results support the formation of a sub-surface channel that pumped heat from the subtropics into the high latitude North Atlantic. We further suggest that this mechanism could have delayed ice sheet growth at the end of the Pliocene.
Sediment archives in the terrestrial and marine realm are regularly analyzed to infer changes in climate, tectonic, or anthropogenic boundary conditions of the past. However, contradictory observations have been made regarding whether short period events are faithfully preserved in stratigraphic archives; for instance, in marine sediments offshore large river systems. On the one hand, short period events are hypothesized to be non-detectable in the signature of terrestrially derived sediments due to buffering during sediment transport along large river systems. On the other hand, several studies have detected signals of short period events in marine records offshore large river systems. We propose that this apparent discrepancy is related to the lack of a differentiation between different types of signals and the lack of distinction between river response times and signal propagation times. In this review, we (1) expand the definition of the term ‘signal’ and group signals in sub-categories related to hydraulic grain size characteristics, (2) clarify the different types of ‘times’ and suggest a precise and consistent terminology for future use, and (3) compile and discuss factors influencing the times of signal transfer along sediment routing systems and how those times vary with hydraulic grain size characteristics. Unraveling different types of signals and distinctive time periods related to signal propagation addresses the discrepancies mentioned above and allows a more comprehensive exploration of event preservation in stratigraphy – a prerequisite for reliable environmental reconstructions from terrestrially derived sedimentary records.
The Aral Sea desiccation and related changes in hydroclimatic conditions on a regional level is a hot topic for past decades. The key problem of scientific research projects devoted to an investigation of modern Aral Sea basin hydrological regime is its discontinuous nature - the only limited amount of papers takes into account the complex runoff formation system entirely. Addressing this challenge we have developed a continuous prediction system for assessing freshwater inflow into the Small Aral Sea based on coupling stack of hydrological and data-driven models. Results show a good prediction skill and approve the possibility to develop a valuable water assessment tool which utilizes the power of classical physically based and modern machine learning models both for territories with complex water management system and strong water-related data scarcity. The source code and data of the proposed system is available on a Github page (https://github.com/SMASHIproject/IWRM2018).
Identifying abrupt transitions is a key question in various disciplines. Existing transition detection methods, however, do not rigorously account for time series uncertainties, often neglecting them altogether or assuming them to be independent and qualitatively similar. Here, we introduce a novel approach suited to handle uncertainties by representing the time series as a time-ordered sequence of probability density functions. We show how to detect abrupt transitions in such a sequence using the community structure of networks representing probabilities of recurrence. Using our approach, we detect transitions in global stock indices related to well-known periods of politico-economic volatility. We further uncover transitions in the El Niño-Southern Oscillation which coincide with periods of phase locking with the Pacific Decadal Oscillation. Finally, we provide for the first time an ‘uncertainty-aware’ framework which validates the hypothesis that ice-rafting events in the North Atlantic during the Holocene were synchronous with a weakened Asian summer monsoon.
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.
This study analyzes the influence of local and regional climatic factors on the stable isotopic composition of rainfall in the Vietnamese Mekong Delta (VMD) as part of the Asian monsoon region. It is based on 1.5 years of weekly rainfall samples. In the first step, the isotopic composition of the samples is analyzed by local meteoric water lines (LMWLs) and single-factor linear correlations. Additionally, the contribution of several regional and local factors is quantified by multiple linear regression (MLR) of all possible factor combinations and by relative importance analysis. This approach is novel for the interpretation of isotopic records and enables an objective quantification of the explained variance in isotopic records for individual factors. In this study, the local factors are extracted from local climate records, while the regional factors are derived from atmospheric backward trajectories of water particles. The regional factors, i.e., precipitation, temperature, relative humidity and the length of backward trajectories, are combined with equivalent local climatic parameters to explain the response variables delta O-18, delta H-2, and d-excess of precipitation at the station of measurement. The results indicate that (i) MLR can better explain the isotopic variation in precipitation (R-2 = 0.8) compared to single-factor linear regression (R-2 = 0.3); (ii) the isotopic variation in precipitation is controlled dominantly by regional moisture regimes (similar to 70 %) compared to local climatic conditions (similar to 30 %); (iii) the most important climatic parameter during the rainy season is the precipitation amount along the trajectories of air mass movement; (iv) the influence of local precipitation amount and temperature is not sig-nificant during the rainy season, unlike the regional precipitation amount effect; (v) secondary fractionation processes (e.g., sub-cloud evaporation) can be identified through the d-excess and take place mainly in the dry season, either locally for delta O-18 and delta H-2, or along the air mass trajectories for d-excess. The analysis shows that regional and local factors vary in importance over the seasons and that the source regions and transport pathways, and particularly the climatic conditions along the pathways, have a large influence on the isotopic composition of rainfall. Although the general results have been reported qualitatively in previous studies (proving the validity of the approach), the proposed method provides quantitative estimates of the controlling factors, both for the whole data set and for distinct seasons. Therefore, it is argued that the approach constitutes an advancement in the statistical analysis of isotopic records in rainfall that can supplement or precede more complex studies utilizing atmospheric models. Due to its relative simplicity, the method can be easily transferred to other regions, or extended with other factors. The results illustrate that the interpretation of the isotopic composition of precipitation as a recorder of local climatic conditions, as for example performed for paleorecords of water isotopes, may not be adequate in the southern part of the Indochinese Peninsula, and likely neither in other regions affected by monsoon processes. However, the presented approach could open a pathway towards better and seasonally differentiated reconstruction of paleoclimates based on isotopic records.
Ice-rich yedoma-dominated landscapes store considerable amounts of organic carbon (C) and nitrogen (N) and are vulnerable to degradation under climate warming. We investigate the C and N pools in two thermokarst-affected yedoma landscapes - on Sobo-Sise Island and on Bykovsky Peninsula in the north of eastern Siberia. Soil cores up to 3m depth were collected along geomorphic gradients and analysed for organic C and N contents. A high vertical sampling density in the profiles allowed the calculation of C and N stocks for short soil column intervals and enhanced understanding of within-core parameter variability. Profile-level C and N stocks were scaled to the landscape level based on landform classifications from 5 m resolution, multispectral RapidEye satellite imagery. Mean landscape C and N storage in the first metre of soil for Sobo-Sise Island is estimated to be 20.2 kg C m(-2) and 1.8 kg N m(-2) and for Bykovsky Peninsula 25.9 kg C m(-2) and 2.2 kg N m(-2). Radiocarbon dating demonstrates the Holocene age of thermokarst basin deposits but also suggests the presence of thick Holoceneage cover layers which can reach up to 2 m on top of intact yedoma landforms. Reconstructed sedimentation rates of 0.10-0.57 mm yr(-1) suggest sustained mineral soil accumulation across all investigated landforms. Both yedoma and thermokarst landforms are characterized by limited accumulation of organic soil layers (peat). We further estimate that an active layer deepening of about 100 cm will increase organic C availability in a seasonally thawed state in the two study areas by similar to 5.8 Tg (13.2 kg C m(-2)). Our study demonstrates the importance of increasing the number of C and N storage inventories in ice-rich yedoma and thermokarst environments in order to account for high variability of permafrost and thermokarst environments in pan-permafrost soil C and N pool estimates.
Background. Metal recycling factories (MRFs) have developed rapidly in Nigeria as recycling policies have been increasingly embraced. These MRFs are point sources for introducing potentially toxic elements (PTEs) into environmental media. Objectives. The aim of this study was to determine the constituents (elemental and mineralogy) of the wastes (slag and particulate matter, (PM)) and soils around the MRFs and to determine the level of pollution within the area. Methods. Sixty samples (30 slag samples, 15 soil samples and 15 PM samples) were collected for this study. The soils, slag and PM samples were analyzed for elemental constituents using inductively coupled plasma optical emission spectrometry. Mineralogy of the PM was determined using scanning electron microscope-energy dispersive spectroscopy (SEM-EDS), and soil mineralogy was determined by an X-ray diffractometer (XRD). Results. The results of the soil analyses revealed the following concentrations for the selected metals in mg/kg include lead (Pb) (21.0-2399.0), zinc (Zn) (56.0-4188.0), copper (Cu) (10.0-1470.0), nickel (Ni) (6.0-215.0), chromium (Cr) (921.0-1737.0) and cadmium (Cd) (below detectable limit (Bdl)-18.1). For the slags the results were Pb (68.0-.333.0), Zn (1364.0-3062), Cu (119.0-1470.0), Ni (12.0-675.0), Cr (297-1737) and Cd (Bdl-15.8). The results in mu g/g for the metal analysis in PM were Pb (4.6-160.0), Zn (18.0-471.0), Cu (2.5-11.0), Ni (0.8-4.2), and Cr (2.5-11.0), while Cd was undetected. The slags are currently utilized for filling the foundations of buildings and roads, providing additional pathways for the introduction of PTEs into the environment from the suspended materials generated from mechanical breakdown of the slags. Conclusions. The MRFs were found to have impacted the quality of environmental media through the introduction of PTEs, impairing soil quality, in addition to PM, which can have detrimental health consequences. Further studies on the health implications of these pollutants and their impacts on human health are needed. Competing Interests. The authors declare no competing financial interests
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.
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.