Two-staged uplift of the southeast margin of the Tibetan plateau revealed by river longitudinal profiles

The Jinsha River, upstream of the Yangtze River, flows through the southeastern margin of the plateau, providing an opportunity to study the plateau uplift relative to the Sichuan basin through fluvial geomorphology. Eight tributaries of the Jinsha River located at the transition boundary from the plateau to the basin were analyzed for the feature of longitudinal profile and its evolution. The spatial distribution of slope-break knickpoints in the catchment was determined by slope-area analysis and integration analysis. Further, the magnitude and timing of the base level lowering (i.e., uplift of the Tibetan Plateau) were estimated by reconstructing the river longitudinal profile. The results show that there are two levels of slope-break knickpoints, with average elevations of 2865 m and 1523 m, corresponding to base levels relatively lowering of 1400 m and 1100 m, respectively. With reference to the long-term downcutting rate of the Jinsha River, the timings of these two phases base level lowering and uplifts are estimated to be ∼21–24 and ∼7–13 Ma, respectively. These are close to the timing of two uplifts revealed by inversion modeling of river longitudinal profiles. Therefore, the uplift of the studied region on the Tibetan Plateau can be divided into two stages, the early stage (∼20 Ma) uplift from 1500 to 2900 m and the late stage (∼10 Ma) uplift from 2900 to 4000 m, which further supports the staged outside uplift of the Tibetan Plateau.


Introduction
In the Cenozoic, a collision of the Eurasian and Indian plates formed the Tibetan Plateau (Figure 1), known as the third pole of the Earth, whose uplift has shaped the current geomorphic pattern of Asia, influencing the material cycle among the Earth's spheres, and having a profound impact on regional and global climate and environmental changes (Berner et al., 1983;An et al., 2001;Chen et al., 2021). Reconstructing the uplift process of the Tibetan Plateau facilitates elucidating the uplift and deformation mechanisms of the plateau and its environmental effects. Previous studies have shown that the uplift of the plateau is a multi-staged process, rather than a simple event (Ding et al., 2017;Cao et al., 2022). The timing and magnitude of the plateau uplift are still highly controversial, which is generally associated with different methods and study areas. For example, studies based on ancient OPEN ACCESS EDITED BY Daniele Spatola, Faculty of Mathematics, Physics, and Natural Sciences, Sapienza University of Rome, Italy pollens, palaeontological fossils, and a series of paleolandform from the plateau margins suggest a cumulative uplift of 3500-4000 m since the late Pliocene (Xu et al., 1973;Li et al., 1979;Deng et al., 2012;Chen et al., 2020), whereas other studies based on tectonic and thermal history suggest that the plateau initiated to uplift significantly earlier in the central plateau during the Eocene (Royden et al., 2008;Ding et al., 2017), and margins mainly started in Miocene (Wang et al., 2012a;Li et al., 2019). In addition, the debate on the uplift mechanism of the Tibetan Plateau has also mainly focused on two models: tectonic extrusion (Tapponnier et al., 1982;Tapponnier et al., 2001) and the lower crustal flow (Clark and Royden, 2000).
As a zone of steep topographic change in the transition from the plateau to the surrounding basin, the southeastern margin of the Tibetan Plateau was suggested to be the key for the study of plateau uplift. A range of evidence from tectonic deformation, topographic relief, sedimentary records, and isotope palaeoaltimetries indicate that the eastern edge of the plateau reached to high elevation during the Eocene (~45-40 Ma) (Wang et al., 2012a;Cao et al., 2019;Su et al., 2019;Xiong et al., 2020), while deformation at the eastern margin of the plateau shows another more widely recorded period (~15-10 Ma) of its uplift (Clark et al., 2005;Ouimet et al., 2010;Tao et al., 2020). On the southeastern edge of the Tibetan Plateau, the palaeo-elevation reconstructed by stable isotopes shows that the plateau existed in the area during the Eocene-Oligocene, while the southernmost part of the region was uplifted to its present elevation during the Miocene (Cook et al., 2013;Hoke et al., 2014;Li et al., 2015;Tang et al., 2021). Furthermore, some researchers consider the time when large rivers incised deeply into the plateau as the timing of plateau uplift. For instance, thermochronology studies on the Dadu River and Yalong River constrain the uplift of the southeastern plateau at~9-13 Ma by river incision histories (Clark et al., 2005;Ahmed, 2019).
Therefore, the results obtained by different methods in this area are inconsistent, probably due to the limitations and uncertainties of these methods. For example, low-temperature thermochronological methods have been widely applied for constraining the timing of the plateau uplift by the denudation rate, but cannot estimate the detailed magnitude of uplift. Stable isotope palaeoaltimetries can calculate the magnitude of uplift, but they are always influenced by the "inverse altitude effect" due to temperature and other effects (Jing et al., 2022). The deformation analysis can also be probably affected by activities in the local domains. Hence, the uplift process on the southeastern edge of the plateau requires to be further investigated through diversified evidence.
In recent years, fluvial geomorphology has been successfully applied to study tectonic uplift, with channel steepness becoming a good proxy for differential uplift at the active edge of the plateau (Kirby et al., 2003). In the southeastern Tibetan Plateau, the Yangtze River and its tributaries, one of the largest drainage systems in the world, originate and flow through the region.
It is an ideal place for studying structural geomorphology concerning the uplift and expansion of the Tibetan Plateau since the Cenozoic, especially the Miocene. Therefore, we target the upper reaches of the Yangtze River (Jinsha River and its tributaries) within the plateau-basin transition zone of the southeastern Tibetan Plateau ( Figure 1) and utilize the fluvial geomorphology analysis methods (e.g., mobile knickpoint analysis, longitudinal river profile reconstruction and inversion simulation of the river bedrock channel) to constrain the uplift history of the southeastern margin of the Tibetan Plateau.
2 Background of the study area

