550 Geowissenschaften
Refine
Year of publication
- 2020 (261) (remove)
Document Type
- Article (182)
- Postprint (48)
- Doctoral Thesis (26)
- Other (2)
- Monograph/Edited Volume (1)
- Report (1)
- Review (1)
Is part of the Bibliography
- yes (261)
Keywords
- climate change (12)
- remote sensing (7)
- climate-change (6)
- dynamics (6)
- Europe (5)
- diversity (5)
- model (5)
- models (5)
- numerical modeling (5)
- Andes (4)
Institute
- Institut für Geowissenschaften (188)
- Institut für Umweltwissenschaften und Geographie (66)
- Institut für Biochemie und Biologie (8)
- Institut für Physik und Astronomie (7)
- Institut für Mathematik (3)
- Extern (2)
- Fachgruppe Volkswirtschaftslehre (1)
- Hasso-Plattner-Institut für Digital Engineering gGmbH (1)
- Potsdam Institute for Climate Impact Research (PIK) e. V. (1)
- Wirtschaftswissenschaften (1)
In the comment on "Varves of the Dead Sea sedimentary record." Quaternary Science Reviews 215 (Ben Dor et al., 2019): 173-184. by R. Bookman, two recently published papers are suggested to prove that the interpretation of the laminated sedimentary sequence of the Dead Sea, deposited mostly during MIS2 and Holocene pluvials, as annual deposits (i.e., varves) is wrong. In the following response, we delineate several lines of evidence which coalesce to demonstrate that based on the vast majority of evidence, including some of the evidence provided in the comment itself, the interpretation of these sediments as varves is the more likely scientific conclusion. We further discuss the evidence brought up in the comment and its irrelevance and lack of robustness for addressing the question under discussion.
Bayesian geomorphology
(2020)
The rapidly growing amount and diversity of data are confronting us more than ever with the need to make informed predictions under uncertainty. The adverse impacts of climate change and natural hazards also motivate our search for reliable predictions. The range of statistical techniques that geomorphologists use to tackle this challenge has been growing, but rarely involves Bayesian methods. Instead, many geomorphic models rely on estimated averages that largely miss out on the variability of form and process. Yet seemingly fixed estimates of channel heads, sediment rating curves or glacier equilibrium lines, for example, are all prone to uncertainties. Neighbouring scientific disciplines such as physics, hydrology or ecology have readily embraced Bayesian methods to fully capture and better explain such uncertainties, as the necessary computational tools have advanced greatly. The aim of this article is to introduce the Bayesian toolkit to scientists concerned with Earth surface processes and landforms, and to show how geomorphic models might benefit from probabilistic concepts. I briefly review the use of Bayesian reasoning in geomorphology, and outline the corresponding variants of regression and classification in several worked examples.
Strong hydroclimatic controls on vulnerability to subsurface nitrate contamination across Europe
(2020)
Subsurface contamination due to excessive nutrient surpluses is a persistent and widespread problem in agricultural areas across Europe. The vulnerability of a particular location to pollution from reactive solutes, such as nitrate, is determined by the interplay between hydrologic transport and biogeochemical transformations. Current studies on the controls of subsurface vulnerability do not consider the transient behaviour of transport dynamics in the root zone. Here, using state-of-the-art hydrologic simulations driven by observed hydroclimatic forcing, we demonstrate the strong spatiotemporal heterogeneity of hydrologic transport dynamics and reveal that these dynamics are primarily controlled by the hydroclimatic gradient of the aridity index across Europe. Contrasting the space-time dynamics of transport times with reactive timescales of denitrification in soil indicate that similar to 75% of the cultivated areas across Europe are potentially vulnerable to nitrate leaching for at least onethird of the year. We find that neglecting the transient nature of transport and reaction timescale results in a great underestimation of the extent of vulnerable regions by almost 50%. Therefore, future vulnerability and risk assessment studies must account for the transient behaviour of transport and biogeochemical transformation processes.
Strong hydroclimatic controls on vulnerability to subsurface nitrate contamination across Europe
(2020)
Subsurface contamination due to excessive nutrient surpluses is a persistent and widespread problem in agricultural areas across Europe. The vulnerability of a particular location to pollution from reactive solutes, such as nitrate, is determined by the interplay between hydrologic transport and biogeochemical transformations. Current studies on the controls of subsurface vulnerability do not consider the transient behaviour of transport dynamics in the root zone. Here, using state-of-the-art hydrologic simulations driven by observed hydroclimatic forcing, we demonstrate the strong spatiotemporal heterogeneity of hydrologic transport dynamics and reveal that these dynamics are primarily controlled by the hydroclimatic gradient of the aridity index across Europe. Contrasting the space-time dynamics of transport times with reactive timescales of denitrification in soil indicate that similar to 75% of the cultivated areas across Europe are potentially vulnerable to nitrate leaching for at least onethird of the year. We find that neglecting the transient nature of transport and reaction timescale results in a great underestimation of the extent of vulnerable regions by almost 50%. Therefore, future vulnerability and risk assessment studies must account for the transient behaviour of transport and biogeochemical transformation processes.
We present a new set of global and local sea‐level projections at example tide gauge locations under the RCP2.6, RCP4.5, and RCP8.5 emissions scenarios. Compared to the CMIP5‐based sea‐level projections presented in IPCC AR5, we introduce a number of methodological innovations, including (i) more comprehensive treatment of uncertainties, (ii) direct traceability between global and local projections, and (iii) exploratory extended projections to 2300 based on emulation of individual CMIP5 models. Combining the projections with observed tide gauge records, we explore the contribution to total variance that arises from sea‐level variability, different emissions scenarios, and model uncertainty. For the period out to 2300 we further breakdown the model uncertainty by sea‐level component and consider the dependence on geographic location, time horizon, and emissions scenario. Our analysis highlights the importance of local variability for sea‐level change in the coming decades and the potential value of annual‐to‐decadal predictions of local sea‐level change. Projections to 2300 show a substantial degree of committed sea‐level rise under all emissions scenarios considered and highlight the reduced future risk associated with RCP2.6 and RCP4.5 compared to RCP8.5. Tide gauge locations can show large ( > 50%) departures from the global average, in some cases even reversing the sign of the change. While uncertainty in projections of the future Antarctic ice dynamic response tends to dominate post‐2100, we see substantial differences in the breakdown of model variance as a function of location, time scale, and emissions scenario.
We present a new set of global and local sea‐level projections at example tide gauge locations under the RCP2.6, RCP4.5, and RCP8.5 emissions scenarios. Compared to the CMIP5‐based sea‐level projections presented in IPCC AR5, we introduce a number of methodological innovations, including (i) more comprehensive treatment of uncertainties, (ii) direct traceability between global and local projections, and (iii) exploratory extended projections to 2300 based on emulation of individual CMIP5 models. Combining the projections with observed tide gauge records, we explore the contribution to total variance that arises from sea‐level variability, different emissions scenarios, and model uncertainty. For the period out to 2300 we further breakdown the model uncertainty by sea‐level component and consider the dependence on geographic location, time horizon, and emissions scenario. Our analysis highlights the importance of local variability for sea‐level change in the coming decades and the potential value of annual‐to‐decadal predictions of local sea‐level change. Projections to 2300 show a substantial degree of committed sea‐level rise under all emissions scenarios considered and highlight the reduced future risk associated with RCP2.6 and RCP4.5 compared to RCP8.5. Tide gauge locations can show large ( > 50%) departures from the global average, in some cases even reversing the sign of the change. While uncertainty in projections of the future Antarctic ice dynamic response tends to dominate post‐2100, we see substantial differences in the breakdown of model variance as a function of location, time scale, and emissions scenario.
The simulation of broad-band (0.1 to 10 + Hz) ground-shaking over deep and spatially extended sedimentary basins at regional scales is challenging. We evaluate the ground-shaking of a potential M 6.5 earthquake in the southern Lower Rhine Embayment, one of the most important areas of earthquake recurrence north of the Alps, close to the city of Cologne in Germany. In a first step, information from geological investigations, seismic experiments and boreholes is combined for deriving a harmonized 3D velocity and attenuation model of the sedimentary layers. Three alternative approaches are then applied and compared to evaluate the impact of the sedimentary cover on ground-motion amplification. The first approach builds on existing response spectra ground-motion models whose amplification factors empirically take into account the influence of the sedimentary layers through a standard parameterization. In the second approach, site-specific 1D amplification functions are computed from the 3D basin model. Using a random vibration theory approach, we adjust the empirical response spectra predicted for soft rock conditions by local site amplification factors: amplifications and associated ground-motions are predicted both in the Fourier and in the response spectra domain. In the third approach, hybrid physics-based ground-motion simulations are used to predict time histories for soft rock conditions which are subsequently modified using the 1D site-specific amplification functions computed in method 2. For large distances and at short periods, the differences between the three approaches become less notable due to the significant attenuation of the sedimentary layers. At intermediate and long periods, generic empirical ground-motion models provide lower levels of amplification from sedimentary soils compared to methods taking into account site-specific 1D amplification functions. In the near-source region, hybrid physics-based ground-motions models illustrate the potentially large variability of ground-motion due to finite source effects.
Borehole leakage is a common and complex issue. Understanding the fluid flow characteristics of a cemented area inside a borehole is crucial to monitor and quantify the wellbore integrity as well as to find solutions to minimise existing leakages. In order to improve our understanding of the flow behaviour of cemented boreholes, we investigated experimental data of a large-scale borehole leakage tests by means of numerical modelling using three different conceptual models. The experiment was performed with an autoclave system consisting of two vessels bridged by a cement-filled casing. After a partial bleed-off at the well-head, a sustained casing pressure was observed due to fluid flow through the cementsteel composite. The aim of our simulations is to investigate and quantify the permeability of the cement-steel composite. From our model results, we conclude that the flow occurred along a preferential flow path at the cement-steel interface. Thus, the inner part of the cement core was impermeable during the duration of the experiment. The preferential flow path can be described as a highly permeable and highly porous area with an aperture of about 5 mu m and a permeability of 3 . 10(-12) m(2) (3 Darcy). It follows that the fluid flow characteristics of a cemented area inside a borehole cannot be described using one permeability value for the entire cement-steel composite. Furthermore, it can be concluded that the quality of the cement and the filling process regarding the cement-steel interface is crucial to minimize possible well leakages.
The southern Central Andes (SCA) (between 27 degrees S and 40 degrees S) is bordered to the west by the convergent margin between the continental South American Plate and the oceanic Nazca Plate. The subduction angle along this margin is variable, as is the deformation of the upper plate. Between 33 degrees S and 35 degrees S, the subduction angle of the Nazca plate increases from sub-horizontal (< 5 degrees) in the north to relatively steep (similar to 30 degrees) in the south. The SCA contain inherited lithological and structural heterogeneities within the crust that have been reactivated and overprinted since the onset of subduction and associated Cenozoic deformation within the Andean orogen. The distribution of the deformation within the SCA has often been attributed to the variations in the subduction angle and the reactivation of these inherited heterogeneities. However, the possible influence that the thickness and composition of the continental crust have had on both short-term and long-term deformation of the SCA is yet to be thoroughly investigated. For our investigations, we have derived density distributions and thicknesses for various layers that make up the lithosphere and evaluated their relationships with tectonic events that occurred over the history of the Andean orogeny and, in particular, investigated the short- and long-term nature of the present-day deformation processes. We established a 3D model of lithosphere beneath the orogen and its foreland (29 degrees S-39 degrees S) that is consistent with currently available geological and geophysical data, including the gravity data. The modelled crustal configuration and density distribution reveal spatial relationships with different tectonic domains: the crystalline crust in the orogen (the magmatic arc and the main orogenic wedge) is thicker (similar to 55 km) and less dense (similar to 2900 kg/m(3)) than in the forearc (similar to 35 km, similar to 2975 kg/m(3)) and foreland (similar to 30 km, similar to 3000 kg/m(3)). Crustal thickening in the orogen probably occurred as a result of stacking of low-density domains, while density and thickness variations beneath the forearc and foreland most likely reflect differences in the tectonic evolution of each area following crustal accretion. No clear spatial relationship exists between the density distribution within the lithosphere and previously proposed boundaries of crustal terranes accreted during the early Paleozoic. Areas with ongoing deformation show a spatial correlation with those areas that have the highest topographic gradients and where there are abrupt changes in the average crustal-density contrast. This suggests that the short-term deformation within the interior of the Andean orogen and its foreland is fundamentally influenced by the crustal composition and the relative thickness of different crustal layers. A thicker, denser, and potentially stronger lithosphere beneath the northern part of the SCA foreland is interpreted to have favoured a strong coupling between the Nazca and South American plates, facilitating the development of a sub-horizontal slab.
Arctic lakes located in permafrost regions are susceptible to catastrophic drainage. In this study, we reconstructed historical lake drainage events on the western Arctic Coastal Plain of Alaska between 1955 and 2017 using USGS topographic maps, historical aerial photography (1955), and Landsat Imagery (ca. 1975, ca. 2000, and annually since 2000). We identified 98 lakes larger than 10 ha that partially (>25% of area) or completely drained during the 62-year period. Decadal-scale lake drainage rates progressively declined from 2.0 lakes/yr (1955-1975), to 1.6 lakes/yr (1975-2000), and to 1.2 lakes/yr (2000-2017) in the ~30,000-km(2) study area. Detailed Landsat trend analysis between 2000 and 2017 identified two years, 2004 and 2006, with a cluster (five or more) of lake drainages probably associated with bank overtopping or headward erosion. To identify future potential lake drainages, we combined the historical lake drainage observations with a geospatial dataset describing lake elevation, hydrologic connectivity, and adjacent lake margin topographic gradients developed with a 5-m-resolution digital surface model. We identified ~1900 lakes likely to be prone to drainage in the future. Of the 20 lakes that drained in the most recent study period, 85% were identified in this future lake drainage potential dataset. Our assessment of historical lake drainage magnitude, mechanisms and pathways, and identification of potential future lake drainages provides insights into how arctic lowland landscapes may change and evolve in the coming decades to centuries.
Arctic lakes located in permafrost regions are susceptible to catastrophic drainage. In this study, we reconstructed historical lake drainage events on the western Arctic Coastal Plain of Alaska between 1955 and 2017 using USGS topographic maps, historical aerial photography (1955), and Landsat Imagery (ca. 1975, ca. 2000, and annually since 2000). We identified 98 lakes larger than 10 ha that partially (>25% of area) or completely drained during the 62-year period. Decadal-scale lake drainage rates progressively declined from 2.0 lakes/yr (1955-1975), to 1.6 lakes/yr (1975-2000), and to 1.2 lakes/yr (2000-2017) in the ~30,000-km(2) study area. Detailed Landsat trend analysis between 2000 and 2017 identified two years, 2004 and 2006, with a cluster (five or more) of lake drainages probably associated with bank overtopping or headward erosion. To identify future potential lake drainages, we combined the historical lake drainage observations with a geospatial dataset describing lake elevation, hydrologic connectivity, and adjacent lake margin topographic gradients developed with a 5-m-resolution digital surface model. We identified ~1900 lakes likely to be prone to drainage in the future. Of the 20 lakes that drained in the most recent study period, 85% were identified in this future lake drainage potential dataset. Our assessment of historical lake drainage magnitude, mechanisms and pathways, and identification of potential future lake drainages provides insights into how arctic lowland landscapes may change and evolve in the coming decades to centuries.
Most hydrological studies rely on a model calibrated using discharge alone. However, judging the model reliability based on such calibration is problematic, as it does not guarantee the correct representation of internal hydrological processes. This study aims (a) to develop a comprehensive multi-objective calibration framework using remote sensing vegetation data and hydrological signatures (flow duration curve - FDC, and baseflow index) in addition to discharge, and (b) to apply this framework for calibration of the Soil and Water Assessment Tool (SWAT) in a typical Andean catchment. Overall, our calibration approach outperformed traditional discharge-based and FDC signature-based calibration strategies in terms of vegetation, streamflow, and flow partitioning simulation. New hydrological insights for the region are the following: baseflow is the main component of the streamflow sustaining the long dry-season flow, and pasture areas offer higher water yield and baseflow than other land-cover types. The proposed approach could be used in other data-scarce regions with complex topography.
Land degradation and water availability in semi-arid regions are interdependent challenges for management that are influenced by climatic and anthropogenic changes. Erosion and high sediment loads in rivers cause reservoir siltation and decrease storage capacity, which pose risk on water security for citizens, agriculture, and industry. In regions where resources for management are limited, identifying spatial-temporal variability of sediment sources is crucial to decrease siltation. Despite widespread availability of rigorous methods, approaches simplifying spatial and temporal variability of erosion are often inappropriately applied to very data sparse semi-arid regions. In this work, we review existing approaches for mapping erosional hotspots, and provide an example of spatial-temporal mapping approach in two case study regions. The barriers limiting data availability and their effects on erosion mapping methods, their validation, and resulting prioritization of leverage management areas are discussed.
Land degradation and water availability in semi-arid regions are interdependent challenges for management that are influenced by climatic and anthropogenic changes. Erosion and high sediment loads in rivers cause reservoir siltation and decrease storage capacity, which pose risk on water security for citizens, agriculture, and industry. In regions where resources for management are limited, identifying spatial-temporal variability of sediment sources is crucial to decrease siltation. Despite widespread availability of rigorous methods, approaches simplifying spatial and temporal variability of erosion are often inappropriately applied to very data sparse semi-arid regions. In this work, we review existing approaches for mapping erosional hotspots, and provide an example of spatial-temporal mapping approach in two case study regions. The barriers limiting data availability and their effects on erosion mapping methods, their validation, and resulting prioritization of leverage management areas are discussed.
Understanding the hydrologic connectivity between kettle holes and shallow groundwater, particularly in reaction to the highly variable local meteorological conditions, is of paramount importance for tracing water in a hydro(geo)logically complex landscape and thus for integrated water resource management. This article is aimed at identifying the dominant hydrological processes affecting the kettle holes' water balance and their interactions with the shallow groundwater domain in the Uckermark region, located in the north-east of Germany. For this reason, based on the stable isotopes of oxygen (delta O-18) and hydrogen (delta H-2), an isotopic mass balance model was employed to compute the evaporative loss of water from the kettle holes from February to August 2017. Results demonstrated that shallow groundwater inflow may play the pivotal role in the processes taking part in the hydrology of the kettle holes in the Uckermark region. Based on the calculated evaporation/inflow (E/I) ratios, most of the kettle holes (86.7%) were ascertained to have a partially open, flow-through-dominated system. Moreover, we identified an inverse correlation between E/I ratios and the altitudes of the kettle holes. The same holds for electrical conductivity (EC) and the altitudes of the kettle holes. In accordance with the findings obtained from this study, a conceptual model explaining the interaction between the shallow groundwater and the kettle holes of Uckermark was developed. The model exhibited that across the highest altitudes, the recharge kettle holes are dominant, where a lower ratio of E/I and a lower EC was detected. By contrast, the lowest topographical depressions represent the discharge kettle holes, where a higher ratio of E/I and EC could be identified. The kettle holes existing in between were categorized as flow-through kettle holes through which the recharge takes place from one side and discharge from the other side.
Understanding the hydrologic connectivity between kettle holes and shallow groundwater, particularly in reaction to the highly variable local meteorological conditions, is of paramount importance for tracing water in a hydro(geo)logically complex landscape and thus for integrated water resource management. This article is aimed at identifying the dominant hydrological processes affecting the kettle holes' water balance and their interactions with the shallow groundwater domain in the Uckermark region, located in the north-east of Germany. For this reason, based on the stable isotopes of oxygen (delta O-18) and hydrogen (delta H-2), an isotopic mass balance model was employed to compute the evaporative loss of water from the kettle holes from February to August 2017. Results demonstrated that shallow groundwater inflow may play the pivotal role in the processes taking part in the hydrology of the kettle holes in the Uckermark region. Based on the calculated evaporation/inflow (E/I) ratios, most of the kettle holes (86.7%) were ascertained to have a partially open, flow-through-dominated system. Moreover, we identified an inverse correlation between E/I ratios and the altitudes of the kettle holes. The same holds for electrical conductivity (EC) and the altitudes of the kettle holes. In accordance with the findings obtained from this study, a conceptual model explaining the interaction between the shallow groundwater and the kettle holes of Uckermark was developed. The model exhibited that across the highest altitudes, the recharge kettle holes are dominant, where a lower ratio of E/I and a lower EC was detected. By contrast, the lowest topographical depressions represent the discharge kettle holes, where a higher ratio of E/I and EC could be identified. The kettle holes existing in between were categorized as flow-through kettle holes through which the recharge takes place from one side and discharge from the other side.
An exceptionally strong stationary planetary wave with Zonal Wavenumber 1 led to a sudden stratospheric warming (SSW) in the Southern Hemisphere in September 2019. Ionospheric data from European Space Agency's Swarm satellite constellation mission show prominent 6-day variations in the dayside low-latitude region at this time, which can be attributed to forcing from the middle atmosphere by the Rossby normal mode "quasi-6-day wave" (Q6DW). Geopotential height measurements by the Microwave Limb Sounder aboard National Aeronautics and Space Administration's Aura satellite reveal a burst of global Q6DW activity in the mesosphere and lower thermosphere during the SSW, which is one of the strongest in the record. The Q6DW is apparently generated in the polar stratosphere at 30-40 km, where the atmosphere is unstable due to strong vertical wind shear connected with planetary wave breaking. These results suggest that an Antarctic SSW can lead to ionospheric variability through wave forcing from the middle atmosphere.
Plain Language Summary: A sudden stratospheric warming (SSW) is an extreme wintertime polar meteorological phenomenon occurring mostly over the Arctic region. Studies have shown that Arctic SSW can influence the entire atmosphere. In September 2019, a rare SSW event occurred in the Antarctic region, providing an opportunity to investigate its broader impact on the whole atmosphere. We present observations from the middle atmosphere and ionosphere during this event, noting unusually strong wave activity throughout this region. Our results suggest that an Antarctic SSW can have a significant impact on the whole atmosphere system similar to those due to Arctic events.
An exceptionally strong stationary planetary wave with Zonal Wavenumber 1 led to a sudden stratospheric warming (SSW) in the Southern Hemisphere in September 2019. Ionospheric data from European Space Agency's Swarm satellite constellation mission show prominent 6-day variations in the dayside low-latitude region at this time, which can be attributed to forcing from the middle atmosphere by the Rossby normal mode "quasi-6-day wave" (Q6DW). Geopotential height measurements by the Microwave Limb Sounder aboard National Aeronautics and Space Administration's Aura satellite reveal a burst of global Q6DW activity in the mesosphere and lower thermosphere during the SSW, which is one of the strongest in the record. The Q6DW is apparently generated in the polar stratosphere at 30-40 km, where the atmosphere is unstable due to strong vertical wind shear connected with planetary wave breaking. These results suggest that an Antarctic SSW can lead to ionospheric variability through wave forcing from the middle atmosphere.
Plain Language Summary: A sudden stratospheric warming (SSW) is an extreme wintertime polar meteorological phenomenon occurring mostly over the Arctic region. Studies have shown that Arctic SSW can influence the entire atmosphere. In September 2019, a rare SSW event occurred in the Antarctic region, providing an opportunity to investigate its broader impact on the whole atmosphere. We present observations from the middle atmosphere and ionosphere during this event, noting unusually strong wave activity throughout this region. Our results suggest that an Antarctic SSW can have a significant impact on the whole atmosphere system similar to those due to Arctic events.
Previous studies have explored the consequences of flood events for exposed households and companies by focusing on single flood events. Less is known about the consequences of experiencing repeated flood events for the resilience of households and companies. In this paper, we therefore explore how multiple floods experience affects the resilience of exposed households and companies. Resilience was made operational through individual appraisals of households and companies' ability to withstand and recover from material as well as health and psychological impacts of the 2013 flood in Germany. The paper is based on three different datasets including more than 2000 households and 300 companies that were affected by the 2013 flood. The surveys revealed that the resilience of households seems to increase, but only with regard to their subjectively appraised ability to withstand impacts on mobile goods and equipment (e.g., cars, TV, and radios). In regard to the ability of households to withstand overall financial consequences of repetitive floods, evidence for nonlinear (quadratic) trends can be found. With regard to psychological and health-related consequences, the findings are mixed but provide tentative evidence for eroding resilience among households. Companies' resilience increased with respect to material assets but appears to decrease with respect to ability to recover. We conclude by arguing that clear and operational definitions of resilience are required so that evidence-based resilience baselines can be established to assess whether resilience is eroding or improving over time.
Previous studies have explored the consequences of flood events for exposed households and companies by focusing on single flood events. Less is known about the consequences of experiencing repeated flood events for the resilience of households and companies. In this paper, we therefore explore how multiple floods experience affects the resilience of exposed households and companies. Resilience was made operational through individual appraisals of households and companies' ability to withstand and recover from material as well as health and psychological impacts of the 2013 flood in Germany. The paper is based on three different datasets including more than 2000 households and 300 companies that were affected by the 2013 flood. The surveys revealed that the resilience of households seems to increase, but only with regard to their subjectively appraised ability to withstand impacts on mobile goods and equipment (e.g., cars, TV, and radios). In regard to the ability of households to withstand overall financial consequences of repetitive floods, evidence for nonlinear (quadratic) trends can be found. With regard to psychological and health-related consequences, the findings are mixed but provide tentative evidence for eroding resilience among households. Companies' resilience increased with respect to material assets but appears to decrease with respect to ability to recover. We conclude by arguing that clear and operational definitions of resilience are required so that evidence-based resilience baselines can be established to assess whether resilience is eroding or improving over time.
The hyperthermal events of the Cenozoic, including the Paleocene-Eocene Thermal Maximum, provide an opportunity to investigate the potential effects of climate warming on marine ecosystems. Here, we examine the shallow benthic marine communities preserved in the late Cretaceous to Eocene strata on the Gulf Coastal Plain (United States). In stark contrast to the ecological shifts following the end-Cretaceous mass extinction, our data show that the early Cenozoic hyperthermals did not have a long-term impact on the generic diversity nor composition of the Gulf Coastal Plain molluscan communities. We propose that these communities were resilient to climate change because molluscs are better adapted to high temperatures than other taxa, as demonstrated by their physiology and evolutionary history. In terms of resilience, these communities differ from other shallow-water carbonate ecosystems, such as reef communities, which record significant changes during the early Cenozoic hyperthermals. These data highlight the strikingly different responses of community types, i.e., the almost imperceptible response of molluscs versus the marked turnover of foraminifera and reef faunas. The impact on molluscan communities may have been low because detrimental conditions did not devastate the entire Gulf Coastal Plain, allowing molluscs to rapidly recolonise vacated areas once harsh environmental conditions ameliorated.
The hyperthermal events of the Cenozoic, including the Paleocene-Eocene Thermal Maximum, provide an opportunity to investigate the potential effects of climate warming on marine ecosystems. Here, we examine the shallow benthic marine communities preserved in the late Cretaceous to Eocene strata on the Gulf Coastal Plain (United States). In stark contrast to the ecological shifts following the end-Cretaceous mass extinction, our data show that the early Cenozoic hyperthermals did not have a long-term impact on the generic diversity nor composition of the Gulf Coastal Plain molluscan communities. We propose that these communities were resilient to climate change because molluscs are better adapted to high temperatures than other taxa, as demonstrated by their physiology and evolutionary history. In terms of resilience, these communities differ from other shallow-water carbonate ecosystems, such as reef communities, which record significant changes during the early Cenozoic hyperthermals. These data highlight the strikingly different responses of community types, i.e., the almost imperceptible response of molluscs versus the marked turnover of foraminifera and reef faunas. The impact on molluscan communities may have been low because detrimental conditions did not devastate the entire Gulf Coastal Plain, allowing molluscs to rapidly recolonise vacated areas once harsh environmental conditions ameliorated.
Groundwater is the biggest single source of high-quality freshwater worldwide, which is also continuously threatened by the changing climate. In this paper, we investigate the response of the regional groundwater system to climate change under three global warming levels (1.5, 2, and 3 ∘C) in a central German basin (Nägelstedt). This investigation is conducted by deploying an integrated modeling workflow that consists of a mesoscale hydrologic model (mHM) and a fully distributed groundwater model, OpenGeoSys (OGS). mHM is forced with climate simulations of five general circulation models under three representative concentration pathways. The diffuse recharges estimated by mHM are used as boundary forcings to the OGS groundwater model to compute changes in groundwater levels and travel time distributions. Simulation results indicate that groundwater recharges and levels are expected to increase slightly under future climate scenarios. Meanwhile, the mean travel time is expected to decrease compared to the historical average. However, the ensemble simulations do not all agree on the sign of relative change. Changes in mean travel time exhibit a larger variability than those in groundwater levels. The ensemble simulations do not show a systematic relationship between the projected change (in both groundwater levels and travel times) and the warming level, but they indicate an increased variability in projected changes with adjusting the enhanced warming level from 1.5 to 3 ∘C. Correspondingly, it is highly recommended to restrain the trend of global warming.
Groundwater is the biggest single source of high-quality freshwater worldwide, which is also continuously threatened by the changing climate. In this paper, we investigate the response of the regional groundwater system to climate change under three global warming levels (1.5, 2, and 3 ∘C) in a central German basin (Nägelstedt). This investigation is conducted by deploying an integrated modeling workflow that consists of a mesoscale hydrologic model (mHM) and a fully distributed groundwater model, OpenGeoSys (OGS). mHM is forced with climate simulations of five general circulation models under three representative concentration pathways. The diffuse recharges estimated by mHM are used as boundary forcings to the OGS groundwater model to compute changes in groundwater levels and travel time distributions. Simulation results indicate that groundwater recharges and levels are expected to increase slightly under future climate scenarios. Meanwhile, the mean travel time is expected to decrease compared to the historical average. However, the ensemble simulations do not all agree on the sign of relative change. Changes in mean travel time exhibit a larger variability than those in groundwater levels. The ensemble simulations do not show a systematic relationship between the projected change (in both groundwater levels and travel times) and the warming level, but they indicate an increased variability in projected changes with adjusting the enhanced warming level from 1.5 to 3 ∘C. Correspondingly, it is highly recommended to restrain the trend of global warming.
Recent climatic changes have the potential to severely alter river runoff, particularly in snow-dominated river basins. Effects of changing snow covers superimpose with changes in precipitation and anthropogenic modifications of the watershed and river network. In the attempt to identify and disentangle long-term effects of different mechanisms, we employ a set of analytical tools to extract long-term changes in river runoff at high resolution. We combine quantile sampling with moving average trend statistics and empirical mode decomposition and apply these tools to discharge data recorded along rivers with nival, pluvial and mixed flow regimes as well as temperature and precipitation data covering the time frame 1869-2016. With a focus on central Europe, we analyse the long-term impact of snow cover and precipitation changes along with their interaction with reservoir constructions.
Our results show that runoff seasonality of snow-dominated rivers decreases. Runoff increases in winter and spring, while discharge decreases in summer and at the beginning of autumn. We attribute this redistribution of annual flow mainly to reservoir constructions in the Alpine ridge. During the course of the last century, large fractions of the Alpine rivers were dammed to produce hydropower. In recent decades, runoff changes induced by reservoir constructions seem to overlap with changes in snow cover. We suggest that Alpine signals propagate downstream and affect runoff far outside the Alpine area in river segments with mixed flow regimes. Furthermore, our results hint at more (intense) rain-fall in recent decades. Detected increases in high discharge can be traced back to corresponding changes in precipitation.
Submerged sequences of marine terraces potentially provide crucial information of past sea-level positions. However, the distribution and characteristics of drowned marine terrace sequences are poorly known at a global scale. Using bathymetric data and novel mapping and modeling techniques, we studied a submerged sequence of marine terraces in the Bay of Biscay with the objective to identify the distribution and morphologies of submerged marine terraces and the timing and conditions that allowed their formation and preservation. To accomplish the objectives a high-resolution bathymetry (5 m) was analyzed using Geographic Information Systems and TerraceM(R). The successive submerged terraces were identified using a Surface Classification Model, which linearly combines the slope and the roughness of the surface to extract fossil sea-cliffs and fossil rocky shore platforms. For that purpose, contour and hillshaded maps were also analyzed. Then, shoreline angles, a geomorphic marker located at the intersection between the fossil sea-cliff and platform, were mapped analyzing swath profiles perpendicular to the isobaths. Most of the submerged strandlines are irregularly preserved throughout the continental shelf. In summary, 12 submerged terraces with their shoreline angles between approximately: -13 m (T1), -30 and -32 m (T2), -34 and 41 m (T3), -44 and -47 m (T4), -49 and 53 m (T5), -55 and 58 m (T6), -59 and 62 m (T7), -65 and 67 m (T8), -68 and 70 m (T9), -74 and -77 m (T10), -83 and -86 m (T11) and -89 and 92 m (T12). Nevertheless, the ones showing the best lateral continuity and preservation in the central part of the shelf are T3, T4, T5, T7, T8, and T10. The age of the terraces has been estimated using a landscape evolution model. To simulate the formation and preservation of submerged terraces three different scenarios: (i) 20-0 ka; (ii) 128-0 ka; and (iii) 128-20 ka, were compared. The best scenario for terrace generation was between 128 and 20 Ka, where T3, T5, and T7 could have been formed.
Submerged sequences of marine terraces potentially provide crucial information of past sea-level positions. However, the distribution and characteristics of drowned marine terrace sequences are poorly known at a global scale. Using bathymetric data and novel mapping and modeling techniques, we studied a submerged sequence of marine terraces in the Bay of Biscay with the objective to identify the distribution and morphologies of submerged marine terraces and the timing and conditions that allowed their formation and preservation. To accomplish the objectives a high-resolution bathymetry (5 m) was analyzed using Geographic Information Systems and TerraceM(R). The successive submerged terraces were identified using a Surface Classification Model, which linearly combines the slope and the roughness of the surface to extract fossil sea-cliffs and fossil rocky shore platforms. For that purpose, contour and hillshaded maps were also analyzed. Then, shoreline angles, a geomorphic marker located at the intersection between the fossil sea-cliff and platform, were mapped analyzing swath profiles perpendicular to the isobaths. Most of the submerged strandlines are irregularly preserved throughout the continental shelf. In summary, 12 submerged terraces with their shoreline angles between approximately: -13 m (T1), -30 and -32 m (T2), -34 and 41 m (T3), -44 and -47 m (T4), -49 and 53 m (T5), -55 and 58 m (T6), -59 and 62 m (T7), -65 and 67 m (T8), -68 and 70 m (T9), -74 and -77 m (T10), -83 and -86 m (T11) and -89 and 92 m (T12). Nevertheless, the ones showing the best lateral continuity and preservation in the central part of the shelf are T3, T4, T5, T7, T8, and T10. The age of the terraces has been estimated using a landscape evolution model. To simulate the formation and preservation of submerged terraces three different scenarios: (i) 20-0 ka; (ii) 128-0 ka; and (iii) 128-20 ka, were compared. The best scenario for terrace generation was between 128 and 20 Ka, where T3, T5, and T7 could have been formed.
The origin of Asian monsoons
(2020)
The Cenozoic inception and development of the Asian monsoons remain unclear and have generated much debate, as several hypotheses regarding circulation patterns at work in Asia during the Eocene have been proposed in the few last decades. These include (a) the existence of modern-like monsoons since the early Eocene; (b) that of a weak South Asian monsoon (SAM) and little to no East Asian monsoon (EAM); or (c) a prevalence of the Intertropical Convergence Zone (ITCZ) migrations, also referred to as Indonesian-Australian monsoon (I-AM). As SAM and EAM are supposed to have been triggered or enhanced primarily by Asian palaeogeographic changes, their possible inception in the very dynamic Eocene palaeogeographic context remains an open question, both in the modelling and field-based communities. We investigate here Eocene Asian climate conditions using the IPSL-CM5A2 (Sepulchre et al., 2019) earth system model and revised palaeogeographies. Our Eocene climate simulation yields atmospheric circulation patterns in Asia substantially different from modern conditions. A large high-pressure area is simulated over the Tethys ocean, which generates intense low tropospheric winds blowing southward along the western flank of the proto-Himalayan-Tibetan plateau (HTP) system. This low-level wind system blocks, to latitudes lower than 10 degrees N, the migration of humid and warm air masses coming from the Indian Ocean. This strongly contrasts with the modern SAM, during which equatorial air masses reach a latitude of 20-25 degrees N over India and southeastern China. Another specific feature of our Eocene simulation is the widespread subsidence taking place over northern India in the midtroposphere (around 5000 m), preventing deep convective updraught that would transport water vapour up to the condensation level. Both processes lead to the onset of a broad arid region located over northern India and over the HTP. More humid regions of high seasonality in precipitation encircle this arid area, due to the prevalence of the Intertropical Convergence Zone (ITCZ) migrations (or Indonesian-Australian monsoon, I-AM) rather than monsoons. Although the existence of this central arid region may partly result from the specifics of our simulation (model dependence and palaeogeographic uncertainties) and has yet to be confirmed by proxy records, most of the observational evidence for Eocene monsoons are located in the highly seasonal transition zone between the arid area and the more humid surroundings. We thus suggest that a zonal arid climate prevailed over Asia before the initiation of monsoons that most likely occurred following Eocene palaeogeographic changes. Our results also show that precipitation seasonality should be used with caution to infer the presence of a monsoonal circulation and that the collection of new data in this arid area is of paramount importance to allow the debate to move forward.
The origin of Asian monsoons
(2020)
The Cenozoic inception and development of the Asian monsoons remain unclear and have generated much debate, as several hypotheses regarding circulation patterns at work in Asia during the Eocene have been proposed in the few last decades. These include (a) the existence of modern-like monsoons since the early Eocene; (b) that of a weak South Asian monsoon (SAM) and little to no East Asian monsoon (EAM); or (c) a prevalence of the Intertropical Convergence Zone (ITCZ) migrations, also referred to as Indonesian-Australian monsoon (I-AM). As SAM and EAM are supposed to have been triggered or enhanced primarily by Asian palaeogeographic changes, their possible inception in the very dynamic Eocene palaeogeographic context remains an open question, both in the modelling and field-based communities. We investigate here Eocene Asian climate conditions using the IPSL-CM5A2 (Sepulchre et al., 2019) earth system model and revised palaeogeographies. Our Eocene climate simulation yields atmospheric circulation patterns in Asia substantially different from modern conditions. A large high-pressure area is simulated over the Tethys ocean, which generates intense low tropospheric winds blowing southward along the western flank of the proto-Himalayan-Tibetan plateau (HTP) system. This low-level wind system blocks, to latitudes lower than 10 degrees N, the migration of humid and warm air masses coming from the Indian Ocean. This strongly contrasts with the modern SAM, during which equatorial air masses reach a latitude of 20-25 degrees N over India and southeastern China. Another specific feature of our Eocene simulation is the widespread subsidence taking place over northern India in the midtroposphere (around 5000 m), preventing deep convective updraught that would transport water vapour up to the condensation level. Both processes lead to the onset of a broad arid region located over northern India and over the HTP. More humid regions of high seasonality in precipitation encircle this arid area, due to the prevalence of the Intertropical Convergence Zone (ITCZ) migrations (or Indonesian-Australian monsoon, I-AM) rather than monsoons. Although the existence of this central arid region may partly result from the specifics of our simulation (model dependence and palaeogeographic uncertainties) and has yet to be confirmed by proxy records, most of the observational evidence for Eocene monsoons are located in the highly seasonal transition zone between the arid area and the more humid surroundings. We thus suggest that a zonal arid climate prevailed over Asia before the initiation of monsoons that most likely occurred following Eocene palaeogeographic changes. Our results also show that precipitation seasonality should be used with caution to infer the presence of a monsoonal circulation and that the collection of new data in this arid area is of paramount importance to allow the debate to move forward.
Groundwater levels are monitored by environmental agencies to support the sustainable use of groundwater resources. For this purpose continuous and spatially comprehensive monitoring in high spatial and temporal resolution is desired. This leads to large datasets that have to be checked for quality and analysed to distinguish local anthropogenic influences from natural variability of the groundwater level dynamics at each well. Both technical problems with the measurements as well as local anthropogenic influences can lead to local anomalies in the hydrographs. We suggest a fast and efficient screening method for the identification of well-specific peculiarities in hydrographs of groundwater head monitoring networks. The only information required is a set of time series of groundwater heads all measured at the same instants of time. For each well of the monitoring network a reference hydrograph is calculated, describing expected “normal” behaviour at the respective well as is typical for the monitored region. The reference hydrograph is calculated by multiple linear regression of the observed hydrograph with the “stable” principal components (PCs) of a principal component analysis of all groundwater head series of the network as predictor variables. The stable PCs are those PCs which were found in a random subsampling procedure to be rather insensitive to the specific selection of the analysed observation wells, i.e. complete series, and to the specific selection of measurement dates. Hence they can be considered to be representative for the monitored region in the respective period. The residuals of the reference hydrograph describe local deviations from the normal behaviour. Peculiarities in the residuals allow the data to be checked for measurement errors and the wells with a possible anthropogenic influence to be identified. The approach was tested with 141 groundwater head time series from the state authority groundwater monitoring network in northeastern Germany covering the period from 1993 to 2013 at an approximately weekly frequency of measurement.
According to established understanding, deep-water formation in the North Atlantic and Southern Ocean keeps the deep ocean cold, counter-acting the downward mixing of heat from the warmer surface waters in the bulk of the world ocean. Therefore, periods of strong Atlantic meridional overturning circulation (AMOC) are expected to coincide with cooling of the deep ocean and warming of the surface waters. It has recently been proposed that this relation may have reversed due to global warming, and that during the past decades a strong AMOC coincides with warming of the deep ocean and relative cooling of the surface, by transporting increasingly warmer waters downward. Here we present multiple lines of evidence, including a statistical evaluation of the observed global mean temperature, ocean heat content, and different AMOC proxies, that lead to the opposite conclusion: even during the current ongoing global temperature rise a strong AMOC warms the surface. The observed weakening of the AMOC has therefore delayed global surface warming rather than enhancing it.
Social Media Abstract:
The overturning circulation in the Atlantic Ocean has weakened in response to global warming, as predicted by climate models. Since it plays an important role in transporting heat, nutrients and carbon, a slowdown will affect global climate processes and the global mean temperature. Scientists have questioned whether this slowdown has worked to cool or warm global surface temperatures. This study analyses the overturning strength and global mean temperature evolution of the past decades and shows that a slowdown acts to reduce the global mean temperature. This is because a slower overturning means less water sinks into the deep ocean in the subpolar North Atlantic. As the surface waters are cold there, the sinking normally cools the deep ocean and thereby indirectly warms the surface, thus less sinking implies less surface warming and has a cooling effect. For the foreseeable future, this means that the slowing of the overturning will likely continue to slightly reduce the effect of the general warming due to increasing greenhouse gas concentrations.
Pollen records from Siberia are mostly absent in global or Northern Hemisphere synthesis works. Here we present a taxonomically harmonized and temporally standardized pollen dataset that was synthesized using 173 palynological records from Siberia and adjacent areas (northeastern Asia, 42-75 degrees N, 50-180 degrees E). Pollen data were taxonomically harmonized, i.e. the original 437 taxa were assigned to 106 combined pollen taxa. Age-depth models for all records were revised by applying a constant Bayesian age-depth modelling routine. The pollen dataset is available as count data and percentage data in a table format (taxa vs. samples), with age information for each sample. The dataset has relatively few sites covering the last glacial period between 40 and 11.5 ka (calibrated thousands of years before 1950 CE) particularly from the central and western part of the study area. In the Holocene period, the dataset has many sites from most of the area, with the exception of the central part of Siberia. Of the 173 pollen records, 81 % of pollen counts were downloaded from open databases (GPD, EPD, PANGAEA) and 10 % were contributions by the original data gatherers, while a few were digitized from publications. Most of the pollen records originate from peatlands (48 %) and lake sediments (33 %). Most of the records (83 %) have >= 3 dates, allowing the establishment of reliable chronologies. The dataset can be used for various purposes, including pollen data mapping (example maps for Larix at selected time slices are shown) as well as quantitative climate and vegetation reconstructions. The datasets for pollen counts and pollen percentages are available at https://doi.org/10.1594/PANGAEA.898616 (Cao et al., 2019a), also including the site information, data source, original publication, dating data, and the plant functional type for each pollen taxa.
Pollen records from Siberia are mostly absent in global or Northern Hemisphere synthesis works. Here we present a taxonomically harmonized and temporally standardized pollen dataset that was synthesized using 173 palynological records from Siberia and adjacent areas (northeastern Asia, 42-75 degrees N, 50-180 degrees E). Pollen data were taxonomically harmonized, i.e. the original 437 taxa were assigned to 106 combined pollen taxa. Age-depth models for all records were revised by applying a constant Bayesian age-depth modelling routine. The pollen dataset is available as count data and percentage data in a table format (taxa vs. samples), with age information for each sample. The dataset has relatively few sites covering the last glacial period between 40 and 11.5 ka (calibrated thousands of years before 1950 CE) particularly from the central and western part of the study area. In the Holocene period, the dataset has many sites from most of the area, with the exception of the central part of Siberia. Of the 173 pollen records, 81 % of pollen counts were downloaded from open databases (GPD, EPD, PANGAEA) and 10 % were contributions by the original data gatherers, while a few were digitized from publications. Most of the pollen records originate from peatlands (48 %) and lake sediments (33 %). Most of the records (83 %) have >= 3 dates, allowing the establishment of reliable chronologies. The dataset can be used for various purposes, including pollen data mapping (example maps for Larix at selected time slices are shown) as well as quantitative climate and vegetation reconstructions. The datasets for pollen counts and pollen percentages are available at https://doi.org/10.1594/PANGAEA.898616 (Cao et al., 2019a), also including the site information, data source, original publication, dating data, and the plant functional type for each pollen taxa.
New Zealand's Alpine Fault is a large, platebounding strike-slip fault, which ruptures in large (M-w > 8) earthquakes. We conducted field and laboratory analyses of fault rocks to assess its fault zone architecture. Results reveal that the Alpine Fault Zone has a complex geometry, comprising an anastomosing network of multiple slip planes that have accommodated different amounts of displacement. This contrasts with the previous perception of the Alpine Fault Zone, which assumes a single principal slip zone accommodated all displacement. This interpretation is supported by results of drilling projects and geophysical investigations. Furthermore, observations presented here show that the young, largely unconsolidated sediments that constitute the footwall at shallow depths have a significant influence on fault gouge rheological properties and structure.
New Zealand's Alpine Fault is a large, platebounding strike-slip fault, which ruptures in large (M-w > 8) earthquakes. We conducted field and laboratory analyses of fault rocks to assess its fault zone architecture. Results reveal that the Alpine Fault Zone has a complex geometry, comprising an anastomosing network of multiple slip planes that have accommodated different amounts of displacement. This contrasts with the previous perception of the Alpine Fault Zone, which assumes a single principal slip zone accommodated all displacement. This interpretation is supported by results of drilling projects and geophysical investigations. Furthermore, observations presented here show that the young, largely unconsolidated sediments that constitute the footwall at shallow depths have a significant influence on fault gouge rheological properties and structure.
Landscapes in high northern latitudes are assumed to be highly sensitive to future global change, but the rates and long-term trajectories of changes are rather uncertain. In the boreal zone, fires are an important factor in climate-vegetation interactions and biogeochemical cycles. Fire regimes are characterized by small, frequent, low-intensity fires within summergreen boreal forests dominated by larch, whereas evergreen boreal forests dominated by spruce and pine burn large areas less frequently but at higher intensities. Here, we explore the potential of the monosaccharide anhydrides (MA) levoglucosan, mannosan and galactosan to serve as proxies of low-intensity biomass burning in glacial-to-interglacial lake sediments from the high northern latitudes. We use sediments from Lake El'gygytgyn (cores PG 1351 and ICDP 5011-1), located in the far north-east of Russia, and study glacial and interglacial samples of the last 430 kyr (marine isotope stages 5e, 6, 7e, 8, 11c and 12) that had different climate and biome configurations. Combined with pollen and non-pollen palynomorph records from the same samples, we assess how far the modern relationships between fire, climate and vegetation persisted during the past, on orbital to centennial timescales. We find that MAs attached to particulates were well-preserved in up to 430 kyr old sediments with higher influxes from low-intensity biomass burning in interglacials compared to glacials. MA influxes significantly increase when summergreen boreal forest spreads closer to the lake, whereas they decrease when tundra-steppe environments and, especially, Sphagnum peatlands spread. This suggests that low-temperature fires are a typical characteristic of Siberian larch forests also on long timescales. The results also suggest that low-intensity fires would be reduced by vegetation shifts towards very dry environments due to reduced biomass availability, as well as by shifts towards peatlands, which limits fuel dryness. In addition, we observed very low MA ratios, which we interpret as high contributions of galactosan and mannosan from biomass sources other than those currently monitored, such as the moss-lichen mats in the understorey of the summergreen boreal forest. Overall, sedimentary MAs can provide a powerful proxy for fire regime reconstructions and extend our knowledge of long-term natural fire-climate-vegetation feedbacks in the high northern latitudes.
Landscapes in high northern latitudes are assumed to be highly sensitive to future global change, but the rates and long-term trajectories of changes are rather uncertain. In the boreal zone, fires are an important factor in climate-vegetation interactions and biogeochemical cycles. Fire regimes are characterized by small, frequent, low-intensity fires within summergreen boreal forests dominated by larch, whereas evergreen boreal forests dominated by spruce and pine burn large areas less frequently but at higher intensities. Here, we explore the potential of the monosaccharide anhydrides (MA) levoglucosan, mannosan and galactosan to serve as proxies of low-intensity biomass burning in glacial-to-interglacial lake sediments from the high northern latitudes. We use sediments from Lake El'gygytgyn (cores PG 1351 and ICDP 5011-1), located in the far north-east of Russia, and study glacial and interglacial samples of the last 430 kyr (marine isotope stages 5e, 6, 7e, 8, 11c and 12) that had different climate and biome configurations. Combined with pollen and non-pollen palynomorph records from the same samples, we assess how far the modern relationships between fire, climate and vegetation persisted during the past, on orbital to centennial timescales. We find that MAs attached to particulates were well-preserved in up to 430 kyr old sediments with higher influxes from low-intensity biomass burning in interglacials compared to glacials. MA influxes significantly increase when summergreen boreal forest spreads closer to the lake, whereas they decrease when tundra-steppe environments and, especially, Sphagnum peatlands spread. This suggests that low-temperature fires are a typical characteristic of Siberian larch forests also on long timescales. The results also suggest that low-intensity fires would be reduced by vegetation shifts towards very dry environments due to reduced biomass availability, as well as by shifts towards peatlands, which limits fuel dryness. In addition, we observed very low MA ratios, which we interpret as high contributions of galactosan and mannosan from biomass sources other than those currently monitored, such as the moss-lichen mats in the understorey of the summergreen boreal forest. Overall, sedimentary MAs can provide a powerful proxy for fire regime reconstructions and extend our knowledge of long-term natural fire-climate-vegetation feedbacks in the high northern latitudes.
Recent policy changes highlight the need for citizens to take adaptive actions to reduce flood-related impacts. Here, we argue that these changes represent a wider behavioral turn in flood risk management (FRM). The behavioral turn is based on three fundamental assumptions: first, that the motivations of citizens to take adaptive actions can be well understood so that these motivations can be targeted in the practice of FRM; second, that private adaptive measures and actions are effective in reducing flood risk; and third, that individuals have the capacities to implement such measures. We assess the extent to which the assumptions can be supported by empirical evidence. We do this by engaging with three intellectual catchments. We turn to research by psychologists and other behavioral scientists which focus on the sociopsychological factors which influence individual motivations (Assumption 1). We engage with economists, engineers, and quantitative risk analysts who explore the extent to which individuals can reduce flood related impacts by quantifying the effectiveness and efficiency of household-level adaptive measures (Assumption 2). We converse with human geographers and sociologists who explore the types of capacities households require to adapt to and cope with threatening events (Assumption 3). We believe that an investigation of the behavioral turn is important because if the outlined assumptions do not hold, there is a risk of creating and strengthening inequalities in FRM. Therefore, we outline the current intellectual and empirical knowledge as well as future research needs. Generally, we argue that more collaboration across intellectual catchments is needed, that future research should be more theoretically grounded and become methodologically more rigorous and at the same time focus more explicitly on the normative underpinnings of the behavioral turn.
Recent policy changes highlight the need for citizens to take adaptive actions to reduce flood-related impacts. Here, we argue that these changes represent a wider behavioral turn in flood risk management (FRM). The behavioral turn is based on three fundamental assumptions: first, that the motivations of citizens to take adaptive actions can be well understood so that these motivations can be targeted in the practice of FRM; second, that private adaptive measures and actions are effective in reducing flood risk; and third, that individuals have the capacities to implement such measures. We assess the extent to which the assumptions can be supported by empirical evidence. We do this by engaging with three intellectual catchments. We turn to research by psychologists and other behavioral scientists which focus on the sociopsychological factors which influence individual motivations (Assumption 1). We engage with economists, engineers, and quantitative risk analysts who explore the extent to which individuals can reduce flood related impacts by quantifying the effectiveness and efficiency of household-level adaptive measures (Assumption 2). We converse with human geographers and sociologists who explore the types of capacities households require to adapt to and cope with threatening events (Assumption 3). We believe that an investigation of the behavioral turn is important because if the outlined assumptions do not hold, there is a risk of creating and strengthening inequalities in FRM. Therefore, we outline the current intellectual and empirical knowledge as well as future research needs. Generally, we argue that more collaboration across intellectual catchments is needed, that future research should be more theoretically grounded and become methodologically more rigorous and at the same time focus more explicitly on the normative underpinnings of the behavioral turn.
Model-derived relationships between chlorophyll a (Chl-a) and nutrients and temperature have fundamental implications for understanding complex interactions among water quality measures used for lake classification, yet accuracy comparisons of different approaches are scarce. Here, we (1) compared Chl-a model performances across linear and nonlinear statistical approaches; (2) evaluated single and combined effects of nutrients, depth, and temperature as lake surface water temperature (LSWT) or altitude on Chl-a; and (3) investigated the reliability of the best water quality model across 13 lakes from perialpine and central Balkan mountain regions. Chl-a was modelled using in situ water quality data from 157 European lakes; elevation data and LSWT in situ data were complemented by remote sensing measurements. Nonlinear approaches performed better, implying complex relationships between Chl-a and the explanatory variables. Boosted regression trees, as the best performing approach, accommodated interactions among predictor variables. Chl-a-nutrient relationships were characterized by sigmoidal curves, with total phosphorus having the largest explanatory power for our study region. In comparison with LSWT, utilization of altitude, the often-used temperature surrogate, led to different influence directions but similar predictive performances. These results support utilizing altitude in models for Chl-a predictions. Compared to Chl-a observations, Chl-a predictions of the best performing approach for mountain lakes (oligotrophic-eutrophic) led to minor differences in trophic state categorizations. Our findings suggest that both models with LSWT and altitude are appropriate for water quality predictions of lakes in mountain regions and emphasize the importance of incorporating interactions among variables when facing lake management challenges.
Model-derived relationships between chlorophyll a (Chl-a) and nutrients and temperature have fundamental implications for understanding complex interactions among water quality measures used for lake classification, yet accuracy comparisons of different approaches are scarce. Here, we (1) compared Chl-a model performances across linear and nonlinear statistical approaches; (2) evaluated single and combined effects of nutrients, depth, and temperature as lake surface water temperature (LSWT) or altitude on Chl-a; and (3) investigated the reliability of the best water quality model across 13 lakes from perialpine and central Balkan mountain regions. Chl-a was modelled using in situ water quality data from 157 European lakes; elevation data and LSWT in situ data were complemented by remote sensing measurements. Nonlinear approaches performed better, implying complex relationships between Chl-a and the explanatory variables. Boosted regression trees, as the best performing approach, accommodated interactions among predictor variables. Chl-a-nutrient relationships were characterized by sigmoidal curves, with total phosphorus having the largest explanatory power for our study region. In comparison with LSWT, utilization of altitude, the often-used temperature surrogate, led to different influence directions but similar predictive performances. These results support utilizing altitude in models for Chl-a predictions. Compared to Chl-a observations, Chl-a predictions of the best performing approach for mountain lakes (oligotrophic-eutrophic) led to minor differences in trophic state categorizations. Our findings suggest that both models with LSWT and altitude are appropriate for water quality predictions of lakes in mountain regions and emphasize the importance of incorporating interactions among variables when facing lake management challenges.
Much of contemporary landslide research is concerned with predicting and mapping susceptibility to slope failure. Many studies rely on generalised linear models with environmental predictors that are trained with data collected from within and outside of the margins of mapped landslides. Whether and how the performance of these models depends on sample size, location, or time remains largely untested. We address this question by exploring the sensitivity of a multivariate logistic regression-one of the most widely used susceptibility models-to data sampled from different portions of landslides in two independent inventories (i.e. a historic and a multi-temporal) covering parts of the eastern rim of the Fergana Basin, Kyrgyzstan. We find that considering only areas on lower parts of landslides, and hence most likely their deposits, can improve the model performance by >10% over the reference case that uses the entire landslide areas, especially for landslides of intermediate size. Hence, using landslide toe areas may suffice for this particular model and come in useful where landslide scars are vague or hidden in this part of Central Asia. The model performance marginally varied after progressively updating and adding more landslides data through time. We conclude that landslide susceptibility estimates for the study area remain largely insensitive to changes in data over about a decade. Spatial or temporal stratified sampling contributes only minor variations to model performance. Our findings call for more extensive testing of the concept of dynamic susceptibility and its interpretation in data-driven models, especially within the broader framework of landslide risk assessment under environmental and land-use change.
Vegetation has long been hypothesized to influence the nature and rates of surface processes. We test the possible impact of vegetation and climate on denudation rates at orogen scale by taking advantage of a pronounced along-strike gradient in rainfall and vegetation density in the Himalaya. We combine 12 new 10Be denudation rates from the Sutlej Valley and 123 published denudation rates from fluvially- dominated catchments in the Himalaya with remotely-sensed measures of vegetation density and rainfall metrics, and with tectonic and lithologic constraints. In addition, we perform topographic analyses to assess the contribution of vegetation and climate in modulating denudation rates along strike. We observe variations in denudation rates and the relationship between denudation and topography along strike that are most strongly controlled by local rainfall amount and vegetation density, and cannot be explained by along-strike differences in tectonics or lithology. A W–E along-strike decrease in denudation rate variability positively correlates with the seasonality of vegetation density (R = 0.95, p < 0.05), and negatively correlates with mean vegetation density (R = −0.84, p < 0.05). Vegetation density modulates the topographic response to changing denudation rates, such that the functional relationship between denudation rate and topographic steepness becomes increasingly linear as vegetation density increases. We suggest that while tectonic processes locally control the pattern of denudation rates across strike of the Himalaya (i.e., S–N), along strike of the orogen (i.e., E–W) climate exerts a measurable influence on how denudation rates scatter around long-term, tectonically-controlled erosion, and on the functional relationship between topography and denudation
We propose a novel way to measure and analyze networks of drainage divides from digital elevation models. We developed an algorithm that extracts drainage divides based on the drainage basin boundaries defined by a stream network. In contrast to streams, there is no straightforward approach to order and classify divides, although it is intuitive that some divides are more important than others. A meaningful way of ordering divides is the average distance one would have to travel down on either side of a divide to reach a common stream location. However, because measuring these distances is computationally expensive and prone to edge effects, we instead sort divide segments based on their tree-like network structure, starting from endpoints at river confluences. The sorted nature of the network allows for assigning distances to points along the divides, which can be shown to scale with the average distance downslope to the common stream location. Furthermore, because divide segments tend to have characteristic lengths, an ordering scheme in which divide orders increase by 1 at junctions mimics these distances. We applied our new algorithm to the Big Tujunga catchment in the San Gabriel Mountains of southern California and studied the morphology of the drainage divide network. Our results show that topographic metrics, like the downstream flow distance to a stream and hillslope relief, attain characteristic values that depend on the drainage area threshold used to derive the stream network. Portions along the divide network that have lower than average relief or are closer than average to streams are often distinctly asymmetric in shape, suggesting that these divides are unstable. Our new and automated approach thus helps to objectively extract and analyze divide networks from digital elevation models.
A mid-aperture broad-band test array (OBS array DOCTAR) was deployed from June 2011 to April 2012 about 100 km north of the Gloria fault in the Eastern North Atlantic in about 5000 m water depth. In addition arrays were installed on Madeira Island and in western Portugal mainland. For the first time in the Eastern North Atlantic, we recorded a large number of high frequency Po and So waves from local and regional small and moderate earthquakes (M-L < 4). An incoherent beamforming method was adapted to scan continuous data for such Po and So arrivals applying a sliding window waveform migration and frequency-wavenumber technique. We identify about 320 Po and 1550 So arrivals and compare the phase onsets with the ISC catalogue (ISC 2015) for the same time span. Up to a distance of 6 degrees to the DOCTAR stations all events listed in the ISC catalogue could be associated to Po and So phases. Arrivals from events in more than 10 degrees distance could be identified only in some cases. Only few Po and/or So arrivals were detected for earthquakes from the European and African continental area, the continental shelf regions and for earthquakes within or northwest of the Azores plateau. Unexpectedly, earthquake clusters are detected within the oceanic plates north and south of the Gloria fault and far from plate boundaries, indicating active intraplate structures. We also observe and locate numerous small magnitude earthquakes on the segment of the Gloria fault directly south of DOCTAR, which likely coincides with the rupture of the 25 November 1941 event. Local small magnitude earthquakes located beneath DOCTAR show hypocentres up to 30 km depth and strike-slip focal mechanisms. A comparison with detections at temporary mid-aperture arrays on Madeira and in western Portugal shows that the deep ocean array performs much better than the island and the continental array regarding the detection threshold for events in the oceanic plates.
We conclude that sparsely distributed mid-aperture seismic arrays in the deep ocean could decrease the detection and location threshold for seismicity with M-L < 4 in the oceanic plate and might constitute a valuable tool to monitor oceanic plate seismicity.
Detailed organic geochemical and carbon isotopic (delta C-13 and Delta C-14) analyses are performed on permafrost deposits affected by coastal erosion (Herschel Island, Canadian Beaufort Sea) and adjacent marine sediments (Herschel Basin) to understand the fate of organic carbon in Arctic nearshore environments. We use an end-member model based on the carbon isotopic composition of bulk organic matter to identify sources of organic carbon. Monte Carlo simulations are applied to quantify the contribution of coastal permafrost erosion to the sedimentary carbon budget. The models suggest that similar to 40% of all carbon released by local coastal permafrost erosion is efficiently trapped and sequestered in the nearshore zone. This highlights the importance of sedimentary traps in environments such as basins, lagoons, troughs, and canyons for the carbon sequestration in previously poorly investigated, nearshore areas.
Plain Language Summary Increasing air and sea surface temperatures at high latitudes leads to accelerated thaw, destabilization, and erosion of perennially frozen soils (i.e., permafrost), which are often rich in organic carbon. Coastal erosion leads to an increased mobilization of organic carbon into the Arctic Ocean, which there can be converted into greenhouse gases and may therefore contribute to further warming. Carbon decomposition can be limited if organic matter is efficiently deposited on the seafloor, buried in marine sediments, and thus removed from the short-term carbon cycle. Basins, canyons, and troughs near the coastline can serve as sediment traps and potentially accommodate large quantities of organic carbon along the Arctic coast. Here we use biomarkers (source-specific molecules), stable carbon isotopes, and radiocarbon to identify the sources of organic carbon in the nearshore zone of the southern Canadian Beaufort Sea near Herschel Island. We quantify the contribution of coastal permafrost erosion to the sedimentary carbon budget of the area and estimate that more than a third of all carbon released by local permafrost erosion is efficiently trapped in marine sediments. This highlights the importance of regional sediment traps for carbon sequestration.
We study the source properties of the 2005 Kashmir earthquake and its aftershocks to unravel the seismotectonics of the NW Himalayan syntaxis. The mainshock and larger aftershocks have been simultaneously relocated using phase data. We use back-projection of high-frequency energy from multiple teleseismic arrays to model the spatio-temporal evolution of the mainshock rupture. Our analysis reveal a bilateral rupture, which initially propagated SE and then NW of the epicenter, with an average rupture velocity of similar to 2 km s(-1). The area of maximum energy release is parallel to and bound by the surface rupture. Incorporating rupture propagation and velocity, we model the mainshock as a line source using P- and SH-waveform inversion. Our result confirms that the mainshock occurred on a NE dipping (similar to 35 degrees) fault plane, with centroid depth of similar to 10 km. Integrated source time function show that majority of the energy was released in the first similar to 20 s, and was confined above the hypocenter. From waveform inverted fault dimension and seismic moment, we argue that the mainshock had an additional similar to 25 km blind rupture beyond the NW Himalayan syntaxis. Combining this with findings from previous studies, we conjecture that the blind rupture propagated NW of the syntaxis underneath a weak detachment overlain by infra-Cambrian salt layer, and terminated in a wedge thrust. All moderate-to-large aftershocks, NW of the mainshock rupture, are concentrated at the edge of the blind rupture termination. Source modeling of these aftershocks reveal thrust mechanism with centroid depths of 2-10 km, and fault planes oriented subparallel to the mainshock rupture. To study the influence of mainshock rupture on aftershock occurrence, we compute Coulomb failure stress on aftershock faults. All these aftershocks lie in the positive Coulomb stress change region. This suggest that the aftershocks have been triggered by either co-seismic or post-seismic slip on the mainshock fault.
Multiple P-T-d-t paths reveal the evolution of the final Nuna assembly in northeast Australia
(2020)
The final assembly of the Mesoproterozoic supercontinent Nuna was marked by the collision of Laurentia and Australia at 1.60 Ga, which is recorded in the Georgetown Inlier of NE Australia. Here, we decipher the metamorphic evolution of this final Nuna collisional event using petrostructural analysis, major and trace element compositions of key minerals, thermodynamic modelling, and multi-method geochronology. The Georgetown Inlier is characterised by deformed and metamorphosed 1.70-1.62 Ga sedimentary and mafic rocks, which were intruded byc. 1.56 Ga old S-type granites. Garnet Lu-Hf and monazite U-Pb isotopic analyses distinguish two major metamorphic events (M1 atc. 1.60 Ga and M2 atc. 1.55 Ga), which allows at least two composite fabrics to be identified at the regional scale-c. 1.60 Ga S1 (consisting in fabrics S1a and S1b) andc. 1.55 Ga S2 (including fabrics S2a and S2b). Also, three tectono-metamorphic domains are distinguished: (a) the western domain, with S1 defined by low-P(LP) greenschist facies assemblages; (b) the central domain, where S1 fabric is preserved as medium-P(MP) amphibolite facies relicts, and locally as inclusion trails in garnet wrapped by the regionally dominant low-Pamphibolite facies S2 fabric; and (c) the eastern domain dominated by upper amphibolite to granulite facies S2 foliation. In the central domain, 1.60 GaMP-medium-T(MT) metamorphism (M1) developed within the staurolite-garnet stability field, with conditions ranging from 530-550 degrees C at 6-7 kbar (garnet cores) to 620-650 degrees C at 8-9 kbar (garnet rims), and it is associated with S1 fabric. The onset of 1.55 GaLP-high-T(HT) metamorphism (M2) is marked by replacement of staurolite by andalusite (M2a/D2a), which was subsequently pseudomorphed by sillimanite (M2b/D2b) where granite and migmatite are abundant.P-Tconditions ranged from 600 to 680 degrees C and 4-6 kbar for the M2b sillimanite stage. 1.60 Ga garnet relicts within the S2 foliation highlight the progressive obliteration of the S1 fabric by regional S2 in the central zone during peak M2 metamorphism. In the eastern migmatitic complex, partial melting of paragneiss and amphibolite occurred syn- to post-S2, at 730-770 degrees C and 6-8 kbar, and at 750-790 degrees C and 6 kbar, respectively. The pressure-temperature-deformation-time paths reconstructed for the Georgetown Inlier suggest ac. 1.60 Ga M1/D1 event recorded under greenschist facies conditions in the western domain and under medium-Pand medium-Tconditions in the central domain. This event was followed by the regional 1.56-1.54 Ga low-Pand high-Tphase (M2/D2), extensively recorded in the central and eastern domains. Decompression between these two metamorphic events is ascribed to an episode of exhumation. The two-stage evolution supports the previous hypothesis that the Georgetown Inlier preserves continental collisional and subsequent thermal perturbation associated with granite emplacement.
Diffusive transport and sorption processes of uranium in the Swiss Opalinus Clay were investigated as a function of partial pressure of carbon dioxide pCO(2), varying mineralogy in the facies and associated changes in porewater composition. Simulations were conducted in one-dimensional diffusion models on the 100 m-scale for a time of one million years using a bottom-up approach based on mechanistic surface complexation models as well as cation exchange to quantify sorption. Speciation calculations have shown, uranium is mainly present as U(VI) and must therefore be considered as mobile for in-situ conditions. Uranium migrated up to 26 m in both, the sandy and the carbonate-rich facies, whereas in the shaly facies 16 m was the maximum. The main species was the anionic complex CaUO2(CO3)(3)(2-) . Hence, anion exclusion was taken into account and further reduced the migration distances by 30 %. The concentrations of calcium and carbonates reflected by the set pCO(2) determine speciation and activity of uranium and consequently the sorption behaviour. Our simulation results allow for the first time to prioritize on the far-field scale the governing parameters for diffusion and sorption of uranium and hence outline the sensitivity of the system. Sorption processes are controlled in descending priority by the carbonate and calcium concentrations, pH, pe and the clay mineral content. Therefore, the variation in porewater composition resulting from the heterogeneity of the facies in the Opalinus Clay formation needs to be considered in the assessment of uranium migration in the far field of a potential repository.
Barite scales in geothermal installations are a highly unwanted effect of circulating deep saline fluids. They build up in the reservoir if supersaturated fluids are re-injected, leading to irreversible loss of injectivity. A model is presented for calculating the total expected barite precipitation. To determine the related injectivity decline over time, the spatial precipitation distribution in the subsurface near the injection well is assessed by modelling barite growth kinetics in a radially diverging Darcy flow domain. Flow and reservoir properties as well as fluid chemistry are chosen to represent reservoirs subject to geothermal exploration located in the North German Basin (NGB) and the Upper Rhine Graben (URG) in Germany. Fluids encountered at similar depths are hotter in the URG, while they are more saline in the NGB. The associated scaling amount normalised to flow rate is similar for both regions. The predicted injectivity decline after 10 years, on the other hand, is far greater for the NGB (64%) compared to the URG (24%), due to the temperature- and salinity-dependent precipitation rate. The systems in the NGB are at higher risk. Finally, a lightweight score is developed for approximating the injectivity loss using the Damkohler number, flow rate and total barite scaling potential. This formula can be easily applied to geothermal installations without running complex reactive transport simulations.
We present a novel data set of subnational economic output, Gross Regional Product (GRP), for more than 1500 regions in 77 countries that allows us to empirically estimate historic climate impacts at different time scales. Employing annual panel models, long-difference regressions and cross-sectional regressions, we identify effects on productivity levels and productivity growth. We do not find evidence for permanent growth rate impacts but we find robust evidence that temperature affects productivity levels considerably. An increase in global mean surface temperature by about 3.5°C until the end of the century would reduce global output by 7–14% in 2100, with even higher damages in tropical and poor regions. Updating the DICE damage function with our estimates suggests that the social cost of carbon from temperature-induced productivity losses is on the order of 73–142$/tCO2 in 2020, rising to 92–181$/tCO2 in 2030. These numbers exclude non-market damages and damages from extreme weather events or sea-level rise.
Atmospheric dynamics of extreme discharge events from 1979 to 2016 in the southern Central Andes
(2020)
During the South-American Monsoon season, deep convective systems occur at the eastern flank of the Central Andes leading to heavy rainfall and flooding. We investigate the large- and meso-scale atmospheric dynamics associated with extreme discharge events (> 99.9th percentile) observed in two major river catchments meridionally stretching from humid to semi-arid conditions in the southern Central Andes. Based on daily gauge time series and ERA-Interim reanalysis, we made the following three key observations: (1) for the period 1940-2016 daily discharge exhibits more pronounced variability in the southern, semi-arid than in the northern, humid catchments. This is due to a smaller ratio of discharge magnitudes between intermediate (0.2 year return period) and rare events (20 year return period) in the semi-arid compared to the humid areas; (2) The climatological composites of the 40 largest discharge events showed characteristic atmospheric features of cold surges based on 5-day time-lagged sequences of geopotential height at different levels in the troposphere; (3) A subjective classification revealed that 80% of the 40 largest discharge events are mainly associated with the north-northeastward migration of frontal systems and 2/3 of these are cold fronts, i.e. cold surges. This work highlights the importance of cold surges and their related atmospheric processes for the generation of heavy rainfall events and floods in the southern Central Andes.
This study, based on new and high quality in situ observations, quantifies for the first time, the individual contributions of light-absorbing aerosols (black carbon (BC), brown carbon (BrC) and dust) to aerosol absorption over the Indo-Gangetic Plain (IGP) and the Himalayan foothill region, a relatively poorly studied region with several sensitive ecosystems of global importance, as well as highly vulnerable populations. The annual and seasonal average single scattering albedo (SSA) over Kathmandu is the lowest of all the locations. The SSA over Kathmandu is < 0.89 during all seasons, which confirms the dominance of light-absorbing carbonaceous aerosols from local and regional sources over Kathmandu. It is observed here that the SSA decreases with increasing elevation, confirming the dominance of light absorbing carbonaceous aerosols at higher elevations. In contrast, the SSA over the IGP does not exhibit a pronounced spatial variation. BC dominates (>= 75%) the aerosol absorption over the IGP and the Himalayan foothills throughout the year. Higher BC concentration at elevated locations in the Himalayas leads to lower SSA at elevated locations in the Himalayas. The contribution of dust to aerosol absorption is higher throughout the year over the IGP than over the Himalayan foothills. The aerosol absorption over South Asia is very high, exceeding available observations over East Asia, and also exceeds previous model estimates. This quantification will be valuable as observational constraints to help improve regional simulations of climate change, impacts on the glaciers and the hydrological cycle, and will help to direct the focus towards BC as the main contributor to aerosol-induced warming in the region.
We analyze the spatiotemporal evolution of seismicity during a sequence of moderate (an M-w 4.7 foreshock and M-w 5.8 mainshock) earthquakes occurring in September 2019 at the transition between a creeping and a locked segment of the North Anatolian fault in the central Sea of Marmara, northwest Turkey. To investigate in detail the seismicity evolution, we apply a matched-filter technique to continuous waveforms, thus reducing the magnitude threshold for detection. Sequences of foreshocks preceding the two largest events are clearly seen, exhibiting two different behaviors: a long-term activation of the seismicity along the entire fault segment and a short-term concentration around the epicenters of the large events. We suggest a two-scale preparation phase, with aseismic slip preparing the mainshock final rupture a few days before, and a cascade mechanism leading to the nucleation of the mainshock. Thus, our study shows a combination of seismic and aseismic slip during the foreshock sequence changing the strength of the fault, bringing it closer to failure.
During the period 750-600 Ma ago, prior to the final break-up of the supercontinent Rodinia, the crust of both the North American Craton and Baltica was intruded by significant amounts of rift-related magmas originating from the mantle. In the Proterozoic crust of Southern Norway, the 580 Ma old Fen carbonatite-ultramafic complex is a representative of this type of rocks. In this paper, we report the occurrence of an ultramafic lamprophyre dyke which possibly is linked to the Fen complex, although Ar-40/Ar-39 data from phenocrystic phlogopite from the dyke gave an age of 686 +/- 9 Ma. The lamprophyre dyke was recently discovered in one of the Kongsberg silver mines at Vinoren, Norway. Whole rock geochemistry, geochronological and mineralogical data from the ultramafic lamprophyre dyke are presented aiming to elucidate its origin and possible geodynamic setting. From the whole-rock composition of the Vinoren dyke, the rock could be recognized as transitional between carbonatite and kimberlite-II (orangeite). From its diagnostic mineralogy, the rock is classified as aillikite. The compositions and xenocrystic nature of several of the major and accessory minerals from the Vinoren aillikite are characteristic for diamondiferous rocks (kimberlites/lamproites/UML): Phlogopite with kinoshitalite-rich rims, chromite-spinel-ulvospinel series, Mg- and Mn-rich ilmenites, rutile and lucasite-(Ce). We suggest that the aillikite melt formed during partial melting of a MARID (mica-amphibole-rutile-ilmenite-diopside)-like source under CO2 fluxing. The pre-rifting geodynamic setting of the Vinoren aillikite before the Rodinia supercontinent breakup suggests a relatively thick SCLM (Subcontinental Lithospheric Mantle) during this stage and might indicate a diamond-bearing source for the parental melt. This is in contrast to the about 100 Ma younger Fen complex, which were derived from a thin SCLM.
Geopolitical shifts and the changing significance of borders in the EU's neighbourhood are usually understood as a matter of international power politics. Factors that accompany geopolitical impact on borders, such as media coverage of geopolitical change, often appear as secondary or irrelevant. However the recent Ukraine conflict revealed the contrary as pro-EU attitudes were strongly supported by 'western' media. Therefore this paper seeks to clarify the role of news media in creating perspectives and attitudes on geopolitical shifts and the significance of European borders. Empirical evidence on the coverage of the evolving Ukraine crisis by German news sources portrays the media as promoters of biased framings and imaginaries which suggest that the EU be a potential conflict party in the newly evolving geostrategic confrontation in its eastern neighbourhood. The findings indicate that during critical periods of the Ukraine crisis media reports combined rising euphoria about Europe and 'the West', as defenders of the 'good cause', with excessive moral polarising and the discursive normalisation of a rhetoric of escalation. Imaginaries of a bipolar world (The West against Russia) and a new Cold War prepared the ground for a new understanding of European borders and neighbourhood relations as being manipulable at will.
The steady increase of ground-motion data not only allows new possibilities but also comes with new challenges in the development of ground-motion models (GMMs). Data classification techniques (e.g., cluster analysis) do not only produce deterministic classifications but also probabilistic classifications (e.g., probabilities for each datum to belong to a given class or cluster). One challenge is the integration of such continuous classification in regressions for GMM development such as the widely used mixed-effects model. We address this issue by introducing an extension of the mixed-effects model to incorporate data weighting. The parameter estimation of the mixed-effects model, that is, fixed-effects coefficients of the GMMs and the random-effects variances, are based on the weighted likelihood function, which also provides analytic uncertainty estimates. The data weighting permits for earthquake classification beyond the classical, expert-driven, binary classification based, for example, on event depth, distance to trench, style of faulting, and fault dip angle. We apply Angular Classification with Expectation-maximization, an algorithm to identify clusters of nodal planes from focal mechanisms to differentiate between, for example, interface- and intraslab-type events. Classification is continuous, that is, no event belongs completely to one class, which is taken into account in the ground-motion modeling. The theoretical framework described in this article allows for a fully automatic calibration of ground-motion models using large databases with automated classification and processing of earthquake and ground-motion data. As an example, we developed a GMM on the basis of the GMM by Montalva et al. (2017) with data from the strong-motion flat file of Bastias and Montalva (2016) with similar to 2400 records from 319 events in the Chilean subduction zone. Our GMM with the data-driven classification is comparable to the expert-classification-based model. Furthermore, the model shows temporal variations of the between-event residuals before and after large earthquakes in the region.
High-pressure experiments were performed to investigate the effectiveness, rate and mechanism of carbonation of serpentinites by a carbon-saturated COH fluid at 1.5-2.5 GPa and 375-700 degrees C. This allows a better understanding of the fate and redistribution of slab-derived carbonic fluids when they react with the partially hydrated mantle within and above the subducting slab under pressure and temperature conditions corresponding to the forearc mantle. Interactions between carbon-saturated CO2-H2O-CH4 fluids and serpentinite were investigated using natural serpentinite cylinders with natural grain sizes and shapes in piston-cylinder experiments. The volatile composition of post-run fluids was quantified by gas chromatography. Solid phases were examined by Raman spectroscopy, electron microscopy and laser ablation inductively coupled plasma mass spectrometry. Textures, porosity and phase abundances of recovered rock cores were visualized and quantified by three-dimensional, high-resolution computed tomography. We find that carbonation of serpentinites is efficient at sequestering CO2 from the interacting fluid into newly formed magnesite. Time-series experiments demonstrate that carbonation is completed within similar to 96 h at 2 GPa and 600 degrees C. With decreasing CO2, aq antigorite is replaced first by magnesite + quartz followed by magnesite + talc + chlorite in distinct, metasomatic fronts. Above antigorite stability magnesite + enstatite + talc + chlorite occur additionally. The formation of fluid-permeable reaction zones enhances the reaction rate and efficiency of carbonation. Carbonation probably occurs via an interface-coupled replacement process, whereby interconnected porosity is present within reaction zones after the experiment. Consequently, carbonation of serpentinites is self-promoting and efficient even if fluid flow is channelized into veins. We conclude that significant amounts of carbonates may accumulate, over time, in the hydrated forearc mantle.
Aerosol emissions from human activities are extensive and changing rapidly over Asia. Model simulations and satellite observations indicate a dipole pattern in aerosol emissions and loading between South Asia and East Asia, two of the most heavily polluted regions of the world. We examine the previously unexplored diverging trends in the existing dipole pattern of aerosols between East and South Asia using the high quality, two-decade long ground-based time series of observations of aerosol properties from the Aerosol Robotic Network (AERONET), from satellites (Moderate Resolution Imaging Spectroradiometer (MODIS) and Ozone Monitoring Instrument (OMI)), and from model simulations (Modern-Era Retrospective Analysis for Research and Applications, version 2 (MERRA-2). The data cover the period since 2001 for Kanpur (South Asia) and Beijing (East Asia), two locations taken as being broadly representative of the respective regions. Since 2010 a dipole in aerosol optical depth (AOD) is maintained, but the trend is reversed-the decrease in AOD over Beijing (East Asia) is rapid since 2010, being 17% less in current decade compared to first decade of twenty-first century, while the AOD over South Asia increased by 12% during the same period. Furthermore, we find that the aerosol composition is also changing over time. The single scattering albedo (SSA), a measure of aerosol's absorption capacity and related to aerosol composition, is slightly higher over Beijing than Kanpur, and has increased from 0.91 in 2002 to 0.93 in 2017 over Beijing and from 0.89 to 0.92 during the same period over Kanpur, confirming that aerosols in this region have on an average become more scattering in nature. These changes have led to a notable decrease in aerosol-induced atmospheric heating rate (HR) over both regions between the two decades, decreasing considerably more over East Asia (- 31%) than over South Asia (- 9%). The annual mean HR is lower now, it is still large (>= 0.6 K per day), which has significant climate implications. The seasonal trends in AOD, SSA and HR are more pronounced than their respective annual trends over both regions. The seasonal trends are caused mainly by the increase/decrease in anthropogenic aerosol emissions (sulfate, black carbon and organic carbon) while the natural aerosols (dust and sea salt) did not change significantly over South and East Asia during the last two decades. The MERRA-2 model is able to simulate the observed trends in AODs well but not the magnitude, while it also did not simulate the SSA values or trends well. These robust findings based on observations of key aerosol parameters and previously unrecognized diverging trends over South and East Asia need to be accounted for in current state-of-the-art climate models to ensure accurate quantification of the complex and evolving impact of aerosols on the regional climate over Asia.
Forested areas are assumed not to be influenced by erosion processes. However, forest soils of Northern Germany in a hummocky ground moraine landscape can sometimes exhibit a very shallow thickness on crest positions and buried soils on slope positions. The question consequently is: Are these on-going or ancient erosional and depositional processes? Plutonium isotopes act as soil erosion/deposition tracers for recent (last few decades) processes. Here, we quantified the 239+240PU inventories in a small, forested catchment (ancient forest "Melzower Forst", deciduous trees), which is characterised by a hummocky terrain including a kettle hole. Soil development depths (depth to C horizon) and 239+240PU inventories along a catena of sixteen different profiles were determined and correlated to relief parameters. Moreover, we compared different modelling approaches to derive erosion rates from Pu data. <br /> We find a strong relationship between soil development depths, distance-to-sink and topography along the catena. Fully developed Retisols (thicknesses > 1 m) in the colluvium overlay old land surfaces as documented by fossil Ah horizons. However, we found no relationship of Pu-based erosion rates to any relief parameter. Instead, 239+240PU inventories showed a very high local, spatial variability (36-70 Bq m(-2)). Low annual rainfall, spatially distributed interception and stem flow might explain the high variability of the 239+240PU inventories, giving rise to a patchy input pattern. Different models resulted in quite similar erosion and deposition rates (max: -5 t ha(-1) yr(-1) to +7.3 t ha(-1) yr(-1)). Although some rates are rather high, the magnitude of soil erosion and deposition - in terms of soil thickness change - is negligible during the last 55 years. The partially high values are an effect of the patchy Pu deposition on the forest floor. This forest has been protected for at least 240 years. Therefore rather natural events and anthropogenic activities during medieval times or even earlier must have caused the observed soil pattern, which documents strong erosion and deposition processes.
Corona and the climate
(2020)
The Frasnian-Famennian (F-F) transition of Late Devonian was a critical episode in geological history, recording a major mass extinction event. In this study, we focus on an F-F succession from a deep marine context in Bancheng, southern Guangxi, South China, to investigate coeval changes in pelagic environments of the Paleo-Tethys Ocean. The studied succession is exclusively composed of bedded cherts intercalated with multiple siliceous volcanic ash beds. A SIMS zircon U-Pb Concordia age of 367.8 +/- 2.5 Ma is reported for a tuffaceous layer slightly above the F-F boundary. Geochemical ratios of Al/(Al + Fe + Mn), Ce/Ce*, Y/Ho, and Al, Fe contents in bedded cherts indicate that they are of predominantly biogenic/chemical origin with some terrigenous inputs. Negligible enrichment of redox sensitive elements (Mo, U, V) and low V/Cr ratios (<2) suggest persistently oxic conditions existed in the deep pelagic basin at Bancheng, South China during the F-F transition. These findings call into question the widely held hypothesis that marine anoxia was the primary killing mechanism for the F-F crisis. In contrast, multiple tuffaceous layers throughout the F-F boundary succession indicate frequent volcanic activity, which could have released massive amounts of greenhouse gases into the atmosphere, inducing climate warming. This scenario may have increased continental weathering and riverine fluxes into the ocean, reconciling the increases in Al2O3 content and Al/(Al + Fe + Mn) ratio across the F-F boundary. Documentation of persistently oxic conditions and frequent volcanic activitiy provides new perspectives on the inter-relationship between volcanism, climate, and oceanic redox fluctuation during the F-F biotic crisis.
Compound natural hazards likeEl Ninoevents cause high damage to society, which to manage requires reliable risk assessments. Damage modelling is a prerequisite for quantitative risk estimations, yet many procedures still rely on expert knowledge, and empirical studies investigating damage from compound natural hazards hardly exist. A nationwide building survey in Peru after theEl Ninoevent 2017 - which caused intense rainfall, ponding water, flash floods and landslides - enables us to apply data-mining methods for statistical groundwork, using explanatory features generated from remote sensing products and open data. We separate regions of different dominant characteristics through unsupervised clustering, and investigate feature importance rankings for classifying damage via supervised machine learning. Besides the expected effect of precipitation, the classification algorithms select the topographic wetness index as most important feature, especially in low elevation areas. The slope length and steepness factor ranks high for mountains and canyons. Partial dependence plots further hint at amplified vulnerability in rural areas. An example of an empirical damage probability map, developed with a random forest model, is provided to demonstrate the technical feasibility.
Different upper tail indicators exist to characterize heavy tail phenomena, but no comparative study has been carried out so far. We evaluate the shape parameter (GEV), obesity index, Gini index and upper tail ratio (UTR) against a novel benchmark of tail heaviness - the surprise factor. Sensitivity analyses to sample size and changes in scale-to-location ratio are carried out in bootstrap experiments. The UTR replicates the surprise factor best but is most uncertain and only comparable between records of similar length. For samples with symmetric Lorenz curves, shape parameter, obesity and Gini indices provide consistent indications. For asymmetric Lorenz curves, however, the first two tend to overestimate, whereas Gini index tends to underestimate tail heaviness. We suggest the use of a combination of shape parameter, obesity and Gini index to characterize tail heaviness. These indicators should be supported with calculation of the Lorenz asymmetry coefficients and interpreted with caution.
The sensitivity of key hydrologic variables and hydropower generation to climate change in the Lake Malawi and Shire River basins is assessed. The study adapts the mesoscale Hydrological Model (mHM) which is applied separately in the Upper Lake Malawi and Shire River basins. A particular Lake Malawi model, which focuses on reservoir routing and lake water balance, has been developed and is interlinked between the two basins. Climate change projections from 20 Coordinated Regional Climate Downscaling Experiment (CORDEX) models for Africa based on two scenarios (RCP4.5 and RCP8.5) for the periods 2021-2050 and 2071-2100 are used. An annual temperature increase of 1 degrees C decreases mean lake level and outflow by 0.3 m and 17%, respectively, signifying the importance of intensified evaporation for Lake Malawi's water budget. Meanwhile, a +5% (-5%) deviation in annual rainfall changes mean lake level by +0.7 m (-0.6 m). The combined effects of temperature increase and rainfall decrease result in significantly lower flows in the Shire River. The hydrological river regime may change from perennial to seasonal with the combination of annual temperature increase and precipitation decrease beyond 1.5 degrees C (3.5 degrees C) and -20% (-15%). The study further projects a reduction in annual hydropower production between 1% (RCP8.5) and 2.5% (RCP4.5) during 2021-2050 and between 5% (RCP4.5) and 24% (RCP8.5) during 2071-2100. The results show that it is of great importance that a further development of hydro energy on the Shire River should take into account the effects of climate change, e.g., longer low flow periods and/or higher discharge fluctuations, and thus uncertainty in the amount of electricity produced.
The occurrence of refugia beyond the arctic treeline and genetic adaptation therein play a crucial role of largely unknown effect size. While refugia have potential for rapidly colonizing the tundra under global warming, the taxa may be maladapted to the new environmental conditions. Understanding the genetic composition and age of refugia is thus crucial for predicting any migration response.
Here, we genotype 194 larch individuals from an similar to 1.8 km(2)area in northcentral Siberia on the southern Taimyr Peninsula by applying an assay of 16 nuclear microsatellite markers. For estimating the age of clonal individuals, we counted tree rings at sections along branches to establish a lateral growth rate that was then combined with geographic distance.
Findings reveal that the predominant reproduction type is clonal (58.76%) by short distance spreading of ramets. One outlier of clones 1 km apart could have been dispersed by reindeer. In clonal groups and within individuals, we find that somatic mutations accumulate with geographic distance. Clonal groups of two or more individuals are observed. Clonal age estimates regularly suggest individuals as old as 2,200 years, which coincides with a major environmental change that forced a treeline retreat in the region.
We conclude that individuals with clonal growth mode were naturally selected as it lowers the likely risk of extinction under a harsh environment. We discuss this legacy from the past that might now be a maladaptation and hinder expansion under currently strongly increasing temperatures.
Variation of deuterium excess in surface waters across a 5000-m elevation gradient in eastern Nepal
(2020)
The strong elevation gradient of the Himalaya allows for investigation of altitude and orographic impacts on surface water delta O-18 and delta D stable isotope values. This study differentiates the time- and altitude-variable contributions of source waters to the Arun River in eastern Nepal. It provides isotope data along a 5000-m gradient collected from tributaries as well as groundwater, snow, and glacial-sourced surface waters and time-series data from April to October 2016. We find nonlinear trends in delta O-18 and delta D lapse rates with high-elevation lapse rates (4000-6000 masl) 5-7 times more negative than low-elevation lapse rates (1000-3000 masl). A distinct seasonal signal in delta O-18 and delta D lapse rates indicates time-variable source-water contributions from glacial and snow meltwater as well as precipitation transitions between the Indian Summer Monsoon and Winter Westerly Disturbances. Deuterium excess correlates with the extent of snowpack and tracks melt events during the Indian Summer Monsoon season. Our analysis identifies the influence of snow and glacial melt waters on river composition during low-flow conditions before the monsoon (April/May 2016) followed by a 5-week transition to the Indian Summer Monsoon-sourced rainfall around mid-June 2016. In the post-monsoon season, we find continued influence from glacial melt waters as well as ISM-sourced groundwater.
Projected changes in compound flood hazard from riverine and coastal floods in northwestern Europe
(2020)
Compound flooding in coastal regions, that is, the simultaneous or successive occurrence of high sea levels and high river flows, is expected to increase in a warmer world. To date, however, there is no robust evidence on projected changes in compound flooding for northwestern Europe. We combine projected storm surges and river floods with probabilistic, localized relative sea-level rise (SLR) scenarios to assess the future compound flood hazard over northwestern coastal Europe in the high (RCP8.5) emission scenario. We use high-resolution, dynamically downscaled regional climate models (RCM) to drive a storm surge model and a hydrological model, and analyze the joint occurrence of high coastal water levels and associated river peaks in a multivariate copula-based approach. The RCM-forced multimodel mean reasonably represents the observed spatial pattern of the dependence strength between annual maxima surge and peak river discharge, although substantial discrepancies exist between observed and simulated dependence strength. All models overestimate the dependence strength, possibly due to limitations in model parameterizations. This bias affects compound flood hazard estimates and requires further investigation. While our results suggest decreasing compound flood hazard over the majority of sites by 2050s (2040-2069) compared to the reference period (1985-2005), an increase in projected compound flood hazard is limited to around 34% of the sites. Further, we show the substantial role of SLR, a driver of compound floods, which has frequently been neglected. Our findings highlight the need to be aware of the limitations of the current generation of Earth system models in simulating coastal compound floods.
Neodymium isotopic composition (epsilon Nd) has enjoyed widespread use as a palaeotracer, principally because it behaves quasi-conservatively in the modern ocean. However, recent bottom water epsilon Nd reconstructions from the eastern North Atlantic are difficult to interpret under assumptions of conservative behaviour. The observation that this apparent departure from conservative behaviour increases with enhanced ice-rafted debris (IRD) fluxes has resulted in the suggestion that IRD leads to the overprinting of bottom water epsilon Nd through reversible scavenging. In this study, a simple water column model successfully reproduces epsilon Nd reconstructions from the eastern North Atlantic at the Last Glacial Maximum and Heinrich Stadial 1, and demonstrates that the changes in scavenging intensity required for good model-data fit is in good agreement with changes in the observed IRD flux. Although uncertainties in model parameters preclude a more definitive conclusion, the results indicate that the suggestion of IRD as a source of non-conservative behaviour in the epsilon Nd tracer is reasonable and that further research into the fundamental chemistry underlying the marine neodymium cycle is necessary to increase confidence in assumptions of conservative epsilon Nd behaviour in the past.
Lakes cover large parts of the climatically sensitive Arctic landscape and respond rapidly to environmental change. Arctic lakes have different origins and include the predominant thermokarst lakes, which are small, young and highly dynamic, as well as large, old and stable glacial lakes. Freshwater diatoms dominate the primary producer community in these lakes and can be used to detect biotic responses to climate and environmental change. We used specific diatom metabarcoding on sedimentary DNA, combined with next-generation sequencing and diatom morphology, to assess diatom diversity in five glacial and 15 thermokarst lakes within the easternmost expanse of the Siberian treeline ecotone in Chukotka, Russia. We obtained 163 verified diatom sequence types and identified 176 diatom species morphologically. Although there were large differences in taxonomic assignment using the two approaches, they showed similar high abundances and diversity of Fragilariceae and Aulacoseiraceae. In particular, the genetic approach detected hidden within-lake variations of fragilarioids in glacial lakes and dominance of centric Aulacoseira species, whereas Lindavia ocellata was predominant using morphology. In thermokarst lakes, sequence types and valve counts also detected high diversity of Fragilariaceae, which followed the vegetation gradient along the treeline. Ordination analyses of the genetic data from glacial and thermokarst lakes suggest that concentrations of sulfate (SO42-), an indicator of the activity of sulfate-reducing microbes under anoxic conditions, and bicarbonate (HCO3-), which relates to surrounding vegetation, have a significant influence on diatom community composition. For thermokarst lakes, we also identified lake depth as an important variable, but SO42- best explains diatom diversity derived from genetic data, whereas HCO3- best explains the data from valve counts. Higher diatom diversity was detected in glacial lakes, most likely related to greater lake age and different edaphic settings, which gave rise to diversification and endemism. In contrast, small, dynamic thermokarst lakes are inhabited by stress-tolerant fragilarioids and are related to different vegetation types along the treeline ecotone. Our study demonstrated that genetic investigations of lake sediments can be used to interpret climate and environmental responses of diatoms. It also showed how lake type affects diatom diversity, and that such genetic analyses can be used to track diatom community changes under ongoing warming in the Arctic.
Multi-year index-based insurance for adapting Water Utility Companies to hydrological drought
(2020)
The sustainability of water utility companies is threatened by non-stationary drivers, such as climate and anthropogenic changes. To cope with potential economic losses, instruments such as insurance are useful for planning scenarios and mitigating impacts, but data limitations and risk uncertainties affect premium estimation and, consequently, business sustainability. This research estimated the possible economic impacts of business interruption to the Sao Paulo Water Utility Company derived from hydrological drought and how this could be mitigated with an insurance scheme. Multi-year insurance (MYI) was proposed through a set of "change" drivers: the climate driver, through forcing the water evaluation and planning system (WEAP) hydrological tool; the anthropogenic driver, through water demand projections; and the economic driver, associated with recent water price policies adopted by the utility company during water scarcity periods. In our study case, the evaluated indices showed that MYI contracts that cover only longer droughts, regardless of the magnitude, offer better financial performance than contracts that cover all events (in terms of drought duration). Moreover, through MYI contracts, we demonstrate solvency for the insurance fund in the long term and an annual average actuarially fair premium close to the total expected revenue reduction.
Increasing interests in hydrocarbon resources at depths have drawn greater attentions to the deeply-buried carbonate reservoirs in the Tarim Basin in China. In this study, the cyclic dolomite rocks of Upper Cambrian Lower Qiulitag Group from four outcrop sections in northwestern Tarim Basin were selected to investigate and evaluate the petrophysical properties in relation to depositional facies and cyclicity. The Lower Qiulitag Group includes ten lithofacies, which were deposited in intermediate to shallow subtidal, restricted shallow subtidal, intertidal, and supratidal environments on a carbonate ramp system. These lithofacies are vertically stacked into repeated shallowing-upward, meter-scale cycles which are further grouped into six third-order depositional sequences (Sq1 to Sq6). There are variable types of pore spaces in the Lower Qiulitag Group dolomite rocks, including interparticle, intraparticle, and fenestral pores of primary origin, inter crystal, and vuggy pores of late diagenetic modification. The porosity in the dolomites is generally facies-selective as that the microbially-originated thrombolites and stromatolites generally yield a relatively high porosity. In contrast, the high-energy ooidal grainstones generally have very low porosity. In this case, the microbialite-based peritidal cycles and peritidal cycle-dominated highstand (or regressive) successions have relatively high volumes of pore spaces, although highly fluctuating (or vertical inhomogeneous). Accordingly, the grainstone-based subtidal cycles and subtidal cycle-dominated transgressive successions generally yield extremely low porosity. This scenario indicates that porosity development and preservation in the thick dolomite successions are primarily controlled by depositional facies which were influenced by sea-level fluctuations of different orders and later diagenetic overprinting.
Understanding the physical mechanisms governing fluid-induced fault slip is important for improved mitigation of seismic risks associated with large-scale fluid injection. We conducted fluid-induced fault slip experiments in the laboratory on critically stressed saw-cut sandstone samples with high permeability using different fluid pressurization rates. Our experimental results demonstrate that fault slip behavior is governed by fluid pressurization rate rather than injection pressure. Slow stick-slip episodes (peak slip velocity < 4 mu m/s) are induced by fast fluid injection rate, whereas fault creep with slip velocity < 0.4 mu m/s mainly occurs in response to slow fluid injection rate. Fluid-induced fault slip may remain mechanically stable for loading stiffness larger than fault stiffness. Independent of fault slip mode, we observed dynamic frictional weakening of the artificial fault at elevated pore pressure. Our observations highlight that varying fluid injection rates may assist in reducing potential seismic hazards of field-scale fluid injection projects. <br /> Plain Language Summary Human-induced earthquakes from field-scale fluid injection projects including enhanced geothermal system and deep wastewater injection have been documented worldwide. Although it is clear that fluid pressure plays a crucial role in triggering fault slip, the physical mechanism behind induced seismicity still remains poorly understood. We performed laboratory tests, and here we present two fluid-induced slip experiments conducted on permeable Bentheim sandstone samples crosscut by a fault that is critically stressed. Fault slip is then triggered by pumping the water from the bottom end of the sample at different fluid injection rates. Our results show that fault slip is controlled by fluid pressure increase rate rather than by the absolute magnitude of fluid pressure. In contrast to episodes of relatively rapid but stable sliding events caused by a fast fluid injection rate, fault creep is observed during slow fluid injection. Strong weakening of the dynamic friction coefficient of the experimental fault is observed at elevated pore pressure, independent of fault slip mode. These results may provide a better understanding of the complex behavior of fluid-induced fault slip on the field scale.
Crustal seismic attenuation parameters in the western region of the North Anatolian Fault Zone
(2020)
Detailed knowledge of the crustal structure along the North Anatolian Fault Zone can help in understanding past and present tectonic processes in relation to the deformation history. To estimate the frequency-dependent crustal attenuation parameters beneath the western part of the North Anatolian Fault Zone we apply acoustic radiative transfer theory under the assumption of multiple isotropic scattering to generate synthetic seismogram envelopes. The inversion depends on finding an optimal fit between observed and synthetically computed coda wave envelopes in five frequency bands. 2-D lateral variation of intrinsic and scattering attenuation at various frequencies tends to three crustal blocks (i.e., Armutlu-Almacik, Istanbul-Zonguldak and Sakarya Zones) separated by the southern and northern branches of the western part of the North Anatolian Fault Zone. Overall, scattering attenuation appears to be dominant over intrinsic attenuation in the study area at lower frequencies. Relatively low attenuation properties are observed beneath the older Istanbul Zone whereas higher attenuation properties are found for the younger Sakarya Zone. The Armutlu Almacik Zone exhibits more complex lateral variations. Very high attenuation values towards the west characterize the area of the Kuzuluk Basin, a pull-apart basin formed under west-east extension. Our coda-derived moment magnitudes are similar to the local magnitude estimates that were previously calculated for the same earthquakes. For smaller earthquakes (M-L < 2.5), however, the relation between local and moment magnitudes appears to lose its coherency. This may stem from various reasons including the use of seismic data recorded in finite sampling interval, possible biases in local magnitude estimates of earthquake catalogues as well as biases due to wrong assumptions to consider anelastic attenuation terms.
Holocene temperature proxy records are commonly used in quantitative synthesis and model-data comparisons. However, comparing correlations between time series from records collected in proximity to one another with the expected correlations based on climate model simulations indicates either regional or noisy climate signals in Holocene temperature proxy records. In this study, we evaluate the consistency of spatial correlations present in Holocene proxy records with those found in data from the Last Glacial Maximum (LGM). Specifically, we predict correlations expected in LGM proxy records if the only difference to Holocene correlations would be due to more time uncertainty and more climate variability in the LGM. We compare this simple prediction to the actual correlation structure in the LGM proxy records. We found that time series data of ice-core stable isotope records and planktonic foraminifera Mg/Ca ratios were consistent between the Holocene and LGM periods, while time series of Uk'37 proxy records were not as we found no correlation between nearby LGM records. Our results support the finding of highly regional or noisy marine proxy records in the compilation analysed here and suggest the need for further studies on the role of climate proxies and the processes of climate signal recording and preservation.
The Fram Strait is an area with a relatively low and irregular distribution of diatom microfossils in surface sediments, and thus microfossil records are scarce, rarely exceed the Holocene, and contain sparse information about past richness and taxonomic composition. These attributes make the Fram Strait an ideal study site to test the utility of sedimentary ancient DNA (sedaDNA) metabarcoding. Amplifying a short, partial rbcL marker from samples of sediment core MSM05/5-712-2 resulted in 95.7% of our sequences being assigned to diatoms across 18 different families, with 38.6% of them being resolved to species and 25.8% to genus level. Independent replicates show a high similarity of PCR products, especially in the oldest samples. Diatom sedaDNA richness is highest in the Late Weichselian and lowest in Mid- and Late Holocene samples. Taxonomic composition is dominated by cold-water and sea-ice-associated diatoms and suggests several reorganisations - after the Last Glacial Maximum, after the Younger Dryas, and after the Early and after the Mid-Holocene. Different sequences assigned to, amongst others, Chaetoceros socialis indicate the detectability of intra-specific diversity using sedaDNA. We detect no clear pattern between our diatom sedaDNA record and the previously published IP25 record of this core, although proportions of pennate diatoms increase with higher IP25 concentrations and proportions of Nitzschia cf. frigida exceeding 2% of the assemblage point towards past sea-ice presence.
The structural and topographic evolution of orogenic plateaus is an important research topic because of its impact on atmospheric circulation patterns, the amount and distribution of rainfall, and resulting changes in surface processes. The Puna region in the north-western Argentina (between 13 degrees S and 27 degrees S) is part of the Andean Plateau, which is the world's second largest orogenic plateau. In order to investigate the deformational events responsible for the initial growth of this part of the Andean plateau, we carried out structural and stratigraphic investigations within the present-day transition zone between the northern Puna and the adjacent Eastern Cordillera to the east. This transition zone is characterized by ubiquitous exposures of continental middle Eocene redbeds of the Casa Grande Formation. Our structural mapping, together with a sedimentological analysis of these units and their relationships with the adjacent mountain ranges, has revealed growth structures and unconformities that are indicative of syntectonic deposition. These findings support the notion that tectonic shortening in this part of the Central Andes was already active during the middle Paleogene, and that early Cenozoic deformation in the region that now constitutes the Puna occurred in a spatially irregular manner. The patterns of Paleogene deformation and uplift along the eastern margin of the present-day plateau correspond to an approximately north-south oriented swath of reactivated basement heterogeneities (i.e. zones of mechanical weakness) stemming from regional Paleozoic mountain building that may have led to local concentration of deformation belts.
The Indus Molasse records orogenic sedimentation associated with uplift and erosion of the southern margin of Asia in the course of ongoing India-Eurasia collision. Detailed field investigation clarifies the nature and extent of the depositional contact between this molasse and the underlying basement units. We report the first dataset on detrital zircon U-Pb ages, Hf isotopes and apatite U-Pb ages for the autochthonous molasse in the Indus Suture Zone. A latest Oligocene depositional age is proposed on the basis of the youngest detrital zircon U-Pb age peak and is consistent with published biostratigraphic data. Multiple provenance indicators suggest exclusively northerly derivation with no input from India in the lowermost parts of the section. The results provide constraints on the uplift and erosion history of the Ladakh Range following the initial India-Asia collision.
The formation of the supercontinent Pangaea during the Permo-Triassic gave rise to an extreme monsoonal climate (often termed 'mega-monsoon') that has been documented by numerous palaeo-records. However, considerable debate exists about the role of orbital forcing in causing humid intervals in an otherwise arid climate. To shed new light on the forcing of monsoonal variability in subtropical Pangaea, this study focuses on sediment facies and colour variability of playa and alluvial fan deposits in an outcrop from the late Carnian (ca 225 Ma) in the southern Germanic Basin, south-western Germany. The sediments were deposited against a background of increasingly arid conditions following the humid Carnian Pluvial Event (ca 234 to 232 Ma). The ca 2 center dot 4 Myr long sedimentary succession studied shows a tripartite long-term evolution, starting with a distal mud-flat facies deposited under arid conditions. This phase was followed by a highly variable playa-lake environment that documents more humid conditions and finally a regression of the playa-lake due to a return of arid conditions. The red-green (a*) and lightness (L*) records show that this long-term variability was overprinted by alternating wet/dry cycles driven by orbital precession and ca 405 kyr eccentricity, without significant influence of obliquity. The absence of obliquity in this record indicates that high-latitude forcing played only a minor role in the southern Germanic Basin during the late Carnian. This is different from the subsequent Norian when high-latitude signals became more pronounced, potentially related to the northward drift of the Germanic Basin. The recurring pattern of pluvial events during the late Triassic demonstrates that orbital forcing, in particular eccentricity, stimulated the occurrence and intensity of wet phases. It also highlights the possibility that the Carnian Pluvial Event, although most likely triggered by enhanced volcanic activity, may also have been modified by an orbital stimulus.
Cosmic-Ray neutron sensors are widely used to determine soil moisture on the hectare scale. Precise measurements, especially in the case of mobile application, demand for neutron detectors with high counting rates and high signal-to-noise ratios. For a long time Cosmic Ray Neutron Sensing (CRNS) instruments have relied on He-3 as an efficient neutron converter. Its ongoing scarcity demands for technological solutions using alternative converters, which are Li-6 and B-10. Recent developments lead to a modular neutron detector consisting of several B-10-lined proportional counter tubes, which feature high counting rates via its large surface area. The modularity allows for individual shieldings of different segments within the detector featuring the capability of gaining spectral information about the detected neutrons. This opens the possibility for active signal correction, especially useful when applied to mobile measurements, where the influence of constantly changing near-field to the overall signal should be corrected. Furthermore, the signal-to-noise ratio could be increased by combining pulse height and pulse length spectra to discriminate between neutrons and other environmental radiation. This novel detector therefore combines high-selective counting electronics with large-scale instrumentation technology.
The Victoria microplate between the Eastern and Western Branches of the East African Rift System is one of the largest continental microplates on Earth. In striking contrast to its neighboring plates, Victoria rotates counterclockwise with respect to Nubia. The underlying cause of this distinctive rotation has remained elusive so far. Using 3D numerical models, we investigate the role of pre-existing lithospheric heterogeneities in continental microplate rotation. We find that Victoria's rotation is primarily controlled by the distribution of rheologically stronger zones that transmit the drag of the major plates to the microplate and of the mechanically weaker mobile belts surrounding Victoria that facilitate rotation. Our models reproduce Victoria's GPS-derived counterclockwise rotation as well as key complexities of the regional tectonic stress field. These results reconcile competing ideas on the opening of the rift system by highlighting differences in orientation of the far-field divergence, local extension, and the minimum horizontal stress. One of the largest continental microplates on Earth is situated in the center of the East African Rift System, and oddly, the Victoria microplate rotates counterclockwise with respect to the neighboring African tectonic plate. Here, the authors' modelling results suggest that Victoria microplate rotation is caused by edge-driven lithospheric processes related to the specific geometry of rheologically weak and strong regions.
Analysis of past and present stimulation projects reveals that the temporal evolution and growth of maximum observed moment magnitudes may be linked directly to the injected fluid volume and hydraulic energy. Overall evolution of seismic moment seems independent of the tectonic stress regime and is most likely governed by reservoir specific parameters, such as the preexisting structural inventory. Data suggest that magnitudes can grow either in a stable way, indicating the constant propagation of self-arrested ruptures, or unbound, for which the maximum magnitude is only limited by the size of tectonic faults and fault connectivity. Transition between the two states may occur at any time during injection or not at all. Monitoring and traffic light systems used during stimulations need to account for the possibility of unstable rupture propagation from the very beginning of injection by observing the entire seismicity evolution in near-real time and at high resolution for an immediate reaction in injection strategy.
Plain Language Summary Predicting and controlling the size of earthquakes caused by fluid injection is currently the major concern of many projects associated with geothermal energy production. Here, we analyze the magnitude and seismic moment evolution with injection parameters for prominent geothermal and scientific projects to date. Evolution of seismicity seems to be largely independent of the tectonic stress background and seemingly depends on reservoir specific characteristics. We find that the maximum observed magnitudes relate linearly to the injected volume or hydraulic energy. A linear relation suggests stable growth of induced ruptures, as predicted by current models, or rupture growth may no longer depend on the stimulated volume but on tectonics. A system may change between the two states during the course of fluid injection. Close-by and high-resolution monitoring of seismic and hydraulic parameters in near-real time may help identify these fundamental changes in ample time to change injection strategy and manage maximum magnitudes.
Other than commonly assumed in seismology, the phase velocity of Rayleigh waves is not necessarily a single-valued function of frequency. In fact, a single Rayleigh mode can exist with three different values of phase velocity at one frequency. We demonstrate this for the first higher mode on a realistic shallow seismic structure of a homogeneous layer of unconsolidated sediments on top of a half-space of solid rock (LOH). In the case of LOH a significant contrast to the half-space is required to produce the phenomenon. In a simpler structure of a homogeneous layer with fixed (rigid) bottom (LFB) the phenomenon exists for values of Poisson's ratio between 0.19 and 0.5 and is most pronounced for P-wave velocity being three times S-wave velocity (Poisson's ratio of 0.4375). A pavement-like structure (PAV) of two layers on top of a half-space produces the multivaluedness for the fundamental mode. Programs for the computation of synthetic dispersion curves are prone to trouble in such cases. Many of them use mode-follower algorithms which loose track of the dispersion curve and miss the multivalued section. We show results for well established programs. Their inability to properly handle these cases might be one reason why the phenomenon of multivaluedness went unnoticed in seismological Rayleigh wave research for so long. For the very same reason methods of dispersion analysis must fail if they imply wave number k(l)(omega) for the lth Rayleigh mode to be a single-valued function of frequency.. This applies in particular to deconvolution methods like phase-matched filters. We demonstrate that a slant-stack analysis fails in the multivalued section, while a Fourier-Bessel transformation captures the complete Rayleigh-wave signal. Waves of finite bandwidth in the multivalued section propagate with positive group-velocity and negative phase-velocity. Their eigenfunctions appear conventional and contain no conspicuous feature.
Rivers transfer particulate organic carbon (POC) from eroding mountains into geological sinks. Organic carbon source composition and selective mobilization have been shown to affect the type and quantity of POC export, but their combined effects across complex mountain ranges remain underexplored. Here, we examine the variation in organic carbon sourcing and transport in the trans-Himalayan Kali Gandaki River catchment, along strong gradients in precipitation, rock type and vegetation. Combining bulk stable nitrogen, and stable and radioactive organic carbon isotopic composition of bedrock, litter, soil and river sediment samples with a Bayesian end-member mixing approach, we differentiate POC sources along the river and quantify their export. Our analysis shows that POC export from the Tibetan segment of the catchment, where carbon bearing shales are partially covered by aged and modern soils, is dominated by petrogenic POC. Based on our data we re-assess the presence of aged biospheric OC in this part of the catchment, and its contribution to the river load. In the High Himalayan segment, we observed low inputs of petrogenic and biospheric POC, likely due to very low organic carbon concentrations in the metamorphic bedrock, combined with erosion dominated by deep-seated landslides. Our findings show that along the Kali Gandaki River, the sourcing of sediment and organic carbon are decoupled, due to differences in rock organic carbon content, soil and above ground carbon stocks, and geomorphic process activity. While the fast eroding High Himalayas are the principal source of river sediment, the Tibetan headwaters, where erosion rates are lower, are the principal source of organic carbon. To robustly estimate organic carbon export from the Himalayas, the mountain range should be divided into tectono-physiographic zones with distinct organic carbon yields due to differences in substrate and erosion processes and rates.
The regional patterns and timing of the Younger Dryas cooling in the North Atlantic realm were complex and are mechanistically incompletely understood. To enhance understanding of regional climate patterns, we present molecular biomarker records at subannual to annual resolution by mass spectrometry imaging (MSI) of sediments from the Lake Meerfelder Maar covering the Allerod-Younger Dryas transition. These analyses are supported by conventional extraction-based molecular-isotopic analyses, which both validate the imaging results and constrain the sources of the target compounds. The targeted fatty acid biomarkers serve as a gauge of the response of the local aquatic and terrestrial ecosystem to climate change. Based on the comparison of our data with existing data from Meerfelder Maar, we analyse the short-term environmental evolution in Western Europe during the studied time interval and confirm the previously reported delayed hydrological response to Greenland cooling. However, despite a detected delay of Western European environmental change of similar to 135 years, our biomarker data show statistically significant correlation with deuterium excess in Greenland ice core at - annual resolution during this time-transgressive cooling. This suggests a coherent atmospheric forcing across the North Atlantic realm during this transition. We propose that Western European cooling was postponed due to major reorganization of the westerlies that were intermittently forcing warmer and wetter air masses from lower latitudes to Western Europe and thus resulted in delayed cooling relative to Greenland.
(40)A/Ar-39 step-heating of mica and amphibole megacrysts from hauyne-bearing olivine melilitite scoria/tephra from the Zelezna hurka yielded a 435 +/- 108 ka isotope correlation age for phlogopite and a more imprecise 1.55 Ma total gas age of the kaersutite megacryst. The amphibole megacrysts may constitute the first, and the younger phlogopite megacrysts the later phase of mafic, hydrous melilitic magma crystallization. It cannot be ruled out that the amphibole megacrysts are petrogenetically unrelated to tephra and phlogopite megacrysts and were derived from mantle xenoliths or disaggregated older, deep crustal pegmatites. This is in line both with the rarity of amphibole at Zelezna hurka and with the observed signs of magmatic resorption at the edges of amphibole crystals.
Two different approaches are used to assess the impacts associated with natural hazards and climate change in cities. A bottom-up approach uses high resolution data on constituent assets within the urban area. In contrast, a top-down approach uses less detailed information but is consequently more readily transferable. Here, we compare damage curves generated by each approach for coastal flooding in London. To compare them, we fit a log-logistic regression with three parameters to the calculated damage curves. We find that the functions are remarkably similar in their shape, albeit with different inflection points and a maximum damage that differs by 13%-25%. If rescaled, the curves agree almost exactly, which enables damage assessment to be undertaken following the calculation of the three parameters.
Initiation of subduction following the impingement of a hot buoyant mantle plume is one of the few scenarios that allow breaking the lithosphere and recycling a stagnant lid without requiring any preexisting weak zones. Here, we investigate factors controlling the number and shape of retreating subducting slabs formed by plume-lithosphere interaction. Using 3-D thermomechanical models we show that the deformation regime, which defines formation of single-slab or multi-slab subduction, depends on several parameters such as age of oceanic lithosphere, thickness of the crust and large-scale lithospheric extension rate. Our model results indicate that on present-day Earth multi-slab plume-induced subduction is initiated only if the oceanic lithosphere is relatively young (<30-40 Myr, but >10 Myr), and the crust has a typical thickness of 8 km. In turn, development of single-slab subduction is facilitated by older lithosphere and pre-imposed extensional stresses. In early Earth, plume-lithosphere interaction could have led to formation of either episodic short-lived circular subduction when the oceanic lithosphere was young or to multi-slab subduction when the lithosphere was old.
Knickpoints in longitudinal river profiles are proxies for the climatic and tectonic history of active mountains. The analysis of river profiles commonly relies on the assumption that drainage network configurations are stable. Here, we show that this assumption must be made cautiously if changes in contributing area are fast relative to knickpoint migration rates. We studied the Parachute Creek basin in the Roan Plateau, Colorado, United States, where knickpoint retreat occurs in horizontally uniform lithology so that drainage area is the sole governing variable. In this basin, we identified an anomalous catchment in the degree to which a stream power-based model predicted knickpoint locations. The catchment is experiencing area loss as the plateau edge is eroded by cliff migration in proximity to the Colorado River. Model predictions improve if the plateau edge is assumed to have migrated over the time scale of knickpoint retreat. Finally, a Lagrangian model of knickpoint migration enabled us to study the kinematic links between drainage area loss and knickpoint migration and offered constraints on the temporal aspects of area loss. Modeled onset and amount of area loss are consistent with cliff retreat rates along the margin of the Roan Plateau inferred from the incisional history of the upper Colorado River.
The Arctic ecosystem, a reservoir of genetic microbial diversity, represents a virtually unlimited source of microorganisms that could interact with human beings. Despite continuous exploration of Arctic habitats and description of their microbial communities, bacterial phenotypes commonly associated with pathogenicity, such as hemolytic activity, have rarely been reported. In this study, samples of snow, fresh and marine water, soil, and sediment from several habitats in the Arctic archipelago of Svalbard were collected during Summer, 2017. Bacterial isolates were obtained after incubation on oligotrophic media at different temperatures and their hemolytic potential was assessed on sheep blood agar plates. Partial (alpha) or true (beta) hemolysis was observed in 32 out of 78 bacterial species. Genes expressing cytolytic compounds, such as hemolysins, likely increase the general fitness of the producing microorganisms and confer a competitive advantage over the availability of nutrients in natural habitats. In environmental species, the nutrient-acquisition function of these compounds presumably precedes their function as toxins for mammalian erythrocytes. However, in the light of global warming, the presence of hemolytic bacteria in Arctic environments highlights the possible risks associated with these microorganisms in the event of habitat melting/destruction, ecosystem transition, and re-colonization.
Full-waveform-based characterization of acoustic emission activity in a mine-scale experiment
(2020)
Understanding fracturing processes and the hydromechanical relation to induced seismicity is a key question for enhanced geothermal systems (EGS). Commonly massive fluid injection, predominately causing hydroshearing, are used in large-scale EGS but also hydraulic fracturing approaches were discussed. To evaluate the applicability of hydraulic fracturing techniques in EGS, six in situ, multistage hydraulic fracturing experiments with three different injection schemes were performed under controlled conditions in crystalline rock at the Aspo Hard Rock Laboratory (Sweden). During the experiments the near-field ground motion was continuously recorded by 11 piezoelectric borehole sensors with a sampling rate of 1 MHz. The sensor network covered a volume of 30x30x30 m around a horizontal, 28-m-long injection borehole at a depth of 410 m. To extract and characterize massive, induced, high-frequency acoustic emission (AE) activity from continuous recordings, a semi-automated workflow was developed relying on full waveform based detection, classification and location procedures. The approach extended the AE catalogue from 196 triggered events in previous studies to more than 19600 located AEs. The enhanced catalogue, for the first time, allows a detailed analysis of induced seismicity during single hydraulic fracturing experiments, including the individual fracturing stages and the comparison between injection schemes. Beside the detailed study of the spatio-temporal patterns, event clusters and the growth of seismic clouds, we estimate relative magnitudes and b-values of AEs for conventional, cyclic progressive and dynamic pulse injection schemes, the latter two being fatigue hydraulic fracturing techniques. While the conventional fracturing leads to AE patterns clustered in planar regions, indicating the generation of a single main fracture plane, the cyclic progressive injection scheme results in a more diffuse, cloud-like AE distribution, indicating the activation of a more complex fracture network. For a given amount of hydraulic energy (pressure multiplied by injected volume) pumped into the system, the cyclic progressive scheme is characterized by a lower rate of seismicity, lower maximum magnitudes and significantly larger b-values, implying an increased number of small events relative to the large ones. To our knowledge, this is the first direct comparison of high resolution seismicity in a mine-scale experiment induced by different hydraulic fracturing schemes.
The change of the mechanical properties of granular materials with pressure is an important topic associated with many industrial applications. In this paper we investigate the influence of hydrostatic pressure (P-e) on the effective bulk compressibility (C-eff) of a granular material by applying two modified theoretical approaches that are based on contact mechanics and micromechanics, respectively. For a granular material composed of rough grains, an extended contact model is developed to elucidate the effect of roughness of grain surfaces on bulk compressibility. At relatively low pressures, the model predicts that the decrease of bulk compressibility with pressure may be described by a power law with an exponent of -1/2 (i.e., C-eff proportional to P-e(1/2) ), but deviates at intermediate pressures. At elevated pressures beyond full contact, bulk compressibility remains almost unchanged, which may be roughly evaluated by continuum contact mechanics. As an alternative explanation of pressure-dependent bulk compressibility, we suggest a micromechanical model that accounts for effects of different types of pore space present in granular materials. Narrow and compliant inter-granular cracks are approximated by three-dimensional oblate spheroidal cracks with rough surfaces, whereas the equant and stiff pores surrounded by three and four neighboring grains are modeled as tubular pores with cross sections of three and four cusp-like corners, respectively. In this model, bulk compressibility is strongly reduced with increasing pressure by progressive closure of rough-walled cracks. At pressures exceeding crack closure pressure, deformation of the remaining equant pores is largely insensitive to pressure, with almost no further change in bulk compressibility. To validate these models, we performed hydrostatic compression tests on Bentheim sandstone (a granular rock consisting of quartz with high porosity) under a wide range of pressure. The relation between observed microstructures and measured pressure-dependent bulk compressibility is well explained by both suggested models.
Cities will play a key role in the grand challenge of nourishing a growing global population, because, due to their population density, they set the demand. To ensure that food systems are sustainable, as well as nourishing, one solution often suggested is to shorten their supply chains toward a regional rather than a global basis. While such regional systems may have a range of costs and benefits, we investigate the mitigation potential of regionalized urban food systems by examining the greenhouse gas emissions associated with food transport. Using data on food consumption for 7108 urban administrative units (UAUs), we simulate total transport emissions for both regionalized and globalized supply chains. In regionalized systems, the UAUs' demands are fulfilled by peripheral food production, whereas to simulate global supply chains, food demand is met from an international pool (where the origin can be any location globally). We estimate that regionalized systems could reduce current emissions from food transport. However, because longer supply chains benefit from maximizing comparative advantage, this emission reduction would require closing yield gaps, reducing food waste, shifting toward diversified farming, and consuming seasonal produce. Regionalization of food systems will be an essential component to limit global warming to well below 2 degrees C in the future.
According to established understanding, deep-water formation in the North Atlantic and Southern Ocean keeps the deep ocean cold, counter-acting the downward mixing of heat from the warmer surface waters in the bulk of the world ocean. Therefore, periods of strong Atlantic meridional overturning circulation (AMOC) are expected to coincide with cooling of the deep ocean and warming of the surface waters. It has recently been proposed that this relation may have reversed due to global warming, and that during the past decades a strong AMOC coincides with warming of the deep ocean and relative cooling of the surface, by transporting increasingly warmer waters downward. Here we present multiple lines of evidence, including a statistical evaluation of the observed global mean temperature, ocean heat content, and different AMOC proxies, that lead to the opposite conclusion: even during the current ongoing global temperature rise a strong AMOC warms the surface. The observed weakening of the AMOC has therefore delayed global surface warming rather than enhancing it.
Social Media Abstract:
The overturning circulation in the Atlantic Ocean has weakened in response to global warming, as predicted by climate models. Since it plays an important role in transporting heat, nutrients and carbon, a slowdown will affect global climate processes and the global mean temperature. Scientists have questioned whether this slowdown has worked to cool or warm global surface temperatures. This study analyses the overturning strength and global mean temperature evolution of the past decades and shows that a slowdown acts to reduce the global mean temperature. This is because a slower overturning means less water sinks into the deep ocean in the subpolar North Atlantic. As the surface waters are cold there, the sinking normally cools the deep ocean and thereby indirectly warms the surface, thus less sinking implies less surface warming and has a cooling effect. For the foreseeable future, this means that the slowing of the overturning will likely continue to slightly reduce the effect of the general warming due to increasing greenhouse gas concentrations.
Many of the volcanic plateau margins of the eastern, formerly glaciated, foreland of the Patagonian Andes are undermined by giant landslides (>= 10(8) m(3)). One cluster of such landslides extends along the margin of the Meseta del Lago Buenos Aires (MLBA) plateau that is formed mainly by Neogene-Quaternary basalts. The dry climate is at odds with numerous >2-km long earthflows nested within older and larger compound landslides. We present a hydrological analysis, a detailed geomorphic map, interpretations of exposed landslide interiors, and radiocarbon dating of the El Mirador landslide, which is one of the largest and morphologically most representative landslide. We find that the presence of lakes on top of the plateau, causing low infiltration rates, correlates negatively with the abundance of earthflows on compound landslides along the plateau margins. Field outcrops show that the pattern of compound landslides and earthflows is likely controlled by groundwater seepage at the contact between the basalts and underlying soft Miocene molasse. Numerous peat bogs store water and sediment and are more abundant in earthflow-affected areas than in their contributing catchment areas. <br /> Radiocarbon dates indicate that these earthflows displaced metre-thick layers of peat in the late Holocene (<2.5 ka). We conclude that earthflows of the MLBA plateau might be promising proxies of past hydroclimatic conditions in the Patagonian foreland, if strong earthquakes or gradual crustal stress changes due to glacioisostatic rebound can be ruled out.
Drainage divide networks
(2020)
Drainage divides are organized into tree-like networks that may record information about drainage divide mobility. However, views diverge about how to best assess divide mobility. Here, we apply a new approach of automatically extracting and ordering drainage divide networks from digital elevation models to results from landscape evolution model experiments. We compared landscapes perturbed by strike-slip faulting and spatiotemporal variations in erodibility to a reference model to assess which topographic metrics (hillslope relief, flow distance, and chi) are diagnostic of divide mobility. Results show that divide segments that are a minimum distance of similar to 5 km from river confluences strive to attain constant values of hillslope relief and flow distance to the nearest stream. Disruptions of such patterns can be related to mobile divides that are lower than stable divides, closer to streams, and often asymmetric in shape. In general, we observe that drainage divides high up in the network, i.e., at great distances from river confluences, are more susceptible to disruptions than divides closer to these confluences and are thus more likely to record disturbance for a longer time period. We found that across-divide differences in hillslope relief proved more useful for assessing divide migration than other tested metrics. However, even stable drainage divide networks exhibit across-divide differences in any of the studied topographic metrics. Finally, we propose a new metric to quantify the connectivity of divide junctions.
Dispersion-curve inversion of Rayleigh waves to infer subsurface shear-wave velocity is a long-standing problem in seismology. Due to nonlinearity and ill-posedness, sophisticated regularization techniques are required to solve the problem for a stable velocity model. We have formulated the problem as a minimization problem with nonlinear operator constraint and then solve it by using an inexact augmented Lagrangian method, taking advantage of the Haney-Tsai Dix-type relation (a global linear approximation of the nonlinear forward operator). This replaces the original regularized nonlinear problem with iterative minimization of a more tractable regularized linear problem followed by a nonlinear update of the phase velocity (data) in which the update can be performed accurately with any forward modeling engine, for example, the finite-element method. The algorithm allows discretizing the medium with thin layers (for the finite-element method) and thus omitting the layer thicknesses from the unknowns and also allows incorporating arbitrary regularizations to shape the desired velocity model. In this research, we use total variation regularization to retrieve the shear-wave velocity model. We use two synthetic and two real data examples to illustrate the performance of the inversion algorithm with total variation regularization. We find that the method is fast and stable, and it converges to the solution of the original nonlinear problem.
The 2.7-2.9 Ma Ertsberg East Skarn System (EESS) is a world-class Cu-Au skarn that formed within and adjacent to an intrusion within a paleodepth of 0.5 km and > 2.5 km. Its economic mineralisation developed by sustained reaction of magmatic fluid with contact metamorphosed siliciclastic and carbonate rocks at the margin of the adjacent Ertsberg quartz monzodiorite intrusion. Based on high-resolution mineral mapping, chemical analysis and thermodynamic calculations, the multistage formation processes of the exoskarn components of the EESS are examined in the context of changing pressure, temperature, fluid composition and fluid phase. We show that contact metamorphism of dolomitic sediments occurred at 51 +/- 5 MPa, between 700 degrees C and 800 degrees C and in the presence of a H2O-CO2-fluid containing similar to 10 to similar to 70 mol% CO2. This prograde metamorphism formed a forsterite + diopside + calcite + phlogopite + spinel assemblage. Such forsterite-dominated skarns account for similar to 55 vol% of the EESS exoskarns. Rare pargasite (previously unrecognized in this deposit) formed locally in the metamorphosed carbonate sequence where the protolith was composed of supratidal evaporites with dolomitic carbonate and interlayered calc-silicate rocks. The subsequent flux of a lower pressure magmatic gas containing SO2(g) caused sulphate metasomatism. This high temperature gas alteration of the metamorphic assemblage also caused skarn Cu-Fe-sulphide mineralisation. The influx of a SO2 gas through fracture permeability occurred at a temperature between similar to 600 and 700 degrees C and caused calcite to be replaced by anhydrite, with the coupled release of H2S(g). This in-situ release of H2S(g) scavenged trace Cu from the gas phase to deposit Cu-Fe-sulphides, which make the economic value of the distinct. We demonstrate that the formation of metal sulphides within forsterite skarns of the Ertsberg East Skarn System required a minimum flux of similar to 1,050 Mt SO2(g) and show that volcanic degassing may have occurred over a time span of similar to 3,900 years. As the system waned, the ambient fluid resulted in partial retrograde serpentinization of olivine and diopside without carbonation, and at temperatures sufficiently high to preserve anhydrite.
Increasingly available high-frequency data during storm events, when hydrological dynamics most likely activate nitrate storage-flux exchanges, reveal insights into catchment nitrate dynamics. In this study, we explored impacts of seasonality and landscape gradients on nitrate concentration-discharge (C-Q) hysteresis patterns in the Selke catchment, central Germany, which has heterogeneous combinations of meteorological, hydrogeological and land use conditions. Three nested gauging stations established along the main Selke River captured flow and nitrate export dynamics from the uppermost subcatchment (mixed forest and arable land), middle subcatchment (pure steep forest) and lowermost subcatchment (arable and urban land). We collected continuous high-frequency (15-min) discharge and nitrate concentration data from 2012 to 2017 and analyzed the 223 events detected at all three stations. A dominant hysteresis pattern in the uppermost and middle subcatchments was counter-clockwise and combined with an accretion effect, indicating many proximal and mobilized distal nitrate sources. However, 66% of all events at the catchment outlet experienced a dilution effect, possibly due to mechanisms that vary seasonally. During wetting/wet periods (October-March), it was combined mainly with a counter-clockwise pattern due to the dominance of event runoff volume from the uppermost and middle subcatchments. During drying/dry periods (April-September), however, it was combined mainly with a clockwise pattern due to occasional quick surface flows from lowland near-stream urban areas. In addition, the clockwise hysteresis occurred mainly from May-October during mostly drying/dry periods at all three sites, indicating little distal nitrate transport in response to the low terrestrial hydrological connectivity, especially in the lowermost dry and flat sub-catchment. This comprehensive analysis (i.e., clockwise vs. counter-clockwise, accretion vs. dilution) enables in-depth analysis of nitrate export mechanisms during certain periods under different landscape conditions. Specific combination of C-Q relationships could identify target locations for agricultural management actions that decrease nitrate output. Therefore, we strongly encourage long-term multisite and high-frequency monitoring strategies in heterogeneous nested catchment(s), which can help understand process mechanisms, generate data for physical-based water-quality modeling and provide guidance for water and agricultural management.