Refine
Year of publication
Document Type
- Article (2483)
- Doctoral Thesis (471)
- 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 (3274) (remove)
Among the multitude of geomorphological processes, aeolian shaping processes are of special character, Pedogenic dust is one of the most important sources of atmospheric aerosols and therefore regarded as a key player for atmospheric processes. Soil dust emissions, being complex in composition and properties, influence atmospheric processes and air quality and has impacts on other ecosystems. In this because even though their immediate impact can be considered low (exceptions exist), their constant and large-scale force makes them a powerful player in the earth system. dissertation, we unravel a novel scientific understanding of this complex system based on a holistic dataset acquired during a series of field experiments on arable land in La Pampa, Argentina. The field experiments as well as the generated data provide information about topography, various soil parameters, the atmospheric dynamics in the very lower atmosphere (4m height) as well as measurements regarding aeolian particle movement across a wide range of particle size classes between 0.2μm up to the coarse sand.
The investigations focus on three topics: (a) the effects of low-scale landscape structures on aeolian transport processes of the coarse particle fraction, (b) the horizontal and vertical fluxes of the very fine particles and (c) the impact of wind gusts on particle emissions.
Among other considerations presented in this thesis, it could in particular be shown, that even though the small-scale topology does have a clear impact on erosion and deposition patterns, also physical soil parameters need to be taken into account for a robust statistical modelling of the latter. Furthermore, specifically the vertical fluxes of particulate matter have different characteristics for the particle size classes. Finally, a novel statistical measure was introduced to quantify the impact of wind gusts on the particle uptake and its application on the provided data set. The aforementioned measure shows significantly increased particle concentrations during points in time defined as gust event.
With its holistic approach, this thesis further contributes to the fundamental understanding of how atmosphere and pedosphere are intertwined and affect each other.
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.