Geological background
The study area, located in the southeastern Tibetan Plateau, falls within the Kangding-Xichang sub-region of the Qiangtang-Sanjiang tectonic zone, which is one of the major tectonic units of the Tibetan Plateau, tectonically belonging to the north-south extending Sichuan-Yunnan tectonic belt. The northern section of the left strike-slip Xiaojiang Fault ( Figure 2) cuts through the region as the regional main controlling fault, which was strongly active since Miocene (Chen et al., 2011). The Palaeozoic-Cenozoic sedimentary/ metasedimentary rocks, dominated by sandstones, shales, limestones, marls, dolomites, and bioclastic rocks, are exposed in the majority of the study area, as well as several Precambrian basement outcrops. The large Permian basalt (E'meishan basalt) is widely distributed in the area, with a few intrusive rocks only scattered in the western part of the study area. Quaternary sediments are mainly distributed in faulted basins and river valleys (Figure 2).
To the east, the Sichuan Basin, a superimposed basin, has undergone multiple phases of tectonic evolution and sedimentary infilling. The southwestern part of the Sichuan Basin was mainly characterized as an intra-land depressional basin during the Cenozoic (Li et al., 1995;Chen et al., 2008;He et al., 2011). The basin during Paleogene mainly received river-lake deposits, with a

FIGURE 2
Simplified geological map of the study area.

FIGURE 3
Topography of the study area and the range of each catchment.

Frontiers in Earth Science
frontiersin.org 03 gradually decreasing sedimentation. The enhanced tectonic uplift at the end of the Paleogene caused widespread interruption of deposition, and deposition along the piedmont belt began to develop in the Pliocene (Huang and He, 2016).
Large sutures (Wang et al., 2012b) or regional fault zones (Li et al., 2018) along the southwest boundary of the Sichuan basin have been suggested. The north segment was identified as a thrust fault (Wang et al., 2012b). Hence, it is reasonable to deduce the uplift of studied region is caused by the movement of this thrust fault (basin boundary fault).

Morphological setting
The study area, located at the boundary between Sichuan Province and Yunnan Province of China, is situated at the intersection of three geographic units, i.e., the Western Sichuan Plateau, the Yunnan-Guizhou Plateau, and the Sichuan Basin, and exhibits a sharp decline in topography from southwest to northeast ( Figure 3). The geographical coordinates are 26.5°-29.5°N and 102.5°-105.5°E, with a total area of 144,700 km 2 . Within this territory, the highest elevation is 5877 m and the lowest is 219 m, with an elevation difference of 5658 m. The swath profile shows there is a relic plateau surface situated at~2500 m, which is separated by river valleys incised into the plateau. The Jinsha River has the maximum incision depth, up to~1500 m. The transition zone between the plateau and the Sichuan Basin is approximately 100 km wide, with a decreasing altitude from 2500 to~600 m ( Figure 3). This region is dominated by terrain with elevations of 1000-2500 m. The proportion of moderately undulating topography with a relief of 500-1000 m is the highest, accounting for 42.26% of the total area. Daily precipitation data in the region show that rainfall is concentrated from May to October, with a mean annual rainfall of 812 mm . The mean annual temperature is 16°C-18°C and the water-heat combination is generally consistent across the region, belonging to the same subtropical monsoon climate.
In the region, the Jinsha River, the upper reaches of the Yangtze River, spans Sichuan, Tibetan, and Yunnan provinces, with the largest tributary, the Yalong River, and the Min River. The Min River joins the Jinsha River in Yibin, Sichuan Province ( Figure 1 and Figure 2), and the main river starts to be called the Yangtze River. In the study area, there are eight tributaries of the Jinsha River and their catchments distributed between Panzhihua and Yibin, namely, Cha River, Heishui River, Xixi River, Yiutong River, Xisujiao River, Mabian River, Mashu River, and Heng River ( Figure 3). These eight rivers all originate from the plateau surface and share the Sichuan basin as the base level. The Mabian River and Heng River flow into the Sichuan basin and then join the Jinsha River, while the other six rivers join the Jinsha River directly within the plateau. The detailed information on the rivers (Table 1) shows that these eight rivers all have a high relief, up to >2000 m. The length of the rivers varies significantly, with the longest Heng River of~290 km. The Mashu River is the shortest (32 km) one, and the others are between 100 and 200 km long. Further, the Xisujiao River has the smallest basin area of 707 km 2 and the Heng River gives the largest basin of 15016 km 2 . According to the classification scheme of the geomorphic stages of evolution (Strahler, 1952), the eight catchments are mainly in the mature stage. The longitudinal profiles cannot generally be fitted with an exponential function, indicating that the longitudinal profiles have not yet reached the equilibrium state. The hypsometric integral (HI) of eight catchments is found to be 0.31-0.60 with an average of 0.51 based on elevationrelief ratio method, which indicates the most watersheds are at mature stage.
In the upper panel, the purple line represents the southwestnortheast topographic line in the study area, with a 5 km buffer zone to obtain the swath profile shown in orange-pink. The white lines represent the boundaries of the river basin.➀ Cha River; ➁ Heishui River; ➂ Xixi River; ➃ Liutong River; ➄ Xisujiao River; ➅Mabian River; ➆ Mashu River; ➇ Heng River. The lower panel shows the swath profile.

Data and methods
The DEM data were downloaded using the software of 91 Map Assistant. This software implements the download of reprocessed data of ALOS PALSAR 12.5 m resolution DEM, provided by NASA. The initial resolution is 18.6 m, so we resampled it to 20 m resolution to facilitate calculation and simulation.

Migrating knickpoints
Using the stream power incision model can obtain the Slope-area plots where the slope of the fitted line is the concavity of the channel and the intercept of the line is the steepness of the channel (Kirby et al., 2003). Two types of knickpoints can be identified on the Slope-area plots, i.e., vertical-step knickpoint and slope-break knickpoint (Kirby and Whipple, 2012). The channel slope increases sharply near the knickpoint, and the steepness index above and below the point remains unchanged, which is called a vertical-step knickpoint, also known as an immobile knickpoint, and its formation is related to the high erosion resistance of the local bedrock. The slope-breaking knickpoint is also known as a mobile knickpoint, and its formation is related to the decline of the base level. When the base level of a river drops, a slope-breaking knickpoint will form, and migrate upstream as the headward incision of the river, and its vertical migration distance is only time-dependent, resulting in the knickpoint being located at the same elevation in all tributaries of the catchment (Niemann et al., 2001).
Since information may be lost when the original topographic data are smoothed and resampled during the production of the Slope-area plot, Perron and Royden (2013) proposed the integration analysis method, where the catchment area is integrated over the distance from the river source, converting the river longitudinal profile into a straight line with the slope indicating the channel steepness. The uncertainties of the resulting steepness indices are extremely small, due to avoid the errors introduced repeatedly when differentiating the elevation. Slope-area analysis can only determine the distribution of knickpoints for one channel at a time, whereas the integration analysis can determine the distribution of knickpoints for all channels in a catchment simultaneously. Here, we combined both methods for comparison to determine the knickpoints.
Knickpoints form and migrate upstream when the base level drop. The location of the knickpoint records the timing of the relative decline of erosional basis, while the downstream steepness index reflects the lowering rate. Thus, the complete uplift history of the region can be inferred from the migrating knickpoints on the longitudinal profile of the river. At present, the inversion of uplift history is based on the linear hypothesis Rudge et al., 2015;Pedersen et al., 2018).

