Original Research ARTICLE
How Continuous Are the “Relict” Landscapes of Southeastern Tibet?
- 1Department of Earth Sciences, University College London, London, United Kingdom
- 2Department Earth and Planetary Sciences, Birkbeck College, London, United Kingdom
- 3School of Earth Science and Resources, and Research Center for Tibetan Plateau Geology, China University of Geosciences, Beijing, China
Pervasive low-relief, high-elevation surfaces separated by incised canyons are common across the Southeastern margin of the Tibetan Plateau and have been used to define the nature of crustal deformation that drove plateau growth. A common assumption is that these surfaces were once part of a continuous low-relief paleotopography that has undergone surface uplift and dissection. Recent research, however, has questioned this assumption and the derived geodynamic models, which suggests that these surfaces formed in situ through drainage network reorganization and the piracy of upstream drainage area. Here, we test the continuity of the low-relief surfaces across SE Tibet using a new inversion scheme that also illuminates the nature of conflicting hypotheses. Our analysis is based on combining the local information contained in maps of normalized channel steepness with the more distributed and integrated information contained in maps of normalized landscape response time. This allows us to model the formation of a hypothetical landscape prior to rock uplift and dissection. We find that large variations in channel steepness are required along the trunk channels within the inferred paleotopography. This is inconsistent with a low-relief surface prior to surface uplift and indicates that a surface interpolated between remnants cannot be used to robustly measure geodynamic processes in space and time. Furthermore, our inverse framework highlights many different solutions to this ill-posed problem and thus provides an explanation as to why the topography alone cannot be used to provide a unique solution to the debate.
The South Eastern margin of the Tibetan Plateau has grown during the last 50 million years (Myr) through a combination of geodynamic processes: eastwards extrusion driven by motion on strike slip faults (Tapponnier et al., 1982); crustal thickening induced by lower crustal flow (Clark and Royden, 2000; Copley, 2008); slab rollback below its eastern margin (Sternai et al., 2012); and general shortening and crustal thickening due to collision of the Indian plate and Eurasia (Yang et al., 2015; Tan et al., 2019). In order to assess how these mechanisms have contributed to the current topography, it is crucial to determine rock uplift in space and time, commonly, through geochronology and geomorphology. One of the key geomorphic features that has been used to aid in this mapping are pervasive low-relief, high-elevation topography that is separated by incised canyons with high relief. These low-relief surfaces, that have areas typically 10 s of km2, have been identified across SE Tibet and decrease gradually in elevation from ∼5 km to less than 1 km over a lengthscale of approximately 1,000 km (Clark et al., 2006). The relief averaged over 5 km within these low-relief landscapes is typically less than 1 km. This landscape has been argued to result from an increase in rock uplift rate (Clark and Royden, 2000). Rivers respond to rock uplift by steepening and these steepened sections incise backwards into slowly eroding areas that have not yet equilibrated to the new rock uplift rate. These patches of low-relief areas are therefore assumed to be parts of the landscape that have not yet adjusted to the new rock uplift rate. The deeply incised canyons, however, have responded to this change in rock uplift rate.
By mapping the extent of these uplifted surfaces, the region that experienced a change in rock uplift can be estimated. In addition, if the difference between erosion rates on the surfaces and in the canyons can be estimated, the time required for these surfaces to reach their current elevation above the surrounding canyons can be inferred.
For example, Clark and Royden (2000) proposed that a low-viscosity channel in the lower crust accommodates the movement of crustal material from high pressures below the plateau to lower pressures away from the plateau. This mechanism simplifies to a diffusion equation in which elevation diffuses from the plateau leading to spatial and temporal variations in surface uplift (England and Molnar, 1990). Because crustal thickening occurs due to flow in the lower crust, the surface is uplifted but deformation is small and there is limited lateral transport of the surfaces. If the erosion rates on the low-relief surfaces are assumed to be close to zero, mapping these surfaces and dating incision can constrain the uplift rate history and the effective viscosity of the low-viscosity channel (Clark and Royden, 2000; Schoenbohm et al., 2004; Clark et al., 2006; Ouimet et al., 2010; Liu-Zeng et al., 2018). Furthermore, this scenario leads to an overall increase in topographic elevation through time. Throughout this study, we will refer to the idea that the low-relief surfaces are remnants of a relict low-relief landscape that has been dissected with minimal drainage network reorganization as Hypothesis 1 and stress that this is a simplification of the proposed model (Figure 1).
FIGURE 1. Hypothesis 1. (A) Before the onset of rock uplift across SE Tibet, Hypothesis 1 proposes that there was a low-relief landscape with low erosion rates. (B) Due to lower-crustal channel flow, rock uplift rates increased and rivers incised into the low-relief landscape. (C) Remnants of the low-relief landscape are preserved today and have not been removed by erosion. Therefore, the low-relief remnants can potentially be used to map geodynamic processes.
Importantly, the assumption that enables low-relief high-elevations surfaces to map rock uplift rates is that erosion rates on the surfaces were low prior to the change in rock uplift and have not changed since this adjustment, or at least that changes in erosion rate can be quantified. Several factors, however, may lead to temporal changes in erosion rate across low-relief surfaces as they are uplifted and some of these factors may even form these surfaces: the surfaces may experience changing elevation-dependent climatic conditions as they increase in elevation leading to changes in vegetation and changes in weathering processes that will influence erosion rates (Hales and Roering, 2007; Roering et al., 2010; Schaller et al., 2018); climate may change regionally due to the changes in topography resulting from surface uplift leading to changes in the patterns, in space and time, of rainfall (Molnar et al., 2010; Ferrier, et al., 2013; Scherler et al., 2017); Quaternary global climate change may have resulted in enhanced or reduced erosion across some of these surfaces (Egholm et al., 2009; Fox et al., 2015a; Egholm et al., 2017).
Recently, however, a model has been proposed in which the surfaces are formed by changes in erosion rate but no significant change in rock uplift rate driven by dynamic adjustment of the drainage networks due to tectonic strain (Hallet and Molnar, 2001; Yang et al., 2015). In this scenario, low slope rivers with large upstream drainage area and equilibrated channel steepness, lose upstream drainage area, reducing channel steepness and thus erosion rate. Once this erosion rate is reduced, the landscape is out of equilibrium with the regional rock uplift rate and would be advected upwards, while surrounding channels continue to erode at the rock uplift rate. As this low-relief surface is being uplifted, rivers incise into the boundaries of these isolated surfaces. In this scenario, the general features of the topography may stay relatively constant, with new low-relief surfaces forming and being lost to erosion through time. Throughout this study we will refer to this model as Hypothesis 2 (Figure 2). This hypothesis remains controversial and it is currently unclear whether the low-relief surfaces are relicts of a remnant topography and can be used to map rock uplift rates in space and time. In particular, Whipple et al. (2017) argued that the low-relief upland landscape patches are approximately co-planar and decrease in elevation from the northwest to the southeast and that variability in the elevations of bounding knickpoints is expected in natural landscapes. It remains unclear how much variability is expected and how much is due to regional geodynamics or the expected change in elevation along the trunk river profiles. Willett (2017) argued that Whipple et al. (2017a) only analyzed topography at the 100 –1,000 km scale but did not look at individual drainage basins. Furthermore, Willett (2017) suggested that the topography might be consistent with incision into a pre-existing landscape if spatially variable rock uplift is accounted for “pending a detailed study”. In response, Whipple et al. (2017b) agreed that diagnostic criteria in landscape evolution are not easily determined and that it is important to examine multiple criteria in a regional context. Whipple et al. (2017b) also maintained that the low-relief topography was co-planar and defines an upper envelope to the topography and thus the topography of SE Tibet does result from incision into a pre-existing low-relief landscape. Importantly, these studies use a stream-power framework and we will adopt these same assumptions here. For simplicity, we have distilled the two concepts into endmember hypotheses. In reality, it is important to note that neither hypothesis is mutually exclusive nor collectively exhaustive (Whipple et al., 2017a; Whipple et al., 2017b; Willett, 2017).
FIGURE 2. Hypothesis 2. (A) After a topography had developed across SE Tibet, Hypothesis 2 proposes there would have been some low-relief landscapes due to drainage reorganization. (B) During ongoing drainage reorganization driven by shortening perpendicular to major rivers, the tributary circled in (A) captures the upper river. The associated loss in upstream drainage area reduces the erosion rate of the river leading to surface uplift. Hillslope processes continue to erode hillslopes however the river has lost the ability to incise, leading to a reduction in topographic relief. Low-relief surfaces from previous capture events are eaten away by erosion. (C) The basin that lost upstream drainage area continues to have less and less relief and is eaten away by erosion. A capture event upstream of the red circle in (B) leads to the beginning of new low-relief topography. This process is continuous and leads to the lowest-relief landscape being preserved for short amounts of time at the highest elevations. Therefore, the modern topography reveals a snapshot of this process and patches of low-relief landscape at high elevations are not remnant of a pervasive low-relief landscape.
To test these hypotheses and to help diagnose how low-relief, high-elevation topography may form, we utilize linear inverse methods to interpolate a hypothetical topographical surface between the low-relief high-elevation surfaces (Fox, 2019). Importantly, we are not simply testing how co-planar the low-relief high-elevation surfaces are, but ensuring that the topographical surface must follow slope-area relationships with variable channel steepness. Furthermore, we account for the fact that the topographic surface might have been warped by geodynamic processes or offset by brittle faulting. If a topographical surface can be found with “relatively uniform” normalized rock uplift rate, u*, this supports the hypothesis that SE Tibet was a low-relief low-elevation surface that has been uplifted and dissected. Here “relatively uniform” would mean that the inferred normalized rock uplift rate between the low-relief surfaces is similar to the normalized rock uplift rate values observed on the low-relief surfaces. Alternatively, if a surface with large variations in normalized rock uplift rate is found, even without the inclusion of the incised portions of the landscape, this is inconsistent with Hypothesis 1 and instead supports Hypothesis 2. In this way, our approach does not explicitly test Hypothesis 2, but if a low-relief paleotopography can be constructed, the importance of Hypothesis 2 is reduced.
Using Fluvial Metrics to Test the Hypotheses
Fluvial geomorphology, and the associated metrics based on the stream power model, provide insight into how landscapes respond to changes in tectonic forcing (Howard, 1994). The rate of change of river long-profiles, dz/dt, increases with rock uplift rate, u, and decreases with erosion rate, e. In the stream power model, erosion rate is proportional to the upstream drainage area, A and the local channel slope, S, raised to the powers of m and n respectively. The constant of proportionality is the erodibility, K, which encompasses bedrock strength, bedload, hydraulic parameters and climate. Therefore,
and, if dz/dt = 0, u = KAmSn. Because K cannot be inferred directly from a landscape, it is common to measure the normalized channel steepness ksn=Am/nS, where m/n = 0.3–0.8 (Mudd et al., 2018) directly from a digital elevation model (DEM) which provides an estimate of the uplift rate at steady state (Kirby and Whipple, 2012). We build on previous work that has used the stream power model to debate the origin of the topography of SE Tibet and assume n = 1 and discuss the implications of this assumption in Dicussion (Yang et al., 2015; Whipple et al. 2017a; Whipple et al., 2017b; Willett, 2017). A value of m = 0.45 has been determined for parts of SE Tibet and we will use this same value (Kirby et al., 2003; Ouimet et al., 2009; Yang et al., 2015). The normalized channel steepness index approach to asses the topography of SE Tibet assumes steady state rock uplift, a fixed drainage network and a uniform K value in space and time. An advantage of this approach is that normalized channel steepness values can be calculated for every river node in a large digital elevation model, and it is therefore easy to interpret spatial patterns in the data (Kirby et al., 2003). In this respect, this approach is ideal to exploit the large topographic datasets available. A potential limitation is that by calculating slope from a DEM, noise in the dataset is amplified resulting in maps that can be hard to interpret. Therefore, averaging and smoothing is required to interpret normalized channel steepness maps. If low-relief surfaces were part of a continuous and former steady-state surface developed under relatively spatially constant rock uplift rate, similar channel steepness values would be expected across all low-relief surfaces and this has been argued by Whipple et al. (2017a). This would support Hypothesis 1. It is worth noting that Whipple et al. (2017b) do expect some variability in the form of the low-relief surface due to drainage re-arrangement and expected spatial variations in initial conditions, forcing mechanisms, and landscape response.
If the rock uplift rate increases, a fluvial knickpoint will form at the baselevel and propagate upstream at a speed given by KAm, if n = 1. Therefore, the response time at a specific location, x, in the landscape from the baselevel, xb, is given by
As with the normalized channel steepness, K is often unknown and thus it is useful to use a modified response time, or
To provide information on relative channel steepness values,
and the normalized rock uplift rate is u*= u/A0mK, which is also proportional to the normalized channel steepness, i.e., u* = ksn/A0m. We discuss the implications of the n = 1 and highlight how we can account for nonlinearity in the discussion. Furthermore, rivers that have experienced similar rock uplift rate histories should have similar forms as knickpoints will have traveled to similar
Our approach is based on the normalized analytical, steady state stream power model but allows for spatial variations in u* and surface uplift (S.U.) after the formation of the low-relief landscape. In this way, we adopt the same sets of assumptions that are the basis for the ongoing debate about the origin of the topography of SE Tibet and attempt to highlight why there is ongoing debate. If u* and S.U. vary in space, we can write a discrete version of Eq. 1 for a node in a DEM along a channel,
where the ith pixel is upstream of the jth pixel and the lowest most pixel has an elevation of B.L. (baselevel). S.U.i is the surface uplift that the low-relief landscape has experienced following dissection. If S.U. is constant across all of the low-relief landscape, this represents a block uplift scenario but it can also be spatially variable reflecting processes that might lead to surface uplift of the low-relief surface. Long wavelength features might be associated with lower crustal channel flow while shorter features might be driven by local faulting. Importantly, S.U. represents the total amount of surface uplift averaged since the low-relief landscape formed and this surface uplift might have occurred at any point in time. In turn, the elevation of any fluvial node can be predicted using a summation of expressions for downstream pixels, and any u* pattern and S.U. pattern. For a single channel, the ith fluvial node is the sum of differences in
Changes in the u* values along the trunk streams will change the elevations upstream of that specific pixel. If tributaries are utilized in the analysis, these changes will propagate across all the upstream river network. Therefore, the response of the landscape to local changes in S.U. and u* is distinct: increasing u* leads to small increases in elevation of rivers with large upstream drainage areas due to the small changes in
Maps of u* and surface uplift can be determined using inverse methods. For a single channel, Eq. 5 shows that i+1 model parameters are used to describe the single elevation node resulting in an ill-posed inverse problem. However, the node downstream constrains i-1 of the same u* model parameters and thus only an additional two parameters are required to describe the elevation of the ith node: one u* parameter and one S.U. parameter. This is still an ill-posed inverse problem which is exacerbated by the fact that we only use nodes from the low-relief surfaces to constrain the model parameters. It is the branching network of river channels that provides redundant information and smoothness constraints that enables maps of u* to be inferred from limited elevation pixels (Sternai et al., 2012). This is achieved by simplifying the drainage network to reduce the number of nodes, by discretizing space into blocks of constant u* and S.U. values that are 25 km by 25 km large and by introducing smoothness constraints on u* and S.U. (Fox, 2019). This smoothness is important where the data do not resolve the parameters, but rough maps can be produced if required by the data.
The aim of the inversion is to find a topography that represents the pre-incision landscape plus surface uplift and this requires finding the model parameters (u* and S.U. values) that minimize the misfit between predicted and observed elevations of low-relief surfaces, and the roughness of the u* and S.U. maps. Some form of regularization is required for most interpolation algorithms. For example, Clark et al. (2006) used a spline to interpolate between low-relief landscapes with a tension parameter determining how smooth the surface is. Here we use negative Laplacian operators to quantify smoothness. Weighting terms, α and λ, are used to determine how to minimize smoothness constraints compared to the fit to the low-relief surfaces. This results in a linear system of equations:
where G is the 2npixel × nnodes forward model as defined in Eq. 5. Each row of G contains differences in
The dataset used for the inversion is therefore just the elevations and
FIGURE 3. The dataset used for the analysis. (A) The HydroSheds hydrologically conditioned topographic dataset with a resolution of ∼90 m (Lehner et al., 2008). (B) Drainage network topology, upstream drainage area and values of were calculated from this dataset with a baselevel set at 500 m for pixels with an upstream drainage area >5 km2, with an m value of 0.45, and n value of 1 and a scaling area A0 of 1 m2 after (Yang et al., 2015). The box outlines the extent of Figure 4.
FIGURE 4. 4Subsampling the data. (A) The complete dataset was trimmed by randomly removing 90% of the pixels outside of the low-relief surfaces and 75% within the low-relief surfaces. (B) Local normalized channel steepness maps values highlight regional patterns but are very noisy and thus must be filtered or averaged to make interpretations. (C) Zoomed in portion showing data distribution and variability in the normalized steepness data.
The results are presented as two maps, u* and S.U. and the corresponding misfit value (Figure 5). Here we focus on one result out of several models produced with different values of α and λ. We show the other results in Supplementary Movie S1 and this also provides a means to visualize the sensitivity of the results to the smoothing parameters. In general, small values of the weighting parameters lead to rough models that fit the data well, whereas large weighting parameters produce smoother results that may fit the data poorly (Figure 5A; Supplementary Movie S1). With large values for the weighting parameters, the inversion attempts to minimize the misfit associated with the smoothness constraints (that the Laplacian of the u* values and the S.U. values should be close to zero) resulting in smooth maps. With small values for the weighting parameters, the inversion finds results that fit the data well but are more sensitive to anomalous topography and geomorphic noise, and therefore, results may be meaningless. Within this context, geomorphic noise encapsulates artifacts in the DEM, landslides blocking rivers and over-steepening rivers and small-scale variations in erodibility. It is expected that this noise far exceeds the uncertainty of the DEM (Fox et al., 2015b).
FIGURE 5. Predicted normalized rock uplift rates and spatially variable surface uplift. (A) A two-dimensional L-surface showing the trade-off between model misfit and the two damping parameters which control the roughness of the normalized rock uplift (α) and the surface uplift map (λ). However, each damping parameters influences both maps as the u* and S.U. both predict the topography and control the misfit between the observed and predicted topography, which is minimized by the inversion. A model is preferred that has a small misfit value, but also results in smooth maps. This model is shown by the circle in (A) and the maps are presented in (B) and (C). The sensitivity of the inversion to the damping parameters is presented in Supplementary Movie S1. Low u* values are observed where low-relief surfaces are found but higher values are observed along the trunk streams and the confluence of the Jinsha river with the Niulan River (NR) and aspects of the drainage network, such as large parallel rivers in the Three Rivers area (TR) and the large bends in the Jinsha River (JR), and the Yalong River (YR). High values are also observed close to the Dadu River (DR) which is related to the steep Longmen Shan. These anomalies are more apparent in the rougher models.
In order to choose a preferred model, we employ the same principles behind the L-curve (Hansen, 1992) and search for the damping variables that result in a model with low model roughness and low misfit (i.e., close to the corner of a 2-D L-surface, Richards et al., 2016). A misfit value of approximately 100 m is reasonable: misfit values lower than this suggest that we are fitting noise and misfit values higher than this indicate that the model is not fitting the data. We focus our discussion of the results on features that are robust across a range of values for the damping parameters.
Results show that there are clear short wavelength variations in normalized rock uplift rate and that high values are observed at large confluences and large bends in trunk rivers. By contrast, if the network geometry and values were consistent with a landscape that is being dissected by a wave of incision due to increased rock uplift, a map of relatively uniform u* values less than approximately 200 would be expected. This is because values less than 200 are observed across the low-relief surfaces. We refer to areas with higher than expected u* values as anomalies and these indicate that a topographic surface with low u*, or a low-relief paleotopography, cannot be interpolated between low-relief surfaces and that considerable drainage network reorganization has occurred. In particular, we see high anomalies at the prominent bends in two major tributaries to the Yangtze, the Jinsha and Yalong Rivers, at a confluence with the Niuling River and in the upper reaches within the Dadu (Yang et al., 2020; Suhail et al., 2020) catchment (Figure 5B). The timings of the capture events that formed these anomalies are debated, with the First Bend of the Yangtze either forming in the Eocene (Zheng et al., 2020) or the Pleistocene (Deng et al., 2020). High anomalies also follow these rivers, along with the Lancang, up into Tibet. These anomalies are most pronounced when maps of S.U. are least variable (Supplementary Movie S1) and this is because changes in elevation can be produced with changes in S.U. or changes in u*. Therefore, decreasing changes in S.U. leads to increased variability in u*. A high anomaly is also observed in the Three Gorges area. Crucially, the elevations of the large rivers, where unusual drainage has been previously mapped (Clark et al., 2004), are not used in the inversion. The high normalized channel steepness values that we predict are the result of rivers from low-relief surfaces flowing into the trunk stream at elevations that require anomalously steep trunk river segments. High values are also observed at the eastern end of the Nantinghe Fault (Figure 5B).
Across many of the low-relief surfaces, the residuals are coupled (Figure 6A): large positive values are associated with large negative values, which is caused by drainage divide migration (Fox et al., 2014). Importantly, these residuals highlight that the model is not fitting all channels that are over-steepened and under-steepened by drainage divide migration and capture. To explain all these residuals with the model, the normalized channel steepness map would be even rougher, providing further support for drainage network reorganization.
FIGURE 6. Investigating the results of the inversion. (A) The predicted - observed residuals for the low-relief surfaces show that in areas where the misfit is large and positive there are often negative residuals close by. This suggests that drainage divides and river capture are causing pairs of over- and under- steepened
We see evidence that active structures across SE Tibet are influencing our inferred values of S.U. (Figures 6B,C), and this correspondence highlights that short wavelength structures are also contributing to deformation of the low-relief surface. Furthermore, this suggests that damping parameters used in the inversion are not obscuring expected features. In particular, there is a clear change in S.U. across the Lijiang Fault, which is also an area of high strain rates (Copley, 2008; Pan and Shen, 2017; Li et al., 2019). Low values are observed to the east of the Lijiang fault where the Xianschuihe fault splits and to the west of the Lijiang fault, and these areas correspond to areas of GPS-resolved subsidence (Pan and Shen, 2017).
Our predicted elevations can also be displayed along the paths of rivers (Figures 7B,C). This provides an indication of how much incision has occurred at a specific location, if Hypothesis 1 is correct. It is important to point out that here we use incision to indicate the amount of deepening of valleys with respect to the surrounding peaks. In this respect incision is not equal to total erosion. If the erosion rates on the low-relief landscape are assumed to be close to zero, however, this incision is equal to the amount of erosion and could be measured with thermochronometry. The predicted elevations are a function of the pre-incision topography, itself, a result of the
FIGURE 7. Comparison with river profiles. (A) Location map showing the rivers plotted in (B) and (C); (B) River long-profiles extracted from the down sampled DEM, black lines. The crosses are the low-relief nodes used in the inversion, these appear as thick black lines where they are closely spaced. The pink lines are the predicted elevations and account for spatial variability in u* and S.U. and are therefore not monotonic. The stepped nature is due to the fine resolution of the pixels used for the inversion and the fact that the values of u* and S.U. vary over short wavelengths; (C) River long-profiles extracted from Mekong and a tributary of the Salween. Larger amounts of incision are predicted for these rivers.
FIGURE 8. Predicted incision. The amount of predicted incision across SE Tibet assuming that the low-relief topographical surface represents a pre-incision landscape. Incision is focused along the rivers in the Three Rivers region, as would be expected from a simple topographic interpolation. Several plotting artifacts are highlighted with red circles and these reflect areas that are unconstrained by the dataset.
We have chosen to discretize the tectonic history into two stages. The first stage represents a steady state landscape before the onset of enhanced rock uplift. The second stage represents the cumulative uplift and incision of this landscape and does not make any assumptions about how this varies in time. However, in Hypothesis 1 channel flow would lead to spatial and temporal variations in surface uplift (Clark and Royden, 2000; Schoenbohm et al., 2004). A full inversion of the landscape using river profile modeling with rock uplift rate variable in space and time (Roberts et al., 2012; Fox et al., 2014; Rudge et al., 2015) would require inferring maps of u* at many times in the past. Given that one of our goals is to investigate whether the planform geometry of the drainage network has evolved through time by using a simplified interpolation method, a full inversion that requires assuming this is not justified. Nevertheless, it is important to consider how our simplifications may impact our results. If the areas closest to the plateau were uplifted first, fluvial incision would begin earlier here and there would be less of the paleotopography preserved. Because areas outside of the low-relief surfaces are not used in the analysis, this simply means that we have less data to constrain the surface at these locations.
Effects of Non-Linear Relationships Between Channel Steepness and Erosion Rate
Previous models used to support the conflicting hypotheses, are based on the linear stream power model (i.e., n = 1) and if this underlying model is incorrect then the hypotheses might be slightly different (Goren et al., 2014). We do not attempt to quantify the implications of these underlying assumptions and have adopted the simplified linear stream power model. If n > 1, river segments that are steeper will propagate upstream faster than less steep river segments and vice versa. This would lead to the margins of the low-relief landscape having different
Cosmogenic-derived catchment averaged erosion rates from across SE Tibet show a non-linear relationship with catchment averaged normalized channel steepness and this suggests that n = 2 (Ouimet et al., 2009). This is based on the assumption that bedrock erodibility is spatially uniform across the study area, that the catchments analyzed are in steady state, and that erosion rates are set by fluvial erosion. Therefore, we explore how increasing the value of n from 1 to 2 would impact our analysis. Royden and Perron (2013), show that the steady state stream power model can be written as
where u* is a simplification of
FIGURE 9. Nonlinearities in the stream power model. (A) Normalized channel steepness values for n = 0.5. Overall magnitudes are reduced and there is less spatial variability in the u* map; (B) If n = 2, more spatial variability is predicted, and magnitudes increase. The box shows the regional extent of the map in Figure 10.
In an attempt to determine a value for n for a separate study, Ma et al. (2020) selected the catchment wide erosion rates from Ouimet et al. (2009) for only the Dadu catchment. They found a linear relationship between channel steepness and catchment averaged erosion rates could explain the data. It is unclear whether this is simply because the area analyzed was much smaller or whether slightly different methods were used to calculate normalized channel steepness. Therefore, the value of n for SE Tibet remains debated.
Effects of Spatial Variations in Erodibility
Importantly, we assume that over the scale of SE Tibet there is an effective average erodibility, enabling us to calculate
FIGURE 10. The influence of lithology on predicted normalized rock uplift rate. (A) Intrusive rocks colored by the predicted normalized rock uplift rates with the faults included as in Figure 5B. A band of high u* can be seen at approximately 28oN 97oW which corresponds to a dense band of intrusive units. However, high and low values of u* are observed across the area; (B) The opposite of Figure 10A showing the predicted normalized rock uplift rates for the other lithologies, with the outcrops of intrusive rocks shown in black. High and low values of u* are observed across both lithological categories and the correlation is not clear. Geological units are simplified after Pan et al. (2004).
It is also important to recognize that other factors may lead to differences between the erodibility of the trunk rivers and the low-relief surfaces, potentially enhancing the development of these low-relief surfaces by increasing K (decreasing our inferred u*), facilitating trunk river incision. For example, the bedload of the large rivers may contain a significant amount of material from hard lithologies from higher up the catchment that would not make up the bedload of the locally drained low-relief surfaces. This would lead to a difference in the contrast between bedrock and bedload in the trunk rivers compared to the low-relief surfaces, and it is this contrast that may lead to local differences in erodibility (Sklar and Dietrich, 2001). In addition, it is clear that the trunk streams flow parallel to major structures, such as suture zones and strike slip faults. It may be expected that these structures would locally increase the erodibility (Roy et al., 2016) and this may facilitate erosion of the trunk rivers and the preservation of the potentially harder low-relief surfaces. Both of these factors may also change through time as different rocks are exhumed leading to changes in bedload or different structures are revealed for rivers to incise into.
We have also assumed that precipitation is uniform in order to relate calculate
Our approach provides a means to quantitatively assess the continuity of a hypothetical paleo-topography with relatively uniform normalized rock uplift rate. We find that large variability in u* is required if surface uplift is uniform, but that a relatively smooth surface can be found if spatial variability is permitted. However, even with this increased flexibility, anomalously high u* values are found along the main trunk streams. In addition, differences between model-predicted elevations and observed elevations are often positive one side of a drainage divide, but negative on the other side suggesting that these tributaries are out of equilibrium with the regional topography. It should be noted that these residuals only correspond to the low-relief high-elevation parts of the landscape. In order for the u* to be smoother and more uniform, without significantly increasing model misfit, a rougher S.U. map is required, however, this would require short wavelength surface uplift structures, as opposed to a long wavelength mechanism such as lower crustal flow (Clark and Royden, 2000). Ultimately, this choice represents a compromise between roughness and data fit, however, our framework provides a method to ask: what is the relative importance of the two hypotheses? The preferred inversion result (Figure 4) fits almost 75% of the analyzed low-relief nodes to within 100 m, however a model that has smooth u* and S.U. maps, consistent with Hypothesis 1 (α = 102,λ = 101.5) only fits approximately 55% of the data. Therefore, low-relief paleotopography might have existed but the processes required by this model (reduced to its simplest form) only explains about half the data.
While our results lend support to the existence of a low-relief surface that has been uplifted and dissected, this model cannot explain all of the data, and additional factors may have modified erosion rates across the low-relief surfaces, either driven by climatic changes (Zhang et al., 2016; Nie et al., 2018), drainage reorganization (Clark et al., 2004; McPhillips et al., 2016) by drainage divide migration (Yang et al., 2015; Willett, 2017). In addition, the surface is likely to have been deformed by local faulting and lateral movement of low-relief topography as SE Tibet extrudes. However, then the role of a lower crustal low-viscosity channel becomes less clear. What is clear from this analysis is that a topographic surface interpolated between low-relief surfaces found across SE Tibet cannot be used to measure geodynamic processes in time and space. There is limited evidence for low-relief paleotopography in the current topographic data alone and other processes may have entirely shaped modern SE Tibet. Our approach also highlights the limit of utilizing topography alone to interpret landscapes in terms of hypotheses as it is can be an ill-posed inverse problem with multiple scenarios that will lead to indistinguishable features. It is this non-uniqueness that leads to debates about the development of SE Tibet (Yang et al., 2015; Whipple et al., 2017a; Whipple et al., 2017b; Willett, 2017). Future work could use our approach and incorporate diverse datasets and more complexity directly into the analysis. This could include spatially and temporally variable precipitation and erodibility, lateral translation of topography and river networks or hypothesized drainage capture events.
Although our approach provides a means to quantitatively diagnose the continuity of a hypothetical low-relief topographic surface that follows expected slope-area scaling, it does not provide a quantitative test of Hypothesis 1 or 2. In part, this is because the exact smoothness of the two surfaces and degree of data misfit required for Hypothesis 1 have not been specified in the literature. Our approach incorporates the strengths of the different approaches used to extract information from river networks, outlined in Using Fluvial Metrics to Test the Hypotheses, but adds additional flexibility. Furthermore, because we are not extracting individual channels from the DEM our approach potentially uses more of the data and we are not selecting which channels to highlight. In order to really understand the topographic development of SE Tibet, data that resolve temporal changes in erosion through time at the scale of individual patches of low-relief topography are required. Thermochronology would be the ideal tool to measure the exhumation rate history across the low-relief landscapes and the surrounding canyons (Gourbet et al., 2020) however, the amount of relief change is right at the limits of the method (Braun, 2002). What is clear from thermochronology is that exhumation rates have been highly variable in space and time, with the amount of exhumation at some locations being far greater than the amount predicted from valley incision alone (Zhang et al., 2016; Cook et al., 2018; Liu-Zeng et al., 2018; Ge et al., 2020; Replumaz et al., 2020). If this exhumation is associated with local high rates of rock uplift, and if the erosion rates on the surrounding low-relief surfaces remain low, large amounts of S.U. would be observed leading to obvious discontinuities in the low-relief surface. These have not been observed and thus this highlights the importance of processes that create low-relief surfaces in situ, supporting Hypothesis 2.
In order to falsify Hypothesis 2, future work could search for an absence of low-relief landscapes at intermediate elevations. This is challenging however as Hypothesis 2 predicts that landscapes that have lost upstream drainage area progressively lose relief and are thus most prominent at the highest elevations. Furthermore, the landscape is developing in an area of active deformation. This means that tectonics may respond to changes in surface topography in such a way that sets a limit on how high the low-relief landscapes can be uplifted. This would mean that all the low-relief landscapes might attain coherent, gradually varying, elevations set by the strength of the crust and faults. If only the elevations of the low-relief surfaces were used, this would support Hypothesis 1. We have not mapped low-relief surfaces across SE Tibet and have used published surfaces for our analysis (Clark et al., 2004), however, there does appear to be low-relief landscapes at intermediate elevations across the Upper Yangtze River (Liu et al., 2019). Alternatively, in situ erosion rate histories within the low-relief landscapes could be measured using low-temperature thermochronometry. The fingerprint of Hypothesis 2 would be decreases in exhumation rate associated with area loss coupled with increases in exhumation in surrounding canyons.
Data Availability Statement
Publicly available datasets were analyzed in this study. This data can be found here: https://www.hydrosheds.org.
All authors contributed to the design, analysis, and writing of the manuscript.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
We thank G. Hilley, J. P. Avouac and two reviewers for providing helpful feedback and S. D. Willett for stimulating discussions. This study was supported by NERC (NE/N015479/1).
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/feart.2020.587597/full#supplementary-material
Bernard, T., Sinclair, H. D., Gailleton, B., Mudd, S. M., and Ford, M. (2019). Lithological control on the post-orogenic topography and erosion history of the Pyrenees. Earth Planet Sci. Lett. 518, 53–66. doi:10.1016/j.epsl.2019.04.034
Botsyun, S., Sepulchre, P., Donnadieu, Y., Risi, C., Licht, A., and Rugenstein, J. K. C. (2019). Revised paleoaltimetry data show low Tibetan Plateau elevation during the Eocene. Science 363 (6430), eaaq1436. doi:10.1126/science.aaq1436
Clark, M. K., Royden, L. H., Whipple, K. X., Burch el, B. C., Zhang, X., and Tang, W. (2006). Use of a regional, relict landscape to measure vertical deformation of the eastern Tibetan Plateau. J. Geophys. Res. Earth Surface 111, 187-10. doi:10.1029/2005jf000294
Clark, M. K., Schoenbohm, L. M., Royden, L. H., Whipple, K. X., Burchfiel, B. C., Zhang, X., et al. (2004). Surface uplift, tectonics, and erosion of eastern Tibet from large-scale drainage patterns. Tectonics 23 (1). doi:10.1029/2002tc001402
Cook, K. L., Hovius, N., Wittmann, H., Heimsath, A. M., and Lee, Y.-H. (2018). Causes of rapid uplift and exceptional topography of Gongga Shan on the eastern margin of the Tibetan Plateau. Earth Planet Sci. Lett. 481, 328–337. doi:10.1016/j.epsl.2017.10.043
Croissant, T., and Braun, J., (2014). Constraining the stream power law: a novel approach combining a landscape evolution model and an inversion method. Earth Surf. Dynam. Discuss. 2 (1), 155–166. doi:10.5194/esurf-2-155-2014
Deng, B., Chew, D., Mark, C., Liu, S., Cogné, N., Jiang, L., et al. (2020). Late Cenozoic drainage reorganization of the paleo-Yangtze river constrained by multi-proxy provenance analysis of the Paleo-lake Xigeda. Geol. Soc. Am. Bull. doi:10.1130/b35579.1
Egholm, D. L., Jansen, J. D., Brædstrup, C. F., Pedersen, V. K., Andersen, J. L., Ugelvig, S. V., et al. (2017). Formation of plateau landscapes on glaciated continental margins. Nat. Geosci. 10 (8), 592. doi:10.1038/ngeo2980
Fox, M., Bodin, T., and Shuster, D. L. (2015b). Abrupt changes in the rate of Andean Plateau uplift from reversible jump Markov chain Monte Carlo inversion of river profiles. Geomorphology 238, 1–14. doi:10.1016/j.geomorph.2015.02.022
Fox, M., Goren, L., May, D. A., and Willett, S. D. (2014). Inversion of fluvial channels for paleorock uplift rates in Taiwan. J. Geophys. Res. Earth Surf. 119 (9), 1853–1875. doi:10.1002/2014jf003196
Ge, Y., Liu-Zeng, J., Zhang, J., Wang, W., Tian, Y., Fox, M., et al. (2020). Spatio-temporal variation in rock exhumation linked to large-scale shear zones in the southeastern Tibetan Plateau. Sci. China Earth Sci. 63, 512–532. doi:110.1007/s11430-019-9567-y
Goren, L., Fox, M., and Willett, S. D. (2014). Tectonics from fluvial topography using formal linear inversion: theory and applications to the Inyo Mountains, California. J. Geophys. Res. Earth Surf. 119, 1651–1681. doi:10.1002/2014jf003079
Gourbet, L., Yang, R., Fellin, M. G., Paquette, J.-L., Willett, S. D., Gong, J., et al. (2020). Evolution of the Yangtze river network, southeastern Tibet: insights from thermochronology and sedimentology. Lithosphere 12 (1), 3–18. doi:10.1130/l1104.1
Henck, A. C., Huntington, K. W., Stone, J. O., Montgomery, D. R., and Hallet, B. (2011). Spatial controls on erosion in the three rivers region, southeastern Tibet and southwestern China. Earth Planet Sci. Lett. 303 (1-2), 71–83. doi:10.1016/j.epsl.2010.12.038
Kirby, E., Whipple, K. X., Tang, W., and Chen, Z. (2003). Distribution of active rock uplift along the eastern margin of the Tibetan Plateau: inferences from bedrock channel longitudinal profiles. J. Geophys. Res. Solid Earth 108 (B4), 2217. doi:10.1029/2001jb000861
Li, Y., Liu, M., Li, Y., and Chen, L. (2019). Active crustal deformation in southeastern Tibetan Plateau: the kinematics and dynamics. Earth Planet Sci. Lett. 523, 115708. doi:10.1016/j.epsl.2019.07.010
Liu, F., Gao, H., Pan, B., Li, Z., and Su, H. (2019). Quantitative analysis of planation surfaces of the upper Yangtze River in the Sichuan-yunnan region, Southwest China. Front. Earth Sci. 13 (1), 55–74. doi:10.1007/s11707-018-0707-y
Liu-Zeng, J., Zhang, J., McPhillips, D., Reiners, P., Wang, W., Pik, R., et al. (2018). Multiple episodes of fast exhumation since Cretaceous in southeast Tibet, revealed by low-temperature thermochronology. Earth Planet Sci. Lett. 490, 62–76. doi:10.1016/j.epsl.2018.03.011
McPhillips, D., Hoke, G. D., Liu-Zeng, J., Bierman, P. R., Rood, D. H., and Niedermann, S. (2016). Dating the incision of the Yangtze River gorge at the First Bend using three-nuclide burial ages. Geophys. Res. Lett. 43 (1), 101–110. doi:10.1002/2015gl066780
Molnar, P., Boos, W. R., and Battisti, D. S. (2010). Orographic controls on climate and paleoclimate of Asia: thermal and mechanical roles for the Tibetan Plateau. Annu. Rev. Earth Planet Sci. 38, 77–102. doi:10.1146/annurev-earth-040809-152456
Mudd, S. M., Attal, M., Milodowski, D. T., Grieve, S. W. D., and Valters, D. A. (2014). A statistical framework to quantify spatial variation in channel gradients using the integral method of channel profile analysis. J. Geophys. Res. Earth Surf. 119, 138–152. doi:10.1002/2013jf002981
Nie, J., Ruetenik, G., Gallagher, K., Hoke, G., Garzione, C. N., Wang, W., et al. (2018). Rapid incision of the Mekong River in the middle Miocene linked to monsoonal precipitation. Nat. Geosci. 11, 944–948. doi:10.1038/s41561-018-0244-z
Ouimet, W. B., Whipple, K. X., and Granger, D. E. (2009). Beyond threshold hillslopes: channel adjustment to base-level fall in tectonically active mountain ranges. Geology 37 (7), 579–582. doi:10.1130/g30013a.1
Pan, G. T., Ding, J., Yao, D. S., and Wang, L. Q. (2004). Guidebook of 1: 1,500,000 geologic map of the Qinghai–Xizang (Tibet) plateau and adjacent areas. Chengdu, China: Chengdu Cartographic Publishing House, 48.
Replumaz, A., San José, M., Margirier, A., van der Beek, P., Gautheron, C., Leloup, P. H., et al. (2020). Tectonic control on rapid late Miocene—quaternary incision of the Mekong River Knickzone, southeast Tibetan Plateau. Tectonics 39 (2), e2019TC005782. doi:10.1029/2019tc005782
Roering, J. J., Marshall, J., Booth, A. M., Mort, M., and Jin, Q. (2010). Evidence for biotic controls on topography and soil production. Earth Planet Sci. Lett. 298 (1-2), 183–190. doi:10.1016/j.epsl.2010.07.040
Roberts, G. G., Paul, J. D., White, N., and Winterbourne, J. (2012). Temporal and spatial evolution of dynamic support from river profiles: a framework for Madagascar. Geochem. Geophys. Geosyst. 13 (4), 1087. doi:10.1029/2012gc004040
Roy, S. G., Tucker, G. E., Koons, P. O., Smith, S. M., and Upton, P. (2016). A fault runs through it: modeling the influence of rock strength and grain-size distribution in a fault-damaged landscape. J. Geophys. Res. Earth Surf. 121 (10), 1911–1930. doi:10.1002/2015jf003662
Rudge, J. F., Roberts, G. G., White, N. J., and Richardson, C. N. (2015). Uplift histories of Africa and Australia from linear inverse modeling of drainage inventories. J. Geophys. Res. Earth Surf. 120 (5), 894–914. doi:10.1002/2014jf003297
Schaller, M., Ehlers, T. A., Lang, K. A. H., Schmid, M., and Fuentes-Espoz, J. P. (2018). Addressing the contribution of climate and vegetation cover on hillslope denudation, Chilean Coastal Cordillera (26°-38°S). Earth Planet Sci. Lett. 489, 111–122. doi:10.1016/j.epsl.2018.02.026
Scherler, D., DiBiase, R. A., Fisher, G. B., and Avouac, J.-P. (2017). Testing monsoonal controls on bedrock river incision in the Himalaya and Eastern Tibet with a stochastic-threshold stream power model. J. Geophys. Res. Earth Surf. 122 (7), 1389–1429. doi:10.1002/2016jf004011
Schoenbohm, L. M., Whipple, K. X., Burchfiel, B. C., and Chen, L. (2004). Geomorphic constraints on surface uplift, exhumation, and plateau growth in the Red River region, Yunnan Province, China. Geol. Soc. Am. Bull. 116, 895–909. doi:10.1130/b25364.1
Styron, R. H., Garcia, J., and Pagani, M. (2017). “The GEM Global Active Faults Database: the growth and synthesis of a worldwide database of active structures for PSHA, research, and education,” in AGU fall meeting abstracts, Washington, DC, December 11–15, 2017.
Suhail, H. A., Yang, R., Chen, H., and Rao, G. (2020). The impact of river capture on the landscape development of the Dadu River drainage basin, eastern Tibetan Plateau. J. Asian Earth Sci. 198, 104377. doi:10.1016/j.jseaes.2020.104377
Tan, X., Liu, Y., Lee, Y.-H., Lu, R., Xu, X., Suppe, J., et al. (2019). Parallelism between the maximum exhumation belt and the Moho ramp along the eastern Tibetan Plateau margin: coincidence or consequence? Earth Planet Sci. Lett. 507, 73–84. doi:10.1016/j.epsl.2018.12.001
Tapponnier, P., Peltzer, G. L. D. A. Y., Le Dain, A. Y., Armijo, R., and Cobbold, P. (1982). Propagating extrusion tectonics in Asia: new insights from simple experiments with plasticine. Geology 10 (12), 611–616. doi:10.1130/0091-7613(1982)10<611:petian>2.0.co;2
Whipple, K. X., DiBiase, R. A., Ouimet, W. B., and Forte, A. M. (2017a). Preservation or piracy: diagnosing low-relief, high-elevation surface formation mechanisms. Geology 45, 91–94. doi:10.1130/g38490.1
Whipple, K. X., DiBiase, R. A., Ouimet, W. B., and Forte, A. M. (2017b). Preservation or piracy: diagnosing low-relief, high-elevation surface formation mechanisms: reply. Geology 45, e422. doi:10.1130/g39252y.1
Whipple, K. X., Forte, A. M., DiBiase, R. A., Gasparini, N. M., and Ouimet, W. B. (2017c). Timescales of landscape response to divide migration and drainage capture: implications for the role of divide mobility in landscape evolution. J. Geophys. Res. Earth Surf. 122 (1), 248–273. doi:10.1002/2016jf003973
Yang, R., Suhail, H. A., Gourbet, L., Willett, S. D., Fellin, M. G., Lin, X., et al. (2020). Early Pleistocene drainage pattern changes in Eastern Tibet: constraints from provenance analysis, thermochronometry, and numerical modeling. Earth Planet Sci. Lett. 531, 115955. doi:10.1016/j.epsl.2019.115955
Zhang, H., Oskin, M. E., Liu-Zeng, J., Zhang, P., Reiners, P. W., and Xiao, P. (2016). Pulsed exhumation of interior eastern Tibet: implications for relief generation mechanisms and the origin of high-elevation planation surfaces. Earth Planet Sci. Lett. 449, 176–185. doi:10.1016/j.epsl.2016.05.048
Keywords: paleotopographic reconstruction, South East Tibet, geodynamics, landscape evolution, river incision
Citation: Fox M, Carter A and Dai J-G (2020) How Continuous Are the “Relict” Landscapes of Southeastern Tibet?. Front. Earth Sci. 8:587597. doi: 10.3389/feart.2020.587597
Received: 26 July 2020; Accepted: 14 October 2020;
Published: 12 November 2020.
Edited by:Paolo Ballato, Roma Tre University, Italy
Reviewed by:Adam Forte, Louisiana State University, United States
Wolfgang Schwanghart, University of Potsdam, Germany
Copyright © 2020 Fox, Carter and Dai. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Matthew Fox, email@example.com