Stream Power Incision Model Research Articles

Quantitative analysis of the fluvial geomorphology and erosion on Hainan island: implications for the source-to-sink system in the NW South China Sea

The sediment delivery from Hainan island into the NW South China Sea during the Quaternary is less well-defined. An investigation into the uplift, exhumation and fluvial geomorphology of the Hainan island is crucial for improving our understanding on the source-to-sink system in this region. In this study, we employed the digital elevation analysis, the stream power incision model and the cosmogenic nuclide isotope analysis to unravel how and why the Hainan island provide sediment to the NW South China Sea. The results show that the average HI values of the main catchments on the Hainan island are below 0.35, which indicates that the rivers on the Hainan island tend to be stable. After mapping the channel steepness index of the Hainan island, a high channel steepness index is constrained in the central and western Hainan island. The χ analysis shows minor variations in χ values across the Hainan island, except of the watershed between the Wanquanhe and Nandujiang rivers, as well as that between the Changhuajiang and Nandujiang rivers, which indicates the potential migration of the river watersheds here. In addition, we carried out the cosmogenic nuclide 10Be analysis on the river sand from the Hainan island. The spatial distribution of 10Be concentrations is identified to be high in the west and low in the east. The catchment-averaged erosion rates are quantified to be 35 m/Myr, 42.5 m/Myr and 69.9 m/Myr for the Changhuajiang, Wanquanhe and Lingshuihe catchments, respectively. We also estimated the sediment deliveries from the Changhuajiang, Wanquanhe and Lingshuihe catchments to the Qiongdongnan and Yinggehai Basins to be 0.46 Mt/yr, 0.45 Mt/yr and 0.21 Mt/yr, respectively. We compared our results of the fluvial geomorphology analysis and erosion rate estimates of the Hainan island with those of the Taiwan island. We concluded that the difference of the erosion pattern and river evolution of the Hainan and Taiwan islands may be associated with the variation of the tectonics. Our investigation on the evolution of the Hainan island will improve our understanding on the source-to-sink systems in the NW South China Sea.

When Does the Concavity Index Constrain Stream Power Parameters?

AbstractBy defining the attributes of river networks, we can quantitatively extract records of climatic and tectonic changes from them. The stream power incision model (SPIM) provides a framework within which this can be achieved, as it facilitates the calculation of the relative rock uplift from river characteristics. One parameter that has been widely employed in tectonic and fluvial geomorphology is the channel steepness index, a metric that can represent the normalized rock uplift rate experienced by a river. However, to accurately infer the channel steepness index, we must accurately estimate m/n, the ratio between the two positive exponents of the SPIM. Present methodologies to constrain m/n rely on an assumption that rock uplift and erodibility are spatially invariant. These conditions are rarely present on Earth. In this study, we use a synthetic example and examples from the Siwalik Hills and Olympic Mountains to demonstrate how existing methodologies to constrain m/n produce systematic errors when there is spatial variation, and particularly spatial gradients, in the processes driving landscape evolution. To solve this problem, we present a methodology to estimate m/n based on a large river network inversion that accounts for spatial variation in landscapes. After demonstrating that the methodology can accurately recover m/n in our synthetic landscape, we show that our methodology can reconcile contrasting observations in the Siwaliks, and is critical to inferring accurate values of channel steepness index in the Olympic Mountains. This highlights the utility of large topographic inversions for investigating landscape dynamics.

Stochastic in Space and Time: 2. Effects of Simulating Orographic Gradients in Daily Runoff Variability on Bedrock River Incision

AbstractThe extent to which climate and tectonics can be coupled rests on the degree to which topography and erosion rates scales linearly. The stream power incision model (SPIM) is commonly used to interpret such relationships, but is limited in probing mechanisms. A promising modification to stream power models are stochastic‐threshold incision models (STIM) which incorporate both variability in discharge and a threshold to erosion. In this family of models, the form of the topography erosion rate relationship is largely controlled by runoff variability. Applications of STIM typically assume temporally variable, but spatially uniform and synchronous runoff generating events, an assumption that is likely broken in regions with complicated orography. To address this limitation, we develop a new 1D STIM model, which we refer to as spatial‐STIM. This modified version of STIM allows for stochasticity in both time and space and is driven by empirical relationships between topography and runoff statistics. Coupling between mean runoff and runoff variability via topography in spatial‐STIM generates highly nonlinear relationships between steady‐state topography and erosion rates. We find that whether the daily statistics of runoff are spatially linked or unlinked, which sets whether there is spatial synchronicity in the recurrence interval of runoff generating events, is a fundamental control on landscape evolution. Many empirical topography—erosion rate data sets are based on data that span across the endmember scenarios of linked versus unlinked behavior. It is thus questionable whether singular SPIM relationships fit to those data can be meaningfully related to their associated hydroclimatic conditions.

