Refine
Year of publication
Document Type
- Article (2817) (remove)
Keywords
- Holocene (35)
- climate change (30)
- Earthquake source observations (21)
- Climate change (18)
- Pollen (18)
- erosion (18)
- permafrost (18)
- Himalaya (17)
- Tibetan Plateau (17)
- Seismicity and tectonics (15)
Institute
- Institut für Geowissenschaften (2817) (remove)
Quantifying the resilience of vegetated ecosystems is key to constraining both present-day and future global impacts of anthropogenic climate change. Here we apply both empirical and theoretical resilience metrics to remotely-sensed vegetation data in order to examine the role of water availability and variability in controlling vegetation resilience at the global scale. We find a concise global relationship where vegetation resilience is greater in regions with higher water availability. We also reveal that resilience is lower in regions with more pronounced inter-annual precipitation variability, but find less concise relationships between vegetation resilience and intra-annual precipitation variability. Our results thus imply that the resilience of vegetation responds differently to water deficits at varying time scales. In view of projected increases in precipitation variability, our findings highlight the risk of ecosystem degradation under ongoing climate change.
Vegetation dynamics depend on both the amount of precipitation and its variability over time. Here, the authors show that vegetation resilience is greater where water availability is higher and where precipitation is more stable from year to year.
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.
During the last 5 Ma the Earth's ocean-atmosphere system passed through several major transitions, many of which are discussed as possible triggers for human evolution. A classic in this context is the possible influence of the closure of the Panama Strait, the intensification of Northern Hemisphere Glaciation, a stepwise increase in aridity in Africa, and the first appearance of the genus Homo about 2.5 - 2.7 Ma ago. Apart from the fact that the correlation between these events does not necessarily imply causality, many attempts to establish a relationship between climate and evolution fail due to the challenge of precisely localizing an a priori unknown number of changes potentially underlying complex climate records. The kernel-based Bayesian inference approach applied here allows inferring the location, generic shape, and temporal scale of multiple transitions in established records of Plio-Pleistocene African climate. By defining a transparent probabilistic analysis strategy, we are able to identify conjoint changes occurring across the investigated terrigenous dust records from Ocean Drilling Programme (ODP) sites in the Atlantic Ocean (ODP 659), Arabian (ODP 721/722) and Mediterranean Sea (ODP 967). The study indicates a two-step transition in the African climate proxy records at (2.35-2.10) Ma and (1.70 - 1.50) Ma, that may be associated with the reorganization of the Hadley-Walker Circulation. .
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.
A volcanic eruption is usually preceded by seismic precursors, but their interpretation and use for forecasting the eruption onset time remain a challenge. A part of the eruptive processes in open conduits of volcanoes may be similar to those encountered in geysers. Since geysers erupt more often, they are useful sites for testing new forecasting methods. We tested the application of Permutation Entropy (PE) as a robust method to assess the complexity in seismic recordings of the Strokkur geyser, Iceland. Strokkur features several minute-long eruptive cycles, enabling us to verify in 63 recorded cycles whether PE behaves consistently from one eruption to the next one. We performed synthetic tests to understand the effect of different parameter settings in the PE calculation. Our application to Strokkur shows a distinct, repeating PE pattern consistent with previously identified phases in the eruptive cycle. We find a systematic increase in PE within the last 15 s before the eruption, indicating that an eruption will occur. We quantified the predictive power of PE, showing that PE performs better than seismic signal strength or quiescence when it comes to forecasting eruptions.
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 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.
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.
Neoarchean (similar to 2.73-2.70 Ga) accretionary history of the eastern Dharwar Craton, India
(2022)
Cratonic mid-crustal plutons may contain supracrustal enclaves that preserve evidence of an earlier growth history. The Eastern Dharwar craton records Neoarchean two-stage accretionary sequential growth (2.70 and 2.55 Ga) and a chronology of their enclaves could refine orogenic models. To test whether the metamorphic history of their enclaves was related to any of these stages, phase equilibria modelling and combined Lu-Hf and Sm-Nd geochronology on garnet were conducted on metapsammite, now preserved as garnet-orthopyroxene-cordierite gneiss. Phase equilibria modelling indicates peak metamorphic conditions, similar to 850 degrees C and similar to 8.5 kbar (M1a), were followed by near isothermal decompression to 5-6 kbar (M1b) and isobaric cooling to similar to 800 degrees C (M1c). The thermobaric gradient related to peak metamorphic conditions, similar to 30 degrees C kbar(-1), is typical of collisional orogens. Regression of the whole-rock and garnet, for sample S17b, yield Lu-Hf isochron ages of 2733 +/- 29 Ma, and for sample S18, 2724 +/- 13 Ma. A Lu-Hf weighted mean age for the porphyroblastic garnet suggests growth at 2725.5 +/- 11.9 Ma during the M1a-M1b stages. In contrast, the whole-rock sample S17b and the garnet fractions yield a Sm-Nd isochron age of 2696 +/- 10 Ma. From sample S18 the whole rock, garnet fractions, and orthopyroxene yield an isochron age of 2683 +/- 15 Ma. The garnet Sm-Nd weighted mean age at 2692.0 +/- 8.3 Ma constrains the M1b-M1c stages. We suggest that the protoliths to these supracrustal enclaves were deposited in an arc tectonic setting and underwent thickening followed by heating during peeled-back lithospheric convergence. Therefore, the earliest of the craton-forming accretionary stages is preserved as the similar to 2.73 Ga granulite-facies enclaves, marginally older than the 2.70-2.65 Ga cratonic greenstone volcanism. Tectonic exhumation of these mid-crustal granulite enclaves was in response to the late-Proterozoic (similar to 1.7 Ga) Bhopalpatnam orogeny.
The geometry of carbonate platforms reflects the interaction of several factors. However, the impact of carbonate-producing organisms has been poorly investigated so far. This study applies stratigraphic forward modelling (SFM) and sensitivity analysis to examine, referenced to the Miocene Llucmajor Platform, the effect of changes of dominant biotic production in the oligophotic and euphotic zones on platform geometry. Our results show that the complex interplay of carbonate production rates, bathymetry and variations in accommodation space control the platform geometry. The main driver of progradation is the oligophotic production of rhodalgal sediments during the lowstands. This study demonstrates that platform geometry and internal architecture varies significantly according to the interaction of the predominant carbonate-producing biotas. The input parameters for this study are based on well-understood Miocene carbonate biotas with characteristic euphotic, oligophotic and photo-independent carbonate production in which it is crucial that each carbonate-producing class is modelled explicitly within the simulation run and not averaged with a single carbonate production-depth profile. This is important in subsurface exploration studies based on stratigraphic forward models where the overall platform geometry may be approximated through calibration runs, and constrained by seismic surveys and wellbores. However, the internal architecture is likely to be oversimplified without an in-depth understanding of the target carbonate system and a transfer to forward modelling parameters.
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.
The Mackenzie Delta (MD) is a permafrost-bearing region along the coasts of the Canadian Arctic which exhibits high sub-permafrost gas hydrate (GH) reserves. The GH occurring at the Mallik site in the MD is dominated by thermogenic methane (CH4), which migrated from deep conventional hydrocarbon reservoirs, very likely through the present fault systems. Therefore, it is assumed that fluid flow transports dissolved CH4 upward and out of the deeper overpressurized reservoirs via the existing polygonal fault system and then forms the GH accumulations in the Kugmallit-Mackenzie Bay Sequences. We investigate the feasibility of this mechanism with a thermo-hydraulic-chemical numerical model, representing a cross section of the Mallik site. We present the first simulations that consider permafrost formation and thawing, as well as the formation of GH accumulations sourced from the upward migrating CH4-rich formation fluid. The simulation results show that temperature distribution, as well as the thickness and base of the ice-bearing permafrost are consistent with corresponding field observations. The primary driver for the spatial GH distribution is the permeability of the host sediments. Thus, the hypothesis on GH formation by dissolved CH4 originating from deeper geological reservoirs is successfully validated. Furthermore, our results demonstrate that the permafrost has been substantially heated to 0.8-1.3 degrees C, triggered by the global temperature increase of about 0.44 degrees C and further enhanced by the Arctic Amplification effect at the Mallik site from the early 1970s to the mid-2000s.
The 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.
On 7 January 2020, an M-w 6.4 earthquake occurred in the northeastern Caribbean, a few kilometers offshore of the island of Puerto Rico. It was the mainshock of a complex seismic sequence, characterized by a large number of energetic earthquakes illuminating an east-west elongated area along the southwestern coast of Puerto Rico. Deformation fields constrained by Interferometric Synthetic Aperture Radar and Global Navigation Satellite System data indicate that the coseismic movements affected only the western part of the island. To assess the mainshock's source fault parameters, we combined the geodetically derived coseismic deformation with teleseismic waveforms using Bayesian inference. The results indicate a roughly east-west oriented fault, dipping northward and accommodating similar to 1.4 m of transtensional motion. Besides, the determined location and orientation parameters suggest an offshore continuation of the recently mapped North Boqueron Bay-Punta Montalva fault in southwest Puerto Rico. This highlights the existence of unmapped faults with moderate-to-large earthquake potential within the Puerto Rico region.
Despite more than half a century of hominin fossil discoveries in eastern Africa, the regional environmental context of hominin evolution and dispersal is not well established due to the lack of continuous palaeoenvironmental records from one of the proven habitats of early human populations, particularly for the Pleistocene epoch. Here we present a 620,000-year environmental record from Chew Bahir, southern Ethiopia, which is proximal to key fossil sites. Our record documents the potential influence of different episodes of climatic variability on hominin biological and cultural transformation. The appearance of high anatomical diversity in hominin groups coincides with long-lasting and relatively stable humid conditions from similar to 620,000 to 275,000 years bp (episodes 1-6), interrupted by several abrupt and extreme hydroclimate perturbations. A pattern of pronounced climatic cyclicity transformed habitats during episodes 7-9 (similar to 275,000-60,000 years bp), a crucial phase encompassing the gradual transition from Acheulean to Middle Stone Age technologies, the emergence of Homo sapiens in eastern Africa and key human social and cultural innovations. Those accumulative innovations plus the alignment of humid pulses between northeastern Africa and the eastern Mediterranean during high-frequency climate oscillations of episodes 10-12 (similar to 60,000-10,000 years bp) could have facilitated the global dispersal of H. sapiens.
We use a dense seismic network on the Reykjanes Peninsula, Iceland, to image a group of earthquakes at 10-12 km depth, 2 km north-east of 2021 Fagradalsfjall eruption site. These deep earthquakes have a lower frequency content compared to earthquakes located in the upper, brittle crust and are similar to deep long period (DLP) seismicity observed at other volcanoes in Iceland and around the world. We observed several swarms of DLP earthquakes between the start of the study period (June 2020) and the initiation of the 3-week-long dyke intrusion that preceded the eruption in March 2021. During the eruption, DLP earthquake swarms returned 1 km SW of their original location during periods when the discharge rate or fountaining style of the eruption changed. The DLP seismicity is therefore likely to be linked to the magma plumbing system beneath Fagradalsfjall. However, the DLP seismicity occurred similar to 5 km shallower than where petrological modelling places the near-Moho magma storage region in which the Fagradalsfjall lava was stored. We suggest that the DLP seismicity was triggered by the exsolution of CO2-rich fluids or the movement of magma at a barrier to the transport of melt in the lower crust. Increased flux through the magma plumbing system during the eruption likely adds to the complexity of the melt migration process, thus causing further DLP seismicity, despite a contemporaneous magma channel to the surface.
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.
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.
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.
Fifteen N-butylpyridinium salts - five monometallic [C4Py](2)[MBr4] and ten bimetallic [C4Py](2)[(M0.5M0.5Br4)-M-a-Br-b] (M=Co, Cu, Mn, Ni, Zn) - were synthesized, and their structures and thermal and electrochemical properties were studied. All the compounds are ionic liquids (ILs) with melting points between 64 and 101 degrees C. Powder and single-crystal X-ray diffraction show that all ILs are isostructural. The electrochemical stability windows of the ILs are between 2 and 3 V. The conductivities at room temperature are between 10(-5) and 10(-6) S cm(-1). At elevated temperatures, the conductivities reach up to 10(-4) S cm(-1) at 70 degrees C. The structures and properties of the current bromide-based ILs were also compared with those of previous examples using chloride ligands, which illustrated differences and similarities between the two groups of ILs.
Devolatilization of subducting lithologies liberates COH-fluids. These may become partially sequestered in peridotites in the slab and the overlying forearc mantle, affecting the cycling of volatiles and fluid mobile elements in subduction zones. Here we assess the magnitudes, timescales and mechanism of channelized injection of COH-fluids doped with Ca-aq(2+), Sr-aq(2+) and Ba-aq(2+) into the dry forearc mantle by performing piston cylinder experiments between 1-2.5 GPa and 600-700 degrees C. Cylindrical cores of natural spinel-bearing harzburgites were used as starting materials. Based on mineral assemblage and composition three reaction zones are distinguishable from the rim towards the core of primary olivine and orthopyroxene grains. Zone 1 contains carbonates + quartz +/- kyanite and zone 2 contains carbonates + talc +/- chlorite. Olivine is further replaced in zone 3 by either antigorite+ magnesite or magnesite +talc within or above antigorite stability, respectively. Orthopyroxene is replaced in zone 3 by talc + chlorite. Mineral assemblages and the compositions of secondary minerals depend on fluid composition and the replaced primary silicate. The extent of alteration depends on fluid CO2 content and fluid/rock-ratio, and is further promoted by fluid permeable reaction zones and reaction driven cracking. Our results show that COH-fluid induced metasomatism of the forearc mantle is self-perpetuating and efficient at sequestering Ca-aq(2+), Sr-aq(2+), Ba-aq(2+) and CO2aq into newly formed carbonates. This process is fast with 90% of the available C sequestered and nearly 50% of the initial minerals altered at 650 degrees C, 2 GPa within 55 h. The dissolution of primary silicates under high COH-fluid/rock-ratios, as in channelized fluid flow, enriches SiO2aq in the fluid, while CO2aq is sequestered into carbonates. In an open system, the remaining CO2-depleted, Si-enriched aqueous fluid may cause Si-metasomatism in the forearc further away from the injection of the COH-fluid into peridotite.
ABSTRACT: Structural evolution of cesium triiodide at high pressures has been revealed by synchrotron single-crystal X-ray diffraction. Cesium triiodide undergoes a first-order phase transition above 1.24(3) GPa from an orthorhombic to a trigonal system. This transition is coupled with severe reorganization of the polyiodide network from a layered to three-dimensional architecture. Quantum chemical calculations show that even though the two polymorphic phases are nearly isoenergetic under ambient conditions, the PV term is decisive in stabilizing the trigonal polymorph above the transition point. Phonon calculations using a non-local correlation functional that accounts for dispersion interactions confirm that this polymorph is dynamically unstable under ambient conditions. The high-pressure behavior of crystalline CsI3 can be correlated with other alkali metal trihalides, which undergo a similar sequence of structural changes upon load.
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.
Die Hochwasserkatastrophe im Juli 2021 in Westdeutschland erfordert eine kritische Diskussion über die Abschätzung der Hochwassergefährdung, Aktualisierung von Hochwassergefahrenkarten und Kommunikation von extremen Hochwasserszenarien. In der vorliegenden Arbeit wurde die Extremwertstatistik für die jährlichen maximalen Spitzenabflüsse am Pegel Altenahr im Ahrtal mit und ohne Berücksichtigung historischer Hochwasser berechnet und verglichen. Die Schätzung der Wiederkehrperiode für das aktuelle Hochwasser mittels Generalisierter Extremwertverteilung (GEV) unter Berücksichtigung historischer Hochwasser schwankt zwischen etwa 2.600 und über 58.700 Jahren (90%-Konfidenzintervall) mit einem Median bei etwa 8.600 Jahren, wogegen die Schätzung, die nur auf der systematisch gemessenen Abflusszeitreihe von 74 Jahren basiert, theoretisch eine Wiederkehrperiode von über 100 Millionen Jahren ergeben würde. Die Berücksichtigung der historischen Hochwasser führt zu einer dramatischen Änderung der Hochwasserquan-
tile, die für eine Gefahrenkartierung zugrunde gelegt werden. Die Anpassung der GEV an die Zeitreihe mit historischen Hochwassern zeigt dennoch, dass das GEV-Modell möglicherweise die Grundgesamtheit der Hochwasser im Ahrtal nicht adäquat abbilden kann. Es könnte sich im vorliegenden Fall um eine gemischte Stichprobe handeln, in der die extremen Hochwasser im Vergleich zu kleineren Ereignissen durch besondere Prozesse hervorgerufen werden. Somit könnten die Wahrscheinlichkeiten von extremen Hochwassern deutlich größer sein, als aus dem GEV-Modell hervorgeht. Hier sollte in Zukunft die Anwendung einer prozessbasierten Mischverteilung
untersucht werden. Der Vergleich von amtlichen Gefahrenkarten zu Extremhochwassern (HQextrem) im Ahrtal mit den Überflutungsflächen vom Juli 2021
zeigt eine deutliche Diskrepanz in den betroffenen Gebieten und die Notwendigkeit, die Grundlagen zur Erstellung der Extremszenarien zu überdenken. Die hydrodynamisch-numerischen Simulationen von 1.000-jährlichen Hochwassern (HQ1000) unter Berücksichtigung historischer Ereignisse und des größten historischen Hochwassers 1804 können die Gefährdung des Juli-Hochwassers 2021 deutlich besser widerspiegeln, wenngleich auch diese beiden Szenarien die Überflutungsflächen unterschätzen. Besondere Effekte wie die Verklausung von Brücken und die geomorphologischen Änderungen im Flussschlauch führten zu noch größeren Überflutungs- flächen im Juli 2021, als die Simulationsergebnisse zeigten. Basierend auf dieser Analyse wird eine einheitliche Festlegung von HQextrem bei Hochwassergefahrenkartierungen in Deutschland vorgeschlagen, die sich an höheren Hochwasserquantilen im Bereich von HQ1000 orientiert. Zusätzlich sollen simulationsbasierte Rekonstruktionen von den größten verlässlich dokumentierten historischen Hochwassern und/oder synthetische Worst-Case-Szenarien in den Hochwassergefahrenkarten gesondert dargestellt werden. Damit wird ein wichtiger Beitrag geleistet, um die potenziell betroffene Bevölkerung und das Katastrophenmanagement vor Überraschungen durch sehr seltene und extreme Hochwasser in Zukunft besser zu schützen.
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.
This study deals with the East Beni Suef Basin (Eastern Desert, Egypt) and aims to evaluate the source-generative potential, reconstruct the burial and thermal history, examine the most influential parameters on thermal maturity modeling, and improve on the models already published for the West Beni Suef to ultimately formulate a complete picture of the whole basin evolution.
Source rock evaluation was carried out based on TOC, Rock-Eval pyrolysis, and visual kerogen petrography analyses. Three kerogen types (II, II/III, and III) are distinguished in the East Beni Suef Basin, where the Abu Roash "F" Member acts as the main source rock with good to excellent source potential, oil-prone mainly type II kerogen, and immature to marginal maturity levels.
The burial history shows four depositional and erosional phases linked with the tectonic evolution of the basin. A hiatus (due to erosion or non-deposition) has occurred during the Late Eocene-Oligocene in the East Beni Suef Basin, while the West Beni Suef Basin has continued subsiding.
Sedimentation began later (Middle to Late Albian) with lower rates in the East Beni Suef Basin compared with the West Beni Suef Basin (Early Albian). The Abu Roash "F" source rock exists in the early oil window with a present-day transformation ratio of about 19% and 21% in the East and West Beni Suef Basin, respectively, while the Lower Kharita source rock, which is only recorded in the West Beni Suef Basin, has reached the late oil window with a present-day transformation ratio of about 70%.
The magnitude of erosion and heat flow have proportional and mutual effects on thermal maturity.
We present three possible scenarios of basin modeling in the East Beni Suef Basin concerning the erosion from the Apollonia and Dabaa formations.
Results of this work can serve as a basis for subsequent 2D and/or 3D basin modeling, which are highly recommended to further investigate the petroleum system evolution of the Beni Suef Basin.
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.
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 Gofa Province and the Chew Bahir Basin of southern Ethiopia constitute tectonically active regions, where the Southern Main Ethiopian Rift converges with the Northern Kenya Rift through a wide zone of extensional deformation with several north to northeast-trending, left-stepping en-e & PRIME;chelon basins. This sector of the Southern Main Ethiopian Rift is characterized by a semi-arid climate and a largely uniform lithology, and thus provides ideal conditions for studying the different parameters that define the tectonic and geomorphic features of this complex kinematic transfer zone. In this study, the degree of tectonic activity, spatiotemporal variations in extension, and the nature of kinematic linkage between different fault systems of the transfer zone are constrained by detailed quantitative geomorphic analysis of river catchments and focused field work. We analyzed fluvial and landscape morphometric characteristics in combination with structural, seismicity, and climatic data to better evaluate the tectono-geomorphic history of this transfer zone. Our data reveal significant north-south variations in the degree of extension from the Sawula Basin in the north (mature) to the Chew Bahir Basin in the south (juvenile). First, normalized channel-steepness indices and the spatial arrangement of knickpoints in footwall-draining streams suggest a gradual, southward shift in extensional deformation and recent tectonic activity. Second, based on 1-k(m) radius local relief and mean-hillslope maximum values that are consistent with ksn anomalies, we confirm strain localization within zones of fault interaction. Third, morphometric indices such as hypsometry, basin asymmetry factor, and valley floor width to valley height ratio also indicate a north to south gradient in tectonic activity, highlighting the importance of such a wide transfer zone with diffuse extension linking different rift segments during the break-up of continental crust.
Magmatic continental rifts often constitute nascent plate boundaries, yet long-term extension rates and transient rate changes associated with these early stages of continental breakup remain difficult to determine. Here, we derive a time-averaged minimum extension rate for the inner graben of the Northern Kenya Rift (NKR) of the East African Rift System for the last 0.5 m.y. We use the TanDEM-X science digital elevation model to evaluate fault-scarp geometries and determine fault throws across the volcano-tectonic axis of the inner graben of the NKR. Along rift-perpendicular profiles, amounts of cumulative extension are determined, and by integrating four new Ar-40/Ar-39 radiometric dates for the Silali volcano into the existing geochronology of the faulted volcanic units, time-averaged extension rates are calculated. This study reveals that in the inner graben of the NKR, the long-term extension rate based on mid-Pleistocene to recent brittle deformation has minimum values of 1.0-1.6 mm yr(-1), locally with values up to 2.0 mm yr(-1). A comparison with the decadal, geodetically determined extension rate reveals that at least 65% of the extension must be accommodated within a narrow, 20-km-wide zone of the inner rift. In light of virtually inactive border faults of the NKR, we show that extension is focused in the region of the active volcano-tectonic axis in the inner graben, thus highlighting the maturing of continental rifting in the NKR.
The 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.
The Devonian Las Chacras-Potrerillos batholith comprises six nested monzonitic to granitic intrusions with metaluminous to weakly peraluminous composition and a Sr-Nd isotopic signature indicating a dominantly juvenile mantle-derived source. The chemically most evolved units in the southern batholith contain a large number of intra-granitic, pod-shaped tourmaline-bearing pegmatites. This study uses in situ chemical and boron isotopic analyses of tourmaline from nine of these pegmatites to discuss their relationship to the respective host intrusions and the implications of their B-isotope composition for the source and evolution of the magmas. The tourmalines reveal a diversity in element composition (e.g., FeO, MgO, TiO2, CaO, MnO, F) which distinguishes individual pegmatites from one another. However, all have a narrow 5 11 B range of -13.7 to -10.5%0 (n = 100) which indicates a relatively uniform magmatic system and similar temperature conditions during tourmaline crystallization. The average delta(11) B value of -11.7%0 is typical for S-type granites and is within the range reported for peraluminous granites. pegmatites, and metamorphic units of the Ordovician basement into which the Las Chacras-Potrerillos batholith intruded. The B-isotope evidence argues for a crustal boron source like that of the Ordovician basement, in contrast to the metaluminous to weakly peraluminous composition and juvenile initial Sr and Nd isotope ratios of the Las Chacras-Potrerillos batholith magmas. We propose that the boron was not derived from the magma source region but was incorporated from dehydration melting of elastic metasedimentary rocks higher up in the crustal column.
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.
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.
Beyond CO2 equivalence
(2022)
In this article we review the physical and chemical properties of methane (CH4) relevant to impacts on climate, ecosystems, and air pollution, and examine the extent to which this is reflected in climate and air pollution governance. Although CH4 is governed under the UNFCCC climate regime, its treatment there is limited to the ways in which it acts as a "CO2 equivalent" climate forcer on a 100-year time frame. The UNFCCC framework neglects the impacts that CH4 has on near-term climate, as well its impacts on human health and ecosystems, which are primarily mediated by methane's role as a precursor to tropospheric ozone. Frameworks for air quality governance generally address tropospheric ozone as a pollutant, but do not regulate CH4 itself. Methane's climate and air quality impacts, together with its alarming rise in atmospheric concentrations in recent years, make it clear that mitigation of CH4 emissions needs to be accelerated globally. We examine challenges and opportunities for further progress on CH4 mitigation within the international governance landscapes for climate change and air pollution.
Deep hydrothermal Mo, W, and base metal mineralization at the Sweet Home mine (Detroit City portal) formed in response to magmatic activity during the Oligocene. Microthermometric data of fluid inclusions trapped in greisen quartz and fluorite suggest that the early-stage mineralization at the Sweet Home mine precipitated from low- to medium-salinity (1.5-11.5 wt% equiv. NaCl), CO2-bearing fluids at temperatures between 360 and 415 degrees C and at depths of at least 3.5 km. Stable isotope and noble gas isotope data indicate that greisen formation and base metal mineralization at the Sweet Home mine was related to fluids of different origins. Early magmatic fluids were the principal source for mantle-derived volatiles (CO2, H2S/SO2, noble gases), which subsequently mixed with significant amounts of heated meteoric water. Mixing of magmatic fluids with meteoric water is constrained by delta H-2(w)-delta O-18(w) relationships of fluid inclusions. The deep hydrothermal mineralization at the Sweet Home mine shows features similar to deep hydrothermal vein mineralization at Climax-type Mo deposits or on their periphery. This suggests that fluid migration and the deposition of ore and gangue minerals in the Sweet Home mine was triggered by a deep-seated magmatic intrusion. The findings of this study are in good agreement with the results of previous fluid inclusion studies of the mineralization of the Sweet Home mine and from Climax-type Mo porphyry deposits in the Colorado Mineral Belt.
Continental rifting is responsible for the generation of major sedimentary basins, both during rift inception and during the formation of rifted continental margins. Geophysical and field studies revealed that rifts feature complex networks of normal faults but the factors controlling fault network properties and their evolution are still matter of debate. Here, we employ high-resolution 2D geodynamic models (ASPECT) including two-way coupling to a surface processes (SP) code (FastScape) to conduct 12 models of major rift types that are exposed to various degrees of erosion and sedimentation. We further present a novel quantitative fault analysis toolbox (Fatbox), which allows us to isolate fault growth patterns, the number of faults, and their length and displacement throughout rift history. Our analysis reveals that rift fault networks may evolve through five major phases: (a) distributed deformation and coalescence, (b) fault system growth, (c) fault system decline and basinward localization, (d) rift migration, and (e) breakup. These phases can be correlated to distinct rifted margin domains. Models of asymmetric rifting suggest rift migration is facilitated through both ductile and brittle deformation within a weak exhumation channel that rotates subhorizontally and remains active at low angles. In sedimentation-starved settings, this channel satisfies the conditions for serpentinization. We find that SP are not only able to enhance strain localization and to increase fault longevity but that they also reduce the total length of the fault system, prolong rift phases and delay continental breakup.
Wetlands are dynamic ecosystems that require continuous monitoring and assessment of degradation status to design strategies for their sustainable management. While hydrology provides the primary functional control for the wetland ecosystem, the loss of landscape connectivity influences wetland degradation in a major way as it leads to fragmentation. This article aims to integrate hydrogeomorphic and ecological concepts for the assessment of degradation status and its causal factors for a large wetland in the western Ganga plains, India, the Haiderpur, using a wetlandscape approach. We have used a remote-sensing-based approach, which offers a powerful tool for assessing and linking cross-scale structures, functions, and controls in a wetlandscape. The Haiderpur, a Ramsar site since December 2021, is an artificial wetland located on the right bank of the Ganga River wherein the inflows are controlled by a barrage constructed on the Ganga River apart from smaller tributaries flowing in from the north. A novel aspect of this work is the integration of river dynamics and its connectivity to the wetlandscape to understand the spatiotemporal variability in the waterspread area in the wetland. In this work, we have developed an integrated wetlandscape assessment approach by evaluating wetland's geomorphic and hydrological connectivity status for the period 1993-2019 (25 years) across three different spatial scales - regional, catchment, and wetland. We have highlighted the ecological implications of connectivity and patch dynamics for developing sustainable wetland management plans.
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.
Landslides
(2022)
Erosion by landslides is a common phenomenon in mountain regions around the globe, affecting all climatic zones. Landslides facilitate bedrock weathering, pedogenesis and ecological succession, being key drivers of biodiversity. Landslide chronosequences have long been used for studies of vegetation succession in initial ecosystems, but they further offer ideal model systems for studies of soil development and microbial community succession. In this review we synthesize the state of knowledge on the role of landslides in ecosystems, their influence on element cycles and interactions with biota. Further, we discuss feedback mechanisms between global warming, landslide activity and greenhouse gas emissions. In the view of increasing anthropogenic influence and climate change, soils are becoming a critical resource. Due to their ubiquity, landslide chronosequences have the potential to provide critical insights into soil development under different climates and thereby contribute to future soil restoration efforts.
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.
The stabilizing properties of mineral-organic carbon (OC) interactions have been studied in many soil environments (temperate soils, podzol lateritic soils, and paddy soils). Recently, interest in their role in permafrost regions is increasing as permafrost was identified as a hotspot of change. In thawing ice-rich permafrost regions, such as the Yedoma domain, 327-466 Gt of frozen OC is buried in deep sediments. Interactions between minerals and OC are important because OC is located very near the mineral matrix. Mineral surfaces and elements could mitigate recent and future greenhouse gas emissions through physical and/or physicochemical protection of OC. The dynamic changes in redox and pH conditions associated with thermokarst lake formation and drainage trigger metal-oxide dissolution and precipitation, likely influencing OC stabilization and microbial mineralization. However, the influence of thermokarst processes on mineral-OC interactions remains poorly constrained. In this study, we aim to characterize Fe, Mn, Al, and Ca minerals and their potential protective role for OC. Total and selective extractions were used to assess the crystalline and amorphous oxides or complexed metal pools as well as the organic acids found within these pools. We analyzed four sediment cores from an ice-rich permafrost area in Central Yakutia, which were drilled (i) in undisturbed Yedoma uplands, (ii) beneath a recent lake formed within Yedoma deposits, (iii) in a drained thermokarst lake basin, and (iv) beneath a mature thermokarst lake from the early Holocene period. We find a decrease in the amount of reactive Fe, Mn, Al, and Ca in the deposits on lake formation (promoting reduction reactions), and this was largely balanced by an increase in the amount of reactive metals in the deposits on lake drainage (promoting oxidation reactions). We demonstrate an increase in the metal to C molar ratio on thermokarst process, which may indicate an increase in metal-C bindings and could provide a higher protective role against microbial mineralization of organic matter. Finally, we find that an increase in mineral-OC interactions corresponded to a decrease in CO2 and CH4 gas emissions on thermokarst process. Mineral-OC interactions could mitigate greenhouse gas production from permafrost thaw as soon as lake drainage occurs.
Tropical Lake Sentani in the Indonesian Province Papua consists of four separate basins and is surrounded by a catchment with a very diverse geology. We characterized the surface sediment (upper 5 cm) of the lake's four sub-basins based on multivariate statistical analyses (principal component analysis, hierarchical clustering) of major element compositions obtained by X-ray fluorescence scanning. Three types of sediment are identified based on distinct compositional differences between rivers, shallow/proximal and deep/distal lake sediments. The different sediment types are mainly characterized by the correlation of elements associated with redox processes (S, Mn, Fe), carbonates (Ca), and detrital input (Ti, Al, Si, K) derived by river discharge. The relatively coarse-grained river sediments mainly derive form the mafic catchment geology and contribution of the limestone catchment geology is only limited. Correlation of redox sensitive and detrital elements are used to reveal oxidation conditions, and indicate oxic conditions in river samples and reducing conditions for lake sediments. Organic carbon (TOC) generally correlates with redox sensitive elements, although a correlation between TOC and individual elements change strongly between the three sediment types. Pyrite is the quantitatively dominant reduced sulfur mineral, monosulfides only reach appreciable concentrations in samples from rivers draining mafic and ultramafic catchments. Our study shows large spatial heterogeneity within the lake's sub-basins that is mainly caused by catchment geology and topography, river runoff as well as the bathymetry and the depth of the oxycline. We show that knowledge about lateral heterogeneity is crucial for understanding the geochemical and sedimentological variations recorded by these sediments. The highly variable conditions make Lake Sentani a natural laboratory, with its different sub-basins representing different depositional environments under identical tropical climate conditions.
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.
In the past decade, sediment connectivity has become a widely recognized characteristic of a geomorphic system. However, the quantification of functional connectivity (i.e. connectivity which arises due to the actual occurrence of sediment transport processes) and its variation over space and time is still a challenge. In this context, this study assesses the effects of expected future phenomena in the context of climate change (i.e. glacier retreat, permafrost degradation or meteorological extreme events) on sediment transport dynamics in a glacierised Alpine basin. The study area is the Sulden river basin (drainage area 130 km(2)) in the Italian Alps, which is composed of two geomorphologically diverse sub-basins. Based on graph theory, we evaluated the spatio-temporal variations in functional connectivity in these two sub-basins. The graph-object, obtained by manually mapping sediment transport processes between landforms, was adapted to 6 different hydro-meteorological scenarios, which derive from combining base, heatwave and rainstorm conditions with snowmelt and glacier-melt periods. For each scenario and each sub-basin, the sediment transport network and related catchment characteristics were analysed. To compare the effects of the scenarios on functional connectivity, we introduced a connectivity degree, calculated based on the area of the landforms involved in sediment cascades. Results indicate that the area of the basin connected to its outlet in terms of sediment transport might feature a six-fold increase in case of rainstorm conditions compared to "average " meteorological conditions assumed for the base scenario. Furthermore, markedly different effects of climate change on sediment connectivity are expected between the two sub-catchments due to their contrasting morphological and lithological characteristics, in terms of relative importance of rainfall triggered colluvial processes vs temperature-driven proglacial fluvial dynamics.
Reducing greenhouse gas emissions in food systems is becoming more challenging as food is increasingly consumed away from producer regions, highlighting the need to consider emissions embodied in trade in agricultural emissions accounting.
To address this, our study explores recent trends in trade-adjusted agricultural emissions of food items at the global, regional, and national levels.
We find that emissions are largely dependent on a country’s consumption patterns and their agricultural emission intensities relative to their trading partners’.
The absolute differences between the production-based and trade-adjusted emissions accounting approaches are especially apparent for major agricultural exporters and importers and where large shares of emission-intensive items such as ruminant meat, milk products and rice are involved.
In relative terms, some low-income and emerging and developing economies with consumption of high emission intensity food products show large differences between approaches.
Similar trends are also found under various specifications that account for trade and re-exports differently.
These findings could serve as an important element towards constructing national emissions reduction targets that consider trading partners, leading to more effective emissions reductions overall.
The concepts of CO2 emission, global warming, climate change, and their environmental impacts are of utmost importance for the understanding and protection of the ecosystems.
Among the natural sources of gases into the atmosphere, the contribution of geogenic sources plays a crucial role. However, while subaerial emissions are widely studied, submarine outgassing is not yet well understood.
In this study, we review and catalog 122 literature and unpublished data of submarine emissions distributed in ten coastal areas of the Aegean Sea. This catalog includes descriptions of the degassing vents through in situ observations, their chemical and isotopic compositions, and flux estimations.
Temperatures and pH data of surface seawaters in four areas affected by submarine degassing are also presented.
This overview provides useful information to researchers studying the impact of enhanced seawater CO2 concentrations related either to increasing CO2 levels in the atmosphere or leaking carbon capture and storage systems.
In-depth understanding of the reorganization of the hydrological cycle in response to global climate change is crucial in highly sensitive regions like the eastern Mediterranean, where water availability is a major factor for socioeconomic and political development.
The sediments of Lake Lisan provide a unique record of hydroclimatic change during the last glacial to Holocene transition (ca. 24-11 ka) with its tremendous water level drop of similar to 240 m that finally led to its transition into the present hypersaline water body-the Dead Sea.
Here we utilize high-resolution sedimentological analyses from the marginal terraces and deep lake to reconstruct an unprecedented seasonal record of the last millennia of Lake Lisan. Aragonite varve formation in intercalated intervals of our record demonstrates that a stepwise long-term lake level decline was interrupted by almost one millennium of rising or stable water level.
Even periods of pronounced water level drops indicated by gypsum deposition were interrupted by decades of positive water budgets.
Our results thus highlight that even during major climate change at the end of the last glacial, decadal to millennial periods of relatively stable or positive moisture supply occurred which could have been an important premise for human sedentism.
The region of the Andean back-arc of northwestern Argentina has been struck by several magnitude >= 6 crustal earthquakes since the first historically recorded event in 1692. One of these events corresponds to the Anta earthquake on 25 August 1948, with epicenter in the Santa Barbara System causing three deaths and severe damage in Salta and Jujuy provinces with maximum Modified Mercalli seismic intensities (MMI) of IX. We collected and digitized analog seismograms of this earthquake from worldwide seismic observatories in order to perform first-motion analysis and modeling of long-period teleseismic P-waveforms. Our results indicate a simple seismic source of M0 = 2.85 x 1019 N m consistent with a moment magnitude Mw = 6.9. We have also tested for the focal depth determining a shallow source at 8 km with a reverse focal mechanism solution with a minor dextral strike-slip component (strike 20 degrees, dip 30 degrees, rake 120 degrees) from the best fit of waveforms. Using magnitude size empirical relationships, the comparison of the obtained Mw 6.9 magnitude value and the ca. 10,000 km2 area of MMI >= IX from our seismic intensity map, which was obtained from newspaper and many historical reports, indicates a rupture length of 42 +/- 8 km for the Anta earthquake. We show our results in a 3D geological model around the epicentral area, which integrates modern seismicity, geological data, and information of a previously studied east-west cross section located a few kilometers south of the 1948 epicenter. The integration of all available information provides evidence of the re-activation of the Pie de la Sierra del Gallo fault during the 1948 Mw 6.9 shallow earthquake; this thrust fault bounds the Santa Barbara System along its western foothill.
A novel idea for an optimal time delay state space reconstruction from uni- and multivariate time series is presented. The entire embedding process is considered as a game, in which each move corresponds to an embedding cycle and is subject to an evaluation through an objective function. This way the embedding procedure can be modeled as a tree, in which each leaf holds a specific value of the objective function. By using a Monte Carlo ansatz, the proposed algorithm populates the tree with many leafs by computing different possible embedding paths and the final embedding is chosen as that particular path, which ends at the leaf with the lowest achieved value of the objective function. The method aims to prevent getting stuck in a local minimum of the objective function and can be used in a modular way, enabling practitioners to choose a statistic for possible delays in each embedding cycle as well as a suitable objective function themselves. The proposed method guarantees the optimization of the chosen objective function over the parameter space of the delay embedding as long as the tree is sampled sufficiently. As a proof of concept, we demonstrate the superiority of the proposed method over the classical time delay embedding methods using a variety of application examples. We compare recurrence plot-based statistics inferred from reconstructions of a Lorenz-96 system and highlight an improved forecast accuracy for map-like model data as well as for palaeoclimate isotope time series. Finally, we utilize state space reconstruction for the detection of causality and its strength between observables of a gas turbine type thermoacoustic combustor.
A large landslide (frozen debris avalanche) occurred at Assapaat on the south coast of the Nuussuaq Peninsula in Central West Greenland on June 13, 2021, at 04:04 local time. We present a compilation of available data from field observations, photos, remote sensing, and seismic monitoring to describe the event. Analysis of these data in combination with an analysis of pre- and post-failure digital elevation models results in the first description of this type of landslide. The frozen debris avalanche initiated as a 6.9 * 10(6) m(3) failure of permafrozen talus slope and underlying colluvium and till at 600-880 m elevation. It entrained a large volume of permafrozen colluvium along its 2.4 km path in two subsequent entrainment phases accumulating a total volume between 18.3 * 10(6) and 25.9 * 10(6) m(3). About 3.9 * 10(6) m(3) is estimated to have entered the Vaigat strait; however, no tsunami was reported, or is evident in the field. This is probably because the second stage of entrainment along with a flattening of slope angle reduced the mobility of the frozen debris avalanche. We hypothesise that the initial talus slope failure is dynamically conditioned by warming of the ice matrix that binds the permafrozen talus slope. When the slope ice temperature rises to a critical level, its shear resistance is reduced, resulting in an unstable talus slope prone to failure. Likewise, we attribute the large-scale entrainment to increasing slope temperature and take the frozen debris avalanche as a strong sign that the permafrost in this region is increasingly at a critical state. Global warming is enhanced in the Arctic and frequent landslide events in the past decade in Western Greenland let us hypothesise that continued warming will lead to an increase in the frequency and magnitude of these types of landslides. Essential data for critical arctic slopes such as precipitation, snowmelt, and ground and surface temperature are still missing to further test this hypothesis. It is thus strongly required that research funds are made available to better predict the change of landslide threat in the Arctic.
Measuring the variability of incoming neutrons locally would be usefull for the cosmic-ray neutron sensing (CRNS) method. As the measurement of high energy neutrons is not so easy, alternative particles can be considered for such purpose. Among them, muons are particles created from the same cascade of primary cosmic-ray fluxes that generate neutrons at the ground. In addition, they can be easily detected by small and relatively inexpensive detectors. For these reasons they could provide a suitable local alternative to incoming corrections based on remote neutron monitor data. The reported measurements demonstrated that muon detection system can detect incoming cosmic-ray variations locally. Furthermore the precision of this measurement technique is considered adequate for many CRNS applications.
In this study, 3-D models of P-wave velocity (Vp) and P-wave and S-wave ratio (Vp/Vs) of the crust and upper mantle in the Eastern and eastern Southern Alps (northern Italy and southern Austria) were calculated using local earthquake tomography (LET). The data set includes high-quality arrival times from well-constrained hypocenters observed by the dense, temporary seismic networks of the AlpArray AASN and SWATH-D. The resolution of the LET was checked by synthetic tests and analysis of the model resolution matrix. The small inter-station spacing (average of similar to 15 km within the SWATH-D network) allowed us to image crustal structure at unprecedented resolution across a key part of the Alps. The derived P velocity model revealed a highly heterogeneous crustal structure in the target area. One of the main findings is that the lower crust is thickened, forming a bulge at 30-50 km depth just south of and beneath the Periadriatic Fault and the Tauern Window. This indicates that the lower crust decoupled both from its mantle substratum as well as from its upper crust. The Moho, taken to be the iso-velocity contour of Vp = 7.25 km/s, agrees with the Moho depth from previous studies in the European and Adriatic forelands. It is shallower on the Adriatic side than on the European side. This is interpreted to indicate that the European Plate is subducted beneath the Adriatic Plate in the Eastern and eastern Southern Alps.
Climate Benefits of Cleaner Energy Transitions in East and South Asia Through Black Carbon Reduction
(2022)
The state of air pollution has historically been tightly linked to how we produce and use energy. Air pollutant emissions over Asia are now changing rapidly due to cleaner energy transitions; however, magnitudes of benefits for climate and air quality remain poorly quantified. The associated risks involve adverse health impacts, reduced agricultural yields, reduced freshwater availability, contributions to climate change, and economic costs. We focus particularly on climate benefits of energy transitions by making first-time use of two decades of high quality observations of atmospheric loading of light-absorbing black carbon (BC) over Kanpur (South Asia) and Beijing (East Asia) and relating these observations to changing energy, emissions, and economic trends in India and China. Our analysis reveals that absorption aerosol optical depth (AAOD) due to BC has decreased substantially, by 40% over Kanpur and 60% over Beijing between 2001 and 2017, and thus became decoupled from regional economic growth. Furthermore, the resultant decrease in BC emissions and BC AAOD over Asia is regionally coherent and occurs primarily due to transitions into cleaner energies (both renewables and fossil fuels) and not due to the decrease in primary energy supply or decrease in use of fossil use and biofuels and waste. Model simulations show that BC aerosols alone contribute about half of the surface temperature change (warming) of the total forcing due to greenhouse gases, natural and internal variability, and aerosols, thus clearly revealing the climate benefits due to a reduction in BC emissions, which would significantly reduce global warming. However, this modeling study excludes responses from natural variability, circulation, and sea ice responses, which cause relatively strong temperature fluctuations that may mask signals from BC aerosols. Our findings show additional benefits for climate (beyond benefits of CO2 reduction) and for several other issues of sustainability over South and East Asia, provide motivation for ongoing cleaner energy production, and consumption transitions, especially when they are associated with reduced emissions of air pollutants. Such an analysis connecting the trends in energy transitions and aerosol absorption loading, unavailable so far, is crucial for simulating the aerosol climate impacts over Asia which is quite uncertain.
In this study we present a novel method for the automatic detection of minerals and elements using hyperspectral transmittance imaging microscopy measurements of complete thin sections (HyperTIM).
This is accomplished by using a hyperspectral camera system that operates in the visible and near-infrared (VNIR) range with a specifically designed sample holder, scanning setup, and a microscope lens.
We utilize this method on a monazite ore thin section from Steenkampskraal (South Africa), which we analyzed for the rare earth element (REE)-bearing mineral monazite ((Ce,Nd,La)PO4), with high concentrations of Nd. The transmittance analyses with the hyperspectral VNIR camera can be used to identify REE minerals and Nd in thin sections.
We propose a three-point band depth index, the Nd feature depth index (NdFD), and its related product the Nd band depth index (NdBDI), which enables automatic mineral detection and classification for the Nd-bearing monazites in thin sections. In combination with the average concentration of the relative Nd content, it permits a destruction-free, total concentration calculation for Nd across the entire thin section.
A novel approach for estimating precipitation patterns is developed here and applied to generate a new hydrologically corrected daily precipitation dataset, called RAIN4PE (Rain for Peru and Ecuador), at 0.1 degrees spatial resolution for the period 1981-2015 covering Peru and Ecuador. It is based on the application of 1) the random forest method to merge multisource precipitation estimates (gauge, satellite, and reanalysis) with terrain elevation, and 2) observed and modeled streamflow data to first detect biases and second further adjust gridded precipitation by inversely applying the simulated results of the ecohydrological model SWAT (Soil and Water Assessment Tool). Hydrological results using RAIN4PE as input for the Peruvian and Ecuadorian catchments were compared against the ones when feeding other uncorrected (CHIRP and ERA5) and gauge-corrected (CHIRPS, MSWEP, and PISCO) precipitation datasets into the model. For that, SWAT was calibrated and validated at 72 river sections for each dataset using a range of performance metrics, including hydrograph goodness of fit and flow duration curve signatures. Results showed that gauge-corrected precipitation datasets outperformed uncorrected ones for streamflow simulation. However, CHIRPS, MSWEP, and PISCO showed limitations for streamflow simulation in several catchments draining into the Pacific Ocean and the Amazon River. RAIN4PE provided the best overall performance for streamflow simulation, including flow variability (low, high, and peak flows) and water budget closure. The overall good performance of RAIN4PE as input for hydrological modeling provides a valuable criterion of its applicability for robust countrywide hydrometeorological applications, including hydroclimatic extremes such as droughts and floods. Significance StatementWe developed a novel precipitation dataset RAIN4PE for Peru and Ecuador by merging multisource precipitation data (satellite, reanalysis, and ground-based precipitation) with terrain elevation using the random forest method. Furthermore, RAIN4PE was hydrologically corrected using streamflow data in watersheds with precipitation underestimation through reverse hydrology. The results of a comprehensive hydrological evaluation showed that RAIN4PE outperformed state-of-the-art precipitation datasets such as CHIRP, ERA5, CHIRPS, MSWEP, and PISCO in terms of daily and monthly streamflow simulations, including extremely low and high flows in almost all Peruvian and Ecuadorian catchments. This underlines the suitability of RAIN4PE for hydrometeorological applications in this region. Furthermore, our approach for the generation of RAIN4PE can be used in other data-scarce regions.
Assessment of climate change impact on discharge of the lakhmass catchment (Northwest Tunisia)
(2022)
The Mediterranean region is increasingly recognized as a climate change hotspot but is highly underrepresented in hydrological climate change studies. This study aims to investigate the climate change effects on the hydrology of Lakhmass catchment in Tunisia. Lakhmass catchment is a part of the Medium Valley of Medjerda in northwestern Tunisia that drains an area of 126 km(2). First, the Hydrologiska Byrans Vattenbalansavdelning light (HBV-light) model was calibrated and validated successfully at a daily time step to simulate discharge during the 1981-1986 period. The Nash Sutcliffe Efficiency and Percent bias (NSE, PBIAS) were (0.80, +2.0%) and (0.53, -9.5%) for calibration (September 1982-August 1984) and validation (September 1984-August 1986) periods, respectively. Second, HBV-light model was considered as a predictive tool to simulate discharge in a baseline period (1981-2009) and future projections using data (precipitation and temperature) from thirteen combinations of General Circulation Models (GCMs) and Regional Climatic Models (RCMs). We used two trajectories of Representative Concentration Pathways, RCP4.5 and RCP8.5, suggested by the Intergovernmental Panel on Climate Change (IPCC). Each RCP is divided into three projection periods: near-term (2010-2039), mid-term (2040-2069) and long-term (2070-2099). For both scenarios, a decrease in precipitation and discharge will be expected with an increase in air temperature and a reduction in precipitation with almost 5% for every +1 degrees C of global warming. By long-term (2070-2099) projection period, results suggested an increase in temperature with about 2.7 degrees C and 4 degrees C, and a decrease in precipitation of approximately 7.5% and 15% under RCP4.5 and RCP8.5, respectively. This will likely result in a reduction of discharge of 12.5% and 36.6% under RCP4.5 and RCP8.5, respectively. This situation calls for early climate change adaptation measures under a participatory approach, including multiple stakeholders and water users.
Understanding the influence of climate change and population pressure on human conflict remains a critically important topic in the social sciences. Long-term records that evaluate these dynamics across multiple centuries and outside the range of modern climatic variation are especially capable of elucidating the relative effect of-and the interaction between-climate and demography. This is crucial given that climate change may structure population growth and carrying capacity, while both climate and population influence per capita resource availability. This study couples paleoclimatic and demographic data with osteological evaluations of lethal trauma from 149 directly accelerator mass spectrometry C-14-dated individuals from the Nasca highland region of Peru. Multiple local and supraregional precipitation proxies are combined with a summed probability distribution of 149 C-14 dates to estimate population dynamics during a 700-y study window. Counter to previous findings, our analysis reveals a precipitous increase in violent deaths associated with a period of productive and stable climate, but volatile population dynamics. We conclude that favorable local climate conditions fostered population growth that put pressure on the marginal and highly circumscribed resource base, resulting in violent resource competition that manifested in over 450 y of internecine warfare. These findings help support a general theory of intergroup violence, indicating that relative resource scarcity-whether driven by reduced resource abundance or increased competition-can lead to violence in subsistence societies when the outcome is lower per capita resource availability.
Efforts have been made in the past to enhance building exposure models on a regional scale with increasing spatial resolutions by integrating different data sources. This work follows a similar path and focuses on the downscaling of the existing SARA exposure model that was proposed for the residential building stock of the communes of Valparaiso and Vina del Mar (Chile). Although this model allowed great progress in harmonising building classes and characterising their differential physical vulnerabilities, it is now outdated, and in any case, it is spatially aggregated over large administrative units. Hence, to more accurately consider the impact of future earthquakes on these cities, it is necessary to employ more reliable exposure models. For such a purpose, we propose updating this existing model through a Bayesian approach by integrating ancillary data that has been made increasingly available from Volunteering Geo-Information (VGI) activities. Its spatial representation is also optimised in higher resolution aggregation units that avoid the inconvenience of having incomplete building-by-building footprints. A worst-case earthquake scenario is presented to calculate direct economic losses and highlight the degree of uncertainty imposed by exposure models in comparison with other parameters used to generate the seismic ground motions within a sensitivity analysis. This example study shows the great potential of using increasingly available VGI to update worldwide building exposure models as well as its importance in scenario-based seismic risk assessment.
The release of greenhouse gases from the large organic carbon stock in permafrost deposits in the circumarctic regions may accelerate global warming upon thaw. The extent of this positive climate feedback is thought to be largely controlled by the microbial degradability of the organic matter preserved in these sediments. In addition, weathering and oxidation processes may release inorganic carbon preserved in permafrost sediments as CO2, which is generally not accounted for. We used C-13 and C-14 analysis and isotopic mass balances to differentiate and quantify organic and inorganic carbon released as CO2 in the field from an active retrogressive thaw slump of Pleistocene-age Yedoma and during a 1.5-years incubation experiment. The results reveal that the dominant source of the CO2 released from freshly thawed Yedoma exposed as thaw mound is Pleistocene-age organic matter (48-80%) and to a lesser extent modern organic substrate (3-34%). A significant portion of the CO2 originated from inorganic carbon in the Yedoma (17-26%). The mixing of young, active layer material with Yedoma at a site on the slump floor led to the preferential mineralization of this young organic carbon source. Admixtures of younger organic substrates in the Yedoma thaw mound were small and thus rapidly consumed as shown by lower contributions to the CO2 produced during few weeks of aerobic incubation at 4 degrees C corresponding to approximately one thaw season. Future CO2 fluxes from the freshly thawed Yedoma will contain higher proportions of ancient inorganic (22%) and organic carbon (61-78%) as suggested by the results at the end, after 1.5 years of incubation. The increasing contribution of inorganic carbon during the incubation is favored by the accumulation of organic acids from microbial organic matter degradation resulting in lower pH values and, in consequence, in inorganic carbon dissolution. Because part of the inorganic carbon pool is assumed to be of pedogenic origin, these emissions would ultimately not alter carbon budgets. The results of this study highlight the preferential degradation of younger organic substrates in freshly thawed Yedoma, if available, and a substantial release of CO2 from inorganic sources.
Groundwater is critical in supporting current and future reliable water supply throughout Africa. Although continental maps of groundwater storage and recharge have been developed, we currently lack a clear understanding on how the controls on groundwater recharge vary across the entire continent. Reviewing the existing literature, we synthesize information on reported groundwater recharge controls in Africa. We find that 15 out of 22 of these controls can be characterised using global datasets. We develop 11 descriptors of climatic, topographic, vegetation, soil and geologic properties using global datasets, to characterise groundwater recharge controls in Africa. These descriptors cluster Africa into 15 Recharge Landscape Units for which we expect recharge controls to be similar. Over 80% of the continents land area is organized by just nine of these units. We also find that aggregating the Units by similarity into four broader Recharge Landscapes (Desert, Dryland, Wet tropical and Wet tropical forest) provides a suitable level of landscape organisation to explain differences in ground-based long-term mean annual recharge and recharge ratio (annual recharge / annual precipitation) estimates. Furthermore, wetter Recharge Landscapes are more efficient in converting rainfall to recharge than drier Recharge Landscapes as well as having higher annual recharge rates. In Dryland Recharge Landscapes, we found that annual recharge rates largely varied according to mean annual precipitation, whereas recharge ratio estimates increase with increasing monthly variability in P-PET. However, we were unable to explain why ground based estimates of recharge signatures vary across other Recharge Landscapes, in which there are fewer ground based recharge estimates, using global datasets alone. Even in dryland regions, there is still considerable unexplained variability in the estimates of annual recharge and recharge ratio, stressing the limitations of global datasets for investigating ground-based information.
Despite its high-seismogenic potential, the details of the seismogenic processes of Zagros Simply Folded Belt (SFB) remains debated. Three large earthquakes (M-w 7.3, 5.9 and 6.3) struck in the Lurestan arc of the Zagros SFB in 2017 and 2018. The sequence was recorded by seismic stations at regional, and teleseismic distances. Coseismic surface displacements, measured by Sentinel-1A/B satellites, provide additional data and a unique opportunity to study these earthquakes in detail. Here, we complement previous studies of the coseismic slip distribution of the 12 November 2017 M-w 7.3 Ezgeleh earthquake by a detailed analysis of its aftershocks, and we analysed the rupture process of the two interrelated earthquakes (25 August 2018 M-w 5.9 Tazehabad and the 25 November 2018 M-w 6.3 Sarpol-e Zahab earthquakes). We model the surface displacements obtained from Interferometric Synthetic Aperture Radar (InSAR) measurements and seismic records. We conduct non-linear probabilistic optimizations based on joint InSAR and seismic data to obtain finite-fault rupture of these earthquakes. The Lurestan arc earthquakes were followed by a sustained aftershock activity, with 133 aftershocks exceeding M-n 4.0 until 30 December 2019. We rely on the permanent seismic networks of Iran and Iraq to relocate similar to 700 M-n 3 + events and estimate moment tensor solutions for 85 aftershocks down to M-w 4.0. The 2017 Ezgeleh earthquake has been considered to activate a low-angle (similar to 17 degrees) dextral-thrust fault at the depth of 10-20 km. However, most of its aftershocks have shallow centroid depths (8-12 km). The joint interpretation of finite source models, moment tensor and hypocentral location indicate that the 2018 Tazehabad and Sarpol-e Zahab earthquakes ruptured different strike-slip structures, providing evidence for the activation of the sinistral and dextral strike-slip faults, respectively. The deformation in the Lurestan arc is seismically accommodated by a complex fault system involving both thrust and strike-slip faults. Knowledge about the deformation characteristics is important for the understanding of crustal shortening, faulting and hazard and risk assessment in this region.
The Altiplano-Puna plateau, in Central Andes, is the second-largest continental plateau on Earth, extending between 22 degrees and 27 degrees S at an average altitude of 4400 m. The Puna plateau has been formed in consequence of the subduction of the oceanic Nazca Plate beneath the continental South American plate, which has an average crustal thickness of 50 km at this location. A large seismicity cluster, the Jujuy cluster, is observed at depth of 150-250 km beneath the central region of the Puna plateau. The cluster is seismically very active, with hundreds of earthquakes reported and a peak magnitude MW 6.6 on 25th August 2006. The cluster is situated in one of three band of intermediate-depth focus seismicity, which extend parallel to the trench roughly North to South. It has been hypothesized that the Jujuy cluster could be a seismic nest, a compact seismogenic region characterized by a high stationary activity relative to its surroundings. In this study, we collected more than 40 years of data from different catalogs and proof that the cluster meets the three conditions of a seismic nest. Compared to other known intermediate depth nests at Hindu Kush (Afganisthan) or Bucaramanga (Colombia), the Jujuy nest presents an outstanding seismicity rate, with more than 100 M4+ earthquakes per year. We additionally performed a detailed analysis of the rupture process of some of the largest earthquakes in the nest, by means of moment tensor inversion and directivity analysis. We focused on the time period 2017-2018, where the seismic monitoring was the most extended. Our results show that earthquakes in the nest take place within the eastward subducting oceanic plate, but rupture along sub-horizontal planes dipping westward. We suggest that seismicity at Jujuy nest is controlled by dehydration processes, which are also responsible for the generation of fluids ascending to the crust beneath the Puna volcanic region. We use the rupture plane and nest geometry to provide a constraint to maximal expected magnitude, which we estimate as MW -6.7.
Ground-motion models (GMMs) are often used to predict the random distribution of Spectral accelerations (SAs) at a site due to a nearby earthquake. In probabilistic seismic hazard and risk assessment, large earthquakes occurring close to a site are considered as critical scenarios. GMMs are expected to predict realistic SAs with low within-model uncertainty (sigma(mu)) for such rare scenarios. However, the datasets used to regress GMMs are usually deficient of data from critical scenarios. The (Kotha et al., A Regionally Adaptable Ground-Motion Model for Shallow Crustal Earthquakes in Europe Bulletin of Earthquake Engineering 18:4091-4125, 2020) GMM developed from the Engineering strong motion (ESM) dataset was found to predict decreasing short-period SAs with increasing M-W >= M-h = 6.2, and with large sigma(mu) at near-source distances <= 30km. In this study, we updated the parametrisation of the GMM based on analyses of ESM and the Near source strong motion (NESS) datasets. With M-h = 5.7, we could rectify the M-W scaling issue, while also reducing sigma(mu). at M-W >= M-h. We then evaluated the GMM against NESS data, and found that the SAs from a few large, thrust-faulting events in California, New Zealand, Japan, and Mexico are significantly higher than GMM median predictions. However, recordings from these events were mostly made on soft-soil geology, and contain anisotropic pulse-like effects. A more thorough non-ergodic treatment of NESS was not possible because most sites sampled unique events in very diverse tectonic environments. We provide an updated set of GMM coefficients,sigma(mu), and heteroscedastic variance models; while also cautioning against its application for M-W <= 4 in low-moderate seismicity regions without evaluating the homogeneity of M-W estimates between pan-European ESM and regional datasets.
Megathrust earthquakes impose changes of differential stress and pore pressure in the lithosphere-asthenosphere system that are transiently relaxed during the postseismic period primarily due to afterslip, viscoelastic and poroelastic processes.
Especially during the early postseismic phase, however, the relative contribution of these processes to the observed surface deformation is unclear.
To investigate this, we use geodetic data collected in the first 48 days following the 2010 Maule earthquake and a poro-viscoelastic forward model combined with an afterslip inversion.
This model approach fits the geodetic data 14% better than a pure elastic model. Particularly near the region of maximum coseismic slip, the predicted surface poroelastic uplift pattern explains well the observations.
If poroelasticity is neglected, the spatial afterslip distribution is locally altered by up to +/- 40%.
Moreover, we find that shallow crustal aftershocks mostly occur in regions of increased postseismic pore-pressure changes, indicating that both processes might be mechanically coupled.
Improving nitrogen (N) status in European water bodies is a pressing issue. N levels depend not only on current but also past N inputs to the landscape, that have accumulated through time in legacy stores (e.g., soil, groundwater).
Catchment-scale N models, that are commonly used to investigate in-stream N levels, rarely examine the magnitude and dynamics of legacy components.
This study aims to gain a better understanding of the long-term fate of the N inputs and its uncertainties, using a legacy-driven N model (ELEMeNT) in Germany's largest national river basin (Weser; 38,450 km(2)) over the period 1960-2015.
We estimate the nine model parameters based on a progressive constraining strategy, to assess the value of different observational data sets.
We demonstrate that beyond in-stream N loading, soil N content and in-stream N concentration allow to reduce the equifinality in model parameterizations.
We find that more than 50% of the N surplus denitrifies (1480-2210 kg ha(-1)) and the stream export amounts to around 18% (410-640 kg ha(-1)), leaving behind as much as around 230-780 kg ha(-1) of N in the (soil) source zone and 10-105 kg ha(-1) in the subsurface.
A sensitivity analysis reveals the importance of different factors affecting the residual uncertainties in simulated N legacies, namely hydrologic travel time, denitrification rates, a coefficient characterizing the protection of organic N in source zone and N surplus input.
Our study calls for proper consideration of uncertainties in N legacy characterization, and discusses possible avenues to further reduce the equifinality in water quality modeling.
Marine sedimentary archives are routinely used to reconstruct past environmental changes. In many cases, bioturbation and sedimentary mixing affect the proxy time-series and the age-depth relationship. While idealized models of bioturbation exist, they usually assume homogeneous mixing, thus that a single sample is representative for the sediment layer it is sampled from.
However, it is largely unknown to which extent this assumption holds for sediments used for paleoclimate reconstructions.
To shed light on
1) the age-depth relationship and its full uncertainty,
2) the magnitude of mixing processes affecting the downcore proxy variations, and
3) the representativity of the discrete sample for the sediment layer, we designed and performed a case study on South China Sea sediment material which was collected using a box corer and which covers the last glacial cycle.
Using the radiocarbon content of foraminiferal tests as a tracer of time, we characterize the spatial age-heterogeneity of sediments in a three-dimensional setup. In total, 118 radiocarbon measurements were performed on defined small- and large-volume bulk samples ( similar to 200 specimens each) to investigate the horizontal heterogeneity of the sediment. Additionally, replicated measurements on small numbers of specimens (10 x 5 specimens) were performed to assess the heterogeneity within a sample volume. Visual assessment of X-ray images and a quantitative assessment of the mixing strength show typical mixing from bioturbation corresponding to around 10 cm mixing depth.
Notably, our 3D radiocarbon distribution reveals that the horizontal heterogeneity (up to 1,250 years), contributing to the age uncertainty, is several times larger than the typically assumed radiocarbon based age-model error (single errors up to 250 years). Furthermore, the assumption of a perfectly bioturbated layer with no mixing underneath is not met.
Our analysis further demonstrates that the age-heterogeneity might be a function of sample size; smaller samples might contain single features from the incomplete mixing and are thus less representative than larger samples.
We provide suggestions for future studies, optimal sampling strategies for quantitative paleoclimate reconstructions and realistic uncertainty in age models, as well as discuss possible implications for the interpretation of paleoclimate records.
Soil bacteria play a fundamental role in pedogenesis. However, knowledge about both the impact of climate and slope aspects on microbial communities and the consequences of these items in pedogenesis is lacking. Therefore, soil-bacterial communities from four sites and two different aspects along the climate gradient of the Chilean Coastal Cordillera were investigated. Using a combination of microbiological and physicochemical methods, soils that developed in arid, semi-arid, mediterranean, and humid climates were analyzed. Proteobacteria, Acidobacteria, Chloroflexi, Verrucomicrobia, and Planctomycetes were found to increase in abundance from arid to humid climates, while Actinobacteria and Gemmatimonadetes decreased along the transect. Bacterial-community structure varied with climate and aspect and was influenced by pH, bulk density, plant-available phosphorus, clay, and total organic-matter content. Higher bacterial specialization was found in arid and humid climates and on the south-facing slope and was likely promoted by stable microclimatic conditions. The presence of specialists was associated with ecosystem-functional traits, which shifted from pioneers that accumulated organic matter in arid climates to organic decomposers in humid climates. These findings provide new perspectives on how climate and slope aspects influence the composition and functional capabilities of bacteria, with most of these capabilities being involved in pedogenetic processes.
Shallow earthquakes frequently disturb the hydrological and mechanical state of the subsurface, with consequences for hazard and water management. Transient post-seismic hydrological behavior has been widely reported, suggesting that the recovery of material properties (relaxation) following ground shaking may impact groundwater fluctuations. However, the monitoring of seismic velocity variations associated with earthquake damage and hydrological variations are often done assuming that both effects are independent. In a field site prone to highly variable hydrological conditions, we disentangle the different forcing of the relative seismic velocity variations delta v retrieved from a small dense seismic array in Nepal in the aftermath of the 2015 M-w 7.8 Gorkha earthquake. We successfully model transient damage effects by introducing a universal relaxation function that contains a unique maximum relaxation timescale for the main shock and the aftershocks, independent of the ground shaking levels. Next, we remove the modeled velocity from the raw data and test whether the corresponding residuals agree with a background hydrological behavior we inferred from a previously calibrated groundwater model. The fitting of the delta v data with this model is improved when we introduce transient hydrological properties in the phase immediately following the main shock. This transient behavior, interpreted as an enhanced permeability in the shallow subsurface, lasts for similar to 6 months and is shorter than the damage relaxation (similar to 1 yr). Thus, we demonstrate the capability of seismic interferometry to deconvolve transient hydrological properties after earthquakes from non-linear mechanical recovery.
Rainfall-intense summer monsoon seasons on the Indian subcontinent that are exceeding long-term averages cause widespread floods and landslides.
Here we show that the latest generation of coupled climate models robustly project an intensification of very rainfall-intense seasons (June-September).
Under the shared socioeconomic pathway SSP5-8.5, very wet monsoon seasons as observed in only 5 years in the period 1965-2015 are projected to occur 8 times more often in 2050-2100 in the multi-model average.
Under SSP2-4.5, these seasons become only a factor of 6 times more frequent, showing that even modest efforts to mitigate climate change can have a strong impact on the frequency of very strong rainfall seasons.
Besides, we find that the increasing risk of extreme seasonal rainfall is accompanied by a shift from days with light rainfall to days with moderate or heavy rainfall. Additionally, the number of wet days is projected to increase.
Phase transitions in molecular crystals are often determined by intermolecular interactions. The cage complex of [Co(C12H30N8)](3+) . 3 NO3- is reported to undergo a disorder-order phase transition at T-c1 approximate to 133 K upon cooling. Temperature-dependent neutron and synchrotron diffraction experiments revealed satellite reflections in addition to main reflections in the diffraction patterns below T-c1. The modulation wave vector varies as function of temperature and locks in at T-c3 approximate to 98 K. Here, we demonstrate that the crystal symmetry lowers from hexagonal to monoclinic in the incommensurately modulated phases in T-c1<T<T-c3. Distinctive levels of competitions: trade-off between longer N-H...O and shorter C-H...O hydrogen bonds; steric constraints to dense C-H...O bonds give rise to pronounced modulation of the basic structure. Severely frustrated crystal packing in the incommensurate phase is precursor to optimal balance of intermolecular interactions in the lock-in phase.
A detailed analysis of horizontal and vertical particulate matter (PM) fluxes during wind erosion has been done, based on measurements of PM smaller than 10, 2.5, and 1.0 mu mm, at windward and leeward positions on a measuring field. The three fractions of PM measurement are differently influenced by the increasing wind and shear velocities of the wind. The measured concentrations of the coarser fractions of the fine dust, PM10, and PM2.5, increase with wind and shear velocity, whereas the PM1.0 concentrations show no clear correlation to the shear velocity. The share of PM2.5 on PM10 depends on the measurement height and wind speed and varies between 4 and 12 m/s at the 1 m height ranging from 25% to 7% (average 10%), and at the 4 m height from 39% to 23% (average 30%). Although general relationships between wind speed, PM concentration, and horizontal and vertical fluxes could be found, the contribution of the measuring field was very low, as balances of incoming and outgoing fluxes show. Consequently, the measured PM concentrations are determined from a variety of sources, such as traffic on unpaved roads, cattle drives, tillage operations, and wind erosion, and thus, represent all components of land use and landscape structure in the near and far surroundings of the measuring field. The current results may reflect factors from the landscape scale rather than the influence of field-related variables. The measuring devices used to monitor PM concentrations showed differences of up to 20%, which led to considerable deviations when determining total balances. Differences up to 67% between the calculated fluxes prove the necessity of a previous calibration of the devices used. (c) 2022 International Research and Training Centre on Erosion and Sedimentation/the World Association for Sedimentation and Erosion Research.
The quantification of spatial propagation of extreme precipitation events is vital in water resources planning and disaster mitigation. However, quantifying these extreme events has always been challenging as many traditional methods are insufficient to capture the nonlinear interrelationships between extreme event time series. Therefore, it is crucial to develop suitable methods for analyzing the dynamics of extreme events over a river basin with a diverse climate and complicated topography. Over the last decade, complex network analysis emerged as a powerful tool to study the intricate spatiotemporal relationship between many variables in a compact way. In this study, we employ two nonlinear concepts of event synchronization and edit distance to investigate the extreme precipitation pattern in the Ganga river basin. We use the network degree to understand the spatial synchronization pattern of extreme rainfall and identify essential sites in the river basin with respect to potential prediction skills. The study also attempts to quantify the influence of precipitation seasonality and topography on extreme events. The findings of the study reveal that (1) the network degree is decreased in the southwest to northwest direction, (2) the timing of 50th percentile precipitation within a year influences the spatial distribution of degree, (3) the timing is inversely related to elevation, and (4) the lower elevation greatly influences connectivity of the sites. The study highlights that edit distance could be a promising alternative to analyze event-like data by incorporating event time and amplitude and constructing complex networks of climate extremes.
On their way from inland to the ocean, flowing water bodies, their constituents and their biotic communities are ex-posed to complex transport and transformation processes. However, detailed process knowledge as revealed by La-grangian measurements adjusted to travel time is rare in large rivers, in particular at hydrological extremes. To fill this gap, we investigated autotrophic processes, heterotrophic carbon utilization, and micropollutant concentrations applying a Lagrangian sampling design in a 600 km section of the River Elbe (Germany) at historically low discharge. Under base flow conditions, we expect the maximum intensity of instream processes and of point source impacts. Phy-toplankton biomass and photosynthesis increased from upstream to downstream sites but maximum chlorophyll con-centration was lower than at mean discharge. Concentrations of dissolved macronutrients decreased to almost complete phosphate depletion and low nitrate values. The longitudinal increase of bacterial abundance and production was less pronounced than in wetter years and bacterial community composition changed downstream. Molecular analyses revealed a longitudinal increase of many DOM components due to microbial production, whereas saturated lipid-like DOM, unsaturated aromatics and polyphenols, and some CHOS surfactants declined. In decomposition exper-iments, DOM components with high O/C ratios and high masses decreased whereas those with low O/C ratios, low masses, and high nitrogen content increased at all sites. Radiocarbon age analyses showed that DOC was relatively old (890-1870 years B.P.), whereas the mineralized fraction was much younger suggesting predominant oxidation of algal lysis products and exudates particularly at downstream sites. Micropollutants determining toxicity for algae (terbuthylazine, terbutryn, isoproturon and lenacil), hexachlorocyclohexanes and DDTs showed higher concentrations from the middle towards the downstream part but calculated toxicity was not negatively correlated to phytoplankton. Overall, autotrophic and heterotrophic process rates and micropollutant concentrations increased from up-to down-stream reaches, but their magnitudes were not distinctly different to conditions at medium discharges.
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.
The James Ross archipelago houses numerous lakes and ponds. In this region, a vast diatom and cyanobacterial variety has been reported; however, the prokaryotic diversity in microbial mats from these lakes remains poorly explored.
Here, a high-throughput sequencing of 16S rRNA gene in microbial mats from Lake Bart-Roja in James Ross Island and lakes Pan Negro and North Pan Negro located in Vega Island was performed. Combined with mineralogical and environmental characteristics, we analyzed the diversity and structure of the microbial communities. Sequences assigned to Archaea were extremely low, while Bacteria domain prevailed with the abundance of Proteobacteria (mostly Betaproteobacteriales) followed by Bacteroidetes, Verrucomicrobia, Firmicutes, and Cyanobacteria.
Local environmental conditions, such as conductivity and Eh, provided differential microbial assemblages that might have implications in the oligotrophic status of the lakes. Consequently, a clear segregation at the family level was observed.
In this sense, the assigned diversity was related to taxa recognized as denitrifiers and sulfur oxidizers. Particularly, in Lake Pan Negro sulfur-reducing and methanogenic representatives were also found and positively correlate with alkalinity and water depth.
Moreover, Deinococcus-Thermus was observed in Lake Bart-Roja, while Melainabacteria (Cyanobacteria)-poorly reported in Antarctic mats-was detected in Lake Pan Negro. Epsilonbacteraeota was exclusively found in this lake, suggesting new potential phylotypes. This study contributes to the understanding of the diversity, composition, and structure of Antarctic benthic microbial ecosystems and provides highly valuable information, which can be used as a proxy to evaluate environmental changes affecting Antarctic microbiota.
Large-Scale interseismic strain mapping of the NE Tibetan Plateau from Sentinel-1 Interferometry
(2022)
The launches of the Sentinel-1 synthetic aperture radar satellites in 2014 and 2016 started a new era of high-resolution velocity and strain rate mapping for the continents. However, multiple challenges exist in tying independently processed velocity data sets to a common reference frame and producing high-resolution strain rate fields. We analyze Sentinel-1 data acquired between 2014 and 2019 over the northeast Tibetan Plateau, and develop new methods to derive east and vertical velocities with similar to 100 m resolution and similar to 1 mm/yr accuracy across an area of 440,000 km(2). By implementing a new method of combining horizontal gradients of filtered east and interpolated north velocities, we derive the first similar to 1 km resolution strain rate field for this tectonically active region. The strain rate fields show concentrated shear strain along the Haiyuan and East Kunlun Faults, and local contractional strain on fault junctions, within the Qilianshan thrusts, and around the Longyangxia Reservoir. The Laohushan-Jingtai creeping section of the Haiyuan Fault is highlighted in our data set by extremely rapid strain rates. Strain across unknown portions of the Haiyuan Fault system, including shear on the eastern extension of the Dabanshan Fault and contraction at the western flank of the Quwushan, highlight unmapped tectonic structures. In addition to the uplift across most of the lowlands, the vertical velocities also contain climatic, hydrological or anthropogenic-related deformation signals. We demonstrate the enhanced view of large-scale active tectonic processes provided by high-resolution velocities and strain rates derived from Sentinel-1 data and highlight associated wide-ranging research applications.
In seismic risk assessment, the sources of uncertainty associated with building exposure modelling have not received as much attention as other components related to hazard and vulnerability. Conventional practices such as assuming absolute portfolio compositions (i.e., proportions per building class) from expert-based assumptions over aggregated data crudely disregard the contribution of uncertainty of the exposure upon earthquake loss models. In this work, we introduce the concept that the degree of knowledge of a building stock can be described within a Bayesian probabilistic approach that integrates both expert-based prior distributions and data collection on individual buildings. We investigate the impact of the epistemic uncertainty in the portfolio composition on scenario-based earthquake loss models through an exposure-oriented logic tree arrangement based on synthetic building portfolios. For illustrative purposes, we consider the residential building stock of Valparaiso (Chile) subjected to seismic ground-shaking from one subduction earthquake. We have found that building class reconnaissance, either from prior assumptions by desktop studies with aggregated data (top-down approach), or from building-by-building data collection (bottom-up approach), plays a fundamental role in the statistical modelling of exposure. To model the vulnerability of such a heterogeneous building stock, we require that their associated set of structural fragility functions handle multiple spectral periods. Thereby, we also discuss the relevance and specific uncertainty upon generating either uncorrelated or spatially cross-correlated ground motion fields within this framework. We successively show how various epistemic uncertainties embedded within these probabilistic exposure models are differently propagated throughout the computed direct financial losses. This work calls for further efforts to redesign desktop exposure studies, while also highlighting the importance of exposure data collection with standardized and iterative approaches.
Understanding catchment controls on catchment solute export is a prerequisite for water quality management. StorAge Selection (SAS) functions encapsulate essential information about catchment functioning in terms of discharge selection preference and solute export dynamics. However, they lack information on the spatial origin of solutes when applied at the catchment scale, thereby limiting our understanding of the internal (subcatchment) functioning. Here, we parameterized SAS functions in a spatially explicit way to understand the internal catchment responses and transport dynamics of reactive dissolved nitrate (N-NO3). The model was applied in a nested mesoscale catchment (457 km(2)), consisting of a mountainous partly forested, partly agricultural subcatchment, a middle-reach forested subcatchment, and a lowland agricultural subcatchment. The model captured flow and nitrate concentration dynamics not only at the catchment outlet but also at internal gauging stations. Results reveal disparate subsurface mixing dynamics and nitrate export among headwater and lowland subcatchments. The headwater subcatchment has high seasonal variation in subsurface mixing schemes and younger water in discharge, while the lowland subcatchment has less pronounced seasonality in subsurface mixing and much older water in discharge. Consequently, nitrate concentration in discharge from the headwater subcatchment shows strong seasonality, whereas that from the lowland subcatchment is stable in time. The temporally varying responses of headwater and lowland subcatchments alternate the dominant contribution to nitrate export in high and low-flow periods between subcatchments. Overall, our results demonstrate that the spatially explicit SAS modeling provides useful information about internal catchment functioning, helping to develop or evaluate spatial management practices.
The region of West Bohemia and Upper Palatinate belongs to the West Bohemian Massif. The study area is situated at the junction of three different Variscan tectonic units and hosts the ENE-WSW trending Ohre Rift as well as many different fault systems. The entire region is characterized by ongoing magmatic processes in the intra-continental lithospheric mantle expressed by a series of phenomena, including e.g. the occurrence of repeated earthquake swarms and massive degassing of mantle derived CO2 in form of mineral springs and mofettes. Ongoing active tectonics is mainly manifested by Cenozoic volcanism represented by different Quaternary volcanic structures. All these phenomena make the Ohre Rift a unique target area for European intra-continental geo-scientific research. With magnetotelluric (MT) measurements we image the subsurface distribution of the electrical resistivity and map possible fluid pathways. Two-dimensional (2D) inversion results by Munoz et al. (2018) reveal a conductive channel in the vicinity of the earthquake swarm region that extends from the lower crust to the surface forming a pathway for fluids into the region of the mofettes. A second conductive channel is present in the south of their model; however, their 2D inversions allow ambiguous interpretations of this feature. Therefore, we conducted a large 3D MT field experiment extending the study area towards the south. The 3D inversion result matches well with the known geology imaging different fluid/magma reservoirs at crust-mantle depth and mapping possible fluid pathways from the reservoirs to the surface feeding known mofettes and spas. A comparison of 3D and 2D inversion results suggests that the 2D inversion results are considerably characterized by 3D and off-profile structures. In this context, the new results advocate for the swarm earthquakes being located in the resistive host rock surrounding the conductive channels; a finding in line with observations e.g. at the San Andreas Fault, California.
We construct and examine the prototype of a deep learning-based ground-motion model (GMM) that is both fully data driven and nonergodic. We formulate ground-motion modeling as an image processing task, in which a specific type of neural network, the U-Net, relates continuous, horizontal maps of earthquake predictive parameters to sparse observations of a ground-motion intensity measure (IM). The processing of map-shaped data allows the natural incorporation of absolute earthquake source and observation site coordinates, and is, therefore, well suited to include site-, source-, and path-specific amplification effects in a nonergodic GMM. Data-driven interpolation of the IM between observation points is an inherent feature of the U-Net and requires no a priori assumptions. We evaluate our model using both a synthetic dataset and a subset of observations from the KiK-net strong motion network in the Kanto basin in Japan. We find that the U-Net model is capable of learning the magnitude???distance scaling, as well as site-, source-, and path-specific amplification effects from a strong motion dataset. The interpolation scheme is evaluated using a fivefold cross validation and is found to provide on average unbiased predictions. The magnitude???distance scaling as well as the site amplification of response spectral acceleration at a period of 1 s obtained for the Kanto basin are comparable to previous regional studies.
The main Marmara fault (MMF) extends for 150 km through the Sea of Marmara and forms the only portion of the North Anatolian fault zone that has not ruptured in a large event (Mw >7) for the last 250 yr. Accordingly, this portion is potentially a major source contributing to the seismic hazard of the Istanbul region. On 26 September 2019, a sequence of moderate-sized events started along the MMF only 20 km south of Istanbul and were widely felt by the population. The largest three events, 26 September Mw 5.8 (10:59 UTC), 26 September 2019 Mw 4.1 (11:26 UTC), and 20 January 2020 Mw 4.7 were recorded by numerous strong-motion seismic stations and the resulting ground motions were compared to the predicted means resulting from a set of the most recent ground-motion prediction equations (GMPEs). The estimated residuals were used to investigate the spatial variation of ground motion across the Marmara region. Our results show a strong azimuthal trend in ground-motion residuals, which might indicate systematically repeating directivity effects toward the eastern Marmara region.
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.
Lithium and boron are trace components of magmas, released during exsolution of a gas phase during volcanic activity.
In this study, we determine the diffusivity and isotopic fractionation of Li and B in hydrous silicate melts.
Two glasses were synthesized with the same rhyolitic composition (4.2 wt% water), having different Li and B contents; these were studied in diffusion-couple experiments that were performed using an internally heated pressure vessel, operated at 300 MPa in the temperature range 700-1250 degrees C for durations from 0 s to 24 h. From this we determined activation energies for Li and B diffusion of 57 +/- 4 kJ/mol and 152 +/- 15 kJ/mol with pre-exponential factors of 1.53 x 10(-7) m(2)/s and 3.80 x 10(-8) m(2)/s, respectively.
Lithium isotopic fractionation during diffusion gave beta values between 0.15 and 0.20, whereas B showed no clear isotopic fractionation.
Our Li diffusivities and isotopic fractionation results differ somewhat from earlier published values, but overall confirm that Li diffusivity increases with water content. Our results on B diffusion show that similarly to Li, B mobility increases in the presence of water.
By applying the Eyring relation, we confirm that B diffusivity is limited by viscous flow in silicate melts.
Our results on Li and B diffusion present a new tool for understanding degassing-related processes, offering a potential geospeedometer to measure volcanic ascent rates.
Data recorded by distributed acoustic sensing (DAS) along an optical fibre sample the spatial and temporal properties of seismic wavefields at high spatial density. Often leading to massive amount of data when collected for seismic monitoring along many kilometre long cables. The spatially coherent signals from weak seismic arrivals within the data are often obscured by incoherent noise. We present a flexible and computationally efficient filtering technique, which makes use of the dense spatial and temporal sampling of the data and that can handle the large amount of data. The presented adaptive frequency-wavenumber filter suppresses the incoherent seismic noise while amplifying the coherent wavefield. We analyse the response of the filter in time and spectral domain, and we demonstrate its performance on a noisy data set that was recorded in a vertical borehole observatory showing active and passive seismic phase arrivals. Lastly, we present a performant open-source software implementation enabling real-time filtering of large DAS data sets.
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.
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 creation of building exposure models for seismic risk assessment is frequently challenging due to the lack of availability of detailed information on building structures. Different strategies have been developed in recent years to overcome this, including the use of census data, remote sensing imagery and volunteered graphic information (VGI). This paper presents the development of a building-by-building exposure model based exclusively on openly available datasets, including both VGI and census statistics, which are defined at different levels of spatial resolution and for different moments in time. The initial model stemming purely from building-level data is enriched with statistics aggregated at the neighbourhood and city level by means of a Monte Carlo simulation that enables the generation of full realisations of damage estimates when using the exposure model in the context of an earthquake scenario calculation. Though applicable to any other region of interest where analogous datasets are available, the workflow and approach followed are explained by focusing on the case of the German city of Cologne, for which a scenario earthquake is defined and the potential damage is calculated. The resulting exposure model and damage estimates are presented, and it is shown that the latter are broadly consistent with damage data from the 1978 Albstadt earthquake, notwithstanding the differences in the scenario. Through this real-world application we demonstrate the potential of VGI and open data to be used for exposure modelling for natural risk assessment, when combined with suitable knowledge on building fragility and accounting for the inherent uncertainties.
Volcano-seismic signals such as long-period events and tremor are important indicators for volcanic activity and unrest. However, their wavefield is complex and characterization and location using traditional seismological instrumentation is often difficult.
In 2019 we recorded the full seismic wavefield using a newly developed 3C rotational sensor co-located with a 3C traditional seismometer on Etna, Italy. We compare the performance of the rotational sensor, the seismometer and the Istituto Nazionale di Geofisica e Vulcanologia-Osservatorio Etneo (INGV-OE) seismic network with respect to the analysis of complex volcano-seismic signals. We create event catalogs for volcano-tectonic (VT) and long-period (LP) events combining a STA/LTA algorithm and cross-correlations.
The event detection based on the rotational sensor is as reliable as the seismometer-based detection. The LP events are dominated by SH-type waves. Derived SH phase velocities range from 500 to 1,000 m/s for LP events and 300-400 m/s for volcanic tremor. SH-waves compose the tremor during weak volcanic activity and SH- and SV-waves during sustained strombolian activity.
We derive back azimuths using (a) horizontal rotational components and (b) vertical rotation rate and transverse acceleration. The estimated back azimuths are consistent with the INGV-OE event location for (a) VT events with an epicentral distance larger than 3 km and some closer events, (b) LP events and tremor in the main crater area. Measuring the full wavefield we can reliably analyze the back azimuths, phase velocities and wavefield composition for VT, LP events and tremor in regions that are difficult to access such as volcanoes.
Regional flood frequency analysis (RFFA) is a powerful method for interrogating hydrological series since it combines observational time series from several sites within a region to estimate risk-relevant statistical parameters with higher accuracy than from single-site series. Since RFFA extreme value estimates depend on the shape of the selected distribution of the data-generating stochastic process, there is need for a suitable goodness-of-distributional-fit measure in order to optimally utilize given data. Here we present a novel, least-squares-based measure to select the optimal fit from a set of five distributions, namely Generalized Extreme Value (GEV), Generalized Logistic, Gumbel, Log-Normal Type III and Log-Pearson Type III. The fit metric is applied to annual maximum discharge series from six hydrological stations along the Sava River in South-eastern Europe, spanning the years 1961 to 2020. Results reveal that (1) the Sava River basin can be assessed as hydrologically homogeneous and (2) the GEV distribution provides typically the best fit. We offer hydrological-meteorological insights into the differences among the six stations. For the period studied, almost all stations exhibit statistically insignificant trends, which renders the conclusions about flood risk as relevant for hydrological sciences and the design of regional flood protection infrastructure.
Purpose:
Soil erosion by water yields sediment to surface reservoirs, reducing their storage capacities, changing their geometry, and degrading water quality. Sediment reuse, i.e., fertilization of agricultural soils with the nutrient-enriched sediment from reservoirs, has been proposed as a recovery strategy. However, the sediment needs to meet certain criteria. In this study, we characterize sediments from the densely dammed semiarid Northeast Brazil by VNIR-SWIR spectroscopy and assess the effect of spectral resolution and spatial scale on the accuracy of N, P, K, C, electrical conductivity, and clay prediction models.
Methods
Sediment was collected in 10 empty reservoirs, and physical and chemical laboratory analyses as well as spectral measurements were performed. The spectra, initially measured at 1 nm spectral resolution, were resampled to 5 and 10 nm, and samples were analysed for both high and low spectral resolution at three spatial scales, namely (1) reservoir, (2) catchment, and (3) regional scale.
Results
Partial least square regressions performed from good to very good in the prediction of clay and electrical conductivity from reservoir (<40 km(2)) to regional (82,500 km(2)) scales. Models for C and N performed satisfactorily at the reservoir scale, but degraded to unsatisfactory at the other scales. Models for P and K were more unstable and performed from unsatisfactorily to satisfactorily at all scales. Coarsening spectral resolution by up to 10 nm only slightly degrades the models' performance, indicating the potential of characterizing sediment from spectral data captured at lower resolutions, such as by hyperspectral satellite sensors.
Conclusion:
By reducing the costly and time-consuming laboratory analyses, the method helps to promote the sediment reuse as a practice of soil and water conservation.
The Walker Circulation (WC) is an east-west trending band of atmospheric circulation cells along the equator and the predominant controller of heat and moisture transport in the tropics. Its variability is closely linked to the sea-surface temperature (SST) changes across the Pacific, the Indian and the Atlantic Oceans and can have pronounced effects on the humidity regimes of the adjacent continents. In recent years, the evolution of the WC during the Plioand Pleistocene epochs has been intensely studied in the context of the effectiveness of the tropics in modulating global climate change (e.g., the intensification of Northern Hemisphere glaciation). However, the onset of the modern WC pattern as well as its global impact during the Plioand Pleistocene is controversially assessed in the literature. For its onset, previous studies have suggested dates ranging between 2.4 and 0.8 million years ago (Myr), while its argued impact ranges from crucially influencing the increase of Northern Hemisphere ice sheet growth by channelling heat and moisture from the tropics into the high latitudes to having no effect on global ice volume changes. In order to achieve a comprehensive understanding of the spatiotemporal evolution of the WC during this time frame, we statistically analysed 30 globally distributed SST records covering the low and high latitudes between 3.5 and 1.5 Myr, encompassing the Late Pliocene to Early Pleistocene. We utilized a statistical change-point regression model to determine significant change points in the SST evolution of the (sub)-tropics and high latitudes that potentially relate to changes in the WC. We find that the WC experienced a multifaceted evolution between the Late Pliocene and the Early Pleistocene with significant transitional steps at-2.7 and-2.1 Ma. Our results suggest after the Late Pliocene, a pre-modern WC set in, which was characterized by a progressively strengthened Pacific Walker Cell alongside a weakened Indian Ocean Walker Cell. This change was potentially triggered by the constriction of the Indonesian seaway, an important transmitter between the Pacific and Indian Ocean. The ensuing mode of the WC intensified until-2.1 Myr, when SST values around the global scale signalled a progressive strengthening of the Indian Walker Cell in phase with the progressive strengthening of the Pacific and Atlantic Cells. Our findings indicate that a shift from a pre-modern to a modern-like WC potentially only occurred during the mid-Pleistocene.
The Arctic is rich in aquatic systems and experiences rapid warming due to climate change. The accelerated warming causes permafrost thaw and the mobilization of organic carbon.
When dissolved organic carbon is mobilized, this DOC can be transported to aquatic systems and degraded in the water bodies and further downstream. Here, we analyze the influence of different landscape components on DOC concentrations and export in a small (6.45 km(2)) stream catchment in the Lena River Delta.
The catchment includes lakes and ponds, with the flow path from Pleistocene yedoma deposits across Holocene non-yedoma deposits to the river outlet. In addition to DOC concentrations, we use radiocarbon dating of DOC as well as stable oxygen and hydrogen isotopes (delta O-18 and delta D) to assess the origin of DOC.
We find significantly higher DOC concentrations in the Pleistocene yedoma area of the catchment compared to the Holocene non-yedoma area with medians of 5 and 4.5 mg L-1 (p < 0.05), respectively. When yedoma thaw streams with high DOC concentration reach a large yedoma thermokarst lake, we observe an abrupt decrease in DOC concentration, which we attribute to dilution and lake processes such as mineralization. The DOC ages in the large thermokarst lake (between 3,428 and 3,637 C-14 y BP) can be attributed to a mixing of mobilized old yedoma and Holocene carbon.
Further downstream after the large thermokarst lake, we find progressively younger DOC ages in the stream water to its mouth, paired with decreasing DOC concentrations. This process could result from dilution with leaching water from Holocene deposits and/or emission of ancient yedoma carbon to the atmosphere. Our study shows that thermokarst lakes and ponds may act as DOC filters, predominantly by diluting incoming waters of higher DOC concentrations or by re-mineralizing DOC to CO2 and CH4.
Nevertheless, our results also confirm that the small catchment still contributes DOC on the order of 1.2 kg km(-2) per day from a permafrost landscape with ice-rich yedoma deposits to the Lena River.
Ground motion with strong-velocity pulses can cause significant damage to buildings and structures at certain periods; hence, knowing the period and velocity amplitude of such pulses is critical for earthquake structural engineering.
However, the physical factors relating the scaling of pulse periods with magnitude are poorly understood.
In this study, we investigate moderate but damaging earthquakes (M-w 6-7) and characterize ground- motion pulses using the method of Shahi and Baker (2014) while considering the potential static-offset effects.
We confirm that the within-event variability of the pulses is large. The identified pulses in this study are mostly from strike-slip-like earthquakes. We further perform simulations using the freq uency-wavenumber algorithm to investigate the causes of the variability of the pulse periods within and between events for moderate strike-slip earthquakes.
We test the effect of fault dips, and the impact of the asperity locations and sizes. The simulations reveal that the asperity properties have a high impact on the pulse periods and amplitudes at nearby stations.
Our results emphasize the importance of asperity characteristics, in addition to earthquake magnitudes for the occurrence and properties of pulses produced by the forward directivity effect.
We finally quantify and discuss within- and between-event variabilities of pulse properties at short distances.
It is widely recognized that collisional mountain belt topography is generated by crustal thickening and lowered by river bedrock erosion, linking climate and tectonics(1-4). However, whether surface processes or lithospheric strength control mountain belt height, shape and longevity remains uncertain. Additionally, how to reconcile high erosion rates in some active orogens with long-term survival of mountain belts for hundreds of millions of years remains enigmatic. Here we investigate mountain belt growth and decay using a new coupled surface process(5,6) and mantle-scale tectonic model(7). End-member models and the new non-dimensional Beaumont number, Bm, quantify how surface processes and tectonics control the topographic evolution of mountain belts, and enable the definition of three end-member types of growing orogens: type 1, non-steady state, strength controlled (Bm > 0.5); type 2, flux steady state(8), strength controlled (Bm approximate to 0.4-0.5); and type 3, flux steady state, erosion controlled (Bm < 0.4). Our results indicate that tectonics dominate in Himalaya-Tibet and the Central Andes (both type 1), efficient surface processes balance high convergence rates in Taiwan (probably type 2) and surface processes dominate in the Southern Alps of New Zealand (type 3). Orogenic decay is determined by erosional efficiency and can be subdivided into two phases with variable isostatic rebound characteristics and associated timescales. The results presented here provide a unified framework explaining how surface processes and lithospheric strength control the height, shape, and longevity of mountain belts.
Modern pollen-vegetation-climate relationships underpin palaeovegetation and palaeoclimate reconstructions from fossil pollen records. East Siberia is an ideal area for investigating the relationships between modern pollen assemblages and near natural vegetation under cold continental climate conditions. Reliable pollen-based quantitative vegetation and climate reconstructions are still scarce due to the limited number of modern pollen datasets. Furthermore, differences in pollen representation of samples from lake sediments and soils are not well understood. Here, we present a new pollen dataset of 48 moss/soil and 24 lake surface-sediment samples collected in Chukotka and central Yakutia in East Siberia. The pollen-vegetation-climate relationships were investigated by ordination analyses. Generally, tundra and taiga vegetation types can be well distinguished in the surface pollen assemblages. Moss/soil and lake samples contain generally similar pollen assemblages as revealed by a Procrustes comparison with some exceptions. Overall, modern pollen assemblages reflect the temperature and precipitation gradients in the study areas as revealed by constrained ordination analysis. We estimate the relative pollen productivity (RPP) of major taxa and the relevant source area of pollen (RSAP) for moss/soil samples from Chukotka and central Yakutia using Extended R-Value (ERV) analysis. The RSAP of the tundra-forest transition area in Chukotka and taiga area in central Yakutia are ca. 1300 and 360 m, respectively. For Chukotka, RPPs relative to both Poaceae and Ericaceae were estimated while RPPs for central Yakutia were relative only to Ericaceae. Relative to Ericaceae (reference taxon, RPP = 1), Larix, Betula, Picea, and Pinus are overrepresented while Alnus, Cyperaceae, Poaceae, and Salix are underrepresented in the pollen spectra. Our estimates are in general agreement with previously published values and provide the basis for reliable quantitative reconstructions of East Siberian vegetation.
Urban surface runoff management via best management practices (BMP) and low impact development (LID) has earned significant recognition owing to positive environmental and ecological impacts. However, due to the complexity of the parameters involved, the estimation of LID efficiency in attenuating the urban surface runoff at the watershed scale is challenging. A planning analysis of employing Green Roofs and Infiltration Trenches as BMPs/LIDs practices for urban surface runoff control is presented in this study. A multi-objective optimization decision-making framework is established by coupling SWMM (Storm Water Management Model) with NSGA-II models to check the performance of BMPs/LIDs concerning the cost-benefit analysis of LID at the watershed scale. Two urbanized areas belonging to Central Delhi in India were used as case studies. The results showed that the SWMM model is useful in simulating optimization problems for managing urban surface runoff. The optimum scenarios efficiently minimized the urban runoff volume while maintaining the BMPs/LIDs implementation costs and size. With BMPs/LIDs implementation, the reduction in runoff volume increases as expenses increase initially; however, there is no noticeable reduction in flood volume after a certain threshold. Contrasted with the haphazard arrangement of BMPs/LIDs, the proposed approach demonstrates 22%-24% runoff reductions for the same expenditures in watershed 1 and 23%-26% in watershed 2. The result of the study provides insights into planning and management of the urban surface runoff control with LID practices. The proposed framework assists the hydrologists in optimum selection and placements of BMPs/LIDs practices to acquire the most extreme ecological advantages with the least expenses.
Extracting information about past tectonic or climatic environmental changes from sedimentary records is a key objective of provenance research. Interpreting the imprint of such changes remains challenging as signals might be altered in the sediment-routing system.
We investigate the sedimentary provenance of the Oligocene/Miocene Upper Austrian Northern Alpine Foreland Basin and its response to the tectonically driven exhumation of the Tauern Window metamorphic dome (28 +/- 1 Ma) in the Eastern European Alps by using the unprecedented combination of Nd isotopic composition of bulk-rock clay-sized samples and partly previously published multi-proxy (Nd isotopic composition, trace-element geochemistry, U-Pb dating) sand-sized apatite single-grain analysis.
The basin offers an excellent opportunity to investigate environmental signal propagation into the sedimentary record because comprehensive stratigraphic and seismic datasets can be combined with present research results. The bulk-rock clay-sized fraction epsilon Nd values of well-cutting samples from one well on the northern basin slope remained stable at similar to-9.7 from 27 to 19 Ma but increased after 19 Ma to similar to-9.1. In contrast, apatite single-grain distributions, which were extracted from 22 drill-core samples, changed significantly around 23.3 Ma from apatites dominantly from low-grade (<upper amphibolite-facies) metamorphic sources with Permo-Mesozoic and late Variscan U-Pb ages and epsilon Nd values of -4.4 to dominantly high-grade metamorphic apatites with late Variscan U-Pb ages and epsilon Nd values of -2.2.
The change in apatite single-grain distributions at 23.3 Ma is interpreted to result from the exposure of a new Upper Austroalpine source nappe with less negative epsilon Nd values triggered by the ongoing Tauern Window exhumation. Combining these data with the clay-sized bulk-rock epsilon Nd values reveals that the provenance changed 4-5 Myrs later at 19 Ma in the clay-sized fraction.
Reasons for the delayed provenance-change recording are rooted in the characteristics of the applied methods.
Whereas single-grain distributions of orogen-wide sediment-routing systems can be dominated by geographically small areas with high erosion and mineral fertility rates, bulk-rock methods integrate over the entire drainage basin, thus diminishing extreme values. Hence, by combining these two methods, spatial information are uncovered, enabling a previously unattained understanding of the underlying environmental change.
A 3-D crustal shear wave velocity model and Moho map below the Semail Ophiolite, eastern Arabia
(2022)
The Semail Ophiolite in eastern Arabia is the largest and best-exposed slice of oceanic lithosphere on land. Detailed knowledge of the tectonic evolution of the shallow crust, in particular during and after ophiolite obduction in Late Cretaceous times is contrasted by few constraints on physical and compositional properties of the middle and lower continental crust below the obducted units. The role of inherited, pre-obduction crustal architecture remains therefore unaccounted for in our understanding of crustal evolution and the present-day geology. Based on seismological data acquired during a 27-month campaign in northern Oman, Ambient Seismic Noise Tomography and Receiver Function analysis provide for the first time a 3-D radially anisotropic shear wave velocity (V-S) model and a consistent Moho map below the iconic Semail Ophiolite. The model highlights deep crustal boundaries that segment the eastern Arabian basement in two distinct units. The previously undescribed Western Jabal Akhdar Zone separates Arabian crust with typical continental properties and a thickness of similar to 40-45 km in the northwest from a compositionally different terrane in the southeast that is interpreted as a terrane accreted during the Pan-African orogeny in Neoproterozoic times. East of the Ibra Zone, another deep crustal boundary, crustal thickness decreases to 30-35 km and very high lower crustal V-S suggest large-scale mafic intrusions into, and possible underplating of the Arabian continental crust that occurred most likely during Permian breakup of Pangea. Mafic reworking is sharply bounded by the (upper crustal) Semail Gap Fault Zone, northwest of which no such high velocities are found in the crust. Topography of the Oman Mountains is supported by a mild crustal root and Moho depth below the highest topography, the Jabal Akhdar Dome, is similar to 42 km. Radial anisotropy is robustly resolved in the upper crust and aids in discriminating dipping allochthonous units from autochthonous sedimentary rocks that are indistinguishable by isotropic V-S alone. Lateral thickness variations of the ophiolite highlight the Haylayn Ophiolite Massif on the northern flank of Jabal Akhdar Dome and the Hawasina Window as the deepest reaching unit. Ophiolite thickness is similar to 10 km in the southern and northern massifs, and <= 5 km elsewhere.