Analysis of knickpoints
Slope-area analysis was implemented using the ArcGIS stream profiler toolbox and MATLAB programs, with a smoothing window of 250 m (Whipple et al., 2007). Differently, Chi analysis was performed in Matlab using the Topotoolbox (Schwanghart and Scherier, 2014). By comparing the results calculated with different parameters, a minimum area of the catchment of 10 5 m 2 , and a smoothing window of 100 m in the integration method were adopted.

Reconstruction of longitudinal profile
There is a proportional relationship between the channel slope and the basin area above the knickpoint, and the principle of the channel profile reconstruction is to project this proportional relationship below the knickpoint to obtain the palaeochannel (Kirby and Whipple, 2012). In general, directly using the steepness index (k s ) and concavity index (θ) of the river above the knickpoints will produce a large error in the projected paleochannel, so we adopted the method proposed by Forte and Whipple (2019), in which, the original chi and elevation values were preprocessed by unsupervised learning. The difference in elevation between the paleochannel and the modern channel is the downcutting depth, which is the lowering amplitude of the base level near the start point of headward erosion in this study.

Modelling of river uplift history
Under conditions of homogeneous watershed lithology, uniform rainfall distribution (K is a constant), and a linear relationship between erosion and channel slope (n = 1), from the start point, the elapsed time of a knickpoint migration can be expressed as: (1) The temporal information of knickpoint migration is transformed to be represented by the river Chi(χ), and the solution obtained using the characteristic linear method, given that the rate of watershed uplift is spatially uniform, is where the unit of t* is the same as Chi and U* is the dimensionless uplift rate . The present elevation (t=0) of the river can be expressed as the integral of the uplift rate over the response time: Differentiation of the Chi values gives the river uplift history as: where U* (t*) is the dimensionless uplift history, which can be transformed into the dimensioned uplift history U(t) by Eq. 3. The erosion coefficient K is a key parameter in the inversion of the uplift history of a river, and the lower the erosion rate, the longer the uplift history obtained. Considering the size of the studied watershed, we adopt the measured erosion coefficient (K = 1.11 ± 0.10 10 −6 m 0.1/ a) of the neighboring Dadu River for the calculation (Ma et al., 2020), which has a similar geological, geomorphological and climatic conditions to studied rivers. Further, Goren et al. (2014) suggested that the best inversion results exist at intervals with time nodes of 10, 20, and 40. It was calculated that the uplift signals were similar at different time intervals and a time node of 40 is chosen in this study. We use a set of open-source codes, RiverProAnalysis, which was based on the Matlab platform and integrated all of these functions .

