Refine
Year of publication
Document Type
- Article (2810)
- Doctoral Thesis (507)
- Postprint (139)
- Other (75)
- Review (52)
- Monograph/Edited Volume (34)
- Preprint (17)
- Conference Proceeding (13)
- Habilitation Thesis (12)
- Master's Thesis (6)
Keywords
- climate change (52)
- Holocene (44)
- erosion (28)
- permafrost (27)
- Himalaya (26)
- remote sensing (24)
- Climate change (23)
- Tibetan Plateau (22)
- Earthquake source observations (21)
- Pollen (21)
Institute
- Institut für Geowissenschaften (3673) (remove)
Water research is introduced from the combined perspectives of natural and social science and cases of citizen and stakeholder coproduction of knowledge. Using the overarching notion of transdisciplinarity, we examine how interdisciplinary and participatory water research has taken place and could be developed further. It becomes apparent that water knowledge is produced widely within society, across certified disciplinary experts and noncertified expert stakeholders and citizens. However, understanding and management interventions may remain partial, or even conflicting, as much research across and between traditional disciplines has failed to integrate disciplinary paradigms due to philosophical, methodological, and communication barriers. We argue for more agonistic relationships that challenge both certified and noncertified knowledge productively. These should include examination of how water research itself embeds and is embedded in social context and performs political work. While case studies of the cultural and political economy of water knowledge exist, we need more empirical evidence on how exactly culture, politics, and economics have shaped this knowledge and how and at what junctures this could have turned out differently. We may thus channel the coproductionist critique productively to bring perspectives, alternative knowledges, and implications into water politics where they were not previously considered; in an attempt to counter potential lock-in to particular water policies and technologies that may be inequitable, unsustainable, or unacceptable. While engaging explicitly with politics, transdisciplinary water research should remain attentive to closing down moments in the research process, such as framings, path-dependencies, vested interests, researchers’ positionalities, power, and scale.
A tale of shifting relations
(2021)
Understanding the dynamics between the East Asian summer (EASM) and winter monsoon (EAWM) is needed to predict their variability under future global warming scenarios. Here, we investigate the relationship between EASM and EAWM as well as the mechanisms driving their variability during the last 10,000 years by stacking marine and terrestrial (non-speleothem) proxy records from the East Asian realm. This provides a regional and proxy independent signal for both monsoonal systems. The respective signal was subsequently analysed using a linear regression model. We find that the phase relationship between EASM and EAWM is not time-constant and significantly depends on orbital configuration changes. In addition, changes in the Atlantic Meridional Overturning circulation, Arctic sea-ice coverage, El Niño-Southern Oscillation and Sun Spot numbers contributed to millennial scale changes in the EASM and EAWM during the Holocene. We also argue that the bulk signal of monsoonal activity captured by the stacked non-speleothem proxy records supports the previously argued bias of speleothem climatic archives to moisture source changes and/or seasonality.
A tale of shifting relations
(2021)
Understanding the dynamics between the East Asian summer (EASM) and winter monsoon (EAWM) is needed to predict their variability under future global warming scenarios. Here, we investigate the relationship between EASM and EAWM as well as the mechanisms driving their variability during the last 10,000 years by stacking marine and terrestrial (non-speleothem) proxy records from the East Asian realm. This provides a regional and proxy independent signal for both monsoonal systems. The respective signal was subsequently analysed using a linear regression model. We find that the phase relationship between EASM and EAWM is not time-constant and significantly depends on orbital configuration changes. In addition, changes in the Atlantic Meridional Overturning circulation, Arctic sea-ice coverage, El Niño-Southern Oscillation and Sun Spot numbers contributed to millennial scale changes in the EASM and EAWM during the Holocene. We also argue that the bulk signal of monsoonal activity captured by the stacked non-speleothem proxy records supports the previously argued bias of speleothem climatic archives to moisture source changes and/or seasonality.
Sustainable development goals (SDGs) have set the 2030 agenda to transform our world by tackling multiple challenges humankind is facing to ensure well-being, economic prosperity, and environmental protection. In contrast to conventional development agendas focusing on a restricted set of dimensions, the SDGs provide a holistic and multidimensional view on development. Hence, interactions among the SDGs may cause diverging results. To analyze the SDG interactions we systematize the identification of synergies and trade-offs using official SDG indicator data for 227 countries. A significant positive correlation between a pair of SDG indicators is classified as a synergy while a significant negative correlation is classified as a trade-off. We rank synergies and trade-offs between SDGs pairs on global and country scales in order to identify the most frequent SDG interactions. For a given SDG, positive correlations between indicator pairs were found to outweigh the negative ones in most countries. Among SDGs the positive and negative correlations between indicator pairs allowed for the identification of particular global patterns. SDG 1 (No poverty) has synergetic relationship with most of the other goals, whereas SDG 12 (Responsible consumption and production) is the goal most commonly associated with trade-offs. The attainment of the SDG agenda will greatly depend on whether the identified synergies among the goals can be leveraged. In addition, the highlighted trade-offs, which constitute obstacles in achieving the SDGs, need to be negotiated and made structurally nonobstructive by deeper changes in the current strategies.
Most South Asian countries have challenges in ensuring water, energy, and food (WEF) security, which are often interacting positively or negatively. To address these challenges, the nexus approach provides a framework to identify the interactions of the WEF sectors as an integrated system. However, most nexus studies only qualitatively discuss the interactions between these sectors. This study conducts a systematic analysis of the WEF security nexus in South Asia by using open data sources at the country scale. We analyze interactions between the WEF sectors statistically, defining positive and negative correlations between the WEF security indicators as synergies and trade-offs, respectively. By creating networks of the synergies and trade-offs, we further identify most positively and negatively influencing indicators in the WEF security nexus. We observe a larger share of trade-offs than synergies within the water and energy sectors and a larger share of synergies than trade-offs among the WEF sectors for South Asia. However, these observations vary across the South Asian countries. Our analysis highlights that strategies on promoting sustainable energy and discouraging fossil fuel use could have overall positive effects on the WEF security nexus in the countries. This study provides evidence for considering the WEF security nexus as an integrated system rather than just a combination of three different sectors or securities.
High-pressure is a key feature of deep subsurface environments. High partial pressure of dissolved gasses plays an important role in microbial metabolism, because thermodynamic feasibility of many reactions depends on the concentration of reactants. For gases, this is controlled by their partial pressure, which can exceed 1 MPa at in situ conditions. Therefore, high hydrostatic pressure alone is not sufficient to recreate true deep subsurface in situ conditions, but the partial pressure of dissolved gasses has to be controlled as well. We developed an incubation system that allows for incubations at hydrostatic pressure up to 60 MPa, temperatures up to 120 degrees C, and at high gas partial pressure. The composition and partial pressure of gasses can be manipulated during the experiment. To keep costs low, the system is mainly made from off-the-shelf components with only very few custommade parts. A flexible and inert PVDF (polyvinylidene fluoride) incubator sleeve, which is almost impermeable for gases, holds the sample and separates it from the pressure fluid. The flexibility of the incubator sleeve allows for sub-sampling of the medium without loss of pressure. Experiments can be run in both static and flow-through mode. The incubation system described here is usable for versatile purposes, not only the incubation of microorganisms and determination of growth rates, but also for chemical degradation or extraction experiments under high gas saturation, e.g., fluid-gas-rock-interactions in relation to carbon dioxide sequestration. As an application of the system we extracted organic compounds from sub-bituminous coal using H2O as well as a H2O-CO2 mixture at elevated temperature (90 degrees C) and pressure (5 MPa). Subsamples were taken at different time points during the incubation and analyzed by ion chromatography. Furthermore we demonstrated the applicability of the system for studies of microbial activity, using samples from the Isis mud volcano. We could detect an increase in sulfate reduction rate upon the addition of methane to the sample.
Current evidence suggests that many of the major events in hominin evolution occurred in East Africa. Hence, over the past two decades, there has been intensive work undertaken to understand African palaeoclimate and tectonics in order to put together a coherent picture of how the environment of Africa has varied over the past 10 Myr. A new consensus is emerging that suggests the unusual geology and climate of East Africa created a complex, environmentally very variable setting. This new understanding of East African climate has led to the pulsed climate variability hypothesis that suggests the long-term drying trend in East Africa was punctuated by episodes of short alternating periods of extreme humidity and aridity which may have driven hominin speciation, encephalization and dispersals out of Africa. This hypothesis is unique as it provides a conceptual framework within which other evolutionary theories can be examined: first, at macro-scale comparing phylogenetic gradualism and punctuated equilibrium; second, at a more focused level of human evolution comparing allopatric speciation, aridity hypothesis, turnover pulse hypothesis, variability selection hypothesis, Red Queen hypothesis and sympatric speciation based on sexual selection. It is proposed that each one of these mechanisms may have been acting on hominins during these short periods of climate variability, which then produce a range of different traits that led to the emergence of new species. In the case of Homo erectus (sensu lato), it is not just brain size that changes but life history (shortened inter-birth intervals, delayed development), body size and dimorphism, shoulder morphology to allow thrown projectiles, adaptation to long-distance running, ecological flexibility and social behaviour. The future of evolutionary research should be to create evidence-based meta-narratives, which encompass multiple mechanisms that select for different traits leading ultimately to speciation.
Ambitious climate policies, as well as economic development, education, technological progress and less resource-intensive lifestyles, are crucial elements for progress towards the UN Sustainable Development Goals (SDGs). However, using an integrated modelling framework covering 56 indicators or proxies across all 17 SDGs, we show that they are insufficient to reach the targets. An additional sustainable development package, including international climate finance, progressive redistribution of carbon pricing revenues, sufficient and healthy nutrition and improved access to modern energy, enables a more comprehensive sustainable development pathway. We quantify climate and SDG outcomes, showing that these interventions substantially boost progress towards many aspects of the UN Agenda 2030 and simultaneously facilitate reaching ambitious climate targets. Nonetheless, several important gaps remain; for example, with respect to the eradication of extreme poverty (180 million people remaining in 2030). These gaps can be closed by 2050 for many SDGs while also respecting the 1.5 °C target and several other planetary boundaries.
A submerged pine forest from the early Holocene in the Mecklenburg Lake District, northern Germany
(2018)
For the first time, evidence of a submerged pine forest from the early Holocene can be documented in a central European lake. Subaquatic tree stumps were discovered in Lake Giesenschlagsee at a depth of between 2 and 5m using scuba divers, side-scan sonar and a remotely operated vehicle. Several erect stumps, anchored to the ground by roots, represent an insitu record of this former forest. Botanical determination revealed the stumps to be Scots pine (Pinus sylvestris) with an individual tree age of about 80years. The trees could not be dated by means of dendrochronology, as they are older than the regional reference chronology for pine. Radiocarbon ages from the wood range from 10880 +/- 210 to 10370 +/- 130cal. a BP, which is equivalent to the mid-Preboreal to early Boreal biozones. The trees are rooted in sedge peat, which can be dated to this period as well, using pollen stratigraphical analysis. Tilting of the peat bed by 4m indicates subsidence of the ground due to local dead ice melting, causing the trees to become submerged and preserved for millennia. Together with recently detected Lateglacial insitu tree occurrences in nearby lakes, the submerged pine forest at Giesenschlagsee represents a new and highly promising type of geo-bio-archive for the wider region. Comparable insitu pine remnants occur at some terrestrial (buried setting) and marine (submerged setting) sites in northern central Europe and beyond, but they partly differ in age. In general, the insitu pine finds document shifts of the zonal boreal forest ecosystem during the late Quaternary.
Response spectra are of fundamental importance in earthquake engineering and represent a standard measure in seismic design for the assessment of structural performance. However, unlike Fourier spectral amplitudes, the relationship of response spectral amplitudes to seismological source, path, and site characteristics is not immediately obvious and might even be considered counterintuitive for high oscillator frequencies. The understanding of this relationship is nevertheless important for seismic-hazard analysis. The purpose of the present study is to comprehensively characterize the variation of response spectral amplitudes due to perturbations of the causative seismological parameters. This is done by calculating the absolute parameter sensitivities (sensitivity coefficients) defined as the partial derivatives of the model output with respect to its input parameters. To derive sensitivities, we apply algorithmic differentiation (AD). This powerful approach is extensively used for sensitivity analysis of complex models in meteorology or aerodynamics. To the best of our knowledge, AD has not been explored yet in the seismic-hazard context. Within the present study, AD was successfully implemented for a proven and extensively applied simulation program for response spectra (Stochastic Method SIMulation [SMSIM]) using the TAPENADE AD tool. We assess the effects and importance of input parameter perturbations on the shape of response spectra for different regional stochastic models in a quantitative way. Additionally, we perform sensitivity analysis regarding adjustment issues of groundmotion prediction equations.
Solar wind observations show that geomagnetic storms are mainly driven by interplanetary coronal mass ejections (ICMEs) and corotating or stream interaction regions (C/SIRs). We present a binary classifier that assigns one of these drivers to 7,546 storms between 1930 and 2015 using ground‐based geomagnetic field observations only. The input data consists of the long‐term stable Hourly Magnetospheric Currents index alongside the corresponding midlatitude geomagnetic observatory time series. This data set provides comprehensive information on the global storm time magnetic disturbance field, particularly its spatial variability, over eight solar cycles. For the first time, we use this information statistically with regard to an automated storm driver identification. Our supervised classification model significantly outperforms unskilled baseline models (78% accuracy with 26[19]% misidentified interplanetary coronal mass ejections [corotating or stream interaction regions]) and delivers plausible driver occurrences with regard to storm intensity and solar cycle phase. Our results can readily be used to advance related studies fundamental to space weather research, for example, studies connecting galactic cosmic ray modulation and geomagnetic disturbances. They are fully reproducible by means of the underlying open‐source software (Pick, 2019, http://doi.org/10.5880/GFZ.2.3.2019.003)
Solar wind observations show that geomagnetic storms are mainly driven by interplanetary coronal mass ejections (ICMEs) and corotating or stream interaction regions (C/SIRs). We present a binary classifier that assigns one of these drivers to 7,546 storms between 1930 and 2015 using ground‐based geomagnetic field observations only. The input data consists of the long‐term stable Hourly Magnetospheric Currents index alongside the corresponding midlatitude geomagnetic observatory time series. This data set provides comprehensive information on the global storm time magnetic disturbance field, particularly its spatial variability, over eight solar cycles. For the first time, we use this information statistically with regard to an automated storm driver identification. Our supervised classification model significantly outperforms unskilled baseline models (78% accuracy with 26[19]% misidentified interplanetary coronal mass ejections [corotating or stream interaction regions]) and delivers plausible driver occurrences with regard to storm intensity and solar cycle phase. Our results can readily be used to advance related studies fundamental to space weather research, for example, studies connecting galactic cosmic ray modulation and geomagnetic disturbances. They are fully reproducible by means of the underlying open‐source software (Pick, 2019, http://doi.org/10.5880/GFZ.2.3.2019.003)
During the last years Urban Green Infrastructure (UGI) has evolved as a research focus across Europe. UGI can be understood as a multifunctional network of different urban green spaces and elements contributing to urban benefits. Urban agriculture has gained increasing research interest in this context. While a strong focus has been made on functions and benefits of small scale activities, the question is still open, whether these findings can be up-scaled and transferred to the farmland scale. Furthermore, multifunctionality of urban and peri-urban agriculture is rarely being considered in the landscape context. This research aims to address these gaps and harnesses the question if agricultural landscapes – which in many European metropolitan regions provide significant spatial potential – can contribute to UGI as multifunctional green spaces. This work considers multifunctionality qualitatively based on stakeholder opinion, using a participatory research approach. This study provides new insights in peri-urban farmland potentials for UGI development, resulting into a strategy framework. Furthermore, it reflects on the role of the stakeholder involvement for `multifunctionality planning´. It suggests that it helps to define meaningful bundles of intertwined functions that interact on different scales, helping to deal with non-linearity of multiple functions and to better manage them simultaneously.
A SSHAC Level 3 Probabilistic Seismic Hazard Analysis for a New-Build Nuclear Site in South Africa
(2015)
A probabilistic seismic hazard analysis has been conducted for a potential nuclear power plant site on the coast of South Africa, a country of low-to-moderate seismicity. The hazard study was conducted as a SSHAC Level 3 process, the first application of this approach outside North America. Extensive geological investigations identified five fault sources with a non-zero probability of being seismogenic. Five area sources were defined for distributed seismicity, the least active being the host zone for which the low recurrence rates for earthquakes were substantiated through investigations of historical seismicity. Empirical ground-motion prediction equations were adjusted to a horizon within the bedrock at the site using kappa values inferred from weak-motion analyses. These adjusted models were then scaled to create new equations capturing the range of epistemic uncertainty in this region with no strong motion recordings. Surface motions were obtained by convolving the bedrock motions with site amplification functions calculated using measured shear-wave velocity profiles.
Mineral resource exploration and mining is an essential part of today's high-tech industry. Elements such as rare-earth elements (REEs) and copper are, therefore, in high demand. Modern exploration techniques from multiple platforms (e.g., spaceborne and airborne), to detect and map the spectral characteristics of the materials of interest, require spectral libraries as an essential reference. They include field and laboratory spectral information in combination with geochemical analyses for validation. Here, we present a collection of REE- and copper-related hyperspectral spectra with associated geochemical information. The libraries contain reflectance spectra from rare-earth element oxides, REE-bearing minerals, copper-bearing minerals and mine surface samples from the Apliki copper-gold-pyrite mine in the Republic of Cyprus. The samples were measured with the HySpex imaging spectrometers in the visible and near infrared (VNIR) and shortwave infrared (SWIR) range (400-2500 nm). The geochemical validation of each sample is provided with the reflectance spectra. The spectral libraries are openly available to assist future mineral mapping campaigns and laboratory spectroscopic analyses. The spectral libraries and corresponding geochemistry are published via GFZ Data Services with the following DOIs: https://doi.org/10.5880/GFZ.1.4.2019.004 (13 REE-bearing minerals and 16 oxide powders, Koerting et al., 2019a), https://doi.org/10.5880/GFZ.1.4.2019.003 (20 copper-bearing minerals, Koellner et al., 2019), and https://doi.org/10.5880/GFZ.1.4.2019.005 (37 copper-bearing surface material samples from the Apliki coppergold-pyrite mine in Cyprus, Koerting et al., 2019b). All spectral libraries are united and comparable by the internally consistent method of hyperspectral data acquisition in the laboratory.
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.
Processes driving the production, transformation and transport of methane (CH4 / in wetland ecosystems are highly complex. We present a simple calculation algorithm to separate open-water CH4 fluxes measured with automatic chambers into diffusion-and ebullition-derived components. This helps to reveal underlying dynamics, to identify potential environmental drivers and, thus, to calculate reliable CH4 emission estimates. The flux separation is based on identification of ebullition-related sudden concentration changes during single measurements. Therefore, a variable ebullition filter is applied, using the lower and upper quartile and the interquartile range (IQR). Automation of data processing is achieved by using an established R script, adjusted for the purpose of CH4 flux calculation. The algorithm was validated by performing a laboratory experiment and tested using flux measurement data (July to September 2013) from a former fen grassland site, which converted into a shallow lake as a result of rewetting. Ebullition and diffusion contributed equally (46 and 55 %) to total CH4 emissions, which is comparable to ratios given in the literature. Moreover, the separation algorithm revealed a concealed shift in the diurnal trend of diffusive fluxes throughout the measurement period. The water temperature gradient was identified as one of the major drivers of diffusive CH4 emissions, whereas no significant driver was found in the case of erratic CH4 ebullition events.
The investigation of stresses, faults, structure and seismic hazards requires a good understanding and mapping of earthquake rupture and slip. Constraining the finite source of earthquakes from seismic and geodetic waveforms is challenging because the directional effects of the rupture itself are small and dynamic numerical solutions often include a large number of free parameters. The computational effort is large and therefore difficult to use in an exploratory forward modelling or inversion approach. Here, we use a simplified self-similar fracture model with only a few parameters, where the propagation of the fracture front is decoupled from the calculation of the slip. The approximative method is flexible and computationally efficient. We discuss the strengths and limitations of the model with real-case examples of well-studied earthquakes. These include the M-w 8.3 2015 Illapel, Chile, megathrust earthquake at the plate interface of a subduction zone and examples of continental intraplate strike-slip earthquakes like the M-w 7.1 2016 Kumamoto, Japan, multisegment variable slip event or the M-w 7.5 2018 Palu, Indonesia, supershear earthquake. Despite the simplicity of the model, a large number of observational features ranging from different rupture-front isochrones and slip distributions to directional waveform effects or high slip patches are easy to model. The temporal evolution of slip rate and rise time are derived from the incremental growth of the rupture and the stress drop without imposing other constraints. The new model is fast and implemented in the open-source Python seismology toolbox Pyrocko, ready to study the physics of rupture and to be used in finite source inversions.
We present a new, seismologically consistent expression for the total area and volume of populations of earthquake-triggered landslides. This model builds on a set of scaling relationships between key parameters, such as landslide spatial density, seismic ground acceleration, fault length, earthquake source depth, and seismic moment. To assess the model we have assembled and normalized a catalog of landslide inventories for 40 shallow, continental earthquakes. Low landscape steepness causes systematic overprediction of the total area and volume of landslides. When this effect is accounted for, the model predicts the total landslide volume of 63% of 40 cases to within a factor 2 of the volume estimated from observations (R-2 = 0.76). The prediction of total landslide area is also sensitive to the landscape steepness, but less so than the total volume, and it appears to be sensitive to controls on the landslide size-frequency distribution, and possibly the shaking duration. Some outliers are likely associated with exceptionally strong rock mass in the epicentral area, while others may be related to seismic source complexities ignored by the model. However, the close match between prediction and estimate for about two thirds of cases in our database suggests that rock mass strength is similar in many cases and that our simple seismic model is often adequate, despite the variety of lithologies and tectonic settings covered. This makes our expression suitable for integration into landscape evolution models and application to the anticipation or rapid assessment of secondary hazards associated with earthquakes.
Pseudotachylyte veins frequently associated with mylonites and ultramylonites occur within migmatitic paragneisses, metamonzodiorites, as well as felsic and mafic granulites at the base of the section of the Hercynian lower crust exposed in Calabria (Southern Italy). The crustal section is tectonically superposed on lower grade units. Ultramylonites and pseudotachylytes are particularly well developed in migmatitic paragneisses, whereas sparse fault-related pseudotachylytes and thin mylonite/ultramylonite bands occur in granulite-facies rocks. The presence of sillimanite and clinopyroxene in ultramylonites and mylonites indicates that relatively high-temperature conditions preceded the formation of pseudotachylytes. We have analysed pseudotachylytes from different rock types to ascertain their deep crustal origin and to better understand the relationships between brittle and ductile processes during deformation of the deeper crust. Different protoliths were selected to test how lithology controls pseudotachylyte composition and textures. In migmatites and felsic granulites, euhedral or cauliflower-shaped garnets directly crystallized from pseudotachylyte melts of near andesitic composition. This indicates that pseudotachylytes originated at deep crustal conditions (> 0.75 GPa). In mafic protoliths, quenched needle-to-feather-shaped high-alumina orthopyroxene occurs in contact with newly crystallized plagioclase. The pyroxene crystallizes in garnet-free and garnet-bearing veins. The simultaneous growth of orthopyroxene and plagioclase as well as almandine, suggests lower crustal origin, with pressures in excess of 0.85 GPa. The existence of melts of different composition in the same vein indicates the stepwise, non-equilibrium conditions of frictional melting. Melt formed and intruded into pre-existing anisotropies. In mafic granulites, brittle faulting is localized in a previously formed thin high-temperature mylonite bands. migmatitic gneisses are deformed into ultramylonite domains characterized by s-c fabric. Small grain size and fluids lowered the effective stress on the c planes favouring a seismic event and the consequent melt generation. Microstructures and ductile deformation of pseudotachylytes suggest continuous ductile flow punctuated by episodes of high-strain rate, leading to seismic events and melting.
Volcanic eruptions are often preceded by seismic activity that can be used to quantify the volcanic activity. In order to allow consistent inference of the volcanic activity state from the observed seismicity patterns, objective and time-invariant classification results achievable by automatic systems should be preferred. Most automatic classification approaches need a large preclassified data set for training the system. However, in case of a volcanic crisis, we are often confronted with a lack of training data due to insufficient prior observations. In the worst case (e. g., volcanic crisis related reconfiguration of stations), there are even no prior observations available. Finally, due to the imminent crisis there might be no time for the time-consuming process of preparing a training data set. For this reason, we have developed a novel seismic-event spotting technique in order to be less dependent on previously acquired data bases and classification schemes. We are using a learning-while-recording approach based on a minimum number of reference waveforms, thus allowing for the build-up of a classification scheme as early as interesting events have been identified. First, short-term wave-field parameters (here, polarization and spectral attributes) are extracted from a continuous seismic data stream. The sequence of multidimensional feature vectors is then used to identify a fixed number of clusters in the feature space. Based on this general description of the overall wave field by a mixture of multivariate Gaussians, we are able to learn particular event classifiers (here, hidden Markov models) from a single waveform example. To show the capabilities of this new approach we apply the algorithm to a data set recorded at Soufriere Hills volcano, Montserrat. Supported by very high classification rates, we conclude that the suggested approach provides a valuable tool for volcano monitoring systems.
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.
Cosmic-Ray neutron sensing (CRS) is a unique approach to measure soil moisture at field scale filling the gap of current methodologies. However, CRS signal is affected by all the hydrogen pools on the land surface and understanding their relative importance plays an important role for the application of the method e.g., validation of remote sensing products and data assimilation. In this study, a soil moisture scaling approach is proposed to estimate directly the correct CRS soil moisture based on the soil moisture profile measured at least in one position within the field. The approach has the advantage to avoid the need to introduce one correction for each hydrogen contribution and to estimate indirectly all the related time-varying hydrogen pools. Based on the data collected in three crop seasons, the scaling approach shows its ability to identify and to quantify the seasonal biomass water equivalent. Additionally, the analysis conducted at sub-daily time resolution is able to quantify the daily vertical redistribution of the water biomass and the rainfall interception, showing promising applications of the CRS method also for these types of measurements. Overall, the study underlines how not only soil moisture but all the specific hydrological processes in the soil-plant-atmosphere continuum should be considered for a proper evaluation of the CRS signal. For this scope, the scaling approach reveals to be a simple and pragmatic analysis that can be easily extended to other experimental sites. (C) 2015 Elsevier B.V. All rights reserved.
Groundwater recharge (GWR) is one of the most challenging water fluxes to estimate, as it relies on observed data that are often limited in many developing countries.
This study developed an innovative water budget method using satellite products for estimating the spatially distributed GWR at monthly and annual scales in tropical wet sedimentary regions despite cloudy conditions.
The distinctive features proposed in this study include the capacity to address 1) evapotranspiration estimations in tropical wet regions frequently overlaid by substantial cloud cover; and 2) seasonal root-zone water storage estimations in sedimentary regions prone to monthly variations.
The method also utilises satellite-based information of the precipitation and surface runoff. The GWR was estimated and validated for the hydrologically contrasting years 2016 and 2017 over a tropical wet sedimentary region located in North-eastern Brazil, which has substantial potential for groundwater abstraction.
This study showed that applying a cloud-cleaning procedure based on monthly compositions of biophysical data enables the production of a reasonable proxy for evapotranspiration able to estimate groundwater by the water budget method.
The resulting GWR rates were 219 (2016) and 302 (2017) mm yr(-1), showing good correlations (CC = 0.68 to 0.83) and slight underestimations (PBIAS =-13 to-9%) when compared with the referenced estimates obtained by the water table fluctuation method for 23 monitoring wells. Sensitivity analysis shows that water storage changes account for +19% to-22% of our monthly evaluation.
The satellite-based approach consistently demonstrated that the consideration of cloud-cleaned evapotranspiration and root-zone soil water storage changes are essential for a proper estimation of spatially distributed GWR in tropical wet sedimentary regions because of their weather seasonality and cloudy conditions.
The Canadian database on contemporary crustal stress has not been revised systematically in the past two decades. Here we present the results of our new compilation that contains 514 new data records for the orientation data of maximum compressive horizontal stress and 188 data records that were re-assessed. In total the Canadian stress database has now 1667 data records, which is an increase of about 45%. From these data, a new Canadian Stress map as well as one for the Province of Alberta is presented.
To analyse the stress pattern, we use the quasi median on the circle as a smoothing algorithm that generates a smoothed stress map of the maximum compressive horizontal stress orientation on a regular grid. The newly introduced quasi interquartile range on the circle estimates the spreading of the data and is used as a measure for the wave-length of the stress pattern. The result of the hybrid wavelength analysis confirms that long spatial wavelength stress patterns (>= 1000 km) exist in large areas in Canada. The observed stress pattern is transmitted through the intra-plate regions.
The results reveal that shorter spatial wave length variation of the maximum compressive horizontal stress orientation of less than 200 km, prevails particularly in south-eastern and western Canada. Regional stress sources such as density contrasts, active fault systems, crustal structures, etc. might have a significant impact in these regions. In contrast to these variations, the observed stress pattern in the Alberta Basin is very homogeneous and mainly controlled by plate boundary forces and body forces. The influence of curvature of the Rocky Mountains salient in southern Alberta is minimal. The present-day horizontal stress orientations determined herein have important implications for the production of hydrocarbons and geothermal energy in the Alberta Basin. (C) 2014 Elsevier B.V. All rights reserved.
Like almost all fields of science, hydrology has benefited to a large extent from the tremendous improvements in scientific instruments that are able to collect long-time data series and an increase in available computational power and storage capabilities over the last decades. Many model applications and statistical analyses (e.g., extreme value analysis) are based on these time series. Consequently, the quality and the completeness of these time series are essential. Preprocessing of raw data sets by filling data gaps is thus a necessary procedure. Several interpolation techniques with different complexity are available ranging from rather simple to extremely challenging approaches. In this paper, various imputation methods available to the hydrological researchers are reviewed with regard to their suitability for filling gaps in the context of solving hydrological questions. The methodological approaches include arithmetic mean imputation, principal component analysis, regression-based methods and multiple imputation methods. In particular, autoregressive conditional heteroscedasticity (ARCH) models which originate from finance and econometrics will be discussed regarding their applicability to data series characterized by non-constant volatility and heteroscedasticity in hydrological contexts. The review shows that methodological advances driven by other fields of research bear relevance for a more intensive use of these methods in hydrology. Up to now, the hydrological community has paid little attention to the imputation ability of time series models in general and ARCH models in particular.
A review of source models to further the understanding of the seismicity of the Groningen field
(2022)
The occurrence of felt earthquakes due to gas production in Groningen has initiated numerous studies and model attempts to understand and quantify induced seismicity in this region. The whole bandwidth of available models spans the range from fully deterministic models to purely empirical and stochastic models. In this article, we summarise the most important model approaches, describing their main achievements and limitations. In addition, we discuss remaining open questions and potential future directions of development.
Damage due to floods has increased during the last few decades, and further increases are expected in several regions due to climate change and growing vulnerability. To address the projected increase in flood risk, a combination of structural and non-structural flood risk mitigation measures is considered as a promising adaptation strategy. Such a combination takes into account that flood defence systems may fail, and prepares for unexpected crisis situations via land-use planning and private damage reduction, e.g. via building precautionary measures, and disaster response. However, knowledge about damage-reducing measures is scarce and often fragmented since based on case studies. For instance, it is believed that private precautionary measures, like shielding with water shutters or building fortification, are especially effective in areas with frequent flood events and low flood water levels. However, some of these measures showed a significant damage-reducing effect also during the extreme flood event in 2002 in Germany. This review analyses potentials of land-use planning and private flood precautionary measures as components of adaptation strategies for global change. Focus is on their implementation, their damage-reducing effects and their potential contribution to address projected changes in flood risk, particularly in developed countries.
The Seismic Hazard Inferred from Tectonics based on the Global Strain Rate Map (SHIFT_GSRM) earthquake forecast was designed to provide high-resolution estimates of global shallow seismicity to be used in seismic hazard assessment. This model combines geodetic strain rates with global earthquake parameters to characterize long-term rates of seismic moment and earthquake activity. Although SHIFT_GSRM properly computes seismicity rates in seismically active continental regions, it underestimates earthquake rates in subduction zones by an average factor of approximately 3. We present a complementary method to SHIFT_GSRM to more accurately forecast earthquake rates in 37 subduction segments, based on the conservation of moment principle and the use of regional interface seismicity parameters, such as subduction dip angles, corner magnitudes, and coupled seismogenic thicknesses. In seven progressive steps, we find that SHIFT_GSRM earthquake-rate underpredictions are mainly due to the utilization of a global probability function of seismic moment release that poorly captures the great variability among subduction megathrust interfaces. Retrospective test results show that the forecast is consistent with the observations during the 1 January 1977 to 31 December 2014 period. Moreover, successful pseudoprospective evaluations for the 1 January 2015 to 31 December 2018 period demonstrate the power of the regionalized earthquake model to properly estimate subduction-zone seismicity.
We report on a receiver function study of the crust and upper mantle within DESERT, a multidisciplinary geophysical project to study the lithosphere across the Dead Sea Transform (DST). A temporary seismic network was operated on both sides of the DST between 2000 April and 2001 June. The depth of the Moho increases smoothly from about 30 to 34-38 km towards the east across the DST, with significant north-south variations east of the DST. These Moho depth estimates from receiver functions are consistent with results from steep-and wide-angle controlled-source techniques. Steep-angle reflections and receiver functions reveal an additional discontinuity in the lower crust, but only east of the DST. This leads to the conclusion that the internal crustal structure east and west of the DST is different. The P to S converted phases from both discontinuities at 410 and 660 km are delayed by 2 s with respect to the IASP91 global reference model. This would indicate that the transition zone is consistent with the global average, but the upper mantle above 410 km is 3-4 per cent slower than the standard earth model
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.
We propose a RAndom Interacting Network (RAIN) model to study the interactions between a pair of complex networks. The model involves two major steps: (i) the selection of a pair of nodes, one from each network, based on intra-network node-based characteristics, and (ii) the placement of a link between selected nodes based on the similarity of their relative importance in their respective networks. Node selection is based on a selection fitness function and node linkage is based on a linkage probability defined on the linkage scores of nodes. The model allows us to relate within-network characteristics to between-network structure. We apply the model to the interaction between the USA and Schengen airline transportation networks (ATNs). Our results indicate that two mechanisms: degree-based preferential node selection and degree-assortative link placement are necessary to replicate the observed inter-network degree distributions as well as the observed inter-network assortativity. The RAIN model offers the possibility to test multiple hypotheses regarding the mechanisms underlying network interactions. It can also incorporate complex interaction topologies. Furthermore, the framework of the RAIN model is general and can be potentially adapted to various real-world complex systems.
Sea-surface temperatures (SSTs) based on radiolarian assemblage changes are estimated for the last 160 kyr, from a sediment core (Y9) recovered from Pukaki Saddle, northeast of Campbell Plateau. Site Y9 lies beneath Subantarctic Surface Water (SAW) immediately to the north of the Subantarctic Front (SAF), which in this region is bathymetrically constrained by the edges of Campbell Plateau and defines the northern boundary of the Antarctic Circumpolar Current (ACC). Radiolarian assemblages are characterised by an exceptionally high abundance of the Antarctic to subantarctic species Antarctissa spp. (up to 68%), especially during glacial intervals. SST estimates are derived using Factor Analysis and the Modern Analog Technique. Both methods capture the glacial-interglacial (G-I) pattern. The SST reconstructions show the changing relative influence of distinct water masses during the past G-I cycle, with major temperature variations of the order of 7-9 degrees C at glacial Terminations. Glacials (marine isotope stages (MIS) 6 and 2) are associated with particularly cool SSTs that are indicative of a more vigorous SAF/ACC and an enhancement of the inflow through Pukaki Saddle and/or frequent development of cold-core eddies at the SAF. By contrast, the influence of warmer waters and relaxation of the ACC during interglacials can be inferred from temperatures slightly warmer (e.g., mid-Holocene) and/or comparable to present day (e.g., MIS 5e). During these intervals, relatively warmer temperatures most likely indicate a higher warmcore eddy activity due to a strengthened Subtropical Front and/or a weakened inflow of cool water through Pukaki Saddle and/or an increased stratification in the Campbell Plateau region. Furthermore, the SST record is characterised by an abrupt warming at ca. 10 kyr (i.e., Termination l), the occurrence of a reversal at Termination I, and a warming event at the end of MIS 4, coinciding with the A4 event in the Byrd ice core. These characteristics, together with the pronounced G-I cycle shown by the SST estimates, suggest that Site Y9 is influenced by major oceanographic changes in the SW Pacific and responds to thermal changes at high southern latitudes. (C) 2014 Elsevier B.V. All rights reserved.
Dynamic earthquake rupture modeling provides information on the rupture physics as the rupture velocity, frictions or tractions acting during the rupture process. Nevertheless, as often based on spatial gridded preset geometries, dynamic modeling is depending on many free parameters leading to both a high non-uniqueness of the results and large computation times. That decreases the possibilities of full Bayesian error analysis.
To assess the named problems we developed the quasi-dynamic rupture model which is presented in this work. It combines the kinematic Eikonal rupture model with a boundary element method for quasi-static slip calculation.
The orientation of the modeled rupture plane is defined by a previously performed moment tensor inversion. The simultanously inverted scalar seismic moment allows an estimation of the extension of the rupture. The modeled rupture plane is discretized by a set of rectangular boundary elements. For each boundary element an applied traction vector is defined as the boundary value.
For insights in the dynamic rupture behaviour the rupture front propagation is calculated for incremental time steps based on the 2D Eikonal equation. The needed location-dependent rupture velocity field is assumed to scale linearly with a layered shear wave velocity field.
At each time all boundary elements enclosed within the rupture front are used to calculate the quasi-static slip distribution. Neither friction nor stress propagation are considered. Therefore the algorithm is assumed to be “quasi-static”. A series of the resulting quasi-static slip snapshots can be used as a quasi-dynamic model of the rupture process.
As many a priori information is used from the earth model (shear wave velocity and elastic parameters) and the moment tensor inversion (rupture extension and orientation) our model is depending on few free parameters as the traction field, the linear factor between rupture and shear wave velocity and the nucleation point and time. Hence stable and fast modeling results are obtained as proven from the comparison to different infinite and finite static crack solutions.
First dynamic applications show promissing results. The location-dependent rise time is automatically derived by the model. Different simple kinematic models as the slip-pulse or the penny-shaped crack model can be reproduced as well as their corresponding slip rate functions. A source time function (STF) approximation calculated from the cumulative sum of moment rates of each boundary element gives results similar to theoretical and empirical known STFs.
The model was also applied to the 2015 Illapel earthquake. Using a simple rectangular rupture geometry and a 2-layered traction regime yields good estimates of both the rupture front propagation and the slip patterns which are comparable to literature results. The STF approximation shows a good fit with previously published STFs.
The quasi-dynamic rupture model is hence able to fastly calculate reproducable slip results. That allows to test full Bayesian error analysis in the future. Further work on a full seismic source inversion or even a traction field inversion can also extend the scope of our model.
Investigation of transient soil moisture profiles yields valuable information of near- surface processes. A recently developed reconstruction algorithm based on the telegraph equation allows the inverse estimation of soil moisture profiles along coated, three rod TDR probes. Laboratory experiments were carried out to prove the results of the inversion and to understand the influence of probe rod deformation and solid objects close to the probe in heterogonous media. Differences in rod geometry can lead to serious misinterpretations in the soil moisture profile but have small influence on the average soil moisture along the probe. Solids in the integration volume have almost no effect on average soil moisture but result in locally slightly decreased moisture values. Inverted profiles obtained in a loamy soil with a clay content of about 16% were in good agreement with independent measurements.
Investigation of transient soil moisture profiles yields valuable information of near- surface processes. A recently developed reconstruction algorithm based on the telegraph equation allows the inverse estimation of soil moisture profiles along coated, three rod TDR probes. Laboratory experiments were carried out to prove the results of the inversion and to understand the influence of probe rod deformation and solid objects close to the probe in heterogeneous media. Differences in rod geometry can lead to serious misinterpretations in the soil moisture profile, but have small influence on the average soil moisture along the probe. Solids in the integration volume have almost no effect on average soil moisture, but result in locally slightly decreased moisture values. Inverted profiles obtained in a loamy soil with a clay content of about 16% were in good agreement with independent measurements.
Python is at the forefront of scientific computation for seismologists and therefore should be introduced to students interested in becoming seismologists. On its own, Python is open source and well designed with extensive libraries. However, Python code can also be executed, visualized, and communicated to others with "Jupyter Notebooks". Thus, Jupyter Notebooks are ideal for teaching students Python and scientific computation. In this article, we designed an openly available Python library and collection of Jupyter Notebooks based on defined scientific computation learning goals for seismology students. The Notebooks cover topics from an introduction to Python to organizing data, earthquake catalog statistics, linear regression, and making maps. Our Python library and collection of Jupyter Notebooks are meant to be used as course materials for an upper-division data analysis course in an Earth Science Department, and the materials were tested in a Probabilistic Seismic Hazard course. However, seismologists or anyone else who is interested in Python for data analysis and map making can use these materials.
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.
Extracellular DNA (eDNA) is a ubiquitous biological compound in aquatic sediment and soil. Previous studies suggested that eDNA plays an important role in biogeochemical element cycling, horizontal gene transfer and stabilization of biofilm structures. Previous methods for eDNA extraction were either not suitable for oligotrophic sediments or only allowed quantification but no genetic analyses. Our procedure is based on cell detachment and eDNA liberation from sediment particles by sequential washing with an alkaline sodium phosphate buffer followed by a separation of cells and eDNA. The separated eDNA is then bound onto silica particles and purified, whereas the intracellular DNA from the separated cells is extracted using a commercial kit. The method provides extra- and intracellular DNA of high purity that is suitable for downstream applications like PCR. Extracellular DNA was extracted from organic-rich shallow sediment of the Baltic Sea, glacially influenced sediment of the Barents Sea and from the oligotrophic South Pacific Gyre. The eDNA concentration in these samples varied from 23 to 626 ng g(-1) wet weight sediment. A number of experiments were performed to verify each processing step. Although extraction efficiency is higher than other published methods, it is not fully quantitative. (C) 2014 Elsevier B.V. All rights reserved.
Extreme weather events are likely to occur more often under climate change and the resulting effects on ecosystems could lead to a further acceleration of climate change. But not all extreme weather events lead to extreme ecosystem response. Here, we focus on hazardous ecosystem behaviour and identify coinciding weather conditions. We use a simple probabilistic risk assessment based on time series of ecosystem behaviour and climate conditions. Given the risk assessment terminology, vulnerability and risk for the previously defined hazard are estimated on the basis of observed hazardous ecosystem behaviour.
We apply this approach to extreme responses of terrestrial ecosystems to drought, defining the hazard as a negative net biome productivity over a 12-month period. We show an application for two selected sites using data for 1981-2010 and then apply the method to the pan-European scale for the same period, based on numerical modelling results (LPJmL for ecosystem behaviour; ERA-Interim data for climate).
Our site-specific results demonstrate the applicability of the proposed method, using the SPEI to describe the climate condition. The site in Spain provides an example of vulnerability to drought because the expected value of the SPEI is 0.4 lower for hazardous than for non-hazardous ecosystem behaviour. In northern Germany, on the contrary, the site is not vulnerable to drought because the SPEI expectation values imply wetter conditions in the hazard case than in the non-hazard case.
At the pan-European scale, ecosystem vulnerability to drought is calculated in the Mediterranean and temperate region, whereas Scandinavian ecosystems are vulnerable under conditions without water shortages. These first model- based applications indicate the conceptual advantages of the proposed method by focusing on the identification of critical weather conditions for which we observe hazardous ecosystem behaviour in the analysed data set. Application of the method to empirical time series and to future climate would be important next steps to test the approach.
Extreme weather events are likely to occur more often under climate change and the resulting effects on ecosystems could lead to a further acceleration of climate change. But not all extreme weather events lead to extreme ecosystem response. Here, we focus on hazardous ecosystem behaviour and identify coinciding weather conditions. We use a simple probabilistic risk assessment based on time series of ecosystem behaviour and climate conditions. Given the risk assessment terminology, vulnerability and risk for the previously defined hazard are estimated on the basis of observed hazardous ecosystem behaviour.
We apply this approach to extreme responses of terrestrial ecosystems to drought, defining the hazard as a negative net biome productivity over a 12-month period. We show an application for two selected sites using data for 1981-2010 and then apply the method to the pan-European scale for the same period, based on numerical modelling results (LPJmL for ecosystem behaviour; ERA-Interim data for climate).
Our site-specific results demonstrate the applicability of the proposed method, using the SPEI to describe the climate condition. The site in Spain provides an example of vulnerability to drought because the expected value of the SPEI is 0.4 lower for hazardous than for non-hazardous ecosystem behaviour. In northern Germany, on the contrary, the site is not vulnerable to drought because the SPEI expectation values imply wetter conditions in the hazard case than in the non-hazard case.
At the pan-European scale, ecosystem vulnerability to drought is calculated in the Mediterranean and temperate region, whereas Scandinavian ecosystems are vulnerable under conditions without water shortages. These first model- based applications indicate the conceptual advantages of the proposed method by focusing on the identification of critical weather conditions for which we observe hazardous ecosystem behaviour in the analysed data set. Application of the method to empirical time series and to future climate would be important next steps to test the approach.
Residential assets, comprising buildings and household contents, are a major source of direct flood losses. Existing damage models are mostly deterministic and limited to particular countries or flood types. Here, we compile building-level losses from Germany, Italy and the Netherlands covering a wide range of fluvial and pluvial flood events. Utilizing a Bayesian network (BN) for continuous variables, we find that relative losses (i.e. loss relative to exposure) to building structure and its contents could be estimated with five variables: water depth, flow velocity, event return period, building usable floor space area and regional disposable income per capita. The model's ability to predict flood losses is validated for the 11 flood events contained in the sample. Predictions for the German and Italian fluvial floods were better than for pluvial floods or the 1993 Meuse river flood. Further, a case study of a 2010 coastal flood in France is used to test the BN model's performance for a type of flood not included in the survey dataset. Overall, the BN model achieved better results than any of 10 alternative damage models for reproducing average losses for the 2010 flood. An additional case study of a 2013 fluvial flood has also shown good performance of the model. The study shows that data from many flood events can be combined to derive most important factors driving flood losses across regions and time, and that resulting damage models could be applied in an open data framework.
This study aims to establish, evaluate, and apply a modern pollen-climate transfer function from the transition zone between arctic tundra and light-needled taiga in Arctic Siberia. Lacustrine samples (n = 96) from the northern Siberian lowlands of Yakutia were collected along four north-to-south transects crossing the arctic forest line. Samples span a broad temperature and precipitation gradient (mean July temperature, T-July: 7.5-18.7 degrees C; mean annual precipitation, P-ann: 114-315 mm/yr). Redundancy analyses are used to examine the relationship between the modern pollen signal and corresponding vegetation types and climate. Performance of transfer functions for T-July and P-ann were cross-validated and tested for spatial autocorrelation effects. The root mean square errors of prediction are 1.67 degrees C for T-July and 40 mm/yr for P-ann. A climate reconstruction based on fossil pollen spectra from a Siberian Arctic lake sediment core spanning the Holocene yielded cold conditions for the Late Glacial (1-2 degrees C below present T-July). Warm and moist conditions were reconstructed for the early to mid Holocene (2 degrees C higher T-July than present), and climate conditions similar to modern ones were reconstructed for the last 4000 years. In conclusion, our modern pollen data set fills the gap of existing regional calibration sets with regard to the underrepresented Siberian tundra-taiga transition zone. The Holocene climate reconstruction indicates that the temperature deviation from modern values was only moderate despite the assumed Arctic sensitivity to present climate change.
Extra-tropical circulation systems impede poleward moisture advection by the Indian Summer Monsoon. In this context, the Himalayan range is believed to insulate the south Asian circulation from extra-tropical influences and to delineate the northern extent of the Indian Summer Monsoon in central Asia. Paleoclimatic evidence, however, suggests increased moisture availability in the Early Holocene north of the Himalayan range which is attributed to an intensification of the Indian Summer Monsoon. Nevertheless, mechanisms leading to a surpassing of the Himalayan range and the northern maximum extent of summer monsoonal influence remain unknown. Here we show that the Kunlun barrier on the northern Tibetan Plateau [similar to 36 degrees N] delimits Indian Summer Monsoon precipitation during the Holocene. The presence of the barrier relocates the insulation effect 1,000 km further north, allowing a continental low intensity branch of the Indian Summer Monsoon which is persistent throughout the Holocene. Precipitation intensities at its northern extent seem to be driven by differentiated solar heating of the Northern Hemisphere indicating dependency on energy-gradients rather than absolute radiation intensities. The identified spatial constraints of monsoonal precipitation will facilitate the prediction of future monsoonal precipitation patterns in Central Asia under varying climatic conditions.
A partially non-ergodic ground-motion prediction equation is estimated for Europe and the Middle East. Therefore, a hierarchical model is presented that accounts for regional differences. For this purpose, the scaling of ground-motion intensity measures is assumed to be similar, but not identical in different regions. This is achieved by assuming a hierarchical model, where some coefficients are treated as random variables which are sampled from an underlying global distribution. The coefficients are estimated by Bayesian inference. This allows one to estimate the epistemic uncertainty in the coefficients, and consequently in model predictions, in a rigorous way. The model is estimated based on peak ground acceleration data from nine different European/Middle Eastern regions. There are large differences in the amount of earthquakes and records in the different regions. However, due to the hierarchical nature of the model, regions with only few data points borrow strength from other regions with more data. This makes it possible to estimate a separate set of coefficients for all regions. Different regionalized models are compared, for which different coefficients are assumed to be regionally dependent. Results show that regionalizing the coefficients for magnitude and distance scaling leads to better performance of the models. The models for all regions are physically sound, even if only very few earthquakes comprise one region.
Early Ilerdian (Early Eocene, Shallow Benthic Zones 5 and 6) carbonate systems of the Pyrenees shelf were deposited after a time of severe climatic ('Paleocene-Eocene Thermal Maximum, PETM') and phylogenetic ('Larger Foraminifer Turnover') changes. They reflect the radiation of nummulitid, alveolinid, and orbitolitid larger foraminifera after remarkable biotic changes at the end of the Paleocene, and announce their subsequent flourishing in the Middle Eocene. A paleoenvironmental model for tropical carbonate environments of this particular time interval is provided herein. During the Early Ilerdian, the inner and middle ramp deposits from Minerve, Campo and Serraduy revealed the end-member of a tropical carbonate factory with carbonate production dominated by the end-members of biotically (photo-autotrophic skeletal) controlled and biotically induced carbonate precipitation. Inner platform environments are dominated by alveolinids and in part by orbitolitids, middle platform environments are dominated by nummulitids. Corals are present, but they do not form reefs, which is a typical feature for the Eocene. Nummulite shoal complexes, which are well-known from the Middle Eocene are also absent during the studied Early Ilerdian interval, which may reflect the early evolutionary stage of this group
The spread of shrubs in Namibian savannas raises questions about the resilience of these ecosystems to global change. This makes it necessary to understand the past dynamics of the vegetation, since there is no consensus on whether shrub encroachment is a new phenomenon, nor on its main drivers. However, a lack of long-term vegetation datasets for the region and the scarcity of suitable palaeoecological archives, makes reconstructing past vegetation and land cover of the savannas a challenge.
To help meet this challenge, this study addresses three main research questions: 1) is pollen analysis a suitable tool to reflect the vegetation change associated with shrub encroachment in savanna environments? 2) Does the current encroached landscape correspond to an alternative stable state of savanna vegetation? 3) To what extent do pollen-based quantitative vegetation reconstructions reflect changes in past land cover?
The research focuses on north-central Namibia, where despite being the region most affected by shrub invasion, particularly since the 21st century, little is known about the dynamics of this phenomenon.
Field-based vegetation data were compared with modern pollen data to assess their correspondence in terms of composition and diversity along precipitation and grazing intensity gradients. In addition, two sediment cores from Lake Otjikoto were analysed to reveal changes in vegetation composition that have occurred in the region over the past 170 years and their possible drivers. For this, a multiproxy approach (fossil pollen, sedimentary ancient DNA (sedaDNA), biomarkers, compound specific carbon (δ13C) and deuterium (δD) isotopes, bulk carbon isotopes (δ13Corg), grain size, geochemical properties) was applied at high taxonomic and temporal resolution. REVEALS modelling of the fossil pollen record from Lake Otjikoto was run to quantitatively reconstruct past vegetation cover. For this, we first made pollen productivity estimates (PPE) of the most relevant savanna taxa in the region using the extended R-value model and two pollen dispersal options (Gaussian plume model and Lagrangian stochastic model). The REVEALS-based vegetation reconstruction was then validated using remote sensing-based regional vegetation data.
The results show that modern pollen reflects the composition of the vegetation well, but diversity less well. Interestingly, precipitation and grazing explain a significant amount of the compositional change in the pollen and vegetation spectra. The multiproxy record shows that a state change from open Combretum woodland to encroached Terminalia shrubland can occur over a century, and that the transition between states spans around 80 years and is characterized by a unique vegetation composition. This transition is supported by gradual environmental changes induced by management (i.e. broad-scale logging for the mining industry, selective grazing and reduced fire activity associated with intensified farming) and related land-use change. Derived environmental changes (i.e. reduced soil moisture, reduced grass cover, changes in species composition and competitiveness, reduced fire intensity) may have affected the resilience of Combretum open woodlands, making them more susceptible to change to an encroached state by stochastic events such as consecutive years of precipitation and drought, and by high concentrations of pCO2. We assume that the resulting encroached state was further stabilized by feedback mechanisms that favour the establishment and competitiveness of woody vegetation.
The REVEALS-based quantitative estimates of plant taxa indicate the predominance of a semi-open landscape throughout the 20th century and a reduction in grass cover below 50% since the 21st century associated with the spread of encroacher woody taxa. Cover estimates show a close match with regional vegetation data, providing support for the vegetation dynamics inferred from multiproxy analyses. Reasonable PPEs were made for all woody taxa, but not for Poaceae.
In conclusion, pollen analysis is a suitable tool to reconstruct past vegetation dynamics in savannas. However, because pollen cannot identify grasses beyond family level, a multiproxy approach, particularly the use of sedaDNA, is required. I was able to separate stable encroached states from mere woodland phases, and could identify drivers and speculate about related feedbacks. In addition, the REVEALS-based quantitative vegetation reconstruction clearly reflects the magnitude of the changes in the vegetation cover that occurred during the last 130 years, despite the limitations of some PPEs.
This research provides new insights into pollen-vegetation relationships in savannas and highlights the importance of multiproxy approaches when reconstructing past vegetation dynamics in semi-arid environments. It also provides the first time series with sufficient taxonomic resolution to show changes in vegetation composition during shrub encroachment, as well as the first quantitative reconstruction of past land cover in the region. These results help to identify the different stages in savanna dynamics and can be used to calibrate predictive models of vegetation change, which are highly relevant to land management.
In soils and sediments there is a strong coupling between local biogeochemical processes and the distribution of water, electron acceptors, acids and nutrients. Both sides are closely related and affect each other from small scale to larger scales. Soil structures such as aggregates, roots, layers or macropores enhance the patchiness of these distributions. At the same time it is difficult to access the spatial distribution and temporal dynamics of these parameter. Noninvasive imaging techniques with high spatial and temporal resolution overcome these limitations. And new non-invasive techniques are needed to study the dynamic interaction of plant roots with the surrounding soil, but also the complex physical and chemical processes in structured soils. In this study we developed an efficient non-destructive in-situ method to determine biogeochemical parameters relevant to plant roots growing in soil. This is a quantitative fluorescence imaging method suitable for visualizing the spatial and temporal pH changes around roots. We adapted the fluorescence imaging set-up and coupled it with neutron radiography to study simultaneously root growth, oxygen depletion by respiration activity and root water uptake. The combined set up was subsequently applied to a structured soil system to map the patchy structure of oxic and anoxic zones induced by a chemical oxygen consumption reaction for spatially varying water contents. Moreover, results from a similar fluorescence imaging technique for nitrate detection were complemented by a numerical modeling study where we used imaging data, aiming to simulate biodegradation under anaerobic, nitrate reducing conditions.
Stream restoration aims at an enhancement of ecological habitats, an increase of water retention within a landscape and sometimes even at an improvement of biogeochemical functions of lotic ecosystems. For the latter, good exchange between groundwater and stream water is often considered to be of major importance. In this study hydraulic connectivity between river and aquifer was investigated for a four years period, covering the restoration of an old oxbow after the second year. The oxbow became reconnected to the stream and the clogging layer in the oxbow was excavated. We expected increasing hydraulic connectivity between oxbow and aquifer after restoration of the stream, and decreasing hydraulic connectivity for the former shortcut due to increased clogging. To test that hypothesis, the spatial and temporal characteristics of the coupled groundwater-stream water system before and after the restoration were analysed by principal component analyses of time series of groundwater heads and stream water levels. The first component depicted between 53% and 70% of the total variance in the dataset for the different years. It captured the propagation of the pressure signal induced by stream water level fluctuations throughout the adjacent aquifer. Thus it could be used as a measure of hydraulic connectivity between stream and aquifer. During the first year, the impact of stream water level fluctuations decreased with distance from the regulated river (shortcut), whereas the hydraulic connection of the oxbow to the adjacent aquifer was very low. After restoration of the stream we observed a slight but not significant increase of hydraulic connectivity in the oxbow in the second year after restoration, but no change for the former shortcut. There is some evidence that the pattern of hydraulic connectivity at the study site is by far more determined by the natural heterogeneity of hydraulic conductivities of the floodplain sediments and the initial construction of the shortcut rather than by the clogging layer in the oxbow. (C) 2015 The Authors. Published by Elsevier B.V.
A novel approach for estimating precipitation patterns is developed here and applied to generate a new hydrologically corrected daily precipitation dataset, called RAIN4PE (Rain for Peru and Ecuador), at 0.1 degrees spatial resolution for the period 1981-2015 covering Peru and Ecuador. It is based on the application of 1) the random forest method to merge multisource precipitation estimates (gauge, satellite, and reanalysis) with terrain elevation, and 2) observed and modeled streamflow data to first detect biases and second further adjust gridded precipitation by inversely applying the simulated results of the ecohydrological model SWAT (Soil and Water Assessment Tool). Hydrological results using RAIN4PE as input for the Peruvian and Ecuadorian catchments were compared against the ones when feeding other uncorrected (CHIRP and ERA5) and gauge-corrected (CHIRPS, MSWEP, and PISCO) precipitation datasets into the model. For that, SWAT was calibrated and validated at 72 river sections for each dataset using a range of performance metrics, including hydrograph goodness of fit and flow duration curve signatures. Results showed that gauge-corrected precipitation datasets outperformed uncorrected ones for streamflow simulation. However, CHIRPS, MSWEP, and PISCO showed limitations for streamflow simulation in several catchments draining into the Pacific Ocean and the Amazon River. RAIN4PE provided the best overall performance for streamflow simulation, including flow variability (low, high, and peak flows) and water budget closure. The overall good performance of RAIN4PE as input for hydrological modeling provides a valuable criterion of its applicability for robust countrywide hydrometeorological applications, including hydroclimatic extremes such as droughts and floods. Significance StatementWe developed a novel precipitation dataset RAIN4PE for Peru and Ecuador by merging multisource precipitation data (satellite, reanalysis, and ground-based precipitation) with terrain elevation using the random forest method. Furthermore, RAIN4PE was hydrologically corrected using streamflow data in watersheds with precipitation underestimation through reverse hydrology. The results of a comprehensive hydrological evaluation showed that RAIN4PE outperformed state-of-the-art precipitation datasets such as CHIRP, ERA5, CHIRPS, MSWEP, and PISCO in terms of daily and monthly streamflow simulations, including extremely low and high flows in almost all Peruvian and Ecuadorian catchments. This underlines the suitability of RAIN4PE for hydrometeorological applications in this region. Furthermore, our approach for the generation of RAIN4PE can be used in other data-scarce regions.
Nitraria is a halophytic taxon (i.e., adapted to saline environments) that belongs to the plant family Nitrariaceae and is distributed from the Mediterranean, across Asia into the south-eastern tip of Australia. This taxon is thought to have originated in Asia during the Paleogene (66-23 Ma), alongside the proto-Paratethys epicontinental sea. The evolutionary history of Nitraria might hold important clues on the links between climatic and biotic evolution but limited taxonomic documentation of this taxon has thus far hindered this line of research. Here we investigate if the pollen morphology and the chemical composition of the pollen wall are informative of the evolutionary history of Nitraria and could explain if origination along the proto-Paratethys and dispersal to the Tibetan Plateau was simultaneous or a secondary process. To answer these questions, we applied a novel approach consisting of a combination of Fourier Transform Infrared spectroscopy (FTIR), to determine the chemical composition of the pollen wall, and pollen morphological analyses using Light Microscopy (LM) and Scanning Electron Microscopy (SEM). We analysed our data using ordinations (principal components analysis and non-metric multidimensional scaling), and directly mapped it on the Nitrariaceae phylogeny to produce a phylomorphospace and a phylochemospace. Our LM, SEM and FTIR analyses show clear morphological and chemical differences between the sister groups Peganum and Nitraria. Differences in the morphological and chemical characteristics of highland species (Nitraria schoberi, N. sphaerocarpa, N. sibirica and N. tangutorum) and lowland species (Nitraria billardierei and N. retusa) are very subtle, with phylogenetic history appearing to be a more important control on Nitraria pollen than local environmental conditions. Our approach shows a compelling consistency between the chemical and morphological characteristics of the eight studied Nitrariaceae species, and these traits are in agreement with the phylogenetic tree. Taken together, this demonstrates how novel methods for studying fossil pollen can facilitate the evolutionary investigation of living and extinct taxa, and the environments they represent.
In the past, floods were basically managed by flood control mechanisms. The focus was set on the reduction of flood hazard. The potential consequences were of minor interest. Nowadays river flooding is increasingly seen from the risk perspective, including possible consequences. Moreover, the large-scale picture of flood risk became increasingly important for disaster management planning, national risk developments and the (re-) insurance industry. Therefore, it is widely accepted that risk-orientated flood management ap-proaches at the basin-scale are needed. However, large-scale flood risk assessment methods for areas of several 10,000 km² are still in early stages. Traditional flood risk assessments are performed reach wise, assuming constant probabilities for the entire reach or basin. This might be helpful on a local basis, but where large-scale patterns are important this approach is of limited use. Assuming a T-year flood (e.g. 100 years) for the entire river network is unrealistic and would lead to an overestimation of flood risk at the large scale. Due to the lack of damage data, additionally, the probability of peak discharge or rainfall is usually used as proxy for damage probability to derive flood risk. With a continuous and long term simulation of the entire flood risk chain, the spatial variability of probabilities could be consider and flood risk could be directly derived from damage data in a consistent way.
The objective of this study is the development and application of a full flood risk chain, appropriate for the large scale and based on long term and continuous simulation. The novel approach of ‘derived flood risk based on continuous simulations’ is introduced, where the synthetic discharge time series is used as input into flood impact models and flood risk is directly derived from the resulting synthetic damage time series.
The bottleneck at this scale is the hydrodynamic simu-lation. To find suitable hydrodynamic approaches for the large-scale a benchmark study with simplified 2D hydrodynamic models was performed. A raster-based approach with inertia formulation and a relatively high resolution of 100 m in combination with a fast 1D channel routing model was chosen.
To investigate the suitability of the continuous simulation of a full flood risk chain for the large scale, all model parts were integrated into a new framework, the Regional Flood Model (RFM). RFM consists of the hydrological model SWIM, a 1D hydrodynamic river network model, a 2D raster based inundation model and the flood loss model FELMOps+r. Subsequently, the model chain was applied to the Elbe catchment, one of the largest catchments in Germany. For the proof-of-concept, a continuous simulation was per-formed for the period of 1990-2003. Results were evaluated / validated as far as possible with available observed data in this period. Although each model part introduced its own uncertainties, results and runtime were generally found to be adequate for the purpose of continuous simulation at the large catchment scale.
Finally, RFM was applied to a meso-scale catchment in the east of Germany to firstly perform a flood risk assessment with the novel approach of ‘derived flood risk assessment based on continuous simulations’. Therefore, RFM was driven by long term synthetic meteorological input data generated by a weather generator. Thereby, a virtual time series of climate data of 100 x 100 years was generated and served as input to RFM providing subsequent 100 x 100 years of spatially consistent river discharge series, inundation patterns and damage values. On this basis, flood risk curves and expected annual damage could be derived directly from damage data, providing a large-scale picture of flood risk. In contrast to traditional flood risk analysis, where homogenous return periods are assumed for the entire basin, the presented approach provides a coherent large-scale picture of flood risk. The spatial variability of occurrence probability is respected. Additionally, data and methods are consistent. Catchment and floodplain processes are repre-sented in a holistic way. Antecedent catchment conditions are implicitly taken into account, as well as physical processes like storage effects, flood attenuation or channel–floodplain interactions and related damage influencing effects. Finally, the simulation of a virtual period of 100 x 100 years and consequently large data set on flood loss events enabled the calculation of flood risk directly from damage distributions. Problems associated with the transfer of probabilities in rainfall or peak runoff to probabilities in damage, as often used in traditional approaches, are bypassed.
RFM and the ‘derived flood risk approach based on continuous simulations’ has the potential to provide flood risk statements for national planning, re-insurance aspects or other questions where spatially consistent, large-scale assessments are required.
A new view of Ecuador's complex geodynamics has been developed in the course of modeling seismic source zones for probabilistic seismic hazard analysis. This study focuses on two aspects of the plates' interaction at a continental scale: (a) age-related differences in rheology between Farallon and Nazca plates—marked by the Grijalva rifted margin and its inland projection—as they subduct underneath central Ecuador, and (b) the rapidly changing convergence obliquity resulting from the convex shape of the South American northwestern continental margin. Both conditions satisfactorily explain several characteristics of the observed seismicity and of the interseismic coupling. Intermediate-depth seismicity reveals a severe flexure in the Farallon slab as it dips and contorts at depth, originating the El Puyo seismic cluster. The two slabs position and geometry below continental Ecuador also correlate with surface expressions observable in the local and regional geology and tectonics. The interseismic coupling is weak and shallow south of the Grijalva rifted margin and increases northward, with a heterogeneous pattern locally associated to the Carnegie ridge subduction. High convergence obliquity is responsible for the North Andean Block northeastward movement along localized fault systems. The Cosanga and Pallatanga fault segments of the North Andean Block-South American boundary concentrate most of the seismic moment release in continental Ecuador. Other inner block faults located along the western border of the inter-Andean Depression also show a high rate of moderate-size earthquake production. Finally, a total of 19 seismic source zones were modeled in accordance with the proposed geodynamic and neotectonic scheme.
The age models of fluvio-lacustrine sedimentary sequences are often subject of discussions in paleoclimate research. The techniques employed to build an age model are very diverse, ranging from visual or intuitive estimation of the age-depth relationship over linear or spline interpolations between age control points to sophisticated Bayesian techniques also taking into account the most likely deposition times of the type of sediment within the sequence. All these methods, however, fail in detecting abrupt variations in sedimentation rates, including the possibility of episodes of no deposition (hiatus), which is the strength of the method presented in this work. The new technique simply compares the deposition time of equally thick sediment slices from the differences of subsequent radiometric age dates and the unit deposition times of the various sediment types. The percentage overlap of the distributions of these two sources of information, together with the evidence from the sedimentary record, helps to build an age model of complex sequences including abrupt variations in the rate of deposition including one or many hiatuses. (C) 2014 Elsevier B.V. All rights reserved.
Van Allen Probes measurements revealed the presence of the most unusual structures in the ultra-relativistic radiation belts. Detailed modeling, analysis of pitch angle distributions, analysis of the difference between relativistic and ultra-realistic electron evolution, along with theoretical studies of the scattering and wave growth, all indicate that electromagnetic ion cyclotron (EMIC) waves can produce a very efficient loss of the ultra-relativistic electrons in the heart of the radiation belts. Moreover, a detailed analysis of the profiles of phase space densities provides direct evidence for localized loss by EMIC waves. The evolution of multi-MeV fluxes shows dramatic and very sudden enhancements of electrons for selected storms. Analysis of phase space density profiles reveals that growing peaks at different values of the first invariant are formed at approximately the same radial distance from the Earth and show the sequential formation of the peaks from lower to higher energies, indicating that local energy diffusion is the dominant source of the acceleration from MeV to multi-MeV energies. Further simultaneous analysis of the background density and ultra-relativistic electron fluxes shows that the acceleration to multi-MeV energies only occurs when plasma density is significantly depleted outside of the plasmasphere, which is consistent with the modeling of acceleration due to chorus waves.
Sulfur is an important component in volcanic gases at the Earth surface but also present in the deep Earth in hydrothermal or magmatic fluids. Little is known about the evolution of such fluids during ascent in the crust. A new optical cell was developed for in situ Raman spectroscopic investigations on fluids allowing abrupt or continuous changes of pressure up to 200 MPa at temperatures up to 750 degrees C. The concept is based on a flexible gold bellow, which separates the sample fluid from the pressure medium water. To avoid reactions between aggressive fluids and the pressure cell, steel components in contact with the fluid are shielded by gold foil. The cell was tested to study redox reactions in fluids using aqueous ammonium sulfate solutions as a model system. During heating at constant pressure of 130 MPa, sulfate ions transform first to HSO4- ions and then to molecular units such as H2SO4. Variation of pressure shows that the stability of sulfate species relies on fluid density, i.e., highly charged species are stable only in high-density fluids. Partial decomposition of ammonium was evident above 550 degrees C by the occurrence of a nitrogen peak in the Raman spectra. Reduced sulfur species were observed above 700 degrees C by Raman signals near 2590 cm(-1) assigned to HS- and H2S. No clear evidence for the formation of sulfur dioxide was found in contrary to previous studies on aqueous H2SO4, suggesting very reducing conditions in our experiments. Fluid-mineral interaction was studied by inserting into the cell a small, semi-open capsule filled with a mixture of pyrite and pyrrhotite. Oxidation of the sample assembly was evident by transformation of pyrite to pyrrhotite. As a consequence, sulfide species were observed in the fluid already at temperatures of similar to 600 degrees C.
Quantitative estimates of sea-level rise in the Mediterranean Basin become increasingly accurate thanks to detailed satellite monitoring. However, such measuring campaigns cover several years to decades, while longer-term sea-level records are rare for the Mediterranean. We used a data archeological approach to reanalyze monthly mean sea-level data of the Antalya-I (1935–1977) tide gauge to fill this gap. We checked the accuracy and reliability of these data before merging them with the more recent records of the Antalya-II (1985–2009) tide gauge, accounting for an eight-year hiatus. We obtain a composite time series of monthly and annual mean sea levels spanning some 75 years, providing the longest record for the eastern Mediterranean Basin, and thus an essential tool for studying the region's recent sea-level trends. We estimate a relative mean sea-level rise of 2.2 ± 0.5 mm/year between 1935 and 2008, with an annual variability (expressed here as the standard deviation of the residuals, σresiduals = 41.4 mm) above that at the closest tide gauges (e.g., Thessaloniki, Greece, σresiduals = 29.0 mm). Relative sea-level rise accelerated to 6.0 ± 1.5 mm/year at Antalya-II; we attribute roughly half of this rate (~3.6 mm/year) to tectonic crustal motion and anthropogenic land subsidence. Our study highlights the value of data archeology for recovering and integrating historic tide gauge data for long-term sea-level and climate studies.
With increasing amount of strong motion data, Ground Motion Prediction Equation (GMPE) developers are able to quantify empirical site amplification functions (delta S2S(s)) from GMPE residuals, for use in site-specific Probabilistic Seismic Hazard Assessment. In this study, we first derive a GMPE for 5% damped Pseudo Spectral Acceleration (g) of Active Shallow Crustal earthquakes in Japan with 3.4 <= M-w <= 7.3 and 0 <= R-JB <= 600km. Using k-mean spectral clustering technique, we then classify our estimated delta S2S(s)(T = 0.01 - 2s) of 588 wellcharacterized sites, into 8 site clusters with distinct mean site amplification functions, and within-cluster site-tosite variability similar to 50% smaller than the overall dataset variability (phi(S2S)). Following an evaluation of existing schemes, we propose a revised data-driven site classification characterized by kernel density distributions of V-s30, V-s10, H-800, and predominant period (T-G) of the site clusters.
The complementary advantages of GPS and seismic measurements are well recognized in seismotectonic monitoring studies. Therefore, integrated processing of the two data streams has been proposed recently in an attempt to obtain accurate and reliable information of surface displacements associated with earthquakes. A hitherto still critical issue in the integrated processing is real-time detection and precise estimation of the transient baseline error in the seismic records. Here, we report on a new approach by introducing the seismic acceleration corrected by baseline errors into the state equation system. The correction is performed and regularly updated in short epochs (with increments which may be as short as seconds), so that station position, velocity, and acceleration can be constrained very tightly and baseline error can be estimated as a random-walk process. With the adapted state equation system, our study highlights the use of a new approach developed for integrated processing of GPS and seismic data by means of sequential least-squares adjustment. The efficiency of our approach is demonstrated and validated using simulated, experimental, and real datasets. The latter were collected at collocated GPS and seismic stations around the 4 April 2010, E1 Mayor-Cucapah earthquake (Mw, 7.2). The results have shown that baseline errors of the strong-motion sensors are corrected precisely and high-precision seismic displacements are real-timely obtained by the new approach.
Flow accumulation algorithms estimate the steady state of flow on real or modeled topographic surfaces and are crucial for hydrological and geomorphological assessments, including delineation of river networks, drainage basins, and sediment transport processes. Existing flow accumulation algorithms are typically designed to compute flows on regular grids and are not directly applicable to arbitrarily sampled topographic data such as lidar point clouds. In this study we present a random sampling scheme that generates homogeneous point densities, in combination with a novel flow path tracing approach-the Facet-Flow Network (FFN)-that estimates flow accumulation in terms of specific catchment area (SCA) on triangulated surfaces. The random sampling minimizes biases due to spatial sampling and the FFN allows for direct flow estimation from point clouds. We validate our approach on a Gaussian hill surface and study the convergence of its SCA compared to the analytical solution. Here, our algorithm outperforms the multiple flow direction algorithm, which is optimized for divergent surfaces. We also compute the SCA of a 6-km(2)-steep, vegetated catchment on Santa Cruz Island, California, based on airborne lidar point-cloud data. Point-cloud-based SCA values estimated by our method compare well with those estimated by the D-infinity or multiple flow direction algorithm on gridded data. The advantage of computing SCA from point clouds becomes relevant especially for divergent topography and for small drainage areas: These are depicted with much more detail due to the higher sampling density of point clouds.
The polar and subtropical jet streams are strong upper-level winds with a crucial influence on weather throughout the Northern Hemisphere midlatitudes. In particular, the polar jet is located between cold arctic air to the north and warmer subtropical air to the south. Strongly meandering states therefore often lead to extreme surface weather. The parameter values of the detection scheme are optimized using simulated annealing and a skill function that accounts for the zonal-mean jet stream position (Rikus, 2015). After the successful optimization process, we apply our scheme to reanalysis data covering 1979-2015 and calculate seasonal-mean probabilistic maps and trends in wind strength and position of jet streams. We present longitudinally defined probability distributions of the positions for both jets for all on the Northern Hemisphere seasons. This shows that winter is characterized by two well-separated jets over Europe and Asia (ca. 20 degrees W to 140 degrees E). In contrast, summer normally has a single merged jet over the western hemisphere but can have both merged and separated jet states in the eastern hemisphere.
Frequent and intense rainfall events demand innovative techniques to better predict the extreme rainfall dynamics. This task requires essentially the assessment of the basic types of atmospheric processes that trigger extreme rainfall, and then to examine the differences between those processes, which may help to identify key patterns to improve predictive algorithms. We employ tools from network theory to compare the spatial features of extreme rainfall over the Japanese archipelago and surrounding areas caused by two atmospheric processes: the Baiu front, which occurs mainly in June and July (JJ), and the tropical storms from August to November (ASON). We infer from complex networks of satellite-derived rainfall data, which are based on the nonlinear correlation measure of event synchronization. We compare the spatial scales involved in both systems and identify different regions which receive rainfall due to the large spatial scale of the Baiu and tropical storm systems. We observed that the spatial scales involved in the Baiu driven rainfall extremes, including the synoptic processes behind the frontal development, are larger than tropical storms, which even have long tracks during extratropical transitions. We further delineate regions of coherent rainfall during the two seasons based on network communities, identifying the horizontal (east-west) rainfall bands during JJ over the Japanese archipelago, while during ASON these bands align with the island arc of Japan.
The combined passive and active seismic TRANSALP experiment produced an unprecedented high-resolution crustal image of the Eastern Alps between Munich and Venice. The European and Adriatic Mohos (EM and AM, respectively) are clearly imaged with different seismic techniques: near-vertical incidence reflections and receiver functions (RFs). The European Moho dips gently southward from 35 km beneath the northern foreland to a maximum depth of 55 km beneath the central part of the Eastern Alps, whereas the Adriatic Moho is imaged primarily by receiver functions at a relatively constant depth of about 40 km. In both data sets, we have also detected first-order Alpine shear zones, such as the Helvetic detachment, Inntal fault and SubTauern ramp in the north. Apart from the Valsugana thrust, receiver functions in the southern part of the Eastern Alps have also observed a north dipping interface, which may penetrate the entire Adriatic crust [Adriatic Crust Interface (ACI)]. Deep crustal seismicity may be related to the ACI. We interpret the ACI as the currently active retroshear zone in the doubly vergent Alpine collisional belt. (C) 2004 Elsevier B.V. All rights reserved
Despite advanced seismological techniques, automatic source characterization for microseismic earthquakes remains difficult and challenging since current inversion and modelling of high-frequency signals are complex and time consuming. For real-time applications such as induced seismicity monitoring, the application of standard methods is often not fast enough for true complete real-time information on seismic sources. In this paper, we present an alternative approach based on recent advances in deep learning for rapid source-parameter estimation of microseismic earthquakes. The seismic inversion is represented in compact form by two convolutional neural networks, with individual feature extraction, and a fully connected neural network, for feature aggregation, to simultaneously obtain full moment tensor and spatial location of microseismic sources. Specifically, a multibranch neural network algorithm is trained to encapsulate the information about the relationship between seismic waveforms and underlying point-source mechanisms and locations. The learning-based model allows rapid inversion (within a fraction of second) once input data are available. A key advantage of the algorithm is that it can be trained using synthetic seismic data only, so it is directly applicable to scenarios where there are insufficient real data for training. Moreover, we find that the method is robust with respect to perturbations such as observational noise and data incompleteness (missing stations). We apply the new approach on synthesized and example recorded small magnitude (M <= 1.6) earthquakes at the Hellisheioi geothermal field in the Hengill area, Iceland. For the examined events, the model achieves excellent performance and shows very good agreement with the inverted solutions determined through standard methodology. In this study, we seek to demonstrate that this approach is viable for microseismicity real-time estimation of source parameters and can be integrated into advanced decision-support tools for controlling induced seismicity.
Cosmogenic nuclides are typically used to either constrain an exposure age, a burial age, or an erosion rate. Constraining the landscape history and past erosion rates in previously glaciated terrains is, however, notoriously difficult because it involves a large number of unknowns. The potential use of cosmogenic nuclides in landscapes with a complex history of exposure and erosion is therefore often quite limited. Here, we present a novel multi-nuclide approach to study the landscape evolution and past erosion rates in terrains with a complex exposure history, particularly focusing on regions that were repeatedly covered by glaciers or ice sheets during the Quaternary. The approach, based on the Markov Chain Monte Carlo (MCMC) technique, focuses on mapping the range of landscape histories that are consistent with a given set of measured cosmogenic nuclide concentrations. A fundamental assumption of the model approach is that the exposure history at the site/location can be divided into two distinct regimes: i) interglacial periods characterized by zero shielding due to overlying ice and a uniform interglacial erosion rate, and ii) glacial periods characterized by 100% shielding and a uniform glacial erosion rate. We incorporate the exposure history in the model framework by applying a threshold value to the global marine benthic delta O-18 record and include the threshold value as a free model parameter, hereby taking into account global changes in climate. However, any available information on the glacial-interglacial history at the sampling location, in particular the timing of the last deglaciation event, is readily incorporated in the model to constrain the inverse problem. Based on the MCMC technique, the model delineates the most likely exposure history, including the glacial and interglacial erosion rates, which, in turn, makes it possible to reconstruct an exhumation history at the site. We apply the model to two landscape scenarios based on synthetic data and two landscape scenarios based on paired Be-10/Al-26 data from West Greenland, which makes it possible to quantify the denudation rate at these locations. The model framework, which currently incorporates any combination of the following nuclides Be-10, Al-26, C-14, and Ne-21, is highly flexible and can be adapted to many different landscape settings. The model framework may also be used in combination with physics-based landscape evolution models to predict nuclide concentrations at different locations in the landscape. This may help validate the landscape models via comparison to measured nuclide concentrations or to devise new effective sampling strategies. (C) 2015 The Authors. Published by Elsevier B.V.
Background and Aims Dynamic processes occurring at the soil-root interface crucially influence soil physical, chemical and biological properties at a local scale around the roots, and are technically challenging to capture in situ. This study presents a novel multi-imaging approach combining fluorescence and neutron radiography that is able to simultaneously monitor root growth, water content distribution, root respiration and root exudation.
Methods Germinated seeds of white lupins (Lupinus albus) were planted in boron-free glass rhizotrons. After 11 d, the rhizotrons were wetted from the bottom and time series of fluorescence and neutron images were taken during the subsequent day and night cycles for 13 d. The following day (i.e. 25 d after planting) the rhizotrons were again wetted from the bottom and the measurements were repeated. Fluorescence sensor foils were attached to the inner sides of the glass and measurements of oxygen and pH were made on the basis of fluorescence intensity. The experimental set-up allowed for simultaneous fluorescence imaging and neutron radiography.
Key Results The interrelated patterns of root growth and distribution in the soil, root respiration, exudation and water uptake could all be studied non-destructively and at high temporal and spatial resolution. The older parts of the root system with greater root-length density were associated with fast decreases of water content and rapid changes in oxygen concentration. pH values around the roots located in areas with low soil water content were significantly lower than the rest of the root system.
Conclusions The results suggest that the combined imaging set-up developed here, incorporating fluorescence intensity measurements, is able to map important biogeochemical parameters in the soil around living plants with a spatial resolution that is sufficiently high enough to relate the patterns observed to the root system.
A modern pollen dataset from China and Mongolia (18-52 degrees N, 74-132 degrees E) is investigated for its potential use in climate reconstructions. The dataset includes 2559 samples, 229 terrestrial pollen taxa and four climatic variables - mean annual precipitation (P-ann): 35-2091 mm, mean annual temperature (T-ann): -12.1-25.8 degrees C, mean temperature in the coldest month (Mt(co).): -33.8-21.7 degrees C, and mean temperature in the warmest month (Mt(wa)): 03-29.8 degrees C. Modern pollen-climate relationships are assessed using canonical correspondence analysis (CCA), Huisman-Olff-Fresco (HOF) models, the modern analogue technique (MAT), and weighted averaging partial least squares (WA-PLS). Results indicate that P-ann is the most important climatic determinant of pollen distribution and the most promising climate variable for reconstructions, as assessed by the coefficient of determination between observed and predicted environmental values (r(2)) and root mean square error of prediction (RMSEP). Mt(co) and Mt(wa) may be reconstructed too, but with caution. Samples from different depositional environments influence the performance of cross-validation differently, with samples from lake sediment-surfaces and moss polsters having the best fit with the lowest RMSEP. The better model performances of MAT are most probably caused by spatial autocorrelation. Accordingly, the WA-PLS models of this dataset are deemed most suitable for reconstructing past climate quantitatively because of their more reliable predictive power. (C) 2014 Elsevier B.V. All rights reserved.
AimFossil pollen spectra from lake sediments in central and western Mongolia have been used to interpret past climatic variations, but hitherto no suitable modern pollen-climate calibration set has been available to infer past climate changes quantitatively. We established such a modern pollen dataset and used it to develop a transfer function model that we applied to a fossil pollen record in order to investigate: (1) whether there was a significant moisture response to the Younger Dryas event in north-western Mongolia; and (2) whether the early Holocene was characterized by dry or wet climatic conditions.
LocationCentral and western Mongolia.
MethodsWe analysed pollen data from surface sediments from 90 lakes. A transfer function for mean annual precipitation (P-ann) was developed with weighted averaging partial least squares regression (WA-PLS) and applied to a fossil pollen record from Lake Bayan Nuur (49.98 degrees N, 93.95 degrees E, 932m a.s.l.). Statistical approaches were used to investigate the modern pollen-climate relationships and assess model performance and reconstruction output.
ResultsRedundancy analysis shows that the modern pollen spectra are characteristic of their respective vegetation types and local climate. Spatial autocorrelation and significance tests of environmental variables show that the WA-PLS model for P-ann is the most valid function for our dataset, and possesses the lowest root mean squared error of prediction.
Main conclusionsPrecipitation is the most important predictor of pollen and vegetation distributions in our study area. Our quantitative climate reconstruction indicates a dry Younger Dryas, a relatively dry early Holocene, a wet mid-Holocene and a dry late Holocene.
Aim: Fossil pollen spectra from lake sediments on the Tibetan Plateau have been used for qualitative climate reconstruction, but no modern pollen-climate calibration set based on lake sediments is available to infer past climate quantitatively. This study aims to develop such a dataset and apply it to fossil data. Location: The Tibetan Plateau, between 30 and 40 degrees N and 87 and 103 degrees E. Methods: We collected surface sediments from 112 lakes and analysed them palynologically. The lakes span a wide range of mean annual precipitation (P-ann; 31-1022 mm), mean annual temperature (T-ann; -6.5 to 1 degrees C), and mean July temperature (T-July; 2.6-19.7 degrees C). Redundancy analysis showed that the modern pollen spectra are characteristic of their respective vegetation types and local climate. Transfer functions for P-ann, T-ann and T-July were developed with weighted averaging partial least squares. Model performance was assessed by leave-one-out cross-validation. Results: The root mean square errors of prediction (RMSEP) were 104 mm (P-ann), 1.18 degrees C (T-ann) and 1.17 degrees C (T-July). The RMSEPs, when expressed as percentages of the gradient sampled, were 10.6% (P-ann), 15.7% (T-ann) and 11.9% (T-July). These low values indicate the good performance of our models. An application of the models to fossil pollen spectra covering the last c. 50 kyr yielded realistic results for Luanhaizi Lake in the Qilian Mountains on the north-eastern Tibetan Plateau (modern P-ann 480 mm; T-ann-1 degrees C). T-ann and P-ann values similar to present ones were reconstructed for late Marine Isotope Stage 3, with minimum values for the Last Glacial Maximum (c. 300 mm and 2 degrees C below present), and maximum values for the early Holocene (c. 70 mm and 0.5 degrees C greater than present). Main conclusions: The modern pollen-climate calibration set will potentially be useful for quantitative climate reconstructions from lake-sediment pollen spectra from the Tibetan Plateau, an area of considerable climatic and biogeographical importance.
Atmospheric water vapour content is a key variable that controls the development of deep convective storms and rainfall extremes over the central Andes. Direct measurements of water vapour are challenging; however, recent developments in microwave processing allow the use of phase delays from L-band radar to measure the water vapour content throughout the atmosphere: Global Navigation Satellite System (GNSS)-based integrated water vapour (IWV) monitoring shows promising results to measure vertically integrated water vapour at high temporal resolutions. Previous works also identified convective available potential energy (CAPE) as a key climatic variable for the formation of deep convective storms and rainfall in the central Andes. Our analysis relies on GNSS data from the Argentine Continuous Satellite Monitoring Network, Red Argentina de Monitoreo Satelital Continuo (RAMSAC) network from 1999 to 2013. CAPE is derived from version 2.0 of the ECMWF’s (European Centre for Medium-Range Weather Forecasts) Re-Analysis (ERA-interim) and rainfall from the TRMM (Tropical Rainfall Measuring Mission) product. In this study, we first analyse the rainfall characteristics of two GNSS-IWV stations by comparing their complementary cumulative distribution function (CCDF). Second, we separately derive the relation between rainfall vs. CAPE and GNSS-IWV. Based on our distribution fitting analysis, we observe an exponential relation of rainfall to GNSS-IWV. In contrast, we report a power-law relationship between the daily mean value of rainfall and CAPE at the GNSS-IWV station locations in the eastern central Andes that is close to the theoretical relationship based on parcel theory. Third, we generate a joint regression model through a multivariable regression analysis using CAPE and GNSS-IWV to explain the contribution of both variables in the presence of each other to extreme rainfall during the austral summer season. We found that rainfall can be characterised with a higher statistical significance for higher rainfall quantiles, e.g., the 0.9 quantile based on goodness-of-fit criterion for quantile regression. We observed different contributions of CAPE and GNSS-IWV to rainfall for each station for the 0.9 quantile. Fourth, we identify the temporal relation between extreme rainfall (the 90th, 95th, and 99th percentiles) and both GNSS-IWV and CAPE at 6 h time steps. We observed an increase before the rainfall event and at the time of peak rainfall—both for GNSS-integrated water vapour and CAPE. We show higher values of CAPE and GNSS-IWV for higher rainfall percentiles (99th and 95th percentiles) compared to the 90th percentile at a 6-h temporal scale. Based on our correlation analyses and the dynamics of the time series, we show that both GNSS-IWV and CAPE had comparable magnitudes, and we argue to consider both climatic variables when investigating their effect on rainfall extremes.
Atmospheric water vapour content is a key variable that controls the development of deep convective storms and rainfall extremes over the central Andes. Direct measurements of water vapour are challenging; however, recent developments in microwave processing allow the use of phase delays from L-band radar to measure the water vapour content throughout the atmosphere: Global Navigation Satellite System (GNSS)-based integrated water vapour (IWV) monitoring shows promising results to measure vertically integrated water vapour at high temporal resolutions. Previous works also identified convective available potential energy (CAPE) as a key climatic variable for the formation of deep convective storms and rainfall in the central Andes. Our analysis relies on GNSS data from the Argentine Continuous Satellite Monitoring Network, Red Argentina de Monitoreo Satelital Continuo (RAMSAC) network from 1999 to 2013. CAPE is derived from version 2.0 of the ECMWF’s (European Centre for Medium-Range Weather Forecasts) Re-Analysis (ERA-interim) and rainfall from the TRMM (Tropical Rainfall Measuring Mission) product. In this study, we first analyse the rainfall characteristics of two GNSS-IWV stations by comparing their complementary cumulative distribution function (CCDF). Second, we separately derive the relation between rainfall vs. CAPE and GNSS-IWV. Based on our distribution fitting analysis, we observe an exponential relation of rainfall to GNSS-IWV. In contrast, we report a power-law relationship between the daily mean value of rainfall and CAPE at the GNSS-IWV station locations in the eastern central Andes that is close to the theoretical relationship based on parcel theory. Third, we generate a joint regression model through a multivariable regression analysis using CAPE and GNSS-IWV to explain the contribution of both variables in the presence of each other to extreme rainfall during the austral summer season. We found that rainfall can be characterised with a higher statistical significance for higher rainfall quantiles, e.g., the 0.9 quantile based on goodness-of-fit criterion for quantile regression. We observed different contributions of CAPE and GNSS-IWV to rainfall for each station for the 0.9 quantile. Fourth, we identify the temporal relation between extreme rainfall (the 90th, 95th, and 99th percentiles) and both GNSS-IWV and CAPE at 6 h time steps. We observed an increase before the rainfall event and at the time of peak rainfall—both for GNSS-integrated water vapour and CAPE. We show higher values of CAPE and GNSS-IWV for higher rainfall percentiles (99th and 95th percentiles) compared to the 90th percentile at a 6-h temporal scale. Based on our correlation analyses and the dynamics of the time series, we show that both GNSS-IWV and CAPE had comparable magnitudes, and we argue to consider both climatic variables when investigating their effect on rainfall extremes.
The features of Middle Miocene deposits in the Puna-Eastern Cordillera transition (Valles Calchaquies) indicate that Cenozoic deformation, sedimentation and volcanism follow a complex spatiotemporal relationship. The intense volcanic activity recorded in the eastern Puna border between 14 and 11.5 Ma coincides with the occurrence of one of the most important deformation events of the Neogene tectonic evolution in the region. Studies performed across the Puna-Eastern Cordillera transition show different relationships between volcanic deposits of ca. 13.5-12.1 Ma and the Oligocene-Miocene Angastaco Formation. In this paper we describe the ash-flow tuff deposits which are the first of this type found concordant in the sedimentary fill of Valles Calchaquies. Several analyses performed on these pyroclastic deposits allow a correlation to be made with the Alto de Las Lagunas Ignimbrite (ca. 13.5 Ma) of the Pucarilla-Cerro Tipillas Volcanic Complex located in the Puna. Outcrops of the ca. 13.5 Ma pyroclastic deposits are recognised within the Puna and the Valle Calchaqui. However, in the southern prolongation of the Valle de Hualfin (Tiopampa-Pucarilla depression) that separates the Puna from the Valle Calchaqui at these latitudes, these deposits are partially eroded and buried, and thus their occurrence is recorded only by abundant volcanic clasts included in conglomerates of the Angastaco Formation. The sedimentation of the Angastaco Formation was aborted at ca. 12 Ma in the Tiopampa-Pucarilla depression by the Pucarilla Ignimbrite, which unconformably covers the synorogenic units. On the contrary, in the Valle Calchaqui the sedimentation of the Angastaco Formation continued until the Late Miocene. The different relationships between the Miocene Angastaco Formation and the ignimbrites with ages of ca. 13.5 and ca. 12 Ma reveal that in this short period (-1.5 m.y.) a significant deformation event took place and resulted in marked palaeogeographic changes, as evidenced by stratigraphic-sedimentological and chronological records in the Angastaco Formation. (C) 2017 Elsevier B.V. All rights reserved.
This manuscript proposes a method to assess hydrological drought in semi-arid environments under high impoundment rate and applies it to the semi-arid Jaguaribe River basin in Brazil. It analyzes droughts (1) in the largest reservoir systems; (2) in the Upper Basin, considering 4744 reservoirs, 800 wells and almost 18,000 cisterns; and (3) in reservoirs of different sizes during multiyear droughts. Results show that the water demand is constrained in the basin; hydrological and meteorological droughts are often out of phase; there is a negative correlation between storage level and drought severity; and the small systems cannot cope with long-term droughts.
The Pleistocene archeological record in East Africa has revealed unusual accumulations of Acheulean handaxes at prehistoric sites. In particular, there has been intensive debate concerning whether the artifact accumulation at the Middle Pleistocene Olorgesailie (Southern Kenya Rift) and Kariandusi (Central Kenya Rift) sites were a result of fluvial reworking or of in situ deposition by hominids. We used a two-step approach to test the hypothesis of fluvial reworking. Firstly, the behavior of handaxes in water currents was investigated in a current flume and the flow threshold required to reorientate the handaxes was determined. The results of these experiments suggested that, in relatively high energy and non-steady flow conditions, handaxes will reorientate themselves perpendicular to the current direction. Secondly, an automated image analysis routine was developed and applied to archeological plans from three Acheulean sites, two at Olorgesailie and one at Kariandusi, in order to determine the orientations of the handaxes. A Rayleigh test was then applied to the orientation data to test for a preferred orientation. The results revealed that the handaxes at the Upper Kariandusi Site and the Olorgesailie Main Site Mid Trench had a preferential orientation, suggesting reworking by a paleocurrent. The handaxes from the Olorgesailie Main Site H/6A, however, appeared to be randomly oriented and in situ deposition by the producers therefore remains a possibility.
We apply and evaluate a recent machine learning method for the automatic classification of seismic waveforms. The method relies on Dynamic Bayesian Networks (DBN) and supervised learning to improve the detection capabilities at 3C seismic stations. A time-frequency decomposition provides the basis for the required signal characteristics we need in order to derive the features defining typical "signal" and "noise" patterns. Each pattern class is modeled by a DBN, specifying the interrelationships of the derived features in the time-frequency plane. Subsequently, the models are trained using previously labeled segments of seismic data. The DBN models can now be compared against in order to determine the likelihood of new incoming seismic waveform segments to be either signal or noise. As the noise characteristics of seismic stations varies smoothly in time (seasonal variation as well as anthropogenic influence), we accommodate in our approach for a continuous adaptation of the DBN model that is associated with the noise class. Given the difficulty for obtaining a golden standard for real data (ground truth) the proof of concept and evaluation is shown by conducting experiments based on 3C seismic data from the International Monitoring Stations, BOSA and LPAZ.
Ancient evaporite deposits are geological archives of depositional environments characterized by a long‐term negative precipitation balance and bear evidence for global ocean element mass balance calculations. Here, Cretaceous selenite pseudomorphs from western Anatolia (‘Rosetta Marble’) — characterized by their exceptional morphological preservation — and their ‘marine’ geochemical signatures are described and interpreted in a process‐oriented context. These rocks recorded Late Cretaceous high‐pressure/low‐temperature, subduction‐related metamorphism with peak conditions of 1·0 to 1·2 GPa and 300 to 400°C. Metre‐scale, rock‐forming radiating rods, now present as fibrous calcite marble, clearly point to selenitic gypsum as the precursor mineral. Stratigraphic successions are recorded along a reconstructed proximal to distal transect. The cyclical alternation of selenite beds and radiolarian ribbon‐bedded cherts in the distal portions are interpreted as a two type of seawater system. During arid intervals, shallow marine brines cascaded downward into basinal settings and induced precipitation. During more humid times, upwelling‐induced radiolarian blooms caused the deposition of radiolarite facies. Interestingly, there is no comparable depositional setting known from the Cenozoic world. Meta‐selenite geochemical data (δ13C, δ18O and 87Sr/86Sr) plot within the range of reconstructed middle Cretaceous seawater signatures. Possible sources for the 13C‐enriched (mean 2·2‰) values include methanogenesis, gas hydrates and cold seep fluid exhalation. Spatially resolved component‐specific analysis of a rock slab displays isotopic variances between meta‐selenite crystals (mean δ13C 2·2‰) and host matrix (mean δ13C 1·3‰). The Cretaceous evaporite‐pseudomorphs of Anatolia represent a basin wide event coeval with the Aptian evaporites of the Proto‐Atlantic and the pseudomorphs share many attributes, including lateral distribution of 600 km and stratigraphic thickness of 1·5 to 2·0 km, with the evaporites formed during the younger Messinian salinity crisis. The Rosetta Marble of Anatolia may represent the best‐preserved selenite pseudomorphs worldwide and have a clear potential to act as a template for the study of meta‐selenite in deep time.
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.
Computer-based simulation models are frequently used in hydrological research and engineering but also in other fields of environmental sciences. New case studies often require existing model concepts to be adapted. Extensions may be necessary due to the peculiarities of the studied natural system or subtleties of anthropogenic control. In other cases, simplifications must be made in response to scarce data, incomplete knowledge, or restrictions set by the spatio-temporal scale of application. This paper introduces an open-source modeling framework called ECHSE designed to cope with the above-mentioned challenges. It provides a lightweight infrastructure for the rapid development of new, reusable simulation tools and, more importantly, the safe modification of existing formulations. ECHSE-based models treat the simulated system as a collection of interacting objects. Although feedbacks are generally supported, the majority of the objects' interactions is expected to be of the feed-forward type. Therefore, the ECHSE software is particularly useful in the context of hydrological catchment modeling. Conversely, it is unsuitable, e.g., for fully hydrodynamic simulations and groundwater flow modeling. The focus of the paper is put on a comprehensible outline of the ECHSE's fundamental concepts and limitations. For the purpose of illustration, a specific, ECHSE-based solution for hydrological catchment modeling is presented which has undergone testing in a number of river basins. (C) 2015 Elsevier Ltd. All rights reserved.
A total of 271 pollen records were selected from a large collection of both raw and digitized pollen spectra from eastern continental Asia (70 degrees-135 degrees E and 18 degrees-55 degrees N). Following pollen percentage recalculations, taxonomic homogenization, and age-depth model revision, the pollen spectra were interpolated at a 500-year resolution and a taxonomically harmonized and temporally standardized fossil pollen dataset established with 226 pollen taxa, covering the last 22 cal lea. Of the 271 pollen records, 85% were published since 1990, with reliable chronologies and high temporal resolutions; of these, 50% have raw data with complete pollen assemblages, ensuring the quality of this dataset The pollen records available for each 500-year time slice are well distributed over all main vegetation types and climatic zones of the study area, making their pollen spectra suitable for paleovegetation and paleoclimate research. Such a dataset can be used as an example for the development of similar datasets for other regions of the world.
A large landslide (frozen debris avalanche) occurred at Assapaat on the south coast of the Nuussuaq Peninsula in Central West Greenland on June 13, 2021, at 04:04 local time. We present a compilation of available data from field observations, photos, remote sensing, and seismic monitoring to describe the event. Analysis of these data in combination with an analysis of pre- and post-failure digital elevation models results in the first description of this type of landslide. The frozen debris avalanche initiated as a 6.9 * 10(6) m(3) failure of permafrozen talus slope and underlying colluvium and till at 600-880 m elevation. It entrained a large volume of permafrozen colluvium along its 2.4 km path in two subsequent entrainment phases accumulating a total volume between 18.3 * 10(6) and 25.9 * 10(6) m(3). About 3.9 * 10(6) m(3) is estimated to have entered the Vaigat strait; however, no tsunami was reported, or is evident in the field. This is probably because the second stage of entrainment along with a flattening of slope angle reduced the mobility of the frozen debris avalanche. We hypothesise that the initial talus slope failure is dynamically conditioned by warming of the ice matrix that binds the permafrozen talus slope. When the slope ice temperature rises to a critical level, its shear resistance is reduced, resulting in an unstable talus slope prone to failure. Likewise, we attribute the large-scale entrainment to increasing slope temperature and take the frozen debris avalanche as a strong sign that the permafrost in this region is increasingly at a critical state. Global warming is enhanced in the Arctic and frequent landslide events in the past decade in Western Greenland let us hypothesise that continued warming will lead to an increase in the frequency and magnitude of these types of landslides. Essential data for critical arctic slopes such as precipitation, snowmelt, and ground and surface temperature are still missing to further test this hypothesis. It is thus strongly required that research funds are made available to better predict the change of landslide threat in the Arctic.
Lakes are dominant and diverse landscape features in the Arctic, but conventional land cover classification schemes typically map them as a single uniform class. Here, we present a detailed lake-centric geospatial database for an Arctic watershed in northern Alaska. We developed a GIS dataset consisting of 4362 lakes that provides information on lake morphometry, hydrologic connectivity, surface area dynamics, surrounding terrestrial ecotypes, and other important conditions describing Arctic lakes. Analyzing the geospatial database relative to fish and bird survey data shows relations to lake depth and hydrologic connectivity, which are being used to guide research and aid in the management of aquatic resources in the National Petroleum Reserve in Alaska. Further development of similar geospatial databases is needed to better understand and plan for the impacts of ongoing climate and land-use changes occurring across lake-rich landscapes in the Arctic.
We present a new autoclave that enables in situ characterization of hydrothermal fluids at high pressures and high temperatures at synchrotron x-ray radiation sources. The autoclave has been specifically designed to enable x-ray absorption spectroscopy in fluids with applications to mineral solubility and element speciation analysis in hydrothermal fluids in complex compositions. However, other applications, such as Raman spectroscopy, in high-pressure fluids are also possible with the autoclave. First experiments were run at pressures between 100 and 600 bars and at temperatures between 25 degrees C and 550 degrees C, and preliminary results on scheelite dissolution in fluids of different compositions show that the autoclave is well suited to study the behavior of ore-forming metals at P-T conditions relevant to the Earth's crust.
A hydrochemical approach to quantify the role of return flow in a surface flow-dominated catchment
(2017)
Stormflow generation in headwater catchments dominated by subsurface flow has been studied extensively, yet catchments dominated by surface flow have received less attention. We addressed this by testing whether stormflow chemistry is controlled by either (a) the event-water signature of overland flow, or (b) the pre-event water signature of return flow. We used a high-resolution hydrochemical data set of stormflow and end-members of multiple storms in an end-member mixing analysis to determine the number of end-members needed to explain stormflow, characterize and identify potential end-members, calculate their contributions to stormflow, and develop a conceptual model of stormflow. The arrangement and relative positioning of end-members in stormflow mixing space suggest that saturation excess overland flow (26-48%) and return flow from two different subsurface storage pools (17-53%) are both similarly important for stormflow. These results suggest that pipes and fractures are important flow paths to rapidly release stored water and highlight the value of within-event resolution hydrochemical data to assess the full range and dynamics of flow paths.
Lacustrine sediments have been widely used to investigate past climatic and environmental changes on millennial to seasonal time scales. Sedimentary archives of lakes in mountainous regions may also record non-climatic events such as earthquakes. We argue herein that a set of 64 annual laminae couplets reconciles a stratigraphically inconsistent accelerator mass spectrometry (AMS) C-14 chronology in a similar to 4-m-long sediment core from Lake Mengda, in the north-eastern Tibetan Plateau. The laminations suggest the lake was formed by a large landslide, triggered by the 1927 Gulang earthquake (M = 8.0). The lake sediment sequence can be separated into three units based on lithologic, sedimentary, and isotopic characteristics. Starting from the bottom of the sequence, these are: (1) unweathered, coarse, sandy valley-floor deposits or landslide debris that pre-date the lake, (2) landslide-induced, fine-grained soil or reworked landslide debris with a high organic content, and (3) lacustrine sediments with low organic content and laminations. These annual laminations provide a high-resolution record of anthropogenic and environmental changes during the twentieth century, recording enhanced sediment input associated with two phases of construction activities. The high mean sedimentation rates of up to 4.8 mm year(-1) underscore the potential for reconstructing such distinct sediment pulses in remote, forested, and seemingly undisturbed mountain catchments.
A comprehensive hydrometeorological dataset is presented spanning the period 1 January 201131 December 2014 to improve the understanding of the hydrological processes leading to flash floods and the relation between rainfall, runoff, erosion and sediment transport in a mesoscale catchment (Auzon, 116 km(2)) of the Mediterranean region. Badlands are present in the Auzon catchment and well connected to high-gradient channels of bedrock rivers which promotes the transfer of suspended solids downstream. The number of observed variables, the various sensors involved (both in situ and remote) and the space-time resolution (similar to km(2), similar to min) of this comprehensive dataset make it a unique contribution to research communities focused on hydrometeorology, surface hydrology and erosion. Given that rainfall is highly variable in space and time in this region, the observation system enables assessment of the hydrological response to rainfall fields. Indeed, (i) rainfall data are provided by rain gauges (both a research network of 21 rain gauges with a 5 min time step and an operational network of 10 rain gauges with a 5 min or 1 h time step), S-band Doppler dual-polarization radars (1 km(2), 5 min resolution), disdrometers (16 sensors working at 30 s or 1 min time step) and Micro Rain Radars (5 sensors, 100m height resolution). Additionally, during the special observation period (SOP-1) of the HyMeX (Hydrological Cycle in the Mediterranean Experiment) project, two X-band radars provided precipitation measurements at very fine spatial and temporal scales (1 ha, 5 min). (ii) Other meteorological data are taken from the operational surface weather observation stations of Meteo-France (including 2m air temperature, atmospheric pressure, 2 m relative humidity, 10m wind speed and direction, global radiation) at the hourly time resolution (six stations in the region of interest). (iii) The monitoring of surface hydrology and suspended sediment is multi-scale and based on nested catchments. Three hydrometric stations estimate water discharge at a 2-10 min time resolution. Two of these stations also measure additional physico-chemical variables (turbidity, temperature, conductivity) and water samples are collected automatically during floods, allowing further geochemical characterization of water and suspended solids. Two experimental plots monitor overland flow and erosion at 1 min time resolution on a hillslope with vineyard. A network of 11 sensors installed in the intermittent hydrographic network continuously measures water level and water temperature in headwater subcatchments (from 0.17 to 116 km(2)) at a time resolution of 2-5 min. A network of soil moisture sensors enables the continuous measurement of soil volumetric water content at 20 min time resolution at 9 sites. Additionally, concomitant observations (soil moisture measurements and stream gauging) were performed during floods between 2012 and 2014. Finally, this dataset is considered appropriate for understanding the rainfall variability in time and space at fine scales, improving areal rainfall estimations and progressing in distributed hydrological and erosion modelling.
Reliable hydrological monitoring is the basis for sound water management in drained wetlands. Since statistical methods cannot be employed for unobserved or sparsely monitored areas, the primary design (first set-up) may be arbitrary in most instances. The objective of this paper is therefore to provide a guideline for designing the initial hydrological monitoring network. A scheme is developed that handles different parts of monitoring and hydrometry in wetlands, focusing on the positioning of surface water and groundwater gauges. For placement of the former, control units are used which correspond to areas whose water levels can be regulated separately. The latter are arranged depending on hydrological response units, defined by combinations of soil type and land use, and the chosen surface water monitoring sites. A practical application of the approach is shown for an investigation area in the Spreewald region in north-east Germany. The presented scheme leaves a certain degree of freedom to its user, allowing the inclusion of expert knowledge or special concerns. Based on easily obtainable data, the developed hydrological network serves as a first step in the iterative procedure of monitoring network optimisation. Copyright (c) 2013 John Wiley & Sons, Ltd.
Growing attention to phytoplankton mixotrophy as a trophic strategy has led to significant revisions of traditional pelagic food web models and ecosystem functioning. Although some empirical estimates of mixotrophy do exist, a much broader set of in situ measurements are required to (i) identify which organisms are acting as mixotrophs in real time and to (ii) assess the contribution of their heterotrophy to biogeochemical cycling. Estimates are needed through time and across space to evaluate which environmental conditions or habitats favour mixotrophy: conditions still largely unknown. We review methodologies currently available to plankton ecologists to undertake estimates of plankton mixotrophy, in particular nanophytoplankton phago-mixotrophy. Methods are based largely on fluorescent or isotopic tracers, but also take advantage of genomics to identify phylotypes and function. We also suggest novel methods on the cusp of use for phago-mixotrophy assessment, including single-cell measurements improving our capacity to estimate mixotrophic activity and rates in wild plankton communities down to the single-cell level. Future methods will benefit from advances in nanotechnology, micromanipulation and microscopy combined with stable isotope and genomic methodologies. Improved estimates of mixotrophy will enable more reliable models to predict changes in food web structure and biogeochemical flows in a rapidly changing world.
A ground motion logic tree for seismic hazard analysis in the stable cratonic region of Europe
(2020)
Regions of low seismicity present a particular challenge for probabilistic seismic hazard analysis when identifying suitable ground motion models (GMMs) and quantifying their epistemic uncertainty. The 2020 European Seismic Hazard Model adopts a scaled backbone approach to characterise this uncertainty for shallow seismicity in Europe, incorporating region-to-region source and attenuation variability based on European strong motion data. This approach, however, may not be suited to stable cratonic region of northeastern Europe (encompassing Finland, Sweden and the Baltic countries), where exploration of various global geophysical datasets reveals that its crustal properties are distinctly different from the rest of Europe, and are instead more closely represented by those of the Central and Eastern United States. Building upon the suite of models developed by the recent NGA East project, we construct a new scaled backbone ground motion model and calibrate its corresponding epistemic uncertainties. The resulting logic tree is shown to provide comparable hazard outcomes to the epistemic uncertainty modelling strategy adopted for the Eastern United States, despite the different approaches taken. Comparison with previous GMM selections for northeastern Europe, however, highlights key differences in short period accelerations resulting from new assumptions regarding the characteristics of the reference rock and its influence on site amplification.
Deep hydrothermal Mo, W, and base metal mineralization at the Sweet Home mine (Detroit City portal) formed in response to magmatic activity during the Oligocene. Microthermometric data of fluid inclusions trapped in greisen quartz and fluorite suggest that the early-stage mineralization at the Sweet Home mine precipitated from low- to medium-salinity (1.5-11.5 wt% equiv. NaCl), CO2-bearing fluids at temperatures between 360 and 415 degrees C and at depths of at least 3.5 km. Stable isotope and noble gas isotope data indicate that greisen formation and base metal mineralization at the Sweet Home mine was related to fluids of different origins. Early magmatic fluids were the principal source for mantle-derived volatiles (CO2, H2S/SO2, noble gases), which subsequently mixed with significant amounts of heated meteoric water. Mixing of magmatic fluids with meteoric water is constrained by delta H-2(w)-delta O-18(w) relationships of fluid inclusions. The deep hydrothermal mineralization at the Sweet Home mine shows features similar to deep hydrothermal vein mineralization at Climax-type Mo deposits or on their periphery. This suggests that fluid migration and the deposition of ore and gangue minerals in the Sweet Home mine was triggered by a deep-seated magmatic intrusion. The findings of this study are in good agreement with the results of previous fluid inclusion studies of the mineralization of the Sweet Home mine and from Climax-type Mo porphyry deposits in the Colorado Mineral Belt.