Geotechnical controls on erodibility in fluvial impact erosion

Abstract. Bedrock incision by rivers is commonly driven by the impacts of moving bedload particles. The speed of incision is modulated by rock properties, which is quantified within a parameter known as erodibility that scales the erosion rate to the erosive action of the flow. Although basic models for the geotechnical controls on rock erodibility have been suggested, large scatter and trends in the remaining relationships indicate that they are incompletely understood. Here, we conducted dedicated laboratory experiments measuring erodibility using erosion mills. In parallel, we measured uniaxial compressive strength, tensile strength, Young's modulus, bulk density, and the Poisson's ratio for the tested lithologies. We find that under the same flow conditions, erosion rates of samples from the same lithology can vary by a factor of up to 60. This indicates that rock properties that may vary over short distances within the same rock can exert a strong control on its erosional properties. The geotechnical properties of the tested lithologies are strongly cross-correlated, preventing a purely empirical determination of their controls on erodibility. The currently prevailing model predicts that erosion rates should scale linearly with Young's modulus and inversely with the square of the tensile strength. We extend this model using first-principle physical arguments, taking into account the geotechnical properties of the impactor. The extended model provides a better description of the data than the existing model. Yet, the fit is far from satisfactory. We suggest that the ratio of mineral grain size to the impactor diameter presents a strong control on erodibility that has not been quantified so far. We also discuss how our laboratory results upscale to real landscapes and long timescales. For both a revised stream power incision model and a sediment-flux-dependent incision model, we suggest that long-term erosion rates scale linearly with erodibility and that, within this theoretical framework, relative laboratory measurements of erodibility can be applied at the landscape scale.

Implicit Algorithm for Threshold Stream Power Incision Model

