Refine
Year of publication
Document Type
- Article (18)
- Doctoral Thesis (9)
- Postprint (1)
Keywords
- erosion (28) (remove)
Institute
- Institut für Geowissenschaften (28) (remove)
Rainfall, snow-, and glacial melt throughout the Himalaya control river discharge, which is vital for maintaining agriculture, drinking water and hydropower generation. However, the spatiotemporal contribution of these discharge components to Himalayan rivers is not well understood, mainly because of the scarcity of ground-based observations. Consequently, there is also little known about the triggers and sources of peak sediment flux events, which account for extensive hydropower reservoir filling and turbine abrasion. We therefore lack basic information on the distribution of water resources and controls of erosion processes. In this thesis, I employ various methods to assess and quantify general characteristics of and links between precipitation, river discharge, and sediment flux in the Sutlej Valley. First, I analyze daily precipitation data (1998-2007) from 80 weather stations in the western Himalaya, to decipher the distribution of rain- and snowfall. Rainfall magnitude frequency analyses indicate that 40% of the summer rainfall budget is attributed to monsoonal rainstorms, which show higher variability in the orogenic interior than in frontal regions. Combined analysis of rainstorms and sediment flux data of a major Sutlej River tributary indicate that monsoonal rainfall has a first order control on erosion processes in the orogenic interior, despite the dominance of snowfall in this region. Second, I examine the contribution of rainfall, snow and glacial melt to river discharge in the Sutlej Valley (s55,000 km2), based on a distributed hydrological model, which covers the period 2000-2008. To achieve high spatial and daily resolution despite limited ground-based observations the hydrological model is forced by daily remote sensing data, which I adjusted and calibrated with ground station data. The calibration shows that the Tropical Rainfall Measuring Mission (TRMM) 3B42 rainfall product systematically overestimates rainfall in semi-arid and arid regions, increasing with aridity. The model results indicate that snowmelt-derived discharge (74%) is most important during the pre-monsoon season (April to June) whereas rainfall (56%) and glacial melt (17%) dominate the monsoon season (July-September). Therefore, climate change most likely causes a reduction in river discharge during the pre-monsoon season, which especially affects the orogenic interior. Third, I investigate the controls on suspended sediment flux in different parts of the Sutlej catchments, based on daily gauging data from the past decade. In conjunction with meteorological data, earthquake records, and rock strength measurements I find that rainstorms are the most frequent trigger of high-discharge events with peaks in suspended sediment concentrations (SSC) that account for the bulk of the suspended sediment flux. The suspended sediment flux increases downstream, mainly due to increases in runoff. Pronounced erosion along the Himalayan Front occurs throughout the monsoon season, whereas efficient erosion of the orogenic interior is confined to single extreme events. The results of this thesis highlight the importance of snow and glacially derived melt waters in the western Himalaya, where extensive regions receive only limited amounts of monsoonal rainfall. These regions are therefore particularly susceptible to global warming with major implications on the hydrological cycle. However, the sediment discharge data show that infrequent monsoonal rainstorms that pass the orographic barrier of the Higher Himalaya are still the primary trigger of the highest-impact erosion events, despite being subordinate to snow and glacially–derived discharge. These findings may help to predict peak sediment flux events and could underpin the strategic development of preventative measures for hydropower infrastructures.
Timber harvesting by clear cutting is known to impose environmental impacts, including severe disturbance of the soil hydraulic properties which intensify the frequency and magnitude of surface runoff and soil erosion. However, it remains unanswered if harvest areas act as sources or sinks for runoff and soil erosion and whether such behavior operates in a steady state or evolves through time. For this purpose, 92 small-scale rainfall simulations of different intensities were carried out under pine plantation conditions and on two clear-cut harvest areas of different age. Nonparametrical Random Forest statistical models were set up to quantify the impact of environmental variables on the hydrological and erosion response. Regardless of the applied rainfall intensity, runoff always initiated first and yielded most under plantation cover. Counter to expectations, infiltration rates increased after logging activities. Once a threshold rainfall intensity of 20mm/h was exceeded, the younger harvest area started to act as a source for both runoff and erosion after connectivity was established, whereas it remained a sink under lower applied rainfall intensities. The results suggest that the impact of microtopography on surface runoff connectivity and water-repellent properties of the topsoil act as first-order controls for the hydrological and erosion processes in such environments. Fast rainfall-runoff response, sediment-discharge-hystereses, and enhanced postlogging groundwater recharge at catchment scale support our interpretation. At the end, we show the need to account for nonstationary hydrological and erosional behavior of harvest areas, a fact previously unappreciated in predictive models.
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.
Reconstructing rates and patterns of colluvial soil redistribution in agrarian (hummocky) landscapes
(2019)
Humans have triggered or accelerated erosion processes since prehistoric times through agricultural practices. Optically stimulated luminescence (OSL) is widely used to quantify phases and rates of the corresponding landscape change, by measuring the last moment of daylight exposure of sediments. However, natural and anthropogenic mixing processes, such as bioturbation and tillage, complicate the use of OSL as grains of different depositional ages become mixed, and grains become exposed to light even long after the depositional event of interest. Instead, OSL determines the stabilization age, indicating when sediments were buried below the active mixing zone. These stabilization ages can cause systematic underestimation when calculating deposition rates. Our focus is on colluvial deposition in a kettle hole in the Uckermark region, northeastern Germany. We took 32 samples from five locations in the colluvium filling the kettle hole to study both spatial and temporal patterns in colluviation. We combined OSL dating with advanced age modelling to determine the stabilization age of colluvial sediments. These ages were combined with an archaeological reconstruction of historical ploughing depths to derive the levels of the soil surface at the moment of stabilization; the deposition depths, which were then used to calculate unbiased deposition rates. We identified two phases of colluvial deposition. The oldest deposits (similar to 5 ka) were located at the fringe of the kettle hole and accumulated relatively slowly, whereas the youngest deposits (<0.3 ka) rapidly filled the central kettle hole with rates of two orders of magnitude higher. We suggest that the latter phase is related to artificial drainage, facilitating accessibility in the central depression for agricultural practices. Our results show the need for numerical dating techniques that take archaeological and soil-geomorphological information into account to identify spatiotemporal patterns of landscape change, and to correctly interpret landscape dynamics in anthropogenically influenced hilly landscapes. (c) 2019 The Authors. Earth Surface Processes and Landforms Published by John Wiley & Sons Ltd.
The global carbon cycle is closely linked to Earth’s climate. In the context of continuously unchecked anthropogenic CO₂ emissions, the importance of natural CO₂ bond and carbon storage is increasing. An important biogenic mechanism of natural atmospheric CO₂ drawdown is the photosynthetic carbon fixation in plants and the subsequent longterm deposition of plant detritus in sediments.
The main objective of this thesis is to identify factors that control mobilization and transport of plant organic matter (pOM) through rivers towards sedimentation basins. I investigated this aspect in the eastern Nepalese Arun Valley. The trans-Himalayan Arun River is characterized by a strong elevation gradient (205 − 8848 m asl) that is accompanied by strong changes in ecology and climate ranging from wet tropical conditions in the Himalayan forelad to high alpine tundra on the Tibetan Plateau. Therefore, the Arun is an excellent natural laboratory, allowing the investigation of the effect of vegetation cover, climate, and topography on plant organic matter mobilization and export in tributaries along the gradient.
Based on hydrogen isotope measurements of plant waxes sampled along the Arun River and its tributaries, I first developed a model that allows for an indirect quantification of pOM contributed to the mainsetm by the Arun’s tributaries. In order to determine the role of climatic and topographic parameters of sampled tributary catchments, I looked for significant statistical relations between the amount of tributary pOM export and tributary characteristics (e.g. catchment size, plant cover, annual precipitation or runoff, topographic measures). On one hand, I demonstrated that pOMsourced from the Arun is not uniformly derived from its entire catchment area. On the other, I showed that dense vegetation is a necessary, but not sufficient, criterion for high tributary pOM export. Instead, I identified erosion and rainfall and runoff as key factors controlling pOM sourcing in the Arun Valley. This finding is supported by terrestrial cosmogenic nuclide concentrations measured on river sands along the Arun and its tributaries in order to quantify catchment wide denudation rates. Highest denudation rates corresponded well with maximum pOM mobilization and export also suggesting the link between erosion and pOM sourcing.
The second part of this thesis focusses on the applicability of stable isotope records such as plant wax n-alkanes in sediment archives as qualitative and quantitative proxy for the variability of past Indian Summer Monsoon (ISM) strength. First, I determined how ISM strength affects the hydrogen and oxygen stable isotopic composition (reported as δD and δ18O values vs. Vienna Standard Mean Ocean Water) of precipitation in the Arun Valley and if this amount effect (Dansgaard, 1964) is strong enough to be recorded in potential paleo-ISM isotope proxies. Second, I investigated if potential isotope records across the Arun catchment reflect ISM strength dependent precipitation δD values only, or if the ISM isotope signal is superimposed by winter precipitation or glacial melt. Furthermore, I tested if δD values of plant waxes in fluvial deposits reflect δD values of environmental waters in the respective catchments.
I showed that surface water δD values in the Arun Valley and precipitation δD from south of the Himalaya both changed similarly during two consecutive years (2011 & 2012) with distinct ISM rainfall amounts (~20% less in 2012). In order to evaluate the effect of other water sources (Winter-Westerly precipitation, glacial melt) and evapotranspiration in the Arun Valley, I analysed satellite remote sensing data of rainfall distribution (TRMM 3B42V7), snow cover (MODIS MOD10C1), glacial coverage (GLIMSdatabase, Global Land Ice Measurements from Space), and evapotranspiration (MODIS MOD16A2). In addition to the predominant ISM in the entire catchment I found through stable isotope analysis of surface waters indications for a considerable amount of glacial melt derived from high altitude tributaries and the Tibetan Plateau. Remotely sensed snow cover data revealed that the upper portion of the Arun also receives considerable winter precipitation, but the effect of snow melt on the Arun Valley hydrology could not be evaluated as it takes place in early summer, several months prior to our sampling campaigns. However, I infer that plant wax records and other potential stable isotope proxy archives below the snowline are well-suited for qualitative, and potentially quantitative, reconstructions of past changes of ISM strength.
Sediments in rivers record the dynamics of erosion processes. While bulk sediment fluxes are easily and routinely obtained, sediment caliber remains underexplored when inferring erosion mechanisms. Yet sediment grain size distributions may be the key to discriminating their origin. We have studied grain size-specific suspended sediment fluxes in the Kali Gandaki, a major trans-Himalayan river. Two strategically located gauging stations enable tracing of sediment caliber on either side of the Himalayan orographic barrier. The data show that fine sediment input into the northern headwaters is persistent, while coarse sediment comes from the High Himalayas during the summer monsoon. A temporally matching landslide inventory similarly indicates the prominence of monsoon-driven hillslope mass wasting. Thus, mechanisms of sediment supply can leave strong traces in the fluvial caliber, which could project well beyond the mountain front and add to the variability of the sedimentary record of orogen erosion.
This study presents the development of 1D and 2D Surface Evolution Codes (SECs) and their coupling to any lithospheric-scale (thermo-)mechanical code with a quadrilateral structured surface mesh.
Both SECs involve diffusion as approach for hillslope processes and the stream power law to reflect riverbed incision. The 1D SEC settles sediment that was produced by fluvial incision in the appropriate minimum, while the supply-limited 2D SEC DANSER uses a fast filling algorithm to model sedimantation. It is based on a cellular automaton. A slope-dependent factor in the sediment flux extends the diffusion equation to nonlinear diffusion. The discharge accumulation is achieved with the D8-algorithm and an improved drainage accumulation routine. Lateral incision enhances the incision's modelling. Following empirical laws, it incises channels of several cells width.
The coupling method enables different temporal and spatial resolutions of the SEC and the thermo-mechanical code. It transfers vertical as well as horizontal displacements to the surface model. A weighted smoothing of the 3D surface displacements is implemented. The smoothed displacement vectors transmit the deformation by bilinear interpolation to the surface model. These interpolation methods ensure mass conservation in both directions and prevent the two surfaces from drifting apart.
The presented applications refer to the evolution of the Pamir orogen. A calibration of DANSER's parameters with geomorphological data and a DEM as initial topography highlights the advantage of lateral incision. Preserving the channel width and reflecting incision peaks in narrow channels, this closes the huge gap between current orogen-scale incision models and observed topographies.
River capturing models in a system of fault-bounded block rotations reaffirm the importance of the lateral incision routine for capturing events with channel initiation. The models show a low probability of river capturings with large deflection angles. While the probability of river capturing is directly depending on the uplift rate, the erodibility inside of a dip-slip fault speeds up headward erosion along the fault: The model's capturing speed increases within a fault.
Coupling DANSER with the thermo-mechanical code SLIM 3D emphasizes the versatility of the SEC. While DANSER has minor influence on the lithospheric evolution of an indenter model, the brittle surface deformation is strongly affected by its sedimentation, widening a basin in between two forming orogens and also the southern part of the southern orogen to south, east and west.
In the last decade growing interest has emerged in quantifying the spatial and temporal variations in mountain building. Until recently, insufficient data have been available to attempt such a task at the scale of large orogens such as the Himalaya. The Himalaya accommodates ongoing convergence between India and Eurasia and is a focal point for studying orogen evolution and hypothesized interactions between tectonics and climate. Here we integrate 1126 published bedrock mineral cooling ages with a transient 1D Monte-Carlo thermal-kinematic erosion model to quantify the denudation histories along similar to 2700 km of the Himalaya. The model free parameter is a temporally variable denudation rate from 50 Ma to present. Thermophysical material properties and boundary conditions were tuned to individual study areas. Monte-Carlo simulations were conducted to identify the range of denudation histories that can reproduce the observed cooling ages. Results indicate large temporal and spatial variations in denudation and these are resolvable across different tectonic units of the Himalaya. More specifically, across > 1000 km of the southern Greater Himalaya denudation rates were highest (similar to 1.5-3 mm/yr) between similar to 10 and 2 Ma and lower (0.5-2.6 mm/yr) over the last 2 My. These differences are best determined in the NW-Himalaya. In contrast to this, across the similar to 2500 km length of the northern Greater Himalaya denudation rates vary over length scales of similar to 300-1700 km. Slower denudation (<1 mm/yr) occurred between 10 and 4 Ma followed by a large increase (1.2-2.6 mm/yr) in the last similar to 4 Ma. We find that only the southern Greater Himalayan Sequence clearly supports a continuous co-evolution of tectonics, climate and denudation. Results from the higher elevation northern Greater Himalaya suggest either tectonic driven variations in denudation due to a ramp-flat geometry in the main decollement and/or recent glacially enhanced denudation.
Logging and large earthquakes are disturbances that may significantly affect hydrological and erosional processes and process rates, although in decisively different ways. Despite numerous studies that have documented the impacts of both deforestation and earthquakes on water and sediment fluxes, a number of details regarding the timing and type of de- and reforestation; seismic impacts on subsurface water fluxes; or the overall geomorphic work involved have remained unresolved. The main objective of this thesis is to address these shortcomings and to better understand and compare the hydrological and erosional process responses to such natural and man-made disturbances. To this end, south-central Chile provides an excellent natural laboratory owing to its high seismicity and the ongoing conversion of land into highly productive plantation forests. In this dissertation I combine paired catchment experiments, data analysis techniques, and physics-based modelling to investigate: 1) the effect of plantation forests on water resources, 2) the source and sink behavior of timber harvest areas in terms of overland flow generation and sediment fluxes, 3) geomorphic work and its efficiency as a function of seasonal logging, 4) possible hydrologic responses of the saturated zone to the 2010 Maule earthquake and 5) responses of the vadose zone to this earthquake. Re 1) In order to quantify the hydrologic impact of plantation forests, it is fundamental to first establish their water balances. I show that tree species is not significant in this regard, i.e. Pinus radiata and Eucalyptus globulus do not trigger any decisive different hydrologic response. Instead, water consumption is more sensitive to soil-water supply for the local hydro-climatic conditions. Re 2) Contradictory opinions exist about whether timber harvest areas (THA) generate or capture overland flow and sediment. Although THAs contribute significantly to hydrology and sediment transport because of their spatial extent, little is known about the hydrological and erosional processes occurring on them. I show that THAs may act as both sources and sinks for overland flow, which in turn intensifies surface erosion. Above a rainfall intensity of ~20 mm/h, which corresponds to <10% of all rainfall, THAs may generate runoff whereas below that threshold they remain sinks. The overall contribution of Hortonian runoff is thus secondary considering the local rainfall regime. The bulk of both runoff and sediment is generated by Dunne, saturation excess, overland flow. I also show that logging may increase infiltrability on THAs which may cause an initial decrease in streamflow followed by an increase after the groundwater storage has been refilled. Re 3) I present changes in frequency-magnitude distributions following seasonal logging by applying Quantile Regression Forests at hitherto unprecedented detail. It is clearly the season that controls the hydro-geomorphic work efficiency of clear cutting. Logging, particularly dry seasonal logging, caused a shift of work efficiency towards less flashy and mere but more frequent moderate rainfall-runoff events. The sediment transport is dominated by Dunne overland flow which is consistent with physics-based modelling using WASA-SED. Re 4) It is well accepted that earthquakes may affect hydrological processes in the saturated zone. Assuming such flow conditions, consolidation of saturated saprolitic material is one possible response. Consolidation raises the hydraulic gradients which may explain the observed increase in discharge following earthquakes. By doing so, squeezed water saturates the soil which in turn increases the water accessible for plant transpiration. Post-seismic enhanced transpiration is reflected in the intensification of diurnal cycling. Re 5) Assuming unsaturated conditions, I present the first evidence that the vadose zone may also respond to seismic waves by releasing pore water which in turn feeds groundwater reservoirs. By doing so, water tables along the valley bottoms are elevated thus providing additional water resources to the riparian vegetation. By inverse modelling, the transient increase in transpiration is found to be 30-60%. Based on the data available, both hypotheses, are not testable. Finally, when comparing the hydrological and erosional effects of the Maule earthquake with the impact of planting exotic plantation forests, the overall observed earthquake effects are comparably small, and limited to short time scales.
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.