Frontiers in Earth Science
frontiersin.org 05  Frontiers in Earth Science frontiersin.org 06 4 Results

Quantitative analysis of river headward erosion
Slope-area analysis has been conducted to obtain the slope-break knickpoints for each river and the k sn values for the different river reaches ( Table 2). Most of the rivers can be divided into three reaches for linear fitting. The k sn values tend to increase downstream, indicating that there are slope-break knickpoints in the rivers (Figure 4). For example, two knickpoints in the Heng River locate at 1767 and 2801 m; only one knickpoint in the Mashu River, Mabian River, and Xisujiao River, which are located at 2742, 1543, and 2134 m, respectively ( Table 2).
The slope-break knickpoints of these eight rivers are concentrated at three zones with similar elevations, with six at high elevation (2742-3026 m), two at middle elevation (2134-2369 m), and five at low elevation (1342-1767 m). The average elevation of each zone is 2865, 2251, and 1523 m, respectively (Table 2).
Integration analysis is able to identify the distribution of knickpoints on the mainstem and tributaries of each river. The obvious turning points on the χ-z plot indicate that the slope-break

Frontiers in Earth Science
frontiersin.org knickpoints, which are generated by the base level lower when they concentrate at the same elevation. Using a threshold of 300 m to find slope break knickpoints, it is possible to determine the location of main mobile knickpoints ( Figure 5). The frequency histogram plot shows that the knickpoints of each river are not concentrated, but that of all rivers exhibits a unimodal distribution with a peak at 1750 m, which is close to the clustering at 1525 m given by the slope-area analysis.

