Refine
Year of publication
Document Type
- Article (9)
- Doctoral Thesis (4)
- Postprint (1)
- Review (1)
Language
- English (15)
Is part of the Bibliography
- yes (15)
Keywords
- landscape evolution (15) (remove)
Continental rifts are excellent regions where the interplay between extension, the build-up of topography, erosion and sedimentation can be evaluated in the context of landscape evolution. Rift basins also constitute important archives that potentially record the evolution and migration of species and the change of sedimentary conditions as a result of climatic change. Finally, rifts have increasingly become targets of resource exploration, such as hydrocarbons or geothermal systems. The study of extensional processes and the factors that further modify the mainly climate-driven surface process regime helps to identify changes in past and present tectonic and geomorphic processes that are ultimately recorded in rift landscapes.
The Cenozoic East African Rift System (EARS) is an exemplary continental rift system and ideal natural laboratory to observe such interactions. The eastern and western branches of the EARS constitute first-order tectonic and topographic features in East Africa, which exert a profound influence on the evolution of topography, the distribution and amount of rainfall, and thus the efficiency of surface processes. The Kenya Rift is an integral part of the eastern branch of the EARS and is characterized by high-relief rift escarpments bounded by normal faults, gently tilted rift shoulders, and volcanic centers along the rift axis.
Considering the Cenozoic tectonic processes in the Kenya Rift, the tectonically controlled cooling history of rift shoulders, the subsidence history of rift basins, and the sedimentation along and across the rift, may help to elucidate the morphotectonic evolution of this extensional province. While tectonic forcing of surface processes may play a minor role in the low-strain rift on centennial to millennial timescales, it may be hypothesized that erosion and sedimentation processes impacted by climate shifts associated with pronounced changes in the availability in moisture may have left important imprints in the landscape.
In this thesis I combined thermochronological, geomorphic field observations, and morphometry of digital elevation models to reconstruct exhumation processes and erosion rates, as well as the effects of climate on the erosion processes in different sectors of the rift. I present three sets of results: (1) new thermochronological data from the northern and central parts of the rift to quantitatively constrain the Tertiary exhumation and thermal evolution of the Kenya Rift. (2) 10Be-derived catchment-wide mean denudation rates from the northern, central and southern rift that characterize erosional processes on millennial to present-day timescales; and (3) paleo-denudation rates in the northern rift to constrain climatically controlled shifts in paleoenvironmental conditions during the early Holocene (African Humid Period).
Taken together, my studies show that time-temperature histories derived from apatite fission track (AFT) analysis, zircon (U-Th)/He dating, and thermal modeling bracket the onset of rifting in the Kenya Rift between 65-50 Ma and about 15 Ma to the present. These two episodes are marked by rapid exhumation and, uplift of the rift shoulders. Between 45 and 15 Ma the margins of the rift experienced very slow erosion/exhumation, with the accommodation of sediments in the rift basin.
In addition, I determined that present-day denudation rates in sparsely vegetated parts of the Kenya Rift amount to 0.13 mm/yr, whereas denudation rates in humid and more densely vegetated sectors of the rift flanks reach a maximum of 0.08 mm/yr, despite steeper hillslopes. I inferred that hillslope gradient and vegetation cover control most of the variation in denudation rates across the Kenya Rift today. Importantly, my results support the notion that vegetation cover plays a fundamental role in determining the voracity of erosion of hillslopes through its stabilizing effects on the land surface.
Finally, in a pilot study I highlighted how paleo-denudation rates in climatic threshold areas changed significantly during times of transient hydrologic conditions and involved a sixfold increase in erosion rates during increased humidity. This assessment is based on cosmogenic nuclide (10Be) dating of quartzitic deltaic sands that were deposited in the northern Kenya Rift during a highstand of Lake Suguta, which was associated with the Holocene African Humid Period. Taken together, my new results document the role of climate variability in erosion processes that impact climatic threshold environments, which may provide a template for potential future impacts of climate-driven changes in surface processes in the course of Global Change.
The tectonically driven growth of mountains reflects the characteristics of the underlying fault systems and the applied tectonic forces. Over time, fault networks might be relatively static, but stress conditions could change and result in variations in fault slip orientation. Such a tectonic landscape would transition from a simple to a composite state: the topography of simple landscapes is correlated with a single set of tectonic boundary conditions, while composite landscapes contain inherited topography due to earlier deformation under different boundary conditions. We use fault interaction modeling to compare vertical displacement fields with topographic metrics to differentiate the two types of landscapes. By successively rotating the axis of maximum horizontal stress, we produce a suite of vertical displacement fields for comparison with real landscapes. We apply this model to a transpressional duplex in the south central Alborz Mountains of Iran, where NW oriented compression was superseded by neotectonic NE compression. The consistency between the modeled displacement field and real landforms indicates that the duplex topography is mostly compatible with the modern boundary conditions, but might include a small remnant from the earlier deformation phase. Our approach is applicable for various tectonic settings and represents an approach to identify the changing boundary conditions that produce composite landscapes. It may be particularly useful for identifying changes that occurred in regions where river profiles may no longer record a signal of the change or where the spatial pattern of uplift is complex.
Drainage divide networks
(2020)
Drainage divides are organized into tree-like networks that may record information about drainage divide mobility. However, views diverge about how to best assess divide mobility. Here, we apply a new approach of automatically extracting and ordering drainage divide networks from digital elevation models to results from landscape evolution model experiments. We compared landscapes perturbed by strike-slip faulting and spatiotemporal variations in erodibility to a reference model to assess which topographic metrics (hillslope relief, flow distance, and chi) are diagnostic of divide mobility. Results show that divide segments that are a minimum distance of similar to 5 km from river confluences strive to attain constant values of hillslope relief and flow distance to the nearest stream. Disruptions of such patterns can be related to mobile divides that are lower than stable divides, closer to streams, and often asymmetric in shape. In general, we observe that drainage divides high up in the network, i.e., at great distances from river confluences, are more susceptible to disruptions than divides closer to these confluences and are thus more likely to record disturbance for a longer time period. We found that across-divide differences in hillslope relief proved more useful for assessing divide migration than other tested metrics. However, even stable drainage divide networks exhibit across-divide differences in any of the studied topographic metrics. Finally, we propose a new metric to quantify the connectivity of divide junctions.
Large earthquakes initiate chains of surface processes that last much longer than the brief moments of strong shaking. Most moderate‐ and large‐magnitude earthquakes trigger landslides, ranging from small failures in the soil cover to massive, devastating rock avalanches. Some landslides dam rivers and impound lakes, which can collapse days to centuries later, and flood mountain valleys for hundreds of kilometers downstream. Landslide deposits on slopes can remobilize during heavy rainfall and evolve into debris flows. Cracks and fractures can form and widen on mountain crests and flanks, promoting increased frequency of landslides that lasts for decades. More gradual impacts involve the flushing of excess debris downstream by rivers, which can generate bank erosion and floodplain accretion as well as channel avulsions that affect flooding frequency, settlements, ecosystems, and infrastructure. Ultimately, earthquake sequences and their geomorphic consequences alter mountain landscapes over both human and geologic time scales. Two recent events have attracted intense research into earthquake‐induced landslides and their consequences: the magnitude M 7.6 Chi‐Chi, Taiwan earthquake of 1999, and the M 7.9 Wenchuan, China earthquake of 2008. Using data and insights from these and several other earthquakes, we analyze how such events initiate processes that change mountain landscapes, highlight research gaps, and suggest pathways toward a more complete understanding of the seismic effects on the Earth's surface.
Plate tectonic boundaries constitute the suture zones between tectonic plates. They are shaped by a variety of distinct and interrelated processes and play a key role in geohazards and georesource formation. Many of these processes have been previously studied, while many others remain unaddressed or undiscovered. In this work, the geodynamic numerical modeling software ASPECT is applied to shed light on further process interactions at continental plate boundaries. In contrast to natural data, geodynamic modeling has the advantage that processes can be directly quantified and that all parameters can be analyzed over the entire evolution of a structure. Furthermore, processes and interactions can be singled out from complex settings because the modeler has full control over all of the parameters involved. To account for the simplifying character of models in general, I have chosen to study generic geological settings with a focus on the processes and interactions rather than precisely reconstructing a specific region of the Earth.
In Chapter 2, 2D models of continental rifts with different crustal thicknesses between 20 and 50 km and extension velocities in the range of 0.5-10 mm/yr are used to obtain a speed limit for the thermal steady-state assumption, commonly employed to address the temperature fields of continental rifts worldwide. Because the tectonic deformation from ongoing rifting outpaces heat conduction, the temperature field is not in equilibrium, but is characterized by a transient, tectonically-induced heat flow signal. As a result, I find that isotherm depths of the geodynamic evolution models are shallower than a temperature distribution in equilibrium would suggest. This is particularly important for deep isotherms and narrow rifts. In narrow rifts, the magnitude of the transient temperature signal limits a well-founded applicability of the thermal steady-state assumption to extension velocities of 0.5-2 mm/yr. Estimation of the crustal temperature field affects conclusions on all temperature-dependent processes ranging from mineral assemblages to the feasible exploitation of a geothermal reservoir.
In Chapter 3, I model the interactions of different rheologies with the kinematics of folding and faulting using the example of fault-propagation folds in the Andean fold-and-thrust belt. The evolution of the velocity fields from geodynamic models are compared with those from trishear models of the same structure. While the latter use only geometric and kinematic constraints of the main fault, the geodynamic models capture viscous, plastic, and elastic deformation in the entire model domain. I find that both models work equally well for early, and thus relatively simple stages of folding and faulting, while results differ for more complex situations where off-fault deformation and secondary faulting are present. As fault-propagation folds can play an important role in the formation of reservoirs, knowledge of fluid pathways, for example via fractures and faults, is crucial for their characterization.
Chapter 4 deals with a bending transform fault and the interconnections between tectonics and surface processes. In particular, the tectonic evolution of the Dead Sea Fault is addressed where a releasing bend forms the Dead Sea pull-apart basin, while a restraining bend further to the North resulted in the formation of the Lebanese mountains. I ran 3D coupled geodynamic and surface evolution models that included both types of bends in a single setup. I tested various randomized initial strain distributions, showing that basin asymmetry is a consequence of strain localization. Furthermore, by varying the surface process efficiency, I find that the deposition of sediment in the pull-apart basin not only controls basin depth, but also results in a crustal flow component that increases uplift at the restraining bend.
Finally, in Chapter 5, I present the computational basis for adding further complexity to plate boundary models in ASPECT with the implementation of earthquake-like behavior using the rate-and-state friction framework. Despite earthquakes happening on a relatively small time scale, there are many interactions between the seismic cycle and the long time spans of other geodynamic processes. Amongst others, the crustal state of stress as well as the presence of fluids or changes in temperature may alter the frictional behavior of a fault segment. My work provides the basis for a realistic setup of involved structures and processes, which is therefore important to obtain a meaningful estimate for earthquake hazards.
While these findings improve our understanding of continental plate boundaries, further development of geodynamic software may help to reveal even more processes and interactions in the future.
For over a century, geomorphologists have attempted to unravel information about landscape evolution, and processes that drive it, using river profiles. Many studies have combined new topographic datasets with theoretical models of channel incision to infer erosion rates, identify rock types with different resistance to erosion, and detect potential regions of tectonic activity. The most common metric used to analyse river profile geometry is channel steepness, or k(s). However, the calculation of channel steepness requires the normalisation of channel gradient by drainage area. This normalisation requires a power law exponent that is referred to as the channel concavity index. Despite the concavity index being crucial in determining channel steepness, it is challenging to constrain. In this contribution, we compare both slope-area methods for calculating the concavity index and methods based on integrating drainage area along the length of the channel, using so-called "chi" (chi) analysis. We present a new chi-based method which directly compares chi values of tributary nodes to those on the main stem; this method allows us to constrain the concavity index in transient landscapes without assuming a linear relationship between chi and elevation. Patterns of the concavity index have been linked to the ratio of the area and slope exponents of the stream power incision model (m/n); we therefore construct simple numerical models obeying detachment-limited stream power and test the different methods against simulations with imposed m and n. We find that chi-based methods are better than slope-area methods at reproducing imposed m/n ratios when our numerical landscapes are subject to either transient uplift or spatially varying uplift and fluvial erodibility. We also test our methods on several real landscapes, including sites with both lithological and structural heterogeneity, to provide examples of the methods' performance and limitations. These methods are made available in a new software package so that other workers can explore how the concavity index varies across diverse landscapes, with the aim to improve our understanding of the physics behind bedrock channel incision.
The response of surface processes to climatic forcing is fundamental for understanding the impacts of climate change on landscape evolution. In the Himalaya, most large rivers feature prominent fill terraces that record an imbalance between sediment supply and transport capacity, presumably due to past fluctuations in monsoon precipitation and/or effects of glaciation at high elevation. Here, we present volume estimates, chronological constraints, and Be-10-derived paleo-erosion rates from a prominent valley fill in the Yamuna catchment, Garhwal Himalaya, to elucidate the coupled response of rivers and hillslopes to Pleistocene climate change. Although precise age control is complicated due to methodological problems, the new data support formation of the valley fill during the late Pleistocene and its incision during the Holocene. We interpret this timing to indicate that changes in discharge and river-transport capacity were major controls. Compared to the present day, late Pleistocene hillslope erosion rates were higher by a factor of similar to 2-4, but appear to have decreased during valley aggradation. The higher late Pleistocene erosion rates are largely unrelated to glacial erosion and could be explained by enhanced sediment production on steep hillslopes due to increased periglacial activity that declined as temperatures increased. Alternatively, erosion rates that decrease during valley aggradation are also consistent with reduced landsliding from threshold hillslopes as a result of rising base levels. In that case, the similarity of paleo-erosion rates near the end of the aggradation period with modern erosion rates might imply that channels and hillslopes are not yet fully coupled everywhere and that present-day hillslope erosion rates may underrepresent long-term incision rates. (C) 2015 Elsevier B.V. All rights reserved.
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.
Piggyback basins on the margins of growing orogens commonly serve as sensitive recorders of the onset of thrust deformation and changes in source areas. The Bieertuokuoyi piggyback basin, located in the hanging wall of the Pamir Frontal Thrust, provides an unambiguous record of the outward growth of the northeast Pamir margin in northwest China from the Miocene through the Quaternary. To reconstruct the deformation along the margin, we synthesized structural mapping, stratigraphy, magnetostratigraphy, and cosmogenic burial dating of basin fill and growth strata. The Bieertuokuoyi basin records the initiation of the Pamir Frontal Thrust and the Takegai Thrust similar to 5-6Ma, as well as clast provenance and paleocurrent changes resulting from the Pliocene-to-Recent uplift and exhumation of the Pamir to the south. Our results show that coeval deformation was accommodated on the major structures on the northeast Pamir margin throughout the Miocene to Recent. Furthermore, our data support a change in the regional kinematics around the Miocene-Pliocene boundary (similar to 5-6Ma). Rapid exhumation of NE Pamir extensional domes, coupled with cessation of the Kashgar-Yecheng Transfer System on the eastern margin of the Pamir, accelerated the outward propagation of the northeastern Pamir margin and the southward propagation of the Kashi-Atushi fold-and-thrust belt in the southern Tian Shan. This coeval deformation signifies the coupling of the Pamir and Tarim blocks and the transfer of shortening north to the Pamir frontal faults and across the quasi-rigid Tarim Basin to the southern Tian Shan Kashi-Atushi fold-and-thrust system.
Drainage-divide migration, controlled by rock-uplift and rainfall patterns, may play a major role in the geomorphic evolution of mountain ranges.
However, divide-migration rates over geologic timescales have only been estimated by theoretical studies and remain empirically poorly constrained.
Geomorphological evidence suggests that the Sierra de Aconquija, on the eastern side of the southern Central Andes, northwest Argentina, is undergoing active westward drainage-divide migration. The mountain range has been subjected to steep rock trajectories and pronounced orographic rainfall for the last several million years, presenting an ideal setting for using low-temperature thermochronometric data to explore its topographic evolution.
We perform three-dimensional thermal-kinematic modeling of previously published thermochronometric data spanning the windward and leeward sides of the range to explore the most likely structural and topographic evolution of the range.
We find that the data can be explained by scenarios involving drainage-divide migration alone, or by scenarios that also involve changes in the structures that have accommodated deformation through time.
By combining new Be-10-derived catchment-average denudation rates with geomorphic constraints on probable fault activity, we conclude that the evolution of the range was likely dominated by west-vergent faulting on a high-angle reverse fault underlying the range, together with westward drainage-divide migration at a rate of several km per million years.
Our findings place new constraints on the magnitudes and rates of drainage-divide migration in real landscapes, quantify the effects of orographic rainfall and erosion on the topographic evolution of a mountain range, and highlight the importance of considering drainage-divide migration when interpreting thermochronometer age patterns.