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.
The aim of this paper is to investigate the ability of various site-condition proxies (SCPs) to reduce ground-motion aleatory variability and evaluate how SCPs capture nonlinearity site effects. The SCPs used here are time-averaged shear-wave velocity in the top 30 m (V-S30), the topographical slope (slope), the fundamental resonance frequency (f(0)) and the depth beyond which V-s exceeds 800 m/s (H800). We considered first the performance of each SCP taken alone and then the combined performance of the 6 SCP pairs [V-S30-f(0)], [V-S30-H-800], [f(0)-slope], [H-800-slope], [V-S30-slope] and [f(0)-H-800]. This analysis is performed using a neural network approach including a random effect applied on a KiK-net subset for derivation of ground-motion prediction equations setting the relationship between various ground-motion parameters such as peak ground acceleration, peak ground velocity and pseudo-spectral acceleration PSA (T), and Mw, RJB, focal depth and SCPs. While the choice of SCP is found to have almost no impact on the median groundmotion prediction, it does impact the level of aleatory uncertainty. VS30 is found to perform the best of single proxies at short periods (T < 0.6 s), while f(0) and H-800 perform better at longer periods; considering SCP pairs leads to significant improvements, with particular emphasis on [V-S30-H-800] and [f(0)-slope] pairs. The results also indicate significant nonlinearity on the site terms for soft sites and that the most relevant loading parameter for characterising nonlinear site response is the "stiff" spectral ordinate at the considered period.
In hydrological models, parameters are used to represent the time-invariant characteristics of catchments and to capture different aspects of hydrological response. Hence, model parameters need to be identified based on their role in controlling the hydrological behaviour. For the identification of meaningful parameter values, multiple and complementary performance criteria are used that compare modelled and measured discharge time series. The reliability of the identification of hydrologically meaningful model parameter values depends on how distinctly a model parameter can be assigned to one of the performance criteria.& para;& para;To investigate this, we introduce the new concept of connective strength between model parameters and performance criteria. The connective strength assesses the intensity in the interrelationship between model parameters and performance criteria in a bijective way. In our analysis of connective strength, model simulations are carried out based on a latin hypercube sampling. Ten performance criteria including Nash-Sutcliffe efficiency (NSE), Kling-Gupta efficiency (KGE) and its three components (alpha, beta and r) as well as RSR (the ratio of the root mean square error to the standard deviation) for different segments of the flow duration curve (FDC) are calculated.& para;& para;With a joint analysis of two regression tree (RT) approaches, we derive how a model parameter is connected to different performance criteria. At first, RTs are constructed using each performance criterion as the target variable to detect the most relevant model parameters for each performance criterion. Secondly, RTs are constructed using each parameter as the target variable to detect which performance criteria are impacted by changes in the values of one distinct model parameter. Based on this, appropriate performance criteria are identified for each model parameter.& para;& para;In this study, a high bijective connective strength between model parameters and performance criteria is found for low- and mid-flow conditions. Moreover, the RT analyses emphasise the benefit of an individual analysis of the three components of KGE and of the FDC segments. Furthermore, the RT analyses highlight under which conditions these performance criteria provide insights into precise parameter identification. Our results show that separate performance criteria are required to identify dominant parameters on low- and mid-flow conditions, whilst the number of required performance criteria for high flows increases with increasing process complexity in the catchment. Overall, the analysis of the connective strength between model parameters and performance criteria using RTs contribute to a more realistic handling of parameters and performance criteria in hydrological modelling.
Reliable flood risk analyses, including the estimation of damage, are an important prerequisite for efficient risk management. However, not much is known about flood damage processes affecting companies. Thus, we conduct a flood damage assessment of companies in Germany with regard to two aspects. First, we identify relevant damage-influencing variables. Second, we assess the prediction performance of the developed damage models with respect to the gain by using an increasing amount of training data and a sector-specific evaluation of the data. Random forests are trained with data from two postevent surveys after flood events occurring in the years 2002 and 2013. For a sector-specific consideration, the data set is split into four subsets corresponding to the manufacturing, commercial, financial, and service sectors. Further, separate models are derived for three different company assets: buildings, equipment, and goods and stock. Calculated variable importance values reveal different variable sets relevant for the damage estimation, indicating significant differences in the damage process for various company sectors and assets. With an increasing number of data used to build the models, prediction errors decrease. Yet the effect is rather small and seems to saturate for a data set size of several hundred observations. In contrast, the prediction improvement achieved by a sector-specific consideration is more distinct, especially for damage to equipment and goods and stock. Consequently, sector-specific data acquisition and a consideration of sector-specific company characteristics in future flood damage assessments is expected to improve the model performance more than a mere increase in data.
Voluminous magmatism during the South Atlantic opening has been considered as a classical example for plume related continental breakup. We present a study of the crustal structure around Walvis Ridge, near the intersection with the African margin. Two wide-angle seismic profiles were acquired. One is oriented NNW–SSE, following the continent–ocean transition and crossing Walvis Ridge. A second amphibious profile runs NW–SE from the Angola Basin into continental Namibia. At the continent–ocean boundary (COB) the mafic crust beneath Walvis Ridge is up to 33 km thick, with a pronounced high-velocity lower crustal body. Towards the south there is a smooth transition to 20–25 km thick crust underlying the COB in the Walvis Basin, with a similar velocity structure, indicating a gabbroic lower crust with associated cumulates at the base. The northern boundary of Walvis Ridge towards the Angola Basin shows a sudden change to oceanic crust only 4–6 km thick, coincident with the projection of the Florianopolis Fracture Zone, one of the most prominent tectonic features of the South Atlantic ocean basin. In the amphibious profile the COB is defined by a sharp transition from oceanic to rifted continental crust, with a magmatic overprint landward of the intersection of Walvis Ridge with the Namibian margin. The continental crust beneath the Congo Craton is 40 km thick, shoaling to 35 km further SE. The velocity models show that massive high-velocity gabbroic intrusives are restricted to a narrow zone directly underneath Walvis Ridge and the COB in the south. This distribution of rift-related magmatism is not easily reconciled with models of continental breakup following the establishment of a large, axially symmetric plume in the Earth's mantle. Rift-related lithospheric stretching and associated transform faulting play an overriding role in locating magmatism, dividing the margin in a magma-dominated southern and an essentially amagmatic northern segment.
Vast areas of the terrestrial Subarctic and Arctic are underlain by permafrost. Landscape evolution is therefore largely controlled by climate-driven periglacial processes. The response of the frozen ground to late Quaternary warm and cold stages is preserved in permafrost sequences, and deducible by multi-proxy palaeoenvironmental approaches. Here, we analyse radiocarbon-dated mid-Wisconsin Interstadial and Holocene lacustrine deposits preserved in the Kit-1 pingo permafrost sequence combined with water and surface sediment samples from nine modern water bodies on Seward Peninsula (NW Alaska) to reconstruct thermokarst dynamics and determine major abiotic factors that controlled the aquatic ecosystem variability. Our methods comprise taxonomical diatom analyses as well as Detrended Correspondence Analysis (DCA) and Redundancy Analysis (RDA). Our results show, that the fossil diatom record reflects thermokarst lake succession since about 42 C-14 kyr BP. Different thermolcarst lake stages during the mid-Wisconsin Interstadial, the late Wisconsin and the early Holocene are mirrored by changes in diatom abundance, diversity, and ecology. We interpret the taxonomical changes in the fossil diatom assemblages in combination with both modern diatom data from surrounding ponds and existing micropalaeontological, sedimentological and mineralogical data from the pingo sequence. A diatom based quantitative reconstruction of lake water pH indicates changing lake environments during mid-Wisconsin to early Holocene stages. Mineralogical analyses indicate presence of tephra fallout and its impact on fossil diatom communities. Our comparison of modern and fossil diatom communities shows the highest floristic similarity of modern polygon ponds to the corresponding initial (shallow water) development stages of thermolcarst lakes. We conclude, that mid-Wisconsin thermokarst processes in the study area could establish during relatively warm interstadial climate conditions accompanied by increased precipitation due to approaching coasts, while still high continentality and hence high seasonal temperature gradients led to warm summers in the central part of Beringia. (C) 2017 Elsevier B.V. All rights reserved.
We investigated a well-dated sediment section of a palaeolake situated in the coastal zone of Shikotan Island (Lesser Kurils) for organic sediment-geochemistry and biotic components (diatoms, chironomids, pollen) in order to provide a reconstruction of the palaeoenvironmental changes and palaeo-events (tsunamis, sea-level fluctuations and landslides) in Holocene. During the ca 8000 years of sedimentation the changes in organic sediment-geochemistry and in composition of the diatoms and chironomids as well as the shifts in composition of terrestrial vegetation suggest that the period until ca 5800 cal yr BP was characterized by a warm and humid climate (corresponds to middle Holocene optimum) with climate cooling thereafter. A warm period reconstructed from ca 900 to at least ca 580 cal yr BP corresponds to a transition to a Nara-Heian-Kamakura warm stage and can be correlated to a Medieval Warm Period. After 580 cal yr PB, the lake gradually dried out and climatic signals could not be obtained from the declining lacustrine biological communities, but the increasing role of spruce and disappearance of the oak from the vegetation give evidences of the climate cooling that can be correlated with the LIA. The marine regression stages at the investigated site are identified for ca 6200-5900 (at the end of the middle Holocene transgression), ca 5500-5100 (Middle Jomon regression or Kemigawa regression), and ca 1070-360 cal yr BP (at the end of Heian transgression). The lithological structure of sediments and the diatom compositions give evidences for the multiple tsunami events of different strengths in the Island. Most remarkable of them can be dated at around ca 7000, 6460, 5750, 4800, 950 cal yr BP. The new results help to understand the Holocene environmental history of the Southern Kurils as a part of the Kuril-Kamchatka and Aleutian Marginal Sea-Island Arc Systems in the North-Western Pacific region.
Tree stands in the boreal treeline ecotone are, in addition to climate change, impacted by disturbances such as fire, water-related disturbances and logging. We aim to understand how these disturbances affect growth, age structure, and spatial patterns of larch stands in the north-eastern Siberian treeline ecotone (lower Kolyma River region), an insufficiently researched region. Stand structure of Larix cajanderi Mayr was studied at seven sites impacted by disturbances. Maximum tree age ranged from 44 to 300 years. Young to medium-aged stands had, independent of disturbance type, the highest stand densities with over 4000 larch trees per ha. These sites also had the highest growth rates for tree height and stem diameter. Overall lowest stand densities were found in a polygonal field at the northern end of the study area, with larches growing in distinct " tree islands". At all sites, saplings are significantly clustered. Differences in fire severity led to contrasting stand structures with respect to tree, recruit, and overall stand densities. While a low severity fire resulted in low-density stands with high proportions of small and young larches, high severity fires resulted in high-density stands with high proportions of big trees. At waterdisturbed sites, stand structure varied between waterlogged and drained sites and latitude. These mixed effects of climate and disturbance make it difficult to predict future stand characteristics and the treeline position.
This article presents a review of the current data on the level of paleolimnological knowledge about lakes in the Russian part of the northern Eurasia. The results of investigation of the northwestern European part of Russia as the best paleolimnologically studied sector of the Russian north is presented in detail. The conditions of lacustrine sedimentation at the boundary between the Late Pleistocene and Holocene and the role of different external factors in formation of their chemical composition, including active volcanic activity and possible large meteorite impacts, are also discussed. The results of major paleoclimatic and paleoecological reconstructions in northern Siberia are presented. Particular attention is given to the databases of abiotic and biotic parameters of lake ecosystems as an important basis for quantitative reconstructions of climatic and ecological changes in the Late Pleistocene and Holocene. Keywords: paleolimnology, lakes, bottom sediments, northern.
Adaptation to flood risk
(2017)
As flood impacts are increasing in large parts of the world, understanding the primary drivers of changes in risk is essential for effective adaptation. To gain more knowledge on the basis of empirical case studies, we analyze eight paired floods, that is, consecutive flood events that occurred in the same region, with the second flood causing significantly lower damage. These success stories of risk reduction were selected across different socioeconomic and hydro-climatic contexts. The potential of societies to adapt is uncovered by describing triggered societal changes, as well as formal measures and spontaneous processes that reduced flood risk. This novel approach has the potential to build the basis for an international data collection and analysis effort to better understand and attribute changes in risk due to hydrological extremes in the framework of the IAHSs Panta Rhei initiative. Across all case studies, we find that lower damage caused by the second event was mainly due to significant reductions in vulnerability, for example, via raised risk awareness, preparedness, and improvements of organizational emergency management. Thus, vulnerability reduction plays an essential role for successful adaptation. Our work shows that there is a high potential to adapt, but there remains the challenge to stimulate measures that reduce vulnerability and risk in periods in which extreme events do not occur.
Trends in precipitation over Germany and the Rhine basin related to changes in weather patterns
(2017)
Precipitation as the central meteorological feature for agriculture, water security, and human well-being amongst others, has gained special attention ever since. Lack of precipitation may have devastating effects such as crop failure and water scarcity. Abundance of precipitation, on the other hand, may as well result in hazardous events such as flooding and again crop failure. Thus, great effort has been spent on tracking changes in precipitation and relating them to underlying processes. Particularly in the face of global warming and given the link between temperature and atmospheric water holding capacity, research is needed to understand the effect of climate change on precipitation.
The present work aims at understanding past changes in precipitation and other meteorological variables. Trends were detected for various time periods and related to associated changes in large-scale atmospheric circulation. The results derived in this thesis may be used as the foundation for attributing changes in floods to climate change. Assumptions needed for the downscaling of large-scale circulation model output to local climate stations are tested and verified here.
In a first step, changes in precipitation over Germany were detected, focussing not only on precipitation totals, but also on properties of the statistical distribution, transition probabilities as a measure for wet/dry spells, and extreme precipitation events.
Shifting the spatial focus to the Rhine catchment as one of the major water lifelines of Europe and the largest river basin in Germany, detected trends in precipitation and other meteorological variables were analysed in relation to states of an ``optimal'' weather pattern classification. The weather pattern classification was developed seeking the best skill in explaining the variance of local climate variables.
The last question addressed whether observed changes in local climate variables are attributable to changes in the frequency of weather patterns or rather to changes within the patterns itself. A common assumption for a downscaling approach using weather patterns and a stochastic weather generator is that climate change is expressed only as a changed occurrence of patterns with the pattern properties remaining constant. This assumption was validated and the ability of the latest generation of general circulation models to reproduce the weather patterns was evaluated.
% Paper 1
Precipitation changes in Germany in the period 1951-2006 can be summarised briefly as negative in summer and positive in all other seasons. Different precipitation characteristics confirm the trends in total precipitation: while winter mean and extreme precipitation have increased, wet spells tend to be longer as well (expressed as increased probability for a wet day followed by another wet day). For summer the opposite was observed: reduced total precipitation, supported by decreasing mean and extreme precipitation and reflected in an increasing length of dry spells.
Apart from this general summary for the whole of Germany, the spatial distribution within the country is much more differentiated. Increases in winter precipitation are most pronounced in the north-west and south-east of Germany, while precipitation increases are highest in the west for spring and in the south for autumn. Decreasing summer precipitation was observed in most regions of Germany, with particular focus on the south and west.
The seasonal picture, however, was again differently represented in the contributing months, e.g.\ increasing autumn precipitation in the south of Germany is formed by strong trends in the south-west in October and in the south-east in November. These results emphasise the high spatial and temporal organisation of precipitation changes.
% Paper 2
The next step towards attributing precipitation trends to changes in large-scale atmospheric patterns was the derivation of a weather pattern classification that sufficiently stratifies the local climate variables under investigation. Focussing on temperature, radiation, and humidity in addition to precipitation, a classification based on mean sea level pressure, near-surface temperature, and specific humidity was found to have the best skill in explaining the variance of the local variables. A rather high number of 40 patterns was selected, allowing typical pressure patterns being assigned to specific seasons by the associated temperature patterns. While the skill in explaining precipitation variance is rather low, better skill was achieved for radiation and, of course, temperature.
Most of the recent GCMs from the CMIP5 ensemble were found to reproduce these weather patterns sufficiently well in terms of frequency, seasonality, and persistence.
% Paper 3
Finally, the weather patterns were analysed for trends in pattern frequency, seasonality, persistence, and trends in pattern-specific precipitation and temperature. To overcome uncertainties in trend detection resulting from the selected time period, all possible periods in 1901-2010 with a minimum length of 31 years were considered. Thus, the assumption of a constant link between patterns and local weather was tested rigorously. This assumption was found to hold true only partly. While changes in temperature are mainly attributable to changes in pattern frequency, for precipitation a substantial amount of change was detected within individual patterns.
Magnitude and even sign of trends depend highly on the selected time period. The frequency of certain patterns is related to the long-term variability of large-scale circulation modes.
Changes in precipitation were found to be heterogeneous not only in space, but also in time - statements on trends are only valid for the specific time period under investigation. While some part of the trends can be attributed to changes in the large-scale circulation, distinct changes were found within single weather patterns as well.
The results emphasise the need to analyse multiple periods for thorough trend detection wherever possible and add some note of caution to the application of downscaling approaches based on weather patterns, as they might misinterpret the effect of climate change due to neglecting within-type trends.
The Cotentin Peninsula (Normandy, France) displays sequences of marine terraces and rasas, the latter being wide Late Cenozoic coastal erosion surfaces, that are typical of Western European coasts in Portugal, Spain, France and southern England. Remote sensing imagery and field mapping enabled reappraisal of the Cotentin coastal sequences. From bottom to top, the N Cotentin sequence includes four previously recognized Pleistocene marine terraces (T1 to T4) at elevations <40 m as well as four higher and older rasas (R1 to R4) reaching 200 +/- 5 m in elevation. Low-standing marine terraces are not observed in the central part of the Peninsula and a limited number of terraces are described to the south. The high-standing rasas are widespread all over the peninsula. Such strandline distributions reveal major changes during the Late Cenozoic. Progressive uplift of an irregular sea-floor led to subaerial exposure of bathymetric highs that were carved into rocky platforms, rasas and marine terraces. Eventually, five main islands coalesced and connected to the mainland to the south to form the Cotentin Peninsula. On the basis of previous dating of the last interglacial maximum terrace (i.e. Marine Isotopic Stage, MIS 5e), sequential morphostratigraphy and modelling, we have reappraised uplift rates and derived: (i) mean Upper Pleistocene (i.e. since MIS 5e similar to 122 +/- 6 ka, i.e. kilo annum) apparent uplift rates of 0.04 +/- 0.01 mm/yr, (ii) mean Middle Pleistocene eustasy-corrected uplift rates of 0.09 +/- 0.03 mm/yr, and (iii) low mean Pleistocene uplift rates of 0.01 mm/yr. Extrapolations of these slow rates combined with geological evidence implies that the formation of the sequences from the Cotentin Peninsula occurred between 3 Ma (Pliocene) and 15 Ma (Miocene), which cannot be narrowed down further without additional research. Along the coasts of Western Europe, sequences of marine terraces and rasas are widespread (169 preserve the MIS Se benchmark). In Spain, Portugal, S England and other parts of western France, the sequences morphostratigraphy is very similar to that of Cotentin. The onset of such Western European sequences occurred during the Miocene (e.g. Spain) or Pliocene (e.g. Portugal). We interpret this Neogene-Quaternary coastal uplift as a symptom of the increasing lithospheric compression that accompanies Cenozoic orogenies. (C) 2017 Elsevier B.V. All rights reserved.
Pollen-based quantitative reconstructions of past climate variables is a standard palaeoclimatic approach. Despite knowing that the spatial extent of the calibration-set affects the reconstruction result, guidance is lacking as to how to determine a suitable spatial extent of the pollen-climate calibration-set. In this study, past mean annual precipitation (P-ann) during the Holocene (since 11.5 cal ka BP) is reconstructed repeatedly for pollen records from Qinghai Lake (36.7 degrees N, 100.5 degrees E; north-east Tibetan Plateau), Gonghai Lake (38.9 degrees N, 112.2 degrees E; north China) and Sihailongwan Lake (42.3 degrees N, 126.6 degrees E; north-east China) using calibration-sets of varying spatial extents extracted from the modern pollen dataset of China and Mongolia (2559 sampling sites and 168 pollen taxa in total). Results indicate that the spatial extent of the calibration-set has a strong impact on model performance, analogue quality and reconstruction diagnostics (absolute value, range, trend, optimum). Generally, these effects are stronger with the modern analogue technique (MAT) than with weighted averaging partial least squares (WA-PLS). With respect to fossil spectra from northern China, the spatial extent of calibration-sets should be restricted to radii between ca. 1000 and 1500 km because small-scale calibration-sets (<800 km radius) will likely fail to include enough spatial variation in the modern pollen assemblages to reflect the temporal range shifts during the Holocene, while too broad a scale calibration-set (>1500 km radius) will include taxa with very different pollen-climate relationships. (C) 2017 Elsevier Ltd. All rights reserved.
The Indian Summer Monsoon (ISM) is one of the most important climate systems, whose variability and driving mechanisms are of broad interest for academic and societal communities. Here, we present a well-dated high-resolution pollen analysis from a 4.82-m long sediment core taken from Basomtso, in the southeastern Tibetan Plateau (TP), which depicts the regional climate changes of the past millennium. Our results show that subalpine coniferous forest was dominant around Basomtso from ca. 867 to ca. 750 cal. yr BP, indicating a warm and semi-humid climate. The timberline in the study area significantly decreased from ca. 750 to ca.100 cal. yr BP, and a cold climate, corresponding to the Little Ice Age (LIA) prevailed. Since ca. 100 cal. yr BP, the vegetation type changed to forest-meadow with rising temperatures and moisture. Ordination analysis reveals that the migration of vegetation was dominated by regional temperatures and then by moisture. Further comparisons between the Basomtso pollen record and the regional temperature reconstructions underscore the relevance of the Basomtso record from the southeastern TP for regional and global climatologies. Our pollen based moisture reconstruction demonstrates the strong multicentennial-scale link to ISM variability, providing solid evidence for the increase of monsoonal strengths over the past four centuries. Spectral analysis indicates the potential influence of solar forcing. However, a closer relationship has been observed between multicentennial ISM variations and Indian Ocean sea surface temperature anomalies (SSTs), suggesting that the variations in monsoonal precipitation over the southeastern TP are probably driven by the Indian Ocean Dipole on the multicentennial scale. (C) 2017 Elsevier Ltd. All rights reserved.
Early agriculture can be detected in palaeovegetation records, but quantification of the relative importance of climate and land use in influencing regional vegetation composition since the onset of agriculture is a topic that is rarely addressed. We present a novel approach that combines pollen-based REVEALS estimates of plant cover with climate, anthropogenic land-cover and dynamic vegetation modelling results. This is used to quantify the relative impacts of land use and climate on Holocene vegetation at a sub-continental scale, i.e. northern and western Europe north of the Alps. We use redundancy analysis and variation partitioning to quantify the percentage of variation in vegetation composition explained by the climate and land-use variables, and Monte Carlo permutation tests to assess the statistical significance of each variable. We further use a similarity index to combine pollen based REVEALS estimates with climate-driven dynamic vegetation modelling results. The overall results indicate that climate is the major driver of vegetation when the Holocene is considered as a whole and at the sub-continental scale, although land use is important regionally. Four critical phases of land-use effects on vegetation are identified. The first phase (from 7000 to 6500 BP) corresponds to the early impacts on vegetation of farming and Neolithic forest clearance and to the dominance of climate as a driver of vegetation change. During the second phase (from 4500 to 4000 BP), land use becomes a major control of vegetation. Climate is still the principal driver, although its influence decreases gradually. The third phase (from 2000 to 1500 BP) is characterised by the continued role of climate on vegetation as a consequence of late-Holocene climate shifts and specific climate events that influence vegetation as well as land use. The last phase (from 500 to 350 BP) shows an acceleration of vegetation changes, in particular during the last century, caused by new farming practices and forestry in response to population growth and industrialization. This is a unique signature of anthropogenic impact within the Holocene but European vegetation remains climatically sensitive and thus may continue to respond to ongoing climate change. (C) 2017 Elsevier Ltd. All rights reserved.
Climatic and limnological changes at Lake Karakul (Tajikistan) during the last similar to 29 cal ka
(2017)
We present results of analyses on a sediment core from Lake Karakul, located in the eastern Pamir Mountains, Tajikistan. The core spans the last similar to 29 cal ka. We investigated and assessed processes internal and external to the lake to infer changes in past moisture availability. Among the variables used to infer lake-external processes, high values of grain-size end-member (EM) 3 (wide grain-size distribution that reflects fluvial input) and high Sr/Rb and Zr/Rb ratios (coinciding with coarse grain sizes), are indicative of moister conditions. High values in EM1, EM2 (peaks of small grain sizes that reflect long-distance dust transport or fine, glacially derived clastic input) and TiO2 (terrigenous input) are thought to reflect greater influence of dry air masses, most likely of Westerly origin. High input of dust from distant sources, beginning before the Last Glacial Maximum (LGM) and continuing to the late glacial, reflects the influence of dry Westerlies, whereas peaks in fluvial input suggest increased moisture availability. The early to early-middle Holocene is characterised by coarse mean grain sizes, indicating constant, high fluvial input and moister conditions in the region. A steady increase in terrigenous dust and a decrease in fluvial input from 6.6 cal ka BP onwards points to the Westerlies as the predominant atmospheric circulation through to present, and marks a return to drier and even arid conditions in the area. Proxies for productivity (TOC, TOC/TN, TOCBr), redox potential (Fe/Mn) and changes in the endogenic carbonate precipitation (TIC, delta(18) OCarb) indicate changes within the lake. Low productivity characterised the lake from the late Pleistocene until 6.6 cal ka BP, and increased rapidly afterwards. Lake level remained low until the LGM, but water depth increased to a maximum during the late glacial and remained high into the early Holocene. Subsequently, the water level decreased to its present stage. Today the lake system is mainly climatically controlled, but the depositional regime is also driven by internal limnogeological processes.
Arctic and alpine treelines worldwide differ in their reactions to climate change. A northward advance of or densification within the treeline ecotone will likely influence climate-vegetation feedback mechanisms. In our study, which was conducted in the Taimyr Depression in the North Siberian Lowlands, w present a combined field-and model-based approach helping us to better understand the population processes involved in the responses of the whole treeline ecotone, spanning from closed forest to single-tree tundra, to climate warming. Using information on stand structure, tree age, and seed quality and quantity from seven sites, we investigate effects of intra-specific competition and seed availability on the specific impact of recent climate warming on larch stands. Field data show that tree density is highest in the forest-tundra, and average tree size decreases from closed forest to single-tree tundra. Age-structure analyses indicate that the trees in the closed forest and forest-tundra have been present for at least similar to 240 yr. At all sites except the most southerly ones, past establishment is positively correlated with regional temperature increase. In the single-tree tundra, however, a change in growth form from krummholz to erect trees, beginning similar to 130 yr ago, rather than establishment date has been recorded. Seed mass decreases from south to north, while seed quantity increases. Simulations with LAVESI (Larix Vegetation Simulator) further suggest that relative density changes strongly in response to a warming signal in the forest-tundra while intra-specific competition limits densification in the closed forest and seed limitation hinders densification in the single-tree tundra. We find striking differences in strength and timing of responses to recent climate warming. While forest-tundra stands recently densified, recruitment is almost non-existent at the southern and northern end of the ecotone due to autecological processes. Palaeo-treelines may therefore be inappropriate to infer past temperature changes at a fine scale. Moreover, a lagged treeline response to past warming will, via feedback mechanisms, influence climate change in the future.
The concept of a Global Monsoon (GM) has been proposed based on modern precipitation observations, but its application over a wide range of temporal scales is still under debate. Here, we present a synthesis of 268 continental paleo-moisture records collected from monsoonal systems in the Eastern Hemisphere, including the East Asian Monsoon (EAsM), the Indian Monsoon (IM), the East African Monsoon (EAfM), and the Australian Monsoon (AuM) covering the last 18,000 years. The overall pattern of late Glacial to Holocene moisture change is consistent with those inferred from ice cores and marine records. With respect to the last 10,000 years (10 ka), i.e. a period that has high spatial coverage, a Fuzzy c-Means clustering analysis of the moisture index records together with "Xie-Beni" index reveals four clusters of our data set. The paleoclimatic meaning of each cluster is interpreted considering the temporal evolution and spatial distribution patterns. The major trend in the tropical AuM, EAfM, and IM regions is a gradual decrease in moisture conditions since the early Holocene. Moisture changes in the EAsM regions show maximum index values between 8 and 6 ka. However, records located in nearby subtropical areas, i.e. in regions not influenced by the intertropical convergence zone, show an opposite trend compared to the tropical monsoon regions (AuM, EAfM and IM), i.e. a gradual increase. Analyses of modern meteorological data reveal the same spatial patterns as in the paleoclimate records such that, in times of overall monsoon strengthening, lower precipitation rates are observed in the nearby subtropical areas. We explain this pattern as the effect of a strong monsoon circulation suppressing air uplift in nearby subtropical areas, and hence hindering precipitation. By analogy to the modern system, this would mean that during the early Holocene strong monsoon period, the intensified ascending airflows within the monsoon domains led to relatively weaker ascending or even descending airflows in the adjacent subtropical regions, resulting in a precipitation deficit compared to the late Holocene. Our conceptual model therefore integrates regionally contrasting moisture changes into the Global Monsoon hypothesis. (C) 2017 Elsevier Ltd. All rights reserved.
Retrogressive thaw slumps (RTSs) are among the most active landforms in the Arctic; their number has increased significantly over the past decades. While processes initiating discrete RTSs are well identified, the major terrain controls on the development of coastal RTSs at a regional scale are not yet defined. Our research reveals the main geomorphic factors that determine the development of RTSs along a 238km segment of the Yukon Coast, Canada. We (1) show the current extent of RTSs, (2) ascertain the factors controlling their activity and initiation, and (3) explain the spatial differences in the density and areal coverage of RTSs. We mapped and classified 287 RTSs using high-resolution satellite images acquired in 2011. We highlighted the main terrain controls over their development using univariate regression trees model. Coastal geomorphology influenced both the activity and initiation of RTSs: active RTSs and RTSs initiated after 1972 occurred primarily on terrains with slope angles greater than 3.9 degrees and 5.9 degrees, respectively. The density and areal coverage of RTSs were constrained by the volume and thickness of massive ice bodies. Differences in rates of coastal change along the coast did not affect the model. We infer that rates of coastal change averaged over a 39year period are unable to reflect the complex relationship between RTSs and coastline dynamics. We emphasize the need for large-scale studies of RTSs to evaluate their impact on the ecosystem and to measure their contribution to the global carbon budget. Plain Language Summary Retrogressive thaw slumps, henceforth slumps are a type of landslides that occur when permafrost thaws. Slumps are active landforms: they develop quickly and extend over several hectares. Satellite imagery allows to map such slumps over large areas. Our research shows where slumps develop along a 238 km segment of the Yukon Coast in Canada and explains which environments are most suitable for slump occurrence. We found that active and newly developed slumps were triggered where coastal slopes were greater than 3.9 degrees and 5.9 degrees, respectively. We explain that coastal erosion influences the development of slumps by modifying coastal slopes. We found that the highest density of slumps as well as the largest slumps occurred on terrains with high amounts of ice bodies in the ground. This study provides tools to better identify areas in the Arctic that are prone to slump development.
Radiocarbon and optically stimulated luminescence dating of sediments from Lake Karakul, Tajikistan
(2017)
Lake Karakul in the eastern Pamirs is a large and closed-basin lake in a partly glaciated catchment. Two parallel sediment cores were collected from 12 m water depth. The cores were correlated using XRF analysis and dated using radiocarbon and OSL techniques. The age results of the two dating methods are generally in agreement. The correlated composite core of 12.26 m length represents continuous accumulation of sediments in the lake basin since 31 ka. The lake reservoir effect (LRE) remained relatively constant over this period. High sediment accumulation rates (SedARs) were recorded before 23 ka and after 6.5 ka. The relatively close position of the coring location near the eastern shore of the lake implies that high SedARs resulted from low lake levels. Thus, high SedARs and lower lake levels before 23 ka probably reflect cold and dry climate conditions that inhibited the arrival of moist air at high elevation in the eastern Pamirs. Low lake levels after 6.5 ka were probably caused by declining temperatures after the warmer early Holocene, which had caused a reduction in water resources stored as snow, ice and frozen ground in the catchment. Low SedARs during 23-6.5 ka suggest increased lake levels in Lake Karakul. A short-lived increase of SedARs at 15 ka probably corresponds to the rapid melting of glaciers in the Karakul catchment during the Greenland Interstadial le, shortly after glaciers in the catchment had reached their maximum extents. The sediment cores from Lake Karakul represent an important climate archive with robust chronology for the last glacial interglacial cycle from Central Asia. (C) 2017 Elsevier B.V. All rights reserved.
We analyzed chlorophyll-a and Colored Dissolved Organic Matter (CDOM) dynamics from field measurements and assessed the potential of multispectral satellite data for retrieving water-quality parameters in three small surface reservoirs in the Brazilian semiarid region. More specifically, this work is comprised of: (i) analysis of Chl-a and trophic dynamics; (ii) characterization of CDOM; (iii) estimation of Chl-a and CDOM from OLI/Landsat-8 and RapidEye imagery. The monitoring lasted 20 months within a multi-year drought, which contributed to water-quality deterioration. Chl-a and trophic state analysis showed a highly eutrophic status for the perennial reservoir during the entire study period, while the non-perennial reservoirs ranged from oligotrophic to eutrophic, with changes associated with the first events of the rainy season. CDOM characterization suggests that the perennial reservoir is mostly influenced by autochthonous sources, while allochthonous sources dominate the non-perennial ones. Spectral-group classification assigned the perennial reservoir as a CDOM-moderate and highly eutrophic reservoir, whereas the non-perennial ones were assigned as CDOM-rich and oligotrophic-dystrophic reservoirs. The remote sensing initiative was partially successful: the Chl-a was best modelled using RapidEye for the perennial one; whereas CDOM performed best with Landsat-8 for non-perennial reservoirs. This investigation showed potential for retrieving water quality parameters in dry areas with small reservoirs.
The large variety of atmospheric circulation systems affecting the eastern Asian climate is reflected by the complex Asian vegetation distribution. Particularly in the transition zones of these circulation systems, vegetation is supposed to be very sensitive to climate change. Since proxy records are scarce, hitherto a mechanistic understanding of the past spatio-temporal climate-vegetation relationship is lacking. To assess the Holocene vegetation change and to obtain an ensemble of potential mid-Holocene biome distributions for eastern Asia, we forced the diagnostic biome model BIOME4 with climate anomalies of different transient Holocene climate simulations performed in coupled atmosphere-ocean(-vegetation) models. The simulated biome changes are compared with pollen-based biome records for different key regions. In all simulations, substantial biome shifts during the last 6000 years are confined to the high northern latitudes and the monsoon-westerly wind transition zone, but the temporal evolution and amplitude of change strongly depend on the climate forcing. Large parts of the southern tundra are replaced by taiga during the mid-Holocene due to a warmer growing season and the boreal treeline in northern Asia is shifted northward by approx. 4 degrees in the ensemble mean, ranging from 1.5 to 6 degrees in the individual simulations, respectively. This simulated treeline shift is in agreement with pollen-based reconstructions from northern Siberia. The desert fraction in the transition zone is reduced by 21% during the mid-Holocene compared to pre-industrial due to enhanced precipitation. The desert-steppe margin is shifted westward by 5 degrees (1-9 degrees in the individual simulations). The forest biomes are expanded north-westward by 2 degrees, ranging from 0 to 4 degrees in the single simulations. These results corroborate pollen-based reconstructions indicating an extended forest area in north-central China during the mid-Holocene. According to the model, the forest-to-non-forest and steppe-to-desert changes in the climate transition zones are spatially not uniform and not linear since the mid-Holocene.
Temporal and spatial stability of the vegetation climate relationship is a basic ecological assumption for pollen-based quantitative inferences of past climate change and for predicting future vegetation. We explore this assumption for the Holocene in eastern continental Asia (China, Mongolia). Boosted regression trees (BRT) between fossil pollen taxa percentages (Abies, Artemisia, Betula, Chenopodiaceae, Cyperaceae, Ephedra, Picea, Pinus, Poaceae and Quercus) and climate model outputs of mean annual precipitation (P-ann) and mean temperature of the warmest month (Mt(wa)) for 9 and 6 ka (ka = thousand years before present) were set up and results compared to those obtained from relating modern pollen to modern climate. Overall, our results reveal only slight temporal differences in the pollen climate relationships. Our analyses suggest that the importance of P-ann compared with Mt(wa) for taxa distribution is higher today than it was at 6 ka and 9 ka. In particular, the relevance of P-ann for Picea and Pinus increases and has become the main determinant. This change in the climate tree pollen relationship parallels a widespread tree pollen decrease in north-central China and the eastern Tibetan Plateau. We assume that this is at least partly related to vegetation climate disequilibrium originating from human impact. Increased atmospheric CO2 concentration may have permitted the expansion of moisture-loving herb taxa (Cyperaceae and Poaceae) during the late Holocene into arid/semi-arid areas. We furthermore find that the pollen climate relationship between north-central China and the eastern Tibetan Plateau is generally similar, but that regional differences are larger than temporal differences. In summary, vegetation climate relationships in China are generally stable in space and time, and pollen-based climate reconstructions can be applied to the Holocene. Regional differences imply the calibration-set should be restricted spatially.
We have analyzed the recently developed pan-European strong motion database, RESORCE-2012: spectral parameters, such as stress drop (stress parameter, Delta sigma), anelastic attenuation (Q), near surface attenuation (kappa(0)) and site amplification have been estimated from observed strong motion recordings. The selected dataset exhibits a bilinear distance-dependent Q model with average kappa(0) value 0.0308 s. Strong regional variations in inelastic attenuation were also observed: frequency-independent Q(0) of 1462 and 601 were estimated for Turkish and Italian data respectively. Due to the strong coupling between Q and kappa(0), the regional variations in Q have strong impact on the estimation of near surface attenuation kappa(0). kappa(0) was estimated as 0.0457 and 0.0261 s for Turkey and Italy respectively. Furthermore, a detailed analysis of the variability in estimated kappa(0) revealed significant within-station variability. The linear site amplification factors were constrained from residual analysis at each station and site-class type. Using the regional Q(0) model and a site-class specific kappa(0), seismic moments (M-0) and source corner frequencies f (c) were estimated from the site corrected empirical Fourier spectra. Delta sigma did not exhibit magnitude dependence. The median Delta sigma value was obtained as 5.75 and 5.65 MPa from inverted and database magnitudes respectively. A comparison of response spectra from the stochastic model (derived herein) with that from (regional) ground motion prediction equations (GMPEs) suggests that the presented seismological parameters can be used to represent the corresponding seismological attributes of the regional GMPEs in a host-to-target adjustment framework. The analysis presented herein can be considered as an update of that undertaken for the previous Euro-Mediterranean strong motion database presented by Edwards and Fah (Geophys J Int 194(2):1190-1202, 2013a).
Receiver functions (RF) have been used for several decades to study structures beneath seismic stations. Although most available stations are deployed on shore, the number of ocean bottom station (OBS) experiments has increased in recent years. Almost all OBSs have to deal with higher noise levels and a limited deployment time (approximate to 1year), resulting in a small number of usable records of teleseismic earthquakes. Here we use OBSs deployed as midaperture array in the deep ocean (4.5-5.5km water depth) of the eastern mid-Atlantic. We use evaluation criteria for OBS data and beamforming to enhance the quality of the RFs. Although some stations show reverberations caused by sedimentary cover, we are able to identify the Moho signal, indicating a normal thickness (5-8km) of oceanic crust. Observations at single stations with thin sediments (300-400m) indicate that a probable sharp lithosphere-asthenosphere boundary (LAB) might exist at a depth of approximate to 70-80km which is in line with LAB depth estimates for similar lithospheric ages in the Pacific. The mantle discontinuities at approximate to 410km and approximate to 660km are clearly identifiable. Their delay times are in agreement with PREM. Overall the usage of beam-formed earthquake recordings for OBS RF analysis is an excellent way to increase the signal quality and the number of usable events.
Two nuclear explosions were carried out by the Democratic People’s Republic of North Korea in January and September 2016. Epicenters were located close to those of the 2006, 2009, and 2013 previous explosions. We perform a seismological analysis of the 2016 events combining the analysis of full waveforms at regional distances and seismic array beams at teleseismic distances. We estimate the most relevant source parameters, such as source depth, moment release, and full moment tensor (MT). The best MT solution can be decomposed into an isotropic source, directly related with the explosion and an additional deviatoric term, likely due to near‐source interactions with topographic and/or underground facilities features. We additionally perform an accurate resolution test to assess source parameters uncertainties and trade‐offs. This analysis sheds light on source parameters inconsistencies among studies on previous shallow explosive sources. The resolution of the true MT is hindered by strong source parameters trade‐offs, so that a broad range of well‐fitting MT solutions can be found, spanning from a dominant positive isotropic term to a dominant negative vertical compensated linear vector dipole. The true mechanism can be discriminated by additionally modeling first‐motion polarities at seismic arrays at teleseismic distances. A comparative assessment of the 2016 explosion with earlier nuclear tests documents similar vertical waveforms but a significant increase of amplitude for the 2016 explosions, which proves that the 9 September 2016 was the largest nuclear explosion ever performed in North Korea with a magnitude Mw 4.9 and a shallow depth of less than 2 km, although there are no proofs of a fusion explosion. Modeling transversal component waveforms suggests variable size and orientation of the double‐couple components of the 2009, 2013, and 2016 sources.
The increasing numbers of recordings at individual sites allows quantification of empirical linear site-response adjustment factors (delta S2S(s)) from the ground motion prediction equation (GMPE) residuals. The delta S2S(s) are then used to linearly scale the ergodic GMPE predictions to obtain site-specific ground motion predictions in a partially non-ergodic Probabilistic Seismic Hazard Assessment (PSHA). To address key statistical and conceptual issues in the current practice, we introduce a novel empirical region-and site-specific PSHA methodology wherein, (1) site-to-site variability (phi(S2S)) is first estimated as a random-variance in a mixed-effects GMPE regression, (2) delta S2S(s) at new sites with strong motion are estimated using the a priori phi(S2S), and (3) the GMPE site-specific single-site aleatory variability sigma(ss,s) is replaced with a generic site-corrected aleatory variability sigma(0). Comparison of region- and site-specific hazard curves from our method against the traditional ergodic estimates at 225 sites in Europe and Middle East shows an approximate 50% difference in predicted ground motions over a range of hazard levels-a strong motivation to increase seismological monitoring of critical facilities and enrich regional ground motion data sets.
Site-Corrected Magnitude- and Region-Dependent Correlations of Horizontal Peak Spectral Amplitudes
(2017)
Empirical correlations of horizontal peak spectral amplitudes (PSA) are modeled using the total-residuals obtained in a ground motion prediction equation (GMPE) regression. Recent GMPEs moved toward partially non-ergodic region-and site-specific predictions, while the residual correlation models remained largely ergodic. Using mixed-effects regression, we decompose the total-residuals of a pan-European GMPE into between-event, between-site, and event-and-site corrected residuals to investigate the ergodicity in empirical PSA correlations. We first observed that the between-event correlations are magnitude-dependent, partially due to the differences in source spectra, and influence of stress-drop parameter on small and large events. Next, removing the between-site residuals from within-event residuals yields the event-and-site corrected residuals which are found to be region-dependent, possibly due to the regional differences in distance-decay of short period PSAs. Using our site-corrected magnitude- and region-dependent correlations, and the between-site residuals as empirical site-specific ground motion adjustments, we compute partially non-ergodic conditional mean spectra at four well-recorded sites in Europe and Middle Eastern regions.
Most of the deformation associated with the seismic cycle in subduction zones occurs offshore and has been therefore difficult to quantify with direct observations at millennial timescales. Here we study millennial deformation associated with an active splay-fault system in the Arauco Bay area off south central Chile. We describe hitherto unrecognized drowned shorelines using high-resolution multibeam bathymetry, geomorphic, sedimentologic, and paleontologic observations and quantify uplift rates using a Landscape Evolution Model. Along a margin-normal profile, uplift rates are 1.3m/ka near the edge of the continental shelf, 1.5m/ka at the emerged Santa Maria Island, -0.1m/ka at the center of the Arauco Bay, and 0.3m/ka in the mainland. The bathymetry images a complex pattern of folds and faults representing the surface expression of the crustal-scale Santa Maria splay-fault system. We modeled surface deformation using two different structural scenarios: deep-reaching normal faults and deep-reaching reverse faults with shallow extensional structures. Our preferred model comprises a blind reverse fault extending from 3km depth down to the plate interface at 16km that slips at a rate between 3.0 and 3.7m/ka. If all the splay-fault slip occurs during every great megathrust earthquake, with a recurrence of similar to 150-200years, the fault would slip similar to 0.5m per event, equivalent to a magnitude similar to 6.4 earthquake. However, if the splay-fault slips only with a megathrust earthquake every similar to 1000years, the fault would slip similar to 3.7m per event, equivalent to a magnitude similar to 7.5 earthquake.
Aim: Across the planet, grass-dominated biomes are experiencing shrub encroachment driven by atmospheric CO2 enrichment and land-use change. By altering resource structure and availability, shrub encroachment may have important impacts on vertebrate communities. We sought to determine the magnitude and variability of these effects across climatic gradients, continents, and taxa, and to learn whether shrub thinning restores the structure of vertebrate communities. Location: Worldwide. Time period: Contemporary. Major taxa studied: Terrestrial vertebrates. Methods: We estimated relationships between percentage shrub cover and the structure of terrestrial vertebrate communities (species richness, Shannon diversity and community abundance) in experimentally thinned and unmanipulated shrub-encroached grass-dominated biomes using systematic review and meta-analyses of 43 studies published from 1978 to 2016. We modelled the effects of continent, biome, mean annual precipitation, net primary productivity and the normalized difference vegetation index (NDVI) on the relationship between shrub cover and vertebrate community structure. Results: Species richness, Shannon diversity and total abundance had no consistent relationship with shrub encroachment and experimental thinning did not reverse encroachment effects on vertebrate communities. However, some effects of shrub encroachment on vertebrate communities differed with net primary productivity, amongst vertebrate groups, and across continents. Encroachment had negative effects on vertebrate diversity at low net primary productivity. Mammalian and herpetofaunal diversity decreased with shrub encroachment. Shrub encroachment also had negative effects on species richness and total abundance in Africa but positive effects in North America. Main conclusions: Biodiversity conservation and mitigation efforts responding to shrub encroachment should focus on low-productivity locations, on mammals and herpetofauna, and in Africa. However, targeted research in neglected regions such as central Asia and India will be needed to fill important gaps in our knowledge of shrub encroachment effects on vertebrates. Additionally, our findings provide an impetus for determining the mechanisms associated with changes in vertebrate diversity and abundance in shrub-encroached grass-dominated biomes.
A sediment core from Khatanga-12 Lake (Taimyr Peninsula, Krasnoyarsk krai) has been studied. The 131.5-cm-long core covers ca. 7100 years of sedimentation. Chironomid analysis, a qualitative reconstruction of the paleoenvironment in the region, and a quantitative reconstruction of variations of the mean July air temperature and in the water depth of the lake have been performed using Northern Russia chironomid-inferred mean July temperature models (Nazarova et al., 2008, 2011, 2015). Khatanga-12 Lake was formed during the Middle Holocene warming as a result of thermokarst processes. The development of the lake ecosystem at different stages of its development was influenced by climatic and cryolithogenic factors. The Middle Holocene warming, which occurred around 7100-6250 cal. years BP, activated thermokarst processes and resulted in the formation of the lake basin. Later, between 6250 and 4500 cal. years BP, a period of cooling took place, as is proved by chironomid analysis. The bottom sediments of the lake during this period were formed by erosion processes on the lake shores. The reconstructed conditions were close to the modern after 2500 cal. years BP.
The Himalaya has a major influence on global and regional climate, in particular on the Asian monsoon system. The foreland basin of the Himalaya contains a record of tectonics and paleoclimate since the Miocene. Previous work on the evolution of vegetation and climate has focused on the central and western Himalaya, where a shift from C3 to C4 vegetation has been observed at similar to 7 Ma and linked to increased seasonality, but the climatic evolution of the eastern part of the orogen is less well understood. In order to track vegetation as a marker of monsoon intensity and seasonality, we analyzed delta C-13 and 8180 values of soil carbonate and associated delta C-13 values of bulk organic carbon from previously dated sedimentary sections exposing the syn-orogenic detrital Dharamsala and Siwalik Groups in the west, and, for the first time, the Siwalik Group in the east of the Himalayan foreland basin. Sedimentary records span from 20 to 1 Myr in the west (Joginder Nagar, Jawalamukhi, and Haripur Kolar sections) and from 13 to 1 Myr in the east (Kameng section), respectively. The presence of soil carbonate in the west and its absence in the east is a first indication of long-term lateral climatic variation, as soil carbonate requires seasonally arid conditions to develop. delta C-13 values in soil carbonate show a shift from around -10 parts per thousand to -2 parts per thousand at similar to 7 Ma in the west, which is confirmed by delta C-13 analyses on bulk organic carbon that show a shift from around -23 parts per thousand to -19 parts per thousand at the same time. Such a shift in isotopic values is likely to be associated with a change from C3 to C4 vegetation. In contrast, delta C-13 values of bulk organic carbon remain at 23 parts per thousand o in the east. Thus, our data show that the current east -west variation in climate was established at similar to 7 Ma. We propose that the regional change towards a more seasonal climate in the west is linked to a decrease of the influence of the Westerlies, delivering less winter precipitation to the western Himalaya, while the east remained annually humid due to its proximity to the monsoonal moisture source. (C) 2017 Elsevier B.V. All rights reserved.
Direct current systems employing a kinematic surveying strategy allow to analyze the electrical resistivity of the subsurface for large areas (i.e., several hectares). Typical applications are found in precision agriculture, archaeological prospecting and soil sciences. With the typical survey setting, the collected data sets are often characterized by a rather high level of noise and a rather coarse lateral sampling compared to data acquired with fixed electrodes. We therefore present an efficient one-dimensional inversion approach in which we put special attention on modeling the effects of noise. We apply this method to data recorded with a five-offset equatorial dipole-dipole system employing rolling electrodes. By performing several synthetic tests with realistic noise levels, we found that the considered five-configuration soundings allow for a reliable imaging of two-layer cases in the uppermost two meters of the subsurface, where the subsurface can be assumed to follow a horizontally layered geometry within 3 m around the system. By analyzing the corresponding sensitivity functions, we also show that the equatorial dipole-dipole array is relatively well suited for a 1D inversion approach compared to standard in-line electrode arrays. To illustrate this aspect, we show that our method can provide results similar to those obtained with a 2D Wenner imaging procedure for data recorded across a well-constrained 2D target. We finally apply our method to a large five-offset data set acquired in an agricultural study. The final pseudo-3D model of electrical resistivity is in accordance with borehole data available for the surveyed area. Our results demonstrate the applicability and the versatility of the presented inversion approach for large-scale data sets as they are typically collected with such rolling electrode systems. (C) 2017 Elsevier B.V. All rights reserved.
Electromagnetic induction (EMI) sensors using sufficiently low-frequency harmonic sources and sufficiently small loop separations operate in the low-induction-number (LIN) domain for a relatively wide range of background conductivity. These systems are used in diverse near-surface investigations including applications from soil sciences, hydrology, and archaeology. The special case of portable multiconfiguration EMI sensors operating at frequencies <= 20 kHz offers the possibility of using a fast linear deconvolution method to interpret multichannel data sets in three dimensions. Here, we have developed a fast 3D inversion/deconvolution method regularized with 3D smoothness constraints and formulated in the hybrid spectral-spatial domain. Compared with other linear approaches, the spectral-spatial domain formulation significantly reduces the computational cost of the processing and opens the door for real-time 3D interpretation of large data sets consisting of more than 100,000 data points. First, we test our proposed algorithm on synthetic data sets computed with the full Maxwell theory. Then, we apply our method to a real four-configuration EMI data set acquired to map the thickness of peat layers embedded in a sandy environment. For the synthetic and the field example, we compared our result with the result obtained using a standard point-by-point 1D nonlinear inversion approach. This comparison demonstrates that the proposed methodology provides superior lateral resolution compared with the 1D nonlinear inversion, at the same time significantly reducing the computational cost of the processing.
The phrase form and function was established in architecture and biology and refers to the idea that form and functionality are closely correlated, influence each other, and co-evolve. We suggest transferring this idea to hydrological systems to separate and analyze their two main characteristics: their form, which is equivalent to the spatial structure and static properties, and their function, equivalent to internal responses and hydrological behavior. While this approach is not particularly new to hydrological field research, we want to employ this concept to explicitly pursue the question of what information is most advantageous to understand a hydrological system. We applied this concept to subsurface flow within a hillslope, with a methodological focus on function: we conducted observations during a natural storm event and followed this with a hillslope-scale irrigation experiment. The results are used to infer hydrological processes of the monitored system. Based on these findings, the explanatory power and conclusiveness of the data are discussed. The measurements included basic hydrological monitoring methods, like piezometers, soil moisture, and discharge measurements. These were accompanied by isotope sampling and a novel application of 2-D time-lapse GPR (ground-penetrating radar). The main finding regarding the processes in the hillslope was that preferential flow paths were established quickly, despite unsaturated conditions. These flow paths also caused a detectable signal in the catchment response following a natural rainfall event, showing that these processes are relevant also at the catchment scale. Thus, we conclude that response observations (dynamics and patterns, i.e., indicators of function) were well suited to describing processes at the observational scale. Especially the use of 2-D time-lapse GPR measurements, providing detailed subsurface response patterns, as well as the combination of stream-centered and hillslope-centered approaches, allowed us to link processes and put them in a larger context. Transfer to other scales beyond observational scale and generalizations, however, rely on the knowledge of structures (form) and remain speculative. The complementary approach with a methodological focus on form (i.e., structure exploration) is presented and discussed in the companion paper by Jackisch et al. (2017).
The study deals with the identification and characterization of rapid subsurface flow structures through pedo- and geo-physical measurements and irrigation experiments at the point, plot and hillslope scale. Our investigation of flow-relevant structures and hydrological responses refers to the general interplay of form and function, respectively. To obtain a holistic picture of the subsurface, a large set of different laboratory, exploratory and experimental methods was used at the different scales. For exploration these methods included drilled soil core profiles, in situ measurements of infiltration capacity and saturated hydraulic conductivity, and laboratory analyses of soil water retention and saturated hydraulic conductivity. The irrigation experiments at the plot scale were monitored through a combination of dye tracer, salt tracer, soil moisture dynamics, and 3-D time-lapse ground penetrating radar (GPR) methods. At the hillslope scale the subsurface was explored by a 3-D GPR survey. A natural storm event and an irrigation experiment were monitored by a dense network of soil moisture observations and a cascade of 2-D time-lapse GPR "trenches". We show that the shift between activated and non-activated state of the flow paths is needed to distinguish structures from overall heterogeneity. Pedo-physical analyses of point-scale samples are the basis for sub-scale structure inference. At the plot and hillslope scale 3-D and 2-D time-lapse GPR applications are successfully employed as non-invasive means to image subsurface response patterns and to identify flow-relevant paths. Tracer recovery and soil water responses from irrigation experiments deliver a consistent estimate of response velocities. The combined observation of form and function under active conditions provides the means to localize and characterize the structures (this study) and the hydrological processes (companion study Angermann et al., 2017, this issue).
The African Humid Period (AHP) between similar to 15 and 5.5 cal. kyr BP caused major environmental change in East Africa, including filling of the Suguta Valley in the northern Kenya Rift with an extensive (similar to 2150 km(2)), deep (similar to 300 m) lake. Interfingering fluvio-lacustrine deposits of the Baragoi paleo-delta provide insights into the lake-level history and how erosion rates changed during this time, as revealed by delta-volume estimates and the concentration of cosmogenic Be-10 in fluvial sand. Erosion rates derived from delta-volume estimates range from 0.019 to 0.03 mm yr(-1). Be-10-derived paleo-erosion rates at similar to 11.8 cal. kyr BP ranged from 0.035 to 0.086 mm yr(-1), and were 2.7 to 6.6 times faster than at present. In contrast, at similar to 8.7 cal. kyr BP, erosion rates were only 1.8 times faster than at present. Because Be-10-derived erosion rates integrate over several millennia; we modeled the erosion-rate history that best explains the 10Be data using established non-linear equations that describe in situ cosmogenic isotope production and decay. Two models with different temporal constraints (15-6.7 and 12-6.7 kyr) suggest erosion rates that were 25 to 300 times higher than the initial erosion rate (pre-delta formation). That pulse of high erosion rates was short (similar to 4 kyr or less) and must have been followed by a rapid decrease in rates while climate remained humid to reach the modern Be-10-based erosion rate of,similar to 0.013 mm yr(-1). Our simulations also flag the two highest Be-10-derived erosion rates at 11.8 kyr BP related to nonuniform catchment erosion. These changes in erosion rates and processes during the AHP may reflect a strong increase in precipitation, runoff, and erosivity at the arid-to-humid transition either at 15 or similar to 12 cal. kyr BP, before the landscape stabilized again, possibly due to increased soil production and denser vegetation.
Near-surface geophysical techniques are extensively used in a variety of engineering, environmental, geologic, and hydrologic applications. While many of these applications ask for detailed, quantitative models of selected material properties, geophysical data are increasingly used to estimate such target properties. Typically, this estimation procedure relies on a two-step workflow including (1) the inversion of geophysical data and (2) the petrophysical translation of the inverted parameter models into the target properties. Standard deterministic implementations of such a quantitative interpretation result in a single best-estimate model, often without considering and propagating the uncertainties related to the two steps. We address this problem by using a rather novel, particle-swarm-based global joint strategy for data inversion and by implementing Monte Carlo procedures for petrophysical property estimation. We apply our proposed workflow to crosshole ground-penetrating radar, P-, and S-wave data sets collected at a well-constrained test site for a detailed geotechnical characterization of unconsolidated sands. For joint traveltime inversion, the chosen global approach results in ensembles of acceptable velocity models, which are analyzed to appraise inversion-related uncertainties. Subsequently, the entire ensembles of inverted velocity models are considered to estimate selected petrophysical properties including porosity, bulk density, and elastic moduli via well-established petrophysical relations implemented in a Monte Carlo framework. Our results illustrate the potential benefit of such an advanced interpretation strategy; i.e., the proposed workflow allows to study how uncertainties propagate into the finally estimated property models, while concurrently we are able to study the impact of uncertainties in the used petrophysical relations (e.g., the influence of uncertain, user-specified parameters). We conclude that such statistical approaches for the quantitative interpretation of geophysical data can be easily extended and adapted to other applications and geophysical methods and might be an important step toward increasing the popularity and acceptance of geophysical tools in engineering practice.
In the last few years the method of cosmic-ray neutron sensing (CRNS) has gained popularity among hydrologists, physicists, and land-surface modelers. The sensor provides continuous soil moisture data, averaged over several hectares and tens of decimeters in depth. However, the signal still may contain unidentified features of hydrological processes, and many calibration datasets are often required in order to find reliable relations between neutron intensity and water dynamics. Recent insights into environmental neutrons accurately described the spatial sensitivity of the sensor and thus allowed one to quantify the contribution of individual sample locations to the CRNS signal. Consequently, data points of calibration and validation datasets are suggested to be averaged using a more physically based weighting approach. In this work, a revised sensitivity function is used to calculate weighted averages of point data. The function is different from the simple exponential convention by the extraordinary sensitivity to the first few meters around the probe, and by dependencies on air pressure, air humidity, soil moisture, and vegetation. The approach is extensively tested at six distinct monitoring sites: two sites with multiple calibration datasets and four sites with continuous time series datasets. In all cases, the revised averaging method improved the performance of the CRNS products. The revised approach further helped to reveal hidden hydrological processes which otherwise remained unexplained in the data or were lost in the process of overcalibration. The presented weighting approach increases the overall accuracy of CRNS products and will have an impact on all their applications in agriculture, hydrology, and modeling.
The formation of the Pamir is a key component of the India-Asia collision with major implications for lithospheric processes, plateau formation, land-sea configurations and associated climate changes. Although the formation of the Pamir is traditionally linked to Cenozoic processes associated with the India-Asia collision, the contribution of the Mesozoic tectonic evolution remains poorly understood. The Pamir was formed by the suturing of Gondwanan terranes to the south margin of Eurasia, however, the timing and tectonic mechanisms associated with this Mesozoic accretion remain poorly constrained. These processes are recorded by several igneous belts within these terranes, which are not well studied. Within the Southern Pamir, the Albian-Turonian volcanic rocks and comagmatic plutons of the Kyzylrabat Igneous Complex (KIC) provide an important and still unconstrained record of the Pamir evolution. Here we provide the age, origin and the geodynamic setting of the KIC volcanics by studying their petrology, zircon U-Pb geochronology, geochemistry and isotope composition.17 samples from the KIC volcanics yield U-Pb ages spanning from 92 to 110 Ma. The volcanics are intermediate to acidic in composition (SiO2 = 56-69 wt%) and exhibit high-K calc-alkaline and shoshonitic affinity (K2O/Na2O = 12.2 wt%). They show enrichment in LILE and LREE and depletion in HFSE and HREE with negative Ta, Ti and Nb anomalies, suggesting an arc-related tectonic setting for their formation. Low sNd(t) values (from 9.1 to 4.7), relatively high Sr-87/Sr-86(i) ratios (0.7069-0.7096) and broad range of zircon stif values (from 22.6 to 1.5) suggest a mixture of different magma sources. These features suggest that volcanics were derived by crustal under- or intraplating of an enriched subduction-related mantle shoshonitic magmas, by heating and partial melting of the lower crust, and by mixing of both magma components. Our results further imply that the KIC volcanics represent a shoshonitic suite typical of an evolution from active continental arc to post-collisional setting with a steepening of the Benioff zone and thickening of the crust toward the back-arc. This setting is best explained by the subduction- collision transition along the Shyok suture due to accretion of the Kohistan island arc to the Karakoram. This suggests that a significant part of the crustal shortening and thickening accommodated in the Pamir occurred in the Mesozoic before the India-Asia collision with implications for regional tectonic models. This further suggests the Pamir was already a major topographic feature with potentially important paleoclimate forcing such as the monsoonal circulation. (C) 2017 International Association for Gondwana Research. Published by Elsevier B.V. All rights reserved.
The Acheulian site of Gesher Benot Ya‘aqov (GBY) in the Upper Jordan Valley revealed important data on environment and material culture, as well as evidence for hominin behavioural and cognitive patterns documented at the margins of the Hula Palaeo-lake. A 50 m long core (GBY#2) drilled at the archaeological site has provided a long Pleistocene geological, environmental and climatological record, which expands the existing knowledge of hominin-habitat relationships. Bracketed by two basalt flows dated by 40Ar/39Ar and based on the identification of the Matuyama-Brunhes Boundary (MBB) and correlation with the GBY excavation site, the sedimentary sequence provides the climatic history around the MBB. Multi-proxy data including pollen and non-pollen palynomorphs, macro-botanical remains, molluscs and ostracods provide evidence for lake and lake-margin environments during Marine Isotope Stages (MIS) 20 and 19. Semi-moist conditions were followed by a pronounced dry phase during MIS 20, and warm and moist conditions with Quercus-Pistacia woodlands prevailed during MIS 19. In contrast to the reconstructed climate change from relatively dry to moister conditions, the depositional environment developed from an open-water lake during MIS 20 to a lake margin environment in MIS 19. Generally shallower conditions at the core site in MIS 19 resulted from the progradation of the lake shore due to the filling of the basin. Micro-charcoal analysis suggests a likelihood of human-induced fire in some parts of the core, which can be correlated with artefact-containing layers of the GBY excavation site. The Hula Palaeo-lake region provided an ideal niche for hominins and other vertebrates during global glacial-interglacial climate fluctuations at the end of the Early Pleistocene.
One of the major challenges in understanding the evolution of our own species is identifying the role climate change has played in the evolution of hominin species. To clarify the influence of climate, we need long and continuous high-resolution paleoclimate records, preferably obtained from hominin-bearing sediments, that are well-dated by tephro- and magnetostratigraphy and other methods. This is hindered, however, by the fact that fossil-bearing outcrop sediments are often discontinuous, and subject to weathering, which may lead to oxidation and remagnetization. To obtain fresh, unweathered sediments, the Hominin Sites and Paleolakes Drilling Project (HSPDP) collected a ∼216-meter core (WTK13) in 2013 from Early Pleistocene Paleolake Lorenyang deposits in the western Turkana Basin (Kenya). Here, we present the magnetostratigraphy of the WTK13 core, providing a first age model for upcoming HSPDP paleoclimate and paleoenvrionmental studies on the core sediments. Rock magnetic analyses reveal the presence of iron sulfides carrying the remanent magnetizations. To recover polarity orientation from the near-equatorial WTK13 core drilled at 5°N, we developed and successfully applied two independent drill-core reorientation methods taking advantage of (1) the sedimentary fabric as expressed in the Anisotropy of Magnetic Susceptibility (AMS) and (2) the occurrence of a viscous component oriented in the present day field. The reoriented directions reveal a normal to reversed polarity reversal identified as the top of the Olduvai Subchron. From this excellent record, we find no evidence for the ‘Vrica Subchron’ previously reported in the area. We suggest that outcrop-based interpretations supporting the presence of the Vrica Subchron have been affected by the oxidation of iron sulfides initially present in the sediments -as evident in the core record- and by subsequent remagnetization. We discuss the implications of the observed geomagnetic record for human evolution studies.
Sedimentary basins in the interior of orogenic plateaus can provide unique insights into the early history of plateau evolution and related geodynamic processes. The northern sectors of the Iranian Plateau of the Arabia-Eurasia collision zone offer the unique possibility to study middle-late Miocene terrestrial clastic and volcaniclastic sediments that allow assessing the nascent stages of collisional plateau formation. In particular, these sedimentary archives allow investigating several debated and poorly understood issues associated with the long-term evolution of the Iranian Plateau, including the regional spatio-temporal characteristics of sedimentation and deformation and the mechanisms of plateau growth. We document that middle-late Miocene crustal shortening and thickening processes led to the growth of a basement-cored range (Takab Range Complex) in the interior of the plateau. This triggered the development of a foreland-basin (Great Pari Basin) to the east between 16.5 and 10.7Ma. By 10.7Ma, a fast progradation of conglomerates over the foreland strata occurred, most likely during a decrease in flexural subsidence triggered by rock uplift along an intraforeland basement-cored range (Mahneshan Range Complex). This was in turn followed by the final incorporation of the foreland deposits into the orogenic system and ensuing compartmentalization of the formerly contiguous foreland into several intermontane basins. Overall, our data suggest that shortening and thickening processes led to the outward and vertical growth of the northern sectors of the Iranian Plateau starting from the middle Miocene. This implies that mantle-flow processes may have had a limited contribution toward building the Iranian Plateau in NW Iran.
Located along the Izmir-Ankara-Erzincan Suture (IAES), the Maastrichtian - Paleogene Orhaniye Basin has yielded a highly enigmatic-yet poorly dated- Paleogene mammal fauna, the endemic character of which has suggested high faunal provincialism associated with paleogeographic isolation of the Anatolian landmass during the early Cenozoic. Despite its biogeographic significance, the tectono-stratigraphic history of the Orhaniye Basin has been poorly documented; Here, we combine sedimentary, magnetostratigraphic, and geochronological data to infer the chronology and depositional history of the Orhaniye Basin. We then assess how our new data and interpretations for the Orhaniye Basin impact (1) the timing and mechanisms of seaway closure along the IAES and (2) the biogeographic evolution of Anatolia. Our results show that the Orhaniye Basin initially developed as a forearc basin during the Maastrichtian, before shifting to a retroarc foreland basin setting sometime between the early Paleocene and 44 Ma. This chronology supports a two-step scenario for the assemblage of the central Anatolian landmass, with incipient collision during the Paleocene - Early Eocene and final seaway retreat along the IAES during the earliest Late Eocene after the last marine incursion into the foreland basin. Our dating for the Orhaniye mammal fauna (44-43 Ma) indicates the persistence of faunal endemism in northern Anatolia until at least the late Lutetian despite the advanced stage of IAES closure. The tectonic evolution of dispersal corridors linking northern Anatolia with adjacent parts of Eurasia was not directly associated with IAES closure and consecutive uplifts, but rather with the build-up of continental bridges on the margins of Anatolia, in the Alpine and Tibetan-Himalayan orogens.