Amplitude and timing of base level lowering
We choose the Mabian River and Heng River, which take the Sichuan Basin as the base level, to reconstruct the longitudinal profile, and use the bottom edge of the Sichuan Basin as the reference point to calculate the elevation difference between the reconstructed and modern channels as the amplitude of the base level lowering. The results show that the elevation difference below the knickpoint increases rapidly with distance, reaching a maximum at the reference point. Below the reference point, the elevation difference of the Mabian River decreases slightly with distance, while the elevation difference of the Heng River remains almost constant ( Figure 6).
The reconstructed longitudinal profile based on the upstream above a slope-break knickpoint (1543 m) in the Mabian River indicates an 800 m base level lowering. There are two slope-break knickpoints in the Heng River (2801 and 1767 m), and the results indicate corresponding 2500 ± 150 and 1400 m base level lowering, respectively (Table 3). However, the former datum may be subject to some error due to the short river length above the 2801 m knickpoint.
Alluvial deposits are present in caves at different elevations on the valley slope near Shigu of the Jinsha River, and burial ages indicate that the river incised~1 km between 18 and 9 Ma, corresponding to a downcutting rate of~0.11 mm/a (McPhillips et al., 2016). In addition, near Weixi on the right bank of the Jinsha River, Liu-Zeng et al. (2018) collected thermochronology samples along the valley slope and obtained a denudation rate of~0.08 mm/a since~18 Ma in the area by river incision. These two rates are generally consistent. Hence, assuming a long-term river incision rate of~0.11 mm/a in the lower Jinsha River allows calculating the age of base level lowering based on the estimated amplitude of base level lowering. The base level lowering of 2500 ± 150 m yields an approximate age of 21-24 Ma, while the lowering of 800-1400 m gives an approximate age of 7-13 Ma.

River uplift history inversion
Our inversion modeling result (Figure 7) shows there are differences in uplift histories between the Mabian River and Heng River. The uplift rate of the Mabian River decreases with

FIGURE 6
Reconstructed longitudinal profiles of the Mabian River and Heng River (A,B) and estimated river downcutting and base level lowering (C,D). The black lines represent the river longitudinal profiles, the red lines represent the reconstructed longitudinal profiles, dotted lines represent the elevation difference between the actual and predicted river longitudinal profiles. The blue circles indicate the location of the edge of the basin bottom. The blue arrow points to the line is calculated using projection from the low knickpoint, while the red arrow points to, that is, calculated using projection from the high knickpoint. The uncertainty is given in a 95% confidence interval.
Frontiers in Earth Science frontiersin.org fluctuation over time, with the larger uplift rate (~0.34 mm/a) occurring at~17-20 Ma, then dropping to~0.18 mm/a at 16.3 Ma, with insignificant fluctuations thereafter, and then declines again to~0.07 mm/a at~6.8 Ma. Finally, this low uplift rate has been maintained to the present ( Figure 7A). The uplift rate of the Heng River fluctuates over the whole period, with the first peak (~0.28 mm/a) occurring in the Early Miocene (~19-22.7 Ma), reaching a second peak (~0.42 mm/a) in the Late Miocene (~5.8-7 Ma) ( Figure 7B), indicating that two significant uplifts occurred in the region. Both rivers display a quick uplift at 20 Ma. However, the Mabian River does not exhibit the quick uplift at~5.8-7 Ma, which is seen for the Heng River.

