Active Tectonics Revealed by River Profiles along the Puqu Fault

The Puqu Fault is situated in Southern Tibet. It is influenced by the eastward extrusion of Northern Tibet and carries the clockwise rotation followed by the southward extrusion. Thus, the Puqu Fault is bounded by the principal dynamic zones and the tectonic evolution remains active alongside. This study intends to understand the tectonic activity in the Puqu Fault Region from the river profiles obtained from the remotely sensed satellite imagery. A medium resolution Digital Elevation Model (DEM, 20 m) was generated from an Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) stereo pair of images and the stream network in this region was extracted from this DEM. The indices of slope and drainage area were subsequently calculated from this ASTER DEM. Based on the stream power law, the area-slope plots of the streams were delineated to derive the indices of channel concavity and steepness, which are closely related to tectonic activity. The results show the active tectonics varying significantly along the Puqu Fault, although the potential influence of glaciations may exist. These results are expected to be useful for a better understanding of tectonic evolution in Southeastern Tibet.


Introduction
Empirical observations have indicated that steep topography may be closely related to rapid uplift [1].In particular, several studies have revealed that landscape surfaces are efficient indicators related to tectonic uplift and deformation processes [2][3][4].In addition, some studies argued that hillslopes may be influenced by potential balance between erosion and soil production rates [1,[5][6][7].These tectonic activities may alter river network mainly through changes of fluvial incision rates and stream diversion by modifying erosion base level [8][9][10].The river network, on the other hand, is continuously responding to tectonic forcing and accordingly provides potential meaningful information about spatial variation of uplift across landscape [1,[11][12][13][14].Furthermore, the drainage network may continuously adjust to changes of surface slope, and consequently it shows potential in recording information related to evolution of faults and folds [15,16].
Empirical observations from river systems across the earth in past decades also denoted a consistent power-law relationship between channel slope and drainage area.Theoretical arguments for both detachment-limited and transport-limited erosion regimes strongly suggested that rock uplift rate exerts first-order control on this relationship [17,18]: S indicates local channel slope, A is upstream drainage area and ks and θ are channel steepness and concavity indices, respectively.This power-law relationship has been widely used in past studies, and its usefulness has been verified in both detachment-limited and transport-limited river systems.For example, Snyder et al. [19] have analyzed landscape responses to tectonic forcing in the Mendocino Triple Junction region in Northern California; Lague et al. [20] have estimated uplift rate and erodibility using the area-slope relationship in Brittany, France; Kirby and Whipple [21] have successfully quantified differential rock uplift rates through the stream profile analysis in the area of Siwalik Hills in Central Nepal; Schoenbohm et al. [4] have calculated the surface uplift in the Red River region in Yunnan Province of China; Shahzad et al. [14] employed Shuttle Radar Topography Mission (SRTM) data to extract draining network for regional tectonic analysis in Potwar Plateau in Northern Pakistan; a similar approach was utilized by Shahzad et al. [22] to understand tectonic geomorphology in Kaghan Valley in Northern Pakistan; Anoop et al. [23] suggested the existence of large-scale tectonic deformation from the calculated geomorphic indices across the Kaurik-Chango Fault region in the Upper Spiti Valley, Northwestern Himalaya; Fuchs et al. [24] estimated the impacts of tectonics and climate change on the Panj River in Pamir Plateau through an analysis of fluvial incision and stream profile of the Panj river network.These studies have revealed a direct power-law relationship between the steepness index ks and the rock uplift rate U.
Since the channel concavity θ and the steepness ks indicate spatial variation of rock uplift rate, river profile parameters derived from regressions on logarithm slope-area data allow to extract information about regional tectonics.In particular, the analysis of these channel gradients is often on the basis of Digital Elevation Model (DEM) [25].Several methods have been proposed to extract river profiles and channel parameters from DEM.For example, Snyder et al. [19] and Kirby and Whipple [21] have developed a group of built-in functions in ARC/INFO to create flow accumulation arrays and delineate drainage basins; Mahmood and Gloaguen [26] calculated the indices of relative active tectonics (IRAT) in the Geography Information System (GIS) from a DEM in Hindu Kush region.In addition, scientific programming languages such as Matlab and SAS scripts were integrated in GIS to color-code stream profiles by their steepness and concavity indices: Shahzad et al. [22] have developed the package TecDEM using MATLAB for drainage network partition and river profile analysis from DEM and this package was also used by Vezzoli et al. [27] to extract the fluvial network in Western Greater Caucasus.Other widely-used packages able to extract river profiles from DEM include LANDLORD [28] and Terrain Analysis System (TAS, [29]).A recent study attempted to use LiDAR (Light Detection and Ranging) to generate a DEM for river profile extraction [30]; however, the cost of LiDAR is very high, which may decrease its practical use in the field.
This study intends to understand the tectonic activities in the Puqu Fault Region by using the power-law relationship between channel slope and drainage area.The topographic indices of river profile shape and character were derived from a freely available Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) DEM.Past studies have reported that in the boundary of the Southeastern Tibet Plateau, rivers are responding to tectonic uplift and deformation, recording temporal and spatial differences in rock uplift rates [4,31].The river network in the Puqu Fault area, which is situated in the southeastern margin of the Tibet Plateau with numerous tributaries flowing into it, is expected to be a sensitive recorder of the spatial and temporal variations in active tectonics in the eastern Karakoram-Jiali Fault Zone (KJFZ).This study chiefly concentrates on extracting tectonics from topography by analyzing river elevation profiles for these tributaries flowing into the Puqu Fault on the basis of the detachment-limited stream power model, in particular, to understand the tectonic activities in the Puqu Fault region from a series of river profiles that were extracted from the ASTER DEM.

Study Area
The Puqu Fault region is located in the KJFZ in Southeastern Tibet.This region is situated on the boundary of China, India and Myanmar (Figure 1).Due to the topography and remoteness, this area has been so far seldom investigated.

The Karakoram-Jiali Fault Zone in Southeastern Tibet
Southeastern Tibet is an area very limitedly studied especially for tectonic evolution, despite that its importance has already been indicated by several studies.Booth et al. [32] claimed that an important element of understanding the tectonics of Southeastern Asia is by tracing the assembly of Tibet itself.To a larger extent, some studies indicated that Tibet is generally composed of terraces accreted successively to the southern boundary of Eurasia (e.g., [33][34][35]).Significant features of tectonics can be readily observed in Southeastern Tibet.These observations generally suggested a complicated series of tectono-magmatic events shaping this region though more studies are fundamentally needed to further clarify these tectonic activities.
As indicated in Figure 2, a mixture of normal and strike-slip faulting may summarize those active tectonics in Southern Tibet [34,36,37].The normal faults are oriented in the direction of North-South and accordingly may reflect an East-West extrusion.Those North-South faults are widespread despite that the active Indian Plate is moving northward with reference to the relatively dormant Eurasia Plate at a rate about 4 cm/year [38].A typical strike-slip system in Tibet is located in the KJFZ.It is composed of a group of Northwest-Southeast aligned, right-lateral faults which appear to form the normal fault system in Southern Tibet [36].In the East, the KJFZ is primarily bounded by the Jiali Fault Zone.The movement of the Jiali Fault is tightly associated with the clockwise rotation and the extruded body of Eastern Tibet, and is closely connected with the northeastward indentation of the India Plate, and such activity results from the indentation of the India Plate into the Eurasia Plate as well as the related erosion in the Eastern Himalaya Syntax [39].

The Puqu Fault Region
Figure 3 indicates the drainage pattern of the Puqu Fault region.In the North of the Eastern Himalayan Syntaxis, the Jiali Fault extends to Southeast and separates into two main easternmost branches, namely the Puqu Fault in the South and the Parlung Fault in the North, which extends southeastward to the Gaoligong Fault [40].This study mainly focuses on the southern branch, the Puqu Fault.As an area of the easternmost Jiali Fault, the Puqu Fault is located in a principal active fault zone and is influenced by the eastward extrusion of Northern Tibet as well as the clockwise rotation followed by the southward extrusion.Literature is very limited about the tectonic evolution in the Puqu Fault region.Chung et al. [39] concluded that during 23-19 Ma, the northern branch of the Jiali Fault, the Parlung Fault (Figure 2), had undergone intensive strike-slip activities.In addition, Chung et al. [39] argued that an uplift that occurred during 10-6 Ma was related to an event which happened in the region of Himalaya as well as Eastern Tibet.The uplifting is then considered to be a regional event that may result from the indentation of the India Plate into the Eurasia Plate.In a more detailed study conducted by Lee et al. [36], rock samples from the Puqu Fault were collected via a Northeast-Southwest transverse, along which a large amount of sheared granitoids was discovered.Lee et al. [36] showed that these rocks always contain garnet or sillimanite, and corresponding peak metamorphic temperatures and pressures are considered to be higher than 550 °C and around 3-5 kbar.The triangle-shaped Indochina continent is featured by the Cenozoic convergence between the Indian and Asian Plates that developed the Indian-Asian continental collision [41][42][43].
So far, the Puqu Fault region has been narrowly studied due to the extreme inconvenience of field work, for instance, difficult to perform dating technology such as fission track in this remote area.As a result, the remote sensing is one important approach to investigate this area.In this study, the ASTER image was utilized to extract the tectonic indicators in the Puqu Fault region, particularly by analyzing river profiles extracted from the derived DEM that was generated from a stereo pair of images.

DEM Generation
The 30 m ASTER Global Digital Elevation Model (GDEM) is freely available, however in the Puqu Fault region its quality is not good due to some unfilled depressions.In this study the ASTER image was utilized to generate DEM, considering its low cost with decent spatial resolution (15 m).In addition, ASTER offers a back-looking stereo visible and near-infrared (VNIR) telescope with 15 m resolution in Band 3.This viewing geometry enables to generate a DEM with horizontal spatial resolution of 15 m.A DEM was created from a stereo pair of ASTER images acquired on 3 January 2005.Using the image correlation, it is possible to match pixels in these two images and then use the sensor geometry from the computed model to derive x, y, and z positions.The original ASTER images are Level-1A data products which are formally defined as reconstructed, unprocessed instrument data at full resolution.The scene center of the data is 28°47'47" N and 96°48'11" E with the sun elevation angle of 35° and the sun azimuth angle of 158°.
One hundred and eighty-seven pairs of tie points distributed throughout the whole image were identified in order to create an epipolar image pair of the channel 3N and 3B.The total Root Mean Square Error (RMSE) in x direction is 0.53 pixel whereas that in y direction is 0.19 pixel.The result is acceptable considering the difficulty in collecting tie points in this steep mountainous area of the Puqu Fault region.With the input of the ASTER orbital information in the software of PCI Geomatica 10, the epipolar images of a stereo pair of 3N and 3B were extracted.After obtaining the epipolar images, the DEM was extracted from the overlap extended between the epipolar image pairs.The final step was to geocode the epipolar DEM and the output was a DEM reprojected to the ground coordinate system.In this study, a DEM with a sampling interval of 20 m was generated.The final derived DEM is shown in Figure 4. Furthermore, on the basis of the derived DEM, the shaded relief map (Figure 5) was displayed in order to get a preliminary geomorphologic understanding in the Puqu Fault region.

River Network Extraction
The streamlines and river network were extracted from the generated DEM for further calculation of channel parameters.Streamlines define the flowing direction at each point on the topographic surface, starting on drainage divides, tracing the downhill water flowing path, and crossing contour lines.In terms of streamlines tracing, the D8 algorithm [44] was used to approximate the flow directions on the topographic surface and to track the flowing way from each pixel to its eight neighboring pixels.The direction from each pixel to its eight neighboring pixels was defined on the basis of steepest descent.This D8 algorithm is suitable for identifying individual channels, river networks and basin boundaries based on two assumptions: (1) the use of eight discrete flow angles; and (2) each pixel has a single flow direction.To simply put, these two assumptions restrict single stream direction into downstream neighboring pixel.In order to establish a D8 flow grid, the imposed gradients method was utilized to resolve the flow direction in flats that contain pixels not having a neighboring pixel with lower elevation.This method mainly focuses on centering flows within flat valleys and reducing parallel flows.Moreover, the D8 flow grid was created from raster into vector.All pixels that have a flow code of zero were treated as basin outlets, including four edges of the DEM and pixels containing no data within the DEM.The derived vector file stores data for a single basin or several disjoint subbasins.Every pixel in a particular basin is the outlet pixel for a subbasin that is contained in it.Each of these subbasins records several attributes, such as contributing area and relief.This vector file stores all of these attributes for all of the pixels/subbasins.Furthermore, two vector files were created in the software River Tools 3.0.One was the link file that stores attributes for every channel link in the river network, and the other was a stream file that records attributes for every Horton-Strahler stream in the river network.The Strahler Stream Order was employed to define the threshold when extracting the river network.In this study, the Strahler Order of the main river along the Puqu Fault was set to 5. The extraction threshold was defined as the Strahler Order ≥ 3. The derived river network is shown in Figure 6.However, when subsequently calculating the concavity and steepness index of the extracted streams, channels with Strahler Order ≥ 4 were selected to extract the main stem and its tributaries.Finally, all streams with the Strahler Stream Order ≥ 4 were exported into dBASE Table files (.dbf).The dBASE Table files were then imported into EXCEL and MATLAB for further river profile analysis.

River Profile Analysis
River incision is often modeled as competition between uplift and erosion activities with change rates of channel elevation [19,21,45].In the model, river profile evolution is simulated as an equation: dz/dt represents the temporal rate of changes of channel elevation, U indicates the rock uplift rate with respect to a fixed base level and E refers to erosion rate.The erosion process can be described as an empirical relation with the expression: S and A represent the local channel gradient and the upstream drainage area, respectively.K is identified as a dimensional coefficient of erosion whereas m and n are positive constants related to basin hydrology, erosion process and hydraulic geometry [19,21,31,45,46].Combining Equation (2) with Equation ( 3), the incision model can be derived as: In the cases of steady-state (dz/dt = 0), with uniform U and K as well as constants m and n, Equation ( 4) can be further expressed as the equilibrium channel gradient: Equation ( 5) describes the relationship between channel gradient and drainage area.It is often simplified as a power-law function as θ s S k A − = previously mentioned in Equation ( 1) with the coefficient ks representing the channel steepness, equivalent to (when combining Equations ( 1) and ( 5)): Equation ( 6) is often recognized as a strong empirical support for a positive correlation between the steepness index ks and the rock uplift rate U (e.g., [19,21]).The coefficient θ, which represents the intrinsic channel concavity, replaces the ratio m/n, written as: Equation (7) shows that channel concavity θ is related to the constants m and n.What is more important, as for m and n, it is very difficult to obtain them from natural field experiments.Channel concavity θ, namely, the ratio of m/n, is easier to be estimated [31,47].Although the concavity value is theoretically expected to range from 0.35 to 0.6, the observed concavity may be higher considering non-uniform uplift, non-uniform lithology or variations in sediment flux and grain size [4,31].
After extracting the river network in the Puqu Fault region, the next step is to derive the channel parameters such as concavity and steepness from stream elevation profiles.In this study the concavity and steepness were calculated directly from regressions on logarithm slope-area plot.Based on Equation (1), when taking logarithm on both sides of the equation, the equation can be written as: and after reforming the Equation ( 8), the equation turns into: Equation ( 9) indicates that it can get the information of concavity and steepness by putting slope S and drainage area A in the log-log plot and then drawing the linear regression fit line.
Using the generated 20 m resolution ASTER DEM, stream elevation profiles were extracted for 54 tributaries to the Puqu Fault in Southeastern Tibet, with the defined Strahler Order ≥4.Among them, six streams were eliminated because of large pitfalls due to the quality of DEM which cannot be filled.By analyzing the elevation profile and the area-slope plot of each stream, it shows that the channels can be generally defined as two morphologic types (Figure 7).The first one is shown in Figure 7a.It represents the stream containing only one lower concave segment in the river profile and the upper segment is absent or hardly recognized.The downwards linear fit can be identified in the area-slope plot.The second type, as shown in Figure 7b, indicates the stream consisting of two concave channel segments.The upper segment and lower segment are separated by convex transition zones, namely knickzones.The corresponding two downwards fits can be found in area-slope plot separated by data with upwards trend.The channel concavity and steepness were calculated for each upper segment with minimum length of 15 pixels (300 m) and each lower segment with minimum length of 25 pixels (500 m) based on the regression fit.The Root Mean Square Error (RMSE) of each regression was calculated and regressions with RMSE greater than 0.1 were not taken into consideration.The overall regression result for all streams is shown in Table 1.The corresponding location of the streams can be found in Figure 8.The streams 13, 29, 30, 32, 46 and 52 which contain depression areas were eliminated, marked as N/A.Among the remaining 48 rivers, 32 of them were detected with both upper and lower concave segments.The remaining 16 streams were only found with the lower concave segment.The channel parameters of convex segments were not calculated because they may indicate transiting knickzones and have not formed the new equilibrium response yet.Alternatively, convexities of river profiles may be interpreted as a transient response of the river system to the change in tectonic uplift.For each linear regression fitting line, based on the Equation ( 9), the slope of the regression fit is −θ whereas the intercept of y axis is s k log .The results of concavity and steepness are shown in Table 1.

Channel Equilibrium Hypothesis
The concavity and the steepness were calculated based on the stream power law (Equation ( 1)).In spite of the strong empirical support, the stream power-law incision model, as well as the interpretation of model parameters ks and θ in Equations ( 6) and (7), is on the basis of the equilibrium channel hypothesis.Hence, the state of activity of analyzed streams can be evaluated and the usefulness of utilizing the stream power-law function in the Puqu Fault region can be anticipated.First, the uplift rates in the Tibet Plateau have been reported constantly over a long period of time.Although Spicer et al. [48] argued that the elevation of the Southern Tibet Plateau has remained unchanged since 15 Ma ago, other studies reported a rapid but steady uplift in the Tibet Plateau from the late Tertiary and Quaternary (e.g., [49,50]).Even for the southern Tibet fault system, no modern activity was detected [51].Such timescale as the later Tertiary and Quaternary, is long enough in comparison with the landscape response time, and this sufficiency makes the hypothesis of steady condition plausible [52].Secondly, the equilibrium process is continuously self-adjusting to accord with the migration of knickzones.In consequence, the new condition is always being created for the maintenance of steady state, thus forming a new equilibrium slope as the response.Because this elaboration may not be sufficient to ensure the feasibility of stream power-law function, in order to maximally avoid the misuse of the stream power-law model, all convex parts of the stream profiles were excluded, considering that under certain circumstances the large scale convexity is always related to the transition knickzone.Limited focus on those concave profiles may maximally eliminate the possibility of transient state.Thirdly, although for some stream profiles knickzones may exist, the overall profiles of the Puqu Fault are relatively smooth and without abrupt knickpoints which are always recognized as the obvious indications of transient state.That is, the knickzone separates two equilibrated profiles: upper segments adjusted to initial conditions and lower segments adjusted to new conditions.This is in accordance with the theory that fully-attained equilibrium stream is always without major knickpoints because they generally appear with changes of uplift rate or erosion rate, and knickpoints can also result from non-uniform uplift or changes in lithology.Fourthly, despite the assumption of existing steady state, it is not plausible to assume that the model parameters can only be derived by analyzing equilibrium state data.The reason is that some data may contain both steady and transient information.And for high resolution data, by analyzing elevation profile and area-slope log plot, it is possible to differentiate the steady and transient states by analyzing knickpoints and knickzones.Sometimes it is even more important to analyze such data because it provides the possibility to detect two steady-state profiles [1].Based on the above four reasons, to conclude, in this study it is practical to utilize the power-law stream incision model as indicated in Equation (1).

Potential Effects of Glaciations
Hobley et al. [53] indicated that in Indian Himalaya, glacial modification of the upper reaches of a catchment may have strong influence on the hydraulic scaling of the channel downstream, which may significantly increase the channel concavity indices of downstream reaches exceeding the expected value range of 0.3-0.6.Whether the Puqu Fault is under influence of large glaciations as appeared in other high relief places in Tibet remains an important research consideration.To attempt to understand this question, the glacier map from the Chinese Glacier Information System (CGIS) established in the year 2004 [54] was acquired.The glacier distribution in the Puqu Fault region is shown in Figure 9.As the inventory map indicates, the glaciers existing along the Puqu Fault are not as much as those in the north in the area of the Parlung Fault and these glaciers are distributed in the higher elevation area, namely the upper segment part.Although the existence of glaciers is supposed to smooth the valley and to erode the valley walls-They can also reduce relief and form plateaus if they are ice caps rather than valley glaciers-Only relatively higher concavity and steepness were observed for the stream 11.The derived concavity value is 2.70, a bit lower than the lower segment value 2.86 and much higher than the average concavity value (0.84) of those upper streams in the north of the bend.The calculated s k log value for the upper segment is 19.27, also slightly lower than the lower segment value (21.35) but much higher than the average value (5.75) of the north bend area.Other than the stream 11, no significantly higher concavity and steepness values was noticed for those streams affected by glaciations in the upper concave segment.Thus, it may be argued that the influence of glaciations is simply limited along the Puqu Fault.Furthermore, it could be argued that the effects of past climate changes may facilitate the mountain uplift, thus resulting in the high concavity and steepness indices in the Puqu Fault region because of the post-glacial rebound effect related to the retreat of existed glaciers.However, although the theory of the post-glacial rebound is proven to be true and was verified by some studies in Europe and North America several decades ago (e.g., [55,56]), such argument in the Tibet Plateau remains to be further investigated.Although some studies proposed the hypothesis of "ice sheet" suggesting that the former ice sheet may result in the glacial uplift in the Tibet Plateau [57], it is contradictory to the new findings that the Tibet Plateau is a young plateau rather than an ice covered plateau [58].Also, such hypothesis has been denied by the zero-order test of the purely atmospheric response to Tibetan ice sheet forcing [59].Even without isostatic rebound, the ice sheet may have created the low relief topography which can be depleted by incision, thus potentially causing the two segments: the top segment with low concavity and slope on the relict plateau and the lower segment with high concavity and slope in the incised part of the landscape.Therefore, it may be argued that the uplift in the Puqu Fault region is not due to the rise of land masses which were previously depressed by the huge weight of ice sheets during the past time.That is, the high concavity and steepness values of the land in Southeastern Tibet are not likely to be the consequence of post-glacial rebound.Therefore, river morphologies and profiles tend to be excluded from the significant influence of the potential glaciations.

Open Questions
In the Puqu Fault region, the topography of rivers is usually very steep.As documented above, most of the channel data have slopes in excess of 0.1 and commonly in excess of 0.3.The transition from fluvial channels to steep debris-flow channels is usually reported at slopes between 0.1 and 0.3 [1].In consequence, it may be argued that fluvial processes are not responsible for the observed shape of the profiles.In addition, the analysis in this study was performed on relatively short segments, though the potential use of the stream power law θ s S k A − = is usually performed segments over a range of drainage areas.Some studies use a reference concavity index to calculate a normalized steepness index at every point along their river profiles (e.g., [25,60,61]).Although normalized steepness is useful for short channel segments, concavity is more meaningful to characterize the region where uplift rate may have systematic change along the profile.

Conclusions
The Puqu Fault is located in a tectonically dynamic region in Southern Tibet.This study attempted to understand the tectonic activities in the Puqu Fault region from a series of river profiles that were extracted from an ASTER DEM.The main work was composed of four parts: (1) DEM generation from two ASTER stereo pair images; (2) river network extraction from the ASTER DEM; (3) calculation of channel parameters; and (4) interpretation and discussion of the derived results.First, in terms of DEM generation for the Puqu Fault region, a total of 187 GCPs were manually chosen with a decent RMSE: 0.53 pixel along the x direction and 0.19 pixel along the y direction.The epipolar images were subsequently produced.The DEM was finalized with 20 m resolution after being resampled from the original 15 m resolution for the purpose of noise reduction.Secondly, the river network was extracted for the Puqu Fault region on the basis of the D8 algorithm with a pre-defined threshold of the Strahler Order ≥ 4. Thirdly, the calculation of concavity θ and steepness ks was based on the stream power law θ s S k A − = .A series of 48 stream elevation profiles on both sides of the Puqu Fault were extracted and then the channel slope S and drainage area A of each profile were calculated.
The study further put slope and drainage area in the logarithm plot and derived the concavity and steepness index through the linear regression.Among these 48 streams, 32 of them were detected with both upper and lower concave segments separated by intermediate transition zones, whereas the remaining 16 streams were detected only with the lower segments.For each segment, the channel parameters of concavity θ and logarithm steepness s k log were calculated.This paper further discussed the river profile analysis which may be influenced by the following factors: (1) the usefulness of the stream power law; (2) the potential effects of glaciations; and (3) the steep topography and relatively short channel segments in this region.Currently the understanding of the tectonic activities in the Puqu Fault region is largely limited by the difficulty in accessing this area.In the future, a field trip in Southeastern Tibet area would be very beneficial for geological surveying and result validation across the Puqu Fault region.

Figure 1 .
Figure 1.Location of the Puqu Fault Region.The red box shown in the upper map (information available in Google Earth) shows that the study area is located on the boundary of China, India and Myanmar.The lower map is from the ASTER image with the 3-2-1 (red-green-blue) color combination, showing the main coverage of the Puqu Fault and its surrounding area.The ASTER Level 1A data (raw digital numbers quantized as 8-bit unsigned integers) was acquired and the date of image acquisition was 3 January 2005.

Figure 2 .
Figure 2. The topography and the principal active faults in Tibet and its adjacent regions.Here white arrows indicate the current movements of India and Northern Tibet and the ongoing clockwise rotation followed by the southward extrusion.(a) The simplified geologic map of the study area; (b) The tectonic map showing the activities of the Himalayas and Southern Tibet.Reprint with the permission of [36], 2003 copyright Elsevier).

Figure 3 .
Figure 3.The drainage system in the Puqu Fault Region.Different colors indicate diverse subbasins extracted.

Figure 4 .
Figure 4.The digital elevation model (DEM) derived from the ASTER stereo pair images for the Puqu Fault region.

Figure 5 .
Figure 5.The shaded relief map of the Puqu Fault region with the simulated sun compass angle of 135° and the zenith angle of 40°.The red arrow indicates the bend of the Puqu Fault.

Figure 6 .
Figure 6.The channels were vectorized and exported in 3D view with DEM.The streams with Strahler Order < 3 were excluded.The elevation range is from 2150 to 5040 m.The approximate size of the domain is about 1400 km 2 .

Figure 7 .
Figure 7. River elevation profiles and slope-area plots extracted from a 20 m ASTER DEM for two representative tributary types.(a) The example of a stream with one concave channel segment; (b) The example of a stream with two concave channel segments, one in the upper segment and another in the lower segment.The lines represent the linear regression fit.The arrows indicate the range of each concave segment.

Figure 8 .
Figure 8. Location of the streams marked with the numbers from 1 to 54.

Figure 9 .
Figure 9. Distribution of glaciers in the Puqu Fault region as shown in the glacier map acquired from the Chinese Glacier Information System (CGIS).