Refine
Year of publication
- 2017 (298) (remove)
Document Type
- Article (229)
- Doctoral Thesis (28)
- Postprint (20)
- Other (11)
- Review (5)
- Part of a Book (2)
- Bachelor Thesis (1)
- Conference Proceeding (1)
- Master's Thesis (1)
Keywords
- Holocene (9)
- climate change (5)
- Arctic (3)
- Chile (3)
- Climate change (3)
- Erosion (3)
- Palaeoclimate (3)
- Tibetan Plateau (3)
- change detection (3)
- earthquake (3)
Institute
- Institut für Geowissenschaften (298) (remove)
Normally in paleoseismology, the study of the tectonic slip-rate is performed in trenches on the fault scarp, or by the estimation of fault movements from the geomorphic features. In this work, we have carried out a paleoseismic analysis of the Benis Fault, located in southeast Spain, combined with a geothermal analysis inside a deep cave related to the fault (-350 m). Thus, we have estimated the last earthquake magnitude and time of occurrence from evidence of ceiling collapse and displaced carbonate blocks inside a cave, which is developed across the fault. The magnitude was obtained from the application of the empirical relationship of the fault parameters and coseismic vertical displacement, yielding a value ranging between M 5.9 and M 6.5. Moreover, we dated this paleoearthquake by the paleontological record of a "Lynx pardinus spelaea", with an age of 65 +/- 18 ka BP. Additionally, we have measured the thermal profile of the Benis Cave (-350 m of depth), from single rock point temperature measurements during 2 years. The temperature profile shows three different parts inside the cave, the shallow heterogeneous thermal zone till 50 m depth; the homogeneous thermal zone 150 m till with constant temperature and the hetero-thermal deep zone, deeper than 150 m and till the deepest zone (350 m). Furthermore, we have estimated the Vertical Geothermal Gradient, 1.85 degrees C/100 m for the deepest zone (-150; -290 m). The temperature increases with depth, showing a reverse thermal profile in comparison with normal gradients in deep caves. Finally, we have calculated the heat flux of 0.46 mWm(2). (C) 2016 Elsevier Ltd and INQUA. All rights reserved.
Ecosystem services (ESs) are defined as the contributions that ecosystems make to human wellbeing and are increasingly being used as an approach to explore the importance of ecosystems for humans through their valuation. Although value plurality has been recognised long before the mainstreaming of ESs research, socio-cultural valuation is still underrepresented in ESs assessments. It is the central goal of this PhD dissertation to explore the ability of socio-cultural valuation methods for the operationalisation of ESs research in land management. To address this, I formulated three research objectives that are briefly outlined below and relate to the three studies conducted during this dissertation.
The first objective relates to the assessment of the current role of socio-cultural valuation in ESs research. Human values are central to ESs research yet non-monetary socio-cultural valuation methods have been found underrepresented in the field of ESs science. In regard to the unbalanced consideration of value domains and conceptual uncertainties, I perform a systematic literature review aiming to answer the research question: To what extent have socio-cultural values been addressed in ESs assessments.
The second objective aims to test socio-cultural valuation methods of ESs and their relevance for land use preferences by exploring their methodological opportunities and limitations. Socio-cultural valuation methods have only recently become a focus in ESs research and therefore bear various uncertainties in regard to their methodological implications. To overcome these uncertainties, I analysed responses to a visitor survey. The research questions related to the second objective were: What are the implications of different valuation methods for ESs values? To what extent are land use preferences explained by socio-cultural values of ESs?
The third objective addressed in this dissertation is the implementation of ESs research into land management through socio-cultural valuation. Though it is emphasised that the ESs approach can assist decision making, there is little empirical evidence of the effect of ESs knowledge on land management. I proposed a way to implement transdisciplinary, spatially explicit research on ESs by answering the following research questions: Which landscape features underpinning ESs supply are considered in land management? How can participatory approaches accounting for ESs be operationalised in land management?
The empirical research resulted in five main findings that provide answers to the research questions. First, this dissertation provides evidence that socio-cultural values are an integral part of ESs research. I found that they can be assessed for provisioning, regulating, and cultural services though they are linked to cultural services to a greater degree. Socio-cultural values have been assessed by monetary and non-monetary methods and their assessment is effectively facilitated by stakeholder participation. Second, I found that different methods of socio-cultural valuation revealed different information. Whereas rating revealed a general value of ESs, weighting was found more suitable to identify priorities across ESs. Value intentions likewise differed in the distribution of values, generally implying a higher value for others than for respondents themselves. Third, I showed that ESs values were distributed similarly across groups with differing land use preferences. Thus, I provided empirical evidence that ESs values and landscape values should not be used interchangeably. Fourth, I showed which landscape features important for ESs supply in a Scottish regional park are not sufficiently accounted for in the current management strategy. This knowledge is useful for the identification of priority sites for land management. Finally, I provide an approach to explore how ESs knowledge elicited by participatory mapping can be operationalised in land management. I demonstrate how stakeholder knowledge and values can be used for the identification of ESs hotspots and how these hotspots can be compared to current management priorities.
This dissertation helps to bridge current gaps of ESs science by advancing the understanding of the current role of socio-cultural values in ESs research, testing different methods and their relevance for land use preferences, and implementing ESs knowledge into land management. If and to what extent ESs and their values are implemented into ecosystem management is mainly the choice of the management. An advanced understanding of socio-cultural valuation methods contributes to the normative basis of this management, while the proposal for the implementation of ESs in land management presents a practical approach of how to transfer this type of knowledge into practice. The proposed methods for socio-cultural valuation can support guiding land management towards a balanced consideration of ESs and conservation goals.
Estuarine marshes are ecosystems that are situated at the transition zone between land and water and are thus controlled by physical and biological interactions. Marsh vegetation offers important ecosystem services by filtrating solid and dissolved substances from the water and providing habitat. By buffering a large part of the arriving flow velocity, attenuating wave energy and serving as erosion control for riverbanks, tidal marshes furthermore reduce the destructive effects of storm surges and storm waves and thus contribute to ecosystem-based shore protection. However, in many estuaries, extensive embankments, artificial bank protection, river dredging and agriculture threaten tidal marshes. Global warming might entail additional risks, such as changes in water levels, an increase of the tidal amplitude and a resulting shift of the salinity zones. This can affect the dynamics of the shore and foreland vegetation, and vegetation belts can be narrowed or fragmented. Against this background, it is crucial to gain a better understanding of the processes underlying the spatio temporal vegetation dynamics in brackish marshes. Furthermore, a better understanding of how plant-habitat relationships generate patterns in tidal marsh vegetation is vital to maintain ecosystem functions and assess the response of marshes to environmental change as well as the success of engineering and restoration projects.
For this purpose, three research objectives were addressed within this thesis: (1) to explore the possibility of vegetation serving as self-adaptive shore protection by quantifying the reduction of current velocity in the vegetation belt and the morphologic plasticity of a brackish marsh pioneer, (2) to disentangle the roles of abiotic factors and interspecific competition on species distribution and stand characteristics in brackish marshes, and (3) to develop a mechanistic vegetation model that helps analysing the influence of habitat conditions on the spatio-temporal dynamic of tidal marsh vegetation. These aspects were investigated using a combination of field studies and statistical as well as process-based modelling.
To explore the possibility of vegetation serving as self-adaptive coastal protection, in the first study, we measured current velocity with and without living vegetation, recorded ramet density and plant thickness during two growing periods at two locations in the Elbe estuary and assessed the adaptive value of a larger stem diameter of plants at locations with higher mechanical stress by biomechanical measurements. The results of this study show that under non-storm conditions, the vegetation belt of the marsh pioneer Bolboschoenus maritimus is able to buffer a large proportion of the flow velocity. We were furthermore able to show that morphological traits of plant species are adapted to hydrodynamic forces by demonstrating a positive correlation between ramet thickness and cross-shore current. In addition, our measurements revealed that thicker ramets growing at the front of the vegetation belt have a significantly higher stability than ramets inside the vegetation belt. This self-adaptive effect improves the ability of B. maritimus to grow and persist in the pioneer zone and could provide an adaptive value in habitats with high mechanical stress.
In the second study, we assessed the distribution of the two marsh species and a set of stand characteristics, namely aboveground and belowground biomass, ramet density, ramet height and the percentage of flowering ramets. Furthermore, we collected information on several abiotic habitat factors to test their effect on plant growth and zonation with generalised linear models (GLMs). Our results demonstrate that flow velocity is the main factor controlling the distribution of Bolboschoenus maritimus and Phragmites australis. Additionally, inundation height and duration, as well as intraspecific competition affect distribution patterns. This study furthermore shows that cross-shore flow velocity does not only directly influence the distribution of the two marsh species, but also alters the plants’ occurrence relative to inun-dation height and duration. This suggests an effect of cross-shore flow velocity on their tolerance to inundation. The analysis of the measured stand characteristics revealed a negative effect of total flow velocity on all measured parameters of B. maritimus and thus confirmed our expectation that flow velocity is a decisive stressor which influences the growth of this species.
To gain a better understanding of the processes and habitat factors influencing the spatio-temporal vegetation dynamics in brackish marshes, I built a spatially explicit, mechanistic model applying a pattern-oriented modelling approach. A sensitivity analysis of the para-meters of this dynamic habitat-macrophyte model HaMac suggests that rhizome growth is the key process for the lateral dynamics of brackish marshes. From the analysed habitat factors, P. australis patterns were mainly influenced by flow velocity. The competition with P. australis was of key importance for the belowground biomass of B. maritimus. Concerning vegetation dynamics, the model results emphasise that without the effect of flow velocity the B. maritimus vegetation belt would expand into the tidal flat at locations with present vegetation recession, suggesting that flow velocity is the main reason for vegetation recession at exposed locations.
Overall, the results of this thesis demonstrate that brackish marsh vegetation considerably contributes to flow reduction under average flow conditions and can hence be a valuable component of shore-protection schemes. At the same time, the distribution, growth and expansion of tidal marsh vegetation is substantially influenced by flow. Altogether, this thesis provides a clear step forward in understanding plant-habitat interactions in tidal marshes. Future research should integrate studies of vertical marsh accretion with research on the factors that control the lateral position of marshes.
In general, a moderate drying trend is observed in mid-latitude arid Central Asia since the Mid-Holocene, attributed to the progressively weakening influence of the mid-latitude Westerlies on regional climate. However, as the spatio-temporal pattern of this development and the underlying climatic mechanisms are yet not fully understood, new high-resolution paleoclimate records from this region are needed. Within this study, a sediment core from Lake Son Kol (Central Kyrgyzstan) was investigated using sedimentological, (bio) geochemical, isotopic, and palynological analyses, aiming at reconstructing regional climate development during the last 6000 years. Biogeochemical data, mainly reflecting summer moisture conditions, indicate predominantly wet conditions until 4950 cal. yr BP, succeeded by a pronounced dry interval between 4950 and 3900 cal. yr BP. In the following, a return to wet conditions and a subsequent moderate drying trend until present times are observed. This is consistent with other regional paleoclimate records and likely reflects the gradual Late Holocene diminishment of the amount of summer moisture provided by the mid-latitude Westerlies. However, climate impact of the Westerlies was apparently not only restricted to the summer season but also significant during winter as indicated by recurrent episodes of enhanced allochthonous input through snowmelt, occurring before 6000 cal. yr BP and at 5100-4350, 3450-2850, and 1900-1500 cal. yr BP. The distinct similar to 1500year periodicity of these episodes of increased winter precipitation in Central Kyrgyzstan resembles similar cyclicities observed in paleoclimate records around the North Atlantic, likely indicating a hemispheric-scale climatic teleconnection and an impact of North Atlantic Oscillation (NAO) variability in Central Asia.
A new sedimentary sequence from Lago di Venere on Pantelleria Island, located in the Strait of Sicily between Tunisia and Sicily was recovered. The lake is located in the coastal infra-Mediterranean vegetation belt at 2 m a.s.l. Pollen, charcoal and sedimentological analyses are used to explore linkages among vegetation, fire and climate at a decadal scale over the past 1200 years. A dry period from ad 800 to 1000 that corresponds to the Medieval Warm Period' (WMP) is inferred from sedimentological analysis. The high content of carbonate recorded in this period suggests a dry phase, when the ratio of evaporation/precipitation was high. During this period the island was dominated by thermophilous and drought-tolerant taxa, such as Quercus ilex, Olea, Pistacia and Juniperus. A marked shift in the sediment properties is recorded at ad 1000, when carbonate content became very low suggesting wetter conditions until ad 1850-1900. Broadly, this period coincides with the Little Ice Age' (LIA), which was characterized by wetter and colder conditions in Europe. During this time rather mesic conifers (i.e. Pinus pinaster), shrubs and herbs (e.g. Erica arborea and Selaginella denticulata) expanded, whereas more drought-adapted species (e.g. Q. ilex) declined. Charcoal data suggest enhanced fire activity during the LIA probably as a consequence of anthropogenic burning and/or more flammable fuel (e.g. resinous Pinus biomass). The last century was characterized by a shift to high carbonate content, indicating a change towards drier conditions, and re-expansion of Q. ilex and Olea. The post-LIA warming is in agreement with historical documents and meteorological time series. Vegetation dynamics were co-determined by agricultural activities on the island. Anthropogenic indicators (e.g. Cerealia-type, Sporormiella) reveal the importance of crops and grazing on the island. Our pollen data suggest that extensive logging caused the local extinction of deciduous Quercus pubescens around ad1750.
The central aim of this thesis is to demonstrate the benefits of innovative frequency-based methods to better explain the variability observed in lake ecosystems. Freshwater ecosystems may be the most threatened part of the hydrosphere. Lake ecosystems are particularly sensitive to changes in climate and land use because they integrate disturbances across their entire catchment. This makes understanding the dynamics of lake ecosystems an intriguing and important research priority. This thesis adds new findings to the baseline knowledge regarding variability in lake ecosystems. It provides a literature-based, data-driven and methodological framework for the investigation of variability and patterns in environmental parameters in the time frequency domain.
Observational data often show considerable variability in the environmental parameters of lake ecosystems. This variability is mostly driven by a plethora of periodic and stochastic processes inside and outside the ecosystems. These run in parallel and may operate at vastly different time scales, ranging from seconds to decades. In measured data, all of these signals are superimposed, and dominant processes may obscure the signals of other processes, particularly when analyzing mean values over long time scales. Dominant signals are often caused by phenomena at long time scales like seasonal cycles, and most of these are well understood in the limnological literature. The variability injected by biological, chemical and physical processes operating at smaller time scales is less well understood. However, variability affects the state and health of lake ecosystems at all time scales. Besides measuring time series at sufficiently high temporal resolution, the investigation of the full spectrum of variability requires innovative methods of analysis.
Analyzing observational data in the time frequency domain allows to identify variability at different time scales and facilitates their attribution to specific processes. The merit of this approach is subsequently demonstrated in three case studies. The first study uses a conceptual analysis to demonstrate the importance of time scales for the detection of ecosystem responses to climate change. These responses often occur during critical time windows in the year, may exhibit a time lag and can be driven by the exceedance of thresholds in their drivers. This can only be detected if the temporal resolution of the data is high enough. The second study applies Fast Fourier Transform spectral analysis to two decades of daily water temperature measurements to show how temporal and spatial scales of water temperature variability can serve as an indicator for mixing in a shallow, polymictic lake. The final study uses wavelet coherence as a diagnostic tool for limnology on a multivariate high-frequency data set recorded between the onset of ice cover and a cyanobacteria summer bloom in the year 2009 in a polymictic lake. Synchronicities among limnological and meteorological time series in narrow frequency bands were used to identify and disentangle prevailing limnological processes.
Beyond the novel empirical findings reported in the three case studies, this thesis aims to more generally be of interest to researchers dealing with now increasingly available time series data at high temporal resolution. A set of innovative methods to attribute patterns to processes, their drivers and constraints is provided to help make more efficient use of this kind of data.
The ionosphere, which is strongly influenced by the Sun, is known to be also affected by meteorological processes. These processes, despite having their origin in the troposphere and stratosphere, interact with the upper atmosphere. Such an interaction between atmospheric layers is known as vertical coupling. During geomagnetically quiet times, when near-Earth space is not under the influence of solar storms, these processes become important drivers for ionospheric variability. Studying the link between these processes in the lower atmosphere and the ionospheric variability is important for our understanding of fundamental mechanisms in ionospheric and meteorological research.
A prominent example of vertical coupling between the stratosphere and the ionosphere are the so-called stratospheric sudden warming (SSW) events that occur usually during northern winters and result in an increase in the polar stratospheric temperature and a reversal of the circumpolar winds. While the phenomenon of SSW is confined to the northern polar stratosphere, its influence on the ionosphere can be observed even at equatorial latitudes. During SSW events, the connection between the polar stratosphere and the equatorial ionosphere is believed to be through the modulation of global atmospheric tides. These tides are fundamental for the ionospheric E-region wind dynamo that generates electric fields and currents in the ionosphere. Observations of ionospheric currents indicate a large enhancement of the semidiurnal lunar tide in response to SSW events. Thus, the semidiurnal lunar tide becomes an important driver of ionospheric variability during SSW events.
In this thesis, the ionospheric effect of SSW events is investigated in the equatorial region, where a narrow but an intense E-region current known as the equatorial electrojet (EEJ) flows above the dip equator during the daytime. The day-to-day variability of the EEJ can be determined from magnetic field records at geomagnetic observatories close to the dip equator. Such magnetic data are available for several decades and allows to investigate the impact of SSW events on the EEJ and, even more importantly, helps in understanding the effects of SSW events on the equatorial ionosphere. An excellent long-term record of the geomagnetic field at the equator from 1922 onwards is available for the observatory Huancayo in Peru and is extensively utilized in this study.
The central subject of this thesis is the investigation of lunar tides in the EEJ during SSW events by analyzing long time series. This is done by estimating the lunar tidal amplitude in the EEJ from the magnetic records at Huancayo and by comparing them to measurements of the polar stratospheric wind and temperature, which led to the identification of the known SSW events from 1952 onwards. One goal of this thesis is to identify SSW events that predate 1952. To this end, superposed epoch analysis (SEA) is employed to establish a relationship between the lunar tidal power and the wind and temperature conditions in the lower atmosphere. A threshold value for the lunar tidal power is identified that is discriminative for the known SSW events. This threshold is then used to identify lunar tidal enhancements, which are indicative for any historic SSW events prior to 1952. It can be shown, that the number of lunar tidal enhancements and thus the occurrence frequency of historic SSW events between 1926 and 1952 is similar to the occurrence frequency of the known SSW events from 1952 onwards.
Next to the classic SSW definition, the concept of polar vortex weakening (PVW) is utilized in this thesis. PVW is defined for higher latitudes and altitudes (≈ 40km) than the classical SSW definition (≈ 32km). The correlation between the timing and magnitude of lunar tidal enhancements in the EEJ and the timing and magnitude of PVW is found to be better than for the classic SSW definition. This suggests that the lunar tidal enhancements in the EEJ are closely linked to the state of the middle atmosphere.
Geomagnetic observatories located in different longitudes at the dip equator allow investigating the longitudinally dependent variability of the EEJ during SSW events. For this purpose, the lunar tidal enhancements in the EEJ are determined for the Peruvian and Indian sectors during the major SSW events of the years 2006 and 2009. It is found that the lunar tidal amplitude shows similar enhancements in the Peruvian sector during both SSW events, while the enhancements are notably different for the two events in the Indian sector.
In summary, this thesis shows that lunar tidal enhancements in the EEJ are indeed correlated to the occurrence of SSW events and they should be considered a prominent driver of low latitude ionospheric variability. Secondly, lunar tidal enhancements are found to be longitudinally variable. This suggests that regional effects, such as ionospheric conductivity and the geometry and strength of the geomagnetic field, also play an important role and have to be considered when investigating the mechanisms behind vertical coupling.
In the arable soil landscape of hummocky ground moraines, an erosion-affected spatial differentiation of soils can be observed. Man-made erosion leads to soil profile modifications along slopes with changed solum thickness and modified properties of soil horizons due to water erosion in combination with tillage operations. Soil erosion creates, thereby, spatial patterns of soil properties (e.g., texture and organic matter content) and differences in crop development. However, little is known about the manner in which water fluxes are affected by soil-crop interactions depending on contrasting properties of differently-developed soil horizons and how water fluxes influence the carbon transport in an eroded landscape. To identify such feedbacks between erosion-induced soil profile modifications and the 1D-water and solute balance, high-precision weighing lysimeters equipped with a wide range of sensor technique were filled with undisturbed soil monoliths that differed in the degree of past soil erosion. Furthermore, lysimeter effluent concentrations were analyzed for dissolved carbon fractions in bi-weekly intervals.
The water balance components measured by high precision lysimeters varied from the most eroded to the less eroded monolith up to 83 % (deep drainage) primarily caused due to varying amounts of precipitation and evapotranspiration for a 3-years period. Here, interactions between crop development and contrasting rainfall interception by above ground biomass could explain differences in water balance components. Concentrations of dissolved carbon in soil water samples were relatively constant in time, suggesting carbon leaching was mainly affected by water fluxes in this observation period. For the lysimeter-based water balance analysis, a filtering scheme was developed considering temporal autocorrelation. The minute-based autocorrelation analysis of mass changes from lysimeter time series revealed characteristic autocorrelation lengths ranging from 23 to 76 minutes. Thereby, temporal autocorrelation provided an optimal approximation of precipitation quantities. However, the high temporal resolution in lysimeter time series is restricted by the lengths of autocorrelation.
Erosion-induced but also gradual changes in soil properties were reflected by dynamics of soil water retention properties in the lysimeter soils. Short-term and long-term hysteretic water retention data suggested seasonal wettability problems of soils increasingly limited rewetting of previously dried pore regions. Differences in water retention were assigned to soil tillage operations and the erosion history at different slope positions. The threedimensional spatial pattern of soil types that result from erosional soil profile modifications were also reflected in differences of crop root development at different landscape positions. Contrasting root densities revealed positive relations of root and aboveground plant characteristics. Differences in the spatially-distributed root growth between different eroded soil types provided indications that root development was affected by the erosion-induced soil evolution processes.
Overall, the current thesis corroborated the hypothesis that erosion-induced soil profile modifications affect the soil water balance, carbon leaching and soil hydraulic properties, but also the crop root system is influenced by erosion-induced spatial patterns of soil properties in the arable hummocky post glacial soil landscape. The results will help to improve model predictions of water and solute movement in arable soils and to understand interactions between soil erosion and carbon pathways regarding sink-or-source terms in landscapes.
Modeling loss-propagation in the global supply network: The dynamic agent-based model acclimate
(2017)
World markets are highly interlinked and local economies extensively rely on global supply and value chains. Consequently, local production disruptions, for instance caused by extreme weather events, are likely to induce indirect losses along supply chains with potentially global repercussions. These complex loss dynamics represent a challenge for comprehensive disaster risk assessments. Here, we introduce the numerical agent-based model acclimate designed to analyze the cascading of economic losses in the global supply network. Using national sectors as agents, we apply the model to study the global propagation of losses induced by stylized disasters. We find that indirect losses can become comparable in size to direct ones, but can be efficiently mitigated by warehousing and idle capacities. Consequently, a comprehensive risk assessment cannot focus solely on first-tier suppliers, but has to take the whole supply chain into account. To render the supply network climate-proof, national adaptation policies have to be complemented by international adaptation efforts. In that regard, our model can be employed to assess reasonable leverage points and to identify dynamic bottlenecks inaccessible to static analyses. (C) 2017 Elsevier B.V. All rights reserved.
Model-Based attribution of high-resolution streamflow trends in two alpine basins of Western Austria
(2017)
Several trend studies have shown that hydrological conditions are changing considerably in the Alpine region. However, the reasons for these changes are only partially understood and trend analyses alone are not able to shed much light. Hydrological modelling is one possible way to identify the trend drivers, i.e., to attribute the detected streamflow trends, given that the model captures all important processes causing the trends. We modelled the hydrological conditions for two alpine catchments in western Austria (a large, mostly lower-altitude catchment with wide valley plains and a nested high-altitude, glaciated headwater catchment) with the distributed, physically-oriented WaSiM-ETH model, which includes a dynamical glacier module. The model was calibrated in a transient mode, i.e., not only on several standard goodness measures and glacier extents, but also in such a way that the simulated streamflow trends fit with the observed ones during the investigation period 1980 to 2007. With this approach, it was possible to separate streamflow components, identify the trends of flow components, and study their relation to trends in atmospheric variables. In addition to trends in annual averages, highly resolved trends for each Julian day were derived, since they proved powerful in an earlier, data-based attribution study. We were able to show that annual and highly resolved trends can be modelled sufficiently well. The results provide a holistic, year-round picture of the drivers of alpine streamflow changes: Higher-altitude catchments are strongly affected by earlier firn melt and snowmelt in spring and increased ice melt throughout the ablation season. Changes in lower-altitude areas are mostly caused by earlier and lower snowmelt volumes. All highly resolved trends in streamflow and its components show an explicit similarity to the local temperature trends. Finally, results indicate that evapotranspiration has been increasing in the lower altitudes during the study period.
The eastern flanks of the Central Andes are characterized by deep convection, exposing them to hydrometeorological extreme events, often resulting in floods and a variety of mass movements. We assessed the spatiotemporal pattern of rainfall trends and the changes in the magnitude and frequency of extreme events (ae<yen>95th percentile) along an E-W traverse across the southern Central Andes using rain-gauge and high-resolution gridded datasets (CPC-uni and TRMM 3B42 V7). We generated different climate indices and made three key observations: (1) an increase of the annual rainfall has occurred at the transition between low (< 0.5 km) and intermediate (0.5-3 km) elevations between 1950 and 2014. Also, rainfall increases during the wet season and, to a lesser degree, decreases during the dry season. Increasing trends in annual total amounts characterize the period 1979-2014 in the arid, high-elevation southern Andean Plateau, whereas trend reversals with decreasing annual total amounts were found at low elevations. (2) For all analyzed periods, we observed small or no changes in the median values of the rainfall-frequency distribution, but significant trends with intensification or attenuation in the 95th percentile. (3) In the southern Andean Plateau, extreme rainfall events exhibit trends towards increasing magnitude and, to a lesser degree, frequency during the wet season, at least since 1979. Our analysis revealed that low (< 0.5 km), intermediate (0.5-3 km), and high-elevation (> 3 km) areas respond differently to changing climate conditions, and the transition zone between low and intermediate elevations is characterized by the most significant changes.
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.
Central Anatolia is a low-relief, high-elevation region where decadal-scale deformation rates estimated from space geodesy suggest low strain rates within a stiff microplate. However, numerous Quaternary faults have been mapped within this low-strain region and estimating their slip rate and seismic potential is important for hazard assessments in an area of increasing infrastructural development. Here we focus on the Sultanhani Fault (SF), which constitutes an integral part of the Eskisehir-Cihanbeyli Fault System, and use deformed maximum highstand shorelines of palaeo-lake Konya to estimate tectonic slip rates at millennial scale. Some of these shorelines were previously interpreted as fault scarps, but we provide conclusive evidence for their erosional origin. We found that shoreline-angle elevations estimated from differential GPS profiles record vertical displacements of 10.2 m across the SF. New radiocarbon ages of lacustrine molluscs suggest 22.4 m of relative lake-level fall between 22.1 +/- 0.3 and 21.7 +/- 0.4 cal. kaBP, constraining the timing of abrupt abandonment of the highstand shoreline. Models of lithospheric rebound associated with regressions of the Tuz Golu and Konya palaeolakes predict only similar to 1 m of regional-scale uplift across the Konya Basin. Dislocation models of displaced shorelines suggest fault-slip rates of 1.5 and 1.8 mm yr(-1) for planar and listric fault geometries, respectively, providing reasonable results for the latter. We found fault scarps in the Nasuhpinar mudflat that likely represent the most recent ground-breaking rupture of the SF, with an average vertical displacement of 1.2 +/- 0.5 m estimated from 54 topographic profiles, equivalent to a M similar to 6.5-6.9 earthquake based on empirical scaling laws. If such events were characteristic during the ultimate 21 ka, a relatively short recurrence time of similar to 800-900 yr would be needed to account for the millennial slip rate. Alternatively, the fault scarp at Nasuhpinar might represent a larger earthquake requiring more frequent smaller events to account for the millennial rate. The relatively fast slip rate of the SF over the past 21 ka is unlikely to have persisted over longer timescales and might reflect spatiotemporal variations in deformation rates within kinematically-linked fault systems within Central Anatolia, or a transient perturbation to the local stress field or fault strength. Such perturbation might have been related to climatically controlled changes in surface and near-surface loads and by interactions among the different tectonic processes that have been proposed to drive the overall slow uplift and associated extension in the Central Anatolian Plateau.
Understanding how Earth-surface processes respond to past climatic perturbations is crucial for making informed predictions about future impacts of climate change on sediment "uxes. Sedimentary records provide the archives for inferring these processes, but their interpretation is compromised by our incomplete understanding of how sediment-routing systems respond to millennial-scale climate cycles. We analyzed seven sediment cores recovered from marine turbidite depositional sites along the Chile continental margin. The sites span a pronounced arid-to-humid gradient with variable relief and related sediment connectivity of terrestrial and marine environments. These sites allowed us to study event related depositional processes in different climatic and geomorphic settings from the Last Glacial Maximum to the present day. The three sites reveal a steep decline of turbidite deposition during deglaciation. High rates of sea-level rise postdate the decline in turbidite deposition. Comparison with paleoclimate proxies documents that the spatio-temporal sedimentary pattern rather mirrors the deglacial humidity decrease and concomitant warming with no resolvable lag times. Our results let us infer that declining deglacial humidity decreased "uvial sediment supply. This signal propagated rapidly through the highly connected systems into the marine sink in north-central Chile. In contrast, in south-central Chile, connectivity between the Andean erosional zone and the "uvial transfer zone probably decreased abruptly by sediment trapping in piedmont lakes related to deglaciation, resulting in a sudden decrease of sediment supply to the ocean. Additionally, reduced moisture supply may have contributed to the rapid decline of turbidite deposition. These different causes result in similar depositional patterns in the marine sinks. We conclude that turbiditic strata may constitute reliable recorders of climate change across a wide range of climatic zones and geomorphic conditions. However, the underlying causes for similar signal manifestations in the sinks may differ, ranging from maintained high system connectivity to abrupt connectivity loss. (C) 2017 Elsevier B.V. All rights reserved.
The northern edge of the western central Tien Shan range is bounded by the Issyk-Ata fault situated south of Bishkek, the capital of Kyrgyzstan. Contraction in this thick-skinned orogen occurs with low-strain accumulation and long earthquake recurrence intervals. In the nineteenth to twentieth centuries, a sequence of large earthquakes with magnitudes between 6.9 and 8 affected the northern Tien Shan but left nearly the entire extent of the Issyk-Ata fault unruptured. Here, the only known historic earthquake ruptured in A.D. 1885 (M6.9) along the western end of the Issyk-Ata fault. Because earthquakes in low-strain regions often tend to cluster in time and may promote failure along nearby structures, the earthquake history of the northern Tien Shan represents an exceptional structural setting for studying fault behavior affected by an intraplate earthquake sequence. We present a paleoseismological study from one site (Belek) along the Issyk-Ata fault located east of the A.D. 1885 epicentral area. Our analysis combines a range of tools, including photogrammetry, differential Global Positioning System, 3D visualization, and age modeling with different dating methods (infrared stimulated luminescence, radiocarbon, U-series) to improve the reliability of an event chronology for the trench stratigraphy and fault geometry. We were able to distinguish three different surfacerupturing paleoearthquakes; these affected the area before 10.5 +/- 1.1 cal ka B.P., at similar to 5.6 +/- 1.0 cal ka B.P., and at similar to 630 +/- 100 cal B.P., respectively. Associated paleomagnitudes for the last two earthquakes range between M6.7 and 7.4, with a cumulative slip rate of 0.7 +/- 0.32 mm/yr. We did not find evidence for the A.D. 1885 event at Belek. Our study yielded two main overall results: first, it extends the regional historic and paleoseismic record; second, the documented rupture events along the Issyk-Ata fault suggest that this fault was not affected in its entirety; instead, these events indicate segmented rupture behavior.
This study analyzes the discharge variability of small to medium drainage basins (10(2)-10(4) km(2)) in the southern Central Andes of NW Argentina. The Hilbert-Huang Transform (HHT) was applied to evaluate non-stationary oscillatory modes of variability and trends, based on four time series of monthly normalized discharge anomaly between 1940 and 2015. Statistically significant trends reveal increasing discharge during the past decades and document an intensification of the hydrological cycle during this period. An Ensemble Empirical Mode Decomposition (EEMD) analysis revealed that discharge variability in this region can be best described by five quasi-periodic statistically significant oscillatory modes, with mean periods varying from 1 to 20 y. Moreover, we show that discharge variability is most likely linked to the phases of the Pacific Decadal Oscillation (PDO) at multi-decadal timescales (similar to 20 y) and, to a lesser degree, to the Tropical South Atlantic SST anomaly (TSA) variability at shorter timescales (similar to 2-5 y). Previous studies highlighted a rapid increase in discharge in the southern Central Andes during the 1970s, inferred to have been associated with the global 1976-77 climate shift. Our results suggest that the rapid discharge increase in the NW Argentine Andes coincides with the periodic enhancement of discharge, which is mainly linked to a negative to positive transition of the PDO phase and TSA variability associated with a long-term increasing trend. We therefore suggest that variations in discharge in this region are largely driven by both natural variability and the effects of global climate change. We furthermore posit that the links between atmospheric and hydrologic processes result from a combination of forcings that operate on different spatiotemporal scales. (C) 2017 Elsevier B.V. All rights reserved.
Water infiltration in soil is not only affected by the inherent heterogeneities of soil, but even more by the interaction with plant roots and their water uptake. Neutron tomography is a unique non-invasive 3D tool to visualize plant root systems together with the soil water distribution in situ. So far, acquisition times in the range of hours have been the major limitation for imaging 3D water dynamics. Implementing an alternative acquisition procedure we boosted the speed of acquisition capturing an entire tomogram within 10 s. This allows, for the first time, tracking of a water front ascending in a rooted soil column upon infiltration of deuterated water time-resolved in 3D. Image quality and resolution could be sustained to a level allowing for capturing the root system in high detail. Good signal-to-noise ratio and contrast were the key to visualize dynamic changes in water content and to localize the root uptake. We demonstrated the ability of ultra-fast tomography to quantitatively image quick changes of water content in the rhizosphere and outlined the value of such imaging data for 3D water uptake modelling. The presented method paves the way for time-resolved studies of various 3D flow and transport phenomena in porous systems
The temporal dynamics of climate processes are spread across different timescales and, as such, the study of these processes at only one selected timescale might not reveal the complete mechanisms and interactions within and between the (sub-) processes. To capture the non-linear interactions between climatic events, the method of event synchronization has found increasing attention recently. The main drawback with the present estimation of event synchronization is its restriction to analysing the time series at one reference timescale only. The study of event synchronization at multiple scales would be of great interest to comprehend the dynamics of the investigated climate processes. In this paper, the wavelet-based multi-scale event synchronization (MSES) method is proposed by combining the wavelet transform and event synchronization. Wavelets are used extensively to comprehend multi-scale processes and the dynamics of processes across various timescales. The proposed method allows the study of spatio-temporal patterns across different timescales. The method is tested on synthetic and real-world time series in order to check its replicability and applicability. The results indicate that MSES is able to capture relationships that exist between processes at different timescales.
Fluvial fill terraces in intermontane basins are valuable geomorphic archives that can record tectonically and/or climatically driven changes of the Earth-surface process system. However, often the preservation of fill terrace sequences is incomplete and/or they may form far away from their source areas, complicating the identification of causal links between forcing mechanisms and landscape response, especially over multi-millennial timescales. The intermontane Toro Basin in the southern Central Andes exhibits at least five generations of fluvial terraces that have been sculpted into several-hundred-meter-thick Quaternary valley-fill conglomerates. New surface-exposure dating using nine cosmogenic Be-10 depth profiles reveals the successive abandonment of these terraces with a 100 kyr cyclicity between 75 +/- 7 and 487 +/- 34 ka. Depositional ages of the conglomerates, determined by four Al-26/Be-10 burial samples and U-Pb zircon ages of three intercalated volcanic ash beds, range from 18 +/- 141 to 936 +/- 170 ka, indicating that there were multiple cut-and-fill episodes. Although the initial onset of aggradation at similar to 1 Ma and the overall net incision since ca. 500 ka can be linked to tectonic processes at the narrow basin outlet, the superimposed 100 kyr cycles of aggradation and incision are best explained by eccentricity-driven climate change. Within these cycles, the onset of river incision can be correlated with global cold periods and enhanced humid phases recorded in paleoclimate archives on the adjacent Bolivian Altiplano, whereas deposition occurred mainly during more arid phases on the Altiplano and global interglacial periods. We suggest that enhanced runoff during global cold phases - due to increased regional precipitation rates, reduced evapotranspiration, or both - resulted in an increased sediment-transport capacity in the Toro Basin, which outweighed any possible increases in upstream sediment supply and thus triggered incision. Compared with two nearby basins that record precessional (21-kyr) and long-eccentricity (400-kyr) forcing within sedimentary and geomorphic archives, the recorded cyclicity scales with the square of the drainage basin length. (C) 2017 Elsevier B.V. All rights reserved.
Water infiltration in soil is not only affected by the inherent heterogeneities of soil, but even more by the interaction with plant roots and their water uptake. Neutron tomography is a unique non-invasive 3D tool to visualize plant root systems together with the soil water distribution in situ. So far, acquisition times in the range of hours have been the major limitation for imaging 3D water dynamics. Implementing an alternative acquisition procedure we boosted the speed of acquisition capturing an entire tomogram within 10 s. This allows, for the first time, tracking of a water front ascending in a rooted soil column upon infiltration of deuterated water time-resolved in 3D. Image quality and resolution could be sustained to a level allowing for capturing the root system in high detail. Good signal-to-noise ratio and contrast were the key to visualize dynamic changes in water content and to localize the root uptake. We demonstrated the ability of ultra-fast tomography to quantitatively image quick changes of water content in the rhizosphere and outlined the value of such imaging data for 3D water uptake modelling. The presented method paves the way for time-resolved studies of various 3D flow and transport phenomena in porous systems.
Scenario-neutral response surfaces illustrate the sensitivity of a simulated natural system, represented by a specific impact variable, to systematic perturbations of climatic parameters. This type of approach has recently been developed as an alternative to top-down approaches for the assessment of climate change impacts. A major limitation of this approach is the underrepresentation of changes in the temporal structure of the climate input data (i.e., the seasonal and day-to-day variability) since this is not altered by the perturbation. This paper presents a framework that aims to examine this limitation by perturbing both observed and projected climate data time series for a future period, which both serve as input into a hydrological model (the HBV model). The resulting multiple response surfaces are compared at a common domain, the standardized runoff response surface (SRRS). We apply this approach in a case study catchment in Norway to (i) analyze possible changes in mean and extreme runoff and (ii) quantify the influence of changes in the temporal structure represented by 17 different climate input sets using linear mixed-effect models. Results suggest that climate change induced increases in mean and peak flow runoff and only small changes in low flow. They further suggest that the effect of the different temporal structures of the climate input data considerably affects low flows and floods (at least 21% influence), while it is negligible for mean runoff.
The characteristics of a landscape pose essential factors for hydrological processes. Therefore, an adequate representation of the landscape of a catchment in hydrological models is vital. However, many of such models exist differing, amongst others, in spatial concept and discretisation. The latter constitutes an essential pre-processing step, for which many different algorithms along with numerous software implementations exist. In that context, existing solutions are often model specific, commercial, or depend on commercial back-end software, and allow only a limited or no workflow automation at all. Consequently, a new package for the scientific software and scripting environment R, called lumpR, was developed. lumpR employs an algorithm for hillslope-based landscape discretisation directed to large-scale application via a hierarchical multi-scale approach. The package addresses existing limitations as it is free and open source, easily extendible to other hydrological models, and the workflow can be fully automated. Moreover, it is user-friendly as the direct coupling to a GIS allows for immediate visual inspection and manual adjustment. Sufficient control is furthermore retained via parameter specification and the option to include expert knowledge. Conversely, completely automatic operation also allows for extensive analysis of aspects related to landscape discretisation. In a case study, the application of the package is presented. A sensitivity analysis of the most important discretisation parameters demonstrates its efficient workflow automation. Considering multiple streamflow metrics, the employed model proved reasonably robust to the discretisation parameters. However, parameters determining the sizes of subbasins and hillslopes proved to be more important than the others, including the number of representative hillslopes, the number of attributes employed for the lumping algorithm, and the number of sub-discretisations of the representative hillslopes.
The information about climate change impact on river discharge is vitally important for planning adaptation measures. The future changes can affect different water-related sectors. The main goal of this study was to investigate the potential water resource changes in Ukraine, focusing on three mesoscale river catchments (Teteriv, UpperWestern Bug, and Samara) characteristic for different geographical zones. The catchment scale watershed model—Soil and Water Integrated Model (SWIM)—was setup, calibrated, and validated for the three catchments under consideration. A set of seven GCM-RCM (General Circulation Model-Regional Climate Model) coupled climate scenarios corresponding to RCPs (Representative Concentration Pathways) 4.5 and 8.5 were used to drive the hydrological catchment model. The climate projections, used in the study, were considered as three combinations of low, intermediate, and high end scenarios. Our results indicate the shifts in the seasonal distribution of runoff in all three catchments. The spring high flow occurs earlier as a result of temperature increases and earlier snowmelt. The fairly robust trend is an increase in river discharge in the winter season, and most of the scenarios show a potential decrease in river discharge in the spring.
The information about climate change impact on river discharge is vitally important for planning adaptation measures. The future changes can affect different water-related sectors. The main goal of this study was to investigate the potential water resource changes in Ukraine, focusing on three mesoscale river catchments (Teteriv, UpperWestern Bug, and Samara) characteristic for different geographical zones. The catchment scale watershed model—Soil and Water Integrated Model (SWIM)—was setup, calibrated, and validated for the three catchments under consideration. A set of seven GCM-RCM (General Circulation Model-Regional Climate Model) coupled climate scenarios corresponding to RCPs (Representative Concentration Pathways) 4.5 and 8.5 were used to drive the hydrological catchment model. The climate projections, used in the study, were considered as three combinations of low, intermediate, and high end scenarios. Our results indicate the shifts in the seasonal distribution of runoff in all three catchments. The spring high flow occurs earlier as a result of temperature increases and earlier snowmelt. The fairly robust trend is an increase in river discharge in the winter season, and most of the scenarios show a potential decrease in river discharge in the spring.
The characteristics of a landscape pose essential factors for hydrological processes. Therefore, an adequate representation of the landscape of a catchment in hydrological models is vital. However, many of such models exist differing, amongst others, in spatial concept and discretisation. The latter constitutes an essential pre-processing step, for which many different algorithms along with numerous software implementations exist. In that context, existing solutions are often model specific, commercial, or depend on commercial back-end software, and allow only a limited or no workflow automation at all.
Consequently, a new package for the scientific software and scripting environment R, called lumpR, was developed. lumpR employs an algorithm for hillslope-based landscape discretisation directed to large-scale application via a hierarchical multi-scale approach. The package addresses existing limitations as it is free and open source, easily extendible to other hydrological models, and the workflow can be fully automated. Moreover, it is user-friendly as the direct coupling to a GIS allows for immediate visual inspection and manual adjustment. Sufficient control is furthermore retained via parameter specification and the option to include expert knowledge. Conversely, completely automatic operation also allows for extensive analysis of aspects related to landscape discretisation.
In a case study, the application of the package is presented. A sensitivity analysis of the most important discretisation parameters demonstrates its efficient workflow automation. Considering multiple streamflow metrics, the employed model proved reasonably robust to the discretisation parameters. However, parameters determining the sizes of subbasins and hillslopes proved to be more important than the others, including the number of representative hillslopes, the number of attributes employed for the lumping algorithm, and the number of sub-discretisations of the representative hillslopes.
Global earthquake locations are often associated with very large systematic travel-time residuals even for clear arrivals, especially for regional and near-regional stations in subduction zones because of their strongly heterogeneous velocity structure. Travel-time corrections can drastically reduce travel-time residuals at regional stations and, in consequence, improve the relative location accuracy. We have extended the shrinking-box source-specific station terms technique to regional and teleseismic distances and adopted the algorithm for probabilistic, nonlinear, global-search location. We evaluated the potential of the method to compute precise relative hypocentre locations on a global scale. The method has been applied to two specific test regions using existing P- and pP-phase picks. The first data set consists of 3103 events along the Chilean margin and the second one comprises 1680 earthquakes in the Tonga-Fiji subduction zone. Pick data were obtained from the GEOFON earthquake bulletin, produced using data from all available, global station networks. A set of timing corrections varying as a function of source position was calculated for each seismic station. In this way, we could correct the systematic errors introduced into the locations by the inaccuracies in the assumed velocity structure without explicitly solving for a velocity model. Residual statistics show that the median absolute deviation of the travel-time residuals is reduced by 40-60 per cent at regional distances, where the velocity anomalies are strong. Moreover, the spread of the travel-time residuals decreased by similar to 20 per cent at teleseismic distances (>28 degrees). Furthermore, strong variations in initial residuals as a function of recording distance are smoothed out in the final residuals. The relocated catalogues exhibit less scattered locations in depth and sharper images of the seismicity associated with the subducting slabs. Comparison with a high-resolution local catalogue reveals that our relocation process significantly improves the hypocentre locations compared to standard locations.
The mechanisms leading to large earthquakes are poorly understood and documented. Here we characterize the long-term precursory phase of the 1 April 2014 M(w)8.1 North Chile megathrust. We show that a group of coastal GPS stations accelerated westward 8months before the main shock, corresponding to a M(w)6.5 slow slip event on the subduction interface, 80% of which was aseismic. Concurrent interface foreshocks underwent a diminution of their radiation at high frequency, as shown by the temporal evolution of Fourier spectra and residuals with respect to ground motions predicted by recent subduction models. Such ground motions change suggests that in response to the slow sliding of the subduction interface, seismic ruptures are progressively becoming smoother and/or slower. The gradual propagation of seismic ruptures beyond seismic asperities into surrounding metastable areas could explain these observations and might be the precursory mechanism eventually leading to the main shock.
Steep mountain channels are an important component of the fluvial system. On geological timescales, they shape mountain belts and counteract tectonic uplift by erosion. Their channels are strongly coupled to hillslopes and they are often the main source of sediment transported downstream to low-gradient rivers and to alluvial fans, where commonly settlements in mountainous areas are located. Hence, mountain streams are the cause for one of the main natural hazards in these regions. Due to climate change and a pronounced populating of mountainous regions the attention given to this threat is even growing. Although quantitative studies on sediment transport have significantly advanced our knowledge on measuring and calibration techniques we still lack studies of the processes within mountain catchments. Studies examining the mechanisms of energy and mass exchange on small temporal and spatial scales in steep streams remain sparse in comparison to low-gradient alluvial channels.
In the beginning of this doctoral project, a vast amount of experience and knowledge of a steep stream in the Swiss Prealps had to be consolidated in order to shape the principal aim of this research effort. It became obvious, that observations from within the catchment are underrepresented in comparison to experiments performed at the catchment’s outlet measuring fluxes and the effects of the transported material. To counteract this imbalance, an examination of mass fluxes within the catchment on the process scale was intended. Hence, this thesis is heavily based on direct field observations, which are generally rare in these environments in quantity and quality. The first objective was to investigate the coupling of the channel with surrounding hillslopes, the major sources of sediment. This research, which involved the monitoring of the channel and adjacent hillslopes, revealed that alluvial channel steps play a key role in coupling of channel and hillslopes. The observations showed that hillslope stability is strongly associated with the step presence and an understanding of step morphology and stability is therefore crucial in understanding sediment mobilization. This finding refined the way we think about the sediment dynamics in steep channels and motivated continued research of the step dynamics. However, soon it became obvious that the technological basis for developing field tests and analyzing the high resolution geometry measured in the field was not available. Moreover, for many geometrical quantities in mountain channels definitions and a clear scientific standard was not available. For example, these streams are characterized by a high spatial variability of the channel banks, preventing straightforward calculations of the channel width without a defined reference. Thus, the second and inevitable part of this thesis became the development and evaluation of scientific tools in order to investigate the geometrical content of the study reach thoroughly. The developed framework allowed the derivation of various metrics of step and channel geometry which facilitated research on the a large data set of observations of channel steps. In the third part, innovative, physically-based metrics have been developed and compared to current knowledge on step formation, suggested in the literature. With this analyses it could be demonstrated that the formation of channel steps follow a wide range of hydraulic controls. Due to the wide range of tested parameters channel steps observed in a natural stream were attributed to different mechanisms of step formation, including those based on jamming and those based on key-stones. This study extended our knowledge on step formation in a steep stream and harmonized different, often time seen as competing, processes of step formation. This study was based on observations collected at one point in time. In the fourth part of this project, the findings of the snap-shot observations were extended in the temporal dimension and the derived concepts have been utilized to investigate reach-scale step patterns in response to large, exceptional flood events. The preliminary results of this work based on the long-term analyses of 7 years of long profile surveys showed that the previously observed channel-hillslope mechanism is the responsible for the short-term response of step formation.
The findings of the long-term analyses of step patterns drew a bow to the initial observations of a channel-hillslope system which allowed to join the dots in the dynamics of steep stream. Thus, in this thesis a broad approach has been chosen to gain insights into the complex system of steep mountain rivers. The effort includes in situ field observations (article I), the development of quantitative scientific tools (article II), the reach-scale analyses of step-pool morphology (article III) and its temporal evolution (article IV). With this work our view on the processes within the catchment has been advanced towards a better mechanistic understanding of these fluvial system relevant to improve applied scientific work.
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.
Northeastern (NE) India experiences extraordinarily pronounced seasonal climate, governed by the Indian summer monsoon (ISM). The vulnerability of this region to floods and droughts calls for detailed and highly resolved paleoclimate reconstructions to assess the recurrence rate and driving factors of ISM changes. We use stable oxygen and carbon isotope ratios (delta O-18 and delta C-13) from stalagmite MAW-6 from Mawmluh Cave to infer climate and environmental conditions in NE India over the last deglaciation (16-6ka). We interpret stalagmite delta O-18 as reflecting ISM strength, whereas delta C-13 appears to be driven by local hydroclimate conditions. Pronounced shifts in ISM strength over the deglaciation are apparent from the delta O-18 record, similarly to other records from monsoonal Asia. The ISM is weaker during the late glacial (LG) period and the Younger Dryas, and stronger during the BOlling-Allerod and Holocene. Local conditions inferred from the delta C-13 record appear to have changed less substantially over time, possibly related to the masking effect of changing precipitation seasonality. Time series analysis of the delta O-18 record reveals more chaotic conditions during the late glacial and higher predictability during the Holocene, likely related to the strengthening of the seasonal recurrence of the ISM with the onset of the Holocene.
Rhizosphere processes are highly dynamic in time and space and strongly depend on each other. Key factors influencing pH changes in the rhizosphere are root exudation, respiration, and nutrient supply, which are influenced by soil water content levels. In this study, we measured the real-time distribution of soil water, pH changes, and oxygen distribution in the rhizosphere of young maize plants using a recently developed imaging approach. Neutron radiography was used to capture the root system and soil water distribution, while fluorescence imaging was employed to map soil pH and soil oxygen changes. Germinated seeds of maize (Zea mays L.) were planted in glass rhizotrons equipped with pH and oxygen-sensitive sensor foils. After 20 d, the rhizotrons were wetted from the bottom and time-lapsed images via fluorescence and neutron imaging were taken during the subsequent day and night cycles for 5 d. We found higher water content and stronger acidification in the first 0.5 mm from the root surface compared to the bulk soil, which could be a consequence of root exudation. While lateral roots only slightly acidified their rhizosphere, crown roots induced stronger acidification of up to 1 pH unit. We observed changing oxygen patterns at different soil moisture conditions and increasing towards lateral as well as crown roots while extending laterally with ongoing water logging. Our work indicates that plants alter the rhizosphere pH and oxygen also depending on root type, which may indirectly arise also from differences in age and water content changes. The results presented here were possible only by combining different imaging techniques to examine profiles at the root-soil interface in a comprehensive way during wetting and drying.
Flow patterns in conjunction with seasonal and diurnal temperature variations control ecological and biogeochemical conditions in hyporheic sediments. In particular, hyporheic temperatures have a great impact on many temperature-sensitive microbial processes. In this study, we used 3-D coupled water flow and heat transport simulations applying the HydroGeoSphere code in combination with high-resolution observations of hydraulic heads and temperatures to quantify reach-scale water and heat flux across the river-groundwater interface and hyporheic temperature dynamics of a lowland gravel bed river. The model was calibrated in order to constrain estimates of the most sensitive model parameters. The magnitude and variations of the simulated temperatures matched the observed ones, with an average mean absolute error of 0.7 degrees C and an average Nash Sutcliffe efficiency of 0.87. Our results indicate that nonsubmerged streambed structures such as gravel bars cause substantial thermal heterogeneity within the saturated sediment at the reach scale. Individual hyporheic flow path temperatures strongly depend on the flow path residence time, flow path depth, river, and groundwater temperature. Variations in individual hyporheic flow path temperatures were up to 7.9 degrees C, significantly higher than the daily average (2.8 degrees C), but still lower than the average seasonal hyporheic temperature difference (19.2 degrees C). The distribution between flow path temperatures and residence times follows a power law relationship with exponent of about 0.37. Based on this empirical relation, we further estimated the influence of hyporheic flow path residence time and temperature on oxygen consumption which was found to partly increase by up to 29% in simulations.
Water infiltration in soil is not only affected by the inherent heterogeneities of soil, but even more by the interaction with plant roots and their water uptake. Neutron tomography is a unique non-invasive 3D tool to visualize plant root systems together with the soil water distribution in situ. So far, acquisition times in the range of hours have been the major limitation for imaging 3D water dynamics. Implementing an alternative acquisition procedure we boosted the speed of acquisition capturing an entire tomogram within 10 s. This allows, for the first time, tracking of a water front ascending in a rooted soil column upon infiltration of deuterated water time-resolved in 3D. Image quality and resolution could be sustained to a level allowing for capturing the root system in high detail. Good signal-to-noise ratio and contrast were the key to visualize dynamic changes in water content and to localize the root uptake. We demonstrated the ability of ultra-fast tomography to quantitatively image quick changes of water content in the rhizosphere and outlined the value of such imaging data for 3D water uptake modelling. The presented method paves the way for time-resolved studies of various 3D flow and transport phenomena in porous systems.
Anthropogenically amplified erosion leads to increased fine-grained sediment input into the fluvial system in the 15.000 km2 Kharaa River catchment in northern Mongolia and constitutes a major stressing factor for the aquatic ecosystem. This study uniquely combines the application of intensive monitoring, source fingerprinting and catchment modelling techniques to allow for the comparison of the credibility and accuracy of each single method. High-resolution discharge data were used in combination with daily suspended solid measurements to calculate the suspended sediment budget and compare it with estimations of the sediment budget model SedNet. The comparison of both techniques showed that the development of an overall sediment budget with SedNet was possible, yielding results in the same order of magnitude (20.3 kt a- 1 and 16.2 kt a- 1).
Radionuclide sediment tracing, using Be-7, Cs-137 and Pb-210 was applied to differentiate sediment sources for particles < 10μm from hillslope and riverbank erosion and showed that riverbank erosion generates 74.5% of the suspended sediment load, whereas surface erosion contributes 21.7% and gully erosion only 3.8%. The contribution of the single subcatchments of the Kharaa to the suspended sediment load was assessed based on their variation in geochemical composition (e.g. in Ti, Sn, Mo, Mn, As, Sr, B, U, Ca and Sb). These variations were used for sediment source discrimination with geochemical composite fingerprints based on Genetic Algorithm driven Discriminant Function Analysis, the Kruskal–Wallis H-test and Principal Component Analysis. The contributions of the individual sub-catchment varied from 6.4% to 36.2%, generally showing higher contributions from the sub-catchments in the middle, rather than the upstream portions of the study area.
The results indicate that river bank erosion generated by existing grazing practices of livestock is the main cause for elevated fine sediment input. Actions towards the protection of the headwaters and the stabilization of the river banks within the middle reaches were identified as the highest priority. Deforestation and by lodging and forest fires should be prevented to avoid increased hillslope erosion in the mountainous areas. Mining activities are of minor importance for the overall catchment sediment load but can constitute locally important point sources for particular heavy metals in the fluvial system.
The characteristics of an aboveground cosmic-ray neutron sensor (CRNS) are evaluated for monitoring a mountain snowpack in the Austrian Alps from March 2014 to June 2016. Neutron counts were compared to continuous point-scale snow depth (SD) and snow-water-equivalent (SWE) measurements from an automatic weather station with a maximum SWE of 600 mm (April 2014). Several spatially distributed Terrestrial Laser Scanning (TLS)-based SD and SWE maps were additionally used. A strong nonlinear correlation is found for both SD and SWE. The representative footprint of the CRNS is in the range of 230-270 m. In contrast to previous studies suggesting signal saturation at around 100 mm of SWE, no complete signal saturation was observed. These results imply that CRNS could be transferred into an unprecedented method for continuous detection of spatially averaged SD and SWE for alpine snowpacks, though with sensitivity decreasing with increasing SWE. While initially different functions were found for accumulation and melting season conditions, this could be resolved by accounting for a limited measurement depth. This depth limit is in the range of 200 mm of SWE for dense snowpacks with high liquid water contents and associated snow density values around 450 kg m(-3) and above. In contrast to prior studies with shallow snowpacks, interannual transferability of the results is very high regardless of presnowfall soil moisture conditions. This underlines the unexpectedly high potential of CRNS to close the gap between point-scale measurements, hydrological models, and remote sensing of the cryosphere in alpine terrain.
Both Alpine and Mediterranean areas are considered sensitive to so-called global change, considered as the combination of climate and land use changes. All panels on climate evolution predict future scenarios of increasing frequency and magnitude of floods which are likely to lead to huge geomorphic adjustments of river channels so major metamorphosis of fluvial systems is expected as a result of global change. Such pressures are likely to give rise to major ecological and economic changes and challenges that governments need to address as a matter of priority. Changes in river flow regimes associated with global change are therefore ushering in a new era, where there is a critical need to evaluate hydro-geomorphological hazards from headwaters to lowland areas (flooding can be not just a problem related to being under the water). A key question is how our understanding of these hazards associated with global change can be improved; improvement has to come from integrated research which includes the climatological and physical conditions that could influence the hydrology and sediment generation and hence the conveyance of water and sediments (including the river’s capacity, i.e. amount of sediment, and competence, i.e. channel deformation) and the vulnerabilities and economic repercussions of changing hydrological hazards (including the evaluation of the hydro-geomorphological risks too).
Within this framework, the purpose of this international symposium is to bring together researchers from several disciplines as hydrology, fluvial geomorphology, hydraulic engineering, environmental science, geography, economy (and any other related discipline) to discuss the effects of global change over the river system in relation with floods. The symposium is organized by means of invited talks given by prominent experts, oral lectures, poster sessions and discussion sessions for each individual topic; it will try to improve our understanding of how rivers are likely to evolve as a result of global change and hence address the associated hazards of that fluvial environmental change concerning flooding.
Four main topics are going to be addressed:
- Modelling global change (i.e. climate and land-use) at relevant spatial (regional, local) and temporal (from the long-term to the single-event) scales.
- Measuring and modelling river floods from the hydrological, sediment transport (both suspended and bedload) and channel morphology points of view at different spatial (from the catchment to the reach) and temporal (from the long-term to the single-event) scales.
- Evaluation and assessment of current and future river flooding hazards and risks in a global change perspective.
- Catchment management to face river floods in a changing world.
We are very pleased to welcome you to Potsdam. We hope you will enjoy your participation at the International Symposium on the Effects of Global Change on Floods, Fluvial Geomorphology and Related Hazards in Mountainous Rivers and have an exciting and profitable experience. Finally, we would like to thank all speakers, participants, supporters, and sponsors for their contributions that for sure will make of this event a very remarkable and fruitful meeting. We acknowledge the valuable support of the European Commission (Marie Curie Intra-European Fellowship, Project ‘‘Floodhazards’’, PIEF-GA-2013-622468, Seventh EU Framework Programme) and the Deutschen Forschungsgemeinschaft (Research Training Group “Natural Hazards and Risks in a Changing World” (NatRiskChange; GRK 2043/1) as the symposium would not have been possible without their help. Without your cooperation, this symposium would not be either possible or successful.
Climate or land use?
(2017)
This study intends to contribute to the ongoing discussion on whether land use and land cover changes (LULC) or climate trends have the major influence on the observed increase of flood magnitudes in the Sahel. A simulation-based approach is used for attributing the observed trends to the postulated drivers. For this purpose, the ecohydrological model SWIM (Soil and Water Integrated Model) with a new, dynamic LULC module was set up for the Sahelian part of the Niger River until Niamey, including the main tributaries Sirba and Goroul. The model was driven with observed, reanalyzed climate and LULC data for the years 1950–2009. In order to quantify the shares of influence, one simulation was carried out with constant land cover as of 1950, and one including LULC. As quantitative measure, the gradients of the simulated trends were compared to the observed trend. The modeling studies showed that for the Sirba River only the simulation which included LULC was able to reproduce the observed trend. The simulation without LULC showed a positive trend for flood magnitudes, but underestimated the trend significantly. For the Goroul River and the local flood of the Niger River at Niamey, the simulations were only partly able to reproduce the observed trend. In conclusion, the new LULC module enabled some first quantitative insights into the relative influence of LULC and climatic changes. For the Sirba catchment, the results imply that LULC and climatic changes contribute in roughly equal shares to the observed increase in flooding. For the other parts of the subcatchment, the results are less clear but show, that climatic changes and LULC are drivers for the flood increase; however their shares cannot be quantified. Based on these modeling results, we argue for a two-pillar adaptation strategy to reduce current and future flood risk: Flood mitigation for reducing LULC-induced flood increase, and flood adaptation for a general reduction of flood vulnerability.
The characteristics of a landscape pose essential factors for hydrological processes. Therefore, an adequate representation of the landscape of a catchment in hydrological models is vital. However, many of such models exist differing, amongst others, in spatial concept and discretisation. The latter constitutes an essential pre-processing step, for which many different algorithms along with numerous software implementations exist. In that context, existing solutions are often model specific, commercial, or depend on commercial back-end software, and allow only a limited or no workflow automation at all.
Consequently, a new package for the scientific software and scripting environment R, called lumpR, was developed. lumpR employs an algorithm for hillslope-based landscape discretisation directed to large-scale application via a hierarchical multi-scale approach. The package addresses existing limitations as it is free and open source, easily extendible to other hydrological models, and the workflow can be fully automated. Moreover, it is user-friendly as the direct coupling to a GIS allows for immediate visual inspection and manual adjustment. Sufficient control is furthermore retained via parameter specification and the option to include expert knowledge. Conversely, completely automatic operation also allows for extensive analysis of aspects related to landscape discretisation.
In a case study, the application of the package is presented. A sensitivity analysis of the most important discretisation parameters demonstrates its efficient workflow automation. Considering multiple streamflow metrics, the employed model proved reasonably robust to the discretisation parameters. However, parameters determining the sizes of subbasins and hillslopes proved to be more important than the others, including the number of representative hillslopes, the number of attributes employed for the lumping algorithm, and the number of sub-discretisations of the representative hillslopes.
The Limpopo Basin in southern Africa is prone to droughts which affect the livelihood of millions of people in South Africa, Botswana, Zimbabwe and Mozambique. Seasonal drought early warning is thus vital for the whole region. In this study, the predictability of hydrological droughts during the main runoff period from December to May is assessed using statistical approaches. Three methods (multiple linear models, artificial neural networks, random forest regression trees) are compared in terms of their ability to forecast streamflow with up to 12 months of lead time. The following four main findings result from the study. 1. There are stations in the basin at which standardised streamflow is predictable with lead times up to 12 months. The results show high inter-station differences of forecast skill but reach a coefficient of determination as high as 0.73 (cross validated). 2. A large range of potential predictors is considered in this study, comprising well-established climate indices, customised teleconnection indices derived from sea surface temperatures and antecedent streamflow as a proxy of catchment conditions. El Nino and customised indices, representing sea surface temperature in the Atlantic and Indian oceans, prove to be important teleconnection predictors for the region. Antecedent streamflow is a strong predictor in small catchments (with median 42% explained variance), whereas teleconnections exert a stronger influence in large catchments. 3. Multiple linear models show the best forecast skill in this study and the greatest robustness compared to artificial neural networks and random forest regression trees, despite their capabilities to represent nonlinear relationships. 4. Employed in early warning, the models can be used to forecast a specific drought level. Even if the coefficient of determination is low, the forecast models have a skill better than a climatological forecast, which is shown by analysis of receiver operating characteristics (ROCs). Seasonal statistical forecasts in the Limpopo show promising results, and thus it is recommended to employ them as complementary to existing forecasts in order to strengthen preparedness for droughts.
Flood events can be expressed by a variety of characteristics such as flood magnitude and extent, event duration or incurred loss. Flood estimation and management may benefit from understanding how the different flood characteristics relate to the hydrological catchment conditions preceding the event and to the meteorological conditions throughout the event. In this study, we therefore propose a methodology to investigate the hydro-meteorological controls on different flood characteristics, based on the simulation of the complete flood risk chain from the flood triggering precipitation event, through runoff generation in the catchment, flood routing and possible inundation in the river system and floodplains to flood loss. Conditional cumulative distribution functions and regression tree analysis delineate the seasonal varying flood processes and indicate that the effect of the hydrological pre-conditions, i.e. soil moisture patterns, and of the meteorological conditions, i.e. weather patterns, depends on the considered flood characteristic. The methodology is exemplified for the Elbe catchment. In this catchment, the length of the build-up period, the event duration and the number of gauges undergoing at least a 10-year flood are governed by weather patterns. The affected length and the number of gauges undergoing at least a 2-year flood are however governed by soil moisture patterns. In case of flood severity and loss, the controlling factor is less pronounced. Severity is slightly governed by soil moisture patterns whereas loss is slightly governed by weather patterns. The study highlights that flood magnitude and extent arise from different flood generation processes and concludes that soil moisture patterns as well as weather patterns are not only beneficial to inform on possible flood occurrence but also on the involved flood processes and resulting flood characteristics.
Long-term estimates of the dissolution of silicate rock are generally derived from a range of isotopic proxies, such as the radiogenic strontium isotope ratio (Sr-87/Sr-86), which are preserved in sediment archives. For these systems to fairly represent silicate weathering, the changes in isotopic ratios in terrestrial surface waters should correspond to changes in the overall silicate dissolution. This assumes that the silicate mineral phases that act as sources of a given isotope dissolve at a rate that is proportional to the overall silicate weathering. Bedrock landsliding exhumes large quantities of fresh rock for weathering in transient storage, and rapid weathering in these deposits is controlled primarily by dissolution of the most reactive phases. In this study, we test the hypothesis that preferential weathering of these labile minerals can decouple the dissolution of strontium sources from the actual silicate weathering rates in the rapidly eroding Western Southern Alps (WSA) of New Zealand. We find that rapid dissolution of relatively radiogenic calcite and biotite in landslides leads to high local fluxes in strontium with isotopic ratios that offer no clear discrimination between sources. These higher fluxes of radiogenic strontium are in contrast to silicate weathering rates in landslides that are not systematically elevated. On a mountain belt scale, radiogenic strontium fluxes are not coupled to volumes of recent landslides in large (>100 km(2)) catchments, but silicate weathering fluxes are. Such decoupling is likely due first to the broad variability in the strontium content of carbonate minerals, and second to the combination of radiogenic strontium released from both biotite and carbonate in recent landslides. This study supports previous work suggesting the limited utility of strontium isotopes as a system to study silicate weathering in the WSA. Crucially however, in settings where bedrock landsliding is a dominant erosive process there is potential for both random and systematic bias in isotope proxies if the most reactive phases exposed for dissolution by landslides disproportionately contribute to the proxy of choice. This clearly suggests that the isotopic composition of marine Sr is a proxy for periods of rapid mountain uplift and erosion rather than for the associated enhanced silicate weathering. (C) 2017 Elsevier Ltd. All rights reserved.
Erosion and tectonic uplift are widely thought to be coupled through feedbacks involving orographic precipitation, relief development, and crustal weakening. In many orogenic systems, it can be difficult to distinguish whether true feedbacks exist, or whether observed features are a consequence of tectonic forcing. To help elucidate these interactions, we examine Gongga Shan, a 7556 m peak on the eastern margin of the Tibetan Plateau where cosmogenic Be-10 basin-wide erosion rates reach >5 mm/yr, defining a region of localized rapid erosion associated with a restraining bend in the left-lateral Xianshuihe Fault. Erosion rates are consistent with topography, thermochronometry, and geodetic data, suggesting a stable pattern of uplift and exhumation over at least the past 2-3 My. Transpression along the Xianshuihe Fault, orographically enhanced precipitation, thermally weakened crust, and substantial local relief all developed independently in the Gongga region and existed there prior to the uplift of Gongga Shan. However, only where all of these conditions are present do the observed topographic and erosional extremes exist, and their relative timing indicates that these conditions are not a consequence of rapid uplift. We conclude that their collocation at 3-4 Ma set into motion a series of feedbacks between erosion and uplift that has resulted in the exceptionally high topography and rapid erosion rates observed today.
Mountain channels can be strongly coupled with adjacent hillslopes, exchanging both mass and energy. However, hypotheses of the underlying cause and effect relations are based on indirect observations that do not resolve the mechanics of channel-hillslope coupling at the process scale. Here we present direct observational data of a coupled channel-hillslope system in the catchment area of the Erlenbach, a mountain stream in Switzerland. A slow-moving landslide flanking this alpine stream failed after a flood had eroded an alluvial step in the channel at its base, representing evidence for an upsystem link in channel-hillslope coupling. Progressive accumulation of landslide debris in the channel resulted in a renewed step, stabilizing the hillslope and restoring the channel step in a downsystem link. Thus, upsystem and downsystem coupling mechanisms are joined in a negative feedback cycle. In this cycle, debuttressing and rebuttressing due to channel bed erosion and alluviation are the dominant controls on hillslope stability. Based on an order of magnitude estimate it is plausible that the observed feedback mechanism is a relevant process in the production of coarse (>2 mm) sediment in the Erlenbach.
Permafrost is a distinct feature of the terrestrial Arctic and is vulnerable to climate warming. Permafrost degrades in different ways, including deepening of a seasonally unfrozen surface and localized but rapid development of deep thaw features. Pleistocene ice-rich permafrost with syngenetic ice-wedges, termed Yedoma deposits, are widespread in Siberia, Alaska, and Yukon, Canada and may be especially prone to rapid-thaw processes. Freeze-locked organic matter in such deposits can be re-mobilized on short time-scales and contribute to a carbon-cycle climate feedback. Here we synthesize the characteristics and vulnerability of Yedoma deposits by synthesizing studies on the Yedoma origin and the associated organic carbon pool. We suggest that Yedoma deposits accumulated under periglacial weathering, transport, and deposition dynamics in non-glaciated regions during the late Pleistocene until the beginning of late glacial warming. The deposits formed due to a combination of aeolian, colluvial, nival, and alluvial deposition and simultaneous ground ice accumulation. We found up to 130 gigatons organic carbon in Yedoma, parts of which are well-preserved and available for fast decomposition after thaw. Based on incubation experiments, up to 10% of the Yedoma carbon is considered especially decomposable and may be released upon thaw. The substantial amount of ground ice in Yedoma makes it highly vulnerable to disturbances such as thermokarst and thermo-erosion processes. Mobilization of permafrost carbon is expected to increase under future climate warming. Our synthesis results underline the need of accounting for Yedoma carbon stocks in next generation Earth-System-Models for a more complete representation of the permafrost-carbon feedback.
Presentation and Analysis of a Worldwide Database of Earthquake-Induced Landslide Inventories
(2017)
Earthquake-induced landslide (EQIL) inventories are essential tools to extend our knowledge of the relationship between earthquakes and the landslides they can trigger. Regrettably, such inventories are difficult to generate and therefore scarce, and the available ones differ in terms of their quality and level of completeness. Moreover, access to existing EQIL inventories is currently difficult because there is no centralized database. To address these issues, we compiled EQIL inventories from around the globe based on an extensive literature study. The database contains information on 363 landslide-triggering earthquakes and includes 66 digital landslide inventories. To make these data openly available, we created a repository to host the digital inventories that we have permission to redistribute through the U.S. Geological Survey ScienceBase platform. It can grow over time as more authors contribute their inventories. We analyze the distribution of EQIL events by time period and location, more specifically breaking down the distribution by continent, country, and mountain region. Additionally, we analyze frequency distributions of EQIL characteristics, such as the approximate area affected by landslides, total number of landslides, maximum distance from fault rupture zone, and distance from epicenter when the fault plane location is unknown. For the available digital EQIL inventories, we examine the underlying characteristics of landslide size, topographic slope, roughness, local relief, distance to streams, peak ground acceleration, peak ground velocity, and Modified Mercalli Intensity. Also, we present an evaluation system to help users assess the suitability of the available inventories for different types of EQIL studies and model development.
Prediction of the area affected by earthquake-induced landsliding based on seismological parameters
(2017)
We present an analytical, seismologically consistent expression for the surface area of the region within which most landslides triggered by an earthquake are located (landslide distribution area). This expression is based on scaling laws relating seismic moment, source depth, and focal mechanism with ground shaking and fault rupture length and assumes a globally constant threshold of acceleration for onset of systematic mass wasting. The seismological assumptions are identical to those recently used to propose a seismologically consistent expression for the total volume and area of landslides triggered by an earthquake. To test the accuracy of the model we gathered geophysical information and estimates of the landslide distribution area for 83 earthquakes. To reduce uncertainties and inconsistencies in the estimation of the landslide distribution area, we propose an objective definition based on the shortest distance from the seismic wave emission line containing 95% of the total landslide area. Without any empirical calibration the model explains 56% of the variance in our dataset, and predicts 35 to 49 out of 83 cases within a factor of 2, depending on how we account for uncertainties on the seismic source depth. For most cases with comprehensive landslide inventories we show that our prediction compares well with the smallest region around the fault containing 95% of the total landslide area. Aspects ignored by the model that could explain the residuals include local variations of the threshold of acceleration and processes modulating the surface ground shaking, such as the distribution of seismic energy release on the fault plane, the dynamic stress drop, and rupture directivity. Nevertheless, its simplicity and first-order accuracy suggest that the model can yield plausible and useful estimates of the landslide distribution area in near-real time, with earthquake parameters issued by standard detection routines.
Soil properties show high heterogeneity at different spatial scales and their correct characterization remains a crucial challenge over large areas. The aim of the study is to quantify the impact of different types of uncertainties that arise from the unresolved soil spatial variability on simulated hydrological states and fluxes. Three perturbation methods are presented for the characterization of uncertainties in soil properties. The methods are applied on the soil map of the upper Neckar catchment (Germany), as an example. The uncertainties are propagated through the distributed mesoscale hydrological model (mHM) to assess the impact on the simulated states and fluxes. The model outputs are analysed by aggregating the results at different spatial and temporal scales. These results show that the impact of the different uncertainties introduced in the original soil map is equivalent when the simulated model outputs are analysed at the model grid resolution (i.e. 500 m). However, several differences are identified by aggregating states and fluxes at different spatial scales (by subcatchments of different sizes or coarsening the grid resolution). Streamflow is only sensitive to the perturbation of long spatial structures while distributed states and fluxes (e.g. soil moisture and groundwater recharge) are only sensitive to the local noise introduced to the original soil properties. A clear identification of the temporal and spatial scale for which finer-resolution soil information is (or is not) relevant is unlikely to be universal. However, the comparison of the impacts on the different hydrological components can be used to prioritize the model improvements in specific applications, either by collecting new measurements or by calibration and data assimilation approaches. In conclusion, the study underlines the importance of a correct characterization of uncertainty in soil properties. With that, soil maps with additional information regarding the unresolved soil spatial variability would provide strong support to hydrological modelling applications.
Seismic-hazard assessment is of great importance within the field of engineering seismology. Nowadays, it is common practice to define future seismic demands using probabilistic seismic-hazard analysis (PSHA). Often it is neither obvious nor transparent how PSHA responds to changes in its inputs. In addition, PSHA relies on many uncertain inputs. Sensitivity analysis (SA) is concerned with the assessment and quantification of how changes in the model inputs affect the model response and how input uncertainties influence the distribution of the model response. Sensitivity studies are challenging primarily for computational reasons; hence, the development of efficient methods is of major importance. Powerful local (deterministic) methods widely used in other fields can make SA feasible, even for complex models with a large number of inputs; for example, automatic/algorithmic differentiation (AD)-based adjoint methods. Recently developed derivative-based global sensitivity measures can combine the advantages of such local SA methods with efficient sampling strategies facilitating quantitative global sensitivity analysis (GSA) for complex models. In our study, we propose and implement exactly this combination. It allows an upper bounding of the sensitivities involved in PSHA globally and, therefore, an identification of the noninfluential and the most important uncertain inputs. To the best of our knowledge, it is the first time that derivative-based GSA measures are combined with AD in practice. In addition, we show that first-order uncertainty propagation using the delta method can give satisfactory approximations of global sensitivity measures and allow a rough characterization of the model output distribution in the case of PSHA. An illustrative example is shown for the suggested derivative-based GSA of a PSHA that uses stochastic ground-motion simulations.