Geological influence
The studied eight rivers flow across the southeastern margin of the plateau where Paleozoic and Mesozoic sedimentary rocks are widely exposed in a masaic pattern (Figure 2), which implies this region has a relative homogenous bedrock on a large scale. Generally, the slopebreak knickpoints in the first-order tributaries of the lower Jinsha River obtained from the stream power incision model are mainly distributed in sediments (Figure 2), which may be conducive to the preservation of mobile knickpoints.
However, the erodibilities of sandstones, shales, limestones, marls, dolomites, and bioclastic rocks are usually different, hence the lithologic difference may produce lithologically controlled knickpoints. Except the strike-slip faults, no normal or reverse fault in a large scale has been reported . Even so, the undefined faults or undiscovered faults may produce tectonically controlled knickpoints.
Lithologically controlled and tectonically controlled knickpoints will blend with those knickpoints originated from base level lowering. The observed knickpoints that do not cluster at an altitude may be influenced by local factors. Due to only two slope-break knickpoints at the middle elevation are obtained by the slope-area analysis, they are probably not related to the base level lowering. Projecting knickpoints on the geological map, it can be found that most of the slope-break knickpoints do not have a clear correspondence with stratigraphic boundaries or faults (Figure 2), further suggesting that the formation of these knickpoints is related to the base level lowering.
The fitting of upstream data points of the Xisujiao River and Mabian River is good, while the fitting of downstream data points which exhibit a large scatter is poor (Figure 4). The main reason why the downstream data points do not show an ideal linear distribution is probably that these two rivers have a smaller catchment area and a weaker stream power than other rivers. After the potential slopebreak knickpoints migrate into these tributaries, they are more affected by the river substrate.

Two-stage plateau uplift revealed by changes in base level
The base level reconstructed on mobile knickpoints indicates that the lowering of the Sichuan Basin relative to the southeastern margin of the Tibetan Plateau can be divided into two stages, with the late stage (7-13 Ma) showing a relative lowering of~1100 m and the early stage (21-24 Ma) showing a relative lowering of~1400 m. The changes in the base level may be attributed to 1) the uplift of the Tibetan Plateau, 2) the graben subsidence of the Sichuan Basin, and 3) the intense denudation of the Sichuan Basin. It has been shown that the Sichuan Basin is mainly surrounded by reverse faults, rather than a graben basin (Xiong et al., 2015). Therefore, the base level lowering cannot be attributed to basin subsidence. Quantitative thermal history simulations of fission tracks in apatite samples indicate that the Sichuan Basin has suffered denudation with a thickness generally greater than 1 km during the Neogene (Deng

FIGURE 7
Uplift history simulated by the inversion method for the Mabian River (A) and the Hengjiang River (B). The grey bars indicate the period with an increasing uplift rate.
Frontiers in Earth Science frontiersin.org et al., 2009). This indicates that Miocene denudation did occur, with a denudation rate of approximately 0.056 mm/a, which is significantly lower than the long-term downcutting rate of the river valley (0.11 mm/a). Further, based on a simple isostatic calculation, denudation only results in an elevation loss of only about 1/5 of the denudation thickness. It is difficult to form a large knickpoint due to pure denudation. Considering that the denudation also occurs in the plateau, we infer that the basin denudation is not the main cause of the base level lowering. There is also the possibility that the Sichuan Basin did not subside relative to the plateau, but that a mobile knickpoint crossed the Sichuan Basin and entered the plateau. For example, it is shown that the tributary (Yuan River) flowing into the Jianghan Basin has maintained a low erosion rate (<0.04 km/Ma) since~30 Ma, whereas the tributary (Wu River) flowing into the Sichuan Basin has experienced a sudden increase in downcutting rate to 0.07 km/ Ma at 21-16 Ma. This difference is thought to reflect a decrease in the base level of erosion due to the Yangtze River cutting through the Three Gorges (Figure 1; Jiao et al., 2022). However, the altitude difference between the Sichuan Basin and Jianghan Basin is approximately 300 m, meaning that the base level change due to the Three Gorges cutting-through is approximately equal to the altitude difference, which is obviously smaller than the approximate 1 km base level lowering estimated for the Jinsha River. This suggests that the observed large lowering of the Sichuan Basin relative to the plateau is not related to cutting through the Three Gorges.
Further, there is another possibility that the sudden increase in the incision rate of the Wu River at 21-16 Ma indicates the uplift of the Yungui Plateau, further southeast to the Tibetan Plateau, where the catchment of the Wu River is located (Figure 1). Because the timing of increase incision rate of the Wu River is close to the early change stage (21-24 Ma) of the Mabian and Heng rivers in this study. The Tibetan Plateau gradually declines in height towards the southeast, and gradually transits to the Yungui Plateau. When a tectonic uplift occurs on the Tibetan Plateau, it is logical that the YunGui Plateau also rises correspondingly.
Therefore, the two stages of relative lowering of the base level (Sichuan Basin), as derived from the mobile knickpoints, should be a result of plateau uplift. At least two tectonic uplifts should have occurred in the studied region on the southeastern margin of the Tibetan Plateau, with~1400 m uplift in the early stage (~21-24 Ma) and~1100 m uplift in the late stage (7)(8)(9)(10)(11)(12). This is also generally supported by the results of inversion modeling of river longitudinal profiles. The change in channel uplift rate by inversion modeling occurred at 20 and~5.8-7 Ma, respectively, which is comparable with those of the base level lowering method. The change in base level is slightly earlier than the change in channel uplift rate, which may be related to the uncertainty in the erosion coefficient (1.11 × 10 −6 m 0.1/a) or the long-term incision rate of the valley (0.11 mm/a). Generally, the similarity of the two results indicates that they are controlled by the same geological process, that is, the process of regional uplift.
Furthermore, although the modeled uplift histories of Heng River and Mabian River differ somewhat, both show high uplift rates at~20 Ma, which is consistent with the early stage of uplift (21-24 Ma) of the plateau (Figure 8). The different uplift histories of the two rivers in the late stage may be due to the influence of local fault activity, requiring further work.

