Refine
Year of publication
Document Type
- Article (2483)
- Doctoral Thesis (470)
- Postprint (112)
- Other (72)
- Review (51)
- Monograph/Edited Volume (33)
- Preprint (17)
- Conference Proceeding (12)
- Habilitation Thesis (11)
- Master's Thesis (5)
Keywords
- Holocene (43)
- climate change (38)
- erosion (24)
- Himalaya (22)
- Climate change (21)
- Pollen (20)
- Tibetan Plateau (20)
- remote sensing (20)
- climate (18)
- Erdbeben (17)
Institute
- Institut für Geowissenschaften (3273) (remove)
Glaciated high-alpine areas are fundamentally altered by climate change, with well-known implications for hydrology, e.g., due to glacier retreat, longer snow-free periods, and more frequent and intense summer rainstorms. While knowledge on how these hydrological changes will propagate to suspended sediment dynamics is still scarce, it is needed to inform mitigation and adaptation strategies. To understand the processes and source areas most relevant to sediment dynamics, we analyzed discharge and sediment dynamics in high temporal resolution as well as their patterns on several spatial scales, which to date few studies have done.
We used a nested catchment setup in the Upper Ötztal in Tyrol, Austria, where high-resolution (15 min) time series of discharge and suspended sediment concentrations are available for up to 15 years (2006–2020). The catchments of the gauges in Vent, Sölden and Tumpen range from 100 to almost 800 km2 with 10 % to 30 % glacier cover and span an elevation range of 930 to 3772 m a.s.l. We analyzed discharge and suspended sediment yields (SSY), their distribution in space, their seasonality and spatial differences therein, and the relative importance of short-term events. We complemented our analysis by linking the observations to satellite-based snow cover maps, glacier inventories, mass balances and precipitation data.
Our results indicate that the areas above 2500 m a.s.l., characterized by glacier tongues and the most recently deglaciated areas, are crucial for sediment generation in all sub-catchments. This notion is supported by the synchronous spring onset of sediment export at the three gauges, which coincides with snowmelt above 2500 m but lags behind spring discharge onsets. This points at a limitation of suspended sediment supply as long as the areas above 2500 m are snow-covered. The positive correlation of annual SSY with glacier cover (among catchments) and glacier mass balances (within a catchment) further supports the importance of the glacier-dominated areas. The analysis of short-term events showed that summer precipitation events were associated with peak sediment concentrations and yields but on average accounted for only 21 % of the annual SSY in the headwaters. These results indicate that under current conditions, thermally induced sediment export (through snow and glacier melt) is dominant in the study area.
Our results extend the scientific knowledge on current hydro-sedimentological conditions in glaciated high-alpine areas and provide a baseline for studies on projected future changes in hydro-sedimentological system dynamics.
Glaciated high-alpine areas are fundamentally altered by climate change, with well-known implications for hydrology, e.g., due to glacier retreat, longer snow-free periods, and more frequent and intense summer rainstorms. While knowledge on how these hydrological changes will propagate to suspended sediment dynamics is still scarce, it is needed to inform mitigation and adaptation strategies. To understand the processes and source areas most relevant to sediment dynamics, we analyzed discharge and sediment dynamics in high temporal resolution as well as their patterns on several spatial scales, which to date few studies have done.
We used a nested catchment setup in the Upper Ötztal in Tyrol, Austria, where high-resolution (15 min) time series of discharge and suspended sediment concentrations are available for up to 15 years (2006–2020). The catchments of the gauges in Vent, Sölden and Tumpen range from 100 to almost 800 km2 with 10 % to 30 % glacier cover and span an elevation range of 930 to 3772 m a.s.l. We analyzed discharge and suspended sediment yields (SSY), their distribution in space, their seasonality and spatial differences therein, and the relative importance of short-term events. We complemented our analysis by linking the observations to satellite-based snow cover maps, glacier inventories, mass balances and precipitation data.
Our results indicate that the areas above 2500 m a.s.l., characterized by glacier tongues and the most recently deglaciated areas, are crucial for sediment generation in all sub-catchments. This notion is supported by the synchronous spring onset of sediment export at the three gauges, which coincides with snowmelt above 2500 m but lags behind spring discharge onsets. This points at a limitation of suspended sediment supply as long as the areas above 2500 m are snow-covered. The positive correlation of annual SSY with glacier cover (among catchments) and glacier mass balances (within a catchment) further supports the importance of the glacier-dominated areas. The analysis of short-term events showed that summer precipitation events were associated with peak sediment concentrations and yields but on average accounted for only 21 % of the annual SSY in the headwaters. These results indicate that under current conditions, thermally induced sediment export (through snow and glacier melt) is dominant in the study area.
Our results extend the scientific knowledge on current hydro-sedimentological conditions in glaciated high-alpine areas and provide a baseline for studies on projected future changes in hydro-sedimentological system dynamics.
Earthquakes often rupture across more than one fault segment. If such rupture segmentation occurs on a significant scale, a simple point-source or one-fault model may not represent the rupture process well. As a consequence earthquake characteristics inferred, based on one-source assumptions, may become systematically wrong. This might have effects on follow-up analyses, for example regional stress field inversions and seismic hazard assessments. While rupture segmentation is evident for most M-w > 7 earthquakes, also smaller ones with 5.5 < M-w < 7 can be segmented. We investigate the sensitivity of globally available data sets to rupture segmentation and their resolution to reliably estimate the mechanisms in presence of segmentation. We focus on the sensitivity of InSAR (Interferometric Synthetic Aperture Radar) data in the static near-field and seismic waveforms in the far-field of the rupture and carry out non-linear and Bayesian optimizations of single-source and two-sources kinematic models (double-couple point sources and finite, rectangular sources) using InSAR and teleseismic waveforms separately. Our case studies comprises of four M-w 6-7 earthquakes: the 2009 L'Aquila and 2016 Amatrice (Italy) and the 2005 and 2008 Zhongba (Tibet) earthquakes. We contrast the data misfits of different source complexity by using the Akaike informational criterion (AIC). We find that the AIC method is well suited for data-driven inferences on significant rupture segmentation for the given data sets. This is based on our observation that an AIC-stated significant improvement of data fit for two-segment models over one-segment models correlates with significantly different mechanisms of the two source segments and their average compared to the single-segment mechanism. We attribute these modelled differences to a sufficient sensitivity of the data to resolve rupture segmentation. Our results show that near-field data are generally more sensitive to rupture segmentation of shallow earthquakes than far-field data but that also teleseismic data can resolve rupture segmentation in the studied magnitude range. We further conclude that a significant difference in the modelled source mechanisms for different segmentations shows that an appropriate choice of model segmentation matters for a robust estimation of source mechanisms. It reduces systematic biases and trade-off and thereby improves the knowledge on the rupture. Our study presents a strategy and method to detect significant rupture segmentation such that an appropriate model complexity can be used in the source mechanism inference. A similar, systematic investigation of earthquakes in the range of M-w 5.5-7 could provide important hazard-relevant statistics on rupture segmentation. In these cases single-source models introduce a systematic bias. Consideration of rupture segmentation therefore matters for a robust estimation of source mechanisms of the studied earthquakes.
A comprehensive workflow to analyze ensembles of globally inverted 2D electrical resistivity models
(2022)
Electrical resistivity tomography (ERT) aims at imaging the subsurface resistivity distribution and provides valuable information for different geological, engineering, and hydrological applications. To obtain a subsurface resistivity model from measured apparent resistivities, stochastic or deterministic inversion procedures may be employed. Typically, the inversion of ERT data results in non-unique solutions; i.e., an ensemble of different models explains the measured data equally well. In this study, we perform inference analysis of model ensembles generated using a well-established global inversion approach to assess uncertainties related to the nonuniqueness of the inverse problem. Our interpretation strategy starts by establishing model selection criteria based on different statistical descriptors calculated from the data residuals. Then, we perform cluster analysis considering the inverted resistivity models and the corresponding data residuals. Finally, we evaluate model uncertainties and residual distributions for each cluster. To illustrate the potential of our approach, we use a particle swarm optimization (PSO) algorithm to obtain an ensemble of 2D layer-based resistivity models from a synthetic data example and a field data set collected in Loon-Plage, France. Our strategy performs well for both synthetic and field data and allows us to extract different plausible model scenarios with their associated uncertainties and data residual distributions. Although we demonstrate our workflow using 2D ERT data and a PSObased inversion approach, the proposed strategy is general and can be adapted to analyze model ensembles generated from other kinds of geophysical data and using different global inversion approaches.
The simulation of broad-band (0.1 to 10 + Hz) ground-shaking over deep and spatially extended sedimentary basins at regional scales is challenging. We evaluate the ground-shaking of a potential M 6.5 earthquake in the southern Lower Rhine Embayment, one of the most important areas of earthquake recurrence north of the Alps, close to the city of Cologne in Germany. In a first step, information from geological investigations, seismic experiments and boreholes is combined for deriving a harmonized 3D velocity and attenuation model of the sedimentary layers. Three alternative approaches are then applied and compared to evaluate the impact of the sedimentary cover on ground-motion amplification. The first approach builds on existing response spectra ground-motion models whose amplification factors empirically take into account the influence of the sedimentary layers through a standard parameterization. In the second approach, site-specific 1D amplification functions are computed from the 3D basin model. Using a random vibration theory approach, we adjust the empirical response spectra predicted for soft rock conditions by local site amplification factors: amplifications and associated ground-motions are predicted both in the Fourier and in the response spectra domain. In the third approach, hybrid physics-based ground-motion simulations are used to predict time histories for soft rock conditions which are subsequently modified using the 1D site-specific amplification functions computed in method 2. For large distances and at short periods, the differences between the three approaches become less notable due to the significant attenuation of the sedimentary layers. At intermediate and long periods, generic empirical ground-motion models provide lower levels of amplification from sedimentary soils compared to methods taking into account site-specific 1D amplification functions. In the near-source region, hybrid physics-based ground-motions models illustrate the potentially large variability of ground-motion due to finite source effects.
Thirteen N-butylpyridinium salts, including three monometallic [C4Py](2)[MCl4], nine bimetallic [C4Py](2)[(M1-xMxCl4)-M-a-Cl-b] and one trimetallic compound [C4Py](2)[(M1-y-zMyMz (c) Cl4)-M-a-M-b] (M=Co, Cu, Mn; x=0.25, 0.50 or 0.75 and y=z=0.33), were synthesized and their structure and thermal and electrochemical properties were studied. All compounds are ionic liquids (ILs) with melting points between 69 and 93 degrees C. X-ray diffraction proves that all ILs are isostructural. The conductivity at room temperature is between 10(-4) and 10(-8) S cm(-1). Some Cu-based ILs reach conductivities of 10(-2) S cm(-1), which is, however, probably due to IL dec. This correlates with the optical bandgap measurements indicating the formation of large bandgap semiconductors. At elevated temperatures approaching the melting points, the conductivities reach up to 1.47x10(-1) S cm(-1) at 70 degrees C. The electrochemical stability windows of the ILs are between 2.5 and 3.0 V.
Borehole leakage is a common and complex issue. Understanding the fluid flow characteristics of a cemented area inside a borehole is crucial to monitor and quantify the wellbore integrity as well as to find solutions to minimise existing leakages. In order to improve our understanding of the flow behaviour of cemented boreholes, we investigated experimental data of a large-scale borehole leakage tests by means of numerical modelling using three different conceptual models. The experiment was performed with an autoclave system consisting of two vessels bridged by a cement-filled casing. After a partial bleed-off at the well-head, a sustained casing pressure was observed due to fluid flow through the cementsteel composite. The aim of our simulations is to investigate and quantify the permeability of the cement-steel composite. From our model results, we conclude that the flow occurred along a preferential flow path at the cement-steel interface. Thus, the inner part of the cement core was impermeable during the duration of the experiment. The preferential flow path can be described as a highly permeable and highly porous area with an aperture of about 5 mu m and a permeability of 3 . 10(-12) m(2) (3 Darcy). It follows that the fluid flow characteristics of a cemented area inside a borehole cannot be described using one permeability value for the entire cement-steel composite. Furthermore, it can be concluded that the quality of the cement and the filling process regarding the cement-steel interface is crucial to minimize possible well leakages.
Flood loss data collection and modeling are not standardized, and previous work has indicated that losses from different flood types (e.g., riverine and groundwater) may follow different driving forces. However, different flood types may occur within a single flood event, which is known as a compound flood event. Therefore, we aimed to identify statistical similarities between loss-driving factors across flood types and test whether the corresponding losses should be modeled separately. In this study, we used empirical data from 4,418 respondents from four survey campaigns studying households in Germany that experienced flooding. These surveys sought to investigate several features of the impact process (hazard, socioeconomic, preparedness, and building characteristics, as well as flood type). While the level of most of these features differed across flood type subsamples (e.g., degree of preparedness), they did so in a nonregular pattern. A variable selection process indicates that besides hazard and building characteristics, information on property-level preparedness was also selected as a relevant predictor of the loss ratio. These variables represent information, which is rarely adopted in loss modeling. Models shall be refined with further data collection and other statistical methods. To save costs, data collection efforts should be steered toward the most relevant predictors to enhance data availability and increase the statistical power of results. Understanding that losses from different flood types are driven by different factors is a crucial step toward targeted data collection and model development and will finally clarify conditions that allow us to transfer loss models in space and time. <br /> Key Points <br /> Survey data of flood-affected households show different concurrent flood types, undermining the use of a single-flood-type loss model Thirteen variables addressing flood hazard, the building, and property level preparedness are significant predictors of the building loss ratio Flood type-specific models show varying significance across the predictor variables, indicating a hindrance to model transferability
Una extraordinaria edición actualizada de Observaciones sobre el clima de Lima y sus influencias en los seres organizados, en especial el hombre de José Hipólito Unanue y Pavón ha sido publicado en el 2018 y merece ser anunciado a la comunidad académica peruana e internacional. Este libro incluye un exhaustivo estudio in troductorio del destacado historiador Lizardo Seiner Lizárraga. Esta introducción está conformada de tres partes encabezadas por los siguientes títulos: a) Unanue: una historiografía de amplio espectro b) Textos, registros y referencias: una indagación en la biblioteca de Unanue y c) Observaciones sobre el clima de Lima: una obra en dos tiempos.
The first step towards assessing hazards in seismically active regions involves mapping capable faults and estimating their recurrence times. While the mapping of active faults is commonly based on distinct geologic and geomorphic features evident at the surface, mapping blind seismogenic faults is complicated by the absence of on-fault diagnostic features. Here we investigated the Pichilemu Fault in coastal Chile, unknown until it generated a Mw 7.0 earthquake in 2010. The lack of evident surface faulting suggests activity along a partly-hidden blind fault. We used off-fault deformed marine terraces to estimate a fault-slip rate of 0.52 ± 0.04 m/ka, which, when integrated with satellite geodesy suggests a 2.12 ± 0.2 ka recurrence time for Mw~7.0 normal-faulting earthquakes. We propose that extension in the Pichilemu region is associated with stress changes during megathrust earthquakes and accommodated by sporadic slip during upper-plate earthquakes, which has implications for assessing the seismic potential of cryptic faults along convergent margins and elsewhere.
The first step towards assessing hazards in seismically active regions involves mapping capable faults and estimating their recurrence times. While the mapping of active faults is commonly based on distinct geologic and geomorphic features evident at the surface, mapping blind seismogenic faults is complicated by the absence of on-fault diagnostic features. Here we investigated the Pichilemu Fault in coastal Chile, unknown until it generated a Mw 7.0 earthquake in 2010. The lack of evident surface faulting suggests activity along a partly-hidden blind fault. We used off-fault deformed marine terraces to estimate a fault-slip rate of 0.52 ± 0.04 m/ka, which, when integrated with satellite geodesy suggests a 2.12 ± 0.2 ka recurrence time for Mw~7.0 normal-faulting earthquakes. We propose that extension in the Pichilemu region is associated with stress changes during megathrust earthquakes and accommodated by sporadic slip during upper-plate earthquakes, which has implications for assessing the seismic potential of cryptic faults along convergent margins and elsewhere.
Understanding the relation between injection-induced seismic moment release and operational parameters is crucial for early identification of possible seismic hazards associated with fluid-injection projects. We conducted laboratory fluid-injection experiments on permeable sandstone samples containing a critically stressed fault at different fluid pressurization rates. The observed fluid-induced fault deformation is dominantly aseismic. Fluid-induced stick-slip and fault creep reveal that total seismic moment release of acoustic emission (AE) events is related to total injected volume, independent of respective fault slip behavior. Seismic moment release rate of AE scales with measured fault slip velocity. For injection-induced fault slip in a homogeneous pressurized region, released moment shows a linear scaling with injected volume for stable slip (steady slip and fault creep), while we find a cubic relation for dynamic slip. Our results highlight that monitoring evolution of seismic moment release with injected volume in some cases may assist in discriminating between stable slip and unstable runaway ruptures.
The Big Naryn Complex (BNC) in the East Djetim-Too Range of the Kyrgyz Middle Tianshan block is a tectonized, at least 2 km thick sequence of predominantly felsic to intermediate volcanic rocks intruded by porphyric rhyolite sills. It overlies a basement of metamorphic rocks and is overlain by late Neoproterozoic Djetim-Too Formation sediments; these also occur as tectonic intercalations in the BNC. The up to ca. 1100 m thick Lower Member is composed of predominantly rhyolites-to-dacites and minor basalts, while the at least 900 m thick pyroclastic Upper Member is dominated by rhyolitic-to-dacitic ignimbrites. Porphyric rhyolite sills are concentrated at the top of the Lower Member. A Lower Member rhyolite and a sill sample have LA-ICP-MS U-Pb zircon crystallization ages of 726.1 +/- 2.2 Ma and 720.3 +/- 6.5 Ma, respectively, showing that most of the magmatism occurred within a short time span in the late Tonian-early Cryogenian. Inherited zircons in the sill sample have Neoarchean (2.63, 2.64 Ga), Paleo- (2.33-1.81 Ga), Meso- (1.55 Ga), and Neoproterozoic (ca. 815 Ma) ages, and were derived from a heterogeneous Kuilyu Complex basement. A 1751 +/- 7 Ma Ar-40/Ar-39 age for amphibole from metagabbro is the age of cooling subsequent to Paleoproterozoic metamorphism of the Kuilyu Complex. The large amount of pyroclastic rocks, and their major and trace element compositions, the presence of Neoarchean to Neoproterozoic inherited zircons and a depositional basement of metamorphic rocks point to formation of the BNC in a continental magmatic arc setting.
Meteorological and hydrological drought assessment in Lake Malawi and Shire River basins (1970-2013)
(2020)
The study assesses the variability and trends of both meteorological and hydrological droughts from 1970 to 2013 in Lake Malawi and Shire River basins using the standardized precipitation index (SPI) and standardized precipitation and evaporation index (SPEI) for meteorological droughts and the lake level change index (LLCI) for hydrological droughts. Trends and slopes in droughts and drought drivers are estimated using Mann-Kendall test and Sen's slope, respectively. Results suggest that meteorological droughts are increasing due to a decrease in precipitation which is exacerbated by an increase in temperature (potential evapotranspiration). The hydrological system of Lake Malawi seems to have a >24-month memory towards meteorological conditions, since the 36-month SPEI can predict hydrological droughts 10 months in advance. The study has found the critical lake level that would trigger hydrological drought to be 474.1 m a.s.l. The increase in drought is a concern as this will have serious impacts on water resources and hydropower supply in Malawi.
Filling in for Governments?
(2020)
The 2015 Paris Agreement on climate change abandons the Kyoto Protocol's paradigm of binding emissions targets and relies instead on countries' voluntary contributions. However, the Paris Agreement encourages not only governments but also sub-national governments, corporations and civil society to contribute to reaching ambitious climate goals. In a transition from the regulated architecture of the Kyoto Protocol to the open system of the Paris Agreement, the Agreement seeks to integrate non-state actors into the treaty-based climate regime. In 2014 the secretariat of the United Nations Framework Convention on Climate Change Peru and France created the Non-State Actor Zone for Climate Action (and launched the Global Climate Action portal). In December 2019, this portal recorded more than twenty thousand climate-commitments of private and public non-state entities, making the non-state venues of international climate meetings decisively more exciting than the formal negotiation space. This level engagement and governments' response to it raises a flurry of questions in relation to the evolving nature of the climate regime and climate change governance, including the role of private actors as standard setters and the lack of accountability mechanisms for non-state actions. This paper takes these developments as occasion to discuss the changing role of private actors in the climate regime.
Ground-penetrating radar (GPR) is an established geophysical tool to explore a wide range of near-surface environments. Today, the use of synthetic GPR data is largely limited to 2D because 3D modeling is computationally more expensive. In fact, only recent developments of modeling tools and powerful hardware allow for a time-efficient computation of extensive 3D data sets. Thus, 3D subsurface models and resulting GPR data sets, which are of great interest to develop and evaluate novel approaches in data analysis and interpretation, have not been made publicly available up to now. <br /> We use a published hydrofacies data set of an aquifer-analog study within fluvio-glacial deposits to infer a realistic 3D porosity model showing heterogeneities at multiple spatial scales. Assuming fresh-water saturated sediments, we generate synthetic 3D GPR data across this model using novel GPU-acceleration included in the open-source software gprMax. We present a numerical approach to examine 3D wave-propagation effects in modeled GPR data. Using the results of this examination study, we conduct a spatial model decomposition to enable a computationally efficient 3D simulation of a typical GPR reflection data set across the entire model surface. We process the resulting GPR data set using a standard 3D structural imaging sequence and compare the results to selected input data to demonstrate the feasibility and potential of the presented modeling studies. We conclude on conceivable applications of our 3D GPR reflection data set and the underlying porosity model, which are both publicly available and, thus, can support future methodological developments in GPR and other near-surface geophysical techniques.
The affordability of property-level adaptation measures against flooding is crucial due to the movement toward integrated flood risk management, which requires the individuals threatened by flooding to actively manage flooding. It is surprising to find that affordability is not often discussed, given the important roles that affordability and social justice play regarding flood risk management. This article provides a starting point for investigating the potential rate of unaffordability of flood risk property-level adaptation measures across Europe using two definitions of affordability, which are combined with two different affordability thresholds from within flood risk research. It uses concepts of investment and payment affordability, with affordability thresholds based on residual income and expenditure definitions of unaffordability. These concepts, in turn, are linked with social justice through fairness concerns, in that, all should have equal capability to act, of which affordability is one avenue. In doing so, it was found that, for a large proportion of Europe, property owners generally cannot afford to make one-time payment of the cost of protective measures. These can be made affordable with installment payment mechanisms or similar mechanisms that spread costs over time. Therefore, the movement toward greater obligations for flood-prone residents to actively adapt to flooding should be accompanied by socially accessible financing mechanisms.
Many institutions struggle to tap into the potential of their large archives of radar reflectivity: these data are often affected by miscalibration, yet the bias is typically unknown and temporally volatile. Still, relative calibration techniques can be used to correct the measurements a posteriori. For that purpose, the usage of spaceborne reflectivity observations from the Tropical Rainfall Measuring Mission (TRMM) and Global Precipitation Measurement (GPM) platforms has become increasingly popular: the calibration bias of a ground radar (GR) is estimated from its average reflectivity difference to the spaceborne radar (SR). Recently, Crisologo et al. (2018) introduced a formal procedure to enhance the reliability of such estimates: each match between SR and GR observations is assigned a quality index, and the calibration bias is inferred as a quality-weighted average of the differences between SR and GR. The relevance of quality was exemplified for the Subic S-band radar in the Philippines, which is greatly affected by partial beam blockage. The present study extends the concept of quality-weighted averaging by accounting for path-integrated attenuation (PIA) in addition to beam blockage. This extension becomes vital for radars that operate at the C or X band. Correspondingly, the study setup includes a C-band radar that substantially overlaps with the S-band radar. Based on the extended quality-weighting approach, we retrieve, for each of the two ground radars, a time series of calibration bias estimates from suitable SR overpasses. As a result of applying these estimates to correct the ground radar observations, the consistency between the ground radars in the region of overlap increased substantially. Furthermore, we investigated if the bias estimates can be interpolated in time, so that ground radar observations can be corrected even in the absence of prompt SR overpasses. We found that a moving average approach was most suitable for that purpose, although limited by the absence of explicit records of radar maintenance operations.
River-valley morphology preserves information on tectonic and climatic conditions that shape landscapes. Observations suggest that river discharge and valley-wall lithology are the main controls on valley width. Yet, current models based on these observations fail to explain the full range of cross-sectional valley shapes in nature, suggesting hitherto unquantified controls on valley width. In particular, current models cannot explain the existence of paired terrace sequences that form under cyclic climate forcing. Paired river terraces are staircases of abandoned floodplains on both valley sides, and hence preserve past valley widths. Their formation requires alternating phases of predominantly river incision and predominantly lateral planation, plus progressive valley narrowing. While cyclic Quaternary climate changes can explain shifts between incision and lateral erosion, the driving mechanism of valley narrowing is unknown. Here, we extract valley geometries from climatically formed, alluvial river-terrace sequences and show that across our dataset, the total cumulative terrace height (here: total valley height) explains 90%–99% of the variance in valley width at the terrace sites. This finding suggests that valley height, or a parameter that scales linearly with valley height, controls valley width in addition to river discharge and lithology. To explain this valley-width-height relationship, we reformulate existing valley-width models and suggest that, when adjusting to new boundary conditions, alluvial valleys evolve to a width at which sediment removal from valley walls matches lateral sediment supply from hillslope erosion. Such a hillslope-channel coupling is not captured in current valley-evolution models. Our model can explain the existence of paired terrace sequences under cyclic climate forcing and relates valley width to measurable field parameters. Therefore, it facilitates the reconstruction of past climatic and tectonic conditions from valley topography.
River-valley morphology preserves information on tectonic and climatic conditions that shape landscapes. Observations suggest that river discharge and valley-wall lithology are the main controls on valley width. Yet, current models based on these observations fail to explain the full range of cross-sectional valley shapes in nature, suggesting hitherto unquantified controls on valley width. In particular, current models cannot explain the existence of paired terrace sequences that form under cyclic climate forcing. Paired river terraces are staircases of abandoned floodplains on both valley sides, and hence preserve past valley widths. Their formation requires alternating phases of predominantly river incision and predominantly lateral planation, plus progressive valley narrowing. While cyclic Quaternary climate changes can explain shifts between incision and lateral erosion, the driving mechanism of valley narrowing is unknown. Here, we extract valley geometries from climatically formed, alluvial river-terrace sequences and show that across our dataset, the total cumulative terrace height (here: total valley height) explains 90%–99% of the variance in valley width at the terrace sites. This finding suggests that valley height, or a parameter that scales linearly with valley height, controls valley width in addition to river discharge and lithology. To explain this valley-width-height relationship, we reformulate existing valley-width models and suggest that, when adjusting to new boundary conditions, alluvial valleys evolve to a width at which sediment removal from valley walls matches lateral sediment supply from hillslope erosion. Such a hillslope-channel coupling is not captured in current valley-evolution models. Our model can explain the existence of paired terrace sequences under cyclic climate forcing and relates valley width to measurable field parameters. Therefore, it facilitates the reconstruction of past climatic and tectonic conditions from valley topography.
Sea level rise and coastal erosion have inundated large areas of Arctic permafrost. Submergence by warm and saline waters increases the rate of inundated permafrost thaw compared to sub-aerial thawing on land. Studying the contact between the unfrozen and frozen sediments below the seabed, also known as the ice-bearing permafrost table (IBPT), provides valuable information to understand the evolution of sub-aquatic permafrost, which is key to improving and understanding coastal erosion prediction models and potential greenhouse gas emissions. In this study, we use data from 2D electrical resistivity tomography (ERT) collected in the nearshore coastal zone of two Arctic regions that differ in their environmental conditions (e.g., seawater depth and resistivity) to image and study the subsea permafrost. The inversion of 2D ERT data sets is commonly performed using deterministic approaches that favor smoothed solutions, which are typically interpreted using a user-specified resistivity threshold to identify the IBPT position. In contrast, to target the IBPT position directly during inversion, we use a layer-based model parameterization and a global optimization approach to invert our ERT data. This approach results in ensembles of layered 2D model solutions, which we use to identify the IBPT and estimate the resistivity of the unfrozen and frozen sediments, including estimates of uncertainties. Additionally, we globally invert 1D synthetic resistivity data and perform sensitivity analyses to study, in a simpler way, the correlations and influences of our model parameters. The set of methods provided in this study may help to further exploit ERT data collected in such permafrost environments as well as for the design of future field experiments.
Sea level rise and coastal erosion have inundated large areas of Arctic permafrost. Submergence by warm and saline waters increases the rate of inundated permafrost thaw compared to sub-aerial thawing on land. Studying the contact between the unfrozen and frozen sediments below the seabed, also known as the ice-bearing permafrost table (IBPT), provides valuable information to understand the evolution of sub-aquatic permafrost, which is key to improving and understanding coastal erosion prediction models and potential greenhouse gas emissions. In this study, we use data from 2D electrical resistivity tomography (ERT) collected in the nearshore coastal zone of two Arctic regions that differ in their environmental conditions (e.g., seawater depth and resistivity) to image and study the subsea permafrost. The inversion of 2D ERT data sets is commonly performed using deterministic approaches that favor smoothed solutions, which are typically interpreted using a user-specified resistivity threshold to identify the IBPT position. In contrast, to target the IBPT position directly during inversion, we use a layer-based model parameterization and a global optimization approach to invert our ERT data. This approach results in ensembles of layered 2D model solutions, which we use to identify the IBPT and estimate the resistivity of the unfrozen and frozen sediments, including estimates of uncertainties. Additionally, we globally invert 1D synthetic resistivity data and perform sensitivity analyses to study, in a simpler way, the correlations and influences of our model parameters. The set of methods provided in this study may help to further exploit ERT data collected in such permafrost environments as well as for the design of future field experiments.
The Ice, Cloud, and Land Elevation Satellite-2 (ICESat-2) with its land and vegetation height data product (ATL08), and Global Ecosystem Dynamics Investigation (GEDI) with its terrain elevation and height metrics data product (GEDI Level 2A) missions have great potential to globally map ground and canopy heights. Canopy height is a key factor in estimating above-ground biomass and its seasonal changes; these satellite missions can also improve estimated above-ground carbon stocks. This study presents a novel Sparse Vegetation Detection Algorithm (SVDA) which uses ICESat-2 (ATL03, geolocated photons) data to map tree and vegetation heights in a sparsely vegetated savanna ecosystem. The SVDA consists of three main steps: First, noise photons are filtered using the signal confidence flag from ATL03 data and local point statistics. Second, we classify ground photons based on photon height percentiles. Third, tree and grass photons are classified based on the number of neighbors. We validated tree heights with field measurements (n = 55), finding a root-mean-square error (RMSE) of 1.82 m using SVDA, GEDI Level 2A (Geolocated Elevation and Height Metrics product): 1.33 m, and ATL08: 5.59 m. Our results indicate that the SVDA is effective in identifying canopy photons in savanna ecosystems, where ATL08 performs poorly. We further identify seasonal vegetation height changes with an emphasis on vegetation below 3 m; widespread height changes in this class from two wet-dry cycles show maximum seasonal changes of 1 m, possibly related to seasonal grass-height differences. Our study shows the difficulties of vegetation measurements in savanna ecosystems but provides the first estimates of seasonal biomass changes.
The Ice, Cloud, and Land Elevation Satellite-2 (ICESat-2) with its land and vegetation height data product (ATL08), and Global Ecosystem Dynamics Investigation (GEDI) with its terrain elevation and height metrics data product (GEDI Level 2A) missions have great potential to globally map ground and canopy heights. Canopy height is a key factor in estimating above-ground biomass and its seasonal changes; these satellite missions can also improve estimated above-ground carbon stocks. This study presents a novel Sparse Vegetation Detection Algorithm (SVDA) which uses ICESat-2 (ATL03, geolocated photons) data to map tree and vegetation heights in a sparsely vegetated savanna ecosystem. The SVDA consists of three main steps: First, noise photons are filtered using the signal confidence flag from ATL03 data and local point statistics. Second, we classify ground photons based on photon height percentiles. Third, tree and grass photons are classified based on the number of neighbors. We validated tree heights with field measurements (n = 55), finding a root-mean-square error (RMSE) of 1.82 m using SVDA, GEDI Level 2A (Geolocated Elevation and Height Metrics product): 1.33 m, and ATL08: 5.59 m. Our results indicate that the SVDA is effective in identifying canopy photons in savanna ecosystems, where ATL08 performs poorly. We further identify seasonal vegetation height changes with an emphasis on vegetation below 3 m; widespread height changes in this class from two wet-dry cycles show maximum seasonal changes of 1 m, possibly related to seasonal grass-height differences. Our study shows the difficulties of vegetation measurements in savanna ecosystems but provides the first estimates of seasonal biomass changes.
Seismology, like many scientific fields, e.g., music information retrieval and speech signal pro- cessing, is experiencing exponential growth in the amount of data acquired by modern seismo- logical networks. In this thesis, I take advantage of the opportunities offered by "big data" and by the methods developed in the areas of music information retrieval and machine learning to predict better the ground motion generated by earthquakes and to study the properties of the surface layers of the Earth. In order to better predict seismic ground motions, I propose two approaches based on unsupervised deep learning methods, an autoencoder network and Generative Adversarial Networks. The autoencoder technique explores a massive amount of ground motion data, evaluates the required parameters, and generates synthetic ground motion data in the Fourier amplitude spectra (FAS) domain. This method is tested on two synthetic datasets and one real dataset. The application on the real dataset shows that the substantial information contained within the FAS data can be encoded to a four to the five-dimensional manifold. Consequently, only a few independent parameters are required for efficient ground motion prediction. I also propose a method based on Conditional Generative Adversarial Networks (CGAN) for simulating ground motion records in the time-frequency and time domains. CGAN generates the time-frequency domains based on the parameters: magnitude, distance, and shear wave velocities to 30 m depth (VS30). After generating the amplitude of the time-frequency domains using the CGAN model, instead of classical conventional methods that assume the amplitude spectra with a random phase spectrum, the phase of the time-frequency domains is recovered by minimizing the observed and reconstructed spectrograms. In the second part of this dissertation, I propose two methods for the monitoring and characterization of near-surface materials and site effect analyses. I implement an autocorrelation function and an interferometry method to monitor the velocity changes of near-surface materials resulting from the Kumamoto earthquake sequence (Japan, 2016). The observed seismic velocity changes during the strong shaking are due to the non-linear response of the near-surface materials. The results show that the velocity changes lasted for about two months after the Kumamoto mainshock. Furthermore, I used the velocity changes to evaluate the in-situ strain-stress relationship. I also propose a method for assessing the site proxy "VS30" using non-invasive analysis. In the proposed method, a dispersion curve of surface waves is inverted to estimate the shear wave velocity of the subsurface. This method is based on the Dix-like linear operators, which relate the shear wave velocity to the phase velocity. The proposed method is fast, efficient, and stable. All of the methods presented in this work can be used for processing "big data" in seismology and for the analysis of weak and strong ground motion data, to predict ground shaking, and to analyze site responses by considering potential time dependencies and nonlinearities.
A local and flexible definition of the monsoon season based on hydrological evidence is important for the understanding and management of Himalayan water resources. Here, we present an objective statistical method to retrieve seasonal hydrometeorological transitions. Applied to daily rainfall data (1951-2015), this method shows an average longitudinal delay of similar to 15 days, with later monsoon onset and earlier withdrawal in the western Himalaya, consistent with the continental progression of wet air masses. This delay leads to seasons of different length along the Himalaya and biased precipitation amounts when using uniform calendric monsoon boundaries. In the Central Himalaya annual precipitation has increased, due primarily to an increase of premonsoon precipitation. These findings highlight issues associated with a static definition of monsoon boundaries and call for a deeper understanding of nonmonsoonal precipitation over the Himalayan water tower. <br /> Plain Language Summary Precipitation in the Himalayas determines water availability for the Indian foreland with large socioeconomic implications. Despite its importance, spatial and temporal patterns of precipitation are poorly understood. Here, we estimate the long-term average and trends of seasonal precipitation at the scale of individual catchments draining the Himalayas. We apply a statistical method to detect the timing of hydrometeorological seasons from local precipitation measurements, focusing on monsoon onset and withdrawal. We identify longitudinal and latitudinal delays, resulting in seasons of different length along and across the Himalayas. These spatial patterns and the annual variability of the monsoon boundaries mean that oft-used, fixed calendric dates, for example, 1 June to 30 September, may be inadequate for retrieving monsoon rainfall totals. Moreover, we find that, despite its prominent contribution to annual rainfall totals, the Indian summer monsoon cannot explain the increase of the annual precipitation over the Central Himalayas. Instead, this appears to be mostly driven by changes in premonsoon and winter rainfall. So far, little attention has been paid to premonsoon precipitation, but governed by evaporative processes and surface water availability, it may be enhanced by irrigation and changed land use in the Gangetic foreland.
Most hydrological studies rely on a model calibrated using discharge alone. However, judging the model reliability based on such calibration is problematic, as it does not guarantee the correct representation of internal hydrological processes. This study aims (a) to develop a comprehensive multi-objective calibration framework using remote sensing vegetation data and hydrological signatures (flow duration curve - FDC, and baseflow index) in addition to discharge, and (b) to apply this framework for calibration of the Soil and Water Assessment Tool (SWAT) in a typical Andean catchment. Overall, our calibration approach outperformed traditional discharge-based and FDC signature-based calibration strategies in terms of vegetation, streamflow, and flow partitioning simulation. New hydrological insights for the region are the following: baseflow is the main component of the streamflow sustaining the long dry-season flow, and pasture areas offer higher water yield and baseflow than other land-cover types. The proposed approach could be used in other data-scarce regions with complex topography.
Sediment transport domains in mountain landscapes are characterized by fundamentally different processes and rates depending on several factors, including geology, climate, and biota. Accurately identifying where transitions between transport domains occur is an important step to quantify the past, present, and future contribution of varying erosion and sedimentation processes and enhance our predictive capabilities. We propose a new methodology based on time series of synthetic aperture radar (SAR) interferometric coherence images to map sediment transport regimes across arid and semiarid landscapes. Using 4 years of Sentinel-1 data, we analyze sediment transport regimes for the south-central Andes in northwestern Argentina characterized by steep topographic and climatic gradients. We observe seasonally low coherence during the regional wet season, particularly on hillslopes and in alluvial channels. The spatial distribution of coherence is compared to drainage areas extracted from digital topography to identify two distinct transitions within watersheds: (a) a hillslope-to-fluvial and (b) a fluvial-to-alluvial transition. While transitions within a given basin can be well-constrained, the relative role of each sediment transport domain varies widely over the climatic and topographic gradients. In semiarid regions, we observe larger relative contributions from hillslopes compared to arid regions. Across regional gradients, the range of coherence within basins positively correlates to previously published millennial catchment-wide erosion rates and to topographic metrics used to indicate long-term uplift. Our study suggests that a dense time series of interferometric coherence can be used as a proxy for surface sediment movement and landscape stability in vegetation-free settings at event to decadal timescales.
Large rock slope failures play a pivotal role in long-term landscape evolution and are a major concern in land use planning and hazard aspects. While the failure phase and the time immediately prior to failure are increasingly well studied, the nature of the preparation phase remains enigmatic. This knowledge gap is due, to a large degree, to difficulties associated with instrumenting high mountain terrain and the local nature of classic monitoring methods, which does not allow integral observation of large rock volumes. Here, we analyse data from a small network of up to seven seismic sensors installed during July-October 2018 (with 43 days of data loss) at the summit of the Hochvogel, a 2592 m high Alpine peak. We develop proxy time series indicative of cyclic and progressive changes of the summit. Modal analysis, horizontal-to-vertical spectral ratio data and end-member modelling analysis reveal diurnal cycles of increasing and decreasing coupling stiffness of a 260,000 m(3) large, instable rock volume, due to thermal forcing. Relative seismic wave velocity changes also indicate diurnal accumulation and release of stress within the rock mass. At longer time scales, there is a systematic superimposed pattern of stress increased over multiple days and episodic stress release within a few days, expressed in an increased emission of short seismic pulses indicative of rock cracking. Our data provide essential first order information on the development of large-scale slope instabilities towards catastrophic failure. (c) 2020 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd
Deep geological repositories represent a promising solution for the final disposal of nuclear waste. Due to its low permeability, high sorption capacity and self-sealing potential, Opalinus Clay (OPA) is considered a suitable host rock formation for the long-term storage of nuclear waste in Switzerland and Germany. However, the clay formation is characterized by compositional and structural variabilities including the occurrence of carbonate- and quartz-rich layers, pronounced bedding planes as well as tectonic elements such as pre-existing fault zones and fractures, suggesting heterogeneous rock mass properties.
Characterizing the heterogeneity of host rock properties is therefore essential for safety predictions of future repositories. This includes a detailed understanding of the mechanical and hydraulic properties, deformation behavior and the underlying deformation processes for an improved assessment of the sealing integrity and long-term safety of a deep repository in OPA. Against this background, this thesis presents the results of deformation experiments performed on intact and artificially fractured specimens of the quartz-rich, sandy and clay-rich, shaly facies of OPA. The experiments focus on the influence of mineralogical composition on the deformation behavior as well as the reactivation and sealing properties of pre-existing faults and fractures at different boundary conditions (e.g., pressure, temperature, strain rate).
The anisotropic mechanical properties of the sandy facies of OPA are presented in the first section, which were determined from triaxial deformation experiments using dried and resaturated samples loaded at 0°, 45° and 90° to the bedding plane orientation. A Paterson-type deformation apparatus was used that allowed to investigate how the deformation behavior is influenced by the variation of confining pressure (50 – 100 MPa), temperature (25 – 200 °C), and strain rate (1 × 10-3 – 5 × 10-6 s-1). Constant strain rate experiments revealed brittle to semi-brittle deformation behavior of the sandy facies at the applied conditions. Deformation behavior showed a strong dependence on confining pressure, degree of water saturation as well as bedding orientation, whereas the variation of temperature and strain rate had no significant effect on deformation. Furthermore, the sandy facies displays higher strength and stiffness compared to the clay-rich shaly facies deformed at similar conditions by Nüesch (1991). From the obtained results it can be concluded that cataclastic mechanisms dominate the short-term deformation behavior of dried samples from both facies up to elevated pressure (<200 MPa) and temperature (<200 °C) conditions.
The second part presents triaxial deformation tests that were performed to investigate how structural discontinuities affect the deformation behavior of OPA and how the reactivation of preexisting faults is influenced by mineral composition and confining pressure. To this end, dried cylindrical samples of the sandy and shaly facies of OPA were used, which contained a saw-cut fracture oriented at 30° to the long axis. After hydrostatic pre-compaction at 50 MPa, constant strain rate deformation tests were performed at confining pressures of 5, 20 or 35 MPa. With increasing confinement, a gradual transition from brittle, highly localized fault slip including a stress drop at fault reactivation to semi-brittle deformation behavior, characterized by increasing delocalization and non-linear strain hardening without dynamic fault reactivation, can be observed. Brittle localization was limited by the confining pressure at which the fault strength exceeded the matrix yield strength, above which strain partitioning between localized fault slip and distributed matrix deformation occurred. The sandy facies displayed a slightly higher friction coefficient (≈0.48) compared to the shaly facies (≈0.4). In addition, slide-hold-slide tests were conducted, revealing negative or negligible frictional strengthening, which suggests stable creep and long-term weakness of faults in both facies of OPA. The conducted experiments demonstrate that dilatant brittle fault reactivation in OPA may be favored at high overconsolidation ratios and shallow depths, increasing the risk of seismic hazard and the creation of fluid pathways.
The final section illustrates how the sealing capacity of fractures in OPA is affected by mineral composition. Triaxial flow-through experiments using Argon-gas were performed with dried samples from the sandy and shaly facies of OPA containing a roughened, artificial fracture. Slate, graywacke, quartzite, natural fault gouge, and granite samples were also tested to highlight the influence of normal stress, mineralogy and diagenesis on the sustainability of fracture transmissivity. With increasing normal stress, a non-linear decrease of fracture transmissivity can be observed that resulted in a permanent reduction of transmissivity after stress release. The transmissivity of rocks with a high portion of strong minerals (e.g., quartz) and high unconfined compressive strength was less sensitive to stress changes. In accordance with this, the sandy facies of OPA displayed a higher initial transmissivity that was less sensitive to stress changes compared to the shaly facies. However, transmissivity of rigid slate was less sensitive to stress changes than the sandy facies of OPA, although the slate is characterized by a higher phyllosilicate content. This demonstrates that in addition to mineral composition, other factors such as the degree of metamorphism, cementation and consolidation have to be considered when evaluating the sealing capacity of phyllosilicate-rich rocks.
The results of this thesis highlighted the role of confining pressure on the failure behavior of intact and artificially fractured OPA. Although the quartz-rich sandy facies may be considered as being more favorable for underground constructions due to its higher shear strength and stiffness than the shaly facies, the results indicate that when fractures develop in the sandy facies, they are more conductive and remain more permeable compared to fractures in the clay-dominated shaly facies at a given stress. The results may provide the basis for constitutive models to predict the integrity and evolution of a future repository. Clearly, the influence of composition and consolidation, e.g., by geological burial and uplift, on the mechanical sealing behavior of OPA highlights the need for a detailed site-specific material characterization for a future repository.
Bayesian geomorphology
(2020)
The rapidly growing amount and diversity of data are confronting us more than ever with the need to make informed predictions under uncertainty. The adverse impacts of climate change and natural hazards also motivate our search for reliable predictions. The range of statistical techniques that geomorphologists use to tackle this challenge has been growing, but rarely involves Bayesian methods. Instead, many geomorphic models rely on estimated averages that largely miss out on the variability of form and process. Yet seemingly fixed estimates of channel heads, sediment rating curves or glacier equilibrium lines, for example, are all prone to uncertainties. Neighbouring scientific disciplines such as physics, hydrology or ecology have readily embraced Bayesian methods to fully capture and better explain such uncertainties, as the necessary computational tools have advanced greatly. The aim of this article is to introduce the Bayesian toolkit to scientists concerned with Earth surface processes and landforms, and to show how geomorphic models might benefit from probabilistic concepts. I briefly review the use of Bayesian reasoning in geomorphology, and outline the corresponding variants of regression and classification in several worked examples.
Resolving Earth's surface at the meter scale is essential for an improved understanding of the dynamics of mass-movement processes. In this study, we explore the applicability and potential of digital elevation models (DEMs) derived from stereophotogrammetry to detect debris-flow channels in the Quebrada del Toro in the northwestern Argentine Andes. Our analysis relies on a high-resolution (3 m) DEM created from SPOT-7 tri-stereo satellite data. We carefully validated DEM quality with ∼6,000 differential GPS points and identified optimal parameters for DEM generation in high-relief terrain. After multiple processing steps, we achieved an accuracy of 0.051 ± 1.915 m (1σ) using n = 3,139 control points with cm precision. Previous studies have used the drainage area and slope framework to identify topographic signatures of debris flows within a catchment. We built upon this and investigated individual river-channel segments using connected-component (CC) analysis on meter-scale topographic data. We define CC as segments of similar slope along the channel profile. Based on seven manually identified debris-flow catchments, we developed a debris-flow similarity index using component length and mean channel-segment slope and identified channel segments that have likely been shaped by debris flows. The presented approach has the potential to resolve intra-catchment variability of transport processes, allows to constrain the extent of debris-flow channels more precisely than slope-area analysis, and highlights the versatility of combined space- and field-based observations for natural-hazard assessments.
Spectral analysis is a technique of time-series analysis that decomposes signals into linear combinations of harmonic components. Rooted in the 19th century, spectral analysis gained popularity in palaeoclimatology since the early 1980s. This was partly due to the availability of long time series of past climates, but also the development of new, partly adapted methods and the increasing spread of affordable personal computers. This paper reviews the most important methods of spectral analysis for palaeoclimate time series and discusses the prerequisites for their application as well as advantages and disadvantages. The paper also offers an overview of suitable software, as well as computer code for using the methods on synthetic examples.
The 1911 Chon-Kemin (Kebin) earthquake culminated c. 30 years of remarkable earthquakes in the northern Tien Shan (Kyrgyzstan and Kazakhstan). Building on prior mapping of the event, we traced its rupture in the field and measured more than 50 offset landforms. Cumulative fault rupture length is >155-195 km along 13 fault patches comprising six sections. The patches are separated by changes of dip magnitude or dip direction, or by 4-10 km-wide stepovers. One <40 km section overlaps and is parallel to the main north-dipping rupture but is 7 km north and dips opposite (south). Both ends of the rupture are along mountain front thrust faults demonstrating late Quaternary activity. We computed the moment from each fault patch using the surface fault traces, dip inferred from the traces, 20 km seismogenic thickness, rigidity of 3.3 x 10(10) N m(-2) and dip slip converted from our observations of the largely reverse sense of motion vertical offsets. The discontinuous patches with c. 3-4 m average slip and peak slip of <14 m yield a seismic moment of 4.6 x 10(20) Nm (M-w 7.78) to 7.4 x 10(20) Nm (M-w 7.91). The majority of moment was released along the inner eastern rupture segments. This geological moment is lower by a factor of 1.5 from that determined from teleseismic data.
The Lower Rhine Graben (Central Europe) is a prime example of a seismically active low-strain rift zone characterized by pronounced anthropogenic and climatic overprint of structures, and long recurrence intervals of large earthquakes. These factors render the identification of active faults and surface ruptures difficult. We investigated two fault scarps in the Lower Rhine Graben, to decipher their structural character, offset and potential seismogenic origin. Both scarps were modified by anthropogenic activity. The Hemmerich site lies c. 20 km SW of Cologne, along the Erft Fault. The Untermaubach site lies SW of Duren, where the Schafberg Fault projects into the Rur River valley. At the Hemmerich site, geomorphic and geophysical data, as well as exploratory coring reveal evidence of repeated normal faulting. Geophysical analysis and palaeoseismological excavation at the Untermaubach site reveal a complex fault zone in Holocene gravels characterized by subtle gravel deformation. Differentiation of tectonic and fluvial features was only possible with trenching, because fault structures and grain sizes of the sediments were below the resolution of the geophysical data. Despite these issues, our investigation demonstrates that valuable insight into past earthquakes and seismogenic deformation in a low-strain environment can be revealed using a multidisciplinary approach.
A series of large-magnitude earthquakes above 6.9 occurred in the northern Tien-Shan between 1885 and 1911. The Chilik earthquake of 11 July 1889, has been listed with a magnitude of 8.3, based on sparse macroseismic intensities, constrained by reported damage. Despite the existence of several juvenile fault scarps in the epicentral region, that are possibly associated with the 1889 earthquake, no through-going surface rupture having the dimensions expected for a magnitude 8.3 earthquake has been located - a puzzling dilemma. Could the magnitude have been overestimated? This would have major implications not only for the understanding of the earthquake series, but also for regional hazard estimates. Fortunately, a fragmentary record from an early Rebeur-Paschwitz seismometer exists for the Chilik event, recorded in Wilhelmshaven (Germany). To constrain the magnitude, we compare the late coda waves of this record with those of recent events from Central Asia, recorded on modern instruments in Germany and filtered with Rebeur-Paschwitz instrument characteristics. Additional constraints come from disturbances of historic magnetograms that exist from the Chilik and the 1911 Chon-Kemin earthquakes. Scaling of these historic records confirm a magnitude of about 8 for the 1889 Chilik earthquake, pointing towards a lower crustal contribution to the fault area.
Bayesian geomorphology
(2020)
The rapidly growing amount and diversity of data are confronting us more than ever with the need to make informed predictions under uncertainty. The adverse impacts of climate change and natural hazards also motivate our search for reliable predictions. The range of statistical techniques that geomorphologists use to tackle this challenge has been growing, but rarely involves Bayesian methods. Instead, many geomorphic models rely on estimated averages that largely miss out on the variability of form and process. Yet seemingly fixed estimates of channel heads, sediment rating curves or glacier equilibrium lines, for example, are all prone to uncertainties. Neighbouring scientific disciplines such as physics, hydrology or ecology have readily embraced Bayesian methods to fully capture and better explain such uncertainties, as the necessary computational tools have advanced greatly. The aim of this article is to introduce the Bayesian toolkit to scientists concerned with Earth surface processes and landforms, and to show how geomorphic models might benefit from probabilistic concepts. I briefly review the use of Bayesian reasoning in geomorphology, and outline the corresponding variants of regression and classification in several worked examples.
Cosmic-ray neutron sensing (CRNS) has become an effective method to measure soil moisture at a horizontal scale of hundreds of metres and a depth of decimetres. Recent studies proposed operating CRNS in a network with overlapping footprints in order to cover root-zone water dynamics at the small catchment scale and, at the same time, to represent spatial heterogeneity. In a joint field campaign from September to November 2020 (JFC-2020), five German research institutions deployed 15 CRNS sensors in the 0.4 km2 Wüstebach catchment (Eifel mountains, Germany). The catchment is dominantly forested (but includes a substantial fraction of open vegetation) and features a topographically distinct catchment boundary. In addition to the dense CRNS coverage, the campaign featured a unique combination of additional instruments and techniques: hydro-gravimetry (to detect water storage dynamics also below the root zone); ground-based and, for the first time, airborne CRNS roving; an extensive wireless soil sensor network, supplemented by manual measurements; and six weighable lysimeters. Together with comprehensive data from the long-term local research infrastructure, the published data set (available at https://doi.org/10.23728/b2share.756ca0485800474e9dc7f5949c63b872; Heistermann et al., 2022) will be a valuable asset in various research contexts: to advance the retrieval of landscape water storage from CRNS, wireless soil sensor networks, or hydrogravimetry; to identify scale-specific combinations of sensors and methods to represent soil moisture variability; to improve the understanding and simulation of land–atmosphere exchange as well as hydrological and hydrogeological processes at the hillslope and the catchment scale; and to support the retrieval of soil water content from airborne and spaceborne remote sensing platforms.
Cosmic-ray neutron sensing (CRNS) has become an effective method to measure soil moisture at a horizontal scale of hundreds of metres and a depth of decimetres. Recent studies proposed operating CRNS in a network with overlapping footprints in order to cover root-zone water dynamics at the small catchment scale and, at the same time, to represent spatial heterogeneity. In a joint field campaign from September to November 2020 (JFC-2020), five German research institutions deployed 15 CRNS sensors in the 0.4 km2 Wüstebach catchment (Eifel mountains, Germany). The catchment is dominantly forested (but includes a substantial fraction of open vegetation) and features a topographically distinct catchment boundary. In addition to the dense CRNS coverage, the campaign featured a unique combination of additional instruments and techniques: hydro-gravimetry (to detect water storage dynamics also below the root zone); ground-based and, for the first time, airborne CRNS roving; an extensive wireless soil sensor network, supplemented by manual measurements; and six weighable lysimeters. Together with comprehensive data from the long-term local research infrastructure, the published data set (available at https://doi.org/10.23728/b2share.756ca0485800474e9dc7f5949c63b872; Heistermann et al., 2022) will be a valuable asset in various research contexts: to advance the retrieval of landscape water storage from CRNS, wireless soil sensor networks, or hydrogravimetry; to identify scale-specific combinations of sensors and methods to represent soil moisture variability; to improve the understanding and simulation of land–atmosphere exchange as well as hydrological and hydrogeological processes at the hillslope and the catchment scale; and to support the retrieval of soil water content from airborne and spaceborne remote sensing platforms.
Early onset and late acceleration of rapid exhumation in the Namche Barwa syntaxis, eastern Himalaya
(2020)
The Himalayan syntaxes, characterized by extreme rates of rock exhumation co-located with major trans-orogenic rivers, figure prominently in the debate on tectonic versus erosional forcing of exhumation. Both the mechanism and timing of rapid exhumation of the Namche Barwa massif in the eastern syntaxis remain controversial. It has been argued that coupling between crustal rock advection and surface erosion initiated in the late Miocene (8-10 Ma). Recent studies, in contrast, suggest a Quaternary onset of rapid exhumation linked to a purely tectonic mechanism. We report new multisystem detrital thermochronology data from the most proximal Neogene clastic sediments downstream of Namche Barwa and use a thermo-kinematic model constrained by new and published data to explore its exhumation history. Modeling results show that exhumation accelerated to similar to 4 km/m.y. at ca. 8 Ma and to similar to 9 km/m.y. after ca. 2 Ma. This three-stage history reconciles apparently contradictory evidence for early and late onset of rapid exhumation and suggests efficient coupling between tectonics and erosion since the late Miocene. Quaternary acceleration of exhumation is consistent with river-profile evolution and may be linked to a Quaternary river-capture event.
Python is used in a wide range of geoscientific applications, such as in processing images for remote sensing, in generating and processing digital elevation models, and in analyzing time series. This book introduces methods of data analysis in the geosciences using Python that include basic statistics for univariate, bivariate, and multivariate data sets, time series analysis, and signal processing; the analysis of spatial and directional data; and image analysis. The text includes numerous examples that demonstrate how Python can be used on data sets from the earth sciences. The supplementary electronic material (available online through Springer Link) contains the example data as well as recipes that include all the Python commands featured in the book.
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.
The southern margin of the Central Anatolian Plateau (CAP) records a strong uplift phase after the early Middle Pleistocene, which has been related to the slab break-off of the subducting Arabian plate beneath the Anatolian microplate. During the last 450 kyr the area underwent an uplift phase at a mean rate of similar to 3.2 m/kyr, as suggested by Middle Pleistocene marine sediments exposed at similar to 1,500 m above sea level. These values are significantly higher than the 1.0-1.5 m/kyr estimated since the Late Pleistocene, suggesting temporal variations in uplift rate. To estimate changes in uplift rate during the Pleistocene we studied the marine terraces along the CAP southern margin, mapping the remnants of the platforms and their associated deposits in the field, and used the TerraceM software to identify the position and elevation of associated shoreline angles. We used shoreline angles and the timing of Quaternary marine sedimentation as constrains for a Landscape Evolution Model that simulates wave erosion of an uplifting coast. We applied random optimization algorithms and minimization statistics to find the input parameters that better reproduce the morphology of CAP marine terraces. The best-fitting uplift rate history suggests a significative increase from 1.9 to 3.5 m/kyr between 500 and 200 kyr, followed by an abrupt decrease to 1.4 m/kyr until the present. Our results agree with slab break-off models, which suggest a strong uplift pulse during slab rupture followed by a smoother decrease.
The Arctic is changing rapidly and permafrost is thawing. Especially ice-rich permafrost, such as the late Pleistocene Yedoma, is vulnerable to rapid and deep thaw processes such as surface subsidence after the melting of ground ice. Due to permafrost thaw, the permafrost carbon pool is becoming increasingly accessible to microbes, leading to increased greenhouse gas emissions, which enhances the climate warming.
The assessment of the molecular structure and biodegradability of permafrost organic matter (OM) is highly needed. My research revolves around the question “how does permafrost thaw affect its OM storage?” More specifically, I assessed (1) how molecular biomarkers can be applied to characterize permafrost OM, (2) greenhouse gas production rates from thawing permafrost, and (3) the quality of OM of frozen and (previously) thawed sediments.
I studied deep (max. 55 m) Yedoma and thawed Yedoma permafrost sediments from Yakutia (Sakha Republic). I analyzed sediment cores taken below thermokarst lakes on the Bykovsky Peninsula (southeast of the Lena Delta) and in the Yukechi Alas (Central Yakutia), and headwall samples from the permafrost cliff Sobo-Sise (Lena Delta) and the retrogressive thaw slump Batagay (Yana Uplands). I measured biomarker concentrations of all sediment samples. Furthermore, I carried out incubation experiments to quantify greenhouse gas production in thawing permafrost.
I showed that the biomarker proxies are useful to assess the source of the OM and to distinguish between OM derived from terrestrial higher plants, aquatic plants and microbial activity. In addition, I showed that some proxies help to assess the degree of degradation of permafrost OM, especially when combined with sedimentological data in a multi-proxy approach. The OM of Yedoma is generally better preserved than that of thawed Yedoma sediments. The greenhouse gas production was highest in the permafrost sediments that thawed for the first time, meaning that the frozen Yedoma sediments contained most labile OM. Furthermore, I showed that the methanogenic communities had established in the recently thawed sediments, but not yet in the still-frozen sediments.
My research provided the first molecular biomarker distributions and organic carbon turnover data as well as insights in the state and processes in deep frozen and thawed Yedoma sediments. These findings show the relevance of studying OM in deep permafrost sediments.
Applying conservation of energy to estimate earthquake frequencies from strain rates and stresses
(2020)
Estimating earthquake occurrence rates from the accumulation rate of seismic moment is an established tool of seismic hazard analysis. We propose an alternative, fault-agnostic approach based on the conservation of energy: the Energy-Conserving Seismicity Framework (ENCOS). Working in energy space has the advantage that the radiated energy is a better predictor of the damage potential of earthquake waves than the seismic moment release. In a region, ENCOS balances the stationary power available to cause earthquakes with the long-term seismic energy release represented by the energy-frequency distribution's first moment. Accumulation and release are connected through the average seismic efficiency, by which we mean the fraction of released energy that is converted into seismic waves. Besides measuring earthquakes in energy, ENCOS differs from moment balance essentially in that the energy accumulation rate depends on the total stress in addition to the strain rate tensor. To validate ENCOS, we exemplarily model the energy-frequency distribution around Southern California. We estimate the energy accumulation rate due to tectonic loading assuming poroelasticity and hydrostasis. Using data from the World Stress Map and assuming the frictional limit to estimate the stress tensor, we obtain a power of 0.8 GW. The uncertainty range, 0.3-2.0GW, originates mainly from the thickness of the seismogenic crust, the friction coefficient on preexisting faults, and models of Global Positioning System (GPS) derived strain rates. Based on a Gutenberg-Richter magnitude-frequency distribution, this power can be distributed over a range of energies consistent with historical earthquake rates and reasonable bounds on the seismic efficiency.
Temperature changes and variations in pore fluid salinity may negatively affect the permeability of clay-bearing sandstones with implications for natural fluid flow and geotechnical applications alike. In this study these factors are investigated for a sandstone dominated by illite as the clay phase. Systematic long-term flow-through experiments were conducted and complemented with comprehensive microstructural investigations and the application of Derjaguin-Landau-Verwey-Overbeek (DLVO) theory to explain mechanistically the observed permeability changes. Initially, sample permeability was not affected by low pore fluid salinity indicating strong attraction of the illite particles to the pore walls as supported by electron microprobe analysis (EMPA). Increasing temperature up to 145 degrees C resulted in an irreversible permeability decrease by 1.5 orders of magnitude regardless of the pore fluid composition (i.e., deionized water and 2 M NaCl solution). Subsequently diluting the high salinity pore fluid to below 0.5 M yielded an additional permeability decline by 1.5 orders of magnitude, both at 145 degrees C and after cooling to room temperature. By applying scanning electron microscopy (SEM) and mercury intrusion porosimetry (MIP) thermo-mechanical pore throat closure and illite particle migration were identified as independently operating mechanisms responsible for observed permeability changes during heating and dilution, respectively. These observations indicate that permeability of illite-bearing sandstones will be impaired by heating and exposure to low salinity pore fluids. In addition, chemically induced permeability variations proved to be path dependent with respect to the applied succession of fluid salinity changes.
In the field of hydrological modeling, many alternative representations of natural processes exist. Choosing specific process formulations when building a hydrological model is therefore associated with a high degree of ambiguity and subjectivity. In addition, the numerical integration of the underlying differential equations and parametrization of model structures influence model performance. Identifiability analysis may provide guidance by constraining the a priori range of alternatives based on observations. In this work, a flexible simulation environment is used to build an ensemble of semidistributed, process-based hydrological model configurations with alternative process representations, numerical integration schemes, and model parametrizations in an integrated manner. The flexible simulation environment is coupled with an approach for dynamic identifiability analysis. The objective is to investigate the applicability of the framework to identify the most adequate model. While an optimal model configuration could not be clearly distinguished, interesting results were obtained when relating model identifiability with hydro-meteorological boundary conditions. For instance, we tested the Penman-Monteith and Shuttleworth & Wallace evapotranspiration models and found that the former performs better under wet and the latter under dry conditions. Parametrization of model structures plays a dominant role as it can compensate for inadequate process representations and poor numerical solvers. Therefore, it was found that numerical solvers of high order of accuracy do often, though not necessarily, lead to better model performance. The proposed coupled framework proved to be a straightforward diagnostic tool for model building and hypotheses testing and shows potential for more in-depth analysis of process implementations and catchment functioning.
The Quaternary volcanic fields of the Eifel (Rhineland-Palatinate, Germany) had their last eruptions less than 13,000 years ago. Recently, deep low-frequency (DLF) earthquakes were detected beneath one of the volcanic fields showing evidence of ongoing magmatic activity in the lower crust and upper mantle. In this work, seismic wide- and steep-angle experiments from 1978/1979 and 1987/1988 are compiled, partially reprocessed and interpreted, together with other data to better determine the location, size, shape, and state of magmatic reservoirs in the Eifel region near the crust-mantle boundary. We discuss seismic evidence for a low-velocity gradient layer from 30-36 km depth, which has developed over a large region under all Quaternary volcanic fields of the Rhenish Massif and can be explained by the presence of partial melts. We show that the DLF earthquakes connect the postulated upper mantle reservoir with the upper crust at a depth of about 8 km, directly below one of the youngest phonolitic volcanic centers in the Eifel, where CO(2)originating from the mantle is massively outgassing. A bright spot in the West Eifel between 6 and 10 km depth represents a Tertiary magma reservoir and is seen as a model for a differentiated reservoir beneath the young phonolitic center today. We find that the distribution of volcanic fields is controlled by the Variscan lithospheric structures and terrane boundaries as a whole, which is reflected by an offset of the Moho depth, a wedge-shaped transparent zone in the lower crust and the system of thrusts over about 120 km length.
It has recently been demonstrated that the interaction of a mantle plume with sufficiently old oceanic lithosphere can initiate subduction. However, the existence of large lithospheric heterogeneities, such as a buoyant plateau, in proximity to a rising plume head may potentially hinder the formation of a new subduction zone. Here, we investigate this scenario by means of 3-D numerical thermomechanical modeling. We explore how plume-lithosphere interaction is affected by lithospheric age, relative location of plume head and plateau border, and the strength of the oceanic crust. Our numerical experiments suggest four different geodynamic regimes: (a) oceanic trench formation, (b) circular oceanic-plateau trench formation, (c) plateau trench formation, and (d) no trench formation. We show that regardless of the age and crustal strength of the oceanic lithosphere, subduction can initiate when the plume head is either below the plateau border or at a distance less than the plume radius from the plateau edge. Crustal heterogeneity facilitates subduction initiation of old oceanic lithosphere. High crustal strength hampers the formation of a new subduction zone when the plume head is located below a young lithosphere containing a thick and strong plateau. We suggest that plume-plateau interaction in the western margin of the Caribbean could have resulted in subduction initiation when the plume head impinged onto the oceanic lithosphere close to the border between plateau and oceanic crust.
Flood loss modeling is a central component of flood risk analysis. Conventionally, this involves univariable and deterministic stage-damage functions. Recent advancements in the field promote the use of multivariable and probabilistic loss models, which consider variables beyond inundation depth and account for prediction uncertainty. Although companies contribute significantly to total loss figures, novel modeling approaches for companies are lacking. Scarce data and the heterogeneity among companies impede the development of company flood loss models. We present three multivariable flood loss models for companies from the manufacturing, commercial, financial, and service sector that intrinsically quantify prediction uncertainty. Based on object-level loss data (n = 1,306), we comparatively evaluate the predictive capacity of Bayesian networks, Bayesian regression, and random forest in relation to deterministic and probabilistic stage-damage functions, serving as benchmarks. The company loss data stem from four postevent surveys in Germany between 2002 and 2013 and include information on flood intensity, company characteristics, emergency response, private precaution, and resulting loss to building, equipment, and goods and stock. We find that the multivariable probabilistic models successfully identify and reproduce essential relationships of flood damage processes in the data. The assessment of model skill focuses on the precision of the probabilistic predictions and reveals that the candidate models outperform the stage-damage functions, while differences among the proposed models are negligible. Although the combination of multivariable and probabilistic loss estimation improves predictive accuracy over the entire data set, wide predictive distributions stress the necessity for the quantification of uncertainty.
Giant earthquakes with magnitudes above 8.5 occur only in subduction zones. Despite the developments made in observing large subduction zone earthquakes with geophysical instruments, the factors controlling the maximum size of these earthquakes are still poorly understood. Previous studies have suggested the importance of slab shape, roughness of the plate interface contact, state of the strain in the upper plate, thickness of sediments filling the trenches, and subduction rate. Here, we present 2-D cross-scale numerical models of seismic cycles for subduction zones with various geometries, subduction channel friction configurations, and subduction rates. We found that low-angle subduction and thick sediments in the subduction channel are the necessary conditions for generating giant earthquakes, while the subduction rate has a negligible effect. We suggest that these key parameters determine the maximum magnitude of a subduction earthquake by controlling the seismogenic zone width and smoothness of the subduction interface. This interpretation supports previous studies that are based upon observations and scaling laws. Our modeling results also suggest that low static friction in the sediment-filled subduction channel results in neutral or moderate compressive deformation in the overriding plate for low-angle subduction zones hosting giant earthquakes. These modeling results agree well with observations for the largest earthquakes. Based on our models we predict maximum magnitudes of subduction earthquakes worldwide, demonstrating the fit to magnitudes of all giant earthquakes of the 20th and 21st centuries and good agreement with the predictions based on statistical analyses of observations.
Birds migrating across the Himalayan region fly over the highest peaks in the world, facing immense physiological and climatic challenges. The authors show the different strategies used by birds to cope with these challenges. Many wetland avian species are seen in the high-altitude lakes of the Himalayas and the adjoining Tibetan Plateau, such as Bar-Headed Geese. Ringing programmes have generated information about origins and destinations, and this book is the first to present information on the bird's exact migratory paths. Capitalising on knowledge generated through satellite telemetry, the authors describe the migratory routes of a multitude of birds flying over or skirting the Himalayas. The myriad of threats to migratory birds and the wetland system in the Central Asian Flyway are discussed, with ways to mitigate them. This volume will inform and persuade policy-makers and conservation practitioners to take appropriate measures for the long-term survival of this unique migration
The Coulomb failure stress (CFS) criterion is the most commonly used method for predicting spatial distributions of aftershocks following large earthquakes. However, large uncertainties are always associated with the calculation of Coulomb stress change. The uncertainties mainly arise due to nonunique slip inversions and unknown receiver faults; especially for the latter, results are highly dependent on the choice of the assumed receiver mechanism. Based on binary tests (aftershocks yes/no), recent studies suggest that alternative stress quantities, a distance-slip probabilistic model as well as deep neural network (DNN) approaches, all are superior to CFS with predefined receiver mechanism. To challenge this conclusion, which might have large implications, we use 289 slip inversions from SRCMOD database to calculate more realistic CFS values for a layered half-space and variable receiver mechanisms. We also analyze the effect of the magnitude cutoff, grid size variation, and aftershock duration to verify the use of receiver operating characteristic (ROC) analysis for the ranking of stress metrics. The observations suggest that introducing a layered half-space does not improve the stress maps and ROC curves. However, results significantly improve for larger aftershocks and shorter time periods but without changing the ranking. We also go beyond binary testing and apply alternative statistics to test the ability to estimate aftershock numbers, which confirm that simple stress metrics perform better than the classic Coulomb failure stress calculations and are also better than the distance-slip probabilistic model.
The timing and mechanisms of the Cretaceous sea incursions into Central Asia are still poorly constrained. We provide a new chronostratigraphic framework based on biostratigraphy and magnetostratigraphy together with detailed paleoenvironmental analyses of Cretaceous records of the proto-Paratethys Sea fluctuations in the Tajik and Tarim basins. The Early Cretaceous marine incursion in the western Tajik Basin was followed by major marine incursions during the Cenomanian (ca. 100 Ma) and Santonian (ca. 86 Ma) that reached far into the eastern Tajik and Tarim basins. These marine incursions were separated by a Turonian-Coniacian (ca. 92-86 Ma) regression. Basin-wide tectonic subsidence analyses imply that the Early Cretaceous sea incursion into the Tajik Basin was related to increased Pamir tectonism. We find that thrusting along the northern edge of the Pamir at ca. 130-90 Ma resulted in increased subsidence in a retro-arc basin setting. This tectonic event and coeval eustatic highstand resulted in the maximum observed geographic extent of the sea during the Cenomanian (ca. 100 Ma). The following Turonian-Coniacian (ca. 92-86 Ma) major regression, driven by eustasy, coincides with a sharp slowdown in tectonic subsidence during the late orogenic unloading period with limited thrusting. The Santonian (ca. 86 Ma) major sea incursion was likely controlled by eustasy as evidenced by the coeval fluctuations in the west Siberian Basin. An early Maastrichtian cooling (ca. 71-70 Ma), potentially connected to global Late Cretaceous trends, is inferred from the replacement of mollusk-rich limestones by bryozoan- and echinoderm-rich limestones.
In this paper, we present the new frequency spectrum recurrence analysis technique by means of electro-encephalon signals (EES) analyses. The technique is suitable for time series analysis with noise and disturbances. EES were collected, and alpha waves of the occipital region were analysed by comparing the signals from participants in two states, eyes open and eyes closed. Firstly, EES were characterized and analysed by means of techniques already known to compare with the results of the innovative technique that we present here. We verified that, standard recurrence quantification analysis by means of EES time series cannot statistically distinguish the two states. However, the new frequency spectrum recurrence quantification exhibit quantitatively whether the participants have their eyes open or closed. In sequence, new quantifiers are created for analysing the recurrence concentration on frequency bands. These analyses show that EES with similar frequency spectrum have different recurrence levels revealing different behaviours of the nervous system. The technique can be used to deepen the study on depression, stress, concentration level and other neurological issues and also can be used in any complex system.
The magnetic equator in the Brazilian region has moved over 1,100 km northward since 1957, passing the geomagnetic observatory Tatuoca (TTB), in northern Brazil, around 2013. We recovered and processed TTB hourly mean values of the geomagnetic field horizontal (H) component from 1957 until 2019, allowing the investigation of long-term changes in the daily variation due to the influence of secular variation, solar activity, season, and lunar phase. The H day-to-day variability and the occurrence of the counter electrojet at TTB were also investigated. Until the 1990s, ionospheric solar quiet currents dominated the quiet-time daily variation at TTB. After 2000, the magnitude of the daily variation became appreciably greater due to the equatorial electrojet (EEJ) contribution. The H seasonal and day-to-day variability increased as the magnetic equator approached, but their amplitudes normalized to the average daily variation remained at similar levels. Meanwhile, the amplitude of the lunar variation, normalized in the same way, increased from 5% to 12%. Within the EEJ region, the occurrence rate of the morning counter electrojet (MCEJ) increased with proximity to the magnetic equator, while the afternoon counter electrojet (ACEJ) did not. EEJ currents derived from CHAMP and Swarm satellite data revealed that the MCEJ rate varies with magnetic latitude within the EEJ region while the ACEJ rate is largely constant. Simulations with the Thermosphere-Ionosphere-Electrodynamics General Circulation Model based on different geomagnetic main field configurations suggest that long-term changes in the geomagnetic daily variation at TTB can be attributed to the main field secular variation.
The sediment profile from Lake Goscia(z) over dot in central Poland comprises a continuous, seasonally resolved and exceptionally well-preserved archive of the Younger Dryas (YD) climate variation. This provides a unique opportunity for detailed investigation of lake system responses during periods of rapid climate cooling (YD onset) and warming (YD termination). The new varve record of Lake Goscia(z) over dot presented here spans 1662 years from the late Allerod (AL) to the early Preboreal (PB). Microscopic varve counting provides an independent chronology with a YD duration of 1149+14/-22 years, which confirms previous results of 1140 +/- 40 years. We link stable oxygen isotopes and chironomid-based air temperature reconstructions with the response of various geochemical and varve microfacies proxies especially focusing on the onset and termination of the YD. Cooling at the YD onset lasted similar to 180 years, which is about a century longer than the terminal warming that was completed in similar to 70 years. During the AL/YD transition, environmental proxy data lagged the onset of cooling by similar to 90 years and revealed an increase of lake productivity and internal lake re-suspension as well as slightly higher detrital sediment input. In contrast, rapid warming and environmental changes during the YD/PB transition occurred simultaneously. However, initial changes such as declining diatom deposition and detrital input occurred already a few centuries before the rapid warming at the YD/PB transition. These environmental changes likely reflect a gradual increase in summer air temperatures already during the YD. Our data indicate complex and differing environmental responses to the major climate changes related to the YD, which involve different proxy sensitivities and threshold processes.
Li and B in ascending magmas: an experimental study on their mobility and isotopic fractionation
(2022)
This research study focuses on the behaviour of Li and B during magmatic ascent, and decompression-driven degassing related to volcanic systems. The main objective of this dissertation is to determine whether it is possible to use the diffusion properties of the two trace elements as a tool to trace magmatic ascent rate. With this objective, diffusion-couple and decompression experiments have been performed in order to study Li and B mobility in intra-melt conditions first, and then in an evolving system during decompression-driven degassing.
Synthetic glasses were prepared with rhyolitic composition and an initial water content of 4.2 wt%, and all the experiments were performed using an internally heated pressure vessel, in order to ensure a precise control on the experimental parameters such as temperature and pressure.
Diffusion-couple experiments were performed with a fix pressure 300 MPa. The temperature was varied in the range of 700-1250 °C with durations between 0 seconds and 24 hours. The diffusion-couple results show that Li diffusivity is very fast and starts already at very low temperature. Significant isotopic fractionation occurs due to the faster mobility of 6Li compared to 7Li. Boron diffusion is also accelerated by the presence of water, but the results of the isotopic ratios are unclear, and further investigation would be necessary to well constrain the isotopic fractionation process of boron in hydrous silicate melts. The isotopic ratios results show that boron isotopic fractionation might be affected by the speciation of boron in the silicate melt structure, as 10B and 11B tend to have tetrahedral and trigonal coordination, respectively.
Several decompression experiments were performed at 900 °C and 1000 °C, with pressures going from 300 MPa to 71-77 MPa and durations of 30 minutes, two, five and ten hours, in order to trigger water exsolution and the formation of vesicles in the sample. Textural observations and the calculation of the bubble number density confirmed that the bubble size and distribution after decompression is directly proportional to the decompression rate.
The overall SIMS results of Li and B show that the two trace elements tend to progressively decrease their concentration with decreasing decompression rates. This is explained because for longer decompression times, the diffusion of Li and B into the bubbles has more time to progress and the melt continuously loses volatiles as the bubbles expand their volumes.
For fast decompression, Li and B results show a concentration increase with a δ7Li and δ11B decrease close to the bubble interface, related to the sudden formation of the gas bubble, and the occurrence of a diffusion process in the opposite direction, from the bubble meniscus to the unaltered melt. When the bubble growth becomes dominant and Li and B start to exsolve into the gas phase, the silicate melt close to the bubble gets depleted in Li and B, because of a stronger diffusion of the trace elements into the bubble.
Our data are being applied to different models, aiming to combine the dynamics of bubble nucleation and growth with the evolution of trace elements concentration and isotopic ratios. Here, first considerations on these models will be presented, giving concluding remarks on this research study. All in all, the final remarks constitute a good starting point for further investigations. These results are a promising base to continue to study this process, and Li and B can indeed show clear dependences on decompression-related magma ascent rates in volcanic systems.
The southern Central Andes (SCA) (between 27 degrees S and 40 degrees S) is bordered to the west by the convergent margin between the continental South American Plate and the oceanic Nazca Plate. The subduction angle along this margin is variable, as is the deformation of the upper plate. Between 33 degrees S and 35 degrees S, the subduction angle of the Nazca plate increases from sub-horizontal (< 5 degrees) in the north to relatively steep (similar to 30 degrees) in the south. The SCA contain inherited lithological and structural heterogeneities within the crust that have been reactivated and overprinted since the onset of subduction and associated Cenozoic deformation within the Andean orogen. The distribution of the deformation within the SCA has often been attributed to the variations in the subduction angle and the reactivation of these inherited heterogeneities. However, the possible influence that the thickness and composition of the continental crust have had on both short-term and long-term deformation of the SCA is yet to be thoroughly investigated. For our investigations, we have derived density distributions and thicknesses for various layers that make up the lithosphere and evaluated their relationships with tectonic events that occurred over the history of the Andean orogeny and, in particular, investigated the short- and long-term nature of the present-day deformation processes. We established a 3D model of lithosphere beneath the orogen and its foreland (29 degrees S-39 degrees S) that is consistent with currently available geological and geophysical data, including the gravity data. The modelled crustal configuration and density distribution reveal spatial relationships with different tectonic domains: the crystalline crust in the orogen (the magmatic arc and the main orogenic wedge) is thicker (similar to 55 km) and less dense (similar to 2900 kg/m(3)) than in the forearc (similar to 35 km, similar to 2975 kg/m(3)) and foreland (similar to 30 km, similar to 3000 kg/m(3)). Crustal thickening in the orogen probably occurred as a result of stacking of low-density domains, while density and thickness variations beneath the forearc and foreland most likely reflect differences in the tectonic evolution of each area following crustal accretion. No clear spatial relationship exists between the density distribution within the lithosphere and previously proposed boundaries of crustal terranes accreted during the early Paleozoic. Areas with ongoing deformation show a spatial correlation with those areas that have the highest topographic gradients and where there are abrupt changes in the average crustal-density contrast. This suggests that the short-term deformation within the interior of the Andean orogen and its foreland is fundamentally influenced by the crustal composition and the relative thickness of different crustal layers. A thicker, denser, and potentially stronger lithosphere beneath the northern part of the SCA foreland is interpreted to have favoured a strong coupling between the Nazca and South American plates, facilitating the development of a sub-horizontal slab.
The Salt Range in Pakistan exposes Precambrian to Pleistocene strata outcropping along the Salt Range Thrust (SRT). To better understand the in-situ Cambrian and Pliocene tectonic evolution of the Pakistan Subhimalaya, we have conducted low-temperature thermochronological analysis using apatite (U-Th-Sm)/He and fission track dating. We combine cooling ages from different samples located along the thrust front of the SRT into a thermal model that shows two major cooling events associated with rifting and regional erosion in the Late Palaeozoic and SRT activity since the Pliocene. Our results suggest that the SRT maintained a long-term average shortening rate of similar to 5-6 mm/yr and a high exhumation rate above the SRT ramp since similar to 4 Ma.
AM(W)= 5.1 earthquake on January 21st, 2016 marked the beginning of a significant seismic sequence in the southern Alboran Sea, culminating in aM(W)= 6.3 earthquake on January 25th, and continuing with further moderate magnitude earthquakes until March. We use data from 35 seismic broadband stations in Spain, Morocco and Portugal to relocate the seismicity, estimate seismic moment tensors, and isolate regional apparent source time functions for the main earthquake. Relocation and regional moment tensor inversion consistently yield very shallow depths for the majority of events. We obtain 50 moment tensors for the sequence, showing a mixture of strike-slip faulting for the foreshock and the main event and reverse faulting for the major aftershocks. The leading role of reverse focal mechanisms among the aftershocks may be explained by the geometry of the fault network. The mainshock nucleates at a bend along the left-lateral Al-Idrisi fault, introducing local transpression within the transtensional Alboran Basin. The shallow depths of the 2016 Alboran Sea earthquakes may favor slip-partitioning on the involved faults. Apparent source durations for the main event suggest a similar to 21 km long, asymmetric rupture that propagates primarily toward NE into the restraining fault segment, with fast rupture speed of similar to 3.0 km/s. Consistently, the inversion for laterally variable fault displacement situates the main slip in the restraining segment. The partitioning into strike-slip rupture and dip-slip aftershocks confirms a non-optimal orientation of this segment, and suggests that the 2016 event settled a slip deficit from previous ruptures that could not propagate into the stronger restraining segment.
Due to the major role of greenhouse gas emissions in global climate change, the development of non-fossil energy technologies is essential. Deep geothermal energy represents such an alternative, which offers promising properties such as a high base load capability and a large untapped potential. The present work addresses barite precipitation within geothermal systems and the associated reduction in rock permeability, which is a major obstacle to maintaining high efficiency. In this context, hydro-geochemical models are essential to quantify and predict the effects of precipitation on the efficiency of a system.
The objective of the present work is to quantify the induced injectivity loss using numerical and analytical reactive transport simulations. For the calculations, the fractured-porous reservoirs of the German geothermal regions North German Basin (NGB) and Upper Rhine Graben (URG) are considered.
Similar depth-dependent precipitation potentials could be determined for both investigated regions (2.8-20.2 g/m3 fluid). However, the reservoir simulations indicate that the injectivity loss due to barite deposition in the NGB is significant (1.8%-6.4% per year) and the longevity of the system is affected as a result; this is especially true for deeper reservoirs (3000 m). In contrast, simulations of URG sites indicate a minor role of barite (< 0.1%-1.2% injectivity loss per year). The key differences between the investigated regions are reservoir thicknesses and the presence of fractures in the rock, as well as the ionic strength of the fluids. The URG generally has fractured-porous reservoirs with much higher thicknesses, resulting in a greater distribution of precipitates in the subsurface. Furthermore, ionic strengths are higher in the NGB, which accelerates barite precipitation, causing it to occur more concentrated around the wellbore. The more concentrated the precipitates occur around the wellbore, the higher the injectivity loss.
In this work, a workflow was developed within which numerical and analytical models can be used to estimate and quantify the risk of barite precipitation within the reservoir of geothermal systems. A key element is a newly developed analytical scaling score that provides a reliable estimate of induced injectivity loss. The key advantage of the presented approach compared to fully coupled reservoir simulations is its simplicity, which makes it more accessible to plant operators and decision makers. Thus, in particular, the scaling score can find wide application within geothermal energy, e.g., in the search for potential plant sites and the estimation of long-term efficiency.
As the Arctic coast erodes, it drains thermokarst lakes, transforming them into lagoons, and, eventually, integrates them into subsea permafrost. Lagoons represent the first stage of a thermokarst lake transition to a marine setting and possibly more saline and colder upper boundary conditions. In this research, borehole data, electrical resistivity surveying, and modeling of heat and salt diffusion were carried out at Polar Fox Lagoon on the Bykovsky Peninsula, Siberia. Polar Fox Lagoon is a seasonally isolated water body connected to Tiksi Bay through a channel, leading to hypersaline waters under the ice cover. The boreholes in the center of the lagoon revealed floating ice and a saline cryotic bed underlain by a saline cryotic talik, a thin ice-bearing permafrost layer, and unfrozen ground. The bathymetry showed that most of the lagoon had bedfast ice in spring. In bedfast ice areas, the electrical resistivity profiles suggested that an unfrozen saline layer was underlain by a thick layer of refrozen talik. The modeling showed that thermokarst lake taliks can refreeze when submerged in saltwater with mean annual bottom water temperatures below or slightly above 0 degrees C. This occurs, because the top-down chemical degradation of newly formed ice-bearing permafrost is slower than the refreezing of the talik. Hence, lagoons may precondition taliks with a layer of ice-bearing permafrost before encroachment by the sea, and this frozen layer may act as a cap on gas migration out of the underlying talik.
Reported magnitudes ranging between 7.8 and 8.7 highlight a confusion about the true size of the 1920 Haiyuan earthquake, the largest earthquake recorded in the northeast Tibetan Plateau. We compiled a global data set of previously unlooked-at historical seismograms and performed modern computational analyses on the digitized seismic records. We found the epicenter to be near Haiyuan town and obtained a moment magnitude of M-W=7.90.2. Following traditional approaches, we obtained m(B)=7.90.3 with data from 21 stations and M-S(20)=8.10.2 with data from three stations. Geomorphologically, we mapped the surface rupture and horizontal offsets on high-resolution Pleiades satellite and drone imagery that covered the entire western and middle sections of the 1920 Haiyuan earthquake rupture and compiled offsets reported on the eastern section from field measurements in the 1980s. Careful discrimination between single-event and cumulative offsets suggests average horizontal slips of 3.01.0m on the western section, 4.51.5m on the middle section, and 3.5 +/- 0.5m on the eastern section, indicating a total moment magnitude of M-W=7.8 +/- 0.1. Thus, the seismological and geomorphological results agree within the uncertainties, a weighted average giving a moment magnitude of M-W=7.9 +/- 0.2 for the 1920 Haiyuan earthquake. It is likely that earthquake magnitudes based on the historical M were systematically overestimated. <br /> Plain Language Summary Earthquakes are the main mechanism by which elastic energy accumulating due to tectonic motion is released. As the earthquake magnitude scale is logarithmic, major earthquakes control the bulk of this energy budget and are often the most destructive, like the 1920 Haiyuan earthquake with similar to 230,000 casualties. However, major earthquakes tend to have recurrence periods of several hundred years, longer than our instrumental records. To obtain knowledge of historic major earthquakes, paleoseismologists measure geomorphic offsets and map surface ruptures left by past events and estimate the shaking intensity from historical writings. However, in the case of the Haiyuan earthquake, which happened in the late historic, early instrumental period, the magnitudes reported from these two communities differed significantly. In order to constrain the magnitude of this earthquake for seismic hazard assessment and to reconcile the differences between published magnitudes, we reestimated its magnitude from both newly compiled and digitized seismological records and modern satellite and drone imagery. The results show that the early seismological magnitudes were overestimated, which may affect historical magnitudes systematically. The 1920 Haiyuan earthquake was of a similar magnitude to the 2001 Kokoxili and 2008 Wenchuan earthquakes that also occurred in and around the Tibetan Plateau, instead of more than half a magnitude larger.
The Brazilian Earth System Model (BESM-OA2.5), while simulating the historical period proposed by the fifth phase of the Coupled Model Intercomparison Project (CMIP5), detects an increasing trend in the sea surface height (SSH) on the southern hemisphere oceans relative to that of the pre-industrial era. The increasing trend is accentuated in the CMIP5 RCP4.5 and RCP8.5 future scenarios with higher concentrations of greenhouse gases in the atmosphere. This study sheds light on the sources of such trends in these regions. The results suggest an association with the thermal expansion of the oceans in the upper 700 m due to a gradual warming inflicted by those future scenarios. BESM-OA2.5 presents a surface height increase of 0.11 m in the historical period of 1850-2005. Concerning future projections, BESM-OA2.5 projects SSH increases of 0.14 and 0.23 m (relative to the historical 2005 value) for RCP4.5 and RCP8.5, respectively, by the end of 2100. These increases are predominantly in a band of latitude within 35-60 degrees S in the Atlantic and Indian oceans. The reproducibility of the trend signal detected in the BESM-OA2.5 simulations is confirmed by the results of three other CMIP5 models.
Forecasting and early warning systems are important investments to protect lives, properties, and livelihood. While early warning systems are frequently used to predict the magnitude, location, and timing of potentially damaging events, these systems rarely provide impact estimates, such as the expected amount and distribution of physical damage, human consequences, disruption of services, or financial loss. Complementing early warning systems with impact forecasts has a twofold advantage: It would provide decision makers with richer information to take informed decisions about emergency measures and focus the attention of different disciplines on a common target. This would allow capitalizing on synergies between different disciplines and boosting the development of multihazard early warning systems. This review discusses the state of the art in impact forecasting for a wide range of natural hazards. We outline the added value of impact-based warnings compared to hazard forecasting for the emergency phase, indicate challenges and pitfalls, and synthesize the review results across hazard types most relevant for Europe.
This review provides a synthesis of current knowledge on the morphological and functional traits of testate amoebae, a polyphyletic group of protists commonly used as proxies of past hydrological changes in paleoecological investigations from peatland, lake sediment and soil archives. A trait-based approach to understanding testate amoebae ecology and paleoecology has gained in popularity in recent years, with research showing that morphological characteristics provide complementary information to the commonly used environmental inferences based on testate amoeba (morpho-)species data. We provide a broad overview of testate amoeba morphological and functional traits and trait-environment relationships in the context of ecology, evolution, genetics, biogeography, and paleoecology. As examples we report upon previous ecological and paleoecological studies that used trait-based approaches, and describe key testate amoebae traits that can be used to improve the interpretation of environmental studies. We also highlight knowledge gaps and speculate on potential future directions for the application of trait-based approaches in testate amoeba research.
The origins and development of the arid and highly seasonal steppe-desert biome in Central Asia, the largest of its kind in the world, remain largely unconstrained by existing records. It is unclear how Cenozoic climatic, geological, and biological forces, acting at diverse spatial and temporal scales, shaped Central Asian ecosystems through time. Our synthesis shows that the Central Asian steppe-desert has existed since at least Eocene times but experienced no less than two regime shifts, one at the Eocene-Oligocene Transition and one in the mid-Miocene. These shifts separated three successive "stable states," each characterized by unique floral and faunal structures. Past responses to disturbance in the Asian steppe-desert imply that modern ecosystems are unlikely to recover their present structures and diversity if forced into a new regime. This is of concern for Asian steppes today, which are being modified for human use and lost to desertification at unprecedented rates.
Over the past few decades, azimuthal seismic anisotropy measurements have been widely used proxy to study past and present-day deformation of the lithosphere and to characterize convection in the mantle. Beneath continental regions, distinguishing between shallow and deep sources of anisotropy remains difficult due to poor depth constraints of measurements and a lack of regional-scale geodynamic modeling. Here, we constrain the sources of seismic anisotropy beneath Madagascar where a complex pattern cannot be explained by a single process such as absolute plate motion, global mantle flow, or geology. We test the hypotheses that either Edge-Driven Convection (EDC) or mantle flow derived from mantle wind interactions with lithospheric topography is the dominant source of anisotropy beneath Madagascar. We, therefore, simulate two sets of mantle convection models using regional-scale 3-D computational modeling. We then calculate Lattice Preferred Orientation that develops along pathlines of the mantle flow models and use them to calculate synthetic splitting parameters. Comparison of predicted with observed seismic anisotropy shows a good fit in northern and southern Madagascar for the EDC model, but the mantle wind case only fits well in northern Madagascar. This result suggests the dominant control of the measured anisotropy may be from EDC, but the role of localized fossil anisotropy in narrow shear zones cannot be ruled out in southern Madagascar. Our results suggest that the asthenosphere beneath northern and southern Madagascar is dominated by dislocation creep. Dislocation creep rheology may be dominant in the upper asthenosphere beneath other regions of continental lithosphere.
The Arctic is greatly affected by climate change. Increasing air temperatures drive permafrost thaw and an increase in coastal erosion and river discharge. This results in a greater input of sediment and organic matter into nearshore waters, impacting ecosystems by reducing light transmission through the water column and altering biogeochemistry. This potentially results in impacts on the subsistence economy of local people as well as the climate due to the transformation of suspended organic matter into greenhouse gases. Even though the impacts of increased suspended sediment concentrations and turbidity in the Arctic nearshore zone are well-studied, the mechanisms underpinning this increase are largely unknown. Wave energy and tides drive the level of turbidity in the temperate and tropical parts of the world, and this is generally assumed to also be the case in the Arctic. However, the tidal range is considerably lower in the Arctic, and processes related to the occurrence of permafrost have the potential to greatly contribute to nearshore turbidity. In this study, we use high-resolution satellite imagery alongside in situ and ERA5 reanalysis data of ocean and climate variables in order to identify the drivers of nearshore turbidity, along with its seasonality in the nearshore waters of Herschel Island Qikiqtaruk, in the western Canadian Arctic. Nearshore turbidity correlates well to wind direction, wind speed, significant wave height, and wave period. Nearshore turbidity is superiorly correlated to wind speed at the Beaufort Shelf compared to in situ measurements at Herschel Island Qikiqtaruk, showing that nearshore turbidity, albeit being of limited spatial extent, is influenced by large-scale weather and ocean phenomenons. We show that, in contrast to the temperate and tropical ocean, freshly eroded material is the predominant driver of nearshore turbidity in the Arctic, rather than resuspension, which is caused by the vulnerability of permafrost coasts to thermo-erosion.
The Arctic is greatly affected by climate change. Increasing air temperatures drive permafrost thaw and an increase in coastal erosion and river discharge. This results in a greater input of sediment and organic matter into nearshore waters, impacting ecosystems by reducing light transmission through the water column and altering biogeochemistry. This potentially results in impacts on the subsistence economy of local people as well as the climate due to the transformation of suspended organic matter into greenhouse gases. Even though the impacts of increased suspended sediment concentrations and turbidity in the Arctic nearshore zone are well-studied, the mechanisms underpinning this increase are largely unknown. Wave energy and tides drive the level of turbidity in the temperate and tropical parts of the world, and this is generally assumed to also be the case in the Arctic. However, the tidal range is considerably lower in the Arctic, and processes related to the occurrence of permafrost have the potential to greatly contribute to nearshore turbidity. In this study, we use high-resolution satellite imagery alongside in situ and ERA5 reanalysis data of ocean and climate variables in order to identify the drivers of nearshore turbidity, along with its seasonality in the nearshore waters of Herschel Island Qikiqtaruk, in the western Canadian Arctic. Nearshore turbidity correlates well to wind direction, wind speed, significant wave height, and wave period. Nearshore turbidity is superiorly correlated to wind speed at the Beaufort Shelf compared to in situ measurements at Herschel Island Qikiqtaruk, showing that nearshore turbidity, albeit being of limited spatial extent, is influenced by large-scale weather and ocean phenomenons. We show that, in contrast to the temperate and tropical ocean, freshly eroded material is the predominant driver of nearshore turbidity in the Arctic, rather than resuspension, which is caused by the vulnerability of permafrost coasts to thermo-erosion.
The Bajgan Complex, one of the basement constituents of the arc massif in Iranian Makran forms a rugged, deeply incised terrain. The complex consists of pelitic schists with minor psammitic and basic schists, calc silicate rocks, amphibolites, marbles, metavolcanosediments, mafic and felsic intrusives as well as ultramafic rocks. Metapelitic rocks show an amphibolite facies regional metamorphism and contain garnet, biotite, white mica, quartz, albite ± rutile ± apatite. Thermobarometry of garnet schist yields pressure of more than 9 kbar and temperatures between 560 and 675 °C. The geothermal gradient obtained for the peak of regional metamorphism is 19 °C/km, corresponding to a depth of ca. 31 km. Replacement of garnet by chlorite and epidote suggest greenschist facies metamorphism due to a decrease in temperature and pressure through exhumation and retrograde metamorphism (370–450 °C and 3–6 kbar). The metapelitic rocks followed a ‘clockwise’ P–T path during metamorphism, consistent with thermal decline following tectonic thickening. The formation of medium-pressure metamorphic rocks is related to presence of active subduction of the Neotethys Oceanic lithosphere beneath Eurasia in the Makran.
Mixed sand- and gravel-bed rivers record erosion, transport, and fining signals in their bedload size distributions. Thus, grain-size data are imperative for studying these processes. However, collecting hundreds to thousands of pebble measurements in steep and dynamic high-mountain river settings remains challenging. Using the recently published digital grain-sizing algorithm PebbleCounts, we were able to survey seven large (>= 1,000 m2) channel cross-sections and measure thousands to tens-of-thousands of grains per survey along a 100-km stretch of the trunk stream of the Toro Basin in Northwest Argentina. The study region traverses a steep topographic and environmental gradient on the eastern margin of the Central Andean Plateau. Careful counting and validation allows us to identify measurement errors and constrain percentile uncertainties using large sample sizes. In the coarse >= 2.5 cm fraction of bedload, only the uppermost size percentiles (>= 95th) vary significantly downstream, whereas the 50th and 84th percentiles show less variability. We note a relation between increases in these upper percentiles and along-channel junctions with large, steep tributaries. This signal is strongly influenced by lithology and geologic structures, and mixed with local hillslope input. In steep catchments like the Toro Basin, we suggest nonlinear relationships between geomorphic metrics and grain size, whereby the steepest parts of the landscape exert primary control on the upper grain-size percentiles. Thus, average or median metrics that do not apply weights or thresholds to steeper topography may be less predictive of grain-size distributions in such settings.
During the second phase of the Alpine Fault, Deep Fault Drilling Project (DFDP) in the Whataroa River, South Westland, New Zealand, bedrock was encountered in the DFDP-2B borehole from 238.5–893.2 m Measured Depth (MD). Continuous sampling and meso- to microscale characterisation of whole rock cuttings established that, in sequence, the borehole sampled amphibolite facies, Torlesse Composite Terrane-derived schists, protomylonites and mylonites, terminating 200–400 m above an Alpine Fault Principal Slip Zone (PSZ) with a maximum dip of 62°. The most diagnostic structural features of increasing PSZ proximity were the occurrence of shear bands and reduction in mean quartz grain sizes. A change in composition to greater mica:quartz + feldspar, most markedly below c. 700 m MD, is inferred to result from either heterogeneous sampling or a change in lithology related to alteration. Major oxide variations suggest the fault-proximal Alpine Fault alteration zone, as previously defined in DFDP-1 core, was not sampled.
Water management tools are necessary to guarantee the preservation of natural resources while ensuring optimum utilization. Linear regression models are a simple and quick solution for creating prognostic capabilities. Multivariate models show higher precision than univariate models. In the case of Waiwera, implementation of individual production rates is more accurate than applying just the total production rate. A maximum of approximately 1,075 m3/day can be pumped to ensure a water level of at least 0.5 m a.s.l. in the monitoring well. The model should be renewed annually to implement new data and current water level trends to keep the quality.
Predicting macroscopic elastic rock properties requires detailed information on microstructure
(2017)
Predicting variations in macroscopic mechanical rock behaviour due to microstructural changes, driven by mineral precipitation and dissolution is necessary to couple chemo-mechanical processes in geological subsurface simulations. We apply 3D numerical homogenization models to estimate Young’s moduli for five synthetic microstructures, and successfully validate our results for comparable geometries with the analytical Mori-Tanaka approach. Further, we demonstrate that considering specific rock microstructures is of paramount importance, since calculated elastic properties may deviate by up to 230 % for the same mineral composition. Moreover, agreement between simulated and experimentally determined Young’s moduli is significantly improved, when detailed spatial information are employed.
Integration and development of the energy supply in China and worldwide is a challenge for the years to come. The innovative idea presented here is based on an extension of the “power-to-gas-to-power” technology by establishing a closed carbon cycle. It is an implementation of a low-carbon energy system based on carbon dioxide capture and storage (CCS) to store and reuse wind and solar energy. The Chenjiacun storage project in China compares well with the German case study for the towns Potsdam and Brandenburg/Havel in the Federal State of Brandenburg based on the Ketzin pilot site for CCS.
Dissolved CO2 storage in geological formations with low pressure, low risk and large capacities
(2017)
Geological CO2 storage is a mitigation technology to reduce CO2 emissions from fossil fuel combustion. However, major concerns are the pressure increase and saltwater displacement in the mainly targeted deep groundwater aquifers due to injection of supercritical CO2. The suggested solution is storage of CO2 exclusively in the dissolved state. In our exemplary regional case study of the North East German Basin based on a highly resolved temperature and pressure distribution model and a newly developed reactive transport coupling, we have quantified that 4.7 Gt of CO2 can be stored in solution compared to 1.5 Gt in the supercritical state.
Tree stands in the boreal treeline ecotone are, in addition to climate change, impacted by disturbances such as fire, water-related disturbances and logging. We aim to understand how these disturbances affect growth, age structure, and spatial patterns of larch stands in the north-eastern Siberian treeline ecotone (lower Kolyma River region), an insufficiently researched region. Stand structure of Larix cajanderi Mayr was studied at seven sites impacted by disturbances. Maximum tree age ranged from 44 to 300 years. Young to medium-aged stands had, independent of disturbance type, the highest stand densities with over 4000 larch trees per ha. These sites also had the highest growth rates for tree height and stem diameter. Overall lowest stand densities were found in a polygonal field at the northern end of the study area, with larches growing in distinct " tree islands". At all sites, saplings are significantly clustered. Differences in fire severity led to contrasting stand structures with respect to tree, recruit, and overall stand densities. While a low severity fire resulted in low-density stands with high proportions of small and young larches, high severity fires resulted in high-density stands with high proportions of big trees. At waterdisturbed sites, stand structure varied between waterlogged and drained sites and latitude. These mixed effects of climate and disturbance make it difficult to predict future stand characteristics and the treeline position.
Different lake systems might reflect different climate elements of climate changes, while the responses of lake systems are also divers, and are not completely understood so far. Therefore, a comparison of lakes in different climate zones, during the high-amplitude and abrupt climate fluctuations of the Last Glacial to Holocene transition provides an exceptional opportunity to investigate distinct natural lake system responses to different abrupt climate changes. The aim of this doctoral thesis was to reconstruct climatic and environmental fluctuations down to (sub-) annual resolution from two different lake systems during the Last Glacial-Interglacial transition (~17 and 11 ka). Lake Gościąż, situated in the temperate central Poland, developed in the Allerød after recession of the Last Glacial ice sheets. The Dead Sea is located in the Levant (eastern Mediterranean) within a steep gradient from sub-humid to hyper-arid climate, and formed in the mid-Miocene. Despite their differences in sedimentation processes, both lakes form annual laminations (varves), which are crucial for studies of abrupt climate fluctuations. This doctoral thesis was carried out within the DFG project PALEX-II (Paleohydrology and Extreme Floods from the Dead Sea ICDP Core) that investigates extreme hydro-meteorological events in the ICDP core in relation to climate changes, and ICLEA (Virtual Institute of Integrated Climate and Landscape Evolution Analyses) that intends to better the understanding of climate dynamics and landscape evolutions in north-central Europe since the Last Glacial. Further, it contributes to the Helmholtz Climate Initiative REKLIM (Regional Climate Change and Humans) Research Theme 3 “Extreme events across temporal and spatial scales” that investigates extreme events using climate data, paleo-records and model-based simulations. The three main aims were to (1) establish robust chronologies of the lakes, (2) investigate how major and abrupt climate changes affect the lake systems, and (3) to compare the responses of the two varved lakes to these hemispheric-scale climate changes.
Robust chronologies are a prerequisite for high-resolved climate and environmental reconstructions, as well as for archive comparisons. Thus, addressing the first aim, the novel chronology of Lake Gościąż was established by microscopic varve counting and Bayesian age-depth modelling in Bacon for a non-varved section, and was corroborated by independent age constrains from 137Cs activity concentration measurements, AMS radiocarbon dating and pollen analysis. The varve chronology reaches from the late Allerød until AD 2015, revealing more Holocene varves than a previous study of Lake Gościąż suggested. Varve formation throughout the complete Younger Dryas (YD) even allowed the identification of annually- to decadal-resolved leads and lags in proxy responses at the YD transitions.
The lateglacial chronology of the Dead Sea (DS) was thus far mainly based on radiocarbon and U/Th-dating. In the unique ICDP core from the deep lake centre, continuous search for cryptotephra has been carried out in lateglacial sediments between two prominent gypsum deposits – the Upper and Additional Gypsum Units (UGU and AGU, respectively). Two cryptotephras were identified with glass analyses that correlate with tephra deposits from the Süphan and Nemrut volcanoes indicating that the AGU is ~1000 years younger than previously assumed, shifting it into the YD, and the underlying varved interval into the Bølling/Allerød, contradicting previous assumptions.
Using microfacies analyses, stable isotopes and temperature reconstructions, the second aim was achieved at Lake Gościąż. The YD lake system was dynamic, characterized by higher aquatic bioproductivity, more re-suspended material and less anoxia than during the Allerød and Early Holocene, mainly influenced by stronger water circulation and catchment erosion due to stronger westerly winds and less lake sheltering. Cooling at the YD onset was ~100 years longer than the final warming, while environmental proxies lagged the onset of cooling by ~90 years, but occurred contemporaneously during the termination of the YD. Chironomid-based temperature reconstructions support recent studies indicating mild YD summer temperatures. Such a comparison of annually-resolved proxy responses to both abrupt YD transitions is rare, because most European lake archives do not preserve varves during the YD.
To accomplish the second aim at the DS, microfacies analyses were performed between the UGU (~17 ka) and Holocene onset (~11 ka) in shallow- (Masada) and deep-water (ICDP core) environments. This time interval is marked by a huge but fluctuating lake level drop and therefore the complete transition into the Holocene is only recorded in the deep-basin ICDP core. In this thesis, this transition was investigated for the first time continuously and in detail. The final two pronounced lake level drops recorded by deposition of the UGU and AGU, were interrupted by one millennium of relative depositional stability and a positive water budget as recorded by aragonite varve deposition interrupted by only a few event layers. Further, intercalation of aragonite varves between the gypsum beds of the UGU and AGU shows that these generally dry intervals were also marked by decadal- to centennial-long rises in lake level. While continuous aragonite varves indicate decadal-long stable phases, the occurrence of thicker and more frequent event layers suggests general more instability during the gypsum units. These results suggest a pattern of complex and variable hydroclimate at different time scales during the Lateglacial at the DS.
The third aim was accomplished based on the individual studies above that jointly provide an integrated picture of different lake responses to different climate elements of hemispheric-scale abrupt climate changes during the Last Glacial-Interglacial transition. In general, climatically-driven facies changes are more dramatic in the DS than at Lake Gościąż. Further, Lake Gościąż is characterized by continuous varve formation nearly throughout the complete profile, whereas the DS record is widely characterized by extreme event layers, hampering the establishment of a continuous varve chronology. The lateglacial sedimentation in Lake Gościąż is mainly influenced by westerly winds and minor by changes in catchment vegetation, whereas the DS is primarily influenced by changes in winter precipitation, which are caused by temperature variations in the Mediterranean. Interestingly, sedimentation in both archives is more stable during the Bølling/Allerød and more dynamic during the YD, even when sedimentation processes are different.
In summary, this doctoral thesis presents seasonally-resolved records from two lake archives during the Lateglacial (ca 17-11 ka) to investigate the impact of abrupt climate changes in different lake systems. New age constrains from the identification of volcanic glass shards in the lateglacial sediments of the DS allowed the first lithology-based interpretation of the YD in the DS record and its comparison to Lake Gościąż. This highlights the importance of the construction of a robust chronology, and provides a first step for synchronization of the DS with other eastern Mediterranean archives. Further, climate reconstructions from the lake sediments showed variability on different time scales in the different archives, i.e. decadal- to millennial fluctuations in the lateglacial DS, and even annual variations and sub-decadal leads and lags in proxy responses during the rapid YD transitions in Lake Gościąż. This showed the importance of a comparison of different lake archives to better understand the regional and local impacts of hemispheric-scale climate variability. An unprecedented example is demonstrated here of how different lake systems show different lake responses and also react to different climate elements of abrupt climate changes. This further highlights the importance of the understanding of the respective lake system for climate reconstructions.
Gas hydrates are ice-like crystalline compounds made of water cavities that retain various types of guest molecules. Natural gas hydrates are CH4-rich but also contain higher hydrocarbons as well as CO2, H2S, etc. They are highly dependent of local pressure and temperature conditions. Considering the high energy content, natural gas hydrates are artificially dissociated for the production of methane gas. Besides, they may also dissociate in response to global warming. It is therefore crucial to investigate the hydrate nucleation and growth process at a molecular level. The understanding of how guest molecules in the hydrate cavities respond to warming climate or gas injection is also of great importance.
This thesis is concerned with a systematic investigation of simple and mixed gas hydrates at conditions relevant to the natural hydrate reservoir in Qilian Mountain permafrost, China. A high-pressure cell that integrated into the confocal Raman spectroscopy ensured a precise and continuous characterization of the hydrate phase during formation/dissociation/transformation processes with a high special and spectral resolution. By applying laboratory experiments, the formation of mixed gas hydrates containing other hydrocarbons besides methane was simulated in consideration of the effects from gas supply conditions and sediments. The results revealed a preferential enclathration of different guest molecules in hydrate cavities and further refute the common hypothesis of the coexistence of hydrate phases due to a changing feed gas phase. However, the presence of specific minerals and organic compounds in sediments may have significant impacts on the coexisting solid phases. With regard to the dissociation, the formation damage caused by fines mobilization and migration during hydrate decomposition was reported for the first time, illustrating the complex interactions between fine grains and hydrate particles. Gas hydrates, starting from simple CH4 hydrates to binary CH4—C3H8 hydrates and multi-component mixed hydrates were decomposed by thermal stimulation mimicking global warming. The mechanisms of guest substitution in hydrate structures were studied through the experimental data obtained from CH4—CO2, CH4—mixed gas hydrates and mixed gas hydrates—CO2 systems. For the first time, a second transformation behavior was documented during the transformation process from CH4 hydrates to CO2-rich mixed hydrates. Most of the crystals grew or maintained when exposed to CO2 gas while some others decreased in sizes and even disappeared over time. The highlight of the two last experimental simulations was to visualize and characterize the hydrate crystals which were at different structural transition stages. These experimental simulations enhanced our knowledge about the mixed gas hydrates in natural reservoirs and improved our capability to assess the response to global warming.
Technological progress allows for producing ever more complex predictive models on the basis of increasingly big datasets. For risk management of natural hazards, a multitude of models is needed as basis for decision-making, e.g. in the evaluation of observational data, for the prediction of hazard scenarios, or for statistical estimates of expected damage. The question arises, how modern modelling approaches like machine learning or data-mining can be meaningfully deployed in this thematic field. In addition, with respect to data availability and accessibility, the trend is towards open data. Topic of this thesis is therefore to investigate the possibilities and limitations of machine learning and open geospatial data in the field of flood risk modelling in the broad sense. As this overarching topic is broad in scope, individual relevant aspects are identified and inspected in detail.
A prominent data source in the flood context is satellite-based mapping of inundated areas, for example made openly available by the Copernicus service of the European Union. Great expectations are directed towards these products in scientific literature, both for acute support of relief forces during emergency response action, and for modelling via hydrodynamic models or for damage estimation. Therefore, a focus of this work was set on evaluating these flood masks. From the observation that the quality of these products is insufficient in forested and built-up areas, a procedure for subsequent improvement via machine learning was developed. This procedure is based on a classification algorithm that only requires training data from a particular class to be predicted, in this specific case data of flooded areas, but not of the negative class (dry areas). The application for hurricane Harvey in Houston shows the high potential of this method, which depends on the quality of the initial flood mask.
Next, it is investigated how much the predicted statistical risk from a process-based model chain is dependent on implemented physical process details. Thereby it is demonstrated what a risk study based on established models can deliver. Even for fluvial flooding, such model chains are already quite complex, though, and are hardly available for compound or cascading events comprising torrential rainfall, flash floods, and other processes. In the fourth chapter of this thesis it is therefore tested whether machine learning based on comprehensive damage data can offer a more direct path towards damage modelling, that avoids explicit conception of such a model chain. For that purpose, a state-collected dataset of damaged buildings from the severe El Niño event 2017 in Peru is used. In this context, the possibilities of data-mining for extracting process knowledge are explored as well. It can be shown that various openly available geodata sources contain useful information for flood hazard and damage modelling for complex events, e.g. satellite-based rainfall measurements, topographic and hydrographic information, mapped settlement areas, as well as indicators from spectral data. Further, insights on damaging processes are discovered, which mainly are in line with prior expectations. The maximum intensity of rainfall, for example, acts stronger in cities and steep canyons, while the sum of rain was found more informative in low-lying river catchments and forested areas. Rural areas of Peru exhibited higher vulnerability in the presented study compared to urban areas. However, the general limitations of the methods and the dependence on specific datasets and algorithms also become obvious.
In the overarching discussion, the different methods – process-based modelling, predictive machine learning, and data-mining – are evaluated with respect to the overall research questions. In the case of hazard observation it seems that a focus on novel algorithms makes sense for future research. In the subtopic of hazard modelling, especially for river floods, the improvement of physical models and the integration of process-based and statistical procedures is suggested. For damage modelling the large and representative datasets necessary for the broad application of machine learning are still lacking. Therefore, the improvement of the data basis in the field of damage is currently regarded as more important than the selection of algorithms.
Enhanced geothermal systems (EGS) are considered a cornerstone of future sustainable energy production. In such systems, high-pressure fluid injections break the rock to provide pathways for water to circulate in and heat up. This approach inherently induces small seismic events that, in rare cases, are felt or can even cause damage. Controlling and reducing the seismic impact of EGS is crucial for a broader public acceptance. To evaluate the applicability of hydraulic fracturing (HF) in EGS and to improve the understanding of fracturing processes and the hydromechanical relation to induced seismicity, six in-situ, meter-scale HF experiments with different injection schemes were performed under controlled conditions in crystalline rock in a depth of 410 m at the Äspö Hard Rock Laboratory (Sweden).
I developed a semi-automated, full-waveform-based detection, classification, and location workflow to extract and characterize the acoustic emission (AE) activity from the continuous recordings of 11 piezoelectric AE sensors. Based on the resulting catalog of 20,000 AEs, with rupture sizes of cm to dm, I mapped and characterized the fracture growth in great detail. The injection using a novel cyclic injection scheme (HF3) had a lower seismic impact than the conventional injections. HF3 induced fewer AEs with a reduced maximum magnitude and significantly larger b-values, implying a decreased number of large events relative to the number of small ones. Furthermore, HF3 showed an increased fracture complexity with multiple fractures or a fracture network. In contrast, the conventional injections developed single, planar fracture zones (Publication 1).
An independent, complementary approach based on a comparison of modeled and observed tilt exploits transient long-period signals recorded at the horizontal components of two broad-band seismometers a few tens of meters apart from the injections. It validated the efficient creation of hydraulic fractures and verified the AE-based fracture geometries. The innovative joint analysis of AEs and tilt signals revealed different phases of the fracturing process, including the (re-)opening, growth, and aftergrowth of fractures, and provided evidence for the reactivation of a preexisting fault in one of the experiments (Publication 2). A newly developed network-based waveform-similarity analysis applied to the massive AE activity supports the latter finding.
To validate whether the reduction of the seismic impact as observed for the cyclic injection schemes during the Äspö mine-scale experiments is transferable to other scales, I additionally calculated energy budgets for injection experiments from previously conducted laboratory tests and from a field application. Across all three scales, the cyclic injections reduce the seismic impact, as depicted by smaller maximum magnitudes, larger b-values, and decreased injection efficiencies (Publication 3).
Localisation of deformation is a ubiquitous feature in continental rift dynamics and observed across drastically different time and length scales. This thesis comprises one experimental and two numerical modelling studies investigating strain localisation in (1) a ductile shear zone induced by a material heterogeneity and (2) in an active continental rift setting. The studies are related by the fact that the weakening mechanisms on the crystallographic and grain size scale enable bulk rock weakening, which fundamentally enables the formation of shear zones, continental rifts and hence plate tectonics. Aiming to investigate the controlling mechanisms on initiation and evolution of a shear zone, the torsion experiments of the experimental study were conducted in a Patterson type apparatus with strong Carrara marble cylinders with a weak, planar Solnhofen limestone inclusion. Using state-of-the-art numerical modelling software, the torsion experiments were simulated to answer questions regarding localisation procedure like stress distribution or the impact of rheological weakening. 2D numerical models were also employed to integrate geophysical and geological data to explain characteristic tectonic evolution of the Southern and Central Kenya Rift. Key elements of the numerical tools are a randomized initial strain distribution and the usage of strain softening. During the torsion experiments, deformation begins to localise at the limestone inclusion tips in a process zone, which propagates into the marble matrix with increasing deformation until a ductile shear zone is established. Minor indicators for coexisting brittle deformation are found close to the inclusion tip and presumed to slightly facilitate strain localisation besides the dominant ductile deformation processes. The 2D numerical model of the torsion experiment successfully predicts local stress concentration and strain rate amplification ahead of the inclusion in first order agreement with the experimental results. A simple linear parametrization of strain weaking enables high accuracy reproduction of phenomenological aspects of the observed weakening. The torsion experiments suggest that loading conditions do not affect strain localisation during high temperature deformation of multiphase material with high viscosity contrasts. A numerical simulation can provide a way of analysing the process zone evolution virtually and extend the examinable frame. Furthermore, the nested structure and anastomosing shape of an ultramylonite band was mimicked with an additional second softening step. Rheological weakening is necessary to establish a shear zone in a strong matrix around a weak inclusion and for ultramylonite formation.
Such strain weakening laws are also incorporated into the numerical models of the
Southern and Central Kenya Rift that capture the characteristic tectonic evolution. A three-stage early rift evolution is suggested that starts with (1) the accommodation of strain by a single border fault and flexure of the hanging-wall crust, after which (2) faulting in the hanging-wall and the basin centre increases before (3) the early-stage asymmetry is lost and basinward localisation of deformation occurs. Along-strike variability of rifts can be produced by modifying the initial random noise distribution. In summary, the three studies address selected aspects of the broad range of mechanisms and processes that fundamentally enable the deformation of rock and govern the localisation patterns across the scales. In addition to the aforementioned results, the first and second manuscripts combined, demonstrate a procedure to find new or improve on existing numerical formulations for specific rheologies and their dynamic weakening. These formulations are essential in addressing rock deformation from the grain to the global scale. As within the third study of this thesis, where geodynamic controls on the evolution of a rift were examined and acquired by the integration of geological and geophysical data into a numerical model.
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.
We investigate source processes of fluid-induced seismicity from The Geysers geothermal reservoir in California to determine their relation with hydraulic operations and improve the corresponding seismic hazard estimates. Analysis of 869 well-constrained full moment tensors (M-w 0.8-3.5) reveals significant non-double-couple components (>25%) for about 65% of the events. Volumetric deformation is governed by cumulative injection rates with larger non-double-couple components observed near the wells and during high injection periods. Source mechanisms are magnitude dependent and vary significantly between faulting regimes. Normal faulting events (M-w<2) reveal substantial volumetric components indicating dilatancy in contrast to strike-slip events that have a dominant double-couple source. Volumetric components indicating closure of cracks in the source region are mostly found for reverse faulting events with M-w>2.5. Our results imply that source processes and magnitudes of fluid-induced seismic events are strongly affected by the hydraulic operations, the reservoir stress state, and the faulting regime.
Fragmentation of landscapes creates a transition zone in between natural habitats or different kinds of land use. In forested and agricultural landscapes with transition zones, microclimate and matter cycling are markedly altered. This probably accelerates and is intensified by global warming. However, there is no consensus on defining transition zones and quantifying relevant variables for microclimate and matter cycling across disciplines. This article is an attempt to a) revise definitions and offer a framework for quantitative ecologists, b) review the literature on microclimate and matter cycling in transition zones and c) summarise this information using meta-analysis to better understand bio-geochemical and bio-geophysical processes and their spatial extent in transition zones. We expect altered conditions in soils of transition zones to be 10-20 m with a maximum of 50 m, and 25-50 m for above-ground space with a maximum of 125 m.
We present a high-resolution study of the upper mantle structure of Central Europe, including the western part of the East European Platform, based on S-receiver functions of 345 stations. A distinct contrast is found between Phanerozoic Europe and the East European Craton across the Trans-European Suture Zone. To the west, a pronounced velocity reduction with depth interpreted as lithosphere-asthenosphere boundary (LAB) is found at an average depth of 90 km. Beneath the craton, no strong and continuous LAB conversion is observed. Instead we find a distinct velocity reduction within the lithosphere, at 80-120 km depth. This mid-lithospheric discontinuity (MLD) is attributed to a compositional boundary between depleted and more fertile lithosphere created by late Proterozoic metasomatism. A potential LAB phase beneath the craton is very weak and varies in depth between 180 and 250 km, consistent with a reduced velocity contrast between the lower lithosphere and the asthenosphere. Within the Trans-European Suture Zone, lithospheric structure is characterized by strong heterogeneity. A dipping or step-wise increase to LAB depth of 150 km is imaged from Phanerozoic Europe to 20-22 degrees E, whereas no direct connection to the cratonic LAB or MLD to the east is apparent. At larger depths, a positive conversion associated with the lower boundary of the asthenosphere is imaged at 210-250 km depth beneath Phanerozoic Europe, continuing down to 300 km depth beneath the craton. Conversions from both 410 km and 660 km discontinuities are found at their nominal depth beneath Phanerozoic Europe, and the discontinuity at 410 km depth can also be traced into the craton. A potential negative conversion on top of the 410 km discontinuity found in migrated images is analyzed by modeling and attributed to interference with other converted phases.
The Kenya rift revisited
(2017)
We present three-dimensional (3-D) models that describe the present-day thermal and rheological state of the lithosphere of the greater Kenya rift region aiming at a better understanding of the rift evolution, with a particular focus on plume-lithosphere interactions. The key methodology applied is the 3-D integration of diverse geological and geophysical observations using gravity modelling. Accordingly, the resulting lithospheric-scale 3-D density model is consistent with (i) reviewed descriptions of lithological variations in the sedimentary and volcanic cover, (ii) known trends in crust and mantle seismic velocities as revealed by seismic and seismological data and (iii) the observed gravity field. This data-based model is the first to image a 3-D density configuration of the crystalline crust for the entire region of Kenya and northern Tanzania. An upper and a basal crustal layer are differentiated, each composed of several domains of different average densities. We interpret these domains to trace back to the Precambrian terrane amalgamation associated with the East African Orogeny and to magmatic processes during Mesozoic and Cenozoic rifting phases. In combination with seismic velocities, the densities of these crustal domains indicate compositional differences. The derived lithological trends have been used to parameterise steady-state thermal and rheological models. These models indicate that crustal and mantle temperatures decrease from the Kenya rift in the west to eastern Kenya, while the integrated strength of the lithosphere increases. Thereby, the detailed strength configuration appears strongly controlled by the complex inherited crustal structure, which may have been decisive for the onset, localisation and propagation of rifting.
Moderate to large earthquakes can increase the amount of water feeding stream flows, mobilizing excess water from deep groundwater, shallow groundwater, or the vadose zone. Here we examine the regional pattern of streamflow response to the Maule M8.8 earthquake across Chile's diverse topographic and hydro-climatic gradients. We combine streamflow analyses with groundwater flow modeling and a random forest classifier, and find that, after the earthquake, at least 85 streams had a change in flow. Discharge mostly increased () shortly after the earthquake, liberating an excess water volume of >1.1 km3, which is the largest ever reported following an earthquake. Several catchments had increased discharge of >50 mm, locally exceeding seasonal streamflow discharge under undisturbed conditions. Our modeling results favor enhanced vertical permeability induced by dynamic strain as the most probable process explaining the observed changes at the regional scale. Supporting this interpretation, our random forest classification identifies peak ground velocity and elevation extremes as most important for predicting streamflow response. Given the mean recurrence interval of ∼25 yr for >M8.0 earthquakes along the Peru–Chile Trench, our observations highlight the role of earthquakes in the regional water cycle, especially in arid environments.
Temporal and spatial stability of the vegetation climate relationship is a basic ecological assumption for pollen-based quantitative inferences of past climate change and for predicting future vegetation. We explore this assumption for the Holocene in eastern continental Asia (China, Mongolia). Boosted regression trees (BRT) between fossil pollen taxa percentages (Abies, Artemisia, Betula, Chenopodiaceae, Cyperaceae, Ephedra, Picea, Pinus, Poaceae and Quercus) and climate model outputs of mean annual precipitation (P-ann) and mean temperature of the warmest month (Mt(wa)) for 9 and 6 ka (ka = thousand years before present) were set up and results compared to those obtained from relating modern pollen to modern climate. Overall, our results reveal only slight temporal differences in the pollen climate relationships. Our analyses suggest that the importance of P-ann compared with Mt(wa) for taxa distribution is higher today than it was at 6 ka and 9 ka. In particular, the relevance of P-ann for Picea and Pinus increases and has become the main determinant. This change in the climate tree pollen relationship parallels a widespread tree pollen decrease in north-central China and the eastern Tibetan Plateau. We assume that this is at least partly related to vegetation climate disequilibrium originating from human impact. Increased atmospheric CO2 concentration may have permitted the expansion of moisture-loving herb taxa (Cyperaceae and Poaceae) during the late Holocene into arid/semi-arid areas. We furthermore find that the pollen climate relationship between north-central China and the eastern Tibetan Plateau is generally similar, but that regional differences are larger than temporal differences. In summary, vegetation climate relationships in China are generally stable in space and time, and pollen-based climate reconstructions can be applied to the Holocene. Regional differences imply the calibration-set should be restricted spatially.
Landscape evolution models (LEMs) allow the study of earth surface responses to changing climatic and tectonic forcings. While much effort has been devoted to the development of LEMs that simulate a wide range of processes, the numerical accuracy of these models has received less attention. Most LEMs use first-order accurate numerical methods that suffer from substantial numerical diffusion. Numerical diffusion particularly affects the solution of the advection equation and thus the simulation of retreating landforms such as cliffs and river knickpoints. This has potential consequences for the integrated response of the simulated landscape. Here we test a higher-order flux-limiting finite volume method that is total variation diminishing (TVD-FVM) to solve the partial differential equations of river incision and tectonic displacement. We show that using the TVD-FVM to simulate river incision significantly influences the evolution of simulated landscapes and the spatial and temporal variability of catchment-wide erosion rates. Furthermore, a two-dimensional TVD-FVM accurately simulates the evolution of landscapes affected by lateral tectonic displacement, a process whose simulation was hitherto largely limited to LEMs with flexible spatial discretization. We implement the scheme in TTLEM (TopoToolbox Landscape Evolution Model), a spatially explicit, raster-based LEM for the study of fluvially eroding landscapes in TopoToolbox 2.
Along a subduction zone, great megathrust earthquakes recur either after long seismic gaps lasting several decades to centuries or over much shorter periods lasting hours to a few years when cascading successions of earthquakes rupture nearby segments of the fault. We analyze a decade of continuous Global Positioning System observations along the South American continent to estimate changes in deformation rates between the 2010 Maule (M8.8) and 2015 Illapel (M8.3) Chilean earthquakes. We find that surface velocities increased after the 2010 earthquake, in response to continental-scale viscoelastic mantle relaxation and to regional-scale increased degree of interplate locking. We propose that increased locking occurs transiently during a super-interseismic phase in segments adjacent to a megathrust rupture, responding to bending of both plates caused by coseismic slip and subsequent afterslip. Enhanced strain rates during a super-interseismic phase may therefore bring a megathrust segment closer to failure and possibly triggered the 2015 event.
Atmospheric transport is an understudied mechanism for leaf wax hydrogen isotope applications that contributes to mobilizing and depositing these compounds on the surface of the Earth. While previous efforts have identified the importance of atmospheric leaf wax deposition in remote marine locations, the processes are not well constrained on land in temperate latitudes where lakes are common and sedimentary leaf wax hydrogen isotope values are an attractive tool for understanding past precipitation changes. This work presents results from a field study that was conducted in 2010 and 2011 at Hainich National Park, Germany in order to evaluate the quantity and sources of leaf waxes in the atmosphere. Aerosols were sampled at approximately weekly intervals inside the forest canopy, and n-alkane distributions and hydrogen isotope values were compared with those from major tree species surrounding the sampling site. Despite sampling in what was expected to be a major production center, the distribution and hydrogen isotope values of atmospheric n-alkanes bore little resemblance to those of the local vegetation. Comparison with local meteorological data and to 10-day and 36-h back air mass trajectories indicated shifting effects of winds and temperature, and that mesoscale transport processes were more important than longrange mechanisms. Back trajectories also highlighted source effects, with easterly winds coinciding with relatively lower leaf wax hydrogen isotope values from more continental regions. These results suggest that leaf wax aerosols average over spatial scales that exceed typical surface catchment areas for small lake systems, even in forested areas, yet that the area over which these compounds are derived is still relatively regional. Depositional fluxes were also estimated in order to assess the potential importance of atmospheric transport to sedimentary archives. Although difficult to constrain, these estimates suggest that atmospheric deposition may be non-negligible for lake systems in cases where inputs from rivers or surface runoff are limited. Together, these observations provide new insights on how leaf waxes from different sources are integrated during aeolian transport and the spatial scales over which these processes occur.
Sensitivity and identifiability analyses are common diagnostic tools to address over-parametrization in complex environmental models, but a combined application of the two analyses is rarely conducted. In this study, we performed a temporal global sensitivity analysis using the variance-based method of Sobol’ and a temporal identifiability analysis of model parameters using the dynamic identifiability method (DYNIA). We discuss the relationship between the two analyses with a focus on parameter identification and output uncertainty reduction. The hydrological model HydroGeoSphere was used to simulate daily evapotranspiration, water content, and seepage at the lysimeter scale. We found that identifiability of a parameter does not necessarily reduce output uncertainty. It was also found that the information from the main and total effects (main Sobol' sensitivity indices) is required to allow uncertainty reduction in the model output. Overall, the study highlights the role of combined temporal diagnostic tools for improving our understanding of model behavior.
The present-day vegetation in the tropics is mainly characterized by forests worldwide except in tropical East Africa, where forests only occur as patches at the coast and in the uplands. These forest patches result from the peculiar aridity that is linked to the uplift of the region during the Late Cenozoic. The Late Cenozoic vegetation history of East Africa is of particular interest as it has set the scene for the contemporary events in mammal and hominin evolution. In this study, we investigate the conditions under which these forest patches could have been connected, and a previous continuous forest belt could have extended and fragmented. We apply a dynamic vegetation model with a set of climatic scenarios in which we systematically alter the present-day environmental conditions such that they would be more favourable for a continuous forest belt in tropical East Africa. We consider varying environmental factors, namely temperature, precipitation and atmospheric CO2 concentrations. Our results show that all of these variables play a significant role in supporting the forest biomes and a continuous forest belt could have occurred under certain combinations of these settings. With our current knowledge of the palaeoenvironmental history of East Africa, it is likely that the region hosted these conditions during the Late Cenozoic. Recent improvements on environmental hypotheses of hominin evolution highlight the role of periods of short and extreme climate variability during the Late Cenozoic specific to East Africa in driving evolution. Our results elucidate how the forest biomes of East Africa can appear and disappear under fluctuating environmental conditions and demonstrate how this climate variability might be recognized on the biosphere level.