Refine
Year of publication
Document Type
- Article (19)
- Doctoral Thesis (10)
Language
- English (29)
Is part of the Bibliography
- yes (29)
Keywords
- Himalaya (29) (remove)
Institute
Fluvial terraces, floodplains, and alluvial fans are the main landforms to store sediments and to decouple hillslopes from eroding mountain rivers. Such low-relief landforms are also preferred locations for humans to settle in otherwise steep and poorly accessible terrain. Abundant water and sediment as essential sources for buildings and infrastructure make these areas amenable places to live at. Yet valley floors are also prone to rare and catastrophic sedimentation that can overload river systems by abruptly increasing the volume of sediment supply, thus causing massive floodplain aggradation, lateral channel instability, and increased flooding. Some valley-fill sediments should thus record these catastrophic sediment pulses, allowing insights into their timing, magnitude, and consequences.
This thesis pursues this theme and focuses on a prominent ~150 km2 valley fill in the Pokhara Valley just south of the Annapurna Massif in central Nepal. The Pokhara Valley is conspicuously broad and gentle compared to the surrounding dissected mountain terrain,
and is filled with locally more than 70 m of clastic debris. The area’s main river, Seti Khola, descends from the Annapurna Sabche Cirque at 3500-4500 m asl down to 900 m asl where it incises into this valley fill. Humans began to settle on this extensive
fan surface in the 1750’s when the Trans-Himalayan trade route connected the Higher Himalayas, passing Pokhara city, with the subtropical lowlands of the Terai. High and unstable river terraces and steep gorges undermined by fast flowing rivers with highly seasonal (monsoon-driven) discharge, a high earthquake risk, and a growing population make the Pokhara Valley an ideal place to study the recent geological and geomorphic history of its sediments and the implication for natural hazard appraisals.
The objective of this thesis is to quantify the timing, the sedimentologic and geomorphic processes as well as the fluvial response to a series of strong sediment pulses. I report
diagnostic sedimentary archives, lithofacies of the fan terraces, their geochemical provenance, radiocarbon-age dating and the stratigraphic relationship between them. All these various and independent lines of evidence show consistently that multiple sediment pulses filled the Pokhara Valley in medieval times, most likely in connection with, if not triggered by, strong seismic ground shaking. The geomorphic and sedimentary evidence is
consistent with catastrophic fluvial aggradation tied to the timing of three medieval Himalayan earthquakes in ~1100, 1255, and 1344 AD. Sediment provenance and calibrated radiocarbon-age data are the key to distinguish three individual sediment pulses, as these are not evident from their sedimentology alone. I explore various measures of adjustment and fluvial response of the river system following these massive aggradation pulses. By using proxies such as net volumetric erosion, incision and erosion rates, clast provenance on active river banks, geomorphic markers such as re-exhumed tree trunks in growth position, and knickpoint locations in tributary valleys, I estimate the response of the river network in the Pokhara Valley to earthquake disturbance over several centuries. Estimates of the removed volumes since catastrophic valley filling began, require average net sediment
yields of up to 4200 t km−2 yr−1 since, rates that are consistent with those reported for Himalayan rivers. The lithological composition of active channel-bed load differs from that of local bedrock material, confirming that rivers have adjusted 30-50% depending on data of different tributary catchments, locally incising with rates of 160-220 mm yr−1. In many tributaries to the Seti Khola, most of the contemporary river loads come from a Higher Himalayan source, thus excluding local hillslopes as sources. This imbalance in sediment provenance emphasizes how the medieval sediment pulses must have rapidly traversed up to 70 km downstream to invade the downstream reaches of the tributaries
up to 8 km upstream, thereby blocking the local drainage and thus reinforcing, or locally creating new, floodplain lakes still visible in the landscape today.
Understanding the formation, origin, mechanism and geomorphic processes of this valley fill is crucial to understand the landscape evolution and response to catastrophic sediment pulses. Several earthquake-triggered long-runout rock-ice avalanches or catastrophic dam burst in the Higher Himalayas are the only plausible mechanisms to explain both the geomorphic and sedimentary legacy that I document here. In any case, the Pokhara Valley was most likely hit by a cascade of extremely rare processes over some two centuries starting in the early 11th century. Nowhere in the Himalayas do we find valley fills of
comparable size and equally well documented depositional history, making the Pokhara Valley one of the most extensively dated valley fill in the Himalayas to date. Judging from the growing record of historic Himalayan earthquakes in Nepal that were traced and
dated in fault trenches, this thesis shows that sedimentary archives can be used to directly aid reconstructions and predictions of both earthquake triggers and impacts from a sedimentary-response perspective. The knowledge about the timing, evolution, and response of the Pokhara Valley and its river system to earthquake triggered sediment pulses is important to address the seismic and geomorphic risk for the city of Pokhara. This
thesis demonstrates how geomorphic evidence on catastrophic valley infill can help to independently verify paleoseismological fault-trench records and may initiate re-thinking on post-seismic hazard assessments in active mountain regions.
In the high mountains of Asia, glaciers cover an area of approximately 115,000 km² and constitute one of the largest continental ice accumulations outside Greenland and Antarctica. Their sensitivity to climate change makes them valuable palaeoclimate archives, but also vulnerable to current and predicted Global Warming. This is a pressing problem as snow and glacial melt waters are important sources for agriculture and power supply of densely populated regions in south, east, and central Asia. Successful prediction of the glacial response to climate change in Asia and mitigation of the socioeconomic impacts requires profound knowledge of the climatic controls and the dynamics of Asian glaciers. However, due to their remoteness and difficult accessibility, ground-based studies are rare, as well as temporally and spatially limited. We therefore lack basic information on the vast majority of these glaciers. In this thesis, I employ different methods to assess the dynamics of Asian glaciers on multiple time scales. First, I tested a method for precise satellite-based measurement of glacier-surface velocities and conducted a comprehensive and regional survey of glacial flow and terminus dynamics of Asian glaciers between 2000 and 2008. This novel and unprecedented dataset provides unique insights into the contrasting topographic and climatic controls of glacial flow velocities across the Asian highlands. The data document disparate recent glacial behavior between the Karakoram and the Himalaya, which I attribute to the competing influence of the mid-latitude westerlies during winter and the Indian monsoon during summer. Second, I tested whether such climate-related longitudinal differences in glacial behavior also prevail on longer time scales, and potentially account for observed regionally asynchronous glacial advances. I used cosmogenic nuclide surface exposure dating of erratic boulders on moraines to obtain a glacial chronology for the upper Tons Valley, situated in the headwaters of the Ganges River. This area is located in the transition zone from monsoonal to westerly moisture supply and therefore ideal to examine the influence of these two atmospheric circulation regimes on glacial advances. The new glacial chronology documents multiple glacial oscillations during the last glacial termination and during the Holocene, suggesting largely synchronous glacial changes in the western Himalayan region that are related to gradual glacial-interglacial temperature oscillations with superimposed monsoonal precipitation changes of higher frequency. In a third step, I combine results from short-term satellite-based climate records and surface velocity-derived ice-flux estimates, with topographic analyses to deduce the erosional impact of glaciations on long-term landscape evolution in the Himalayan-Tibetan realm. The results provide evidence for the long-term effects of pronounced east-west differences in glaciation and glacial erosion, depending on climatic and topographic factors. Contrary to common belief the data suggest that monsoonal climate in the central Himalaya weakens glacial erosion at high elevations, helping to maintain a steep southern orographic barrier that protects the Tibetan Plateau from lateral destruction. The results of this thesis highlight how climatic and topographic gradients across the high mountains of Asia affect glacier dynamics on time scales ranging from 10^0 to 10^6 years. Glacial response times to climate changes are tightly linked to properties such as debris cover and surface slope, which are controlled by the topographic setting, and which need to be taken into account when reconstructing mountainous palaeoclimate from glacial histories or assessing the future evolution of Asian glaciers. Conversely, the regional topographic differences of glacial landscapes in Asia are partly controlled by climatic gradients and the long-term influence of glaciers on the topographic evolution of the orogenic system.
Anthropogenic climate change alters the hydrological cycle. While certain areas experience more intense precipitation events, others will experience droughts and increased evaporation, affecting water storage in long-term reservoirs, groundwater, snow, and glaciers. High elevation environments are especially vulnerable to climate change, which will impact the water supply for people living downstream. The Himalaya has been identified as a particularly vulnerable system, with nearly one billion people depending on the runoff in this system as their main water resource. As such, a more refined understanding of spatial and temporal changes in the water cycle in high altitude systems is essential to assess variations in water budgets under different climate change scenarios.
However, not only anthropogenic influences have an impact on the hydrological cycle, but changes to the hydrological cycle can occur over geological timescales, which are connected to the interplay between orogenic uplift and climate change. However, their temporal evolution and causes are often difficult to constrain. Using proxies that reflect hydrological changes with an increase in elevation, we can unravel the history of orogenic uplift in mountain ranges and its effect on the climate.
In this thesis, stable isotope ratios (expressed as δ2H and δ18O values) of meteoric waters and organic material are combined as tracers of atmospheric and hydrologic processes with remote sensing products to better understand water sources in the Himalayas. In addition, the record of modern climatological conditions based on the compound specific stable isotopes of leaf waxes (δ2Hwax) and brGDGTs (branched Glycerol dialkyl glycerol tetraethers) in modern soils in four Himalayan river catchments was assessed as proxies of the paleoclimate and (paleo-) elevation. Ultimately, hydrological variations over geological timescales were examined using δ13C and δ18O values of soil carbonates and bulk organic matter originating from sedimentological sections from the pre-Siwalik and Siwalik groups to track the response of vegetation and monsoon intensity and seasonality on a timescale of 20 Myr.
I find that Rayleigh distillation, with an ISM moisture source, mainly controls the isotopic composition of surface waters in the studied Himalayan catchments. An increase in d-excess in the spring, verified by remote sensing data products, shows the significant impact of runoff from snow-covered and glaciated areas on the surface water isotopic values in the timeseries.
In addition, I show that biomarker records such as brGDGTs and δ2Hwax have the potential to record (paleo-) elevation by yielding a significant correlation with the temperature and surface water δ2H values, respectively, as well as with elevation. Comparing the elevation inferred from both brGDGT and δ2Hwax, large differences were found in arid sections of the elevation transects due to an additional effect of evapotranspiration on δ2Hwax. A combined study of these proxies can improve paleoelevation estimates and provide recommendations based on the results found in this study.
Ultimately, I infer that the expansion of C4 vegetation between 20 and 1 Myr was not solely dependent on atmospheric pCO2, but also on regional changes in aridity and seasonality from to the stable isotopic signature of the two sedimentary sections in the Himalaya (east and west).
This thesis shows that the stable isotope chemistry of surface waters can be applied as a tool to monitor the changing Himalayan water budget under projected increasing temperatures. Minimizing the uncertainties associated with the paleo-elevation reconstructions were assessed by the combination of organic proxies (δ2Hwax and brGDGTs) in Himalayan soil. Stable isotope ratios in bulk soil and soil carbonates showed the evolution of vegetation influenced by the monsoon during the late Miocene, proving that these proxies can be used to record monsoon intensity, seasonality, and the response of vegetation. In conclusion, the use of organic proxies and stable isotope chemistry in the Himalayas has proven to successfully record changes in climate with increasing elevation. The combination of δ2Hwax and brGDGTs as a new proxy provides a more refined understanding of (paleo-)elevation and the influence of climate.
In the Himalaya of Chamba, NW India, a major orographic barrier in front of the Greater Himalayan Range extracts a high proportion of the monsoonal rainfall along its southern slopes and effectively shields the orogen interior from moisture-bearing winds. Along a similar to 100-km-long orogen perpendicular transect, 28 new apatite fission track (AFT) and 30 new zircon (U-Th)/He (ZHe) cooling ages reveal marked variations in age distributions and long-term exhumation rates between the humid frontal range and the semi-arid orogen interior. On the southern topographic front, very young, elevation-invariant AFT ages of <4 Ma have been obtained that are concentrated in a similar to 30-km-wide zone; 1-D-thermal modeling suggests a Plio-Pleistocene mean erosion rate of 0.8-1.9 mm yr(-1). In contrast, AFT and ZHe ages within the orogen interior are older (4-9 and 7-18 Ma, respectively), are positively correlated with sample elevation, and yield lower mean erosion rates (0.3-0.9 mm yr(-1)). Protracted low exhumation rates within the orogen interior over the last similar to 15 Myr prevailed contemporaneously with overall humid conditions and an effective erosional regime within the southern Himalaya. This suggests that the frontal Dhauladar Range was sufficiently high during this time to form an orographic barrier, focusing climatically enhanced erosional processes and tectonic deformation there. Thrusting along the two frontal range-bounding thrust, the Main Central Thrust and the Main Boundary Thrusts, was initiated at least similar to 15 Ma ago and has remained localized since then. The lack of evidence for localized uplift farther north indicates either a rather flat decollement with no ramp or the absence of active duplex systems beneath the interior of Chamba. Exhumational variability within Chamba is best explained as the result of continuous thrusting along a major basal decollement, with a flat beneath the slowly exhuming internal compartments and a steep frontal ramp at the rapidly exhuming frontal range. The pattern in Chamba contrasts with what is observed elsewhere along the Himalaya, where exhumation is focused in a zone similar to 150 km north of the orogenic front. In the NW Himalaya, preserved High Himalayan Crystalline nappes and Lesser Himalayan windows alternate on a relatively small scale of <100 km; these alternations are closely correlated with the pattern of exhumation. Although the spatial distribution of high-exhumation zones varies considerably between individual Himalayan sectors, all of these zones are closely correlated with locally higher rock-uplift rates, sharp topographic discontinuities, and focused orographic precipitation, suggesting strong feedbacks between tectonically driven rock uplift, orographically enhanced precipitation, and erosional processes.
Formation of a Rain Shadow
(2018)
We measure the oxygen and hydrogen stable isotope composition of authigenic clays from Himalayan foreland sediments (Siwalik Group), and from present day small stream waters in eastern Bhutan to explore the impact of uplift of the Shillong Plateau on rain shadow formation over the Himalayan foothills. Stable isotope data from authigenic clay minerals (<2 μm) suggest the presence of three paleoclimatic periods during deposition of the Siwalik Group, between ∼7 and ∼1 Ma. The mean δ18O value in paleometeoric waters, which were in equilibrium with clay minerals, is ∼2.5‰ lower than in modern meteoric and stream waters at the elevation of the foreland basin. We discuss the factors that could have changed the isotopic composition of water over time and we conclude that (a) the most likely and significant cause for the increase in meteoric water δ18O values over time is the “amount effect,” specifically, a decrease in mean annual precipitation. (b) The change in mean annual precipitation over the foreland basin and foothills of the Himalaya is the result of orographic effect caused by the Shillong Plateau's uplift. The critical elevation of the Shillong Plateau required to induce significant orographic precipitation was attained after ∼1.2 Ma. (c) By applying scale analysis, we estimate that the mean annual precipitation over the foreland basin of the eastern Bhutan Himalayas has decreased by a factor of 1.7–2.5 over the last 1–3 million years.
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.
The Himalayan arc stretches >2500 km from east to west at the southern edge of the Tibetan Plateau, representing one of the most important Cenozoic continent-continent collisional orogens. Internal deformation processes and climatic factors, which drive weathering, denudation, and transport, influence the growth and erosion of the orogen. During glacial times wet-based glaciers sculpted the mountain range and left overdeepend and U-shaped valleys, which were backfilled during interglacial times with paraglacial sediments over several cycles. These sediments partially still remain within the valleys because of insufficient evacuation capabilities into the foreland. Climatic processes overlay long-term tectonic processes responsible for uplift and exhumation caused by convergence. Possible processes accommodating convergence within the orogenic wedge along the main Himalayan faults, which divide the range into four major lithologic units, are debated. In this context, the identification of processes shaping the Earth’s surface on short- and on long-term are crucial to understand the growth of the orogen and implications for landscape development in various sectors along the arc. This thesis focuses on both surface and tectonic processes that shape the landscape in the western Indian Himalaya since late Miocene.
In my first study, I dated well-preserved glacially polished bedrock on high-elevated ridges and valley walls in the upper of the Chandra Valley the by means of 10Be terrestrial cosmogenic radionuclides (TCN). I used these ages and mapped glacial features to reconstruct the extent and timing of Pleistocene glaciation at the southern front of the Himalaya. I was able to reconstruct an extensive valley glacier of ~200 km length and >1000 m thickness. Deglaciation of the Chandra Valley glacier started subsequently to insolation increase on the Northern Hemisphere and thus responded to temperature increase. I showed that the timing this deglaciation onset was coeval with retreat of further midlatitude glaciers on the Northern and Southern Hemispheres. These comparisons also showed that the post-LGM deglaciation very rapid, occurred within a few thousand years, and was nearly finished prior to the Bølling/Allerød interstadial.
A second study (co-authorship) investigates how glacial advances and retreats in high mountain environments impact the landscape. By 10Be TCN dating and geomorphic mapping, we obtained maximal length and height of the Siachen Glacier within the Nubra Valley. Today the Shyok and Nubra confluence is backfilled with sedimentary deposits, which are attributed to the valley blocking of the Siachen Glacier 900 m above the present day river level. A glacial dam of the Siachen Glacier blocked the Shyok River and lead to the evolution of a more than 20 km long lake. Fluvial and lacustrine deposits in the valley document alternating draining and filling cycles of the lake dammed by the Siachen Glacier. In this study, we can show that glacial incision was outpacing fluvial incision.
In the third study, which spans the million-year timescale, I focus on exhumation and erosion within the Chandra and Beas valleys. In this study the position and discussed possible reasons of rapidly exhuming rocks, several 100-km away from one of the main Himalayan faults (MFT) using Apatite Fission Track (AFT) thermochronometry. The newly gained AFT ages indicate rapid exhumation and confirm earlier studies in the Chandra Valley. I assume that the rapid exhumation is most likely related to uplift over subsurface structures. I tested this hypothesis by combining further low-temperature thermochronometers from areas east and west of my study area. By comparing two transects, each parallel to the Beas/Chandra Valley transect, I demonstrate similarities in the exhumation pattern to transects across the Sutlej region, and strong dissimilarities in the transect crossing the Dhauladar Range. I conclude that the belt of rapid exhumation terminates at the western end of the Kullu-Rampur window. Therewith, I corroborate earlier studies suggesting changes in exhumation behavior in the western Himalaya. Furthermore, I discussed several causes responsible for the pronounced change in exhumation patterns along strike: 1) the role of inherited pre-collisional features such as the Proterozoic sedimentary cover of the Indian basement, former ridges and geological structures, and 2) the variability of convergence rates along the Himalayan arc due to an increased oblique component towards the syntaxis.
The combination of field observations (geological and geomorphological mapping) and methods to constrain short- and long-term processes (10Be, AFT) help to understand the role of the individual contributors to exhumation and erosion in the western Indian Himalaya. With the results of this thesis, I emphasize the importance of glacial and tectonic processes in shaping the landscape by driving exhumation and erosion in the studied areas.
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.
The India-Eurasia continental collision zone provides a spectacular example of active mountain building and climatic forcing. In order to quantify the critically important process of mass removal, I analyzed spatial and temporal precipitation patterns of the oscillating monsoon system and their geomorphic imprints. I processed passive microwave satellite data to derive high-resolution rainfall estimates for the last decade and identified an abnormal monsoon year in 2002. During this year, precipitation migrated far into the Sutlej Valley in the northwestern part of the Himalaya and reached regions behind orographic barriers that are normally arid. There, sediment flux, mean basin denudation rates, and channel-forming processes such as erosion by debris-flows increased significantly. Similarly, during the late Pleistocene and early Holocene, solar forcing increased the strength of the Indian summer monsoon for several millennia and presumably lead to analogous precipitation distribution as were observed during 2002. However, the persistent humid conditions in the steep, high-elevation parts of the Sutlej River resulted in deep-seated landsliding. Landslides were exceptionally large, mainly due to two processes that I infer for this time: At the onset of the intensified monsoon at 9.7 ka BP heavy rainfall and high river discharge removed material stored along the river, and lowered the baselevel. Second, enhanced discharge, sediment flux, and increased pore-water pressures along the hillslopes eventually lead to exceptionally large landslides that have not been observed in other periods. The excess sediments that were removed from the upstream parts of the Sutlej Valley were rapidly deposited in the low-gradient sectors of the lower Sutlej River. Timing of downcutting correlates with centennial-long weaker monsoon periods that were characterized by lower rainfall. I explain this relationship by taking sediment flux and rainfall dynamics into account: High sediment flux derived from the upstream parts of the Sutlej River during strong monsoon phases prevents fluvial incision due to oversaturation the fluvial sediment-transport capacity. In contrast, weaker monsoons result in a lower sediment flux that allows incision in the low-elevation parts of the Sutlej River.
Late Quaternary valley infill and dissection in the Indus River, western Tibetan Plateau margin
(2014)
The Indus, one of Earth's major rivers, drains large parts of the NW Himalaya and the Transhimalayan ranges that form part of the western Tibetan Plateau margin. In the western Himalayan syntaxis, where local topographic relief exceeds 7 km, the Indus has incised a steep bedrock gorge at rates of several mm yr(-1). Upstream, however, the upper Indus and its tributaries alternate between bedrock gorges and broad alluvial flats flanked by the Ladakh and Zanskar ranges. We review the late Quaternary valley history in this region with a focus on the confluence of the Indus and Zanskar Rivers, where vast alluvial terrace staircases and lake sediments record major episodes of aggradation and incision. New absolute dating of high-level fluvial terrace remnants using cosmogenic Be-10, optically and infrared stimulated luminescence (OSL, IRSL) indicates at least two phases of late Quaternary valley infilling. These phases commenced before similar to 200 ka and similar to 50-20 ka, judging from terrace treads stranded >150 m and similar to 30-40 m above modern river levels, respectively. Numerous stacks of lacustrine sediments that straddle the Indus River >200 km between the city of Leh and the confluence with the Shyok River share a distinct horizontal alignment. Constraints from IRSL samples of lacustrine sequences from the Leh-Spituk area reveal a protracted lake phase from >177 ka to 72 ka, locally accumulating >50-m thick deposits. In the absence of tectonic faulting, major lithological differences, and stream capture, we attribute the formation of this and other large lakes in the region to natural damming by large landslides, glaciers, and alluvial fans. The overall patchy landform age constraints from earlier studies can be reconciled by postulating a major deglacial control on sediment flux, valley infilling, and subsequent incision that has been modulated locally by backwater effects of natural damming. While comparison with Pleistocene monsoon proxies reveals no obvious correlation, a lateor post-glacial sediment pulse seems a more likely source of this widespread sedimentation that has partly buried the dissected bedrock topography. Overall, the long residence times of fluvial, alluvial and lacustrine deposits in the region (>500 ka) support previous studies, but remain striking given the dominantly steep slopes and deeply carved valleys that characterise this high-altitude mountain desert. Recalculated late Quaternary rates of fluvial bedrock incision in the Indus and Zanskar of 1.5 +/- 0.2 mm yr(-1) are at odds with the longevity of juxtaposed valley-fill deposits, unless a lack of decisive lateral fluvial erosion helps to preserve these late Pleistocene sedimentary archives. We conclude that alternating, similar to 10(4)-yr long, phases of massive infilling and incision have dominated the late Quaternary history of the Indus valley below the western Tibetan Plateau margin. (C) 2014 Elsevier Ltd. All rights reserved.
The Himalaya has a major influence on global and regional climate, in particular on the Asian monsoon system. The foreland basin of the Himalaya contains a record of tectonics and paleoclimate since the Miocene. Previous work on the evolution of vegetation and climate has focused on the central and western Himalaya, where a shift from C3 to C4 vegetation has been observed at similar to 7 Ma and linked to increased seasonality, but the climatic evolution of the eastern part of the orogen is less well understood. In order to track vegetation as a marker of monsoon intensity and seasonality, we analyzed delta C-13 and 8180 values of soil carbonate and associated delta C-13 values of bulk organic carbon from previously dated sedimentary sections exposing the syn-orogenic detrital Dharamsala and Siwalik Groups in the west, and, for the first time, the Siwalik Group in the east of the Himalayan foreland basin. Sedimentary records span from 20 to 1 Myr in the west (Joginder Nagar, Jawalamukhi, and Haripur Kolar sections) and from 13 to 1 Myr in the east (Kameng section), respectively. The presence of soil carbonate in the west and its absence in the east is a first indication of long-term lateral climatic variation, as soil carbonate requires seasonally arid conditions to develop. delta C-13 values in soil carbonate show a shift from around -10 parts per thousand to -2 parts per thousand at similar to 7 Ma in the west, which is confirmed by delta C-13 analyses on bulk organic carbon that show a shift from around -23 parts per thousand to -19 parts per thousand at the same time. Such a shift in isotopic values is likely to be associated with a change from C3 to C4 vegetation. In contrast, delta C-13 values of bulk organic carbon remain at 23 parts per thousand o in the east. Thus, our data show that the current east -west variation in climate was established at similar to 7 Ma. We propose that the regional change towards a more seasonal climate in the west is linked to a decrease of the influence of the Westerlies, delivering less winter precipitation to the western Himalaya, while the east remained annually humid due to its proximity to the monsoonal moisture source. (C) 2017 Elsevier B.V. All rights reserved.
The transfer of particulate organic carbon from continents to the ocean is an important component of the global carbon cycle. Transfer to and burial of photosynthetically fixed biospheric organic carbon in marine sediments can effectively sequester atmospheric carbon dioxide over geological timescales. The exhumation and erosion of fossil organic carbon contained in sedimentary rocks, i.e. petrogenic carbon, can result in remineralization, releasing carbon to the atmosphere. In contrast, eroded petrogenic organic carbon that gets transferred back to the ocean and reburied does not affect atmospheric carbon content.
Mountain ranges play a key role in this transfer since they can source vast amounts of sediment including particulate organic carbon. Globally, the export of both, biospheric and petrogenic organic carbon has been linked to sediment export. Additionally, short transfer times from mountains to the ocean and high sediment concentrations have been shown to increase the likelihood of organic carbon burial. While the importance of mountain ranges in the organic carbon cycle is now widely recognized, the processes acting within mountain ranges to influence the storage, cycling and mobilization of organic carbon, as well as carbon fluxes from mountain ranges remain poorly constrained.
In this thesis, I employ different methods to assess the nature and fate of particulate organic carbon in mountain belts, ranging from the molecular to regional landscape scale. These studies are located along the Trans-Himalayan Kali Gandaki River in Central Nepal. This river traverses all major geological and climatic zones of the Himalaya, from the dry northern Tibetan plateau to the high-relief, monsoon dominated steep High Himalaya and the lower relief and abundant vegetation of the Lesser Himalayan region.
First, I document how biospheric organic matter has accumulated during the Holocene in the headwaters of the Kali Gandaki River valley, by combining compound specific isotope measurements with different dating methods and grain size data, and investigate the stability of this organic carbon reservoir on millennial timescales. I show, that around 1.6 ka an eco-geomorphic tipping point occurred leading to a destabilization of the landscape resulting in today’s high erosion rates and the excavation of the aged organic carbon reservoir. This study highlights the climatic and geomorphic controls on biospheric organic carbon storage and release from mountain ranges.
Second, I systematically investigate the spatial variation of particulate organic carbon fluxes across the Himalaya along the Kali Gandaki River, using bulk stable and radioactive isotopes combined with a new Bayesian modeling approach. The detailed dataset allows the distinction of aged and modern biospheric organic carbon as well as petrogenic organic carbon across the Himalayan mountain range and the investigation of the role of climatic and geomorphic factors in their riverine export. The data suggest a decoupling of the particulate organic carbon from the sediment yield along the Kali Gandaki River, partially driven by climatic and geomorphic processes. In contrast to the suspended sediment, a large part of the particulate organic carbon exported by the river originates from the Tibetan part of the catchment and is dominated by petrogenic organic carbon derived from Jurassic shales with only minor contributions of modern and aged biospheric organic carbon. These findings emphasize the importance of organic carbon source distribution and erosion mechanisms in determining the organic carbon export from mountain ranges.
In a third step, I explore the potential of ultra-high resolution mass spectrometry for particulate organic carbon transport studies. I have generated a novel and unprecedented high-resolution molecular dataset, which contains up to 103 molecular formulas of the lipid fraction of particulate organic matter for modern and aged biospheric carbon, petrogenic organic carbon and river sediments. First, I test if this dataset can be used to better resolve different organic carbon sources and to identify new geochemical tracers. Using multivariate statistics, I identify up to 10² characteristic molecular formulas for the major organic carbon sources in the upper part of the Kali Gandaki catchment, and trace their transfer from the surrounding landscape into the river sediment. Second, I test the potential of the molecular dataset to trace molecular transformations along source-to-sink pathways. I identify changes in molecular metrics derived from the dataset, which are characteristic of transformation processes during incorporation of litter into soil, the aging of soil material, and the mobilization of the organic carbon into the river. These two studies demonstrate that high-resolution molecular datasets open a promising analytical window on particulate organic carbon and can provide novel insights into the composition, sourcing and transformation of riverine particulate organic carbon.
Collectively, these studies advance our understanding of the processes contributing to the storage and mobilization of organic carbon in the Central Himalaya, the mountain belt that dominates global erosional fluxes. They do so by identifying the major sources of particulate organic carbon to the Trans-Himalayan Kali Gandaki River, by elucidating their sensitivity to climate and geomorphic processes, and by identifying some of the transformations of this material on the molecular scale. As a result, the thesis demonstrates that the amount and composition of organic carbon routed from mountain belts is a function of the dynamic interactions of geologic, biologic, geomorphic and climatic processes within the mountain belt. This understanding will ultimately help in answering whether the build-up and erosion of mountain ranges over geological time represents a net carbon source or sink to the atmosphere. Beyond this, the thesis contributes to our technical ability to characterize organic matter and attribute it to sources by scoping the potential of high-end molecular analysis.
The Himalayas are a region that is most dependent, but also frequently prone to hazards from changing meltwater resources. This mountain belt hosts the highest mountain peaks on earth, has the largest reserve of ice outside the polar regions, and is home to a rapidly growing population in recent decades. One source of hazard has attracted scientific research in particular in the past two decades: glacial lake outburst floods (GLOFs) occurred rarely, but mostly with fatal and catastrophic consequences for downstream communities and infrastructure. Such GLOFs can suddenly release several million cubic meters of water from naturally impounded meltwater lakes. Glacial lakes have grown in number and size by ongoing glacial mass losses in the Himalayas. Theory holds that enhanced meltwater production may increase GLOF frequency, but has never been tested so far. The key challenge to test this notion are the high altitudes of >4000 m, at which lakes occur, making field work impractical. Moreover, flood waves can attenuate rapidly in mountain channels downstream, so that many GLOFs have likely gone unnoticed in past decades. Our knowledge on GLOFs is hence likely biased towards larger, destructive cases, which challenges a detailed quantification of their frequency and their response to atmospheric warming. Robustly quantifying the magnitude and frequency of GLOFs is essential for risk assessment and management along mountain rivers, not least to implement their return periods in building design codes.
Motivated by this limited knowledge of GLOF frequency and hazard, I developed an algorithm that efficiently detects GLOFs from satellite images. In essence, this algorithm classifies land cover in 30 years (~1988–2017) of continuously recorded Landsat images over the Himalayas, and calculates likelihoods for rapidly shrinking water bodies in the stack of land cover images. I visually assessed such detected tell-tale sites for sediment fans in the river channel downstream, a second key diagnostic of GLOFs. Rigorous tests and validation with known cases from roughly 10% of the Himalayas suggested that this algorithm is robust against frequent image noise, and hence capable to identify previously unknown GLOFs. Extending the search radius to the entire Himalayan mountain range revealed some 22 newly detected GLOFs. I thus more than doubled the existing GLOF count from 16 previously known cases since 1988, and found a dominant cluster of GLOFs in the Central and Eastern Himalayas (Bhutan and Eastern Nepal), compared to the rarer affected ranges in the North. Yet, the total of 38 GLOFs showed no change in the annual frequency, so that the activity of GLOFs per unit glacial lake area has decreased in the past 30 years. I discussed possible drivers for this finding, but left a further attribution to distinct GLOF-triggering mechanisms open to future research.
This updated GLOF frequency was the key input for assessing GLOF hazard for the entire Himalayan mountain belt and several subregions. I used standard definitions in flood hydrology, describing hazard as the annual exceedance probability of a given flood peak discharge [m3 s-1] or larger at the breach location. I coupled the empirical frequency of GLOFs per region to simulations of physically plausible peak discharges from all existing ~5,000 lakes in the Himalayas. Using an extreme-value model, I could hence calculate flood return periods. I found that the contemporary 100-year GLOF discharge (the flood level that is reached or exceeded on average once in 100 years) is 20,600+2,200/–2,300 m3 s-1 for the entire Himalayas. Given the spatial and temporal distribution of historic GLOFs, contemporary GLOF hazard is highest in the Eastern Himalayas, and lower for regions with rarer GLOF abundance. I also calculated GLOF hazard for some 9,500 overdeepenings, which could expose and fill with water, if all Himalayan glaciers have melted eventually. Assuming that the current GLOF rate remains unchanged, the 100-year GLOF discharge could double (41,700+5,500/–4,700 m3 s-1), while the regional GLOF hazard may increase largest in the Karakoram.
To conclude, these three stages–from GLOF detection, to analysing their frequency and estimating regional GLOF hazard–provide a framework for modern GLOF hazard assessment. Given the rapidly growing population, infrastructure, and hydropower projects in the Himalayas, this thesis assists in quantifying the purely climate-driven contribution to hazard and risk from GLOFs.
Oxygen isotope compositions are reported for the first time for the Himalayan metabasites of the Kaghan Valley, Pakistan in this study. The highest metamorphic grades are recorded in the north of the valley, near the India-Asia collision boundary, in the form of high-pressure (HP: Group I) and ultrahigh-pressure (UHP: Group II) eclogites. The rocks show a step-wise decrease in grade from the UHP to HP eclogites and amphibolites. The protoliths of these metabasites were the Permian Panjal Trap basalts (ca. 267 +/- 2.4 Ma), which were emplaced along the northern margin of India when it was part of Gondwana. After the break-up of Gondwana, India drifted northward, subducted beneath Asia and underwent UHP metamorphism during the Eocene (ca. 45 +/- 1.2 Ma). At the regional scale, amphibolites, Group I and II eclogites yielded delta O-18 values of +5.84 and +5.91 parts per thousand, +1.66 to +424 parts per thousand, and -2.25 to +0.76 parts per thousand, respectively, relative to VSMOW. On a more local scale, within a single eclogite body, the delta O-18 values were the lowest (-2.25 to-1.44%.) in the central, the best preserved (least retrograded) parts, and show a systematic increase outward into more retrograded rocks, reaching up to +0.12 parts per thousand. These values are significantly lower than the typical mantle values for basalts of + 5.7 +/- 0.3 parts per thousand. The unusually low or negative delta O-18 values in Group II eclogites potentially resulted from hydrothermal alteration of the protoliths by interactions with meteoric water when the Indian plate was at southern high latitudes (similar to 60 degrees S). The stepwise increase in delta O-18 values, among different eclogite bodies in general and at single outcrop-scales in particular, reflects differing degrees of resetting of the oxygen isotope compositions during exhumation-related retrogression. (C) 2014 Elsevier B.V. All rights reserved.
The Himalayan thrust belt comprises three in-sequence foreland-propagating orogen-scale faults, the Main Central thrust, the Main Boundary thrust, and the Main Frontal thrust. Recently, the Munsiari–Ramgarh–Shumar thrust system has been recognized as an additional, potentially orogen-scale shear zone in the proximal footwall of the Main Central thrust. The timing of the Munsiari, Ramgarh, and Shumar thrusts and their role in Himalayan tectonics are disputed. We present 31 new zircon (U–Th)/He ages from a profile across the central Himachal Himalaya in the Beas River area. Within a ∼40 km wide belt northeast of the Kullu–Larji–Rampur window, ages ranging from to constrain a distinct episode of rapid Pliocene to Present exhumation; north and south of this belt, zircon (U–Th)/He ages are older ( to ). We attribute the Pliocene rapid exhumation episode to basal accretion to the Himalayan thrust belt and duplex formation in the Lesser Himalayan sequence including initiation of the Munsiari thrust. Pecube thermokinematic modelling suggests exhumation rates of ∼2–3 mm/yr from 4–7 to 0 Ma above the duplex contrasting with lower (<0.3 mm/yr) middle-late Miocene exhumation rates. The Munsiari thrust terminates laterally in central Himachal Pradesh. In the NW Indian Himalaya, the Main Central thrust zone comprises the sheared basal sections of the Greater Himalayan sequence and the mylonitic ‘Bajaura nappe’ of Lesser Himalayan affinity. We correlate the Bajaura unit with the Ramgarh thrust sheet in Nepal based on similar lithologies and the middle Miocene age of deformation. The Munsiari thrust in the central Himachal Himalaya is several Myr younger than deformation in the Bajaura and Ramgarh thrust sheets. Our results illustrate the complex and segmented nature of the Munsiari–Ramgarh–Shumar thrust system.
Mountain rivers respond to strong earthquakes by rapidly aggrading to accommodate excess sediment delivered by co-seismic landslides. Detailed sediment budgets indicate that rivers need several years to decades to recover from seismic disturbances, depending on how recovery is defined. We examine three principal proxies of river recovery after earthquake-induced sediment pulses around Pokhara, Nepal's second largest city. Freshly exhumed cohorts of floodplain trees in growth position indicate rapid and pulsed sedimentation that formed a fan covering 150 km2 in a Lesser Himalayan basin with tens of metres of debris between the 11th and 15th centuries AD. Radiocarbon dates of buried trees are consistent with those of nearby valley deposits linked to major medieval earthquakes, such that we can estimate average rates of re-incision since. We combine high-resolution digital elevation data, geodetic field surveys, aerial photos, and dated tree trunks to reconstruct geomorphic marker surfaces. The volumes of sediment relative to these surfaces require average net sediment yields of up to 4200 t km–2 yr–1 for the 650 years since the last inferred earthquake-triggered sediment pulse. The lithological composition of channel bedload differs from that of local bedrock, confirming that rivers are still mostly evacuating medieval valley fills, locally incising at rates of up to 0.2 m yr–1. Pronounced knickpoints and epigenetic gorges at tributary junctions further illustrate the protracted fluvial response; only the distal portions of the earthquake-derived sediment wedges have been cut to near their base. Our results challenge the notion that mountain rivers recover speedily from earthquakes within years to decades. The valley fills around Pokhara show that even highly erosive Himalayan rivers may need more than several centuries to adjust to catastrophic perturbations. Our results motivate some rethinking of post-seismic hazard appraisals and infrastructural planning in active mountain regions.
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.
Despite remarkable tectonostratigraphic similarities along the Himalayan arc, pronounced topographic and exhumational variability exists in different morphotectonic segments. The processes responsible for this segmentation are debated. Of particular interest is a 30- to 40-km-wide orogen-parallel belt of rapid exhumation that extends from central Nepal to the western Himalaya and its possible linkage to a midcrustal ramp in the basal decollement, and the related growth of Lesser Himalayan duplex structures. Here we present 26 new apatite fission track cooling ages from the Beas-Lahul region, at the transition from the Central to the Western Himalaya (77 degrees-78 degrees E) to investigate segmentation in the Himalayan arc from a thermochronologic perspective. Together with previously published data from this part of the orogen, we document significant lateral changes in exhumation between the Dhauladar Range to the west, the Beas-Lahul region, and the Sutlej area to the east of the study area. In contrast to the Himalayan front farther east, exhumation in the far western sectors is focused at the frontal parts of the mountain range and associated with the hanging wall of the Main Boundary Thrust fault ramp. Our results allow us to spatially correlate the termination of the rapid exhumation belt with a midcrustal ramp to the west. We suggest that a plunging anticline at the northwestern edge of the Larji-Kullu-Rampur window represents the termination of the Central Himalayan segment, which is related to the evolution of the Lesser Himalayan duplex. Key Points
We present an integrated study of LA-ICP-MS U-Pb age, Hf isotopes, and trace element geochemistry of zircons from the Himalayan eclogites (mafic rocks) and their host gneisses (felsic rocks) from the Kaghan Valley in Pakistan in order to understand the source and mode of their magmatic protoliths and the effect of metamorphism. Zircons from the so-called Group I (high-pressure) eclogites yielded U-Pb mean ages of 259 +/- 10 Ma (MSWD = 0.74), whereas those of Group II (ultrahigh-pressure) eclogites yielded 48 3 Ma (MSWD = 0.71). In felsic gneisses the central or core domains of zircons yielded ages similar to those from Group I edogites but zircon overgrowth domains yielded 47 +/- 1 Ma (MSWD = 1.9). Trace element data suggest a magmatic origin for Group I-derived (having Th/U ratios: >0.5) and metamorphic origin for Group II -derived (Th/U < 0.07) zircons, respectively. Zircon Hf isotope data, obtained from the same dated spots, show positive initial Hf-176/Hf-177 isotopic ratios referred to as "epsilon(Hf)(t)" of around +10 in Group I eclogites; +7 in Group II eclogites; and +8 in felsic gneisses zircons, respectively, thus indicate a juvenile mantle source for the protolith rocks (Panjal Traps) with almost no contribution from the ancient crustal material. The similar epsilon(Hf)(t) values, identical protolith ages and trace element compositions of zircons in felsic (granites or rhyolites) and mafic (basalt and dolerite) rocks attest to a bimodal magmatism accounting for the Panjal Traps during the Permian. Later, during India-Asia collision in Eocene times, both the felsic and mafic lithologies were subducted to mantle-depths (>90 km; coesite-stable) and experienced ultrahigh-pressure metamorphism before their final exhumation. (C) 2016 Elsevier B.V. All rights reserved.
Fold and thrust belts are characteristic features of collisional orogen that grow laterally through time by deforming the upper crust in response to stresses caused by convergence. The deformation propagation in the upper crust is accommodated by shortening along major folds and thrusts. The formation of these structures is influenced by the mechanical strength of décollements, basement architecture, presence of preexisting structures and taper of the wedge. These factors control not only the sequence of deformation but also cause differences in the structural style.
The Himalayan fold and thrust belt exhibits significant differences in the structural style from east to west. The external zone of the Himalayan fold and thrust belt, also called the Subhimalaya, has been extensively studied to understand the temporal development and differences in the structural style in Bhutan, Nepal and India; however, the Subhimalaya in Pakistan remains poorly studied. The Kohat and Potwar fold and thrust belts (herein called Kohat and Potwar) represent the Subhimalaya in Pakistan. The Main Boundary Thrust (MBT) marks the northern boundary of both Kohat and Potwar, showing that these belts are genetically linked to foreland-vergent deformation within the Himalayan orogen, despite the pronounced contrast in structural style. This contrast becomes more pronounced toward south, where the active strike-slip Kalabagh Fault Zone links with the Kohat and Potwar range fronts, known as the Surghar Range and the Salt Range, respectively. The Surghar and Salt Ranges developed above the Surghar Thrust (SGT) and Main Frontal Thrust (MFT). In order to understand the structural style and spatiotemporal development of the major structures in Kohat and Potwar, I have used structural modeling and low temperature thermochronolgy methods in this study. The structural modeling is based on construction of balanced cross-sections by integrating surface geology, seismic reflection profiles and well data. In order to constrain the timing and magnitude of exhumation, I used apatite (U-Th-Sm)/He (AHe) and apatite fission track (AFT) dating. The results obtained from both methods are combined to document the Paleozoic to Recent history of Kohat and Potwar.
The results of this research suggest two major events in the deformation history. The first major deformation event is related to Late Paleozoic rifting associated with the development of the Neo-Tethys Ocean. The second major deformation event is related to the Late Miocene to Pliocene development of the Himalayan fold and thrust belt in the Kohat and Potwar. The Late Paleozoic rifting is deciphered by inverse thermal modelling of detrital AFT and AHe ages from the Salt Range. The process of rifting in this area created normal faulting that resulted in the exhumation/erosion of Early to Middle Paleozoic strata, forming a major unconformity between Cambrian and Permian strata that is exposed today in the Salt Range. The normal faults formed in Late Paleozoic time played an important role in localizing the Miocene-Pliocene deformation in this area. The combination of structural reconstructions and thermochronologic data suggest that deformation initiated at 15±2 Ma on the SGT ramp in the southern part of Kohat. The early movement on the SGT accreted the foreland into the Kohat deforming wedge, forming the range front. The development of the MBT at 12±2 Ma formed the northern boundary of Kohat and Potwar. Deformation propagated south of the MBT in the Kohat on double décollements and in the Potwar on a single basal décollement. The double décollement in the Kohat adopted an active roof-thrust deformation style that resulted in the disharmonic structural style in the upper and lower parts of the stratigraphic section. Incremental shortening resulted in the development of duplexes in the subsurface between two décollements and imbrication above the roof thrust. Tectonic thickening caused by duplexes resulted in cooling and exhumation above the roof thrust by removal of a thick sequence of molasse strata. The structural modelling shows that the ramps on which duplexes formed in Kohat continue as tip lines of fault propagation folds in the Potwar. The absence of a double décollement in the Potwar resulted in the preservation of a thick sequence of molasse strata there. The temporal data suggest that deformation propagated in-sequence from ~ 8 to 3 Ma in the northern part of Kohat and Potwar; however, internal deformation in the Kohat was more intense, probably required for maintaining a critical taper after a significant load was removed above the upper décollement. In the southern part of Potwar, a steeper basement slope (β≥3°) and the presence of salt at the base of the stratigraphic section allowed for the complete preservation of the stratigraphic wedge, showcased by very little internal deformation. Activation of the MFT at ~4 Ma allowed the Salt Range to become the range front of the Potwar. The removal of a large amount of molasse strata above the MFT ramp enhanced the role of salt in shaping the structural style of the Salt Range and Kalabagh Fault Zone. Salt accumulation and migration resulted in the formation of normal faults in both areas. Salt migration in the Kalabagh fault zone has triggered out-of-sequence movement on ramps in the Kohat.
The amount of shortening calculated between the MBT and the SGT in Kohat is 75±5 km and between the MBT and the MFT in Potwar is 65±5 km. A comparable amount of shortening is accommodated in the Kohat and Potwar despite their different widths: 70 km Kohat and 150 km Potwar. In summary, this research suggests that deformation switched between different structures during the last ~15 Ma through different modes of fault propagation, resulting in different structural styles and the out-of-sequence development of Kohat and Potwar.