Comparison with other uplift records
Our study on the river longitudinal profile shows that the studied region on southeastern Tibetan Plateau underwent two rapid and large uplifts phases relative to the Sichuan Basin during the Miocene, which are generally consistent with previous studies on the eastern margin of the Tibetan Plateau. For example, in the Longmen Mountains, a transition region between the Sichuan Basin and the Tibetan Plateau, Kirby et al. (2003) first used the biotite 40 Ar/ 39 Ar of and zircon/apatite (U-Th)/He to confine the rapid denudation of the central Longmen Mountain at~5-12 Ma. Subsequently, Godard et al. (2009) used thermochronological forward modeling of vertical section data to update the age to 8-11 Ma. In the southern of Longmen Mountain, thermochronological studies suggest that the timing of regional rapid denudation events was confined to 8-10 Ma (Guenthner et al., 2014;Tian et al., 2014). Based on thermal data from the western margin of the Sichuan basin and from along the Xianshuihe fault, the counterclockwise bending/rotation of the Frontiers in Earth Science frontiersin.org Sichuan basin was initiated in late Cenozoic time (~13 Ma) . In addition, Wang et al. (2012a) identified two events of rapid denudation in the central Longmen Mountain, one occurring at the Late Oligocene to Early Miocene (20)(21)(22)(23)(24)(25)(26)(27)(28)(29)(30) and the other occurring at 10-15 Ma; both events in this area were mainly controlled by the activity of Longmen Mountain Fault. Similarly, Zhang et al. (2016) reported two rapid exhumation events in the adjacent Jiulong area during the Early Oligocene (~30 Ma) and Late Miocene (~7-8 Ma), indicating a correlation with the regional activity of the Jiulong Fault. Thermochronological data from the Deqin area of northern Yunnan show two rapid cooling during the Early Cenozoic and Miocene (<~20 Ma), which was mainly influenced by fault activity and river incision in the area (Liu-Zeng et al., 2018). A recent study suggests that rises in the east and west parts of the Northern TP during 15 to 7 Ma reached a high elevation close to that of today, which can be correlated to the second-stage uplift in the lower reach of Jinsha River (Miao et al., 2022). Additionally, Clark et al. (2005) suggest that the timing of major rivers incising into the plateau at the southeastern margin can be considered as the onset of plateau uplift. Rapid erosion of the Dadu River occurred at~10 Ma, with an erosion rate of 0.33 ± 0.04 km/Ma. Results from the Yalong River's transect indicate that rapid erosion occurred earlier than 14 Ma, with an erosion rate of 0.34 ± 0.02 km/Ma. All these correspond with the Rapid erosion of the Yangtze River beginning at~10 Ma (Ouimet et al., 2010). The riverbed profile of the middle Mekong River also records a rapid denudation event that also occurred at 17 Ma, although this incision signal was interpreted as a climate-driven event (Nie et al., 2018).
The above findings all point to rapid regional denudation and river incision since the Early Miocene. The timing is generally comparable with our results. Based on our knickpoint studies, we believe that this event should be more tectonically related, i.e., the tectonic uplift at the southeastern edge of the plateau in the Miocene led to strong river incision due to the relative lowering of base level of erosion. Clark et al. (2005) recovered the palaeo-topography of the southeastern margin of the plateau using planation surfaces and suggest the plateau surface had uplifted~3 km since the late Miocene by the isostatic rebound.
Furthermore, this Miocene tectonic uplift on the southeastern plateau margin is also widespread through the Tibetan Plateau margin. For example, previous studies suggest that the uplift of the northeastern, northwestern and central Tibetan Plateau and the Pamir Plateau occurred at 25-20 Ma (George et al., 2001) and the West Kunlun since Early Miocene (Jiang and Li, 2014;Li et al., 2019). The plateau expanded to the east and north and reached a height and width close to its current state due to a significant uplift at 15-10 Ma (Hough et al., 2011;Wang et al., 2012a). The low-temperature thermochronology, sedimentary, and tectonic deformation records show that main tectonic uplifts and denudations in different parts of the plateau are quasi-simultaneous and can be divided into four stages, respectively 60-35, 25-17, 12-8 (17-12 Ma in southern Tibetan), and 5 Ma . Another comprehensive study also shows that the uplift of the Tibetan Plateau was episodic and can be divided into five stages, reaching an average altitude of 500 m at 60-50 Ma,~1000 m at 45-34 Ma,~3500 m at 25-17 Ma, and a maximum height at 12-8 Ma (Zhang et al., 2013;Cao et al., 2022). All of the above studies indicate that the plateau margins have experienced two-staged tectonic uplifts since the Miocene Zhang et al., 2013), roughly corresponding to the early uplift (21-24 Ma) and late uplift (7-13 Ma) revealed by our study. The magnitude of the two uplifts estimated by previous studies is about 3000 m (Zhang et al., 2013), which is close to that estimated in this study (2500 m). The above comparison suggests that the two-staged uplifts on the southeastern margin of the Tibetan Plateau can roughly represent the uplift process of the entire plateau margins.
This result is also corroborated by paleoclimate studies of the evolution of the East Asian summer monsoon which is sensitive to the uplift of the Tibetan Plateau. The monsoon did not exist in northern East Asia until the plateau reached half of its present height, and was continually strengthened with plateau uplifting, and the aridity to the north of the plateau was simultaneously intensified (Liu, 2000). Intensified aridity will lead to an increase in dust release and a corresponding increase in dust deposition in the downwind region. Studies show that dust deposition experienced a stepwise expansion in the Loess Plateau, from the west at~25 Ma to the center at~20 Ma, and to the east at~7 Ma (Meijer et al., 2021). The episodic expansion can generally be correlated with the two-staged uplifts found in our previous study. The~20 Ma expansion was once associated with the uplift of the Tibetan Plateau (Guo et al., 2008), and the~7 Ma expansion was thought to be due to the establishment of winter monsoon circulation and inland aridification, both of which are associated with the uplift of the Tibetan Plateau .

