Refine
Year of publication
- 2016 (291) (remove)
Document Type
- Article (221)
- Doctoral Thesis (37)
- Other (14)
- Postprint (12)
- Review (5)
- Habilitation Thesis (1)
- Report (1)
Is part of the Bibliography
- yes (291) (remove)
Keywords
- climate change (5)
- erosion (5)
- water balance (5)
- Pollen (4)
- August 2002 flood (3)
- Body waves (3)
- Central Europe (3)
- Floods Directive (3)
- Himalaya (3)
- Himalayas (3)
Institute
- Institut für Geowissenschaften (291) (remove)
Intermontane valley fills
(2016)
Sedimentary valley fills are a widespread characteristic of mountain belts around the world. They transiently store material over time spans ranging from thousands to millions of years and therefore play an important role in modulating the sediment flux from the orogen to the foreland and to oceanic depocenters. In most cases, their formation can be attributed to specific fluvial conditions, which are closely related to climatic and tectonic processes. Hence, valley-fill deposits constitute valuable archives that offer fundamental insight into landscape evolution, and their study may help to assess the impact of future climate change on sediment dynamics.
In this thesis I analyzed intermontane valley-fill deposits to constrain different aspects of the climatic and tectonic history of mountain belts over multiple timescales. First, I developed a method to estimate the thickness distribution of valley fills using artificial neural networks (ANNs). Based on the assumption of geometrical similarity between exposed and buried parts of the landscape, this novel and highly automated technique allows reconstructing fill thickness and bedrock topography on the scale of catchments to entire mountain belts.
Second, I used the new method for estimating the spatial distribution of post-glacial sediments that are stored in the entire European Alps. A comparison with data from exploratory drillings and from geophysical surveys revealed that the model reproduces the measurements with a root mean squared error (RMSE) of 70m and a coefficient of determination (R2) of 0.81. I used the derived sediment thickness estimates in combination with a model of the Last Glacial Maximum (LGM) icecap to infer the lithospheric response to deglaciation, erosion and deposition, and deduce their relative contribution to the present-day rock-uplift rate. For a range of different lithospheric and upper mantle-material properties, the results suggest that the long-wavelength uplift signal can be explained by glacial isostatic adjustment with a small erosional contribution and a substantial but localized tectonic component exceeding 50% in parts of the Eastern Alps and in the Swiss Rhône Valley. Furthermore, this study reveals the particular importance of deconvolving the potential components of rock uplift when interpreting recent movements along active orogens and how this can be used to constrain physical properties of the Earth’s interior.
In a third study, I used the ANN approach to estimate the sediment thickness of alluviated reaches of the Yarlung Tsangpo River, upstream of the rapidly uplifting Namche Barwa massif. This allowed my colleagues and me to reconstruct the ancient river profile of the Yarlung Tsangpo, and to show that in the past, the river had already been deeply incised into the eastern margin of the Tibetan Plateau. Dating of basal sediments from drill cores that reached the paleo-river bed to 2–2.5 Ma are consistent with mineral cooling ages from the Namche Barwa massif, which indicate initiation of rapid uplift at ~4 Ma. Hence, formation of the Tsangpo gorge and aggradation of the voluminous valley fill was most probably a consequence of rapid uplift of the Namche Barwa massif and thus tectonic activity.
The fourth and last study focuses on the interaction of fluvial and glacial processes at the southeastern edge of the Karakoram. Paleo-ice-extent indicators and remnants of a more than 400-m-thick fluvio-lacustrine valley fill point to blockage of the Shyok River, a main tributary of the upper Indus, by the Siachen Glacier, which is the largest glacier in the Karakoram Range. Field observations and 10Be exposure dating attest to a period of recurring lake formation and outburst flooding during the penultimate glaciation prior to ~110 ka. The interaction of Rivers and Glaciers all along the Karakorum is considered a key factor in landscape evolution and presumably promoted headward erosion of the Indus-Shyok drainage system into the western margin of the Tibetan Plateau.
The results of this thesis highlight the strong influence of glaciation and tectonics on valley-fill formation and how this has affected the evolution of different mountain belts. In the Alps valley-fill deposition influenced the magnitude and pattern of rock uplift since ice retreat approximately 17,000 years ago. Conversely, the analyzed valley fills in the Himalaya are much older and reflect environmental conditions that prevailed at ~110 ka and ~2.5 Ma, respectively. Thus, the newly developed method has proven useful for inferring the role of sedimentary valley-fill deposits in landscape evolution on timescales ranging from 1,000 to 10,000,000 years.
Variations in the distribution of mass within an orogen may lead to transient sediment storage, which in turn might affect the state of stress and the level of fault activity. Distinguishing between different forcing mechanisms causing variations of sediment flux and tectonic activity, is therefore one of the most challenging tasks in understanding the spatiotemporal evolution of active mountain belts.
The Himalayan mountain belt is one of the most significant Cenozoic collisional mountain belt, formed due to collision between northward-bound Indian Plate and the Eurasian Plate during the last 55-50 Ma. Ongoing convergence of these two tectonic plates is accommodated by faulting and folding within the Himalayan arc-shaped orogen and the continued lateral and vertical growth of the Tibetan Plateau and mountain belts adjacent to the plateau as well as regions farther north. Growth of the Himalayan orogen is manifested by the development of successive south-vergent thrust systems. These thrust systems divide the orogen into different morphotectonic domains. From north to south these thrusts are the Main Central Thrust (MCT), the Main Boundary Thrust (MBT) and the Main Frontal Thrust (MFT). The growing topography interacts with moisture-bearing monsoonal winds, which results in pronounced gradients in rainfall, weathering, erosion and sediment transport toward the foreland and beyond. However, a fraction of this sediment is trapped and transiently stored within the intermontane valleys or ‘dun’s within the lower-elevation foothills of the range. Improved understanding of the spatiotemporal evolution of these sediment archives could provide a unique opportunity to decipher the triggers of variations in sediment production, delivery and storage in an actively deforming mountain belt and support efforts to test linkages between sediment volumes in intermontane basins and changes in the shallow crustal stress field. As sediment redistribution in mountain belts on timescales of 102-104 years can effect cultural characteristics and infrastructure in the intermontane valleys and may even impact the seismotectonics of a mountain belt, there is a heightened interest in understanding sediment-routing processes and causal relationships between tectonism, climate and topography. It is here at the intersection between tectonic processes and superposed climatic and sedimentary processes in the Himalayan orogenic wedge, where my investigation is focused on. The study area is the intermontane Kangra Basin in the northwestern Sub-Himalaya, because the characteristics of the different Himalayan morphotectonic provinces are well developed, the area is part of a region strongly influenced by monsoonal forcing, and the existence of numerous fluvial terraces provides excellent strain markers to assess deformation processes within the Himalayan orogenic wedge. In addition, being located in front of the Dhauladhar Range the region is characterized by pronounced gradients in past and present-day erosion and sediment processes associated with repeatedly changing climatic conditions. In light of these conditions I analysed climate-driven late Pleistocene-Holocene sediment cycles in this tectonically active region, which may be responsible for triggering the tectonic re-organization within the Himalayan orogenic wedge, leading to out-of-sequence thrusting, at least since early Holocene.
The Kangra Basin is bounded by the MBT and the Sub-Himalayan Jwalamukhi Thrust (JMT) in the north and south, respectively and transiently stores sediments derived from the Dhauladhar Range. The Basin contains ~200-m-thick conglomerates reflecting two distinct aggradation phases; following aggradation, several fluvial terraces were sculpted into these fan deposits. 10Be CRN surface exposure dating of these terrace levels provides an age of 53.4±3.2 ka for the highest-preserved terrace (AF1); subsequently, this surface was incised until ~15 ka, when the second fan (AF2) began to form. AF2 fan aggradation was superseded by episodic Holocene incision, creating at least four terrace levels. We find a correlation between variations in sediment transport and ∂18O records from regions affected by the Indian Summer Monsoon (ISM). During strengthened ISMs sand post-LGM glacial retreat, aggradation occurred in the Kangra Basin, likely due to high sediment flux, whereas periods of a weakened ISM coupled with lower sediment supply coincided with renewed re-incision.
However, the evolution of fluvial terraces along Sub-Himalayan streams in the Kangra sector is also forced by tectonic processes. Back-tilted, folded terraces clearly document tectonic activity of the JMT. Offset of one of the terrace levels indicates a shortening rate of 5.6±0.8 to 7.5±1.0 mm.a-1 over the last ~10 ka. Importantly, my study reveals that late Pleistocene/Holocene out-of-sequence thrusting accommodates 40-60% of the total 14±2 mm.a-1 shortening partitioned throughout the Sub-Himalaya. Importantly, the JMT records shortening at a lower rate over longer timescales hints towards out-of-sequence activity within the Sub-Himalaya. Re-activation of the JMT could be related to changes in the tectonic stress field caused by large-scale sediment removal from the basin. I speculate that the deformation processes of the Sub-Himalaya behave according to the predictions of critical wedge model and assume the following: While >200m of sediment aggradation would trigger foreland-ward propagation of the deformation front, re-incision and removal of most of the stored sediments (nearly 80-85% of the optimum basin-fill) would again create a sub-critical condition of the wedge taper and trigger the retreat of the deformation front.
While tectonism is responsible for the longer-term processes of erosion associated with steepening hillslopes, sediment cycles in this environment are mainly the result of climatic forcing. My new 10Be cosmogenic nuclide exposure dates and a synopsis of previous studies show the late Pleistocene to Holocene alluvial fills and fluvial terraces studied here record periodic fluctuations of sediment supply and transport capacity on timescales of 1000-100000 years. To further evaluate the potential influence of climate change on these fluctuations, I compared the timing of aggradation and incision phases recorded within remnant alluvial fans and terraces with continental climate archives such as speleothems in neighboring regions affected by monsoonal precipitation. Together with previously published OSL ages yielding the timing of aggradation, I find a correlation between variations in sediment transport with oxygen-isotope records from regions affected by the Indian Summer Monsoon (ISM). Accordingly, during periods of increased monsoon intensity (transitions from dry and cold to wet and warm periods – MIS4 to MIS3 and MIS2 to MIS1) (MIS=marine isotope stage) and post-Last Glacial Maximum glacial retreat, aggradation occurred in the Kangra Basin, likely due to high sediment flux. Conversely, periods of weakened monsoon intensity or lower sediment supply coincide with re-incision of the existing basin-fill.
Finally, my study entails part of a low-temperature thermochronology study to assess the youngest exhumation history of the Dhauladhar Range. Zircon helium (ZHe) ages and existing low-temperature data sets (ZHe, apatite fission track (AFT)) across this range, together with 3D thermokinematic modeling (PECUBE) reveals constraints on exhumation and activity of the range-bounding Main Boundary Thrust (MBT) since at least mid-Miocene time. The modeling results indicate mean slip rates on the MBT-fault ramp of ~2 – 3 mm.a-1 since its activation. This has lead to the growth of the >5-km-high frontal Dhauladhar Range and continuous deep-seated exhumation and erosion. The obtained results also provide interesting constraints of deformation patterns and their variation along strike. The results point towards the absence of the time-transient ‘mid-crustal ramp’ in the basal decollement and
duplexing of the Lesser Himalayan sequence, unlike the nearby regions or even the central Nepal domain. A fraction of convergence (~10-15%) is accommodated along the deep-seated MBT-ramp, most likely merging into the MHT. This finding is crucial for a rigorous assessment of the overall level of tectonic activity in the Himalayan morphotectonic provinces as it contradicts recently-published geodetic shortening estimates. In these studies, it has been proposed that the total Himalayan shortening in the NW Himalaya is accommodated within the Sub-Himalaya whereas no tectonic activity is assigned to the MBT.
Understanding the rates and processes of denudation is key to unraveling the dynamic processes that shape active orogens. This includes decoding the roles of tectonic and climate-driven processes in the long-term evolution of high- mountain landscapes in regions with pronounced tectonic activity and steep climatic and surface-process gradients. Well-constrained denudation rates can be used to address a wide range of geologic problems. In steady-state landscapes, denudation rates are argued to be proportional to tectonic or isostatic uplift rates and provide valuable insight into the tectonic regimes underlying surface denudation. The use of denudation rates based on terrestrial cosmogenic nuclide (TCN) such as 10Beryllium has become a widely-used method to quantify catchment-mean denudation rates. Because such measurements are averaged over timescales of 102 to 105 years, they are not as susceptible to stochastic changes as shorter-term denudation rate estimates (e.g., from suspended sediment measurements) and are therefore considered more reliable for a comparison to long-term processes that operate on geologic timescales. However, the impact of various climatic, biotic, and surface processes on 10Be concentrations and the resultant denudation rates remains unclear and is subject to ongoing discussion. In this thesis, I explore the interaction of climate, the biosphere, topography, and geology in forcing and modulating denudation rates on catchment to orogen scales.
There are many processes in highly dynamic active orogens that may effect 10Be concentrations in modern river sands and therefore impact 10Be-derived denudation rates. The calculation of denudation rates from 10Be concentrations, however, requires a suite of simplifying assumptions that may not be valid or applicable in many orogens. I investigate how these processes affect 10Be concentrations in the Arun Valley of Eastern Nepal using 34 new 10Be measurements from the main stem Arun River and its tributaries. The Arun Valley is characterized by steep gradients in climate and topography, with elevations ranging from <100 m asl in the foreland basin to >8,000 asl in the high sectors to the north. This is coupled with a five-fold increase in mean annual rainfall across strike of the orogen. Denudation rates from tributary samples increase toward the core of the orogen, from <0.2 to >5 mm/yr from the Lesser to Higher Himalaya. Very high denudation rates (>2 mm/yr), however, are likely the result of 10Be TCN dilution by surface and climatic processes, such as large landsliding and glaciation, and thus may not be representative of long-term denudation rates. Mainstem Arun denudation rates increase downstream from ~0.2 mm/yr at the border with Tibet to 0.91 mm/yr at its outlet into the Sapt Kosi. However, the downstream 10Be concentrations may not be representative of the entire upstream catchment. Instead, I document evidence for downstream fining of grains from the Tibetan Plateau, resulting in an order-of-magnitude apparent decrease in the measured 10Be concentration.
In the Arun Valley and across the Himalaya, topography, climate, and vegetation are strongly interrelated. The observed increase in denudation rates at the transition from the Lesser to Higher Himalaya corresponds to abrupt increases in elevation, hillslope gradient, and mean annual rainfall. Thus, across strike (N-S), it is difficult to decipher the potential impacts of climate and vegetation cover on denudation rates. To further evaluate these relationships I instead took advantage of an along-strike west-to-east increase of mean annual rainfall and vegetation density in the Himalaya. An analysis of 136 published 10Be denudation rates from along strike of the revealed that median denudation rates do not vary considerably along strike of the Himalaya, ~1500 km E-W. However, the range of denudation rates generally decreases from west to east, with more variable denudation rates in the northwestern regions of the orogen than in the eastern regions. This denudation rate variability decreases as vegetation density increases (R=- 0.90), and increases proportionately to the annual seasonality of vegetation (R=0.99). Moreover, rainfall and vegetation modulate the relationship between topographic steepness and denudation rates such that in the wet, densely vegetated regions of the Himalaya, topography responds more linearly to changes in denudation rates than in dry, sparsely vegetated regions, where the response of topographic steepness to denudation rates is highly nonlinear. Understanding the relationships between denudation rates, topography, and climate is also critical for interpreting sedimentary archives. However, there is a lack of understanding of how terrestrial organic matter is transported out of orogens and into sedimentary archives. Plant wax lipid biomarkers derived from terrestrial and marine sedimentary records are commonly used as paleo- hydrologic proxy to help elucidate these problems. I address the issue of how to interpret the biomarker record by using the plant wax isotopic composition of modern suspended and riverbank organic matter to identify and quantify organic matter source regions in the Arun Valley. Topographic and geomorphic analysis, provided by the 10Be catchment-mean denudation rates, reveals that a combination of topographic steepness (as a proxy for denudation) and vegetation density is required to capture organic matter sourcing in the Arun River.
My studies highlight the importance of a rigorous and careful interpretation of denudation rates in tectonically active orogens that are furthermore characterized by strong climatic and biotic gradients. Unambiguous information about these issues is critical for correctly decoding and interpreting the possible tectonic and climatic forces that drive erosion and denudation, and the manifestation of the erosion products in sedimentary archives.
The lakes in the Kenyan Rift Valley offer the unique opportunity to study a wide range of hydrochemical environmental conditions, ranging from freshwater to highly saline and alkaline lakes. Because little is known about the hydro- and biogeochemical conditions in the underlying lake sediments, it was the aim of this study to extend the already existing data sets with data from porewater and biomarker analyses. Additionally, reduced sulphur compounds and sulphate reduction rates in the sediment were determined. The new data was used to examine the anthropogenic and microbial influence on the lakes sediments as well as the influence of the water chemistry on the degradation and preservation of organic matter in the sediment column. The lakes discussed in this study are: Logipi, Eight (a small crater lake in the region of Kangirinyang), Baringo, Bogoria, Naivasha, Oloiden, and Sonachi.
The biomarker compositions were similar in all studied lake sediments; nevertheless, there were some differences between the saline and freshwater lakes. One of those differences is the occurrence of a molecule related to β-carotene, which was only found in the saline lakes. This molecule most likely originates from cyanobacteria, single-celled organisms which are commonly found in saline lakes. In the two freshwater lakes, stigmasterol, a sterol characteristic for freshwater algae, was found. In this study, it was shown that Lakes Bogoria and Sonachi can be used for environmental reconstructions with biomarkers, because the absence of oxygen at the lake bottoms slowed the degradation process. Other lakes, like for example Lake Naivasha, cannot be used for such reconstructions, because of the large anthropogenic influence. But the biomarkers proved to be a useful tool to study those anthropogenic influences. Additionally, it was observed that horizons with a high concentration of elemental sulphur can be used as temporal markers. Those horizons were deposited during times when the lake levels were very low. The sulphur was deposited by microorganisms which are capable of anoxygenic photosynthesis or sulphide oxidation.
Despite a large number of dated glacial landforms in the Himalaya, the ice extent during the global Last Glacial Maximum (LGM) from 19 to 23 ka is only known to first order. New cosmogenic Be-10 exposure ages from well-preserved glacially polished surfaces, combined with published data, and an improved production rate scaling model allow reconstruction of the LGM ice extent and subsequent deglaciation in the Chandra Valley of NW India. We show that a >1000 m thick valley glacier retreated >150 km within a few thousand years after the onset of LGM deglaciation. By comparing the recession of the Chandra Valley Glacier and other Himalayan glaciers with those of Northern and Southern Hemisphere glaciers, we demonstrate that post-LGM deglaciation was similar and nearly finished prior to the Bolling/Allerod interstadial. Our study supports the view that many Himalayan glaciers advanced during the LGM, likely in response to global variations in temperature.
Elevation-dependent changes in n-alkane delta D and soil GDGTs across the South Central Andes
(2016)
Surface uplift of large plateaus may significantly influence regional climate and more specifically precipitation patterns and temperature, sometimes complicating paleoaltimetry interpretations. Thus, understanding the topographic evolution of tectonically active mountain belts benefits from continued development of reliable proxies to reduce uncertainties in paleoaltimetry reconstructions. Lipid biomarker-based proxies provide a novel approach to stable isotope paleoaltimetry and complement authigenic or pedogenic mineral proxy materials, in particular outside semi-arid climate zones where soil carbonates are not abundant but (soil) organic matter has a high preservation potential. Here we present delta D values of soil-derived n-alkanes and mean annual air temperature (MAT) estimates based on branched glycerol dialkyl glycerol tetraether (brGDGT) distributions to assess their potential for paleoelevation reconstructions in the southern central Andes. We analyzed soil samples across two environmental and hydrological gradients that include a hillslope (26-28 degrees S) and a valley (22-24 degrees S) transect on the windward flanks of Central Andean Eastern Cordillera in NW Argentina. Our results show that present-day n-alkane delta D values and brGDGT-based MAT estimates are both linearly related with elevation and in good agreement with present-day climate conditions. Soil n-alkanes show a delta D lapse rate (A(delta D)) of -1.64 parts per thousand/100 m (R-2 = 0.91, p < 0.01) at the hillslope transect, within the range of delta D lapse rates from precipitation and surface waters in other tropical regions in the Andes like the Eastern Cordillera in Colombia and Bolivia and the Equatorial and Peruvian Andes. BrGDGT-derived soil temperatures are similar to monitored winter temperatures in the region and show a lapse rate of Delta T = -0.51 degrees C/100 m (R-2 = 0.91, p < 0.01), comparable with lapse rates from in situ soil temperature measurements, satellite derived land-surface temperatures at this transect, and weather stations from the Eastern Cordillera at similar latitude. As a result of an increasing leeward sampling position along the valley transect lapse rates are biased towards lower values and display higher scatter (Delta(delta D) = -0.9 parts per thousand/100 m, R-2 = 0.76, p < 0.01 and Delta T = -0.19 degrees C/100 m, R-2 = 0.48, p < 0.05). Despite this higher complexity, they are in line with lapse rates from stream-water samples and in situ soil temperature measurements along the same transect. Our results demonstrate that both soil n-alkane delta D values and MAT reconstructions based on brGDGTs distributions from the hillslope transect (Delta(delta D) = -1.64 parts per thousand/100 m, R-2 = 0.91, p < 0.01 and Delta T = -0.51 degrees C/100 m, R-2 = 0.91, p < 0.01) track the direct effects of orography on precipitation and temperature and hence the combined effects of local and regional hydrology as well as elevation. (C) 2016 Elsevier B.V. All rights reserved.
This paper is mainly based on field work carried out on the Messinian deposits of the Adana Basin ( southern Turkey), as well as on the interpretation of seismic reflection profiles to understand 3D geometries of the basin fill. Chronostratigraphic constraints for the Messinian deposits are from micropaleontological studies on foraminifera, ostracods, and calcareous nannofossils, recently carried out on the Messinian deposits of the Adana Basin. Our results indicate that this basin developed in a marginal area strictly related to the Mediterranean realm. The Messinian deposits of the Adana Basin record all the main steps of the Messinian Salinity Crisis ( MSC) that affected the Mediterranean area at the end of the Miocene. The new stratigraphic model for the Messinian deposits of the Adana Basin provided in this work gives new insights into both the MSC and the Taurus petroleum system. Despite their complete correspondence with the MSC, the Messinian deposits of the Adana Basin show some differences with respect to the current conceptual model for the MSC. For example, in the current conceptual model for the MSC, only one regional erosional surface ( MES) characterizes the MSC deposits. In the Adana Basin, two regional erosional surfaces, named MES1 and MES2, separate the Messinian deposits related to the MSC in Lower Evaporites, Resedimented Lower Evaporites ( RLE), and upper Messinian continental deposits containing a late Lago-Mare ostracod assemblage ( mainly fluvial coarse-grained and fine-grained sediments). In some places, Brecciated Limestones lie just above the MES1 and beneath the RLE. In addition, the RLE are thought to be related to the same step that brought to the Messinian halite deposition throughout the Mediterranean, pointing to a hyperhaline environment. In contrast, the fine-grained deposits of the RLE of the Adana Basin show the occurrence of Parathetyan brackish ostracod fauna ( early Lago-Mare ostracod assemblages), which defines an oligohaline depositional environment for the RLE. In terms of hydrocarbon prospecting, the Messinian evaporites of the Adana Basin have been considered as a perfect seal for the active Taurus petroleum system. Our results show that due to the complex stratigraphic architecture of the basin fill and the occurrence of two regional erosional surfaces ( MES1 and MES2), the Messinian evaporites are discontinuously present both in surface and in the subsurface of the Adana Basin. However, seal properties in the Adana Basin could be found in the Lower Pliocene deep marine clays of the Avadan Formation. This work leads to suggest a new stratigraphical model for the Messinian deposits of the Adana Basin, allowing us to amend the classical scheme with respect to the Messinian, and to officially define some new formations within the stratigraphy of the Adana Basin.
Deciphering the response of sediment routing systems to climatic forcing is fundamental for understanding the impacts of climate change on landscape evolution. In the Kangra Basin (northwest Sub-Himalaya, India), upper Pleistocene to Holocene alluvial fills and fluvial terraces record periodic fluctuations of sediment supply and transport capacity on timescales of 10(3) to 10(5) yr. To evaluate the potential influence of climate change on these fluctuations, we compare the timing of aggradation and incision phases recorded within remnant alluvial fans and terraces with climate archives. New surface-exposure dating of six terrace levels with in-situ cosmogenic Be-10 indicates the onset of incision phases. Two terrace surfaces from the highest level (T1) sculpted into the oldest preserved alluvial fan (AF1) date back to 53.4 +/- 3.2 ka and 43.0 +/- 2.7 ka (1 sigma). T2 surfaces sculpted into the remnants of AF1 have exposure ages of 18.6 +/- 1.2 ka and 15.3 +/- 0.9 ka, while terraces sculpted into the upper Pleistocene-Holocene fan (AF2) provide ages of 9.3 +/- 0.4 ka (T3), 7.1 +/- 0.4 ka (T4), 5.2 +/- 0.4 ka (T5) and 3.6 +/- 0.2 ka (T6). Together with previously published OSL ages yielding the timing of aggradation, we find a correlation between variations in sediment transport with oxygen-isotope records from regions affected by the Indian Summer Monsoon. During periods of increased monsoon intensity and post-Last Glacial Maximum glacial retreat, aggradation occurred in the Kangra Basin, likely due to high sediment flux, whereas periods of weakened monsoon intensity or lower sediment supply coincide with incision. (C) 2016 Elsevier B.V. All rights reserved.
The evolution of through-going normal-fault arrays from initial nucleation to growth and subsequent interaction and mechanical linkage is well documented in many extensional provinces. Over time, these processes. lead to predictable spatial and temporal variations in the amount and rate of displacement accumulated along strike of individual fault segments, which should be manifested in the patterns of footwall exhumation. Here, we investigate the along-strike and vertical distribution of low-temperature apatite (U-Th)/He (AHe) cooling ages along the bounding fault system, the Livingstone fault, of the Karonga Basin of the northern Malawi Rift. The fault evolution and linkage from rift initiation to the present day has been previously constrained through investigations of the hanging wall basin fill. The new cooling ages from the footwall of the Livingstone fault can be related to the adjacent depocentre evolution and across a relay zone between two palaeo-fault segments. Our data are complimented by published apatite fission track (AFT) data and reveal significant variation in rock cooling history along-strike: the centre of the footwall yields younger cooling ages than the former tips of earlier fault segments that are now linked. This suggests that low-temperature thermochronology can detect fault interactions along strike. That these former segment boundaries are preserved within exhumed footwall rocks is a function of the relatively recent linkage of the system. Our study highlights that changes in AHe (and potentially AFT) ages associated with the along-strike displacement profile can occur over relatively short horizontal distances (of a few kilometres). This is fundamentally important in the assessment of the vertical cooling history of footwalls in extensional systems: temporal differences in the rate of tectonically driven exhumation at a given location along fault strike may be of greater importance in controlling changes in rates of vertical exhumation than commonly invoked climatic fluctuations. (C) 2016 Elsevier B.V. All rights reserved.
The southernmost thrust of the Himalayan orogenic wedge that separates the foreland from the orogen, the Main Frontal Thrust, is thought to accommodate most of the ongoing crustal shortening in the Sub-Himalaya. Steepened longitudinal river profile segments, terrace offsets, and back-tilted fluvial terraces within the Kangra reentrant of the NW Sub-Himalaya suggest Holocene activity of the Jwalamukhi Thrust (JMT) and other thrust faults that may be associated with strain partitioning along the toe of the Himalayan wedge. To assess the shortening accommodated by the JMT, we combine morphometric terrain analyses with in situ Be-10-based surface-exposure dating of the deformed terraces. Incision into upper Pleistocene sediments within the Kangra Basin created two late Pleistocene terrace levels (T1 and T2). Subsequent early Holocene aggradation shortly before similar to 10ka was followed by episodic reincision, which created four cut-and-fill terrace levels, the oldest of which (T3) was formed at 10.10.9ka. A vertical offset of 445m of terrace T3 across the JMT indicates a shortening rate of 5.60.8 to 7.51.1mma(-1) over the last similar to 10ka. This result suggests that thrusting along the JMT accommodates 40-60% of the total Sub-Himalayan shortening in the Kangra reentrant over the Holocene. We speculate that this out-of-sequence shortening may have been triggered or at least enhanced by late Pleistocene and Holocene erosion of sediments from the Kangra Basin.
Rainfall in the central Andes associated with the South American Monsoon and the South American Low-Level Jet results from orographic effects on atmospheric circulation exerted by the Andean Plateau and the Eastern Cordillera. However, despite its importance for South American climate, no reliable records exist that allow decoding the evolution of thresholds and interactions between Andean topography and atmospheric circulation, especially regarding the onset of humid conditions in the inherently dry southern central Andes. Here, we employ multi-proxy isotope data of lipid biomarkers, pedogenic carbonates and volcanic glass from the Eastern Cordillera of NW Argentina and present the first long-term evapotranspiration record. We find that regional eco-hydrology and vegetation changes are associated with initiation of moisture transport via the South American Low-Level Jet at 7.6 Ma, and subsequent lateral growth of the orogen at 6.5 Ma. Our results highlight that topographically induced changes in atmospheric circulation patterns, not global climate change, were responsible for late Miocene environmental change in this part of the southern hemisphere. This suggests that mountain building over time fundamentally controlled habitat evolution along the central Andes.
Recent studies have shown that the 1976-77 global climate shift strongly affected the South American climate. In our study, we observed a link between this climate shift and river-discharge variability in the subtropical Southern Central Andes. We analyzed the daily river-discharge time series between 1940 and 1999 from small to medium mountain drainage basins (10(2)-10(4) km(2) ) across a steep climatic and topographic gradient. We document that the discharge frequency distribution changed significantly, with higher percentiles exhibiting more pronounced trends. A change point between 1971 and 1977 marked an intensification of the hydrological cycle, which resulted in increased river discharge. In the upper Rio Bermejo basin of the northernmost Argentine Andes, the mean annual discharge increased by 40% over 7 years. Our findings are important for flood risk management in areas impacted by the 1976-77 climate shift; discharge frequency distribution analysis provides important insights into the variability of the hydrological cycle in the Andean realm.
Extreme hydro-meteorological events, such as severe droughts or heavy rainstorms, constitute primary manifestations of climate variability and exert a critical impact on the natural environment and human society. This is particularly true for high-mountain areas, such as the eastern flank of the southern Central Andes of NW Argentina, a region impacted by deep convection processes that form the basis of extreme events, often resulting in floods, a variety of mass movements, and hillslope processes. This region is characterized by pronounced E-W gradients in topography, precipitation, and vegetation cover, spanning low to medium-elevation, humid and densely vegetated areas to high-elevation, arid and sparsely vegetated environments. This strong E-W gradient is mirrored by differences in the efficiency of surface processes, which mobilize and transport large amounts of sediment through the fluvial system, from the steep hillslopes to the intermontane basins and further to the foreland. In a highly sensitive high-mountain environment like this, even small changes in the spatiotemporal distribution, magnitude and rates of extreme events may strongly impact environmental conditions, anthropogenic activity, and the well-being of mountain communities and beyond. However, although the NW Argentine Andes comprise the catchments for the La Plata river that traverses one of the most populated and economically relevant areas of South America, there are only few detailed investigations of climate variability and extreme hydro-meteorological events.
In this thesis, I focus on deciphering the spatiotemporal variability of rainfall and river discharge, with particular emphasis on extreme hydro-meteorological events in the subtropical southern Central Andes of NW Argentina during the past seven decades. I employ various methods to assess and quantify statistically significant trend patterns of rainfall and river discharge, integrating high-quality daily time series from gauging stations (40 rainfall and 8 river discharge stations) with gridded datasets (CPC-uni and TRMM 3B42 V7), for the period between 1940 and 2015. Evidence for a general intensification of the hydrological cycle at intermediate elevations (~ 0.5 – 3 km asl) at the eastern flank of the southern Central Andes is found both from rainfall and river-discharge time-series analysis during the period from 1940 to 2015. This intensification is associated with the increase of the annual total amount of rainfall and the mean annual discharge. However, most pronounced trends are found at high percentiles, i.e. extreme hydro-meteorological events, particularly during the wet season from December to February.An important outcome of my studies is the recognition of a rapid increase in the amount of river discharge during the period between 1971 and 1977, most likely linked to the 1976-77 global climate shift, which is associated with the North Pacific Ocean sea surface temperature variability. Interestingly, after this rapid increase, both rainfall and river discharge decreased at low and intermediate elevations along the eastern flank of the Andes. In contrast, during the same time interval, at high elevations, extensive areas on the arid Puna de Atacama plateau have recorded increasing annual rainfall totals. This has been associated with more intense extreme hydro-meteorological events from 1979 to 2014. This part of the study reveals that low-, intermediate, and high-elevation sectors in the Andes of NW Argentina respond differently to changing climate conditions.
Possible forcing mechanisms of the pronounced hydro-meteorological variability observed in the study area are also investigated. For the period between 1940 and 2015, I analyzed modes of oscillation of river discharge from small to medium drainage basins (102 to 104 km2), located on the eastern flank of the orogen. First, I decomposed the relevant monthly time series using the Hilbert-Huang Transform, which is particularly appropriate for non-stationary time series that result from non-linear natural processes. I observed that in the study region discharge variability can be described by five quasi-periodic oscillatory modes on timescales varying from 1 to ~20 years. Secondly, I tested the link between river-discharge variations and large-scale climate modes of variability, using different climate indices, such as the BEST ENSO (Bivariate El Niño-Southern Oscillation Time-series) index. This analysis reveals that, although most of the variance on the annual timescale is associated with the South American Monsoon System, a relatively large part of river-discharge variability is linked to Pacific Ocean variability (PDO phases) at multi-decadal timescales (~20 years). To a lesser degree, river discharge variability is also linked to the Tropical South Atlantic (TSA) sea surface temperature anomaly at multi-annual timescales (~2-5 years).
Taken together, these findings exemplify the high degree of sensitivity of high-mountain environments with respect to climatic variability and change. This is particularly true for the topographic transitions between the humid, low-moderate elevations and the semi-arid to arid highlands of the southern Central Andes. Even subtle changes in the hydro-meteorological regime of these areas of the mountain belt react with major impacts on erosional hillslope processes and generate mass movements that fundamentally impact the transport capacity of mountain streams. Despite more severe storms in these areas, the fluvial system is characterized by pronounced variability of the stream power on different timescales, leading to cycles of sediment aggradation, the loss of agriculturally used land and severe impacts on infrastructure.
The collision of bathymetric anomalies, such as oceanic spreading centers, at convergent plate margins can profoundly affect subduction dynamics, magmatism, and the structural and geomorphic evolution of the overriding plate. The Southern Patagonian Andes of South America are a prime example for sustained oceanic ridge collision and the successive formation and widening of an extensive asthenospheric slab window since the Middle Miocene. Several of the predicted upper-plate geologic manifestations of such deep-seated geodynamic processes have been studied in this region, but many topics remain highly debated. One of the main controversial topics is the interpretation of the regional low-temperature thermochronology exhumational record and its relationship with tectonic and/or climate-driven processes, ultimately manifested and recorded in the landscape evolution of the Patagonian Andes. The prominent along-strike variance in the topographic characteristics of the Andes, combined with coupled trends in low-temperature thermochronometer cooling ages have been interpreted in very contrasting ways, considering either purely climatic (i.e. glacial erosion) or geodynamic (slab-window related) controlling factors.
This thesis focuses on two main aspects of these controversial topics. First, based on field observations and bedrock low-temperature thermochronology data, the thesis addresses an existing research gap with respect to the neotectonic activity of the upper plate in response to ridge collision - a mechanism that has been shown to affect the upper plate topography and exhumational patterns in similar tectonic settings. Secondly, the qualitative interpretation of my new and existing thermochronological data from this region is extended by inverse thermal modelling to define thermal histories recorded in the data and evaluate the relative importance of surface vs. geodynamic factors and their possible relationship with the regional cooling record.
My research is centered on the Northern Patagonian Icefield (NPI) region of the Southern Patagonian Andes. This site is located inboard of the present-day location of the Chile Triple Junction - the juncture between the colliding Chile Rise spreading center and the Nazca and Antarctic Plates along the South American convergent margin. As such this study area represents the region of most recent oceanic-ridge collision and associated slab window formation. Importantly, this location also coincides with the abrupt rise in summit elevations and relief characteristics in the Southern Patagonian Andes. Field observations, based on geological, structural and geomorphic mapping, are combined with bedrock apatite (U-Th)/He and apatite fission track (AHe and AFT) cooling ages sampled along elevation transects across the orogen. This new data reveals the existence of hitherto unrecognized neotectonic deformation along the flanks of the range capped by the NPI.
This deformation is associated with the closely spaced oblique collision of successive oceanic-ridge segments in this region over the past 6 Ma. I interpret that this has caused a crustal-scale partitioning of deformation and the decoupling, margin-parallel migration, and localized uplift of a large crustal sliver (the NPI block) along the subduction margin. The location of this uplift coincides with a major increase of summit elevations and relief at the northern edge of the NPI massif. This mechanism is compatible with possible extensional processes along the topographically subdued trailing edge of the NPI block as documented by very recent and possibly still active normal faulting. Taken together, these findings suggest a major structural control on short-wavelength variations in topography in the Southern Patagonian Andes - the region affected by ridge collision and slab window formation.
The second research topic addressed here focuses on using my new and existing bedrock low-temperature cooling ages in forward and inverse thermal modeling. The data was implemented in the HeFTy and QTQt modeling platforms to constrain the late Cenozoic thermal history of the Southern Patagonian Andes in the region of the most recent upper-plate sectors of ridge collision. The data set combines AHe and AFT data from three elevation transects in the region of the Northern Patagonian Icefield. Previous similar studies claimed far-reaching thermal effects of the approaching ridge collision and slab window to affect patterns of Late Miocene reheating in the modelled thermal histories. In contrast, my results show that the currently available data can be explained with a simpler thermal history than previously proposed. Accordingly, a reheating event is not needed to reproduce the observations. Instead, the analyzed ensemble of modelled thermal histories defines a Late Miocene protracted cooling and Pliocene-to-recent stepwise exhumation. These findings agree with the geological record of this region. Specifically, this record indicates an Early Miocene phase of active mountain building associated with surface uplift and an active fold-and-thrust belt, followed by a period of stagnating deformation, peneplanation, and lack of synorogenic deposition in the Patagonian foreland. The subsequent period of stepwise exhumation likely resulted from a combination of pulsed glacial erosion and coeval neotectonic activity. The differences between the present and previously published interpretation of the cooling record can be reconciled with important inconsistencies of previously used model setup. These include mainly the insufficient convergence of the models and improper assumptions regarding the geothermal conditions in the region. This analysis puts a methodological emphasis on the prime importance of the model setup and the need for its thorough examination to evaluate the robustness of the final outcome.
Stable-isotopic and sedimentary records from the orogenic Puna Plateau of NW Argentina and adjacent intermontane basins to the east furnish a unique late Cenozoic record of range uplift and ensuing paleoenvironmental change in the south-central Andes. Today, focused precipitation in this region occurs along the eastern, windward flanks of the Eastern Cordillera and Sierras Pampeanas ranges, while the orogen interior constitutes high-elevation regions with increasingly arid conditions in a westward direction. As in many mountain belts, such hydrologic and topographic gradients are commonly mirrored by a systematic relationship between the oxygen and hydrogen stable isotope ratios of meteoric water and elevation. The glass fraction of isotopically datable volcanic ash intercalated in sedimentary sequences constitutes an environmental proxy that retains a signal of the hydrogen-isotopic composition of ancient precipitation. This isotopic composition thus helps to elucidate the combined climatic and tectonic processes associated with topographic growth, which ultimately controls the spatial patterns of precipitation in mountain belts. However, between 25.5 and 27 degrees S present-day river-based hydrogen isotope lapse rates are very low, possibly due to deep-convective seasonal storms that dominate runoff. If not accounted for, the effects of such conditions on moisture availability in the past may lead to misinterpretations of proxy-records of rainfall. Here, we present hydrogen-isotope data of volcanic glass (delta Dg), extracted from 34 volcanic ash layers in different sedimentary basins of the Eastern Cordillera and the Sierras Pampeanas. Combined with previously published delta Dg records and our refined U-Pb and (U-Th)/He zircon geochronology on 17 tuff samples, we demonstrate hydrogen-isotope variations associated with paleoenvironmental change in the Angastaco Basin, which evolved from a contiguous foreland to a fault-bounded intermontane basin during the late Mio-Pliocene. We unravel the environmental impact of Mio-Pliocene topographic growth and associated orographic effects on long-term hydrogen-isotope records of rainfall in the south-central Andes, and potentially identify temporal variations in regional isotopic lapse rates that may also apply to other regions with similar topographic boundary conditions. (C) 2016 Elsevier B.V. All rights reserved.
The Cenozoic East African Rift System (EARS) is an exceptional example of active continental extension, providing opportunities for furthering our understanding of hydrocarbon plays within rifts. It is divided into structurally distinct western and eastern branches. The western branch comprises deep rift basins separated by transfer zones, commonly localised onto pre-existing structures, offering good regional scale hydrocarbon traps. At a basin-scale, local discrete inherited structures might also play an important role on fault localisation and hydrocarbon distribution. Here, we consider the evolution of the Central basin of the Malawi Rift, in particular the influence of pre-existing structural fabrics. Integrating basin-scale multichannel 2D, and high resolution seismic datasets we constrain the border, Mlowe-Nkhata, fault system (MNF) to the west of the basin and smaller Mbamba fault (MF) to the east and document their evolution. Intra basin structures define a series of horsts, which initiated as convergent transfers, along the basin axis. The horsts are offset along a NE SW striking transfer fault parallel to and along strike of the onshore Karoo (Permo-Triassic) Ruhuhu graben. Discrete pre-existing structures probably determined its location and, oriented obliquely to the extension orientation it accommodated predominantly strike-slip deformation, with more slowly accrued dip-slip. To the north of this transfer fault, the overall basin architecture is asymmetric, thickening to the west throughout; while to the south, an initially symmetric graben architecture became increasingly asymmetric in sediment distribution as strain localised onto the western MNF. The presence of the axial horst increasingly focussed sediment supply to the west. As the transfer fault increased its displacement, so this axial supply was interrupted, effectively starving the south-east while ponding sediments between the western horst margin and the transfer fault. This asymmetric bathymetry and partitioned sedimentation continues to the present-day, overprinting the early basin symmetry and configuration. Sediments deposited earlier become increasingly dissected and fault juxtapositions changed at a small (10-100 m) scale. The observed influence of basin-scale transfer faults on sediment dispersal and fault compartmentalization due to pre-existing structures oblique to the extension orientation is relevant to analogous exploration settings. (C) 2016 Elsevier Ltd. All rights reserved.
Terrigenous sediment supply, marine transport, and depositional processes along tectonically active margins are key to decoding turbidite successions as potential archives of climatic and seismic forcings. Sequence stratigraphic models predict coarse-grained sediment delivery to deep-marine sites mainly during sea-level fall and lowstand. Marine siliciclastic deposition during transgressions and highstands has been attributed to sustained connectivity between terrigenous sources and marine sinks facilitated by narrow shelves. To decipher the controls on Holocene highstand turbidite deposition, we analyzed 12 sediment cores from spatially discrete, coeval turbidite systems along the Chile margin (29 degrees-40 degrees S) with changing climatic and geomorphic characteristics but uniform changes in sea level. Sediment cores from intraslope basins in north-central Chile (29 degrees-33 degrees S) offshore a narrow to absent shelf record a shut-off of turbidite deposition during the Holocene due to postglacial aridification. In contrast, core sites in south-central Chile (36 degrees-40 degrees S) offshore a wide shelf record frequent turbidite deposition during highstand conditions. Two core sites are linked to the Biobio river-canyon system and receive sediment directly from the river mouth. However, intraslope basins are not connected via canyons to fluvial systems but yield even higher turbidite frequencies. High sediment supply combined with a wide shelf and an undercurrent moving sediment toward the shelf edge appear to control Holocene turbidite sedimentation and distribution. Shelf undercurrents may play an important role in lateral sediment transport and supply to the deep sea and need to be accounted for in sediment-mass balances.
The scaling behavior of rainfall has been extensively studied both in terms of event magnitudes and in terms of spatial extents of the events. Different heavy-tailed distributions have been proposed as candidates for both instances, but statistically rigorous treatments are rare. Here we combine the domains of event magnitudes and event area sizes by a spatiotemporal integration of 3-hourly rain rates corresponding to extreme events derived from the quasi-global high-resolution rainfall product Tropical Rainfall Measuring Mission 3B42. A maximum likelihood evaluation reveals that the distribution of spatiotemporally integrated extreme rainfall cluster sizes over the oceans is best described by a truncated power law, calling into question previous statements about scale-free distributions. The observed subpower law behavior of the distribution's tail is evaluated with a simple generative model, which indicates that the exponential truncation of an otherwise scale-free spatiotemporal cluster size distribution over the oceans could be explained by the existence of land masses on the globe.
Coastal fens like the nature reserve "Hutelmoor und Heiligensee" (north-eastern Germany) are important landscape elements along the southern Baltic coast, which exchange fresh water and brackish water with the Baltic Sea. These exchange processes can be understood as experiments with a natural tracer, which may be used to investigate the hydrologic behaviour of these fen systems. With the establishment of coastal protection measures such as dunes and dikes, the installation of surface drainage and, more recently, also nature conservation measures, the hydrologic regime of these coastal wetlands has constantly altered over the last centuries.The rehabilitated wetland "Hutelnnoor und Heiligensee" is suitable for an analysis of hydrologic change as it has been monitored over the time period since nature conservation measures started in the 1990s. Collected data sets included observation of groundwater levels and electrical conductivities, weather data as well as discharge at the outlet of the drainage catchment. In this article, as a second part of the dual publication, processes and quantified process magnitudes have been identified that govern the salt balance of the study area including its variability in space and time. It was detected that - over the period of rehabilitation - salt water entered the catchment with an episodic storm surge by wave overtopping of dunes in 1995. The intruded brackish water was then diluted, which was a slow process extending over decades. It was governed by local groundwater recharge from precipitation and the inflow of relatively fresh groundwater from the hinterland. It is concluded that salt inputs from the Baltic Sea provide a natural tracer of hydrological processes, which can be readily monitored via electrical conductivity measurements.
Hydrological and climatological controls on radiocarbon concentrations in a tropical stalagmite
(2016)
Precisely-dated stalagmites are increasingly important archives for the reconstruction of terrestrial paleoclimate at very high temporal resolution. In-depth understanding of local conditions at the cave site and of the processes driving stalagmite deposition is of paramount importance for interpreting proxy signals incorporated in stalagmite carbonate. Here we present a sub-decadally resolved dead carbon fraction (DCF) record for a stalagmite from Yok Balum Cave (southern Belize). The record is coupled to parallel stable carbon isotope (delta C-13) and U/Ca measurements, as well as radiocarbon (C-14) measurements from soils overlying the cave system. Using a karst carbon cycle model we disentangle the importance of soil and karst processes on stalagmite DCF incorporation, revealing a dominant host rock dissolution control on total DCF. Covariation between DCF, delta C-13, and U/Ca indicates that karst processes are a common driver of all three parameters, suggesting possible use of delta C-13 and trace element ratios to independently quantify DCF variability. A statistically significant multi-decadal lag of variable length exists between DCF and reconstructed solar activity, suggesting that solar activity influenced regional precipitation in Mesoamerica over the past 1500 years, but that the relationship was non-static. Although the precise nature of the observed lag is unclear, solar-induced changes in North Atlantic oceanic and atmospheric dynamics may play a role. (C) 2016 Elsevier Ltd. All rights reserved.
Water scarcity, adaption on climate change, and risk assessment of droughts and floods are critical topics for science and society these days. Monitoring and modeling of the hydrological cycle are a prerequisite to understand and predict the consequences for weather and agriculture. As soil water storage plays a key role for partitioning of water fluxes between the atmosphere, biosphere, and lithosphere, measurement techniques are required to estimate soil moisture states from small to large scales.
The method of cosmic-ray neutron sensing (CRNS) promises to close the gap between point-scale and remote-sensing observations, as its footprint was reported to be 30 ha. However, the methodology is rather young and requires highly interdisciplinary research to understand and interpret the response of neutrons to soil moisture. In this work, the signal of nine detectors has been systematically compared, and correction approaches have been revised to account for meteorological and geomagnetic variations. Neutron transport simulations have been consulted to precisely characterize the sensitive footprint area, which turned out to be 6--18 ha, highly local, and temporally dynamic. These results have been experimentally confirmed by the significant influence of water bodies and dry roads. Furthermore, mobile measurements on agricultural fields and across different land use types were able to accurately capture the various soil moisture states. It has been further demonstrated that the corresponding spatial and temporal neutron data can be beneficial for mesoscale hydrological modeling. Finally, first tests with a gyrocopter have proven the concept of airborne neutron sensing, where increased footprints are able to overcome local effects.
This dissertation not only bridges the gap between scales of soil moisture measurements. It also establishes a close connection between the two worlds of observers and modelers, and further aims to combine the disciplines of particle physics, geophysics, and soil hydrology to thoroughly explore the potential and limits of the CRNS method.
The northern Tien Shan of Kyrgyzstan and Kazakhstan has been affected by a series of major earthquakes in the late 19th and early 20th centuries. To assess the significance of such a pulse of strain release in a continental interior, it is important to analyze and quantify strain release over multiple time scales. We have undertaken paleoseismological investigations at two geomorphically distinct sites (Panfilovkoe and Rot Front) near the Kyrgyz capital Bishkek. Although located near the historic epicenters, both sites were not affected by these earthquakes. Trenching was accompanied by dating stratigraphy and offset surfaces using luminescence, radiocarbon, and Be-10 terrestrial cosmogenic nuclide methods. At Rot Front, trenching of a small scarp did not reveal evidence for surface rupture during the last 5000 years. The scarp rather resembles an extensive debris-flow lobe. At Panfilovkoe, we estimate a Late Pleistocene minimum slip rate of 0.2 +/- 0.1 mm/a, averaged over at least two, probably three earthquake cycles. Dip-slip reverse motion along segmented, moderately steep faults resulted in hanging wall collapse scarps during different events. The most recent earthquake occurred around 3.6 +/- 1.3 kyr ago (1 sigma), with dip-slip offsets between 1.2 and 1.4 m. We calculate a probabilistic paleomagnitude to be between 6.7 and 7.2, which is in agreement with regional data from the Kyrgyz range. The morphotectonic signals in the northern Tien Shan are a prime example of deformation in a tectonically active intracontinental mountain belt and as such can help understand the longer-term coevolution of topography and seismogenic processes in similar structural settings worldwide.
Ecosystems' exposure to climate change - Modeling as support for nature conservation management
(2016)
Storm runoff from the Marikina River Basin frequently causes flood events in the Philippine capital region Metro Manila. This paper presents and evaluates a system to predict short-term runoff from the upper part of that basin (380km(2)). It was designed as a possible component of an operational warning system yet to be installed. For the purpose of forecast verification, hindcasts of streamflow were generated for a period of 15 months with a time-continuous, conceptual hydrological model. The latter was fed with real-time observations of rainfall. Both ground observations and weather radar data were tested as rainfall forcings. The radar-based precipitation estimates clearly outperformed the raingauge-based estimates in the hydrological verification. Nevertheless, the quality of the deterministic short-term runoff forecasts was found to be limited. For the radar-based predictions, the reduction of variance for lead times of 1, 2 and 3hours was 0.61, 0.62 and 0.54, respectively, with reference to a no-forecast scenario, i.e. persistence. The probability of detection for major increases in streamflow was typically less than 0.5. Given the significance of flood events in the Marikina Basin, more effort needs to be put into the reduction of forecast errors and the quantification of remaining uncertainties.
Ongoing climate change is known to cause an increase in the frequency and amplitude of local temperature and precipitation extremes in many regions of the Earth. While gradual changes in the climatological conditions have already been shown to strongly influence plant flowering dates, the question arises if and how extremes specifically impact the timing of this important phenological phase. Studying this question calls for the application of statistical methods that are tailored to the specific properties of event time series. Here, we employ event coincidence analysis, a novel statistical tool that allows assessing whether or not two types of events exhibit similar sequences of occurrences in order to systematically quantify simultaneities between meteorological extremes and the timing of the flowering of four shrub species across Germany. Our study confirms previous findings of experimental studies by highlighting the impact of early spring temperatures on the flowering of the investigated plants. However, previous studies solely based on correlation analysis do not allow deriving explicit estimates of the strength of such interdependencies without further assumptions, a gap that is closed by our analysis. In addition to direct impacts of extremely warm and cold spring temperatures, our analysis reveals statistically significant indications of an influence of temperature extremes in the autumn preceding the flowering.
The 2014 April 1, M-w 8.1 Iquique (Chile) earthquake struck in the Northern Chile seismic gap. With a rupture length of less than 200 km, it left unbroken large segments of the former gap. Early studies were able to model the main rupture features but results are ambiguous with respect to the role of aseismic slip and left open questions on the remaining hazard at the Northern Chile gap. A striking observation of the 2014 earthquake has been its extensive preparation phase, with more than 1300 events with magnitude above M-L 3, occurring during the 15 months preceding the main shock. Increasing seismicity rates and observed peak magnitudes accompanied the last three weeks before the main shock. Thanks to the large data sets of regional recordings, we assess the precursor activity, compare foreshocks and aftershocks and model rupture preparation and rupture effects. To tackle inversion challenges for moderate events with an asymmetric network geometry, we use full waveforms techniques to locate events, map the seismicity rate and derive source parameters, obtaining moment tensors for more than 300 events (magnitudes M-w 4.0-8.1) in the period 2013 January 1-2014 April 30. This unique data set of fore- and aftershocks is investigated to distinguish rupture process models and models of strain and stress rotation during an earthquake. Results indicate that the spatial distributions of foreshocks delineated the shallower part of the rupture areas of the main shock and its largest aftershock, well matching the spatial extension of the aftershocks cloud. Most moment tensors correspond to almost pure double couple thrust mechanisms, consistent with the slab orientation. Whereas no significant differences are observed among thrust mechanisms in different areas, nor among thrust foreshocks and aftershocks, the early aftershock sequence is characterized by the presence of normal fault mechanisms, striking parallel to the trench but dipping westward. These events likely occurred in the shallow wedge structure close to the slab interface and are consequence of the increased extensional stress in this region after the largest events. The overall stress inversion result suggests a minor stress rotation after the main shock, but a significant release of the deviatoric stress. The temporal change in the distribution of focal mechanisms can also be explained in terms of the spatial heterogeneity of the stress field: under such interpretation, the potential of a large megathrust earthquake breaking a larger segment offshore Northern Chile remains high.
This case study evaluates the suitability of radar-based quantitative precipitation estimates (QPEs) for the simulation of streamflow in the Marikina River Basin (MRB), the Philippines. Hourly radar-based QPEs were produced from reflectivity that had been observed by an S-band radar located about 90 km from the MRB. Radar data processing and precipitation estimation were carried out using the open source library wradlib. To assess the added value of the radar-based QPE, we used spatially interpolated rain gauge observations (gauge-only (GO) product) as a benchmark. Rain gauge observations were also used to quantify rainfall estimation errors at the point scale. At the point scale, the radar-based QPE outperformed the GO product in 2012, while for 2013, the performance was similar. For both periods, estimation errors substantially increased from daily to the hourly accumulation intervals. Despite this fact, both rainfall estimation methods allowed for a good representation of observed streamflow when used to force a hydrological simulation model of the MRB. Furthermore, the results of the hydrological simulation were consistent with rainfall verification at the point scale: the radar-based QPE performed better than the GO product in 2012, and equivalently in 2013. Altogether, we could demonstrate that, in terms of streamflow simulation, the radar-based QPE can perform as good as or even better than the GO product - even for a basin such as the MRB which has a comparatively dense rain gauge network. This suggests good prospects for using radar-based QPE to simulate and forecast streamflow in other parts of the Philippines where rain gauge networks are not as dense.
The increase in atmospheric methane concentration, which is determined by an imbalance between its sources and sinks, has led to investigations of the methane cycle in various environments. Aquatic environments are of an exceptional interest due to their active involvement in methane cycling worldwide and in particular in areas sensitive to climate change. Furthermore, being connected with each other aquatic environments form networks that can be spread on vast areas involving marine, freshwater and terrestrial ecosystems. Thus, aquatic systems have a high potential to translate local or regional environmental and subsequently ecosystem changes to a bigger scale. Many studies neglect this connectivity and focus on individual aquatic or terrestrial ecosystems.
The current study focuses on environmental controls of the distribution and aerobic oxidation of methane at the example of two aquatic ecosystems. These ecosystems are Arctic fresh water bodies and the Elbe estuary which represent interfaces between freshwater-terrestrial and freshwater-marine environments, respectively.
Arctic water bodies are significant atmospheric sources of methane. At the same time the methane cycle in Arctic water bodies is strongly affected by the surrounding permafrost environment, which is characterized by high amounts of organic carbon. The results of this thesis indicate that the methane concentrations in Arctic lakes and streams substantially vary between each other being regulated by local landscape features (e.g. floodplain area) and the morphology of the water bodies (lakes, streams and river). The highest methane concentrations were detected in the lake outlets and in a floodplain lake complex. In contrast, the methane concentrations measured at different sites of the Lena River did not vary substantially. The lake complexes in comparison to the Lena River, thus, appear as more individual and heterogeneous systems with a pronounced imprint of the surrounding soil environment. Furthermore, connected with each other Arctic aquatic environments have a large potential to transport methane from methane-rich water bodies such as streams and floodplain lakes to aquatic environments relatively poor in methane such as the Lena River.
Estuaries represent hot spots of oceanic methane emissions. Also, estuaries are intermediate zones between methane-rich river water and methane depleted marine water. Substantiated through this thesis at the example of the Elbe estuary, the methane distribution in estuaries, however, cannot entirely be described by the conservative mixing model i.e. gradual decrease from the freshwater end-member to the marine water end-member. In addition to the methane-rich water from the Elbe River mouth substantial methane input occurs from tidal flats, areas of significant interaction between aquatic and terrestrial environments. Thus, this study demonstrates the complex interactions and their consequences for the methane distribution within estuaries. Also it reveals how important it is to investigate estuaries at larger spatial scales.
Methane oxidation (MOX) rates are commonly correlated with methane concentrations. This was shown in previous research studies and was substantiated by the present thesis. In detail, the highest MOX rates in the Arctic water bodies were detected in methane-rich streams and in the floodplain area while in the Elbe estuary the highest MOX rates were observed at the coastal stations. However, in these bordering environments, MOX rates are affected not only via the regulation through methane concentrations. The MOX rates in the Arctic lakes were shown to be also dependent on the abundance and community composition of methane-oxidising bacteria (MOB), that in turn are controlled by local landscape features (regardless of the methane concentrations) and by the transport of MOB between neighbouring environments. In the Elbe estuary, the MOX rates in addition to the methane concentrations are largely affected by the salinity, which is in turn regulated by the mixing of fresh- and marine waters. The magnitude of the salinity impact on MOX rates thereby depends on the MOB community composition and on the rate of the salinity change.
This study extends our knowledge of environmental controls of methane distribution and aerobic methane oxidation in aquatic environments. It also illustrates how important it is to investigate complex ecosystems rather than individual ecosystems to better understand the functioning of whole biomes.
Large, compressional earthquakes cause surface uplift aswell as widespread mass wasting. Knowledge of their trade-off is fragmentary. Combining a seismologically consistent model of earthquake-triggered landsliding and an analytical solution of coseismic surface displacement, we assess how the mass balance of single earthquakes and earthquake sequences depends on fault size and other geophysical parameters. We find that intermediate size earthquakes (M-w 6-7.3) may cause more erosion than uplift, controlled primarily by seismic source depth and landscape steepness, and less so by fault dip and rake. Such earthquakes can limit topographic growth, but our model indicates that both smaller and larger earthquakes (M-w < 6, M-w > 7.3) systematically cause mountain building. Earthquake sequences with a Gutenberg-Richter distribution have a greater tendency to lead to predominant erosion, than repeating earthquakes of the same magnitude, unless a fault can produce earthquakes with M-w > 8 or more.
Zircon (ZrSiO4), hafnon (HfSiO4) and five intermediate compositions were synthesized from a Pb silicate melt. The resulting crystals were 20-300 mu m in size and displayed sector and growth zoning. Raman spectra were acquired at locations in the sample for which preceding electron microprobe (EMP) analyses revealed sufficient compositional homogeneity. The dataset documents shifts of Raman bands with changing composition. In this study, bands that have previously not been reported were found for the intermediate compositions and for pure hafnon, in particular at wavenumbers less than 200 cm(-1). For these external modes, the dataset provides new insight into the compositional dependence of their frequencies. Density-functional theory calculations support the observations and are used for a detailed interpretation of the spectra. The pitfalls of the EMP analysis along the zircon-hafnon join are highlighted.
Solute concentration variability is of fundamental importance for the chemical and ecological state of streams. It is often closely related to discharge variability and can be characterized in terms of a solute export regime. Previous studies, especially in lowland catchments, report that nitrate is often exported with an accretion pattern of increasing concentrations with increasing discharge. Several modeling approaches exist to predict the export regime of solutes from the spatial relationship of discharge generating zones with solute availability in the catchment. For a small agriculturally managed lowland catchment in central Germany, we show that this relationship is controlled by the depth to groundwater table and its temporal dynamics. Principal component analysis of groundwater level time series from wells distributed throughout the catchment allowed derivation of a representative groundwater level time series that explained most of the discharge variability. Groundwater sampling revealed consistently decreasing nitrate concentrations with an increasing thickness of the unsaturated zone. The relationships of depth to groundwater table to discharge and to nitrate concentration were parameterized and integrated to successfully model catchment discharge and nitrate export on the basis of groundwater level variations alone. This study shows that intensive and uniform agricultural land use likely results in a clear and consistent concentration-depth relationship of nitrate, which can be utilized in simple approaches to predict stream nitrate export dynamics at the catchment scale. (C) 2016 Elsevier Ltd. All rights reserved.
Lake Towuti is a tectonic basin, surrounded by ultramafic rocks. Lateritic soils form through weathering and deliver abundant iron (oxy)hydroxides but very little sulfate to the lake and its sediment. To characterize the sediment biogeochemistry, we collected cores at three sites with increasing water depth and decreasing bottom water oxygen concentrations. Microbial cell densities were highest at the shallow site a feature we attribute to the availability of labile organic matter (OM) and the higher abundance of electron acceptors due to oxic bottom water conditions. At the two other sites, OM degradation and reduction processes below the oxycline led to partial electron acceptor depletion. Genetic information preserved in the sediment as extracellular DNA (eDNA) provided information on aerobic and anaerobic heterotrophs related to Nitrospirae. Chloroflexi, and Therrnoplasmatales. These taxa apparently played a significant role in the degradation of sinking OM. However, eDNA concentrations rapidly decreased with core depth. Despite very low sulfate concentrations, sulfate-reducing bacteria were present and viable in sediments at all three sites, as confirmed by measurement of potential sulfate reduction rates. Microbial community fingerprinting supported the presence of taxa related to Deltaproteobacteria and Firmicutes with demonstrated capacity for iron and sulfate reduction. Concomitantly, sequences of Ruminococcaceae, Clostridiales, and Methanornicrobiales indicated potential for fermentative hydrogen and methane production. Such first insights into ferruginous sediments showed that microbial populations perform successive metabolisms related to sulfur, iron, and methane. In theory, iron reduction could reoxidize reduced sulfur compounds and desorb OM from iron minerals to allow remineralization to methane. Overall, we found that biogeochemical processes in the sediments can be linked to redox differences in the bottom waters of the three sites, like oxidant concentrations and the supply of labile OM. At the scale of the lacustrine record, our geomicrobiological study should provide a means to link the extant subsurface biosphere to past environments.
Earthquakes deform Earth's surface, building long-lasting topographic features and contributing to landscape and mountain formation.
However, seismic waves produced by earthquakes may also destabilize hillslopes, leading to large amounts of soil and bedrock moving downslope. Moreover, static deformation and shaking are suspected to damage the surface bedrock and therefore alter its future properties, affecting hydrological and erosional dynamics. Thus, earthquakes participate both in mountain building and stimulate directly or indirectly their erosion. Moreover, the impact of earthquakes on hillslopes has important implications for the amount of sediment and organic matter delivered to rivers, and ultimately to oceans, during episodic catastrophic seismic crises, the magnitude of life and property losses associated with landsliding, the perturbation and recovery of landscape properties after shaking, and the long term topographic evolution of mountain belts. Several of these aspects have been addressed recently through individual case studies but additional data compilation as well as theoretical or numerical modelling are required to tackle these issues in a more systematic and rigorous manner.
This dissertation combines data compilation of earthquake characteristics, landslide mapping, and seismological data interpretation with physically-based modeling in order to address how earthquakes impact on erosional processes and landscape evolution. Over short time scales (10-100 s) and intermediate length scales (10 km), I have attempted to improve our understanding and ability to predict the amount of landslide debris triggered by seismic shaking in epicentral areas. Over long time scales (1-100 ky) and across a mountain belt (100 km) I have modeled the competition between erosional unloading and building of topography associated with earthquakes. Finally, over intermediate time scales (1-10 y) and at the hillslope scale (0.1-1 km) I have collected geomorphological and seismological data that highlight persistent effects of earthquakes on landscape properties and behaviour.
First, I compiled a database on earthquakes that produced significant landsliding, including an estimate of the total landslide volume and area, and earthquake characteristics such as seismic moment and source depth. A key issue is the accurate conversion of landslide maps into volume estimates. Therefore I also estimated how amalgamation - when mapping errors lead to the bundling of multiple landslide into a single polygon - affects volume estimates from various earthquake-induced landslide inventories and developed an algorithm to automatically detect this artifact. The database was used to test a physically-based prediction of the total landslide area and volume caused by earthquakes, based on seismological scaling relationships and a statistical description of the landscape properties. The model outperforms empirical fits in accuracy, with 25 out of 40 cases well predicted, and allows interpretation of many outliers in physical terms. Apart from seismological complexities neglected by the model I found that exceptional rock strength properties or antecedent conditions may explain most outliers.
Second, I assessed the geomorphic effects of large earthquakes on landscape dynamics by surveying the temporal evolution of precipitation-normalized landslide rate. I found strongly elevated landslide rates following earthquakes that progressively recover over 1 to 4 years, indicating that regolith strength drops and recovers. The relaxation is clearly non-linear for at least one case, and does not seem to correlate with coseismic landslide reactivation, water table level increase or tree root-system recovery. I suggested that shallow bedrock is damaged by the earthquake and then heals on annual timescales. Such variations in ground strength must be translated into shallow subsurface seismic velocities that are increasingly surveyed with ambient seismic noise correlations. With seismic noise autocorrelation I computed the seismic velocity in the epicentral areas of three earthquakes where I constrained a change in landslide rate. We found similar recovery dynamics and timescales, suggesting that seismic noise correlation techniques could be further developed to meaningfully assess ground strength variations for landscape dynamics. These two measurements are also in good agreement with the temporal dynamics of post-seismic surface displacement measured by GPS. This correlation suggests that the surface healing mechanism may be driven by tectonic deformation, and that the surface regolith and fractured bedrock may behave as a granular media that slowly compacts as it is sheared or vibrated.
Last, I compared our model of earthquake-induced landsliding with a standard formulation of surface deformation caused by earthquakes to understand which parameters govern the competition between the building and destruction of topography caused by earthquakes. In contrast with previous studies I found that very large (Mw>8) earthquakes always increase the average topography, whereas only intermediate (Mw ~ 7) earthquakes in steep landscapes may reduce topography. Moreover, I illustrated how the net effect of earthquakes varies with depth or landscape steepness implying a complex and ambivalent role through the life of a mountain belt. Further I showed that faults producing a Gutenberg-Richter distribution of earthquake sizes, will limit topography over a larger range of fault sizes than faults producing repeated earthquakes with a characteristic size.
We present a new, seismologically consistent expression for the total area and volume of populations of earthquake-triggered landslides. This model builds on a set of scaling relationships between key parameters, such as landslide spatial density, seismic ground acceleration, fault length, earthquake source depth, and seismic moment. To assess the model we have assembled and normalized a catalog of landslide inventories for 40 shallow, continental earthquakes. Low landscape steepness causes systematic overprediction of the total area and volume of landslides. When this effect is accounted for, the model predicts the total landslide volume of 63% of 40 cases to within a factor 2 of the volume estimated from observations (R-2 = 0.76). The prediction of total landslide area is also sensitive to the landscape steepness, but less so than the total volume, and it appears to be sensitive to controls on the landslide size-frequency distribution, and possibly the shaking duration. Some outliers are likely associated with exceptionally strong rock mass in the epicentral area, while others may be related to seismic source complexities ignored by the model. However, the close match between prediction and estimate for about two thirds of cases in our database suggests that rock mass strength is similar in many cases and that our simple seismic model is often adequate, despite the variety of lithologies and tectonic settings covered. This makes our expression suitable for integration into landscape evolution models and application to the anticipation or rapid assessment of secondary hazards associated with earthquakes.
Bedrock landslides, by excavating deep below saprolite-rock interfaces, create conditions for weathering in which all mineral phases in a lithology are initially unweathered within landslide deposits. As a result, the most labile phases dominate the weathering immediately after mobilisation and during a transient period of depletion. This mode of dissolution can strongly alter the overall output of solutes from catchments and their contribution to global chemical cycles if landslide-derived material is retained in catchments for extended periods after mass wasting.
The South American Andes are frequently exposed to intense rainfall events with varying moisture sources and precipitation-forming processes. In this study, we assess the spatiotemporal characteristics and geographical origins of rainfall over the South American continent. Using high-spatiotemporal resolution satellite data (TRMM 3B42 V7), we define four different types of rainfall events based on their (1) high magnitude, (2) long temporal extent, (3) large spatial extent, and (4) high magnitude, long temporal and large spatial extent combined. In a first step, we analyze the spatiotemporal characteristics of these events over the entire South American continent and integrate their impact for the main Andean hydrologic catchments. Our results indicate that events of type 1 make the overall highest contributions to total seasonal rainfall (up to 50%). However, each consecutive episode of the infrequent events of type 4 still accounts for up to 20% of total seasonal rainfall in the subtropical Argentinean plains. In a second step, we employ complex network theory to unravel possibly non-linear and long-ranged climatic linkages for these four event types on the high-elevation Altiplano-Puna Plateau as well as in the main river catchments along the foothills of the Andes. Our results suggest that one to two particularly large squall lines per season, originating from northern Brazil, indirectly trigger large, long-lasting thunderstorms on the Altiplano Plateau. In general, we observe that extreme rainfall in the catchments north of approximately 20 degrees S typically originates from the Amazon Basin, while extreme rainfall at the eastern Andean foothills south of 20 degrees S and the Puna Plateau originates from southeastern South America.
Non-linear time series analysis of precipitation events using regional climate networks for Germany
(2016)
Synchronous occurrences of heavy rainfall events and the study of their relation in time and space are of large socio-economical relevance, for instance for the agricultural and insurance sectors, but also for the general well-being of the population. In this study, the spatial synchronization structure is analyzed as a regional climate network constructed from precipitation event series. The similarity between event series is determined by the number of synchronous occurrences. We propose a novel standardization of this number that results in synchronization scores which are not biased by the number of events in the respective time series. Additionally, we introduce a new version of the network measure directionality that measures the spatial directionality of weighted links by also taking account of the effects of the spatial embedding of the network. This measure provides an estimate of heavy precipitation isochrones by pointing out directions along which rainfall events synchronize. We propose a climatological interpretation of this measure in terms of propagating fronts or event traces and confirm it for Germany by comparing our results to known atmospheric circulation patterns.
The analytical evaluation of diurnal temperature variation in riverbed sediments provides detailed information on exchange fluxes between rivers and groundwater. The underlying assumption of the stationary, one-dimensional vertical flow field is frequently violated in natural systems where subsurface water flow often has a significant horizontal component. In this paper, we present a new methodology for identifying the geometry of the subsurface flow field using vertical temperature profiles. The statistical analyses are based on model optimisation and selection and are used to evaluate the shape of vertical amplitude ratio profiles. The method was applied to multiple profiles measured around in-stream geomorphological structures in a losing reach of a gravel bed river. The predominant subsurface flow field was systematically categorised in purely vertical and horizontal (hyporheic, parafluvial) components. The results highlight that river groundwater exchange flux at the head, crest and tail of geomorphological structures significantly deviated from the one-dimensional vertical flow, due to a significant horizontal component. The geometry of the subsurface water flow depended on the position around the geomorphological structures and on the river level. The methodology presented in this paper features great potential for characterising the spatial patterns and temporal dynamics of complex subsurface flow geometries by using measured temperature time series in vertical profiles. (C) 2016 Elsevier B.V. All rights reserved.
Im Graduiertenkolleg NatRiskChange der Universität Potsdam und anderen Forschungseinrichtungen werden beobachtete sowie zukünftig mögliche Veränderungen von Naturgefahren untersucht. Teil des strukturierten Doktorandenprogramms sind sogenannte Task-Force-Einsätze, bei denen die Promovierende zeitlich begrenzt ein aktuelles Ereignis auswerten. Im Zuge dieser Aktivität wurde die Sturzflut vom 29.05.2016 in Braunsbach (Baden-Württemberg) untersucht.
In diesem Bericht werden erste Auswertungen zur Einordnung der Niederschläge, zu den hydrologischen und geomorphologischen Prozessen im Einzugsgebiet des Orlacher Bachs sowie zu den verursachten Schäden beleuchtet.
Die Region war Zentrum extremer Regenfälle in der Größenordnung von 100 mm innerhalb von 2 Stunden. Das 6 km² kleine Einzugsgebiet hat eine sehr schnelle Reaktionszeit, zumal bei vorgesättigtem Boden. Im steilen Bachtal haben mehrere kleinere und größere Hangrutschungen über 8000 m³ Geröll, Schutt und Schwemmholz in das Gewässer eingetragen und möglicherweise kurzzeitige Aufstauungen und Durchbrüche verursacht. Neben den großen Wassermengen mit einer Abflussspitze in einer Größenordnung von 100 m³/s hat gerade die Geschiebefracht zu großen Schäden an den Gebäuden entlang des Bachlaufs in Braunsbach geführt.
Climatic controls on debris-flow activity and sediment aggradation: The Del Medio fan, NW Argentina
(2016)
In the Central Andes, several studies on alluvial terraces and valley fills have linked sediment aggradation to periods of enhanced sediment supply. However, debate continues over whether tectonic or climatic factors are most important in triggering the enhanced supply. The Del Medio catchment in the Humahuaca Basin (Eastern Cordillera, NW Argentina) is located within a transition zone between subhumid and arid climates and hosts the only active debris-flow fan within this intermontane valley. By combining Be-10 analyses of boulder and sediment samples within the Del Medio catchment, with regional morphometric measurements of nearby catchments, we identify the surface processes responsible for aggradation in the Del Medio fan and their likely triggers. We find that the fan surface has been shaped by debris flows and channel avulsions during the last 400 years. Among potential tectonic, climatic, and autogenic factors that might influence deposition, our analyses point to a combination of several favorable factors that drive aggradation. These are in particular the impact of occasional abundant rainfall on steep slopes in rock types prone to failure, located in a region characterized by relatively low rainfall amounts and limited transport capacity. These characteristics are primarily associated with the climatic transition zone between the humid foreland and the arid orogen interior, which creates an imbalance between sediment supply and sediment transfer. The conditions and processes that drive aggradation in the Del Medio catchment today may provide a modern analog for the conditions and processes that drove aggradation in other nearby tributaries in the past.
Fluvial fill terraces preserve sedimentary archives of landscape responses to climate change, typically over millennial timescales. In the Humahuaca Basin of NW Argentina (Eastern Cordillera, southern Central Andes), our 29 new optically stimulated luminescence ages of late Pleistocene fill terrace sediments demonstrate that the timing of past river aggradation occurred over different intervals on the western and eastern sides of the valley, despite their similar bedrock lithology, mean slopes, and precipitation. In the west, aggradation coincided with periods of increasing precipitation, while in the east, aggradation coincided with decreasing precipitation or more variable conditions. Erosion rates and grain size dependencies in our cosmogenic Be-10 analyses of modern and fill terrace sediments reveal an increased importance of landsliding compared to today on the west side during aggradation, but of similar importance during aggradation on the east side. Differences in the timing of aggradation and the Be-10 data likely result from differences in valley geometry, which causes sediment to be temporarily stored in perched basins on the east side. It appears as if periods of increasing precipitation triggered landslides throughout the region, which induced aggradation in the west, but blockage of the narrow bedrock gorges downstream from the perched basins in the east. As such, basin geometry and fluvial connectivity appear to strongly influence the timing of sediment movement through the system. For larger basins that integrate subbasins with differing geometries or degrees of connectivity (like Humahuaca), sedimentary responses to climate forcing are likely attenuated.
A link between chemical weathering and physical erosion exists at the catchment scale over a wide range of erosion rates(1,2). However, in mountain environments, where erosion rates are highest, weathering may be kinetically limited(3-5) and therefore decoupled from erosion. In active mountain belts, erosion is driven by bedrock landsliding(6) at rates that depend strongly on the occurrence of extreme rainfall or seismicity(7). Although landslides affect only a small proportion of the landscape, bedrock landsliding can promote the collection and slow percolation of surface runoff in highly fragmented rock debris and create favourable conditions for weathering. Here we show from analysis of surface water chemistry in the Southern Alps of New Zealand that weathering in bedrock landslides controls the variability in solute load of these mountain rivers. We find that systematic patterns in surface water chemistry are strongly associated with landslide occurrence at scales from a single hillslope to an entire mountain belt, and that landslides boost weathering rates and river solute loads over decades. We conclude that landslides couple erosion and weathering in fast-eroding uplands and, thus, mountain weathering is a stochastic process that is sensitive to climatic and tectonic controls on mass wasting processes.
This manuscript proposes a method to assess hydrological drought in semi-arid environments under high impoundment rate and applies it to the semi-arid Jaguaribe River basin in Brazil. It analyzes droughts (1) in the largest reservoir systems; (2) in the Upper Basin, considering 4744 reservoirs, 800 wells and almost 18,000 cisterns; and (3) in reservoirs of different sizes during multiyear droughts. Results show that the water demand is constrained in the basin; hydrological and meteorological droughts are often out of phase; there is a negative correlation between storage level and drought severity; and the small systems cannot cope with long-term droughts.
The Dead Sea region has faced substantial environmental challenges in recent decades, including water resource scarcity, similar to 1 m annual decreases in the water level, sinkhole development, ascending-brine freshwater pollution, and seismic disturbance risks. Natural processes are significantly affected by human interference as well as by climate change and tectonic developments over the long term. To get a deep understanding of processes and their interactions, innovative scientific approaches that integrate disciplinary research and education are required. The research project DESERVE (Helmholtz Virtual Institute Dead Sea Research Venue) addresses these challenges in an interdisciplinary approach that includes geophysics, hydrology, and meteorology. The project is implemented by a consortium of scientific institutions in neighboring countries of the Dead Sea (Israel, Jordan, Palestine Territories) and participating German Helmholtz Centres (KIT, GFZ, UFZ). A new monitoring network of meteorological, hydrological, and seismic/geodynamic stations has been established, and extensive field research and numerical simulations have been undertaken. For the first time, innovative measurement and modeling techniques have been applied to the extreme conditions of the Dead Sea and its surroundings. The preliminary results show the potential of these methods. First time ever performed eddy covariance measurements give insight into the governing factors of Dead Sea evaporation. High-resolution bathymetric investigations reveal a strong correlation between submarine springs and neo-tectonic patterns. Based on detailed studies of stratigraphy and borehole information, the extension of the subsurface drainage basin of the Dead Sea is now reliably estimated. Originality has been achieved in monitoring flash floods in an arid basin at its outlet and simultaneously in tributaries, supplemented by spatio-temporal rainfall data. Low-altitude, high resolution photogrammetry, allied to satellite image analysis and to geophysical surveys (e.g. shear-wave reflections) has enabled a more detailed characterization of sinkhole morphology and temporal development and the possible subsurface controls thereon. All the above listed efforts and scientific results take place with the interdisciplinary education of young scientists. They are invited to attend joint thematic workshops and winter schools as well as to participate in field experiments. (C) 2015 The Authors. Published by Elsevier B.V.
Oceanic lithospheric S-wave velocities from the analysis of P-wave polarization at the ocean floor
(2016)
Our knowledge of the absolute S-wave velocities of the oceanic lithosphere is mainly based on global surface wave tomography, local active seismic or compliance measurements using oceanic infragravity waves. The results of tomography give a rather smooth picture of the actual S-wave velocity structure and local measurements have limitations regarding the range of elastic parameters or the geometry of the measurement. Here, we use the P-wave polarization (apparent P-wave incidence angle) of teleseismic events to investigate the S-wave velocity structure of the oceanic crust and the upper tens of kilometres of the mantle beneath single stations. In this study, we present an up to our knowledge new relation of the apparent P-wave incidence angle at the ocean bottom dependent on the half-space S-wave velocity. We analyse the angle in different period ranges at ocean bottom stations (OBSs) to derive apparent S-wave velocity profiles. These profiles are dependent on the S-wave velocity as well as on the thickness of the layers in the subsurface. Consequently, their interpretation results in a set of equally valid models. We analyse the apparent P-wave incidence angles of an OBS data set which was collected in the Eastern Mid Atlantic. We are able to determine reasonable S-wave-velocity-depth models by a three-step quantitative modelling after a manual data quality control, although layer resonance sometimes influences the estimated apparent S-wave velocities. The apparent S-wave velocity profiles are well explained by an oceanic PREM model in which the upper part is replaced by four layers consisting of a water column, a sediment, a crust and a layer representing the uppermost mantle. The obtained sediment has a thickness between 0.3 and 0.9 km with S-wave velocities between 0.7 and 1.4 km s(-1). The estimated total crustal thickness varies between 4 and 10 km with S-wave velocities between 3.5 and 4.3 km s(-1). We find a slight increase of the total crustal thickness from similar to 5 to similar to 8 km towards the South in the direction of a major plate boundary, the Gloria Fault. The observed crustal thickening can be related with the known dominant compression in the vicinity of the fault. Furthermore, the resulting mantle S-wave velocities decrease from values around 5.5 to 4.5 km s(-1) towards the fault. This decrease is probably caused by serpentinization and indicates that the oceanic transform fault affects a broad region in the uppermost mantle. Conclusively, the presented method is useful for the estimation of the local S-wave velocity structure beneath ocean bottom seismic stations. It is easy to implement and consists of two main steps: (1) measurement of apparent P-wave incidence angles in different period ranges for real and synthetic data, and (2) comparison of the determined apparent S-wave velocities for real and synthetic data to estimate S-wave velocity-depth models.
Rapidly uplifting coastlines are frequently associated with convergent tectonic boundaries, like subduction zones, which are repeatedly breached by giant megathrust earthquakes. The coastal relief along tectonically active realms is shaped by the effect of sea-level variations and heterogeneous patterns of permanent tectonic deformation, which are accumulated through several cycles of megathrust earthquakes. However, the correlation between earthquake deformation patterns and the sustained long-term segmentation of forearcs, particularly in Chile, remains poorly understood. Furthermore, the methods used to estimate permanent deformation from geomorphic markers, like marine terraces, have remained qualitative and are based on unrepeatable methods. This contrasts with the increasing resolution of digital elevation models, such as Light Detection and Ranging (LiDAR) and high-resolution bathymetric surveys.
Throughout this thesis I study permanent deformation in a holistic manner: from the methods to assess deformation rates, to the processes involved in its accumulation. My research focuses particularly on two aspects: Developing methodologies to assess permanent deformation using marine terraces, and comparing permanent deformation with seismic cycle deformation patterns under different spatial scales along the M8.8 Maule earthquake (2010) rupture zone. Two methods are developed to determine deformation rates from wave-built and wave-cut terraces respectively. I selected an archetypal example of a wave-built terrace at Santa Maria Island studying its stratigraphy and recognizing sequences of reoccupation events tied with eleven radiocarbon sample ages (14C ages). I developed a method to link patterns of reoccupation with sea-level proxies by iterating relative sea level curves for a range of uplift rates. I find the best fit between relative sea-level and the stratigraphic patterns for an uplift rate of 1.5 +- 0.3 m/ka.
A Graphical User Interface named TerraceM® was developed in Matlab®. This novel software tool determines shoreline angles in wave-cut terraces under different geomorphic scenarios. To validate the methods, I select test sites in areas of available high-resolution LiDAR topography along the Maule earthquake rupture zone and in California, USA. The software allows determining the 3D location of the shoreline angle, which is a proxy for the estimation of permanent deformation rates. The method is based on linear interpolations to define the paleo platform and cliff on swath profiles. The shoreline angle is then located by intersecting these interpolations. The
accuracy and precision of TerraceM® was tested by comparing its results with previous assessments, and through an experiment with students in a computer lab setting at the University
of Potsdam.
I combined the methods developed to analyze wave-built and wave-cut terraces to assess regional patterns of permanent deformation along the (2010) Maule earthquake rupture. Wave-built terraces are tied using 12 Infra Red Stimulated luminescence ages (IRSL ages) and shoreline angles in wave-cut terraces are estimated from 170 aligned swath profiles. The comparison of coseismic slip, interseismic coupling, and permanent deformation, leads to three areas of high permanent uplift, terrace warping, and sharp fault offsets. These three areas correlate with regions of high slip and low coupling, as well as with the spatial limit of at least eight historical megathrust ruptures (M8-9.5). I propose that the zones of upwarping at Arauco and Topocalma reflect changes in frictional properties of the megathrust, which result in discrete boundaries for the propagation of mega earthquakes.
To explore the application of geomorphic markers and quantitative morphology in offshore areas I performed a local study of patterns of permanent deformation inferred from hitherto unrecognized drowned shorelines at the Arauco Bay, at the southern part of the (2010) Maule earthquake rupture zone. A multidisciplinary approach, including morphometry, sedimentology, paleontology, 3D morphoscopy, and a landscape Evolution Model is used to recognize, map, and assess local rates and patterns of permanent deformation in submarine environments. Permanent deformation patterns are then reproduced using elastic models to assess deformation rates of an active submarine splay fault defined as Santa Maria Fault System. The best fit suggests a reverse structure with a slip rate of 3.7 m/ka for the last 30 ka. The register of land level changes during the earthquake cycle at Santa Maria Island suggest that most of the deformation may be accrued through splay fault reactivation during mega earthquakes, like the (2010) Maule event. Considering a recurrence time of 150 to 200 years, as determined from historical and geological observations, slip between 0.3 and 0.7 m per event would be required to account for the 3.7 m/ka millennial slip rate. However, if the SMFS slips only every ~1000 years, representing a few megathrust earthquakes, then a slip of ~3.5 m per event would be required to account for the long- term rate. Such event would be equivalent to a magnitude ~6.7 earthquake capable to generate a local tsunami.
The results of this thesis provide novel and fundamental information regarding the amount of permanent deformation accrued in the crust, and the mechanisms responsible for this accumulation at millennial time-scales along the M8.8 Maule earthquake (2010) rupture zone. Furthermore, the results of this thesis highlight the application of quantitative geomorphology and the use of repeatable methods to determine permanent deformation, improve the accuracy of marine terrace assessments, and estimates of vertical deformation rates in tectonically active coastal areas. This is vital information for adequate coastal-hazard assessments and to anticipate realistic earthquake and tsunami scenarios.
The repeated occurrence of exceptional floods within a few years, such as the Rhine floods in 1993 and 1995 and the Elbe and Danube floods in 2002 and 2013, suggests that floods in Central Europe may be organized in flood-rich and flood-poor periods. This hypothesis is studied by testing the significance of temporal clustering in flood occurrence (peak-over-threshold) time series for 68 catchments across Germany for the period 1932-2005. To assess the robustness of the results, different methods are used: Firstly, the index of dispersion, which quantifies the departure from a homogeneous Poisson process, is investigated. Further, the time-variation of the flood occurrence rate is derived by non-parametric kernel implementation and the significance of clustering is evaluated via parametric and non-parametric tests. Although the methods give consistent overall results, the specific results differ considerably. Hence, we recommend applying different methods when investigating flood clustering. For flood estimation and risk management, it is of relevance to understand whether clustering changes with flood severity and time scale. To this end, clustering is assessed for different thresholds and time scales. It is found that the majority of catchments show temporal clustering at the 5% significance level for low thresholds and time scales of one to a few years. However, clustering decreases substantially with increasing threshold and time scale. We hypothesize that flood clustering in Germany is mainly caused by catchment memory effects along with intra- to inter-annual climate variability, and that decadal climate variability plays a minor role. (C) 2016 Elsevier B.V. All rights reserved.
To understand past flood changes in the Rhine catchment and in particular the role of anthropogenic climate change in extreme flows, an attribution study relying on a proper GCM (general circulation model) downscaling is needed. A downscaling based on conditioning a stochastic weather generator on weather patterns is a promising approach. This approach assumes a strong link between weather patterns and local climate, and sufficient GCM skill in reproducing weather pattern climatology. These presuppositions are unprecedentedly evaluated here using 111 years of daily climate data from 490 stations in the Rhine basin and comprehensively testing the number of classification parameters and GCM weather pattern characteristics. A classification based on a combination of mean sea level pressure, temperature, and humidity from the ERA20C reanalysis of atmospheric fields over central Europe with 40 weather types was found to be the most appropriate for stratifying six local climate variables. The corresponding skill is quite diverse though, ranging from good for radiation to poor for precipitation. Especially for the latter it was apparent that pressure fields alone cannot sufficiently stratify local variability. To test the skill of the latest generation of GCMs from the CMIP5 ensemble in reproducing the frequency, seasonality, and persistence of the derived weather patterns, output from 15 GCMs is evaluated. Most GCMs are able to capture these characteristics well, but some models showed consistent deviations in all three evaluation criteria and should be excluded from further attribution analysis.
During expedition 202 aboard the RV Sonne in 2009, 39 seafloor surface sediment sites were sampled over a wide sector of the North Pacific and adjoining Bering Sea. The data served to infer land-ocean linkages of terrigenous sediment supply in terms of major sources and modes of sediment transport within an over-regional context. This is based on an integrated approach dealing with grain-size analysis, bulk mineralogy and clay mineralogy in combination with statistical data evaluation (end-member modelling of grain-size data, fuzzy cluster analysis of mineralogical data). The findings on clay mineralogy served to update those of earlier work extracted from the literature. Today, two processes of terrigenous sediment supply prevail in the study area: far-distance aeolian sediment supply to the pelagic North Pacific, and hemipelagic sediment dispersal from nearby land sources via ocean currents along the continental margins and island arcs. Aeolian particles show the finest grain sizes (clay and fine silt), whereas hemipelagic sediments have high abundances of coarse silt. Exposed sites on seamounts and the continental slope are partly swept by strong currents, leading to residual enrichment of fine sand. Four sediment sources can be distinguished on the basis of distinct index minerals revealed by statistical data analysis: dust plumes from central Asia (quartz, illite), altered materials from the volcanic regions of Kamchatka and the Aleutian Arc (smectite), detritus from the Alaskan Cordillera (chlorite, hornblende), and fluvial detritus from far-eastern Siberia and the Alaska mainland (quartz, feldspar, illite). These findings confirm those of former studies but considerably expand the geographic range of this suite of proxies as far south as 39A degrees N in the open North Pacific. The present integrated methodological approach proved useful in identifying the major modern processes of terrigenous sediment supply to the study region. This aspect deserves attention in the selection of sediment core sites for future palaeoenvironmental reconstructions related to aeolian and glacial dynamics, as well as the recognition of palaeo-ocean circulation patterns in general.
In the past, floods were basically managed by flood control mechanisms. The focus was set on the reduction of flood hazard. The potential consequences were of minor interest. Nowadays river flooding is increasingly seen from the risk perspective, including possible consequences. Moreover, the large-scale picture of flood risk became increasingly important for disaster management planning, national risk developments and the (re-) insurance industry. Therefore, it is widely accepted that risk-orientated flood management ap-proaches at the basin-scale are needed. However, large-scale flood risk assessment methods for areas of several 10,000 km² are still in early stages. Traditional flood risk assessments are performed reach wise, assuming constant probabilities for the entire reach or basin. This might be helpful on a local basis, but where large-scale patterns are important this approach is of limited use. Assuming a T-year flood (e.g. 100 years) for the entire river network is unrealistic and would lead to an overestimation of flood risk at the large scale. Due to the lack of damage data, additionally, the probability of peak discharge or rainfall is usually used as proxy for damage probability to derive flood risk. With a continuous and long term simulation of the entire flood risk chain, the spatial variability of probabilities could be consider and flood risk could be directly derived from damage data in a consistent way.
The objective of this study is the development and application of a full flood risk chain, appropriate for the large scale and based on long term and continuous simulation. The novel approach of ‘derived flood risk based on continuous simulations’ is introduced, where the synthetic discharge time series is used as input into flood impact models and flood risk is directly derived from the resulting synthetic damage time series.
The bottleneck at this scale is the hydrodynamic simu-lation. To find suitable hydrodynamic approaches for the large-scale a benchmark study with simplified 2D hydrodynamic models was performed. A raster-based approach with inertia formulation and a relatively high resolution of 100 m in combination with a fast 1D channel routing model was chosen.
To investigate the suitability of the continuous simulation of a full flood risk chain for the large scale, all model parts were integrated into a new framework, the Regional Flood Model (RFM). RFM consists of the hydrological model SWIM, a 1D hydrodynamic river network model, a 2D raster based inundation model and the flood loss model FELMOps+r. Subsequently, the model chain was applied to the Elbe catchment, one of the largest catchments in Germany. For the proof-of-concept, a continuous simulation was per-formed for the period of 1990-2003. Results were evaluated / validated as far as possible with available observed data in this period. Although each model part introduced its own uncertainties, results and runtime were generally found to be adequate for the purpose of continuous simulation at the large catchment scale.
Finally, RFM was applied to a meso-scale catchment in the east of Germany to firstly perform a flood risk assessment with the novel approach of ‘derived flood risk assessment based on continuous simulations’. Therefore, RFM was driven by long term synthetic meteorological input data generated by a weather generator. Thereby, a virtual time series of climate data of 100 x 100 years was generated and served as input to RFM providing subsequent 100 x 100 years of spatially consistent river discharge series, inundation patterns and damage values. On this basis, flood risk curves and expected annual damage could be derived directly from damage data, providing a large-scale picture of flood risk. In contrast to traditional flood risk analysis, where homogenous return periods are assumed for the entire basin, the presented approach provides a coherent large-scale picture of flood risk. The spatial variability of occurrence probability is respected. Additionally, data and methods are consistent. Catchment and floodplain processes are repre-sented in a holistic way. Antecedent catchment conditions are implicitly taken into account, as well as physical processes like storage effects, flood attenuation or channel–floodplain interactions and related damage influencing effects. Finally, the simulation of a virtual period of 100 x 100 years and consequently large data set on flood loss events enabled the calculation of flood risk directly from damage distributions. Problems associated with the transfer of probabilities in rainfall or peak runoff to probabilities in damage, as often used in traditional approaches, are bypassed.
RFM and the ‘derived flood risk approach based on continuous simulations’ has the potential to provide flood risk statements for national planning, re-insurance aspects or other questions where spatially consistent, large-scale assessments are required.
Flood loss modeling is an important component for risk analyses and decision support in flood risk management. Commonly, flood loss models describe complex damaging processes by simple, deterministic approaches like depth-damage functions and are associated with large uncertainty. To improve flood loss estimation and to provide quantitative information about the uncertainty associated with loss modeling, a probabilistic, multivariable Bagging decision Tree Flood Loss Estimation MOdel (BT-FLEMO) for residential buildings was developed. The application of BT-FLEMO provides a probability distribution of estimated losses to residential buildings per municipality. BT-FLEMO was applied and validated at the mesoscale in 19 municipalities that were affected during the 2002 flood by the River Mulde in Saxony, Germany. Validation was undertaken on the one hand via a comparison with six deterministic loss models, including both depth-damage functions and multivariable models. On the other hand, the results were compared with official loss data. BT-FLEMO outperforms deterministic, univariable, and multivariable models with regard to model accuracy, although the prediction uncertainty remains high. An important advantage of BT-FLEMO is the quantification of prediction uncertainty. The probability distribution of loss estimates by BT-FLEMO well represents the variation range of loss estimates of the other models in the case study.
Flood generation at the scale of large river basins is triggered by the interaction of the hydrological pre-conditions and the meteorological event conditions at different spatial and temporal scales. This interaction controls diverse flood generating processes and results in floods varying in magnitude and extent, duration as well as socio-economic consequences. For a process-based understanding of the underlying cause-effect relationships, systematic approaches are required. These approaches have to cover the complete causal flood chain, including the flood triggering meteorological event in combination with the hydrological (pre-)conditions in the catchment, runoff generation, flood routing, possible floodplain inundation and finally flood losses.
In this thesis, a comprehensive probabilistic process-based understanding of the causes and effects of floods is advanced. The spatial and temporal dynamics of flood events as well as the geophysical processes involved in the causal flood chain are revealed and the systematic interconnections within the flood chain are deciphered by means of the classification of their associated causes and effects. This is achieved by investigating the role of the hydrological pre-conditions and the meteorological event conditions with respect to flood occurrence, flood processes and flood characteristics as well as their interconnections at the river basin scale.
Broadening the knowledge about flood triggers, which up to now has been limited to linking large-scale meteorological conditions to flood occurrence, the influence of large-scale pre-event hydrological conditions on flood initiation is investigated. Using the Elbe River basin as an example, a classification of soil moisture, a key variable of pre-event conditions, is developed and a probabilistic link between patterns of soil moisture and flood occurrence is established. The soil moisture classification is applied to continuously simulated soil moisture data which is generated using the semi-distributed conceptual rainfall-runoff model SWIM. Applying successively a principal component analysis and a cluster analysis, days of similar soil moisture patterns are identified in the period November 1951 to October 2003.
The investigation of flood triggers is complemented by including meteorological conditions described by a common weather pattern classification that represents the main modes of atmospheric state variability. The newly developed soil moisture classification thereby provides the basis to study the combined impact of hydrological pre-conditions and large-scale meteorological event conditions on flood occurrence at the river basin scale.
A process-based understanding of flood generation and its associated probabilities is attained by classifying observed flood events into process-based flood types such as snowmelt floods or long-rain floods. Subsequently, the flood types are linked to the soil moisture and weather patterns. Further understanding of the processes is gained by modeling of the complete causal flood chain, incorporating a rainfall-runoff model, a 1D/2D hydrodynamic model and a flood loss model. A reshuffling approach based on weather patterns and the month of their occurrence is developed to generate synthetic data fields of meteorological conditions, which drive the model chain, in order to increase the flood sample size. From the large number of simulated flood events, the impact of hydro-meteorological conditions on various flood characteristics is detected through the analysis of conditional cumulative distribution functions and regression trees.
The results show the existence of catchment-scale soil moisture patterns, which comprise of large-scale seasonal wetting and drying components as well as of smaller-scale variations related to spatially heterogeneous catchment processes. Soil moisture patterns frequently occurring before the onset of floods are identified. In winter, floods are initiated by catchment-wide high soil moisture, whereas in summer the flood-initiating soil moisture patterns are diverse and the soil moisture conditions are less stable in time. The combined study of both soil moisture and weather patterns shows that the flood favoring hydro-meteorological patterns as well as their interactions vary seasonally. In the analysis period, 18 % of the weather patterns only result in a flood in the case of preceding soil saturation. The classification of 82 past events into flood types reveals seasonally varying flood processes that can be linked to hydro-meteorological patterns. For instance, the highest flood potential for long-rain floods is associated with a weather pattern that is often detected in the presence of so-called ‘Vb’ cyclones. Rain-on-snow and snowmelt floods are associated with westerly and north-westerly wind directions. The flood characteristics vary among the flood types and can be reproduced by the applied model chain. In total, 5970 events are simulated. They reproduce the observed event characteristics between September 1957 and August 2002 and provide information on flood losses. A regression tree analysis relates the flood processes of the simulated events to the hydro-meteorological (pre-)event conditions and highlights the fact that flood magnitude is primarily controlled by the meteorological event, whereas flood extent is primarily controlled by the soil moisture conditions.
Describing flood occurrence, processes and characteristics as a function of hydro-meteorological patterns, this thesis is part of a paradigm shift towards a process-based understanding of floods. The results highlight that soil moisture patterns as well as weather patterns are not only beneficial to a probabilistic conception of flood initiation but also provide information on the involved flood processes and the resulting flood characteristics.
Variational methods are employed in situations where exact Bayesian inference becomes intractable due to the difficulty in performing certain integrals. Typically, variational methods postulate a tractable posterior and formulate a lower bound on the desired integral to be approximated, e.g. marginal likelihood. The lower bound is then optimised with respect to its free parameters, the so-called variational parameters. However, this is not always possible as for certain integrals it is very challenging (or tedious) to come up with a suitable lower bound. Here, we propose a simple scheme that overcomes some of the awkward cases where the usual variational treatment becomes difficult. The scheme relies on a rewriting of the lower bound on the model log-likelihood. We demonstrate the proposed scheme on a number of synthetic and real examples, as well as on a real geophysical model for which the standard variational approaches are inapplicable.
We present a calibration-set based on modern pollen and satellite-based Advanced Very High Resolution Radiometer (AVHRR) observations of woody cover (including needleleaved, broadleaved and total tree cover) in eastern continental Asia, which shows good performance under cross-validation with the modern analogue technique (all the coefficients of determination between observed and predicted values are greater than 0.65). The calibration-set is used to reconstruct woody cover from a taxonomically harmonized and temporally standardized fossil pollen dataset (including 274 cores) with 500-year resolution over the last 22 kyr. The spatial range of forest has not noticeably changed in eastern continental Asia during the last 22 kyr, although woody cover has, especially at the margin of the eastern Tibetan Plateau and in the forest-steppe transition area of north-central China. Vegetation was sparse during the LGM in the present forested regions, but woody cover increased markedly at the beginning of the Bolling/Allerod period (B/A; ca. 14.5 ka BP) and again at the beginning of the Holocene (ca. 11.5 ka BP), and is related to the enhanced strength of the East Asian Summer Monsoon. Forest flourished in the mid Holocene (ca. 8 ka BP) possibly due to favourable climatic conditions. In contrast, cover was stable in southern China (high cover) and arid central Asia (very low cover) throughout the investigated period. Forest cover increased in the north-eastern part of China during the Holocene. Comparisons of these regional pollen-based results with simulated forest cover from runs of a global climate model (for 9, 6 and 0 ka BP (ECHAM5/JSBACH similar to 1.125 degrees spatial resolution)) reveal many similarities in temporal change. The Holocene woody cover history of eastern continental Asia is different from that of other regions, likely controlled by different climatic variables, i.e. moisture in eastern continental Asia; temperature in northern Eurasia and North America. (C) 2016 Elsevier Ltd. All rights reserved.
Late Quaternary paleoenvironmental records from the Chatanika River valley near Fairbanks (Alaska)
(2016)
Perennially-frozen deposits are considered as excellent paleoenvironmental archives similar to lacustrine, deep marine, and glacier records because of the long-term and good preservation of fossil records under stable permafrost conditions. A permafrost tunnel in the Vault Creek Valley (Chatanika River Valley, near Fairbanks) exposes a sequence of frozen deposits and ground ice that provides a comprehensive set of proxies to reconstruct the late Quaternary environmental history of Interior Alaska. The multi-proxy approach includes different dating techniques (radiocarbon-accelerator mass spectrometry [AMS C-14], optically stimulated luminescence [OSL], thorium/uranium radioisotope disequilibria [Th-230/U]), as well as methods of sedimentology, paleoecology, hydrochemistry, and stable isotope geochemistry of ground ice. The studied sequence consists of 36-m-thick late Quaternary deposits above schistose bedrock. Main portions of the sequence accumulated during the early and middle Wisconsin periods. The lowermost unit A consists of about 9-m-thick ice-bonded fluvial gravels with sand and peat lenses. A late Sangamon (MIS 5a) age of unit A is assumed. Spruce forest with birch, larch, and some shrubby alder dominated the vegetation. High presence of Sphagnum spores and Cyperaceae pollen points to mires in the Vault Creek Valley. The overlying unit B consists of 10-m-thick alternating fluvial gravels, loess-like silt, and sand layers, penetrated by small ice wedges. OSL dates support a stadial early Wisconsin (MIS 4) age of unit B. Pollen and plant macrofossil data point to spruce forests with some birch interspersed with wetlands around the site. The following unit C is composed of 15-m-thick ice-rich loess-like and organic-rich silt with fossil bones and large ice wedges. Unit C formed during the interstadial mid-Wisconsin (MIS 3) and stadial late Wisconsin (MIS 2) as indicated by radiocarbon ages. Post-depositional slope processes significantly deformed both, ground ice and sediments of unit C. Pollen data show that spruce forests and wetlands dominated the area. The macrofossil remains of Picea, Larix, and Alnus incana ssp. tenuifolia also prove the existence of boreal coniferous forests during the mid-Wisconsin interstadial, which were replaced by treeless tundra-steppe vegetation during the late Wisconsin stadial. Unit C is discordantly overlain by the 2-m-thick late Holocene deposits of unit D. The pollen record of unit D indicates boreal forest vegetation similar to the modern one. The permafrost record from the Vault Creek tunnel reflects more than 90 ka of periglacial landscape dynamics triggered by fluvial and eolian accumulation, and formation of ice-wedge polygons and post depositional deformation by slope processes. The record represents a typical Wisconsin valley-bottom facies in Central Alaska. (C) 2016 Elsevier Ltd. All rights reserved.
Waterbodies such as lakes and ponds are abundant in vast Arctic landscapes and strongly affect the thermal state of the surrounding permafrost. In order to gain a better understanding of the impact of small-and medium-sized waterbodies on permafrost and the formation of thermokarst, a land surface model was developed that can represent the vertical and lateral thermal interactions between waterbodies and permafrost. The model was validated using temperature measurements from two typical waterbodies located within the Lena River delta in northern Siberia. Impact simulations were performed under current climate conditions as well as under a moderate and a strong climate-warming scenario. The performed simulations demonstrate that small waterbodies can rise the sediment surface temperature by more than 10 degrees C and accelerate permafrost thaw by a factor of between 4 and 5. Up to 70% of this additional heat flux into the ground was found to be dissipated into the surrounding permafrost by lateral ground heat flux in the case of small, shallow, and isolated waterbodies. Under moderate climate warming, the lateral heat flux was found to reduce permafrost degradation underneath waterbodies by a factor of 2. Under stronger climatic warming, however, the lateral heat flux was too small to prevent rapid permafrost degradation. The lateral heat flux was also found to strongly impede the formation of thermokarst. Despite this stabilizing effect, our simulations have demonstrated that underneath shallow waterbodies (<1 m), thermokarst initiation happens 30 to 40 years earlier than in simulations without preexisting waterbody.
We report the first high-resolution (20-50 years) mid- to late Holocene pollen records from Lake Teletskoye, the largest lake in the Altai Mountains, in south-eastern West Siberia. Generally, the mid- to late Holocene (the last 4250 years) vegetation of the north-eastern Altai, as recorded in two studied sediment cores, is characterised by Siberian pine-spruce-fir forests that are similar to those of the present day. A relatively cool and dry interval with July temperatures lower than those of today occurred between 3.9 and 3.6 ka BP. The widespread distribution of open, steppe-like communities with Artemisia, Chenopodiaceae and Cyperaceae reflects maximum deforestation during this interval. After ca. 3.5 ka BP, the coniferous mountain taiga spread significantly, with maximum woody coverage and taiga biome scores between ca. 2.7 and 1.6 ka BP. This coincides well with the highest July temperature (approximately 1 degrees C higher than today) intervals. A short period of cooling about 13-1.4 ka BP could have been triggered by the increased volcanic activity recorded across the Northern Hemisphere. A new period of cooling started around 1100-1150 CE, with the minimum July temperatures occurring between 1450 and 1800 CE. (C) 2016 Elsevier B.V. All rights reserved.
Freshwater ostracods (Crustacea, Ostracoda) are valuable biological indicators. In Arctic environments, their habitat conditions are barely known and the abundance and diversity of ostracods is documented only in scattered records with incomplete ecological characterization. To determine the taxonomic range of ostracod assemblages and their habitat conditions in polygon ponds in the Indigirka Lowland, north-east Siberia, we collected more than 100 living ostracod individuals per site with a plankton net (mesh size 65 mm) and an exhaustor system from 27 water bodies and studied them in the context of substrate and hydrochemical data. During the summer of 2011, a single pond site and its ostracod population was selected for special study. This first record of the ostracod fauna in the Indigirka Lowland comprises eight species and three additional taxa. Fabaeformiscandona krochini and F. groenlandica were documented for the first time in continental Siberia. Repeated sampling of a low-centre polygon pond yielded insights into the population dynamics of F. pedata. We identified air temperature and precipitation as the main external drivers of water temperatures, water levels, ion concentrations and water stable isotope composition on diurnal and seasonal scales.
Glacial legacies on interglacial vegetation at the Pliocene-Pleistocene transition in NE Asia
(2016)
Broad-scale climate control of vegetation is widely assumed. Vegetation-climate lags are generally thought to have lasted no more than a few centuries. Here our palaeoecological study challenges this concept over glacial–interglacial timescales. Through multivariate analyses of pollen assemblages from Lake El’gygytgyn, Russian Far East and other data we show that interglacial vegetation during the Plio-Pleistocene transition mainly reflects conditions of the preceding glacial instead of contemporary interglacial climate. Vegetation–climate disequilibrium may persist for several millennia, related to the combined effects of permafrost persistence, distant glacial refugia and fire. In contrast, no effects from the preceding interglacial on glacial vegetation are detected. We propose that disequilibrium was stronger during the Plio-Pleistocene transition than during the Mid-Pliocene Warm Period when, in addition to climate, herbivory was important. By analogy to the past, we suggest today’s widespread larch ecosystem on permafrost is not in climate equilibrium. Vegetation-based reconstructions of interglacial climates used to assess atmospheric CO2–temperature relationships may thus yield misleading simulations of past global climate sensitivity.
QuestionHow important is the effect of micro-relief and vegetation type on the characteristics of vascular plants and bryophytes in low-centred polygons? LocationSiberian Arctic, Russia. MethodsEight low-centred polygons in northern Siberia were surveyed for vegetation along transects running from the rim to the pond via the rim-pond transition of each polygon and across a vegetation type gradient from open forest to tundra. ResultsThe cover of vascular plants and bryophytes displays no significant differences between the rim and rim-pond transition but is significantly lower in the pond section of the polygons. Alpha-diversity of vascular plants decreases strongly from rim to pond, whereas bryophyte diversity in pond plots is significantly distinct from the rim and the rim-pond transition. There is no clear trend in cover for either plant group along the vegetation type transect and only a weak trend in -diversity. However, both gradients are reflected in the compositional turnover. The applied indicator species analysis identified taxa characteristic of certain environmental conditions. Among others, we found vascular plants primarily characteristic of the rim and bryophyte taxa characteristic of each micro-relief level and vegetation type. ConclusionsThe observed gradual pattern in -diversity and composition of polygonal vegetation suggests that micro-relief is the main driver of changes in the vegetation composition, while vegetation type and the related forest cover change are of subordinate importance for polygonal vegetation patterns along the Siberian tree line.
Abrupt climate changes and fluctuations over short time scales are superimposed on long-term climate changes. Understanding rapid climate fluctuations at the decadal time scale over the past millennium will enhance our understanding of patterns of climate variability and aid in forecasting climate changes in the future. In this study, climate changes on the southeastern Tibetan Plateau over the past millennium were determined from a 4.82-m-long sediment core from Basomtso Lake. At the centennial time scale, the Medieval Climate Anomaly (MCA), Little Ice Age (LIA) and Current Warm Period (CWP) are distinct in the Basomtso region. Rapid climate fluctuations inferred from five episodes with higher sediment input and likely warmer conditions, as well as seven episodes with lower sediment input and likely colder conditions, were well preserved in our record. These episodes with higher and lower sediment input are characterized by abrupt climate changes and short time durations. Spectral analysis indicates that the climate variations at the centennial scale on the southeastern Tibetan Plateau are influenced by solar activity during the past millennium.
We compare the ability of various site-condition proxies (SCPs) to reduce the aleatory variability of ground motion prediction equations (GMPEs). Three SCPs (measured V-S30, inferred V-S30, local topographic slope) and two accelerometric databases (RESORCE and NGA-West2) are considered. An artificial neural network (ANN) approach including a random-effect procedure is used to derive GMPEs setting the relationship between peak ground acceleration (PGA), peak ground velocity (PGV), pseudo-spectral acceleration [PSA(T)], and explanatory variables (M-w, R-JB, and V-S30 or Slope). The analysis is performed using both discrete site classes and continuous proxy values. All "non-measured" SCPs exhibit a rather poor performance in reducing aleatory variability, compared to the better performance of measured V-S30. A new, fully data-driven GMPE based on the NGA-West2 is then derived, with an aleatory variability value depending on the quality of the SCP. It proves very consistent with previous GMPEs built on the same data set. Measuring V-S30 allows for benefit from an aleatory variability reduction up to 15%.
A new view of Ecuador's complex geodynamics has been developed in the course of modeling seismic source zones for probabilistic seismic hazard analysis. This study focuses on two aspects of the plates' interaction at a continental scale: (a) age-related differences in rheology between Farallon and Nazca plates—marked by the Grijalva rifted margin and its inland projection—as they subduct underneath central Ecuador, and (b) the rapidly changing convergence obliquity resulting from the convex shape of the South American northwestern continental margin. Both conditions satisfactorily explain several characteristics of the observed seismicity and of the interseismic coupling. Intermediate-depth seismicity reveals a severe flexure in the Farallon slab as it dips and contorts at depth, originating the El Puyo seismic cluster. The two slabs position and geometry below continental Ecuador also correlate with surface expressions observable in the local and regional geology and tectonics. The interseismic coupling is weak and shallow south of the Grijalva rifted margin and increases northward, with a heterogeneous pattern locally associated to the Carnegie ridge subduction. High convergence obliquity is responsible for the North Andean Block northeastward movement along localized fault systems. The Cosanga and Pallatanga fault segments of the North Andean Block-South American boundary concentrate most of the seismic moment release in continental Ecuador. Other inner block faults located along the western border of the inter-Andean Depression also show a high rate of moderate-size earthquake production. Finally, a total of 19 seismic source zones were modeled in accordance with the proposed geodynamic and neotectonic scheme.
Event coincidence analysis for quantifying statistical interrelationships between event time series
(2016)
Studying event time series is a powerful approach for analyzing the dynamics of complex dynamical systems in many fields of science. In this paper, we describe the method of event coincidence analysis to provide a framework for quantifying the strength, directionality and time lag of statistical interrelationships between event series. Event coincidence analysis allows to formulate and test null hypotheses on the origin of the observed interrelationships including tests based on Poisson processes or, more generally, stochastic point processes with a prescribed inter-event time distribution and other higher-order properties. Applying the framework to country-level observational data yields evidence that flood events have acted as triggers of epidemic outbreaks globally since the 1950s. Facing projected future changes in the statistics of climatic extreme events, statistical techniques such as event coincidence analysis will be relevant for investigating the impacts of anthropogenic climate change on human societies and ecosystems worldwide.
On 16 September 2015, the M-W = 8.2 Illapel megathrust earthquake ruptured the Central Chilean margin. Combining inversions of displacement measurements and seismic waveforms with high frequency (HF) teleseismic backprojection, we derive a comprehensive description of the rupture, which also predicts deep ocean tsunami wave heights. We further determine moment tensors and obtain accurate depth estimates for the aftershock sequence. The earthquake nucleated near the coast but then propagated to the north and updip, attaining a peak slip of 5-6 m. In contrast, HF seismic radiation is mostly emitted downdip of the region of intense slip and arrests earlier than the long period rupture, indicating smooth slip along the shallow plate interface in the final phase. A superficially similar earthquake in 1943 with a similar aftershock zone had a much shorter source time function, which matches the duration of HF seismic radiation in the recent event, indicating that the 1943 event lacked the shallow slip.
High-resolution topographic data greatly facilitate the remote identification of geomorphic features, furnishing valuable information concerning surface processes and characterization of reference markers for quantifying tectonic deformation. Marine terraces have been used as long baseline geodetic markers of relative past sea-level positions, reflecting the interplay between vertical crustal movements and sea-level oscillations. Uplift rates may be determined from the terrace age and the elevation of its shoreline angle, a geomorphic feature that can be correlated with past sea-levels positions. A precise definition of the shoreline angle in time and space is essential to obtain reliable uplift rates with coherent spatial correlation. To improve our ability to rapidly assess and map shoreline angles at regional and local scales, we have developed TerraceM, a MATLAB (R) graphical user interface that allows the shoreline angle and its associated error to be estimated using high-resolution topography. TerraceM uses topographic swath profiles oriented orthogonally to the terrace riser. Four functions are included to analyze the swath profiles and extract the shoreline angle, from both staircase sequences of multiple terraces and rough coasts characterized by eroded remnants of emerged terrace surfaces. The former are measured by outlining the paleocliffs and paieo-platforms and finding their intersection by extrapolating linear regressions, whereas the latter are assessed by automatically detecting peaks of sea-stack tops and back-projecting them to the modern sea cliff. In the absence of rigorous absolute age determinations of marine terraces, their geomorphic age may be estimated using previously published diffusion models. Postprocessing functions are included to obtain first-order statistics of shoreline-angle elevations and their spatial distribution. TerraceM has the ability to process series of profiles from several sites in an efficient and structured workflow. Results may be exported in Google Earth and ESRI shapefile formats. The precision and accuracy of the method have been estimated from a case study at Santa Cruz, California, by comparing TerraceM results with published field measurements. The repeatability was evaluated using multiple measurements made by inexperienced users. TerraceM will improve the efficiency and precision of estimating shoreline-angle elevations in wave-cut terraces in both marine and lacustrine environments.
Although the climate development over the Holocene in the Northern Hemisphere is well known, palaeolimnological climate reconstructions reveal spatiotemporal variability in northern Eurasia. Here we present a multi-proxy study from north-eastern Siberia combining sediment geochemistry, and diatom and pollen data from lake-sediment cores covering the last 38,000 cal. years. Our results show major changes in pyrite content and fragilarioid diatom species distributions, indicating prolonged seasonal lake-ice cover between similar to 13,500 and similar to 8900 cal. years BP and possibly during the 8200 cal. years BP cold event. A pollen-based climate reconstruction generated a mean July temperature of 17.8 degrees C during the Holocene Thermal Maximum (HTM) between similar to 8900 and similar to 4500 cal. years BP. Naviculoid diatoms appear in the late Holocene indicating a shortening of the seasonal ice cover that continues today. Our results reveal a strong correlation between the applied terrestrial and aquatic indicators and natural seasonal climate dynamics in the Holocene. Planktonic diatoms show a strong response to changes in the lake ecosystem due to recent climate warming in the Anthropocene. We assess other palaeolimnological studies to infer the spatiotemporal pattern of the HTM and affirm that the timing of its onset, a difference of up to 3000 years from north to south, can be well explained by climatic teleconnections. The westerlies brought cold air to this part of Siberia until the Laurentide ice sheet vanished 7000 years ago. The apparent delayed ending of the HTM in the central Siberian record can be ascribed to the exceedance of ecological thresholds trailing behind increases in winter temperatures and decreases in contrast in insolation between seasons during the mid to late Holocene as well as lacking differentiation between summer and winter trends in paleolimnological reconstructions. (C) 2015 Elsevier Ltd. All rights reserved.
Carbon and nutrient cycling in kettle hole sediments depending on hydrological dynamics: a review
(2016)
Kettle holes as a specific group of isolated, small lentic freshwater systems (LFS) often are (i) hot spots of biogeochemical cycling and (ii) exposed to frequent sediment desiccation and rewetting. Their ecological functioning is greatly determined by immanent carbon and nutrient transformations. The objective of this review is to elucidate effects of a changing hydrological regime (i.e., dry-wet cycles) on carbon and nutrient cycling in kettle hole sediments. Generally, dry-wet cycles have the potential to increase C and N losses as well as P availability. However, their duration and frequency are important controlling factors regarding direction and intensity of biogeochemical and microbiological responses. To evaluate drought impacts on sediment carbon and nutrient cycling in detail requires the context of the LFS hydrological history. For example, frequent drought events induce physiological adaptation of exposed microbial communities and thus flatten metabolic responses, whereas rare events provoke unbalanced, strong microbial responses. Different potential of microbial resilience to drought stress can irretrievably change microbial communities and functional guilds, gearing cascades of functional responses. Hence, dry-wet events can shift the biogeochemical cycling of organic matter and nutrients to a new equilibrium, thus affecting the dynamic balance between carbon burial and mineralization in kettle holes.
Landscapes can be viewed as spatially heterogeneous areas encompassing terrestrial and aquatic domains. To date, most landscape carbon (C) fluxes have been estimated by accounting for terrestrial ecosystems, while aquatic ecosystems have been largely neglected. However, a robust assessment of C fluxes on the landscape scale requires the estimation of fluxes within and between both landscape components. Here, we compiled data from the literature on C fluxes across the air–water interface from various landscape components. We simulated C emissions and uptake for five different scenarios which represent a gradient of increasing spatial heterogeneity within a temperate young moraine landscape: (I) a homogeneous landscape with only cropland and large lakes; (II) separation of the terrestrial domain into cropland and forest; (III) further separation into cropland, forest, and grassland; (IV) additional division of the aquatic area into large lakes and peatlands; and (V) further separation of the aquatic area into large lakes, peatlands, running waters, and small water bodies These simulations suggest that C fluxes at the landscape scale might depend on spatial heterogeneity and landscape diversity, among other factors. When we consider spatial heterogeneity and diversity alone, small inland waters appear to play a pivotal and previously underestimated role in landscape greenhouse gas emissions that may be regarded as C hot spots. Approaches focusing on the landscape scale will also enable improved projections of ecosystems’ responses to perturbations, e.g., due to global change and anthropogenic activities, and evaluations of the specific role individual landscape components play in regional C fluxes. WIREs Water 2016, 3:601–617. doi: 10.1002/wat2.1147
We review recent progress in studying structural properties of liquids using X-ray absorption spectroscopy coupled with the Paris-Edinburgh press at third-generation synchrotron facilities. This experimental method allows for detecting subtle changes in atomic arrangements of melts over a wide pressure-temperature range. It has been also employed to monitor variations of the local coordination environment of diluted species contained in glasses, liquids and crystalline phases as a function of the pressure and temperature. Such information is of great importance for gaining deeper insights into the physico-chemical properties of liquids at extreme condition, including the understanding of such phenomena as liquid-liquid phase transitions, viscosity drops and various transport properties of geological melts. Here, we describe the experimental approach and discuss its potential in structural characterization on selected scientific highlights. Finally, the current ongoing instrumental developments and future scientific opportunities are discussed.
Ice complex deposits are characteristic, ice-rich formations in northern East Siberia and represent an important part in the arctic carbon pool. Recently, these late Quaternary deposits are the objective of numerous investigations typically relying on outcrop and borehole data. Many of these studies can benefit from a 3D structural model of the subsurface for upscaling their observations or for constraining estimations of inventories, such as the local carbon stock. We have addressed this problem of structural imaging by 3D ground-penetrating radar (GPR), which, in permafrost studies, has been primarily used for 2D profiling. We have used a 3D kinematic GPR surveying strategy at a field site located in the New Siberian Archipelago on top of an ice complex. After applying a 3D GPR processing sequence, we were able to trace two horizons at depths below 20 m. Taking available borehole and outcrop data into account, we have interpreted these two features as interfaces of major lithologic units and derived a 3D cryostratigraphic model of the subsurface. Our data example demonstrated that a 3D surveying and processing strategy was crucial at our field site and showed the potential of 3D GPR to image geologic structures in complex ice-rich permafrost landscapes.
Analysis of time-lapse ground-penetrating radar (GPR) data can provide information regarding subsurface hydrological processes, such as preferential flow. However, the analysis of time-lapse data is often limited by data quality; for example, for noisy input data, the interpretation of difference images is often difficult. Motivated by modern image-processing tools, we have developed two robust GPR attributes, which allow us to distinguish amplitude (contrast similarity) and time-shift (structural similarity) variations related to differences between individual time-lapse GPR data sets. We tested and evaluated our attributes using synthetic data of different complexity. Afterward, we applied them to a field data example, in which subsurface flow was induced by an artificial rainfall event. For all examples, we identified our structural similarity attribute to be a robust measure for highlighting time-lapse changes also in data with low signal-to-noise ratios. We determined that our new attribute-based workflow is a promising tool to analyze time-lapse GPR data, especially for imaging subsurface hydrological processes.
We present an algorithm that performs sequentially one-dimensional inversion of subsurface magnetic permeability and electrical conductivity by using multi-configuration electromagnetic induction sensor data. The presented method is based on the conversion of the in-phase and out-of-phase data into effective magnetic permeability and electrical conductivity of the equivalent homogeneous half-space. In the case of small-offset systems, such as portable electromagnetic induction sensors, for which in-phase and out-of-phase data are moderately coupled, the effective half-space magnetic permeability and electrical conductivity can be inverted sequentially within an iterative scheme. We test and evaluate the proposed inversion strategy using synthetic and field examples. First, we apply it to synthetic data for some highly magnetic environments. Then, the method is tested on real field data acquired in a basaltic environment to image a formation of archaeological interest. These examples demonstrate that a joint interpretation of in-phase and out-of-phase data leads to a better characterisation of the subsurface in magnetic environments such as volcanic areas.
Fossil oyster shells are well-suited to provide palaeotemperature proxies from geologic to seasonal timescales due to their ubiquitous occurrence from Triassic to Quaternary sediments, the seasonal nature of their shell growth and their relative strong resistance to post-mortem alteration. However, the common use to translate calcitic oxygen isotopes into palaeotemperatures is challenged by uncertainties in accounting for past seawater delta O-18, especially in shallow coastal environment where oysters calcify. In principle, the Mg/Ca ratio in oyster shells can provide an alternative palaeothermometer. Several studies provided temperature calibrations for this potential proxy based on modem species, nevertheless their application to palaeo-studies remains hitherto unexplored. Here, we show that past temperature variability in seawater can be obtained from Mg/Ca analyses from selected fossil oyster species and specimens. High-resolution Mg/Ca profiles, combined with delta O-18, were obtained along 41 fossil oyster shells of seven different species from the Palaeogene Proto-Paratethys sea (Central Asia) found in similar as well as different depositional age and environments providing comparison. Suitable Mg/Ca profiles, defined by continuous cyclicity and reproducibility within one shell, are found to be consistent for specimens of the same species but differ systematically between species, implying a dominant species-specific effect on the Mg/Ca signal. Two species studied here (Ostrea (Turkostrea) strictiplicata and Sokolowia buhsii) provide an excellent proxy for palaeoclimate reconstruction from China to Europe in Palaeogene marine sediments. More generally, the protocol developed here can be applied to identify other fossil oyster species suitable for palaeoclimate reconstructions. (C) 2015 Elsevier B.V. All rights reserved.
The loess and paleosol sequences of the Chinese Loess Plateau are composed of Quaternary dust, the origin of which has been the subject of considerable debate. Some recent U-Pb geochronological studies of eolian zircons have proposed the existence of two major wind pathways: from the north and northwest, through the Badan Jaran, Tengger, and Mu Us Deserts during interglacials, and from the west, through the Qaidam Basin during glacials. Others have emphasized the importance of Yellow River sediment supply in the Chinese Loess Plateau sediment budget. However, tracking dust source regions through U-Pb dating in a statistically robust manner is particularly complex given the similar age peaks in the age probability distributions of potential source regions in Central Asia. This paper presents 2410 new U-Pb ages of detrital zircons from wind-eroded strata, Quaternary eolian deposits, and modern river sands in central China in order to increase the robustness and the spatial resolution of zircon age distributions in dust source regions. We then propose a new mixture modeling technique to statistically address the contribution of these different sources to the Chinese Loess Plateau sedimentary budget. Our contribution estimates indicate that eolian supply to the Chinese Loess Plateau is dominated (60%-70%) by reworking of Yellow River sediment. Moreover, evidence of Qaidam Basin-sourced zircons (15%-20%) in both loess (glacial) and paleosol (interglacial) layers corroborates the existence of an erosive wind pathway through the Qaidam Basin during glacials and implies that a substantial portion of the interglacial dust is recycled from older glacial loess. We propose that sediment reworking of Yellow River sediment and older loess deposits by wind on the Chinese Loess Plateau homogenized the eolian zircon populations toward a glacial provenance due to higher (2-20 times) dust accumulation rates during glacials. These findings suggest that the Chinese Loess Plateau has evolved as a more dynamic landform than previous thought, where wind deflation, fluvial input, lateral transport, and accumulation of sediment are equally important. These internal reworking effects would then significantly bias the paleoclimatic interpretations based on eolian dust properties of the Chinese Loess Plateau.
The Pamirs represent the indented westward continuation of the northern margin of the Tibetan Plateau, dividing the Tarim and Tajik basins. Their evolution may be a key factor influencing aridification of the Asian interior, yet the tectonics of the Pamir Salient are poorly understood. We present a provenance study of the Aertashi section, a Paleogene to late Neogene clastic succession deposited in the Tarim basin to the north of the NW margin of Tibet (the West Kunlun) and to the east of the Pamirs. Our detrital zircon U-Pb ages coupled with zircon fission track, bulk rock Sm-Nd, and petrography data document changes in contributing source terranes during the Oligocene to Miocene, which can be correlated to regional tectonics. We propose a model for the evolution of the Pamir and West Kunlun (WKL), in which the WKL formed topography since at least similar to 200 Ma. By similar to 25 Ma, movement along the Pamir-bounding faults such as the Kashgar-Yecheng Transfer System had commenced, marking the onset of Pamir indentation into the Tarim-Tajik basin. This is coincident with basinward expansion of the northern WKL margin, which changed the palaeodrainage pattern within the Kunlun, progressively cutting off the more southerly WKL sources from the Tarim basin. An abrupt change in the provenance and facies of sediments at Aertashi has a maximum age of 14 Ma; this change records when the Pamir indenter had propagated sufficiently far north that the North Pamir was now located proximal to the Aertashi region.
Steppe vegetation represents a key marker of past Asian aridification and is associated with monsoonal intensification. Little is, however, known about the origin of this pre-Oligocene vegetation, its specific composition and how it changed over time and responded to climatic variations. Here, we describe the morphological characters of Ephedraceae pollen in Eocene strata of the Xining Basin and compare the pollen composition with the palynological composition of Late Cretaceous and Paleocene deposits of the Xining Basin and the Quaternary deposits of the Qaidam Basin. We find that the Late Cretaceous steppe was dominated by Gnetaceaepollenites; in the transition from the Cretaceous to the Paleocene, Gnetaceaepollenites became extinct and Ephedripites subgenus Ephedripites dominated the flora with rare occurrences of Ephedripites subgen. Distachyapites; the middle to late Eocene presents a strong increase of Ephedripites subgen. Distachyapites; and the Quaternary/Recent is marked by a significantly lower diversity of Ephedraceae (and Nitrariaceae) compared to the Eocene. In the modern landscape of China, only a fraction of the Paleogene species diversity of Ephedraceae remains and we propose that these alterations in Ephedreaceae composition occurred in response to the climatic changes at least since the Eocene. In particular, the strong Eocene monsoons that enhanced the continental aridification may have played an important role in the evolution of Ephedripites subgen. Distachyapites triggering an evolutionary shift to wind-pollination in this group. Conceivably, the Ephedraceae/Nitrariaceae dominated steppe ended during the Eocene/Oligocene climatic cooling and aridification, which favoured other plant taxa.
The onset of modern central Asian atmospheric circulation is traditionally linked to the interplay of surface uplift of the Mongolian and Tibetan-Himalayan orogens, retreat of the Paratethys sea from central Asia and Cenozoic global cooling. Although the role of these players has not yet been unravelled, the vast dust deposits of central China support the presence of arid conditions and modern atmospheric pathways for the last 25 million years (Myr). Here, we present provenance data from older (42-33 Myr) dust deposits, at a time when the Tibetan Plateau was less developed, the Paratethys sea still present in central Asia and atmospheric pCO(2) much higher. Our results show that dust sources and near-surface atmospheric circulation have changed little since at least 42 Myr. Our findings indicate that the locus of central Asian high pressures and concurrent aridity is a resilient feature only modulated by mountain building, global cooling and sea retreat.
Vegetation has long been hypothesized to influence the nature and rates of surface processes. We test the possible impact of vegetation and climate on denudation rates at orogen scale by taking advantage of a pronounced along-strike gradient in rainfall and vegetation density in the Himalaya. We combine 12 new Be-10 denudation rates from the Sutlej Valley and 123 published denudation rates from fluvially-dominated catchments in the Himalaya with remotely-sensed measures of vegetation density and rainfall metrics, and with tectonic and lithologic constraints. In addition, we perform topographic analyses to assess the contribution of vegetation and climate in modulating denudation rates along strike. We observe variations in denudation rates and the relationship between denudation and topography along strike that are most strongly controlled by local rainfall amount and vegetation density, and cannot be explained by along-strike differences in tectonics or lithology. A W-E along-strike decrease in denudation rate variability positively correlates with the seasonality of vegetation density (R = 0.95, p < 0.05), and negatively correlates with mean vegetation density (R = -0.84, p < 0.05). Vegetation density modulates the topographic response to changing denudation rates, such that the functional relationship between denudation rate and topographic steepness becomes increasingly linear as vegetation density increases. We suggest that while tectonic processes locally control the pattern of denudation rates across strike of the Himalaya (i.e., S-N), along strike of the orogen (i.e., E-W) climate exerts a measurable influence on how denudation rates scatter around long-term, tectonically-controlled erosion, and on the functional relationship between topography and denudation. (C) 2016 Elsevier B.V. All rights reserved.
The role that climate and environmental history may have played in influencing human evolution has been the focus of considerable interest and controversy among paleoanthropologists for decades. Prior attempts to understand the environmental history side of this equation have centered around the study of outcrop sediments and fossils adjacent to where fossil hominins (ancestors or close relatives of modern humans) are found, or from the study of deep sea drill cores. However, outcrop sediments are often highly weathered and thus are unsuitable for some types of paleoclimatic records, and deep sea core records come from long distances away from the actual fossil and stone tool remains. The Hominin Sites and Paleolakes Drilling Project (HSPDP) was developed to address these issues. The project has focused its efforts on the eastern African Rift Valley, where much of the evidence for early hominins has been recovered. We have collected about 2 km of sediment drill core from six basins in Kenya and Ethiopia, in lake deposits immediately adjacent to important fossil hominin and archaeological sites. Collectively these cores cover in time many of the key transitions and critical intervals in human evolutionary history over the last 4 Ma, such as the earliest stone tools, the origin of our own genus Homo, and the earliest anatomically modern Homo sapiens. Here we document the initial field, physical property, and core description results of the 2012-2014 HSPDP coring campaign.