AbstractWe present an O(n) complexity and implicit algorithm for the two‐dimensional solution of the stream power incision model (SPIM) enriched by a discharge threshold term and taking into account variability in rainfall and thus discharge. The algorithm is based on the formulation developed by Deal et al. (2018, https://doi.org/10.1002/2017jf004393) and the generalization of the FastScape algorithm (Braun & Willett, 2013, https://doi.org/10.1016/j.geomorph.2012.10.008) where the slope is approximated by first‐order accurate finite difference. We consider a variety of discharge thresholds that vary in their dependence on channel slope. The algorithm requires finding the root of a non‐linear equation using a Newton‐Raphson iterative scheme. We show that the convergence of this scheme is unconditional, except for a narrow range of model parameters where the threshold increases with the slope and for low discharge variability. We also show that the rate of convergence of the iterative scheme is directly proportional to the slope exponent n in the SPIM. We compare the algorithm to analytical solutions and to numerical solutions obtained using a higher‐order finite difference scheme. We show that the accuracy of the FastScape algorithm and its generalization presented here is comparable to other schemes for values of n > 1. We also confirm that the FastScape algorithm and its generalization to variable discharge + threshold conditions does not need to satisfy the CFL condition and provides an accurate solution for both small and very long time steps. We finally use the new algorithm to quantify how the existence of an erosional threshold strongly affects the length of the post‐orogenic decay of mountain belts.

Morphotectonic assessment of the Gaula river basin, Kumaun lesser Himalaya, Uttarakhand

Morphotectonics is a discipline of geomorphology that explores how landforms are formed or modified by tectonic action. In the present study we uthe sed application of geospatial technology to evaluate the active deformation owing to Main Boundary Thrust (MBT) and Ramgarh Thrust (RT) in Outer part of Kumaun Lesser Himalaya. Active deformation inside river valleys is ascribed downcuttinging and incision of bedrock, offset of river channel, lateral migration of the river, and development of Knickpoints. Therefore, we used quantitative geo-mathematical analysis of The Gaula river basin using ASTER DEM along with survey of India Toposheet (1:50,000) supplemented by field observations. Four morphotectonic metrics, including the basin elongation ratio, the asymmetry factor, the channel sinuosity, and the hypsometric integral, were derived to examine the relative index of active tectonics (RIAT). We applied Stream Power Incision Model (SPIM) by analysing the steepness index (Ksn) using the stream power law. Additionally, a comprehensive evaluation of tributary streams of the Gaula watershed has been carried out to understand relative tectonic activity. Further we used Chi integral (χ) to measure transience dynamics and geometric disequilibrium of Gaula basin. We conclude that most sub-basins are tectonically enhanced, structurally organized and elongated in nature. The landform in the study area is largely controlled by presence of MBT, RT and subsidiary thrusts spread over in the vicinity. Based on our activity classes we concluded that approximately 31% of the area in Gaula river basin falls in very active whereas ∼34% is potentially active. Furthermore, we suggest that the results will be very useful for estimation of hazard potential zones in sub-Himalayan region of foothill zone of Kumaun Himalaya.

Plan‐Form Evolution of Drainage Basins in Response to Tectonic Changes: Insights From Experimental and Numerical Landscapes

AbstractSpatial gradients in rock uplift control the relief and slope distribution in uplifted terrains. Relief and slopes, in turn, promote channelization and fluvial incision. Consequently, the geometry of drainage basins is linked to the spatial pattern of uplift. When the uplift pattern changes, basin geometry is expected to change by migrating water divides. However, the relations between drainage pattern and changing uplift patterns remain elusive. The current study investigates the plan‐view evolution of drainage basins and the reorganization of drainage networks in response to changes in the spatial pattern of uplift, focusing on basin interactions that produce globally observed geometrical scaling relations. We combine landscape evolution experiments and simulations to explore a double‐stage scenario: the emergence of a fluvial network under block uplift conditions followed by tilting that forces drainage reorganization. We find that the globally observed basin spacing ratio and Hack's parameters emerge early in the basin formation and are maintained by differential basin growth. In response to the tilting, main divide migration induces basin size changes. However, basins' scaling relations are mostly preserved within a narrow range of values, assisted by incorporation and disconnection of basins to and from the migrating main divide. Lastly, owing to similarities in landscape dynamics and response rate to uplift pattern changes between experiments and simulations, we conclude that the stream power incision model can represent fluvial erosion processes operating in experimental settings.

Response of drainage to tectonics and PS-InSAR derived deformation study in Bilaspur, northwestern Himalaya, India

The study of drainage patterns in tectonically active regions is conducive to the prediction of regional geomorphology. Subtle subsurface changes can be detected by drainage conditions and manifested in the form of drainage anomalies. The Satluj valley of Bilaspur, which is traversed by numerous faults in northwest Himalayan region, was selected to analyze the effect of active tectonics on drainage evolution. With the Persistent Scatterer Interferometric Synthetic Aperture Radar (PS-InSAR) technique, SENTINEL-1A data were used to estimate the active surface deformation between September 2015 and December 2020. The results show that the region between Barasar Thrust (BrT) and Main Central Thrust (MCT) is undergoing deformation of ±12 mm/yr. The Stream Power Incision Model (SPIM) was used to predict deformation patterns. To validate the tectonic activity generated by the drainage network, seismic b-values were calculated, indicating the accumulating stresses. This study shows the importance of drainage anomalies in tectonically active areas. When used in close combination with other seismotectonic parameters, drainage anomalies can be effective in delineating tectonically active regions.

Delineation of tectonically active zones in the Island Belt Uplift region, Kachchh Basin, western India: A geomorphic and geodetic approach

In the present study, we have identified the tectonically active zones along the Island Belt Uplift (IBU), an upthrown block of E-W trending Island Belt Fault (IBF). Major and minor transverse faults have segmented the IBU into four sub-uplifts, viz. Pachham Uplift (PU), Khadir Uplift (KU), Bela Uplift (BU), and Chorar Uplift (CU) from west to east, respectively. The work includes geodetic and morphotectonic analysis of the IBU region. The Normalized Steepness Index (Ksn) analysis was performed supported by swath profiling to mark the tectonic activity variations along and across the strike of the IBF in the IBU region. Further, the Ksn contour map was prepared by using a power-law function of slope (S) of the longitudinal profiles of the river, concavity index (ө), and upstream area (A), which is an integral part of the Stream Power Incision Model (SPIM). The Ksn analysis suggests that the maximum amount of tectonic activity is restricted along the northern escarpment of the IBU, which falls in the vicinity of IBF and it's associated longitudinal and transverse faults. Based on the results of Ksn analysis, we have marked the tectonically active zones in the IBU region, which shows that amongst four sub-uplifts, the PU and BU were found most active (Zone-I), the KU expresses moderate activity (Zone-II), while CU demonstrated to be least active (Zone-III). The results were also confirmed by documentation of active tectono-geomorphic features like strath terraces, river offsets, straight fault scarps, triangular facets, linear valleys, fault ridges, Quaternary upwarps, laterally displaced crests and alluvial fans within the IBU. The active deformation pattern calculated using the Global Positioning System (GPS) data-set from 2009–2019 along IBU reveal maximum deformation rate of ∼2.5 ± 0.5 mm/yr in the PU, ∼2.0 ± 0.5 mm/yr in the BU while, the KU shows the deformation rate of ∼1.2 ± 0.5 mm/yr, which coincides well with the results of Ksn. The Ksn analysis results were also validated by comparing them with existing Gradient Length Anomaly (GLA) and Differential Interferometric Synthetic Aperture Radar (DInSAR) data. The ongoing tectonic activity in the study area is also manifested by the occurrence of moderate seismic activities in the identified active zones. The active tectono-geomorphic features and fault scarp geometry were analyzed and compared with the existing experimental models of basement fault and overlying sedimentary covers, which allow us to understand the nature of fault components of IBF. The analysis reveals that the IBF has been influenced by both dip-slip and strike-slip fault components of motion, which can be further confirmed by geophysical analysis. The results of Ksn analysis and field evidence were merged to prepare Active Fault Map of the IBU region that can be implemented for future seismic hazard assessment.

Upscaling Sediment‐Flux‐Dependent Fluvial Bedrock Incision to Long Timescales

AbstractFluvial bedrock incision is driven by the impact of moving bedload particles. Mechanistic, sediment‐flux‐dependent incision models have been proposed, but the stream power incision model (SPIM) is frequently used to model landscape evolution over large spatial and temporal scales. This disconnect between the mechanistic understanding of fluvial bedrock incision on the process scale, and the way it is modeled on long time scales presents one of the current challenges in quantitative geomorphology. Here, a mechanistic model of fluvial bedrock incision that is rooted in current process understanding is explicitly upscaled to long time scales by integrating over the distribution of discharge. The model predicts a channel long profile form equivalent to the one yielded by the SPIM, but explicitly resolves the effects of channel width, cross‐sectional shape, bedrock erodibility, and discharge variability. The channel long profile chiefly depends on the mechanics of bedload transport, rather than bedrock incision. In addition to the imposed boundary conditions specifying the upstream supply of water and sediment, and the incision rate, the model includes four free parameters, describing the at‐a‐station hydraulic geometry of channel width, the dependence of bedload transport capacity on channel width, the threshold discharge of bedload motion, and reach‐scale cover dynamics. For certain parameter combinations, no solutions exist. However, by adjusting the free parameters, one or several solutions can usually be found. The controls on and the feedbacks between the free parameters have so far been little studied, but may exert important controls on bedrock channel morphology and dynamics.

Quantifying Normal Fault Evolution from River Profile Analysis in the Northern Basin and Range Province, Southwest Montana, USA

Abstract Over the past few decades, tectonic geomorphology has been widely implemented to constrain spatial and temporal patterns of fault slip, especially where existing geologic or geodetic data are poor. We apply this practice along the eastern margin of Bull Mountain, Southwest Montana, where 15 transient channels are eroding into the flat, upstream relict landscape in response to an ongoing period of increased base level fall along the Western North Boulder fault. We aim to improve constraints on the spatial and temporal slip rates across the Western North Boulder fault zone by applying channel morphometrics, cosmogenic erosion rates, bedrock characteristics, and calibrated reproductions of the modern river profiles using a 1-dimensional stream power incision model that undergoes a change in the rate of base level fall. We perform over 104 base level fall simulations to explore a wide range of fault slip dynamics and stream power parameters. Our best fit simulations suggest that the Western North Boulder fault started as individual fault segments along the middle to southern regions of Bull Mountain that nucleated around 6.2 to 2.5 Ma, respectively. This was followed by the nucleation of fault segments in the northern region around 1.5 to 0.4 Ma. We recreate the evolution of the Western North Boulder fault to show that through time, these individual segments propagate at the fault tips and link together to span over 40 km, with a maximum slip of 462 m in the central portion of the fault. Fault slip rates range from 0.02 to 0.45 mm/yr along strike and are consistent with estimates for other active faults in the region. We find that the timing of fault initiation coincides well with the migration of the Yellowstone hotspot across the nearby Idaho-Montana border and thus attribute the initiation of extension to the crustal bulge from the migrating hotspot. Overall, we provide the first quantitative constraints on fault initiation and evolution of the Western North Boulder fault, perhaps the farthest north basin in the Northern Basin and Range province that such constraints exist. We show that river profiles are powerful tools for documenting the spatial and temporal patterns of normal fault evolution, especially where other geologic/geodetic methods are limited, proving to be a vital tool for accurate tectonic hazard assessments.

A New Data‐Driven Bayesian Inversion of Fluvial Topography Clarifies the Tectonic History of the Corinth Rift and Reveals a Channel Steepness Threshold

AbstractLandscape evolution models that invert topography for rock uplift can improve our understanding of both tectonic and geomorphic processes when properly constrained with data. Here we present a flexible, data‐driven Bayesian approach to invert fluvial topography for tectonic and geomorphic model parameters and apply it to a case study, the uplifting footwall of the Corinth Rift, Greece. We invert transient river profiles and up‐flexed marine terraces to resolve seven unknown parameters in a regional‐to‐flexural uplift tectonic model and the stream power incision model. The best‐fit tectonic parameters are consistent with independent data and predict block uplift rates of ∼0.1 mm yr−1 that changed to flexural uplift rates of ∼1.6 mm yr−1 at ∼0.6 Ma, as the master normal fault initiated. Similarly, the best‐fit geomorphic parameters predict sediment flux consistent with the offshore record and erodibility consistent with previous studies. However, the drainage area exponent, , of ∼2, and slope exponent, , of ∼7, are unusually high, indicating a threshold channel steepness where fluvial topography is largely insensitive to rock uplift rate >0.05 mm yr−1. Analysis indicates channels narrow to accommodate enhanced uplift rates, but channel narrowing only partially explains our results, suggesting that other processes not accounted for in the generic stream power model are also relevant to bedrock river incision in Corinth. Our results help clarify the tectonic and geomorphic evolution of the Corinth Rift, have important implications for studies that invert topography for rock uplift histories, and provide insight into potential limitations of some long‐term river incision models.

The Sliding Ice Incision Model: A New Approach to Understanding Glacial Landscape Evolution

AbstractAnalytical long‐profile models of fluvial incision are fundamental in geomorphology for building intuition about long‐term landscape evolution and interpreting the climatic and tectonic histories of landscapes. However, analogous models are missing for glaciers, despite the fact that glaciers are important and widespread agents of physical erosion. We introduce a new model of glacial erosion in terms of catchment area, ice accumulation rate, and ice surface slope in a form that is similar to the classic stream power incision model for fluvial erosion. We derive an analytical solution for coupled glacial‐fluvial networks under steady‐state climate and uplift, predict an ice surface slope–catchment area power‐law relationship, and predict the scaling of orogen relief with uplift rate. Finally, we note that this theory should enable glacial landscape evolution models that rival modern fluvial models in efficiency.

Temporal profiling of uplift rate along an active fault using river long profile in the Kuchchh region, Western India

River long profiles are widely used as tectonic event markers. Incision/uplift rate estimations based on the dating of river terrace or exposed scarp sediment give long-term averaged uplift rate values with the limitation of not giving the intermediate phases of highs and lows. On the other hand, GPS based estimates provide only decadal-scale uplift rate estimates. In our study, we adopted the long profile inversion methodology to obtain a temporal profile of uplift rate history (over the Late Pleistocene-Holocene time period) along an active fault (the Kuchchh Mainland Fault) in Kuchchh region, western India, which has experienced several large magnitude earthquakes in recent past. We also used long profiles analysis to quantify the transient response of Kuchchh landscape to uplift activity along the Kuchchh Mainland Fault. Steepness index and chi profiles were estimated to display the spatial variability in landscape disequilibrium. Integral method-based chi profiles were used to obtain the optimum value of basin area exponent of the Stream Power Incision Model (SPIM). Our estimates provide a mean uplift rate value of 1.1–1.4 mm/yr along the Kuchchh Mainland Fault during the Late Pleistocene-Holocene period, which matches with the average uplift rate derived from other independent observations in this region. Further, a chronology of major tectonic phases from the estimated uplift rate time series also corroborates with field-based estimates of major incision events in Kuchchh region. We also tested the robustness and sensitivity of our uplift rate estimates to model parameters. Such data can be complementary to the existing sparsely (both spatially and temporally) distributed uplift rate database, which is imperative to seismic hazard assessment.

Calibrating fluvial erosion laws and quantifying river response to faulting in Sardinia, Italy

It is now widely accepted that rivers modify their erosion rates in response to variable rock uplift rates, resulting in changes in channel slope that propagate upstream through time. Therefore, present-day river morphology may contain a record of tectonic history. The simple stream power incision model can, in principle, be used to quantify past uplift rates over a variety of spatial and temporal scales. Nonetheless, the erosional model's exponents of area and slope (m and n respectively) and ‘bedrock erodibility’ (k) remain poorly constrained. In this paper, we will use a geologically and geomorphically well constrained Plio-Pleistocene volcanic landscape in central Sardinia, Italy, to calibrate the stream power erosion equation and to investigate the slip rate of faults that have been seismically quiescent in the historic past. By analysing digital elevation models, geological maps and Landsat imagery, we have identified the geomorphic expression of several volcanic features (eruption centres and basaltic lava flows) and three normal faults with 6 to 8 km fault traces within the outcrop. Downstream, river longitudinal profiles show a similar transient response to relative base level fall, probably as a result of relief inversion at the edge of the volcanic outcrop. From measurements of incision, local slope and upstream catchment area across eight different rivers, we calculate n ≈ 1, m = 0.50 ± 0.02 and, using a landscape age from literature of 2.7 Ma, bedrock erodibility k = 0.10 ± 0.04 m(1−2m) Myr−1. There are also knickpoints on rivers upstream of two normal faults, and we used numerical inverse modelling of the longitudinal profiles to predict the slip rate of these faults since 2.7 Ma. The results from the inverse model show that the erosional parameter values derived in this study can produce theoretical longitudinal profiles that closely resemble observed river profiles upstream of the faults. The lowest misfit theoretical longitudinal profiles were generated by a model of temporally discontinuous footwall uplift with consistently low throw rates (<0.1 mm yr−1). The predicted footwall uplift history is similar for the two faults, both showing periods of fault slip and no fault movement since 2.7 Ma.

GIS and DEM based analysis of incision and drainage reorganization of the Buyuan River basin in the upper Lancang-Mekong of China since the Late Pleistocene

River incision and drainage reorganization have an important impact on the site selection of many major projects including city, road and others, and are the key issues of Quaternary environmental changes. Studies of river incision and river-network adjustment have traditionally been based on extensive field evidence, such as sediment age and beheaded river system. The Buyuan River basin is a large sub-basin of the upper Lancang-Mekong, with high mountains and extremely active erosion. The latter affects the preservation of the Quaternary period sediments leading to difficulties in understanding the main evolution characteristics of the basin. This study investigates differences in the equilibrium state of the longitudinal profile, infers incision rates, and evaluates drainage divide migration timelines using the stream-power incision model, the latest morphological dating, and Chi-plots (χ-z) based on digital elevation models (DEMs) on the GIS software platform. The final results show that two significant erosion base-level decreases occurred in the Late Pleistocene at least. The incision rate of the mainstream might have been 0–2.99 mm/yr since 100 ka BP and 0–3.28 mm/yr since 46 ka BP. The Chi-values across the divides suggest that space limited (or constrained) river reorganization and that there is no severe reorganization in the basin; the imbalance of traceable erosion only exists in local areas. The main driving force for the geomorphologic evolution of the Buyuan River basin is likely climate fluctuations rather than strong tectonic uplift since the Late Pleistocene.

Drainage conformation and transient response of river system in thrust segmentation of Northwest Himachal Himalaya, India

In the northwest Himalayan region, the ongoing out-of-sequence deformation caused emergence of Back Thrust (BT) around 43ka, which attributed the growth of anticlines parallel to the strike of major thrusts. This spatial and temporal geomorphic deformation is accompanied by coexisting river incision in the axial zone of anticlines and base-level fall along with upstream migration of knickpoints in the longitudinal river course. We used the Stream Power Incision model (SPIM) which relates the normalized steepness index (ksn), and concavity index (θref) with the slope (S) of longitudinal river profile as a power-law function of upstream area (A). and estimated these parameters for Beas, Soan, and Sutlej rivers to decipher the spatial distribution of deformation pattern. Using SRTM DEM and based on strike parallel variations of the ksn we identified segmented nature of the NW Himalayan thrusts. Further, the dynamics of river basins and transience along the tributaries and main stem rivers has been analyzed by Chi analysis (χ), which is attributed to geometric disequilibrium of the rivers. The results of ksn and χ analysis together suggests that the crustal shortening of NW Himalaya is caused by segmental reactivation of the thrusts. Further, investigating the change in steepness along longitudinal course aids in understanding the differential uplifts patterns in actively deforming segments of major NW Himalayan thrusts.

Inversion of Dadu River Bedrock Channels for the Late Cenozoic Uplift History of the Eastern Tibetan Plateau

AbstractThe Dadu River in the eastern Tibetan Plateau has the potential to record transient variations of tectonic uplift since the Late Cenozoic. In this study, we extracted bedrock channels of the Dadu River drainage basin from the 30 m ASTER digital elevation model. The longitudinal river profiles were then analyzed, and three generations of knickpoints were recognized via the slope area method. By solving the analytical solution of the linear transient stream power incision model, we determined a Late Cenozoic history of regional tectonic uplift. The results reveal slow pre–middle Miocene uplift rates (&lt;0.1 mm/a) and then increased gradually to 0.33 ± 0.03 mm/a since 12–15 Ma until present. This increase of uplift rate from bedrock channel inversion is consistent with the regional findings of low‐temperature themochronology, emphasizing the potential and significance to utilize channel profiles to infer the tectonic uplift history.

Denudation Rate Changes Along a Fast‐Eroding Mountainous River With Slate Headwaters in Taiwan From 10Be (Meteoric)/9Be Ratios

AbstractThe novel 10Be (meteoric)/9Be system, where 10Be is delivered by precipitation and stable 9Be is released by weathering, provides denudation rates over weathering‐erosion timescales. The new tool is applicable to quartz‐poor lithologies, for example, mafic rock and claystone, which are not readily accessible by the commonly used in situ‐produced 10Be in quartz. We provide a first application of this proxy to a tectonically active mountainous river, the Zhuoshui River in Taiwan. Taiwan Rivers supply a disproportionately high suspended and dissolved flux to the oceans and are often underlain by fine‐grained shale/slate. 10Be (meteoric)/9Be‐derived denudation rates (Dmet) from the Zhuoshui Catchment are highest in the slate‐dominated headwaters (4–8 mm/year), and much lower (1–2 mm/year) along the midlower reaches with mixed lithologies. At the basin‐wide scale, we find a poor correlation between Dmet and basin‐averaged channel steepness despite a small climatic gradient. Because large lithological heterogeneities exist in this basin, we invoke a lithological effect to explain this poor correlation. Relying on a revised stream power incision model that incorporates rock erodibility, the resulting lithology‐ and runoff‐adjusted ksn (kLrsn) can be reconciled with denudation rates with the highest erodibility predicted to prevail in the Miocene slate of low metamorphic grade and high fracture density. This model suggests that the lithological heterogeneity can alter the coupling between surface denudation and channel morphology. On a broader perspective, the successful application of the 10Be (meteoric)/9Be proxy shows its applicability as a tracer for erosion and sediment transport processes in fast‐eroding mountain belts underlain by slate lithologies.

Struggles with stream power: Connecting theory across scales

Two empirically consistent explanations exist for the morphology and dynamics of bedrock rivers. The stream power law appears to explain rates of rock incision at large scale (large drainage basins or landscape scale), but appears to be of limited value in predicting reach-scale morphology and dynamics. The underlying reason why stream power is such an effective tool for understanding landscape evolution is because it is an explicit representation of the inverse power relation between slope and area in rivers incising bedrock. This morphometric relation is nearly universal, subject to exogenous contingencies. Independent testing and validation of the stream power incision model is challenging because of the contingencies that exist in almost all landscapes. Mechanistic models of erosion have been tested and validated at small scales, generally below the reach scale, and opportunities exist for independent testing and validation by using them to predict reach-scale morphologies. Predicting larger scale bedrock river morphology requires the assumption that flow, sediment fluxes and bedrock erosion processes that occur at smaller scale are adequately averaged at larger scales, which is not well supported. Reconciling these two scales of explanation is one of the grand challenges of fluvial geomorphology, unless we accept that theory and explanations will be scale-delimited.