Conclusion
The eight tributaries of the lower Jinsha River are located in the transition zone between the Tibetan Plateau and the Sichuan Basin and take the Sichuan Basin as their erosion base. Slope-area and integration analysis shows that there are two levels of slope-break knickpoints associated with changes in the erosion base on these rivers, with an average elevation of 2865 and 1523 m, respectively. The reconstruction of the longitudinal profiles of the rivers shows a drop in the erosion base of 1400 and 1100 m, respectively. The two-staged uplift events in the studied region on the Tibetan Plateau are estimated to be 21-24 and 7-13 Ma, respectively, with reference to the longterm downcutting rates of the Jinsha River obtained by previous studies. In addition, the numerically inverted river longitudinal profiles also reveal two similar uplift events. Thus, the uplift of the studied region is divided into two stages, the early stage (~20 Ma) when the plateau uplifts from 1500 to 2900 m, and the late stage (~10 Ma) when the plateau uplifts from 2900 to 4000 m. This finding confirms that the stepwise uplift of the Tibetan Plateau proposed by previous studies also occurs at the southeastern edge of the plateau.

Data availability statement
Publicly available datasets were analyzed in this study. This data can be found here: https://search.asf.alaska.edu/.

Author contributions
RP: collected DEM data, made calculation and simulation and prepared the draft. ZH: designed this study and revised the draft. QS: revised the draft and contributed some calculation methods. GL: Frontiers in Earth Science frontiersin.org contributed some ideas and revised the draft. YL: collected the geological data and drew the geological map. XL: contributed some ideas and revised the draft. XW: contributed some ideas and revised the draft.

Funding
This study was funded by the National Natural Science Foundation of China (Grant No. 42021001).