Refine
Has Fulltext
- no (240) (remove)
Year of publication
- 2019 (240) (remove)
Document Type
- Article (216)
- Other (9)
- Review (7)
- Doctoral Thesis (6)
- Habilitation Thesis (2)
Is part of the Bibliography
- yes (240)
Keywords
- climate (4)
- Holocene (3)
- paleoclimate (3)
- remote sensing (3)
- AMOC (2)
- Asia (2)
- Central Andes (2)
- Climate (2)
- Climate change (2)
- Climate change adaptation (2)
Institute
- Institut für Geowissenschaften (240) (remove)
Ongoing development of dynamical atmosphere-ocean general circulation models keep expectations high regarding seasonal predictions of Indian monsoon rainfall. This study compares past and present skill of four currently operating forecasting systems, CFSv2 from NCEP, ENSEMBLES, System 4 and the newest SEAS5 from ECMWF, by analysing correlations of respective hindcasts with observed all-India summer rainfall. For the common time period 1982-2005, only ENSEMBLES and CFSv2 give significantly skilful forecasts. It is shown that skill is highly dependent on the chosen time period. Especially the intense El Nino of 1997 seems to degrade the predictions, most notably for SEAS4 and SEAS5 which seem to be linked to El Nino too strongly. We show that by discarding that year, a regime shift in the 1990s is no longer visible. Overall, we observe a convergence of skill towards the present with correlations of about 0.4 for CFSv2 and of 0.6 for System 4 and SEAS5.
Niche-based species distribution models (SDMs) have become an essential tool in conservation and restoration planning. Given the current threats to freshwater biodiversity, it is of fundamental importance to address scale effects on the performance of niche-based SDMs of freshwater species’ distributions. The scale effects are addressed here in the context of hierarchical catchment ordering, considered as counterpart to coarsening grain-size by increasing grid-cell size. We combine fish occurrence data from the Danube River Basin, the hierarchical catchment ordering and multiple environmental factors representing topographic, climatic and anthropogenic effects to model fish occurrence probability across multiple scales. We focus on 1st to 5th order catchments. The spatial scale (hierarchical catchment order) only marginally influences the mean performance of SDMs, however the uncertainty of the estimates increases with scale. Key predictors and their relative importance are scale and species dependent. Our findings have useful implications for choosing proper species dependent spatial scales for river rehabilitation measures, and for conservation planning in areas where fine grain species data are unavailable.
Lake sediments are increasingly explored as reliable paleoflood archives. In addition to established flood proxies including detrital layer thickness, chemical composition, and grain size, we explore stable oxygen and carbon isotope data as paleoflood proxies for lakes in catchments with carbonate bedrock geology. In a case study from Lake Mondsee (Austria), we integrate high-resolution sediment trapping at a proximal and a distal location and stable isotope analyses of varved lake sediments to investigate flood-triggered detrital sediment flux. First, we demonstrate a relation between runoff, detrital sediment flux, and isotope values in the sediment trap record covering the period 2011-2013 CE including 22 events with daily (hourly) peak runoff ranging from 10 (24) m(3) s(-1) to 79 (110) m(3) s(-1). The three- to ten-fold lower flood-triggered detrital sediment deposition in the distal trap is well reflected by attenuated peaks in the stable isotope values of trapped sediments. Next, we show that all nine flood-triggered detrital layers deposited in a sediment record from 1988 to 2013 have elevated isotope values compared with endogenic calcite. In addition, even two runoff events that did not cause the deposition of visible detrital layers are distinguished by higher isotope values. Empirical thresholds in the isotope data allow estimation of magnitudes of the majority of floods, although in some cases flood magnitudes are overestimated because local effects can result in too-high isotope values. Hence we present a proof of concept for stable isotopes as reliable tool for reconstructing flood frequency and, although with some limitations, even for flood magnitudes.
We show that near-real-time seismic monitoring of fluid injection allowed control of induced earthquakes during the stimulation of a 6.1-km-deep geothermal well near Helsinki, Finland. A total of 18,160 m(3) of fresh water was pumped into crystalline rocks over 49 days in June to July 2018. Seismic monitoring was performed with a 24-station borehole seismometer network. Using near-real-time information on induced-earthquake rates, locations, magnitudes, and evolution of seismic and hydraulic energy, pumping was either stopped or varied-in the latter case, between well-head pressures of 60 and 90 MPa and flow rates of 400 and 800 liters/min. This procedure avoided the nucleation of a project-stopping magnitude M-W 2.0 induced earthquake, a limit set by local authorities. Our results suggest a possible physics-based approach to controlling stimulation-induced seismicity in geothermal projects.
Human mortality shows a pronounced temperature dependence. The minimum mortality temperature (MMT) as a characteristic point of the temperature-mortality relationship is influenced by many factors. As MMT estimates are based on case studies, they are sporadic, limited to data-rich regions, and their drivers have not yet been clearly identified across case studies. This impedes the elaboration of spatially comprehensive impact studies on heat-related mortality and hampers the temporal transfer required to assess climate change impacts. Using 400 MMTs from cities, we systematically establish a generalised model that is able to estimate MMTs (in daily apparent temperature) for cities, based on a set of climatic, topographic and socio-economic drivers. A sigmoid model prevailed against alternative model setups due to having the lowest Akaike Information Criterion (AICc) and the smallest RMSE. We find the long-term climate, the elevation, and the socio-economy to be relevant drivers of our MMT sample within the non-linear parametric regression model. A first model application estimated MMTs for 599 European cities ( >100 000 inhabitants) and reveals a pronounced decrease in MMTs (27.8-16 degrees C) from southern to northern cities. Disruptions of this pattern across regions of similar mean temperatures can be explained by socio-economic standards as noted for central eastern Europe. Our alternative method allows to approximate MMTs independently from the availability of daily mortality records. For the first time, a quantification of climatic and non-climatic MMT drivers has been achieved, which allows to consider changes in socio-economic conditions and climate. This work contributes to the comparability among MMTs beyond location-specific and regional limits and, hence, towards a spatially comprehensive impact assessment for heat-related mortality.
Hungry cities: how local food self-sufficiency relates to climate change, diets, and urbanisation
(2019)
Using a newly developed model approach and combining it with remote sensing, population, and climate data, first insights are provided into how local diets, urbanisation, and climate change relates to local urban food self-sufficiency. In plain terms, by utilizing the global peri-urban (PU) food production potential approximately lbn urban residents (30% of global urban population) can be locally nourished, whereby further urbanisation is by far the largest pressure factor on PU agriculture, followed by a change of diets, and climate change. A simple global food transport model which optimizes transport and neglects differences in local emission intensities indicates that CO2 emissions related to food transport can be reduced by a factor of 10.
Shear-waves are the most energetic body-waves radiated from an earthquake, and are responsible for the destruction of engineered structures. In both short-term emergency response and long-term risk forecasting of disaster-resilient built environment, it is critical to predict spatially accurate distribution of shear-wave amplitudes. Although decades’ old theory proposes a deterministic, highly anisotropic, four-lobed shear-wave radiation pattern, from lack of convincing evidence, most empirical ground-shaking prediction models settled for an oversimplified stochastic radiation pattern that is isotropic on average. Today, using the large datasets of uniformly processed seismograms from several strike, normal, reverse, and oblique-slip earthquakes across the globe, compiled specifically for engineering applications, we could reveal, quantify, and calibrate the frequency-, distance-, and style-of-faulting dependent transition of shear-wave radiation between a stochastic-isotropic and a deterministic-anisotropic phenomenon. Consequent recalibration of empirical ground-shaking models dramatically improved their predictions: with isodistant anisotropic variations of ±40%, and 8% reduction in uncertainty. The outcomes presented here can potentially trigger a reappraisal of several practical issues in engineering seismology, particularly in seismic ground-shaking studies and seismic hazard and risk assessment.
Increased landslide activity on forested hillslopes following two recent volcanic eruptions in Chile
(2019)
Large explosive eruptions can bury landscapes beneath thick layers of tephra. Rivers subsequently overloaded with excess pyroclastic sediments have some of the highest reported specific sediment yields. Much less is known about how hillslopes respond to tephra loads. Here, we report a pulsed and distinctly delayed increase in landslide activity following the eruptions of the Chaiten (2008) and Puyehue-Cordon Caulle (2011) volcanoes in southern Chile. Remote-sensing data reveal that land-slides clustered in densely forested hillslopes mostly two to six years after being covered by tephra. This lagged instability is consistent with a gradual loss of shear strength of decaying tree roots in areas of high tephra loads. Surrounding areas with comparable topography, forest cover, rainfall and lithology maintained landslide rates roughly ten times lower. The landslides eroded the landscape by up to 4.8 mm on average within 30 km of both volcanoes, mobilizing up to 1.6 MtC at rates of about 265 tC km(-2) yr(-1). We suggest that these yields may reinforce the elevated river loads of sediment and organic carbon in the decade after the eruptions. We recommend that studies of post-eruptive mass fluxes and hazards include lagged landslide responses of tephra-covered forested hillslopes, to avoid substantial underestimates.
With the onset of the global food crisis, the discussion about the use and misuse of agricultural market interventions regained academic attention. As a result of economies of scale, centralized policy implementation at the regional level has the potential to reduce the budgetary costs of policies. Borrowing from the literature on international unions and international policy coordination, we develop a conceptual framework to analyze when regional policy implementation makes sense. This is the case whenever spill-overs from centralization are large and policy preferences, driven by country-specific characteristics, are homogeneous. Subsequently, we examine the advantageousness of centralized policy implementation for the West African region regarding the most common food security policies. We show that centralization of trade policies and emergency food reserves is beneficial, while buffer stocks, safety net policies, and producer support policies should be implemented at the national level.
The Sentinel Application Platform (SNAP) architecture facilitates Earth Observation data processing. In this work, we present results from a new Snow Processor for SNAP. We also describe physical principles behind the developed snow property retrieval technique based on the analysis of Ocean and Land Colour Instrument (OLCI) onboard Sentinel-3A/B measurements over clean and polluted snow fields. Using OLCI spectral reflectance measurements in the range 400-1020 nm, we derived important snow properties such as spectral and broadband albedo, snow specific surface area, snow extent and grain size on a spatial grid of 300 m. The algorithm also incorporated cloud screening and atmospheric correction procedures over snow surfaces. We present validation results using ground measurements from Antarctica, the Greenland ice sheet and the French Alps. We find the spectral albedo retrieved with accuracy of better than 3% on average, making our retrievals sufficient for a variety of applications. Broadband albedo is retrieved with the average accuracy of about 5% over snow. Therefore, the uncertainties of satellite retrievals are close to experimental errors of ground measurements. The retrieved surface grain size shows good agreement with ground observations. Snow specific surface area observations are also consistent with our OLCI retrievals. We present snow albedo and grain size mapping over the inland ice sheet of Greenland for areas including dry snow, melted/melting snow and impurity rich bare ice. The algorithm can be applied to OLCI Sentinel-3 measurements providing an opportunity for creation of long-term snow property records essential for climate monitoring and data assimilation studies-especially in the Arctic region, where we face rapid environmental changes including reduction of snow/ice extent and, therefore, planetary albedo.
During lower sea levels in glacial periods, deep permafrost formed on large continental shelf areas of the Arctic Ocean. Subsequent sea level rise and coastal erosion created subsea permafrost, which generally degrades after inundation under the influence of a complex suite of marine, near-shore processes. Global warming is especially pronounced in the Arctic, and will increase the transition to and the degradation of subsea permafrost, with implications for atmospheric climate forcing, offshore infrastructure, and aquatic ecosystems.
This thesis combines new geophysical, borehole observational and modelling approaches to enhance our understanding of subsea permafrost dynamics. Three specific areas for advancement were identified: (I) sparsity of observational data, (II) lacking implementation of salt infiltration mechanisms in models, and (III) poor understanding of the regional differences in key driving parameters. This study tested the combination of spectral ratios of the ambient vibration seismic wavefield, together with estimated shear wave velocity from seismic interferometry analysis, for estimating the thickness of the unfrozen sediment overlying the ice-bonded permafrost offshore. Mesoscale numerical calculations (10^1 to 10^2 m, thousands of years) were employed to develop and solve the coupled heat diffusion and salt transport equations including phase change effects. Model soil parameters were constrained by borehole data, and the impact of a variety of influences during the transgression was tested in modelling studies. In addition, two inversion schemes (particle swarm optimization and a least-square method) were used to reconstruct temperature histories for the past 200-300 years in the Laptev Sea region in Siberia from two permafrost borehole temperature records. These data were evaluated against larger scale reconstructions from the region.
It was found (I) that peaks in spectral ratios modelled for three-layer, one-dimensional systems corresponded with thaw depths. Around Muostakh Island in the central Laptev Sea seismic receivers were deployed on the seabed. Derived depths of the ice-bonded permafrost table were between 3.7-20.7 m ± 15 %, increasing with distance from the coast. (II) Temperatures modelled during the transition to subsea permafrost resembled isothermal conditions after about 2000 years of inundation at Cape Mamontov Klyk, consistent with observations from offshore boreholes. Stratigraphic scenarios showed that salt distribution and infiltration had a large impact on the ice saturation in the sediments. Three key factors were identified that, when changed, shifted the modelled permafrost thaw depth most strongly: bottom water temperatures, shoreline retreat rate and initial temperature before inundation. Salt transport based on diffusion and contribution from arbitrary density-driven mechanisms only accounted for about 50 % of observed thaw depths at offshore sites hundreds to thousands of years after inundation. This bias was found consistently at all three sites in the Laptev Sea region. (III) In the temperature reconstructions, distinct differences in the local temperature histories between the western Laptev Sea and the Lena Delta sites were recognized, such as a transition to warmer temperatures a century later in the western Laptev Sea as well as a peak in warming three decades later. The local permafrost surface temperature history at Sardakh Island in the Lena Delta was reminiscent of the circum-Arctic regional average trends. However, Mamontov Klyk in the western Laptev Sea was consistent to Arctic trends only in the most recent decade and was more similar to northern hemispheric mean trends. Both sites were consistent with a rapid synoptic recent warming.
In conclusion, the consistency between modelled response, expected permafrost distribution, and observational data suggests that the passive seismic method is promising for the determination of the thickness of unfrozen sediment on the continental Arctic shelf. The quantified gap between currently modelled and observed thaw depths means that the impact of degradation on climate forcing, ecosystems, and infrastructure is larger than current models predict. This discrepancy suggests the importance of further mechanisms of salt penetration and thaw that have not been considered – either pre-inundation or post-inundation, or both. In addition, any meaningful modelling of subsea permafrost would have to constrain the identified key factors and their regional differences well. The shallow permafrost boreholes provide missing well-resolved short-scale temperature information in the coastal permafrost tundra of the Arctic. As local differences from circum-Arctic reconstructions, such as later warming and higher warming magnitude, were shown to exist in this region, these results provide a basis for local surface temperature record parameterization of climate and, in particular, permafrost models. The results of this work bring us one step further to understanding the full picture of the transition from terrestrial to subsea permafrost.
The Arctic is directly impacted by climate change. The increase in air temperature drives the thawing of permafrost and an increase in coastal erosion and river discharge. This leads to a greater input of sediment and organic matter into coastal waters, which substantially impacts the ecosystems, the subsistence economy of the local population, and the climate because of the transformation of organic matter into greenhouse gases. Yet, the patterns of sediment dispersal in the nearshore zone are not well known, because ships do not often reach shallow waters and satellite remote sensing is traditionally focused on less dynamic environments. The goal of this study is to use the extensive Landsat archive to investigate sediment dispersal patterns specifically on an exemplary Arctic nearshore environment, where field measurements are often scarce. Multiple Landsat scenes were combined to calculate means of sediment dispersal and sea surface temperature under changing seasonal wind conditions in the nearshore zone of Herschel Island Qikiqtaruk in the western Canadian Arctic since 1982. We use observations in the Landsat red and thermal wavebands, as well as a recently published water turbidity algorithm to relate archive wind data to turbidity and sea surface temperature. We map the spatial patterns of turbidity and water temperature at high spatial resolution in order to resolve transport pathways of water and sediment at the water surface. Our results show that these pathways are clearly related to the prevailing wind conditions, being ESE and NW. During easterly wind conditions, both turbidity and water temperature are significantly higher in the nearshore area. The extent of the Mackenzie River plume and coastal erosion are the main explanatory variables for sediment dispersal and sea surface temperature distributions in the study area. During northwesterly wind conditions, the influence of the Mackenzie River plume is negligible. Our results highlight the potential of high spatial resolution Landsat imagery to detect small-scale hydrodynamic processes, but also show the need to specifically tune optical models for Arctic nearshore environments.
The progress of science is tied to the standardization of measurements, instruments, and data. This is especially true in the Big Data age, where analyzing large data volumes critically hinges on the data being standardized. Accordingly, the lack of community-sanctioned data standards in paleoclimatology has largely precluded the benefits of Big Data advances in the field. Building upon recent efforts to standardize the format and terminology of paleoclimate data, this article describes the Paleoclimate Community reporTing Standard (PaCTS), a crowdsourced reporting standard for such data. PaCTS captures which information should be included when reporting paleoclimate data, with the goal of maximizing the reuse value of paleoclimate data sets, particularly for synthesis work and comparison to climate model simulations. Initiated by the LinkedEarth project, the process to elicit a reporting standard involved an international workshop in 2016, various forms of digital community engagement over the next few years, and grassroots working groups. Participants in this process identified important properties across paleoclimate archives, in addition to the reporting of uncertainties and chronologies; they also identified archive-specific properties and distinguished reporting standards for new versus legacy data sets. This work shows that at least 135 respondents overwhelmingly support a drastic increase in the amount of metadata accompanying paleoclimate data sets. Since such goals are at odds with present practices, we discuss a transparent path toward implementing or revising these recommendations in the near future, using both bottom-up and top-down approaches.
Paleogene evolution and demise of the proto-Paratethys Sea in Central Asia (Tarim and Tajik basins)
(2019)
The proto-Paratethys Sea covered a vast area extending from the Mediterranean Tethys to the Tarim Basin in western China during Cretaceous and early Paleogene. Climate modelling and proxy studies suggest that Asian aridification has been governed by westerly moisture modulated by fluctuations of the proto-Paratethys Sea. Transgressive and regressive episodes of the proto-Paratethys Sea have been previously recognized but their timing, extent and depositional environments remain poorly constrained. This hampers understanding of their driving mechanisms (tectonic and/or eustatic) and their contribution to Asian aridification. Here, we present a new chronostratigraphic framework based on biostratigraphy and magnetostratigraphy as well as a detailed palaeoenvironmental analysis for the Paleogene proto-Paratethys Sea incursions in the Tajik and Tarim basins. This enables us to identify the major drivers of marine fluctuations and their potential consequences on Asian aridification. A major regional restriction event, marked by the exceptionally thick (<= 400 m) shelf evaporites is assigned a Danian-Selandian age (ca. 63-59 Ma) in the Aertashi Formation. This is followed by the largest recorded proto-Paratethys Sea incursion with a transgression estimated as early Thanetian (ca. 59-57 Ma) and a regression within the Ypresian (ca. 53-52 Ma), both within the Qimugen Formation. The transgression of the next incursion in the Kalatar and Wulagen formations is now constrained as early Lutetian (ca. 47-46 Ma), whereas its regression in the Bashibulake Formation is constrained as late Lutetian (ca. 41 Ma) and is associated with a drastic increase in both tectonic subsidence and basin infilling. The age of the final and least pronounced sea incursion restricted to the westernmost margin of the Tarim Basin is assigned as Bartonian-Priabonian (ca. 39.7-36.7 Ma). We interpret the long-term westward retreat of the proto-Paratethys Sea starting at ca. 41 Ma to be associated with far-field tectonic effects of the Indo-Asia collision and Pamir/Tibetan plateau uplift. Short-term eustatic sea level transgressions are superimposed on this long-term regression and seem coeval with the transgression events in the other northern Peri-Tethyan sedimentary provinces for the 1st and 2nd sea incursions. However, the 3rd sea incursion is interpreted as related to tectonism. The transgressive and regressive intervals of the proto-Paratethys Sea correlate well with the reported humid and arid phases, respectively in the Qaidam and Xining basins, thus demonstrating the role of the proto-Paratethys Sea as an important moisture source for the Asian interior and its regression as a contributor to Asian aridification.
Small-scale variations in mineral chemistry, textures, and platinum group element (PGE) mineralization were investigated in the Lower and Middle Group chromitite layers LG6, LG6a, MG1, MG2, and MG2 II from vertical drill core profiles at the Thaba mine in the northwestern limb of the Bushveld Complex. We present detailed geochemical profiles of chromite composition and chromite crystal size distribution curves to shed light on the processes of chromite accumulation and textural modification as well as mineralization. Multiple samples within each layer were assayed for PGE concentrations, and the respective platinum group mineral association was determined by mineral liberation analysis (MLA). There is strong evidence for postcumulus changes in the chromitites. The crystal size distribution curves suggest that the primary chromite texture was coarsened by a combination of adcumulus growth and textural equilibration, while compaction of the crystal mush played only a minor role. Mineral compositions were also modified by postcumulus processes, but because of the very high modal amount of chromite and its local preservation in orthopyroxene oikocrysts, that phase retained much primary information. Vertical variations of chromite composition within chromitite layers and from one layer to another do not support the idea of chromite accumulation from crystal-rich slurries or crystal settling from a large magma chamber. Instead, we favor a successive buildup of chromitite layers by repeated injections of relatively thin layers of chromite-saturated magmas, with in situ crystallization occurring at the crystal mush-magma interface. The adcumulus growth of chromite grains to form massive chromitite required addition of Cr to the layers, which we attribute to downward percolation from the overlying magma. The PGE concentrations are elevated in all chromitite layers compared to adjacent silicate rocks and show a systematic increase upward from LG6 (avg 807 ppb Ir + Ru + Rh + Pt + Pd + Au) to MG2 II (avg 2,062 ppb). There are also significant internal variations in all layers, with enrichments at hanging and/or footwalls. The enriched nature of chromitites in PGEs compared to host pyroxenites is a general feature, independent of the layer thickness. The MLA results distinguish two principal groups of PGE mineral associations: the LG6, LG6, and MG1 are dominated by the malanite series, laurite, and PGE sulfarsenides, while the MG2 and MG2 II layers are characterized by laurite and PGE sulfides as well as Pt-Fe-Sn and PGE-Sb-Bi-Pb alloys. Differences in the PGE associations are attributed to postcumulus alteration of the MG2 and MG2 II layer, while the chromitites below, particularly LG6 and LG6a, contain a more pristine association.
Small-aperture array as a tool to monitor fluid injection- and extraction-induced microseismicity
(2019)
The monitoring of microseismicity during temporary human activities such as fluid injections for hydrofracturing, hydrothermal stimulations or wastewater disposal is a difficult task. The seismic stations often cannot be installed on hard rock, and at quiet places, noise is strongly increased during the operation itself and the installation of sensors in deep wells is costly and often not feasible. The combination of small-aperture seismic arrays with shallow borehole sensors offers a solution. We tested this monitoring approach at two different sites: (1) accompanying a fracking experiment in sedimentary shale at 4km depth and (2) above a gas field under depletion. The small-aperture arrays were planned according to theoretical wavenumber studies combined with simulations considering the local noise conditions. We compared array recordings with recordings available from shallow borehole sensors and give examples of detection and location performance. Although the high-frequency noise on the 50-m-deep borehole sensors was smaller compared to the surface noise before the injection experiment, the signals were highly contaminated during injection by the pumping activities. Therefore, a set of three small-aperture arrays at different azimuths was more suited to detect small events, since noise recorded on these arrays is uncorrelated with each other. Further, we developed recommendations for the adaptation of the monitoring concept to other sites experiencing induced seismicity.
The significance of phytoliths for the control of silicon (Si) fluxes from terrestrial to aquatic ecosystems has been recognized as a key factor. Humankind actively influences Si fluxes by intensified land use, i.e., agriculture and forestry, on a global scale. We hypothesized phytolith distribution and assemblages in soils of agricultural and forestry sites to be controlled by vegetation (which is directed by land use) with direct effects on extractable Si fractions driven mainly by phytolith characteristics, i.e., dissolution status (dissolution signs) and morphology (morphotype proportions). To test our hypothesis we combined different chemical extraction methods (calcium chloride, ammonium oxalate, Tiron) for the quantification of different Si fractions (plant available Si, Si adsorbed to/occluded in pedogenic oxides/hydroxides, amorphous Si) and microscopic techniques (light microscopy, confocal laser scanning microscopy, scanning electron microscopy) for detailed analyses of phytoliths extracted using gravimetric separation (physical extraction) from exemplary loess soils of agricultural (arable land and grassland/meadow) and forestry (beech and pine) sites in Poland. We found differences in dissolution signs, morphotype proportions, and vertical distribution of phytoliths in soil horizons per site. In general, dominant morphotypes of assignable phytoliths in the studied soil profiles were elongate phytoliths and short cells, both of which are typical for grass-dominated vegetation. However, the organic layers of forest soils were dominated by globular phytoliths, which are typical indicators for mosses. As expected soil horizons under different vegetation generally were characterized by differences in extractable Si fractions, especially in the upper soil horizons. However, phytogenic Si pools counter-intuitively showed no correlations with chemically extracted Si fractions and soil pH at all. Our findings indicate that it is necessary to combine microscopic analyses and Si extraction techniques for examinations of Si cycling in biogeosystems, because extractions of Si fractions alone do not allow drawing any conclusions about phytolith characteristics or interactions between phytolith pools and chemically extractable Si fractions and do not necessarily reflect phytogenic Si pool quantities in soils and vice versa.
Lineament mapping, which is an important part of any structural geological investigation, is made more efficient and easier by the availability of optical as well as radar remote sensing data, such as Landsat and Sentinel with medium and high spatial resolutions. However, the results from these multi-resolution data vary due to their difference in spatial resolution and sensitivity to soil occupation. The accuracy and quality of extracted lineaments depend strongly on the spatial resolution of the imagery. Therefore, the aim of this study was to compare the optical Landsat-8, Sentinel-2A, and radar Sentinel-1A satellite data for automatic lineament extraction. The framework of automatic approach includes defining the optimal parameters for automatic lineament extraction with a combination of edge detection and line-linking algorithms and determining suitable bands from optical data suited for lineament mapping in the study area. For the result validation, the extracted lineaments are compared against the manually obtained lineaments through the application of directional filtering and edge enhancement as well as to the lineaments digitized from the existing geological maps of the study area. In addition, a digital elevation model (DEM) has been utilized for an accuracy assessment followed by the field verification. The obtained results show that the best correlation between automatically extracted lineaments, manual interpretation, and the preexisting lineament map is achieved from the radar Sentinel-1A images. The tests indicate that the radar data used in this study, with 5872 and 5865 lineaments extracted from VH and VV polarizations respectively, is more efficient for structural lineament mapping than the Landsat-8 and Sentinel-2A optical imagery, from which 2338 and 4745 lineaments were extracted respectively.
The morphology of marine and lacustrine terraces has been largely used to measure past sea- and lake-level positions and estimate vertical deformation in a wealth of studies focused on climate and tectonic processes. To obtain accurate morphometric assessments of terrace morphology we present TerraceM-2, an improved version of our MatlabR (R) graphic-user interface that provides new methodologies for morphometric analyses as well as landscape evolution and fault-dislocation modeling. The new version includes novel routines to map the elevation and spatial distribution of terraces, to model their formation and evolution, and to estimate fault-slip rates from terrace deformation patterns. TerraceM-2 has significantly improves its processing speed and mapping capabilities, and includes separate functions for developing customized workflows beyond the graphic-user interface. We illustrate these new mapping and modeling capabilities with three examples: mapping lacustrine shorelines in the Dead Sea to estimate deformation across the Dead Sea Fault, landscape evolution modeling to estimate a history of uplift rates in southern Peru, and dislocation modeling of deformed marine terraces in California. These examples also illustrate the need to use topographic data of different resolutions. The new modeling and mapping routines of TerraceM-2 highlight the advantages of an integrated joint mapping and modeling approach to improve the efficiency and precision of coastal terrace metrics in both marine and lacustrine environments.
The paper presents a historical long-distance communication system based on beacon fires in one of the most dynamic and rugged mountain ranges of the world, the Hindu Kush-Karakoram-Himalayas. It was deployed as an early warning system for glacial lake outburst floods, which caused devastating impacts on settlement zones and infrastructure until the middle of the twentieth century. The study revealed that the beacon fire systems were operated in distinct valleys spread over the entire Hindu Kush-Karakoram Region. The remarkable fact is the establishment of fire posts in highly difficult accessible mountain environments with communication distances of several hundred kilometres for individual beacon lines. The warning system was a cooperative natural hazard management, which was operated even across distinct ethnic groups. Distant societies, formerly perceived as isolated villages by physical barriers of the high mountain relief, were in historical times connected not only by challenging trade routes but also by a fast working optical communication system. The findings are discussed in the context of a future sustainable natural hazard management.
We present a new three-dimensional density model of the Central Andes characterizing the structure and composition of the lithosphere together with a geodynamic simulation subjected to continental intraplate shortening. The principal aim of this study is to assess the link between heterogeneities in the lithosphere and different deformation patterns and styles along the orogen-foreland system of the Central Andes. First, we performed a 3D integration of new geological and geophysical data with previous models through forward modelling of Bouguer anomalies. Subsequently, a geodynamic model was set-up and parametrized from the previously obtained 3D structure and composition. We do not find a unambigous correlation between the resulting density configuration and terrane boundaries proposed by other authors. Our models reproduce the observed Bouguer anomaly and deformation patterns in the foreland. We find that thin-skinned deformation in the Subandean fold-and thrust belt is controlled by a thick sedimentary layer and coeval underthrusting of thin crust of the foreland beneath the thick crust of the Andean Plateau. In the adjacent thick-skinned deformation province of the inverted Cretaceous extensional Santa Barbara System sedimentary strata are much thinner and crustal thickness transitions from greater values in the Andean to a more reduced thickness in the foreland. Our results show that deformation processes occur where the highest gradients of lithospheric strength are present between the orogen and the foreland, thus suggesting a spatial correlation between deformation and lithospheric strength.
Foreword
(2019)
The intangible impacts of floods on welfare are not well investigated, even though they are important aspects of welfare. Moreover, flooding has gender based impacts on welfare. These differing impacts create a gender based flood risk resilience gap. We study the intangible impacts of flood risk on the subjective well-being of residents in central Vietnam. The measurement of intangible impacts through subjective well-being is a growing field within flood risk research. We find an initial drop in welfare through subjective well-being across genders when a flood is experienced. Male respondents tended to recover their welfare losses by around 80% within 5 years while female respondents were associated with a welfare recovery of around 70%. A monetization of the impacts floods have on an individual’s subjective well-being shows that for the average female respondent, between 41% to 86% of annual income would be required to compensate subjective well-being losses after 5 years of experiencing a flood. The corresponding value for males is 30% to 57% of annual income. This shows that the intangible impacts of flood risk are important (across genders) and need to be integrated into flood (or climate) risk assessments to develop more socially appropriate risk management strategies.
Flood disasters severely impact human subjective well-being (SWB). Nevertheless, few studies have examined the influence of flood events on individual well-being and how such impacts may be limited by flood protection measures. This study estimates the long term impacts on individual subjective well-being of flood experiences, individual subjective flood risk perceptions, and household flood preparedness decisions. These effects are monetised and placed in context through a comparison with impacts of other adverse events on well-being. We collected data from households in flood-prone areas in France. The results indicate that experiencing a flood has a large negative impact on subjective well-being that is incompletely attenuated over time. Moreover, individuals do not need to be directly affected by floods to suffer SWB losses since subjective well-being is lower for those who expect their flood risk to increase or who have seen a neighbour being flooded. Floodplain inhabitants who prepared for flooding by elevating their home have a higher subjective well-being. A monetisation of the aforementioned well-being impacts shows that a flood requires Euro150,000 in immediate compensation to attenuate SWB losses. The decomposition of the monetised impacts of flood experience into tangible losses and intangible effects on SWB shows that intangible effects are about twice as large as the tangible direct monetary flood losses. Investments in flood protection infrastructure may be under funded if the intangible SWB benefits of flood protection are not taken into account.
Flood risk will increase in many areas around the world due to climate change and increase in economic exposure. This implies that adequate flood insurance schemes are needed to adapt to increasing flood risk and to minimise welfare losses for households in flood-prone areas. Flood insurance markets may need reform to offer sufficient and affordable financial protection and incentives for risk reduction. Here, we present the results of a study that aims to evaluate the ability of flood insurance arrangements in Europe to cope with trends in flood risk, using criteria that encompass common elements of the policy debate on flood insurance reform. We show that the average risk-based flood insurance premium could double between 2015 and 2055 in the absence of more risk reduction by households exposed to flooding. We show that part of the expected future increase in flood risk could be limited by flood insurance mechanisms that better incentivise risk reduction by policyholders, which lowers vulnerability. The affordability of flood insurance can be improved by introducing the key features of public-private partnerships (PPPs), which include public reinsurance, limited premium cross-subsidisation between low- and high-risk households, and incentives for policyholder-level risk reduction. These findings were evaluated in a comprehensive sensitivity analysis and support ongoing reforms in Europe and abroad that move towards risk-based premiums and link insurance with risk reduction, strengthen purchase requirements, and engage in multi-stakeholder partnerships.
One commonly proposed method to limit flood risk is land-use or zoning policies which regulates construction in high-risk areas, in order to reduce economic exposure and its vulnerability to flood events. Although such zoning regulations can be effective in limiting trends in flood risk, they also have adverse impacts on society, for instance by limiting local development of areas near the water. In order to judge whether proposed land-use or zoning policies are a net benefit to society, they should be accepted or rejected based on a societal cost-benefit analysis (CBA). However, conducting a CBA of zoning regulation is complex and comprehensive guidelines of how to do such an analysis are lacking. We offer guidelines for good practice. In order to assess the costs and benefits of zoning as a climate change adaption strategy, they should be assessed at a societal level in order to account for public good features of flood risk reduction strategies, and because costs in one area can be benefits in another region. We propose a multistep process: first, determine the spatial extent of the zoning policy and how interconnected the zoned area is to other locations; second, conduct a CBA using monetary costs and benefits estimated from an integrated hydro-economic model to investigate if total benefits exceed total costs; third, conduct a sensitivity analysis regarding the main assumptions; fourth, conduct a multicriteria analysis (MCA) of the normative outcomes of a zoning policy. A desirable policy is preferred in both the CBA and MCA. This article is categorized under: Engineering Water > Planning Water Human Water > Value of Water Science of Water > Water Extremes Human Water > Methods
Although the Paris Agreement arguably made some progress, interest in supplementary approaches to climate change co-operation persist. This article examines the conditions under which a climate club might emerge and grow. Using agent-based simulations, it shows that even with less than a handful of major actors as initial members, a club can eventually reduce global emissions effectively. To succeed, a club must be initiated by the ‘right’ constellation of enthusiastic actors, offer sufficiently large incentives for reluctant countries and be reasonably unconstrained by conflicts between members over issues beyond climate change. A climate club is particularly likely to persist and grow if initiated by the United States and the European Union. The combination of club-good benefits and conditional commitments can produce broad participation under many conditions.
Large-magnitude fluid-injection induced seismic events are a potential risk for geothermal energy developments worldwide. One potential risk mitigation measure is the application of cyclic injection schemes. After validation at small (laboratory) and meso (mine) scale, the concept has now been applied for the first time at field scale at the Pohang Enhanced Geothermal System (EGS) site in Korea. From 7 August until 14 August 2017 a total of 1756 m(3) of surface water was injected into Pohang well PX-1 at flow rates between 1 and 10 l s(-1), with a maximum wellhead pressure (WHP) of 22.8 MPa, according to a site-specific cyclic soft stimulation schedule and traffic light system. A total of 52 induced microearthquakes were detected in real-time during and shortly after the injection, the largest of M-w 1.9. After that event a total of 1771 m(3) of water was produced back from the well over roughly 1 month, during which time no larger-magnitude seismic event was observed. The hydraulic data set exhibits pressure-dependent injectivity increase with fracture opening between 15 and 17 MPa WHP, but no significant permanent transmissivity increase was observed. The maximum magnitude of the induced seismicity during the stimulation period was below the target threshold of M-w 2.0 and additional knowledge about the stimulated reservoir was gained. Additionally, the technical feasibility of cyclic injection at field scale was evaluated. The major factors that limited the maximum earthquake magnitude are believed to be: limiting the injected net fluid volume, flowback after the occurrence of the largest induced seismic event, using a cyclic injection scheme, the application of a traffic light system, and including a priori information from previous investigations and operations in the treatment design.
In this study, we detect high percentile rainfall events in the eastern central Andes, based on Tropical Rainfall Measuring Mission (TRMM) with a spatial resolution of 0.25 × 0.25°, a temporal resolution of 3 h, and for the duration from 2001 to 2018. We identify three areas with high mean accumulated rainfall and analyze their atmospheric behaviour and rainfall characteristics with specific focus on extreme events. Extreme events are defined by events above the 95th percentile of their daily mean accumulated rainfall. Austral summer (DJF) is the period of the year presenting the most frequent extreme events over these three regions. Daily statistics show that the spatial maxima, as well as their associated extreme events, are produced during the night. For the considered period, ERA-Interim reanalysis data, provided by the European Centre for Medium-Range Weather Forecasts (ECMWF) with 0.75° x0.75° spatial and 6-hourly temporal resolutions, were used for the analysis of the meso- and synoptic-scale atmospheric patterns. Night- and day-time differences indicate a nocturnal overload of northerly and northeasterly low-level humidity flows arriving from tropical South America. Under these conditions, cooling descending air from the mountains may find unstable air at the surface, giving place to the development of strong local convection. Another possible mechanism is presented here: a forced ascent of the low-level flow due to the mountains, disrupting the atmospheric stratification and generating vertical displacement of air trajectories. A Principal Component Analysis (PCA) in T-mode is applied to day- and night-time data during the maximum and extreme events. The results show strong correlation areas over each subregion under study during night-time, whereas during day-time no defined patterns are found. This confirms the observed nocturnal behavior of rainfall within these three hotspots.
Permafrost thaw subjects previously frozen soil organic carbon (SOC) to microbial degradation to the greenhouse gases carbon dioxide (CO2) and methane (CH4). Emission of these gases constitutes a positive feedback to climate warming. Among numerous uncertainties in estimating the strength of this permafrost carbon feedback (PCF), two are: (i) how mineralization of permafrost SOC thawed in saturated anaerobic conditions responds to changes in temperature and (ii) how microbial communities and temperature sensitivities change over time since thaw. To address these uncertainties, we utilized a thermokarst-lake sediment core as a natural chronosequence where SOC thawed and incubated in situ under saturated anaerobic conditions for up to 400 years following permafrost thaw. Initial microbial communities were characterized, and sediments were anaerobically incubated in the lab at four temperatures (0 °C, 3 °C, 10 °C, and 25 °C) bracketing those observed in the lake's talik. Net CH4 production in freshly-thawed sediments near the downward-expanding thaw boundary at the base of the talik were most sensitive to warming at the lower incubation temperatures (0 °C to 3 °C), while the overlying sediments which had been thawed for centuries had initial low abundant methanogenic communities (< 0.02%) and did not experience statistically significant increases in net CH4 production potentials until higher incubation temperatures (10 °C to 25 °C). We propose these observed differences in temperature sensitivities are due to differences in SOM quality and functional microbial community composition that evolve over time; however further research is necessary to better constrain the roles of these factors in determining temperature controls on anaerobic C mineralization.
Issue Despite their rather similar climatic conditions, eastern Eurasia and northern North America are largely covered by different plant functional types (deciduous or evergreen boreal forest) composed of larch or pine, spruce and fir, respectively. I propose that these deciduous and evergreen boreal forests represent alternative quasi-stable states, triggered by their different northern tree refugia that reflect the different environmental conditions experienced during the Last Glacial. Evidence This view is supported by palaeoecological and environmental evidence. Once established, Asian larch forests are likely to have stabilized through a complex vegetation-fire-permafrost soil-climate feedback system. Conclusion With respect to future forest developments, this implies that Asian larch forests are likely to be governed by long-term trajectories and are therefore largely resistant to natural climate variability on time-scales shorter than millennia. The effects of regional human impact and anthropogenic global warming might, however, cause certain stability thresholds to be crossed, meaning that irreversible transitions occur and resulting in marked consequences for ecosystem services on these human-relevant time-scales.
The occurrence of deep low-frequency (DLF) microearthquakes beneath volcanoes is commonly attributed to mass transport in the volcanic plumbing system and used to infer feeding channels from and into magma reservoirs. The key question is how magmas migrate from depth to the shallow crust and whether magma reservoirs are currently being recharged. For the first time since the improvement of the local seismic networks in the East Eifel region (Rhineland-Palatinate, Germany), we detect and locate recurrent DLF earthquakes in the lower crust and upper mantle beneath the Laacher See Volcano (LSV), using a joint data set of permanent sensors and a temporary deployment. So far, eight DLF earthquake sequences were observed in four distinct clusters between 10 and 40 km depth. These clusters of weak events (M-L< 2) align along an approximately 80. southeast dipping line south of the LSV. Moment tensor solutions of these events have large shear components, and the irregular dispersion and long coda of body waves indicate interaction processes between shear cracks and fluids. We find a rotation of P-axes orientation for shallow tectonic earthquakes compared to DLF events, indicating that the stress field in the depth interval of DLF events might favour a vertical migration of magma or magmatic fluids. The caldera of the LSV was formed by the last major eruption of the East Eifel Volcanic Field only 12.9 kyr ago, fed by a shallow magma chamber at 5-8 km depth and erupting a total magma volume of 6.7 km(3). The observed DLF earthquake activity and continuous volcanic gas emissions around the LSV indicate an active magmatic system, possibly connected with an upper mantle melt zone.
The computation of such synthetic GFs is computationally and operationally demanding. As a consequence, the onthe-fly recalculation of synthetic GFs in each iteration of an optimisation is time-consuming and impractical. Therefore, the pre-calculation and efficient storage of synthetic GFs on a dense grid of source to receiver combinations enables the efficient lookup and utilisation of GFs in time-critical scenarios. We present a Python-based framework and toolkit - Pyrocko-GF - that enables the pre-calculation of synthetic GF stores, which are independent of their numerical calculation method and GF transfer function. The framework aids in the creation of such GF stores by interfacing a suite of established numerical forward modelling codes in seismology (computational back ends). So far, interfaces to back ends for layered Earth model cases have been provided; however, the architecture of Pyrocko-GF is designed to cover back ends for other geometries (e.g. full 3-D heterogeneous media) and other physical quantities (e.g. gravity, pressure, tilt). Therefore, Pyrocko-GF defines an extensible GF storage format suitable for a wide range of GF types, especially handling elasticity and wave propagation problems. The framework assists with visualisations, quality control, and the exchange of GF stores, which is supported through an online platform that provides many pre-calculated GF stores for local, regional, and global studies. The Pyrocko-GF toolkit comes with a well-documented application programming interface (API) for the Python programming language to efficiently facilitate forward modelling of geophysical processes, e.g. synthetic waveforms or static displacements for a wide range of source models.
Water stable isotope signatures can provide valuable insights into the catchment internal runoff processes. However, the ability of the water isotope data to constrain the internal apportionments of runoff components in hydrological models for glacierized basins is not well understood. This study developed an approach to simultaneously model the water stable isotopic compositions and runoff processes in a glacierized basin in Central Asia. The fractionation and mixing processes of water stable isotopes in and from the various water sources were integrated into a glacio-hydrological model. The model parameters were calibrated on discharge, snow cover and glacier mass balance data, and additionally isotopic composition of streamflow. We investigated the value of water isotopic compositions for the calibration of model parameters, in comparison to calibration methods without using such measurements. Results indicate that: (1) The proposed isotope-hydrological integrated modeling approach was able to reproduce the isotopic composition of streamflow, and improved the model performance in the evaluation period; (2) Involving water isotopic composition for model calibration reduced the model parameter uncertainty, and helped to reduce the uncertainty in the quantification of runoff components; (3) The isotope-hydrological integrated modeling approach quantified the contributions of runoff components comparably to a three-component tracer-based end-member mixing analysis method for summer peak flows, and required less water tracer data. Our findings demonstrate the value of water isotopic compositions to improve the quantification of runoff components using hydrological models in glacierized basins.
Arboreal epiphytes (plants residing in forest canopies) are present across all major climate zones and play important roles in forest biogeochemistry. The substantial water storage capacity per unit area of the epiphyte "bucket" is a key attribute underlying their capability to influence forest hydrological processes and their related mass and energy flows. It is commonly assumed that the epiphyte bucket remains saturated, or near-saturated, most of the time; thus, epiphytes (particularly vascular epiphytes) can store little precipitation, limiting their impact on the forest canopy water budget. We present evidence that contradicts this common assumption from (i) an examination of past research; (ii) new datasets on vascular epiphyte and epi-soil water relations at a tropical montane cloud forest (Monteverde, Costa Rica); and (iii) a global evaluation of non-vascular epiphyte saturation state using a process-based vegetation model, LiBry. All analyses found that the external and internal water storage capacity of epiphyte communities is highly dynamic and frequently available to intercept precipitation. Globally, non-vascular epiphytes spend <20% of their time near saturation and regionally, including the humid tropics, model results found that non-vascular epiphytes spend similar to 1/3 of their time in the dry state (0-10% of water storage capacity). Even data from Costa Rican cloud forest sites found the epiphyte community was saturated only 1/3 of the time and that internal leaf water storage was temporally dynamic enough to aid in precipitation interception. Analysis of the epi-soils associated with epiphytes further revealed the extent to which the epiphyte bucket emptied-as even the canopy soils were often <50% saturated (29-53% of all days observed). Results clearly show that the epiphyte bucket is more dynamic than currently assumed, meriting further research on epiphyte roles in precipitation interception, redistribution to the surface and chemical composition of "net" precipitation waters reaching the surface.
Nature-based solutions (NBS) have recently received attention due to their potential ability to sustainably reduce hydro-meteorological risks, providing co-benefits for both ecosystems and affected people. Therefore, pioneering research has dedicated efforts to optimize the design of NBS, to evaluate their wider co-benefits and to understand promoting and/or hampering governance conditions for the uptake of NBS. In this article, we aim to complement this research by conducting a comprehensive literature review of factors shaping people’s perceptions of NBS as a means to reduce hydro-meteorological risks. Based on 102 studies, we identified six topics shaping the current discussion in this field of research: (1) valuation of the co-benefits (including those related to ecosystems and society); (2) evaluation of risk reduction efficacy; (3) stakeholder participation; (4) socio-economic and location-specific conditions; (5) environmental attitude, and (6) uncertainty. Our analysis reveals that concerned empirical insights are diverse and even contradictory, they vary in the depth of the insights generated and are often not comparable for a lack of a sound theoretical-methodological grounding. We, therefore, propose a conceptual model outlining avenues for future research by indicating potential inter-linkages between constructs underlying perceptions of NBS to hydro-meteorological risks.
The in-phase response collected by portable loop-loop electromagnetic induction (EMI) sensors operating at low and moderate induction numbers (<= 1) is typically used for sensing the magnetic permeability (or susceptibility) of the subsurface. This is due to the fact that the in-phase response contains a small induction fraction and a preponderant induced magnetization fraction. The magnetization fraction follows the magneto-static equations similarly to the magnetic method but with an active magnetic source. The use of an active source offers the possibility to collect data with several loop-loop configurations, which illuminate the subsurface with different sensitivity patterns. Such multiconfiguration soundings thereby allows the imaging of subsurface magnetic permeability/susceptibility variations through an inversion procedure. This method is not affected by the remnant magnetization and theoretically overcomes the classical depth ambiguity generally encountered with passive geomagnetic data. To invert multiconfiguration in-phase data sets, we propose a novel methodology based on a full-grid 3-D multichannel deconvolution (MCD) procedure. This method allows us to invert large data sets (e.g. consisting of more than a hundred thousand of data points) for a dense voxel-based 3-D model of magnetic susceptibility subject to smoothness constraints. In this study, we first present and discuss synthetic examples of our imaging procedure, which aim at simulating realistic conditions. Finally, we demonstrate the applicability of our method to field data collected across an archaeological site in Auvergne (France) to image the foundations of a Gallo-Roman villa built with basalt rock material. Our synthetic and field data examples demonstrate the potential of the proposed inversion procedure offering new and complementary ways to interpret data sets collected with modern EMI instruments.
Climate science is highly interdisciplinary by nature, so understanding interactions between Earth processes inherently warrants the use of analytical software that can operate across the disciplines of Earth science. Toward this end, we present the Climate Data Toolbox for MATLAB, which contains more than 100 functions that span the major climate-related disciplines of Earth science. The toolbox enables streamlined, entirely scriptable workflows that are intuitive to write and easy to share. Included are functions to evaluate uncertainty, perform matrix operations, calculate climate indices, and generate common data displays. Documentation is presented pedagogically, with thorough explanations of how each function works and tutorials showing how the toolbox can be used to replicate results of published studies. As a well-tested, well-documented platform for interdisciplinary collaborations, the Climate Data Toolbox for MATLAB aims to reduce time spent writing low-level code, let researchers focus on physics rather than coding and encourage more efficacious code sharing. Plain Language Summary This article describes a collection of computer code that has recently been released to help scientists analyze many types of Earth science data. The code in this toolbox makes it easy to investigate things like global warming, El Nino, or other major climate-related processes such as how winds affect ocean circulation. Although the toolbox was designed to be used by expert climate scientists, its instruction manual is well written, and beginners may be able to learn a great deal about coding and Earth science, simply by following along with the provided examples. The toolbox is intended to help scientists save time, help them ensure their analysis is accurate, and make it easy for other scientists to repeat the results of previous studies.
Abstract. The Sea of Marmara, in northwestern Turkey, is a transition zone where the dextral North Anatolian Fault zone (NAFZ) propagates westward from the Anatolian Plate to the Aegean Sea Plate. The area is of interest in the context of seismic hazard of Istanbul, a metropolitan area with about 15 million inhabitants. Geophysical observations indicate that the crust is heterogeneous beneath the Marmara basin, but a detailed characterization of the crustal heterogeneities is still missing. To assess if and how crustal heterogeneities are related to the NAFZ segmentation below the Sea of Marmara, we develop new crustal-scale 3-D density models which integrate geological and seismological data and that are additionally constrained by 3-D gravity modeling. For the latter, we use two different gravity datasets including global satellite data and local marine gravity observation. Considering the two different datasets and the general non-uniqueness in potential field modeling, we suggest three possible “end-member” solutions that are all consistent with the observed gravity field and illustrate the spectrum of possible solutions. These models indicate that the observed gravitational anomalies originate from significant density heterogeneities within the crust. Two layers of sediments, one syn-kinematic and one pre-kinematic with respect to the Sea of Marmara formation are underlain by a heterogeneous crystalline crust. A felsic upper crystalline crust (average density of 2720 kgm⁻³) and an intermediate to mafic lower crystalline crust (average density of 2890 kgm⁻³) appear to be cross-cut by two large, dome-shaped mafic highdensity bodies (density of 2890 to 3150 kgm⁻³) of considerable thickness above a rather uniform lithospheric mantle (3300 kgm⁻³). The spatial correlation between two major bends of the main Marmara fault and the location of the highdensity bodies suggests that the distribution of lithological heterogeneities within the crust controls the rheological behavior along the NAFZ and, consequently, maybe influences fault segmentation and thus the seismic hazard assessment in the region.
Variability of the Cold Season Climate in Central Asia. Part II: Hydroclimatic Predictability
(2019)
Central Asia (CA) is subjected to a large variability of precipitation. This study presents a statistical model, relating precipitation anomalies in three subregions of CA in the cold season (November-March) with various predictors in the preceding October. Promising forecast skill is achieved for two subregions covering 1) Uzbekistan, Turkmenistan, Kyrgyzstan, Tajikistan, and southern Kazakhstan and 2) Iran, Afghanistan, and Pakistan. ENSO in October is identified as the major predictor. Eurasian snow cover and the quasi-biennial oscillation further improve the forecast performance. To understand the physical mechanisms, an analysis of teleconnections between these predictors and the wintertime circulation over CA is conducted. The correlation analysis of predictors and large-scale circulation indices suggests a seasonal persistence of tropical circulation modes and a dynamical forcing of the westerly circulation by snow cover variations over Eurasia. An EOF analysis of pressure and humidity patterns allows separating the circulation variability over CA into westerly and tropical modes and confirms that the identified predictors affect the respective circulation characteristics. Based on the previously established weather type classification for CA, the predictors are investigated with regard to their effect on the regional circulation. The results suggest a modification of the Hadley cell due to ENSO variations, with enhanced moisture supply from the Arabian Gulf during El Nino. They further indicate an influence of Eurasian snow cover on the wintertime Arctic Oscillation (AO) and Northern Hemispheric Rossby wave tracks. Positive anomalies favor weather types associated with dry conditions, while negative anomalies promote the formation of a quasi-stationary trough over CA, which typically occurs during positive AO conditions.
An essential, respected, and critical aspect of the modern practice of science and scientific publishing is peer review. The process of peer review facilitates best practices in scientific conduct and communication, ensuring that manuscripts published are as accurate, valuable, and clearly communicated. The over 216 papers published in Tectonics in 2018 benefit from the time, effort, and expertise of our reviewers who have provided thoughtfully considered advice on each manuscript. This role is critical to advancing our understanding of the evolution of the continents and their margins, as these reviews lead to even clearer and higher-quality papers. In 2018, the over 443 papers submitted to Tectonics were the beneficiaries of more than 1,010 reviews provided by 668 members of the tectonics community and related disciplines. To everyone who has volunteered their time and intellect to peer reviewing, thank you for helping Tectonics and all other AGU Publications provide the best science possible.
Encroachment of shrubs into the unique pastoral grassland ecosystems of the Tibetan Plateau has significant impact on ecosystem services, especially forage production. We developed a process-based ecohydrological model to identify the relative importance of the main drivers of shrub encroachment for the alpine meadows within the Qinghai province. Specifically, we explored the effects of summer livestock grazing (intensity and type of livestock) together with the effects of climate warming, including interactions between herbaceous and woody vegetation and feedback loops between soil, water and vegetation. Under current climatic conditions and a traditional herd composition, an increasing grazing intensity above a threshold value of 0.32 +/- 0.10 large stock units (LSU) ha(-1) day(-1) changes the vegetation composition from herbaceous towards a woody and bare soil dominated system. Very high grazing intensity (above 0.8 LSU ha(-1) day(-1)) leads to a complete loss of any vegetation. Under warmer conditions, the vegetation showed a higher resilience against livestock farming. This resilience is enhanced when the herd has a higher browser : grazer ratio. A cooler climate has a shrub encroaching effect, whereas warmer conditions increase the cover of the herbaceous vegetation. This effect was primarily due to season length and an accompanied competitive loss of slower growing shrubs, rather than evaporative water loss leading to less soil water in deeper soil layers for deeper rooting shrubs. If climate warming is driving current shrub encroachment, we conclude it is only indirectly so. It would be manifest by an advancing shrubline and could be regarded as a climatic escape of specific shrub species such as Potentilla fruticosa. Under the recent high intensity of grazing, only herding by more browsing animals can potentially prevent both shrub encroachment and the complete loss of herbaceous vegetation.
Vermetid reefs and rocky shores are hot spots of biodiversity, often referred to as the subtropical equivalent of coral reefs. The development of the ecosystem depends on the activity of several reef builders, including red crustose coralline algae (CCA) such as Neogoniolithon brassica-florida. Despite its importance, little is known about Neogoniolithon sp. acclimation to rapid changes in light intensity and corresponding photosynthetic activity. To overcome the large spatial variability in the light field (due to location and the porous nature of the rocks) we grew Neogoniolithon sp. on glass slides and characterized its photosynthetic performance in response to various light intensities by following O-2 exchange and fluorescence parameters. This was also performed on rock-inhabiting thalli collected from the east Mediterranean basin. Generally, maximal photosynthetic rate was reached when Neogoniolithon sp. thalli grown under low illumination (such as in protected niches where the light intensity can be as low as 1% of surface illumination) were examined. When exposed to light intensities higher than those experienced during growth, Neogoniolithon sp. activates adaptive/protective mechanisms such as state transition and nonphotochemical fluorescence quenching and increases the dark respiration thereafter. We find that the Fv/Fm parameter (variable/maximal fluorescence) is not suitable to assess photosynthetic performance in Neogoniolithon sp. and propose using instead an alternative parameter recently developed. Our findings help to clarify why Neogoniolithon sp. is usually observed in shaded niches along the reef surfaces.
We analyzed the Lomas de Carabajal area in the intermontane Lerma valley of the Cordillera Oriental to assess the level of neotectonic activity in a densely populated region of northwestern Argentina. In this region, Plio-Pleistocene synorogenic conglomerates are deformed, locally associated with high-angle faults, and NNW-SSE oriented en-echelon folds characterized by wavelengths of < 1 km. The deformed Quaternary units follow the same pattern of deformation as observed in the underlying Neogene deposits; growth-strata geometries are observed near faults. This configuration is compatible with local left-lateral transpressional tectonism driven by ENE-WSW buttressing against the NW-oriented border of a Cretaceous extensional basin (Alemania sub-basin). Optically Stimulated Luminescence analysis of sandy-silty layers interbedded within the folded late Pleistocene conglomeratic sequence helps to determine uplift rates of 0.83-0.87 mm/a during the last 30-40 ka. Nearby the Lomas de Carabajal, a WNW-striking, 3-m-high fault scarp disrupts radiocarbon dated, 10-ka-old loessic deposits providing a Holocene mean uplift rate of 0.30 mm/a. Our data unambiguously show that shallow crustal deformation in the intermontane Lerma valley is ongoing; some of this deformation may be associated with seismicity. Our findings support the notion of temporally and spatially disparate deformation processes in the broken foreland of the northwestern Argentinean Andes.
We analyze trends in compound flooding resulting from high coastal water levels (HCWLs) and peak river discharge over northwestern Europe during 1901-2014. Compound peak discharge associated with 37 stream gauges with at least 70 years of record availability near the North and Baltic Sea coasts is used. Compound flooding is assessed using a newly developed index, compound hazard ratio, that compares the severity of river flooding associated with HCWL with the at-site, T-year (a flood with 1/T chance of being exceeded in any given year) fluvial peak discharge. Our findings suggest a spatially coherent pattern in the dependence between HCWL and river peaks and in compound flood magnitudes and frequency. For higher return levels, we find upward trends in compound hazard ratio frequency at midlatitudes (gauges from 47 degrees N to 60 degrees N) and downward trends along the high latitude (>60 degrees N) regions of northwestern Europe. Plain Language Summary Compound floods in delta areas, that is, the co-occurrence of high coastal water levels (HCWLs) and high river discharge, are a particular challenge for disaster management. Such events are caused by two distinct mechanisms: (1) HCWLs may affect river flows and water levels by backwater effects or by reversing the seaward flow of rivers, particularly in regions with elevation less than 10 m in northwestern Europe. (2) The correlation between HCWL and river flow peaks may also stem from a common meteorological driver. Severe storm periods may be associated with high winds leading to storm surges, and at the same time with high precipitation followed by inland flooding. Understanding the historical trends in compound flooding, owing to changes in relative sea levels, in river flooding and in the dependence between these two drivers, is essential for projecting future changes and disaster management. The risk assessment frameworks are often limited to assessing flood risk from a single driver only. We present a new approach to assess compound flood severity resulting from extreme coastal water level and peak river discharge. We find upward trends in compound flooding for midlatitude regions and downward trends for high latitudes in northwestern Europe.
Compound flooding, such as the co-occurrence of fluvial floods and extreme coastal water levels (CWL), may lead to significant impacts in densely-populated Low Elevation Coastal Zones. They may overstrain disaster management owing to the co-occurrence of inundation from rivers and the sea. Recent studies are limited by analyzing joint dependence between river discharge and either CWL or storm surges, and little is known about return levels of compound flooding, accounting for the covariance between drivers. Here, we assess the compound flood severity and identify hotspots for northwestern Europe during 1970–2014, using a newly developed Compound Hazard Ratio (CHR) that compares the severity of compound flooding associated with extreme CWL with the unconditional T-year fluvial peak discharge. We show that extreme CWL and stronger storms greatly amplify fluvial flood hazards. Our results, based on frequency analyses of observational records during 2013/2014’s winter storm Xaver, reveal that the river discharge of the 50-year compound flood is up to 70% larger, conditioned on the occurrence of extreme CWL, than that of the at-site peak discharge. For this event, nearly half of the stream gauges show increased flood hazards, demonstrating the importance of including the compounding effect of extreme CWL in river flood risk management.
Changes in the steepness of river profiles or abrupt vertical steps (i.e. waterfalls) are thought to be indicative of changes in erosion rates, lithology or other factors that affect landscape evolution. These changes are referred to as knickpoints or knickzones and are pervasive in bedrock river systems. Such features are thought to reveal information about landscape evolution and patterns of erosion, and therefore their locations are often reported in the geomorphic literature. It is imperative that studies reporting knickpoints and knickzones use a reproducible method of quantifying their locations, as their number and spatial distribution play an important role in interpreting tectonically active landscapes. In this contribution we introduce a reproducible knickpoint and knickzone extraction algorithm that uses river profiles transformed by integrating drainage area along channel length (the so-called integral or chi method). The profile is then statistically segmented and the differing slopes and step changes in the elevations of these segments are used to identify knickpoints, knickzones and their relative magnitudes. The output locations of identified knickpoints and knickzones compare favourably with human mapping: we test the method on Santa Cruz Island, CA, using previously reported knickzones and also test the method against a new dataset from the Quadrilatero Ferrifero in Brazil. The algorithm allows for the extraction of varying knickpoint morphologies, including stepped, positive slope-break (concave upward) and negative slope-break knickpoints. We identify parameters that most affect the resulting knickpoint and knickzone locations and provide guidance for both usage and outputs of the method to produce reproducible knickpoint datasets.
Thermokarst lake landscapes are permafrost regions, which are prone to rapid (on seasonal to decadal time scales) changes, affecting carbon and nitrogen cycles. However, there is a high degree of uncertainty related to the balance between carbon and nitrogen cycling and storage. We collected 12 permafrost soil cores from six drained thermokarst lake basins (DTLBs) along a chronosequence north of Teshekpuk Lake in northern Alaska and analyzed them for carbon and nitrogen contents. For comparison, we included three lacustrine cores from an adjacent thermokarst lake and one soil core from a non thermokarst affected remnant upland. This allowed to calculate the carbon and nitrogen stocks of the three primary landscape units (DTLB, lake, and upland), to reconstruct the landscape history, and to analyze the effect of thermokarst lake formation and drainage on carbon and nitrogen stocks. We show that carbon and nitrogen contents and the carbon-nitrogen ratio are considerably lower in sediments of extant lakes than in the DTLB or upland cores indicating degradation of carbon during thermokarst lake formation. However, we found similar amounts of total carbon and nitrogen stocks due to the higher density of lacustrine sediments caused by the lack of ground ice compared to DTLB sediments. In addition, the radiocarbon-based landscape chronology for the past 7,000years reveals five successive lake stages of partially, spatially overlapping DTLBs in the study region, reflecting the dynamic nature of ice-rich permafrost deposits. With this study, we highlight the importance to include these dynamic landscapes in future permafrost carbon feedback models. Plain Language Summary When permanently frozen soils (permafrost) contain ice-rich sediments, the thawing of this permafrost causes the surface to sink, which may result in lake formation. This process, the thaw of ice-rich permafrost and melting of ground ice leads to characteristic landforms-known as thermokarst. Once such a thaw process is initiated in ice-rich sediments, a thaw lake forms and grows by shoreline erosion, eventually expanding until a drainage pathway is encountered and the lake eventually drains, resulting in a drained thermokarst lake basin. In our study, we show that such a thermokarst-affected landscape north of Teshekpuk Lake in northern Alaska is shaped by repeated thaw lake formation and lake drainage events during the past 7,000years, highlighting the dynamic nature of these landscapes. These landscape-scale processes have a big effect on the carbon and nitrogen stored in permafrost soils. We show that large amounts of carbon (>45kg C/m(2)) and nitrogen (>2.6kg N/m(2)) are stored in unfrozen lake sediments and in frozen soil sediments. The findings are important when considering the potential effect that permafrost thaw has for the global climate through releasing carbon and nitrogen, which was frozen and therefore locked away for millennia, from the active carbon cycle.