We investigate segmentation of High Himalayan strain by cross-orogen structures separating western and eastern obliquely convergent sectors from a central orthogonally convergent sector, and evaluate the relationship between the size of regions accumulating strain, their proximity to the toe of the thrust wedge, and recurrence of Mw >7 earthquakes. We present a map of river channel steepness (ksn)—a proxy for rock-uplift rate over 105 yr, for the Himalayan arc—and evaluate the strength of its correlation with Main Himalayan thrust (MHT) coupling (–0.6), earthquake density (0.6), topography (0.6), lithotectonic units (0.5), and precipitation (–0.3) along 40 profiles spanning the Himalaya from 78°E to 92°E. We interpret the ksn map to be foremost a function of recent strain accumulation. This reveals prominent offsets of hinterland strain accumulation collocated with cross-orogen strike-slip and extensional fault systems. Clusters of high-ksn rivers are located near the boundary between the strongly and weakly coupled portions of the MHT, where fault behavior changes from seismogenic to sliding at the rheologic brittle-to-plastic transition (BPT). We propose that the rate at which major MHT earthquakes repeat is related to four parameters: convergence rate (nearly uniform); spatial dimensions of the high-ksn cluster (proxy for volume of material accumulating strain); the high ksn clusters distance from the toe of thrust wedge (fault surface area over which static friction must be overcome); and the degree of obliquity between India-Asia convergence and the local trend of the orogen (proxy for the magnitude of strain partitioning).

The Himalayan thrust wedge is 2500 km long, 250 km wide, and up to 50 km thick below southern Tibet (Nábělek et al., 2009) (Fig. 1). The wedge is floored by an intracontinental thrust fault separating India from Asia known as the Main Himalayan thrust (MHT), whose behavior gradually changes from an aseismically sliding décollement beneath southern Tibet to the Himalayan seismogenic zone across a 50–100-km-wide brittle-plastic transition (BPT) beneath the High Himalaya (Bilham et al., 2001; Ader et al., 2012; Cannon and Murphy, 2014; Stevens and Avouac, 2015) (Fig. 2). The BPT corresponds to the boundary between the strongly coupled portion of the MHT, which only moves during earthquakes, and the weakly coupled portion (Stevens and Avouac, 2015). Clusters of steep river channels are located at the BPT (Fig. 2). Clusters are offset from each other across strain partitioning structures that have been independently proposed to separate the range into orthogonally convergent and obliquely convergent sectors (Drukpa et al., 2006; Velasco et al., 2007; Styron et al., 2011; Murphy et al., 2014) (Fig. 3).

Figure 1

(A) Digital elevation map of the Himalaya with major rivers, peaks 8000 m or higher, and bedrock geology. Geologic map is after DeCelles et al. (1998), Long et al. (2011), Mitra et al. (2010), Yin (2006), Webb et al. (2007), and He et al. (2015). Gray polygons are Indian basement ridges, after Gahalaut and Kundu (2012): DHR—Delhi-Hardwar ridge; FR—Faizabad ridge; MSR—Munger-Saharsa ridge. Tibetan rifts: AD—Ama Drime; CR—Coma rift; GM—Gurla Mandhata; LP—Leo Pargil; TG—Thakkhola graben; YG—Yadong-Gulu rift. SP—Shillong Plateau; K—Karnali river; A—Arun river. (B) Geologic cross sections, after Cannon and Murphy (2014), Robinson et al. (2008), Bollinger et al. (2004), and McQuarrie et al. (2014). There is no vertical exaggeration. See A for cross-section locations. Faults: MHT—Main Himalayan thrust; MFT—Main Frontal thrust; MBT—Main Boundary thrust; MCT—Main Central thrust; STD—South Tibetan detachment.

Figure 1

(A) Digital elevation map of the Himalaya with major rivers, peaks 8000 m or higher, and bedrock geology. Geologic map is after DeCelles et al. (1998), Long et al. (2011), Mitra et al. (2010), Yin (2006), Webb et al. (2007), and He et al. (2015). Gray polygons are Indian basement ridges, after Gahalaut and Kundu (2012): DHR—Delhi-Hardwar ridge; FR—Faizabad ridge; MSR—Munger-Saharsa ridge. Tibetan rifts: AD—Ama Drime; CR—Coma rift; GM—Gurla Mandhata; LP—Leo Pargil; TG—Thakkhola graben; YG—Yadong-Gulu rift. SP—Shillong Plateau; K—Karnali river; A—Arun river. (B) Geologic cross sections, after Cannon and Murphy (2014), Robinson et al. (2008), Bollinger et al. (2004), and McQuarrie et al. (2014). There is no vertical exaggeration. See A for cross-section locations. Faults: MHT—Main Himalayan thrust; MFT—Main Frontal thrust; MBT—Main Boundary thrust; MCT—Main Central thrust; STD—South Tibetan detachment.

Figure 2.

Contour map of Himalayan normalized river channel steepness (ksn), a proxy for rock uplift of 105 yr, overlaid on the degree of Main Himalayan thrust (MHT) coupling (blue shading; darker blue = more coupling from GPS velocities [Stevens and Avouac, 2015]), active Himalayan-Tibetan faults (purple lines) after Styron et al. (2011) and Taylor and Yin (2009), and Mw 2–9 epicenters from A.D. 1966 to 2016 (black dots) from the U.S. Geological Survey National Earthquake Information Center (https://earthquake.usgs.gov/earthquakes/search/). Where rock uplift outpaces river incision the river channel steepens, ksn analysis uses a digital elevation model to calculate longitudinal river steepness normalized to a reference river channel concavity (0.45). For a more complete description see the methods section. Heavy white lines denote continuous clusters of high-ksn river channels referred to in text. Thin white lines are political boundaries.

Figure 2.

Contour map of Himalayan normalized river channel steepness (ksn), a proxy for rock uplift of 105 yr, overlaid on the degree of Main Himalayan thrust (MHT) coupling (blue shading; darker blue = more coupling from GPS velocities [Stevens and Avouac, 2015]), active Himalayan-Tibetan faults (purple lines) after Styron et al. (2011) and Taylor and Yin (2009), and Mw 2–9 epicenters from A.D. 1966 to 2016 (black dots) from the U.S. Geological Survey National Earthquake Information Center (https://earthquake.usgs.gov/earthquakes/search/). Where rock uplift outpaces river incision the river channel steepens, ksn analysis uses a digital elevation model to calculate longitudinal river steepness normalized to a reference river channel concavity (0.45). For a more complete description see the methods section. Heavy white lines denote continuous clusters of high-ksn river channels referred to in text. Thin white lines are political boundaries.

Figure 3.

Map of arc-parallel GPS velocity after Styron et al. (2011), overlaid on a contour map of Himalayan normalized river channel steepness (ksn), and active Himalayan-Tibetan faults (purple lines) after Styron et al. (2010) and Taylor and Yin (2009). Heavy black lines denote strain partitioning structures referred to in text.

Figure 3.

Map of arc-parallel GPS velocity after Styron et al. (2011), overlaid on a contour map of Himalayan normalized river channel steepness (ksn), and active Himalayan-Tibetan faults (purple lines) after Styron et al. (2010) and Taylor and Yin (2009). Heavy black lines denote strain partitioning structures referred to in text.

A record of Himalayan strain accumulation is encoded in the landscape through a competition between erosion and rock uplift, and can be visualized by mapping normalized river channel steepness (ksn), a proxy for rock uplift over 105 yr. Currently, high topography, high-ksn rivers, rapid denudation rates (Godard et al., 2014; Morell et al., 2017), rapid interseismic strain rates (from GPS velocities) (Ader et al., 2012; Bilham et al., 2001; Stevens and Avouac, 2015), orogen-parallel belts of seismicity (Bollinger et al., 2007; Pandey et al., 1995, 1999), and MHT thrust-ramp duplexes are all coincident with the inferred MHT BPT (Fig. 1B) (Pandey et al., 1999; Bollinger et al., 2004a; Herman et al. 2010; Cannon and Murphy, 2014; McQuarrie et al., 2014; Morell et al., 2015; Harvey et al., 2015) tracing the northern limit of the strongly coupled plate boundary (Stevens and Avouac, 2015) (Fig. 2). This suggests that the position of high topography within the thrust wedge is invariant over 105 yr time scales, and that it is a function of thrust wedge geometry, plate boundary thermal structure, and stored elastic strain. As described below, the strongest correlations found in our analysis are between ksn, MHT coupling (–0.6), and earthquake density (0.6) (excluding regions with too few modern earthquakes for statistical significance), and topography (0.6) highlighting that in active orogens, ksn records strain accumulation (modified by substrate erodibility, channel width, sediment load, and climate). All other parameters are correlated with ksn.

The MHT’s surface trace is manifest as a series of thrust-related folds and faults at the toe of the wedge, which form the first spatial, southernmost topographic expression of the Himalaya and are collectively known as the Main Frontal thrust (Fig. 1). The rate of thrusting across the Main Frontal thrust in central Nepal, documented by offset and folded river terraces (Lavé and Avouac, 2001), is within uncertainty of the rate of aseismic sliding recorded by GPS stations in south-central Tibet (Ader et al., 2012; Bilham et al., 2001, 1997; Jackson and Bilham, 1994). This has been interpreted by some researchers to signify that the wedge does not deform internally, implying that rock uplift within the wedge is a passive product of fluxing material over MHT ramps (Cattin and Avouac, 2000; Pandey et al., 1995). However, the suggestion of an entirely passive wedge is at odds with other observations. Recent work has highlighted that the A.D. 1505 earthquake in west Nepal may have ruptured multiple splays within the wedge (West Nepal fault system, Main Boundary thrust) (Hossler et al., 2016; Murphy et al., 2014; Harvey et al., 2015), while aftershocks of the A.D. 2015 Gorkha earthquake appear to have ruptured internal thrust wedge faults at the base of the High Himalaya (Whipple et al., 2016; Hoste-Colomer et al., 2016). Quaternary thrust faulting at the base of the High Himalaya has been proposed on the basis of field relationships and thermochronology (Hodges et al., 2004; Wobus et al., 2005), although evidence of recent surface rupturing thrusts at the base of the High Himalaya has been elusive. Our ksn map, as well as thermokinematic models of late Miocene to recent times and structural restorations of the frontal portion of the wedge are consistent with underplating of Indian crust along a foreland-migrating ramp (Bollinger et al., 2004a; Herman et al., 2010; McQuarrie and Ehlers, 2015; Robinson, 2008).

GPS measurements show that the most-rapid interseismic strain accumulation is located near or above (north of) the 3500 m topographic contour (Ader et al., 2012; Bilham et al., 2001), which has been used to model MHT coupling (Stevens and Avouac, 2015). We use this model to examine the relationship between river channel steepness and interseismic strain (Fig. 2). There is a nonlinear relationship between ksn and rock-uplift rate (Kirby and Whipple, 2001, 2012), where the time required to impart a ksn signal to the landscape is dependent on the long-term (multiple seismic cycles) rock-uplift rate. In the Himalaya, vertical interseismic strain rates are as much as 4–7 mm yr–1 (Jackson and Bilham, 1994; Ader et al., 2012; Grandin et al., 2012; Stevens and Avouac, 2015), however the 25 April 2015 Gorkha earthquake demonstrated that much of this interseismic uplift is elastic and destined to be transferred into horizontal and vertical motion of the strongly coupled portion of the range during plate boundary earthquakes (Elliott et al., 2016; Lindsey et al., 2015).

During the Gorkha earthquake there was a ∼0.5 m lowering of the High Himalaya and 1 m of uplift for the Lesser Himalaya (Elliott et al., 2016; Lindsey et al., 2015), showing that the long-term uplift rate of the High Himalaya may be a fraction of that measured by GPS, interferometric synthetic-aperture radar (InSAR), and spirit leveling. The map presented here documents clusters of high-ksn rivers, most of which straddle the boundary between the strongly coupled and weakly coupled portions of the MHT (Stevens and Avouac, 2015). In regions with a high rate of seismicity, clusters of earthquake epicenters are located south of (foreland to) the steepest river channels (Fig. 2). We follow Cattin and Avouac (2000) and Bollinger et al. (2004b) in interpreting peak seismicity to reflect the leading edge of the BPT and, based on our data, suggest that peak ksn reflects focused elastic strain in the upper crust and plastic deformation of the lower to middle crust (Cannon and Murphy, 2014). Here within the BPT, slip is fed from the plastically deforming lower crust beneath southern Tibet to MHT thrust-ramp duplexes in the mid- to upper crust via underplating or tectonic wedging. The MHT BPT is a mechanical boundary whose location within the thrust wedge is a function temperature along the subduction interface (MHT). The southernmost limit of the BPT, closest to the toe of the thrust wedge, begins where the MHT encounters the 350 °C isotherm (Herman et al., 2010), at which point quartz begins to deform plastically (Hirth and Tullis, 1992). The percentage of deformation accommodated by crystal-plastic creep increases with temperature as the dominant mineral constituents begin to deform plastically (feldspar begins to deform ductily at 600 °C) (Tullis and Yund, 1992). The width of the MHT BPT depends on the geothermal gradient and the local dip of the MHT. The thrust wedge north of the BPT acts as a dashpot (a device which resists motion via viscous friction), storing strain until it overcomes static friction on the strongly coupled portion of the MHT, generating a great earthquake (Bilham et al., 2001; Avouac, 2003; Feldl and Bilham, 2006; Cannon and Murphy, 2014). Either some amount of this strain is permanent or the high-ksn signal present across the High Himalaya represents a purely elastic bulge maintained by ongoing India-Asia convergence as suggested by Meade (2010). If true, then the highest topography of the Himalaya is dynamically maintained by plate convergence. While this idea is unconventional, an elastic strain reservoir equivalent to four Mw 8 earthquakes or a single Mw 9 earthquake has been inferred to exist across the range by comparison of the seismic energy released from all known great Himalayan earthquakes of the past 500 yr with convergence over the same time interval (Bilham and Ambraseys, 2005; Stevens and Avouac, 2016). While the percentage of accumulated strain released in plate boundary earthquakes versus absorbed by internal shortening structures remains uncertain (Bilham and Ambraseys, 2005; Cannon and Murphy, 2014; Meade, 2010; Stevens and Avouac, 2015; Taylor, 2016), the identical rate of frontal anticline growth and convergence of southern Tibet and India (Lavé and Avouac, 2001) indicates that the amount of strain available for internal deformation is small.

Bilham et al. (2001) estimated that up to 10% of the strain accumulated in the High Himalaya is available for crustal thickening. This slow uplift of the High Himalaya is evident in the ksn map, as clusters of high-ksn rivers in the High Himalaya are offset from each other across Himalayan-Tibetan rifts. Himalayan river channel steepness can be used as a proxy for strain accumulation and provides an intermediary time step between the decadal GPS record, the roughly 1000 yr historical record of great earthquakes, and the ±0.5 m.y. record of low-temperature thermochronology.

The surface bedrock geology of the High Himalaya is largely the result of emplacement of the metamorphic core of the orogen between the Main Central thrust and the South Tibetan detachment (Fig. 1) during the early Miocene (Searle and Godin, 2003; Edwards et al., 1999; Murphy and Harrison, 1999; Vannay et al., 2004; Thiede et al., 2006; Kellett et al., 2010). In central Nepal, these structures appear to maintain a low degree of activity today (Hodges et al., 2004; McDermott et al., 2013; Wobus et al., 2005), however no recent surface ruptures have been found, and cooling histories and river channel steepness data sets have been variously interpreted to represent out-of-sequence faulting, duplexing, or passive uplift over MHT ramps (Pandey et al., 1995; Cattin and Avouac, 2000; Avouac, 2003; Bollinger et al., 2004a; Wobus et al., 2005; Herman et al., 2010; McDermott et al., 2013, 2015; Hubbard et al., 2016).

In west Nepal, the South Tibetan detachment is folded (Fig. 1) and defines the Dolpo anticline, a regional structure 400 km long with at least 9 km of structural relief whose frontal limb coincides with a cluster of high-ksn channels interpreted to be within the BPT in west Nepal (Fuchs, 1964; Cannon and Murphy, 2014). Folding initiated after the cessation of South Tibetan detachment motion, regionally constrained to 16–19 Ma (Searle and Godin, 2003; Godin et al., 2006b). The concentration of high-ksn rivers and young zircon (U-Th)/He ages (8.0 ± 1.3 to 2.6 ± 0.7 Ma) (McCallister et al., 2014) at Gurla Mandhata (a core complex marking the western limit of the Dolpo anticline) (Murphy et al., 2002; Murphy and Copeland, 2005) suggests that hinterland folding initiated in the late Miocene and is ongoing today.

In Bhutan, the South Tibetan detachment was active 23–16 Ma (Kellett et al., 2010) or 20.5–14 Ma (Cooper et al., 2015) and appears to be folded into a broad anticline that spans the northern provinces of Punakha, Bumthang, and Kurtoed (Grujic et al., 2002; Long et al., 2011; McQuarrie et al., 2014). The anticline is collocated with the BPT interpreted from our ksn data. Low-temperature thermochronology of the northern Bhutan Himalaya is largely lacking from the literature, but Adams et al. (2013) presented apatite (U-Th)/He ages (2.4 ± 0.4 to 3.6 ± 1.2 Ma) from the southern edge of the northern Bhutan high-ksn cluster, which reflects recent exhumation of the Bhutanese hinterland. The anticlines overlying the BPT of west Nepal and Bhutan, interpreted from ksn data, record recent exhumation, and are bound on one side by major extensional structures, the Gurla Mandhata core complex and the Chomolari detachment (Yadong-Gulu rift) (Murphy et al., 2002; Wu et al., 1998). In both cases, folding appears synchronous with the onset of extension, suggesting that the initiation of hinterland thrust duplexes is kinematically linked with thickening of the lower crust in a constrictional strain field (Murphy and Copeland, 2005; Sundell et al., 2013).

Normalized river-channel steepness (ksn) analysis was conducted using the Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) Global Digital Elevation Model V002 with 30 m resolution from 78°E to 92°E. Although the ASTER digital elevation model has better resolution than previous models, it contains numerous holes and data issues. Deeply incised river canyons sometimes generate “digital dams”, which artificially elevate the river channel. Reaches of rivers with bad data are recognizable on longitudinal profiles (plots of elevation versus downstream distance) by their jagged stair-step appearance. Good ksn fits are not possible for these reaches, hence they were not analyzed. For this reason, deep canyons on along-trunk channels within the High Himalaya were commonly not able to be analyzed. To complete our ksn map in these regions, we analyzed tributary channels instead of trunk channels. This is the advantage of using normalized channel steepness, as it allows the user to quantitatively compare channels with different drainage areas.

Channel steepness analysis is founded on the observation that the longitudinal profiles of graded streams in steady state are characterized by gentle concave-up curves (Hack, 1957). In steady-state landscapes, channel slope, width, and available discharge are balanced by the available stream power (Davis, 1902; Mackin, 1948), implying that the rate of erosion is equal to the rate of rock uplift. This relationship can be expressed using:
graphic

(Kirby and Whipple, 2001; Wobus et al., 2006), in which S is slope, A is upstream drainage area, ksn is normalized channel steepness, and θref is reference concavity (0.45). The digital elevation model was hydrologically corrected in ArcGIS software (https://www.arcgis.com/) using the “fill” function in the hydrology toolbox to fill “digital pits” in the data. Then a flow-direction grid was created by determining the downslope direction at each cell using the “flow direction” function. Finally, using the flow-direction grid as an input, an accumulation grid was created by calculating the accumulated flow at each cell in the digital elevation model using the “flow accumulation” function. Once these hydrologic corrections were completed, the analysis could proceed. In ArcGIS, analysis parameters were selected for each session. In areas where the data quality is good, a 250 m smoothing window was used; in areas with poor data quality, a 500 m smoothing window was used. A 0.45 reference concavity was used for all sessions. Site-specific reference concavity can be calculated, but studies examining its effect on the analysis have all found that the choice of reference concavity changes the values obtained, but not the pattern of high and low channel steepness (Kirby et al., 2003; Whipple, 2004; Tyagi et al., 2009; Olen et al., 2015). Studies of ksn in mountainous regions use reference concavities between 0.4 and 0.5, with 0.45 being the most common choice (Kirby and Whipple, 2001, 2012; Kirby et al., 2003; Wobus et al., 2003, 2006; Duvall et al., 2004; Hodges et al., 2004; Korup, 2006; Clark and Bilham, 2008; Robl et al., 2008; Cyr et al., 2010; Allen et al., 2013; Adams et al., 2016; Harvey et al., 2015). After the parameters were set, a point on a river near its headwaters was selected, and its channel data (elevation, drainage area, and gradient at each point along the channel) were exported to Matlab software (https://www.mathworks.com/products/matlab.html). In Matlab, reaches of rivers with uniform gradient were visually selected on the longitudinal profile, exported to ArcGIS, and plotted on the digital elevation model. Deviations from steady-state longitudinal profiles can be produced by changes in rock-uplift rate, changes in substrate erodibility, and climatic gradients, so the results of the analysis must be evaluated for each of these forcing mechanisms.

Strong range-normal climatic gradients exist across the Himalaya, where the Lesser Himalaya receives two to four times more annual rainfall than the High Himalaya. Significant along-strike climatic gradients also exist across the arc, where the eastern Himalaya is significantly wetter than the western Himalaya. Despite these differences, the High Himalaya can be generalized as being contained within a single relatively dry climate and the Lesser Himalaya as belonging to a wetter climate. Because the emphasis of this study is strain accumulation in the High Himalaya, we have not weighted the flow-accumulation grid for climatic variations. Doing so would lower ksn values slightly in drier regions and raise them slightly in wetter regions, but would not change the overall pattern. Determining appropriate weighting factors for the various arc-normal and arc-parallel climatic variations present across the Himalaya is outside the scope of the current study, and given the low level of correlation between precipitation rates and river channel steepness, it would be unlikely to significantly affect the results of the analysis.

The relationship between rock uplift, upstream drainage area, and river channel morphology in glaciated regions cannot be effectively modeled using the ksn index. Measurements of modern Himalayan glacial equilibrium-line altitudes are broadly bracketed at 4000–5000 m (Owen and Benn, 2005). The Himalaya has a long complicated history of glaciations during which equilibrium-line altitudes have fluctuated considerably. During the last glacial maximum, many Himalayan glaciers reached below 1000 m elevation, and the toe of some nearly reached the Gangetic plain (Owen and Benn, 2005, and references therein). For this reason we, like in all other studies of Himalayan river channel steepness, cannot exclude all landscapes that have had any glacial influence. Instead we tried to exclude active glacial landscapes, while not eliminating more of the High Himalaya than necessary, by clipping all raster data sets excluding elevations above 4500 m. An unavoidable consequence of including areas that experienced past glaciation is the apparent oversteepening of the lowermost portion of tributaries where they join a formerly glaciated trunk channel. Widening of the trunk channel during glaciation produced a knickpoint in the tributary channel at their junction. This is seen in our data as a short (<1 km) anomalously steep segment at many high-elevation trunk-tributary junctions. However, the effect of these glacial artifacts on the orogen-scale map is small due to their limited spatial extent.

To test the dependency of ksn on interseismic versus coseismic uplift and how it relates to the generation of relief, we calculated correlation coefficients for ksn versus MHT coupling, earthquake epicenter density, and topography. We evaluated the contributions of lithology (proxy for substrate erodibility) and precipitation by generating correlation coefficients for ksn versus our geologic map [compiled from DeCelles et al. (1998); Long et al. (2011); Mitra et al. (2010); Yin (2006); Webb et al. (2007); and He et al. (2015)] and a 50 yr interpolation of precipitation data (Hijmans et al., 2005) (Supplemental Figs. S1–S31). The statistical analysis was conducted using ArcGIS and Microsoft Excel software. The ksn data constitute a vector data set; comparing it quantitatively with raster data sets like MHT coupling, topography, and precipitation required converting the individual vector shapefiles to a single raster file. We used the natural neighbors interpolation algorithm with a 1 km grid for the conversion. The final step in preparing the ksn data for display was to contour the data to highlight continuous clusters of high-ksn rivers. We define strain accumulation segments as regions near the BPT whose ksn values are consistently >200 (Fig. 2).

Similarly, earthquake epicenters of Mw >2.5 from A.D. 1966 to 2016 from the U.S. Geological Survey National Earthquake Information Center (NEIC) database (https://earthquake.usgs.gov/earthquakes/search/) were converted to a raster grid using the point density algorithm with an output cell size of 0.05° and circular search radius of 0.5° (Supplemental Figs. S1–S3 [footnote 1]). Data were then extracted from each data set along the same line using the “interpolate line” function on the “3D Analyst” toolbar in ArcGIS, followed by the “profile graph” function.

The resulting data were then exported to Excel, and the magnitudes of the data sets were compared statistically using Pearson correlation coefficients. Correlation coefficients statistically compare how data sets change in relation to one another, i.e., whether they change together (positive correlation) or opposite to each other (negative correlation), or display no systematic relation to one another (no correlation). Coefficients were generated using the data analyst tool in Excel, wherein magnitudes of all data sets being analyzed were arranged in adjacent columns, selected for analysis, and then compared using the formula:
graphic
in which x and y are the sets of values being compared and r is pearson correlation coefficient. After analysis, the magnitudes were normalized to 1 by dividing all values by the greatest value in each individual data set so they could all be displayed at the same scale (Fig. 4).
Figure 4.

(A–D) Plots of swaths showing the spatial relationships between river channel steepness (ksn), earthquake density, Main Himalayan thrust (MHT; intracontinental subduction zone) coupling, and topography. (E–H) Plots of swaths showing ksn and 50 yr interpolation of precipitation data (Hijmans et al., 2005), and bars showing distribution of lithotectonic units. Values normalized to one by dividing every entry in each dataset by the highest number in that dataset.

Figure 4.

(A–D) Plots of swaths showing the spatial relationships between river channel steepness (ksn), earthquake density, Main Himalayan thrust (MHT; intracontinental subduction zone) coupling, and topography. (E–H) Plots of swaths showing ksn and 50 yr interpolation of precipitation data (Hijmans et al., 2005), and bars showing distribution of lithotectonic units. Values normalized to one by dividing every entry in each dataset by the highest number in that dataset.

Comparisons of ksn with MHT coupling, earthquake density, topography, and precipitation are straightforward using existing raster data sets and the methods outlined above. However, there is no existing data set for substrate erodibility, requiring the use of a proxy. We used the material constant mi of the Hoek-Brown failure criterion (Hoek and Brown, 1980, 1997) as a basis for assigning each Himalayan lithotectonic unit an erodibility index. The Hoek-Brown failure criterion is an engineering estimate of rock mass strength; its material constant mi can be expressed as:
graphic
where σt is tensional stress and σc is compressive stress.

Hoek and Brown created their failure criterion as a way to estimate the strength of rock masses based on lithology through testing of uniaxial and tensile strength of different lithologies. We took the mi values of all lithologies within each Himalayan lithotectonic unit and averaged them to reflect average rock mass strength for each (Table 1), and then normalized the values to 1 for comparison with the other data sets. The resulting erodibility index values are: 0.5 for the Siwaliks, 0.7 for the Lesser Himalaya sequence, 1.0 for the Greater Himalaya sequence, and 0.5 for the Tethyan sedimentary sequence, which essentially separates the Himalaya into soft, intermediate, and hard rocks. To characterize the relationship between ksn, MHT coupling, substrate erodibility, and precipitation, we analyzed four orogen-normal profiles spaced at 25 km intervals for ten regions spanning the central Himalaya from Uttarakhand to Bhutan (Supplemental Figs. S1–S3 [footnote 1]). The results reported for the entire range (Table 2) are an average of all ten regions.

TABLE 1.

HOEK-BROWN ROCK MASS STRENGTH VALUES FOR HIMALAYAN LITHOTECTONIC UNITS

TABLE 2.

AVERAGE CORRELATION COEFFICIENTS FOR CORRELATION WITH NORMALIZED RIVER CHANNEL STEEPNESS (ksn), BY HIMALAYAN REGION

The results of the ksn analysis and the statistical analysis of potential forcing mechanisms and their effects are presented in the following sections. First, the role of precipitation and lithology in determining river channel steepness is evaluated. Next, the contribution of strain accumulation (interseismic and coseismic), in the form of MHT coupling and earthquake epicenter density, is examined. Finally, the relationship between topography, ksn, earthquakes, and drainage networks is investigated.

Climate

A map of average rainfall interpolated for the past 50 yr (Hijmans et al., 2005) was used to generate correlation coefficients for ksn and precipitation. The average r (Pearson correlation coefficient) for the entire Himalaya is –0.3 with a p < 0.005 (p-values measure the likelihood that the correlation results from the null hypothesis; any p < 0.05 is considered statistically robust) (Fig. 4; Supplemental Figs. S1–S3 [footnote 1]). No relationship was found between the average rate of rainfall and the strength of the correlation. The low degree of average correlation reveals that either the spatial pattern of Himalayan precipitation has undergone significant changes over the last ten thousand to hundred thousand years, or that erosion rates are decoupled from precipitation rates as suggested by recent studies of the Alpine-Himalayan orogen (Burbank et al., 2003; Godard et al., 2014; Adams et al., 2015, 2016; Forte et al., 2016).

Lithotectonic Unit

This analysis reveals a moderate degree of positive correlation (average r = 0.5, p < 0.05) for ksn and lithology (Fig. 4; Supplemental Figs. S1–S3 [footnote 1]). Regionally, the correlation is strong (0.6–0.8) across east Nepal and Bhutan (see Table 2), demonstrating the effect that the juxtaposition of lithologies with different erodibilities has on river channel steepness. The moderate degree of average correlation between ksn and litholotectonic unit is evidence of the role substrate erodibility plays in modulating the ksn signal; while it can be locally significant, e.g., capstone plateaus protecting more erodible substrates, regionally it appears secondary to rock-uplift rate in determining the pattern of river channel steepness.

MHT Coupling

The MHT coupling model of Stevens and Avouac (2015) is based primarily on GPS data but also uses the northern limit of clusters of earthquake epicenters to aid in locating where the plate boundary becomes uncoupled. MHT coupling displays a moderately strong correlation with ksn (r = –0.6, p < 0.05) (Fig. 4; Supplemental Figs. S1–S3 [footnote 1]). The correlation between the degree of MHT coupling and river channel steepness suggests that in the Himalaya, rock-uplift rate is the primary driver of river channel steepening, and relates it directly to interseismic strain. Tectonic steepening of rivers takes place over hundreds of thousands of years and is the net result of the interaction between interseismic, coseismic, and post-seismic strain accumulation and erosion; the relation between ksn and coseismic strain release will be examined in the next section. Regionally, this correlation is strong across the western Himalaya (average r = –0.7) and becomes weak in the east (average r = –0.3) (see Table 2). Whether this is related to the paucity of GPS stations in easternmost Nepal, Sikkim (northeastern India), Bhutan, and Arunachal Pradesh (northeasternmost India) is unknown. While the GPS data upon which the MHT coupling model is based were recorded over tens of years, it takes far longer for the ksn signal to be recorded in the landscape. While the exact amount of time required to impart a ksn signal to the landscape is unknown, it is certainly over numerous seismic cycles. The strength of correlation between decadal interseismic elastic strain and ksn demonstrates that strain accumulation has been localized in its current location for a period of time on the order of 100,000 yr. Our analysis suggests that of the three forcing mechanisms thought to influence river channel steepness (substrate erodibility, erosion, and rock-uplift rate), in the Himalaya, rock-uplift rate is the primary driver, with secondary influence from lithology, and to a more minor extent by the pattern of modern precipitation. Our correlations agree well with recent studies highlighting the dominance of tectonic processes over precipitation gradients in driving rock uplift (Burbank et al., 2003; Godard et al., 2014; Adams et al., 2015, 2016; Forte et al., 2016).

Earthquake Epicenter Density

Earthquakes occur when accumulated interseismic strain overcomes static friction on the strongly coupled portion of a fault. Stored interseismic (elastic) strain is released, resulting in coseismic displacement of material and inelastic deformation taking the form of crustal thickening, thinning, or lateral offset. The spatial density of earthquake epicenters is used here as a representation of the amount and location of coseismic energy released. We use epicenters from the NEIC catalog from 1 January 1966 to 1 January 2016 because it is the longest continuous record available for the entire Himalayan arc. The 50 yr record of earthquakes used in our analysis is impacted by the presence of the 2015 Gorkha earthquake and its sequence of aftershocks, which dominate the east Nepal segment. Similarly, moderate-magnitude earthquakes and their aftershock sequences create local epicenter “hotspots” across the arc. Along-strike variations in the distribution of modern earthquakes should be viewed as a snapshot of coseismic energy release; we expect that given a sufficiently long record, earthquake clusters near the BPT will span the entire arc. We acknowledge that the 50 yr record used in this analysis is only a meaningful metric of coseismic energy release in regions where there have been modern earthquake sequences, and do not mean to imply that the lack of modern earthquakes in intervening regions suggests a fundamental difference in long-term strain accumulation. Instead we rely on the 105 yr record of strain accumulation provided by ksn mapping to denote regional differences. River channel steepness is impacted by both inter- and coseismic strain and as such should correlate well with any record of strain accumulation and release. This led us to test the relationship of ksn with earthquake epicenter density. In regions with a sufficient number of earthquakes to produce a statistically significant result, epicenter density has a moderately strong correlation with ksn (r = 0.6, p < 0.05) and relates river channel steepness directly to coseismic strain release and concomitant inelastic crustal deformation (Fig. 4; and Supplemental Figs. S1–S3 [footnote 1]). Eight out of the ten regions analyzed met this criterion; the remaining two (Arun river area) and west Bhutan) were not used in calculating the average correlation. The moderately high degree of correlation between ksn and earthquake epicenter density is demonstrative of the link between the steepening of river channels and crustal thickening by internal deformation of the high Himalaya. In regions with sufficient earthquakes for a robust statistical analysis, epicenter density peaks spatially southward of (foreland to) ksn, supporting the idea that the High Himalaya acts as an elastic strain reservoir driving earthquakes on the strongly coupled portion of the plate boundary, and also demonstrates that rock uplift during earthquakes contributes to the ksn signal.

Topography

Topography and ksn display a moderately strong correlation (r = 0.6, p < 0.005), with the strongest correlation in the central and east Nepal segments (average r = 0.7) (Fig. 4; Supplemental Figs. S1–S3 [footnote 1]). This is expected, as the same forces primarily responsible for steepening stream channels (stored elastic strain and crustal thickening) also generate the high topography. The segments with the strongest correlation between ksn and topography are the same segments that host all of the Himalayan peaks >8000 m (excluding the western syntaxis), suggesting that these segments also accumulate strain more rapidly than the rest of the Himalayan range.

Normalized River Channel Steepness

Over its 2500 km length, regions of rapid rock uplift in the High Himalaya are characterized by spatial clusters of high-ksn rivers that are offset from one another across segment boundaries that coincide with linked strike-slip Himalayan-Tibetan rift systems. Nearly all high-ksn clusters overlie weakly coupled portions of the MHT located 60–120 km north of the range front (Figs. 2 and 5). Only west Nepal and Bhutan host high-ksn clusters that overlie strongly coupled portions of the MHT, 30–50 km from the range front (Figs. 2 and 5). The spatial dimensions of the high-ksn clusters vary in width, length, and distance from the range front and are summarized in Table 3.

Figure 5.

Spatial distribution and dimensions of strain accumulation segments (polygons) interpreted from clusters of high normalized river channel steepness (ksn) overlaid on convergence obliquity boundaries and the position of cross-orogen strain partitioning structures. This highlights the size and more foreland position of the east Nepal strain accumulation segment with respect to its neighbors, and its location in the orthogonally convergent sector of the range. This is juxtaposed against strain accumulation in the neighboring west Nepal and Bhutan segments, which are smaller and offset toward the hinterland in a sense consistent with the sense of motion identified for the West Nepal fault system (WNFS) and predicted by the rotation of magnetic remanences across the Yadong-Gulu rift (YGR). MFT—Main Frontal thrust; GM—Gurla Mandhata core complex.

Figure 5.

Spatial distribution and dimensions of strain accumulation segments (polygons) interpreted from clusters of high normalized river channel steepness (ksn) overlaid on convergence obliquity boundaries and the position of cross-orogen strain partitioning structures. This highlights the size and more foreland position of the east Nepal strain accumulation segment with respect to its neighbors, and its location in the orthogonally convergent sector of the range. This is juxtaposed against strain accumulation in the neighboring west Nepal and Bhutan segments, which are smaller and offset toward the hinterland in a sense consistent with the sense of motion identified for the West Nepal fault system (WNFS) and predicted by the rotation of magnetic remanences across the Yadong-Gulu rift (YGR). MFT—Main Frontal thrust; GM—Gurla Mandhata core complex.

TABLE 3.

SPATIAL DIMENSIONS AND DISTANCE FROM THE HIMALAYAN RANGE FRONT OF HIGH-ksn CLUSTERS

High-ksn clusters overlying weakly coupled portions of the MHT are located at variable distances (60–120 km) from the range front, but nearly all are between 70 km and 100 km. Neighboring segments host clusters that are offset by 40–60 km. This indicates that whatever mechanism(s) (passive uplift above MHT ramps, out-of-sequence thrusting, duplexing, tectonic wedging) are responsible for rapid rock uplift within the clusters, it is not purely a function of distance from the range front as would be expected if the geometry of the MHT was uniform along strike. The segment boundaries are collocated with the Himalayan-Tibetan rift systems and strain partitioning structures, suggesting that they may share a common mechanism for their development (Figs. 3 and 5). Two segment boundaries, west Nepal–east Nepal and east Nepal–Bhutan, are located near the northward projection of Indian basement ridges (Gahalaut and Kundu, 2012) (Fig. 1), which may play a role in segmenting the range. The largest cluster of high-ksn river channels is located in the east Nepal segment, which goes from the Thakkhola graben in the west to the Yadong-Gulu rift in the east, just over 500 km (Fig. 2). The high-ksn cluster that defines the east Nepal segment is closer to the toe of the thrust wedge (60–70 km) than any other cluster overlying a weakly coupled portion of the MHT (Fig. 2). Only west Nepal and Bhutan contain high-ksn clusters closer to the range front (30–50 km), but these overlie strongly coupled portions MHT (Fig. 2). In both west Nepal and Bhutan, a second high-ksn cluster parallels the first but overlies the weakly coupled portion of the MHT 80–120 km from the toe of the thrust wedge. In west Nepal, this has been interpreted to reflect an active MHT flat connecting a hinterland ramp with the more foreland ramp (Harvey et al., 2015); our analysis suggests that a similar geometry may be present in Bhutan.

The segments documented here are typically 200–300 km in orogen-parallel length, with the exception of east Nepal which is 500 km long (Figs. 2 and 5; and Table 3). This regular spacing suggests that segment boundaries and Himalayan-Tibetan rifts may be the product of a frequency-dependent mechanism, possibly buckling of the downgoing Indian plate. Buckling of slabs is a function of arc curvature and slab thickness, which, in a homogenous medium, would result in a set wavelength and amplitude (Yin, 2000). Another possibility is that both rifts and ksn segment boundaries are the geometric consequence of the curvature of the Himalayan arc (McCaffrey and Nabelek, 1998).

Structural analysis of the deformed Siwalik Group at the leading edge of the MHT in western Nepal indicates the presence of north-striking lateral ramps in the Main Frontal thrust (Mugnier et al., 1999) that may extend far into the thrust wedge as lateral ramps in the MHT. It is also observed that the intersection of Indian basement ridges and the Main Frontal thrust coincide spatially with the Uttarakhand–west Nepal and west Nepal–east Nepal segment boundaries, and project near the east Nepal–Bhutan boundary (Fig. 1) (Gahalaut and Kundu, 2012). The timing of plastic deformation and metamorphism varies across the Faizabad ridge (west Nepal–east Nepal boundary), suggesting that these differences could be explained by lateral ramps in the MHT produced by subduction of these features (Gibson et al., 2016). Based on our results, these hypothesized lateral ramps may also explain the margins of the strain accumulation segments. Another possibility is that lateral ramps are generated by slip on syncollisional normal faults striking at a high angle to the thrust wedge, such as the Thakkhola graben, Kaurik-Chango (Leo Pargil), Gurla Mandhata–Humla, and Yadong-Gulu rift systems (Fig. 1A). If true, slip on these extensional structures is coeval with thrusting below. Motion on these upper crustal structures and thrusting at lower crustal levels along lateral ramps along the MHT could produce orogen-normal shear zones that are kinematically linked with upper and lower structural levels, as has been previously suggested by Harvey et al. (2015).

Strain Partitioning Boundaries and Himalayan-Tibetan Rifts

Oblique convergence at subduction zones partitions strain into arc-normal and arc-parallel components (McCaffrey, 1992; Tikoff and Teyssier, 1994; Braun and Beaumont, 1995; Ellis et al., 1995; McCaffrey et al., 2000), and is present in the Himalaya as well. The most striking Himalayan example exists in the western Himalaya–Tibetan Plateau where the right-lateral Karakorum fault feeds slip through a major extensional stepover at the Gurla Mandhata core complex to the right-lateral Western Nepal fault system and then to the Main Boundary thrust near the toe of the thrust wedge (Armijo et al., 1986; McCaffrey and Nabalek, 1998; Murphy and Copeland, 2005; Styron et al., 2011; Murphy et al., 2014; Whipp et al., 2014; Silver et al., 2015). This linked system of faults divides the Himalaya into a western obliquely convergent sliver and a central orthogonally convergent region. At the eastern edge of the orthogonally convergent region is the Yadong-Gulu rift, which extends from central Tibet to the Lesser Himalaya and has been suggested to be a strain partitioning boundary (Wu et al., 1998; Velasco et al., 2007; Drukpa et al., 2006; Antolín et al., 2012), although the case for it is less clear than for the Karakorum fault–Western Nepal fault system. The Karakorum fault–Western Nepal fault system and Yadong-Gulu rift separate orthogonally convergent eastern Nepal from its obliquely convergent neighbors. Upon crossing either of these structures, the change is not simply kinematic; the spatiotemporal pattern of seismicity also changes. East Nepal features a continuous band of seismicity interpreted to mark a crustal scale ramp in the MHT (Pandey et al., 1995, 1999; Cattin and Avouac, 2000). Upon crossing from east to west the Western Nepal fault system, this band of seismicity splits, becoming two distinct bands (Cattin and Avouac, 2000; Harvey et al., 2015); similarly, upon crossing from west to east the Yadong-Gulu rift, the band of seismicity disappears, becoming a diffuse scattering of earthquakes across the Kingdom of Bhutan (Velasco et al., 2007; Drukpa et al., 2006) (Fig. 2). The strike-slip nature of the Western Nepal fault system separating west and east Nepal has been documented (Murphy et al., 2014; Silver et al., 2015); however, fault-slip data implicating strike-slip motion along the Yadong-Gulu rift is lacking. Nevertheless, the High Himalaya of Bhutan is offset from that of east Nepal and Sikkim by 70 km in a left-lateral sense (Wu et al., 1998); studies of seismicity find distributed mid-crustal strike-slip deformation across southern Bhutan (Velasco et al., 2007; Drukpa et al., 2006). Magnetic remanence in western Bhutan requires clockwise vertical-axis rotation, consistent with left-lateral slip on the Lingshi fault, an oblique-slip normal fault that runs parallel to the southern Yadong-Gulu rift on the Bhutanese side of the border (Antolín et al., 2012). The pattern of seismicity in Nepal is reflected in river channel steepness as a 500-km-long cluster of high-ksn rivers located immediately north of the band of seismicity. In west Nepal there are two clusters of high-ksn rivers, each associated with a band of seismicity. However, in Bhutan the pattern of river channel steepness is distinct from the diffuse distribution of earthquakes; high-ksn rivers are clustered in two narrow bands, one near the Tibetan Plateau margin and one near the southern border with India. We interpret these two bands to reflect the presence of thrust ramp duplexes or antiformal stacks generated by tectonic wedging, upon which strain accumulates at a relatively slower rate than in east Nepal.

Accumulation and Release of Plate Boundary Strain

The High Himalaya overlies weakly coupled portions of the MHT (Stevens and Avouac, 2015) and separates a stably sliding décollement beneath southern Tibet from strongly coupled portions of the MHT, suggesting that the High Himalaya acts as a dashpot storing interseismic strain (Bilham et al., 1997, 2001; Jouanne et al., 2004; Feldl and Bilham, 2006; Meade, 2010; Ader et al., 2012). In this model, southernmost Tibet and the High Himalaya bend elastically to accommodate India-Asia convergence until the elastically stored energy is sufficient to overcome the static friction on the strongly coupled portion of the MHT. This buildup and release was highlighted by the 25 April 2015 Gorkha earthquake, which resulted in ∼0.5 m of surface lowering of the leading edge of the High Himalaya and 1 m of uplift of the Lesser Himalaya (Lindsey et al., 2015; Elliott et al., 2016). Despite its large magnitude (Mw 7.9), this earthquake did not rupture the strongly coupled MHT. Rather, it propagated along the boundary between weakly and strongly coupled areas, uplifting the intervening moderately coupled region and transferring strain equivalent to a Mw 6.9 earthquake from the High Himalaya to the frontal portion of the range (Grandin et al., 2015).

The spatial distribution of historical earthquake ruptures in the Himalaya is an area of ongoing research (Ambraseys and Jackson, 2003; Kumar et al., 2006, 2010; Mugnier et al., 2013; Murphy et al., 2014; Berthet et al., 2014; Bollinger et al., 2014; Rajendran et al., 2015; Hossler et al., 2016). Our analysis reveals a rough correspondence between ruptures and segment boundaries delineated by high-ksn clusters. While major earthquake ruptures are each largely confined to a single segment, the edges of rupture patches in some cases bleed over into adjacent segments. This indicates that segmented strain accumulation controls where and when major earthquakes initiate, but that segment boundaries are not necessarily barriers to rupture propagation. Rupture patches of large Himalayan earthquakes since A.D. 1900 correspond to segments defined by high-ksn clusters. In our study area, there have been seven major Himalayan earthquakes during the past 500 yr, of which six are confined to individual segments. Only the A.D. 1505 Mw 8.2 earthquake is an exception; it ruptured the entire west Nepal segment, the eastern edge of the Uttarakhand segment, and the western edge of the east Nepal segment.

The most striking correspondence between the ksn analysis and the spatiotemporal distribution of historical earthquakes lies within the Nepal-Bhutan segments. The west Nepal segment has not ruptured since 1505 (Ambraseys and Jackson, 2003) and Bhutan has not ruptured since 1713 (Berthet et al., 2014), while the east Nepal segment has ruptured three times over the last three hundred years in 1833, 1934, and 2015 (Fig. 6). This apparent dichotomy may simply reflect a lack of resolution in our knowledge of the Himalayan seismic cycle, or it may indicate that the central Nepal segment accumulates and releases strain more efficiently than its neighbors. This latter idea is supported by geomorphologic and kinematic considerations.

Figure 6.

Contour map of clusters of high normalized river channel steepness (ksn) overlaid with inferred earthquake ruptures of the past 500 yr, labeled with year of occurrence, after Berthet et al. (2014, and references therein) and Mugnier et al. (2013, and references therein). Heavy white lines denote strain accumulation segments referred to in text. Thin white lines are political boundaries and purple lines are active faults from Taylor and Yin (2009).

Figure 6.

Contour map of clusters of high normalized river channel steepness (ksn) overlaid with inferred earthquake ruptures of the past 500 yr, labeled with year of occurrence, after Berthet et al. (2014, and references therein) and Mugnier et al. (2013, and references therein). Heavy white lines denote strain accumulation segments referred to in text. Thin white lines are political boundaries and purple lines are active faults from Taylor and Yin (2009).

The east Nepal strain accumulation segment occupies the most orthogonally convergent sector of the orogen (Figs. 3 and 5), consistent with pure shear shortening, and hosts a cluster of high-ksn channels that is broader, longer, and closer to the range front than those anywhere else in the Himalaya. The east Nepal segment contains eight out of the ten Himalayan peaks higher than 8000 m, and has well developed frontal anticlines that collect all of the High Himalayan catchments into three outlets along the range front. Taking the highest topography in the range as a measure of accumulated High Himalayan strain, and the most well-developed frontal anticlines in the range as a measure of Sub-Himalayan strain accumulation, suggests that range-normal strain accumulation and release is more efficient here than anywhere else in the Himalaya.

West Nepal and Bhutan lie within the obliquely convergent western and eastern sectors of the Himalaya, respectively, in which some portion of India-Asia convergence is partitioned into orogen-parallel motion (Drukpa et al., 2006; Velasco et al., 2007; Styron et al., 2011; Antolín et al., 2012; Murphy et al., 2014; Whipp et al., 2014). Changes in the size of the hinterland region accumulating strain, its proximity to the toe of the thrust wedge, and lateral strain partitioning provide a plausible explanation for the seemingly disparate seismic cycles of west Nepal, east Nepal, and Bhutan.

Looking at earthquake ruptures even further back, we encounter an apparent change in seismic behavior. There have been four surface-rupturing events identified between A.D. 1000 and A.D. 1500: 1100, 1255, 1344, and 1430, three of which have inferred rupture lengths between 500 and 1000 km requiring magnitudes of Mw 8–9 (Sapkota et al., 2013; Berthet et al., 2014; Mugnier et al., 2013). While this may be the result of conflating two (or more) smaller ruptures spaced closely in time, it may be evidence of a >1000 yr seismic cycle that culminates in an ∼Mw 9.0 earthquake (Kumar et al., 2006; Stevens and Avouac, 2016). Viewing the recent sequence of large earthquakes in central Nepal (in 1833, 1934, and 2015) in this context, only the 1934 surface-rupturing event represents a significant decrease in accumulated strain. The 1833 and 2015 earthquakes instead represent strain transfer from the hinterland to the foreland. This idea is supported by an apparent strain deficit between the magnitude of convergence during the past 500 yr and the magnitude of strain released in great earthquakes, which is estimated to be equivalent to four Mw 8 earthquakes or a single Mw 9 earthquake (Bilham and Ambraseys, 2005; Meade, 2010; Stevens and Avouac, 2016).

Drainage Network

The drainage network of major Himalayan rivers is segmented similarly to channel steepness (Fig. 1), suggesting that channel steepness and drainage network patterns form in response to the same tectonic forcings. The Nepal segments are characterized by large catchments in which numerous parallel drainages are collected into a few pour points at the range front. This process is driven by the location of active anticlines in the Sub-Himalaya, which grow more rapidly than the individual rivers can incise through them, producing lower river reaches that flow parallel to the range front for 70–100 km. The result is a 760-km-long segment of the range hosting only three pour points entering the foreland basin—the Karnali, Narayani, and Arun rivers. In west Nepal, near the Tibetan Plateau boundary, in the upper reaches of the Karnali river, the Mugu Karnali in the west and Dolpo Karnali in the east flow along the axis of the Dolpo syncline, paralleling the range front for 100 km and 150 km, respectively. Axial drainage along the leading edge of the growing Dolpo anticline (Cannon and Murphy, 2014) by the Karnali’s two main trunk channels is funneled into a zig-zag drainage pattern in one of the largest erosional bites in the Himalaya (Fig. 1). The lower reaches of the Karnali contain two hairpin turns, suggesting recent drainage network reorganization by the development of active anticlines at the range front. The 270-km-long Bhutan segment contains five catchment basins, the largest of which hosts a parallel drainage network characteristic of steep, active mountain fronts. Tectonic activity along the range front in Bhutan is expressed differently than in the western and central Himalaya. To the west, the range front in Nepal is characterized by closely spaced anticlines separated from the range to the north by broad flat valleys (duns). This is in stark contrast to Bhutan, where foothill anticlines are intermittent, duns are narrow, the Main Frontal thrust is difficult to recognize, and the first topographic step is commonly associated with the Main Boundary thrust. However, these differences are minor compared to the presence of the Shillong Plateau in the eastern Himalaya foreland basin. The Shillong Plateau is located 100 km south of the range front and has >1 km of relief. It has been characterized as a popup structure bounded to the north and south by oppositely dipping thrusts, however its significance for the thrust wedge and its implications for the interaction between climate and tectonics are debated (Adams et al., 2015; Bilham and England, 2001; Clark and Bilham, 2008; Grujic et al., 2006; Najman et al., 2016). Whatever role the Shillong Plateau plays in accommodating India-Asia convergence, it seems clear that at the range front strain is accumulating at a lower rate than in Nepal, as all drainages flow across the trace of the Main Frontal thrust without disruption.

Despite numerous geologic and geophysical studies of the Himalayan plate boundary, basic questions remain unanswered. Foremost is, what is the geometry of the subduction interface? Seismic reflection surveys have resolved ramps and flats in the MHT (Hauck et al., 1998; Caldwell et al., 2013), although seismic refraction surveys have failed to do so (Nábělek et al., 2009). Geodesists are able to model the pattern of MHT fault coupling using a single inclined plane by assuming that all significant variability in the position and dip of MHT ramps is located within strongly coupled portions (Bollinger et al., 2004b; Jouanne et al., 2004; Stevens and Avouac, 2015). Nevertheless, there is evidence for the existence of MHT ramps located within the weakly coupled portions of the MHT. Balanced cross sections (Schelling and Arita, 1991; DeCelles et al., 1998; Robinson, 2008; McQuarrie and Ehlers, 2015; Hubbard et al., 2016), electrical conductivity (Lemonnier et al., 1999), tectonic geomorphology (Seeber and Gornitz, 1983; Cannon and Murphy, 2014; Harvey et al., 2015; Morell et al., 2015), microseismicity (Fig. 2) (Pandey et al., 1999), and Bayesian modeling of possible MHT geometries (Elliott et al., 2016) all suggest the presence of thrust ramps in essentially the same locations, at the boundary between the strongly and weakly coupled portions of the MHT. Changes in the position of high-ksn clusters are most pronounced between west Nepal and east Nepal, where high-ksn clusters overlying weakly coupled portions of the MHT are offset by 60 km along a sharp boundary coincident with the Gurla Mandhata–Humla core complex (Harvey et al., 2015) (Figs. 2 and 5). This indicates that the geometry of the MHT changes along strike. That the largest offset of high-ksn clusters occurs here has kinematic significance. The Gurla Mandhata–Humla core complex is thought to be an extensional stepover feeding slip from the right-lateral Karakorum fault in southern Tibet into the right-lateral cross-orogen West Nepal fault system (Murphy et al., 2014; Silver et al., 2015).

East Nepal is located in the most orthogonally convergent part of the Himalayan arc (Fig. 5) and during that past 500 yr has had the most rapid recurrence of Mw >7 earthquakes anywhere in the range (Fig. 6). Convergence obliquity and great earthquake recurrence times increase to the west and east. How much convergence obliquity modifies the geometry of the plate margin or modulates its seismic behavior remains unknown, but is increasingly seen as an important seismotectonic variable (McCaffrey and Nabalek, 1998; Jouanne et al., 1999; Banerjee and Bürgmann, 2002; Styron et al., 2011; Murphy et al., 2014; Whipp et al., 2014; Kundu et al., 2014).

Middle Miocene to recent internal deformation of the Himalayan thrust wedge appears to be dominated by thrust duplexing (Bollinger et al., 2004a; Herman et al., 2010; Webb, 2013; He et al., 2015; Larson et al., 2015; Yu et al., 2015). At the surface, duplexing is characterized by a high ratio of crustal thickening to surficial horizontal shortening. This criterion is met in the hinterland of west Nepal where the broad Dolpo anticline records vertical thickening to surficial horizontal shortening ratios of between 2:1 and 5:1, depending on the thickness of the repeated layer (Cannon and Murphy, 2014). Duplexes typically form by underplating, whereby slices of the footwall are accreted to the hanging wall (Fig. 7), however no structural window in the High Himalaya is documented that provides an unambiguous view of the mechanism thickening the hinterland. Alternatively, in the plastically deforming lower crust, repeated tectonic wedging (Webb et al., 2007; Yin, 2006) at the BPT could produce a duplex restricted to the hanging wall of the MHT (Fig. 7). Tectonic wedging focused at the BPT could alleviate some of the slip deficit by allowing the lower crust to accommodate more convergence than the upper crust. In this view, hinterland thickening results from repeated insertion of plastic wedges at a ramp where the MHT BPT is localized, forming an antiformal stack. As the active wedge is inserted at the base of the stack, it forces the preceding wedges up, building topography and steepening river channels at the surface. The brittle upper crust overlying the antiformal stack absorbs elastic strain imparted by the south-southwest–directed sliding of southern Tibet recorded by GPS (Ader et al., 2012). Hinterland lower to mid-crustal duplexing provides a plausible explanation for the coincidence of interseismic strain accumulation in the thrust wedge, spatial clusters of high ksn values, and the highest topography in the range without calling upon a delayed erosional response to ramp migration (Grandin et al., 2012), and diminishes the need for an unbalanced earthquake cycle (Bilham and Ambraseys, 2005; Meade, 2010; Stevens and Avouac, 2016) by allowing the lower crust to accommodate more permanent strain than the upper crust. Recently, several mid-crustal discontinuities have been identified within the Greater Himalaya sequence of Nepal predating the main period of the linked Main Central thrust–South Tibetan detachment motion in the late early to middle Miocene (Carosi et al., 2010; Larson et al., 2015; Larson and Cottle, 2014; Montomoli et al., 2013). These discontinuities have been interpreted as shear zones bounding horses within a hinterland duplex that was exhumed during thrusting on the Main Central thrust. Hinterland duplexes create permanent uplift by thickening the lower crust while storing elastic strain in the upper crust, the net result of which is collocated increased rates of rock uplift and elastic strain accumulation. Currently, we cannot distinguish between these two duplexing models, although the presence of isotopic signatures in anatectic melts of the Gurla Mandhata core complex provides compelling evidence for migration of Indian-affinity material across the MHT (Murphy, 2007).

Figure 7.

Crustal thickening models for the Himalaya highlighting the differences between duplexing (A) and tectonic wedging (B).

Figure 7.

Crustal thickening models for the Himalaya highlighting the differences between duplexing (A) and tectonic wedging (B).

Rapid rates of inferred rock uplift in the Himalaya, recorded by ksn analysis, are predominantly localized within the High Himalaya above weakly coupled portions of the MHT. The time scales required to encode ksn signals in the landscape (105 yr) provide an intermediary time step between GPS data sets and low-temperature thermochronology. Clusters of high-ksn rivers are offset from each other in map view across Himalayan-Tibetan strain partitioning structures and rift systems. This suggests that the location of linked strike-slip and extensional structures is fundamental to the structural segmentation of the orogen. Regions of rapid rock uplift in the High Himalaya separate stably sliding southern Tibet from the Himalayan seismogenic zone, and acts as a dashpot accumulating and storing elastic strain. In addition to acting as a reservoir for accumulating elastic strain, collocation of clusters of high-ksn streams with the highest topography in the Himalaya indicates that they also record crustal thickening. As such, clusters of high-ksn streams represent an accumulation of combined interseismic and coseismic strain. Structural segmentation of the Himalaya based on spatial clustering of high ksn values reveals an association with the spatiotemporal distribution of historical great earthquakes in the central Himalaya. This relationship suggests that the density and spatial dimensions of regions characterized by continuous high values of ksn, their distance from the toe of the thrust wedge, and the degree of convergence obliquity impact the recurrence of major earthquakes in active orogenic systems.

The manuscript benefited from reviews by Dr. Jon Harvey and an anonymous reviewer. The authors would especially like to thank Dr. Harvey for thorough and insightful reviews which led to substantiative improvements to the manuscript.

1Supplemental Figures S1–S3. Locations of the profiles analyzed and the spatial relationships between ksn, MHT coupling, earthquake epicenter density, lithotectonic unit, and climate for each of the 10 regions analyzed in this study. Please visit http://doi.org/10.1130/GES01508.S1 or the full-text article on www.gsapubs.org to view the Supplemental Figures.
Science Editor: Raymond M. Russo
1.
Adams
,
B.A.
,
Hodges
,
K.V.
,
van Soest
,
M.C.
, and
Whipple
,
K.X.
,
2013
,
Evidence for Pliocene–Quaternary normal faulting in the hinterland of the Bhutan Himalaya
:
Lithosphere
 , v.
5
, p.
438
449
, https://doi.org/10.1130/L277.1.
2.
Adams
,
B.A.
,
Hodges
,
K.V.
,
Whipple
,
K.X.
,
Ehlers
,
T.A.
,
Soest
,
M.C.
, and
Wartho
,
J.
,
2015
,
Constraints on the tectonic and landscape evolution of the Bhutan Himalaya from thermochronometry
:
Tectonics
 , v.
34
, p.
1329
1347
, https://doi.org/10.1002/2015TC003853.
3.
Adams
,
B.A.
,
Whipple
,
K.X.
,
Hodges
,
KV.
.
, and
Heimsath
,
A.M.
,
2016
,
In situ development of high-elevation, low-relief landscapes via duplex deformation in the Eastern Himalayan hinterland, Bhutan: Journal of Geophysical Research
:
Earth Surface
 , v.
121
, p.
294
319
, https://doi.org/10.1002/2015JF003508.
4.
Ader
,
T.
,
Avouac
,
J.-P.
,
Liu-Zeng
,
J.
,
Lyon-Caen
,
H.
,
Bollinger
,
L.
,
Galetzka
,
J.
,
Genrich
,
J.
,
Thomas
,
M.
,
Chanard
,
K.
,
Sapkota
,
S.N.
,
Rajaure
,
S.
,
Shrestha
,
P.
,
Ding
,
L.
, and
Flouzat
,
M.
,
2012
,
Convergence rate across the Nepal Himalaya and interseismic coupling on the Main Himalayan Thrust: Implications for seismic hazard
:
Journal of Geophysical Research
 , v.
117
, B04403, https://doi.org/10.1029/2011JB009071.
5.
Allen
,
G.H.
,
Barnes
,
J.B.
,
Pavelsky
,
T.M.
, and
Kirby
,
E.
,
2013
,
Lithologic and tectonic controls on bedrock channel form at the northwest Himalayan front: Journal of Geophysical Research
:
Earth Surface
 , v.
118
, p.
1806
1825
, https://doi.org/10.1002/jgrf.20113.
6.
Ambraseys
,
N.
, and
Jackson
,
D.
,
2003
,
A note on early earthquakes in northern India and southern Tibet
:
Current Science
 , v.
84
, p.
570
582
.
7.
Antolín
,
B.
,
Schill
,
E.
,
Grujic
,
D.
,
Baule
,
S.
,
Quidelleur
,
X.
,
Appel
,
E.
, and
Waldhör
,
M.
,
2012
,
E–W extension and block rotation of the southeastern Tibet: Unravelling late deformation stages in the eastern Himalayas (NW Bhutan) by means of pyrrhotite remanences
:
Journal of Structural Geology
 , v.
42
, p.
19
33
, https://doi.org/10.1016/j.jsg.2012.07.003.
8.
Armijo
,
R.
,
Tapponnier
,
P.
,
Mercier
,
J.L.
, and
Han
,
T.L.
,
1986
,
Quaternary extension in southern Tibet: Field observations and tectonic implications
:
Journal of Geophysical Research
 , v.
91
, p.
13,803
13,872
, https://doi.org/10.1029/JB091iB14p13803.
9.
Avouac
,
J.-P.
,
2003
,
Mountain building, erosion, and the seismic cycle in the Nepal Himalaya
:
Advances in Geophysics
 , v.
46
, p.
1
80
, https://doi.org/10.1016/S0065-2687(03)46001-9.
10.
Banerjee
,
P.
, and
Bürgmann
,
R.
,
2002
,
Convergence across the northwest Himalaya from GPS measurements
:
Geophysical Research Letters
 , v.
29
, no. 13, https://doi.org/10.1029/2002GL015184.
11.
Berthet
,
T.
,
Ritz
,
J.-F.
,
Ferry
,
M.
,
Pelgay
,
P.
,
Cattin
,
R.
,
Drukpa
,
D.
,
Braucher
,
R.
, and
Hetényi
,
G.
,
2014
,
Active tectonics of the eastern Himalaya: New constraints from the first tectonic geomorphology study in southern Bhutan
:
Geology
 , v.
42
, p.
427
430
, https://doi.org/10.1130/G35162.1.
12.
Bilham
,
R.
, and
Ambraseys
,
N.
,
2005
,
Apparent Himalayan slip deficit from the summation of seismic moments for Himalayan earthquakes, 1500–2000
:
Current Science
 , v.
88
, p.
1658
1663
.
13.
Bilham
,
R.
, and
England
,
P.
,
2001
,
Plateau ‘pop-up’ in the great 1897 Assam earthquake
:
Nature
 , v.
410
, p.
806
809
, https://doi.org/10.1038/35071057.
14.
Bilham
,
R.
,
Larson
,
K.
, and
Freymueller
,
J.
,
1997
,
GPS measurements of present-day convergence across the Nepal Himalaya
:
Nature
 , v.
386
, p.
61
64
, https://doi.org/10.1038/386061a0.
15.
Bilham
,
R.
,
Gaur
,
V.K.
, and
Molnar
,
P.
,
2001
,
Himalayan seismic hazard
:
Science
 , v.
293
, p.
1442
1444
, https://doi.org/10.1126/science.1062584.
16.
Bollinger
,
L.
,
Avouac
,
J.P.
,
Beyssac
,
O.
,
Catlos
,
E.J.
,
Harrison
,
T.M.
,
Grove
,
M.
,
Goffé
,
B.
, and
Sapkota
,
S.
,
2004a
,
Thermal structure and exhumation history of the Lesser Himalaya in central Nepal
:
Tectonics
 , v.
23
, TC5015, https://doi.org/10.1029/2003TC001564.
17.
Bollinger
,
L.
,
Avouac
,
J.P.
,
Cattin
,
R.
, and
Pandey
,
M.R.
,
2004b
,
Stress buildup in the Himalaya
:
Journal of Geophysical Research
 , v.
109
, B11405, https://doi.org/10.1029/2003JB002911.
18.
Bollinger
,
L.
,
Perrier
,
F.
,
Avouac
,
J.-P.
,
Sapkota
,
S.
,
Gautam
,
U.
, and
Tiwari
,
D.R.
,
2007
,
Seasonal modulation of seismicity in the Himalaya of Nepal
:
Geophysical Research Letters
 , v.
34
, L08304, https://doi.org/10.1029/2006GL029192.
19.
Bollinger
,
L.
,
Sapkota
,
S.N.
,
Tapponnier
,
P.
,
Klinger
,
Y.
,
Rizza
,
M.
,
Van der Woerd
,
J.
,
Tiwari
,
D.R.
,
Pandey
,
R.
,
Bitri
,
A.
, and
Bes de Berc
,
S.
,
2014
,
Estimating the return times of great Himalayan earthquakes in eastern Nepal: Evidence from the Patu and Bardibas strands of the Main Frontal Thrust: Journal of Geophysical Research
:
Solid Earth
 , v.
119
, p.
7123
7163
, https://doi.org/10.1002/2014JB010970.
20.
Braun
,
J.
, and
Beaumont
,
C.
,
1995
,
Three-dimensional numerical experiments of strain partitioning at oblique plate boundaries: Implications for contrasting tectonic styles in the southern Coast Ranges, California, and central South Island, New Zealand
:
Journal of Geophysical Research
 , v.
100
, p.
18,059
18,074
, https://doi.org/10.1029/95JB01683.
21.
Burbank
,
D.W.
,
Blythe
,
A.E.
,
Putkonen
,
J.
,
Pratt-Sitaula
,
B.
,
Gabet
,
E.
,
Oskin
,
M.
,
Barros
,
A.
, and
Ojha
,
T.P.
,
2003
,
Decoupling of erosion and precipitation in the Himalayas
:
Nature
 , v.
426
, p.
652
655
, https://doi.org/10.1038/nature02187.
22.
Caldwell
,
W.B.
,
Klemperer
,
S.L.
,
Lawrence
,
J.F.
, and
Rai
,
S.S.
,
2013
,
Characterizing the Main Himalayan Thrust in the Garhwal Himalaya, India with receiver function CCP stacking
:
Earth and Planetary Science Letters
 , v.
367
, p.
15
27
, https://doi.org/10.1016/j.epsl.2013.02.009.
23.
Cannon
,
J.M.
, and
Murphy
,
M.A.
,
2014
,
Active lower-crustal deformation and Himalayan seismic hazard revealed by stream channels and regional geology
:
Tectonophysics
 , v.
633
, p.
34
42
, https://doi.org/10.1016/j.tecto.2014.06.031.
24.
Carosi
,
R.
,
Montomoli
,
C.
,
Rubatto
,
D.
, and
Visonà
,
D.
,
2010
,
Late Oligocene high-temperature shear zones in the core of the Higher Himalayan Crystallines (Lower Dolpo, western Nepal)
:
Tectonics
 , v.
29
, TC4029, https://doi.org/10.1029/2008TC002400.
25.
Cattin
,
R.
, and
Avouac
,
J.P.
,
2000
,
Modeling mountain building and the seismic cycle in the Himalaya of Nepal
:
Journal of Geophysical Research
 , v.
105
, p.
13,389
13,407
, https://doi.org/10.1029/2000JB900032.
26.
Clark
,
M.K.
, and
Bilham
,
R.
,
2008
,
Miocene rise of the Shillong Plateau and the beginning of the end for the Eastern Himalaya
:
Earth and Planetary Science Letters
 , v.
269
, p.
337
351
, https://doi.org/10.1016/j.epsl.2008.01.045.
27.
Cooper
,
F.J.
,
Hodges
,
K.V.
,
Parrish
,
R.R.
,
Roberts
,
N.M.W.
, and
Horstwood
,
M.S.A.
,
2015
,
Synchronous N-S and E-W extension at the Tibet-to-Himalaya transition in NW Bhutan
:
Tectonics
 , v.
34
, p.
1375
1395
, https://doi.org/10.1002/2014TC003712.
28.
Cyr
,
A.J.
,
Granger
,
D.E.
,
Olivetti
,
V.
, and
Molin
,
P.
,
2010
,
Quantifying rock uplift rates using channel steepness and cosmogenic nuclide–determined erosion rates: Examples from northern and southern Italy
:
Lithosphere
 , v.
2
, p.
188
198
, https://doi.org/10.1130/L96.1.
29.
Davis
,
W.M.
,
1902
,
Baselevel, grade and peneplain
:
The Journal of Geology
 , v.
10
, p.
77
111
, https://doi.org/10.1086/620982.
30.
DeCelles
,
P.G.
,
Gehrels
,
G.E.
,
Quade
,
J.
,
Ojha
,
T.P.
,
Kapp
,
P.A.
, and
Upreti
,
B.N.
,
1998
,
Neogene foreland basin deposits, erosional unroofing, and the kinematic history of the Himalayan fold-thrust belt, western Nepal
:
Geological Society of America Bulletin
 , v.
110
, p.
2
21
, https://doi.org/10.1130/0016-7606(1998)110<0002:NFBDEU>2.3.CO;2.
31.
Drukpa
,
D.
,
Velasco
,
A.A.
, and
Doser
,
D.I.
,
2006
,
Seismicity in the Kingdom of Bhutan (1937–2003): Evidence for crustal transcurrent deformation
:
Journal of Geophysical Research
 , v.
111
, B06301, https://doi.org/10.1029/2004JB003087.
32.
Duvall
,
A.
,
Kirby
,
E.
, and
Burbank
,
D.
,
2004
,
Tectonic and lithologic controls on bedrock channel profiles and processes in coastal California
:
Journal of Geophysical Research
 , v.
109
, F03002, https://doi.org/10.1029/2003JF000086.
33.
Edwards
,
M.A.
,
Pêcher
,
A.
,
Kidd
,
W.S.F.
,
Burchfiel
,
B.C.
, and
Royden
,
L.H.
,
1999
,
Southern Tibet Detachment System at Khula Kangri, Eastern Himalaya: A large-area, shallow detachment stretching into Bhutan?
:
The Journal of Geology
 , v.
107
, p.
623
631
, https://doi.org/10.1086/314366.
34.
Elliott
,
J.R.
,
Jolivet
,
R.
,
González
,
P.J.
,
Avouac
,
J.-P.
,
Hollingsworth
,
J.
,
Searle
,
M.P.
, and
Stevens
,
V.L.
,
2016
,
Himalayan megathrust geometry and relation to topography revealed by the Gorkha earthquake
:
Nature Geoscience
 , v.
9
, p.
174
180
, https://doi.org/10.1038/ngeo2623.
35.
Ellis
,
S.
,
Fullsack
,
P.
, and
Beaumont
,
C.
,
1995
,
Oblique convergence of the crust driven by basal forcing: Implications for length-scales of deformation and strain partitioning in orogens
:
Geophysical Journal International
 , v.
120
, p.
24
44
, https://doi.org/10.1111/j.1365-246X.1995.tb05909.x.
36.
Feldl
,
N.
, and
Bilham
,
R.
,
2006
,
Great Himalayan earthquakes and the Tibetan plateau
:
Nature
 , v.
444
, p.
165
170
, https://doi.org/10.1038/nature05199.
37.
Forte
,
A.M.
,
Whipple
,
K.X.
,
Bookhagen
,
B.
, and
Rossi
,
M.W.
,
2016
,
Decoupling of modern shortening rates, climate, and topography in the Caucasus
:
Earth and Planetary Science Letters
 , v.
449
, p.
282
294
, https://doi.org/10.1016/j.epsl.2016.06.013.
38.
Fuchs
,
G.R.
,
1964
,
Note on the geology of the Palaeozoics and Mesozoics of the Tibetan zone of the Dolpo region (Nepal–Himalaya)
:
Sonderabdruck aus den Verhandlungen der Geologischen Bundesanstalt
 , v.
1
, p.
6
9
.
39.
Gahalaut
,
V.K.
, and
Kundu
,
B.
,
2012
,
Possible influence of subducting ridges on the Himalayan arc and on the ruptures of great and major Himalayan earthquakes
:
Gondwana Research
 , v.
21
, p.
1080
1088
, https://doi.org/10.1016/j.gr.2011.07.021.
40.
Gibson
,
R.
,
Godin
,
L.
,
Kellett
,
D.A.
,
Cottle
,
J.M.
, and
Archibald
,
D.
,
2016
,
Diachronous deformation along the base of the Himalayan metamorphic core, west-central Nepal
:
Geological Society of America Bulletin
 , v.
128
, p.
860
878
, https://doi.org/10.1130/B31328.1.
41.
Godard
,
V.
,
Bourlès
,
D.L.
,
Spinabella
,
F.
,
Burbank
,
D.W.
,
Bookhagen
,
B.
,
Fisher
,
G.B.
,
Moulin
,
A.
, and
Léanni
,
L.
,
2014
,
Dominance of tectonics over climate in Himalayan denudation
:
Geology
 , v.
42
, p.
243
246
, https://doi.org/10.1130/G35342.1.
42.
Godin
,
L.
,
Grujic
,
D.
,
Law
,
R.D.
, and
Searle
,
M.P.
,
2006b
,
Channel flow, ductile extrusion and exhumation in continental collision zones: An introduction
, in
Law
,
R.D.
,
Searle
,
M.P.
, and
Godin
,
L.
, eds.,
Channel Flow, Ductile Extrusion and Exhumation in Continental Collision Zones: Geological Society of London Special Publication 268
 , p.
1
23
, https://doi.org/10.1144/GSL.SP.2006.268.01.01.
43.
Grandin
,
R.
,
Doin
,
M.-P.
,
Bollinger
,
L.
,
Pinel-Puysségur
,
B.
,
Ducret
,
G.
,
Jolivet
,
R.
, and
Sapkota
,
S.N.
,
2012
,
Long-term growth of the Himalaya inferred from interseismic InSAR measurement
:
Geology
 , v.
40
, p.
1059
1062
, https://doi.org/10.1130/G33154.1.
44.
Grandin
,
R.
,
Vallée
,
M.
,
Satriano
,
C.
,
Lacassin
,
R.
,
Klinger
,
Y.
,
Simoes
,
M.
, and
Bollinger
,
L.
,
2015
,
Rupture process of the Mw = 7.9 2015 Gorkha earthquake (Nepal): Insights into Himalayan megathrust segmentation
:
Geophysical Research Letters
 , v.
42
, p.
8373
8382
, https://doi.org/10.1002/2015GL066044.
45.
Grujic
,
D.
,
Hollister
,
L.S.
, and
Parrish
,
R.R.
,
2002
,
Himalayan metamorphic sequence as an orogenic channel: Insight from Bhutan
:
Earth and Planetary Science Letters
 , v.
198
, p.
177
191
, https://doi.org/10.1016/S0012-821X(02)00482-X.
46.
Grujic
,
D.
,
Coutand
,
I.
,
Bookhagen
,
B.
,
Bonnet
,
S.
,
Blythe
,
A.
, and
Duncan
,
C.
,
2006
,
Climatic forcing of erosion, landscape, and tectonics in the Bhutan Himalayas
:
Geology
 , v.
34
, p.
801
804
, https://doi.org/10.1130/G22648.1.
47.
Hack
,
J.T.
,
1957
,
Studies of longitudinal stream profiles in Virginia and Maryland
:
U.S. Geological Survey Professional Paper 294-B
 ,
53
p.
48.
Harvey
,
J.E.
,
Burbank
,
D.W.
, and
Bookhagen
,
B.
,
2015
,
Along-strike changes in Himalayan thrust geometry: Topographic and tectonic discontinuities in western Nepal
:
Lithosphere
 , v.
7
, p.
511
518
, https://doi.org/10.1130/L444.1.
49.
Hauck
,
M.
,
Nelson
,
K.
,
Brown
,
L.
,
Zhao
,
W.
, and
Ross
,
A.
,
1998
,
Crustal structure of the Himalayan orogen at ∼90° east longitude from Project INDEPTH deep reflection profiles
:
Tectonics
 , v.
17
, p.
481
500
, https://doi.org/10.1029/98TC01314.
50.
He
,
D.
,
Webb
,
A.A.G.
,
Larson
,
K.P.
,
Martin
,
A.J.
, and
Schmitt
,
A.K.
,
2015
,
Extrusion vs. duplexing models of Himalayan mountain building 3: Duplexing dominates from the Oligocene to Present
:
International Geology Review
 , v.
57
, p.
1
27
, https://doi.org/10.1080/00206814.2014.986669.
51.
Herman
,
F.
,
Copeland
,
P.
,
Avouac
,
J.-P.
,
Bollinger
,
L.
,
Mahéo
,
G.
,
Le Fort
,
P.
,
Rai
,
S.
,
Foster
,
D.
,
Pêcher
,
A.
, and
Stüwe
,
K.
,
2010
,
Exhumation, crustal deformation, and thermal structure of the Nepal Himalaya derived from the inversion of thermochronological and thermobarometric data and modeling of the topography
:
Journal of Geophysical Research
 , v.
115
, B06407, https://doi.org/10.1029/2008JB006126.
52.
Hijmans
,
R.J.
,
Cameron
,
S.E.
,
Parra
,
J.L.
,
Jones
,
P.G.
, and
Jarvis
,
A.
,
2005
,
Very high resolution interpolated climate surfaces for global land areas
:
International Journal of Climatology
 , v.
25
, p.
1965
1978
, https://doi.org/10.1002/joc.1276.
53.
Hirth
,
G.
, and
Tullis
,
J.
,
1992
,
Dislocation creep regimes in quartz aggregates
:
Journal of Structural Geology
 , v.
14
, p.
145
159
, https://doi.org/10.1016/0191-8141(92)90053-Y.
54.
Hodges
,
K.V.
,
Wobus
,
C.
,
Ruhl
,
K.
,
Schildgen
,
T.
, and
Whipple
,
K.
,
2004
,
Quaternary deformation, river steepening, and heavy precipitation at the front of the Higher Himalayan ranges
:
Earth and Planetary Science Letters
 , v.
220
, p.
379
389
, https://doi.org/10.1016/S0012-821X(04)00063-9.
55.
Hoek
,
E.
, and
Brown
,
E.T.
,
1980
,
Empirical strength criterion for rock masses, Journal of the Geotechnical Engineering Division (American Society of Civil Engineers)
, v.
106
, p.
1013
1035
.
56.
Hoek
,
E.
, and
Brown
,
E.T.
,
1997
,
Practical estimates of rock mass strength
:
International Journal of Rock Mechanics and Mining Sciences
 , v.
34
, p.
1165
1186
, https://doi.org/10.1016/S1365-1609(97)80069-X.
57.
Hossler
,
T.
,
Bollinger
,
L.
,
Sapkota
,
S.N.
,
Lavé
,
J.
,
Gupta
,
R.M.
, and
Kandel
,
T.P.
,
2016
,
Surface ruptures of large Himalayan earthquakes in Western Nepal: Evidence along a reactivated strand of the Main Boundary Thrust
:
Earth and Planetary Science Letters
 , v.
434
, p.
187
196
, https://doi.org/10.1016/j.epsl.2015.11.042.
58.
Hoste-Colomer
,
R.
,
Bollinger
,
L.
,
Lyon-Caen
,
H.
,
Burtin
,
A.
, and
Adhikari
,
L.B.
,
2016
,
Lateral structure variations and transient swarm revealed by seismicity along the Main Himalayan Thrust north of Kathmandu
:
Tectonophysics
 , v.
714–715
, p.
107
116
, https://doi.org/10.1016/j.tecto.2016.10.004.
59.
Hubbard
,
J.
,
Almeida
,
R.
,
Foster
,
A.
,
Sapkota
,
S.N.
,
Bürgi
,
P.
, and
Tapponnier
,
P.
,
2016
,
Structural segmentation controlled the 2015 Mw 7.8 Gorkha earthquake rupture in Nepal
:
Geology
 , v.
44
, p.
639
642
, https://doi.org/10.1130/G38077.1.
60.
Jackson
,
M.E.
, and
Bilham
,
R.
,
1994
,
1991–1992 GPS measurements across the Nepal Himalaya
:
Geophysical Research Letters
 , v.
21
, p.
1169
1172
, https://doi.org/10.1029/94GL00917.
61.
Jouanne
,
F.
,
Mugnier
,
J.L.
,
Pandey
,
M.R.
,
Gamond
,
J.F.
,
Le Fort
,
P.
,
Serrurier
,
L.
,
Vigny
,
C.
, and
Avouac
,
J.P.
,
1999
,
Oblique convergence in the Himalayas of western Nepal deduced from preliminary results of GPS measurements
:
Geophysical Research Letters
 , v.
26
, p.
1933
1936
, https://doi.org/10.1029/1999GL900416.
62.
Jouanne
,
F.
,
Mugnier
,
J.-L.
,
Gamond
,
J.F.
,
Le Fort
,
P.
,
Pandey
,
M.R.
,
Bollinger
,
L.
,
Flouzat
,
M.
, and
Avouac
,
J.P.
,
2004
,
Current shortening across the Himalayas of Nepal
:
Geophysical Journal International
 , v.
157
, p.
1
14
, https://doi.org/10.1111/j.1365-246X.2004.02180.x.
63.
Kellett
,
D.A
,
Grujic
,
D.
,
Warren
,
C.
,
Cottle
,
J.
,
Jamieson
,
R.
, and
Tenzin
,
T.
,
2010
,
Metamorphic history of a syn-convergent orogen-parallel detachment: The South Tibetan detachment system, Bhutan Himalaya
:
Journal of Metamorphic Geology
 , v.
28
, p.
785
808
, https://doi.org/10.1111/j.1525-1314.2010.00893.x.
64.
Kirby
,
E.
, and
Whipple
,
K.
,
2001
,
Quantifying differential rock-uplift rates via stream profile analysis
:
Geology
 , v.
29
, p.
415
418
, https://doi.org/10.1130/0091-7613(2001)029<0415:QDRURV>2.0.CO;2.
65.
Kirby
,
E.
, and
Whipple
,
K.X.
,
2012
,
Expression of active tectonics in erosional landscapes
:
Journal of Structural Geology
 , v.
44
, p.
54
75
, https://doi.org/10.1016/j.jsg.2012.07.009.
66.
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
:
Journal of Geophysical Research
 , v.
108
, 2217, https://doi.org/10.1029/2001JB000861.
67.
Korup
,
O.
,
2006
,
Rock-slope failure and the river long profile
:
Geology
 , v.
34
, p.
45
48
, https://doi.org/10.1130/G21959.1.
68.
Kumar
,
S.
,
Wesnousky
,
S.G.
,
Rockwell
,
T.K.
,
Briggs
,
R.W.
,
Thakur
,
V.C.
, and
Jayangondaperumal
,
R.
,
2006
,
Paleoseismic evidence of great surface rupture earthquakes along the Indian Himalaya
:
Journal of Geophysical Research
 , v.
111
, B03304, https://doi.org/10.1029/2004JB003309.
69.
Kumar
,
S.
,
Wesnousky
,
S.G.
,
Jayangondaperumal
,
R.
,
Nakata
,
T.
,
Kumahara
,
Y.
, and
Singh
,
V.
,
2010
,
Paleoseismological evidence of surface faulting along the northeastern Himalayan front, India: Timing, size, and spatial extent of great earthquakes
:
Journal of Geophysical Research
 , v.
115
, p.
1
20
, https://doi.org/10.1029/2009JB006789.
70.
Kundu
,
B.
,
Yadav
,
R.K.
,
Bali
,
B.S.
,
Chowdhury
,
S.
, and
Gahalaut
,
V.K.
,
2014
,
Oblique convergence and slip partitioning in the NW Himalaya
:
Implications from GPS measurements: Tectonics
 , v.
33
, p.
2013
2024
, https://doi.org/10.1002/2014TC003633.
71.
Larson
,
K.P.
, and
Cottle
,
J.M.
,
2014
,
Midcrustal discontinuities and the assembly of the Himalayan midcrust
:
Tectonics
 , v.
33
, p.
718
740
, https://doi.org/10.1002/2013TC003452.
72.
Larson
,
K.P.
,
Ambrose
,
T.K.
,
Webb
,
A.G.
,
Cottle
,
J.M.
, and
Shrestha
,
S.
,
2015
,
Reconciling Himalayan midcrustal discontinuities: The Main Central thrust system
:
Earth and Planetary Science Letters
 , v.
429
, p.
139
146
, https://doi.org/10.1016/j.epsl.2015.07.070.
73.
Lavé
,
J.
, and
Avouac
,
J.P.
,
2001
,
Fluvial incision and tectonic uplift across the Himalayas of central Nepal
:
Journal of Geophysical Research
 , v.
106
, p.
26,561
26,591
, https://doi.org/10.1029/2001JB000359.
74.
Lemonnier
,
C.
,
Marquis
,
G.
,
Perrier
,
F.
,
Avouac
,
J.-P.
,
Chitrakar
,
G.
,
Kafle
,
B.
,
Sapkota
,
S.
,
Gautam
,
U.
,
Tiwari
,
D.
, and
Bano
,
M.
,
1999
,
Electrical structure of the Himalaya of central Nepal: High conductivity around the mid-crustal ramp along the MHT
:
Geophysical Research Letters
 , v.
26
, p.
3261
3264
, https://doi.org/10.1029/1999GL008363.
75.
Lindsey
,
E.O.
,
Natsuaki
,
R.
,
Xu
,
X.
,
Shimada
,
M.
,
Hashimoto
,
M.
,
Melgar
,
D.
, and
Sandwell
,
D.T.
,
2015
,
Line-of-sight displacement from ALOS-2 interferometry: Mw 7.8 Gorkha Earthquake and Mw 7.3 aftershock
:
Geophysical Research Letters
 , v.
42
, p.
6655
6661
, https://doi.org/10.1002/2015GL065385.
76.
Long
,
S.
,
McQuarrie
,
N.
,
Tobgay
,
T.
,
Grujic
,
D.
, and
Hollister
,
L.
,
2011
,
Geologic map of Bhutan
:
Journal of Maps
 , v.
7
, p.
184
192
, https://doi.org/10.4113/jom.2011.1159.
77.
Mackin
,
J.H.
,
1948
,
Concept of the graded river
:
Geological Society of America Bulletin
 , v.
59
, p.
463
512
, https://doi.org/10.1130/0016-7606(1948)59[463:COTGR]2.0.CO;2.
78.
McCaffrey
,
R.
,
1992
,
Oblique plate convergence, slip vectors, and forearc deformation
:
Journal of Geophysical Research
 , v.
97
, p.
8905
8915
, https://doi.org/10.1029/92JB00483.
79.
McCaffrey
,
R.
, and
Nabelek
,
J.
,
1998
,
Role of oblique convergence in the active deformation of the Himalayas and southern Tibet plateau
:
Geology
 , v.
26
, p.
691
694
, https://doi.org/10.1130/0091-7613(1998)026<0691:ROOCIT>2.3.CO;2.
80.
McCaffrey
,
R.
,
Zwick
,
P.C.
,
Bock
,
Y.
,
Prawirodirdjo
,
L.
,
Genrich
,
J.F.
,
Stevens
,
C.W.
,
Puntodewo
,
S.S.O.
, and
Subarya
,
C.
,
2000
,
Strain partitioning during oblique plate convergence in northern Sumatra: Geodetic and seismologic constraints and numerical modeling
:
Journal of Geophysical Research
 , v.
105
, p.
28,363
28,376
, https://doi.org/10.1029/1999JB900362.
81.
McCallister
,
A.T.
,
Taylor
,
M.H.
,
Murphy
,
M.A.
,
Styron
,
R.H.
, and
Stockli
,
D.F.
,
2014
,
Thermochronologic constraints on the late Cenozoic exhumation history of the Gurla Mandhata metamorphic core complex, Southwestern Tibet
:
Tectonics
 , v.
33
, p.
27
52
, https://doi.org/10.1002/2013TC003302.
82.
McDermott
,
J.A.
,
Whipple
,
K.X.
,
Hodges
,
K.V.
, and
Soest
,
M.C.
,
2013
,
Evidence for Plio-Pleistocene north-south extension at the southern margin of the Tibetan Plateau, Nyalam region
:
Tectonics
 , v.
32
, p.
317
333
, https://doi.org/10.1002/tect.20018.
83.
McDermott
,
J.A.
,
Hodges
,
K.V.
,
Whipple
,
K.X.
,
Van Soest
,
M.C.
, and
Hurtado
,
J.M.
, Jr
.,
2015
,
Evidence for Pleistocene low-angle normal faulting in the Annapurna-Dhaulagiri region, Nepal
:
The Journal of Geology
 , v.
123
, p.
133
151
, https://doi.org/10.1086/681219.
84.
McQuarrie
,
N.
, and
Ehlers
,
T.A.
,
2015
,
Influence of thrust belt geometry and shortening rate on thermochronometer cooling ages: Insights from thermokinematic and erosion modeling of the Bhutan Himalaya
:
Tectonics
 , v.
34
, p.
1055
1079
, https://doi.org/10.1002/2014TC003783.
85.
McQuarrie
,
N.
,
Tobgay
,
T.
,
Long
,
S.P.
,
Reiners
,
P.W.
, and
Cosca
,
M.A.
,
2014
,
Variable exhumation rates and variable displacement rates: Documenting recent slowing of Himalayan shortening in western Bhutan
:
Earth and Planetary Science Letters
 , v.
386
, p.
161
174
, https://doi.org/10.1016/j.epsl.2013.10.045.
86.
Meade
,
B.J.
,
2010
,
The signature of an unbalanced earthquake cycle in Himalayan topography?
:
Geology
 , v.
38
, p.
987
990
, https://doi.org/10.1130/G31439.1.
87.
Mitra
,
G.
,
Bhattacharyya
,
K.
, and
Mukul
,
M.
,
2010
,
The Lesser Himalayan duplex in Sikkim: Implications for variations in Himalayan shortening
:
Journal of the Geological Society of India
 , v.
75
, p.
289
301
, https://doi.org/10.1007/s12594-010-0016-x.
88.
Montomoli
,
C.
,
Iaccarino
,
S.
,
Carosi
,
R.
,
Langone
,
A.
, and
Visonà
,
D.
,
2013
,
Tectonometamorphic discontinuities within the Greater Himalayan Sequence in Western Nepal (Central Himalaya): Insights on the exhumation of crystalline rocks
:
Tectonophysics
 , v.
608
, p.
1349
1370
, https://doi.org/10.1016/j.tecto.2013.06.006.
89.
Morell
,
K.D.
,
Sandiford
,
M.
,
Rajendran
,
C.P.
,
Rajendran
,
K.
,
Alimanovic
,
A.
,
Fink
,
D.
, and
Sanwal
,
J.
,
2015
,
Geomorphology reveals active décollement geometry in the central Himalayan seismic gap
:
Lithosphere
 , v.
7
, p.
247
256
, https://doi.org/10.1130/L407.1.
90.
Morell
,
K.D.
,
Sandiford
,
M.
,
Kohn
,
B.
,
Codilean
,
A.
,
Fülöp
,
R.H.
, and
Ahmad
,
T.
,
2017
,
Current strain accumulation in the hinterland of the northwest Himalaya constrained by landscape analyses, basin-wide denudation rates, and low temperature thermochronology
:
Tectonophysics
 , v.
721
, p.
70
89
, https://doi.org/10.1016/j.tecto.2017.09.007.
91.
Mugnier
,
J.L.
,
Leturmy
,
P.
,
Mascle
,
G.
,
Huyghe
,
P.
,
Chalaron
,
E.
,
Vidal
,
G.
,
Husson
,
L.
, and
Delcaillau
,
B.
,
1999
,
The Siwaliks of western Nepal: I. Geometry and kinematics
:
Journal of Asian Earth Sciences
 , v.
17
, p.
629
642
, https://doi.org/10.1016/S1367-9120(99)00038-3.
92.
Mugnier
,
J.L.
,
Gajurel
,
A.
,
Huyghe
,
P.
,
Jayangondaperumal
,
R.
,
Jouanne
,
F.
, and
Upreti
,
B.
,
2013
,
Structural interpretation of the great earthquakes of the last millennium in the central Himalaya
:
Earth-Science Reviews
 , v.
127
, p.
30
47
, https://doi.org/10.1016/j.earscirev.2013.09.003.
93.
Murphy
,
M.A.
,
2007
,
Isotopic characteristics of the Gurla Mandhata metamorphic core complex: Implications for the architecture of the Himalayan orogen
:
Geology
 , v.
35
, p.
983
986
, https://doi.org/10.1130/G23774A.1.
94.
Murphy
,
M.A.
, and
Copeland
,
P.
,
2005
,
Transtensional deformation in the central Himalaya and its role in accommodating growth of the Himalayan orogen
:
Tectonics
 , v.
24
, TC4012, https://doi.org/10.1029/2004TC001659.
95.
Murphy
,
M.A.
, and
Harrison
,
T.M.
,
1999
,
Relationship between leucogranites and the Qomolangma detachment in the Rongbuk Valley, south Tibet
:
Geology
 , v.
27
, p.
831
834
, https://doi.org/10.1130/0091-7613(1999)027<0831:RBLATQ>2.3.CO;2.
96.
Murphy
,
M.A.
,
Yin
,
A.
,
Kapp
,
P.
,
Harrison
,
T.M.
,
Manning
,
C.E.
,
Ryerson
,
F.J.
,
Lin
,
D.
, and
Jinghui
,
G.
,
2002
,
Structural evolution of the Gurla Mandhata detachment system, southwest Tibet: Implica0074ions for the eastward extent of the Karakoram fault system
:
Geological Society of America Bulletin
 , v.
114
, p.
428
447
, https://doi.org/10.1130/0016-7606(2002)114<0428:SEOTGM>2.0.CO;2.
97.
Murphy
,
M.A.
,
Taylor
,
M.H.
,
Gosse
,
J.
,
Silver
,
C.R.P.
,
Whipp
,
D.M.
, and
Beaumont
,
C.
,
2014
,
Limit of strain partitioning in the Himalaya marked by large earthquakes in western Nepal
:
Nature Geoscience
 , v.
7
, p.
38
42
, https://doi.org/10.1038/ngeo2017.
98.
Nábělek
,
J.
,
Hetényi
,
G.
,
Vergne
,
J.
,
Sapkota
,
S.
,
Kafle
,
B.
,
Jiang
,
M.
,
Su
,
H.
,
Chen
,
J.
,
Huang
,
B.-S.
, and
the HiCLIMB Team
,
2009
,
Underplating in the Himalaya-Tibet collision zone revealed by the Hi-CLIMB experiment
:
Science
 , v.
325
, p.
1371
1374
, https://doi.org/10.1126/science.1167719.
99.
Najman
,
Y.
,
Bracciali
,
L.
,
Parrish
,
R.R.
,
Chisty
,
E.
, and
Copley
,
A.
,
2016
,
Evolving strain partitioning in the Eastern Himalaya: The growth of the Shillong Plateau
:
Earth and Planetary Science Letters
 , v.
433
, p.
1
9
, https://doi.org/10.1016/j.epsl.2015.10.017.
100.
Olen
,
S.M.
,
Bookhagen
,
B.
,
Hoffmann
,
B.
,
Sachse
,
D.
,
Adhikari
,
D.P.
, and
Strecker
,
M.R.
,
2015
,
Understanding erosion rates in the Himalayan orogen: A case study from the Arun Valley: Journal of Geophysical Research
:
Earth Surface
 , v.
120
, p.
2080
2102
, https://doi.org/10.1002/2014JF003410.
101.
Owen
,
L.A.
, and
Benn
,
D.I.
,
2005
,
Equilibrium-line altitudes of the Last Glacial Maximum for the Himalaya and Tibet: An assessment and evaluation of results
:
Quaternary International
 , v.
138
, p.
55
78
, https://doi.org/10.1016/j.quaint.2005.02.006.
102.
Pandey
,
M.R.
,
Tandukar
,
R.P.
,
Avouac
,
J.P.
,
Lavé
,
J.
, and
Massot
,
J.P.
,
1995
,
Interseismic strain accumulation on the Himalayan crustal ramp (Nepal)
:
Geophysical Research Letters
 , v.
22
, p.
751
754
, https://doi.org/10.1029/94GL02971.
103.
Pandey
,
M.R.
,
Tandukar
,
R.P.
,
Avouac
,
J.P.
,
Vergne
,
J.
, and
Héritier
,
T.
,
1999
,
Seismotectonics of the Nepal Himalaya from a local seismic network
:
Journal of Asian Earth Sciences
 , v.
17
, p.
703
712
, https://doi.org/10.1016/S1367-9120(99)00034-6.
104.
Rajendran
,
C.P.
,
John
,
B.
, and
Rajendran
,
K.
,
2015
,
Medieval pulse of great earthquakes in the central Himalaya: Viewing past activities on the frontal thrust: Journal of Geophysical Research
:
Solid Earth
 , v.
120
, p.
1623
1641
, https://doi.org/10.1002/2014JB011015.
105.
Robinson
,
D.M.
,
2008
,
Forward modeling the kinematic sequence of the central Himalayan thrust belt, western Nepal
:
Geosphere
 , v.
4
, p.
785
801
, https://doi.org/10.1130/GES00163.1.
106.
Robl
,
J.
,
Stüwe
,
K.
, and
Hergarten
,
S.
,
2008
,
Channel profiles around Himalayan river anticlines: Constraints on their formation from digital elevation model analysis
:
Tectonics
 , v.
27
, TC3010, https://doi.org/10.1029/2007TC002215.
107.
Sapkota
,
S.N.
,
Bollinger
,
L.
,
Klinger
,
Y.
,
Tapponnier
,
P.
,
Gaudemer
,
Y.
, and
Tiwari
,
D.
,
2013
,
Primary surface ruptures of the great Himalayan earthquakes in 1934 and 1255
:
Nature Geoscience
 , v.
6
, p.
71
76
, https://doi.org/10.1038/ngeo1669.
108.
Schelling
,
D.
, and
Arita
,
K.
,
1991
,
Thrust tectonics, crustal shortening, and the structure of the far-eastern Nepal Himalaya
:
Tectonics
 , v.
10
, p.
851
862
, https://doi.org/10.1029/91TC01011.
109.
Searle
,
M.P.
, and
Godin
,
L.
,
2003
,
The south Tibetan detachment and the Manaslu leucogranite: A structural reinterpretation and restoration of the Annapurna-Manaslu Himalaya, Nepal
:
The Journal of Geology
 , v.
111
, p.
505
523
, https://doi.org/10.1086/376763.
110.
Seeber
,
L.
, and
Gornitz
,
V.
,
1983
,
River profiles along the Himalayan arc as indicators of active tectonics
:
Tectonophysics
 , v.
92
, p.
335
367
, https://doi.org/10.1016/0040-1951(83)90201-9.
111.
Silver
,
C.R.
,
Murphy
,
M.A.
,
Taylor
,
M.H.
,
Gosse
,
J.
, and
Baltz
,
T.
,
2015
,
Neotectonics of the Western Nepal Fault System: Implications for Himalayan strain partitioning
:
Tectonics
 , v.
34
, p.
2494
2513
, https://doi.org/10.1002/2014TC003730.
112.
Stevens
,
V.L.
, and
Avouac
,
J.P.
,
2015
,
Interseismic coupling on the main Himalayan thrust
:
Geophysical Research Letters
 , v.
42
, p.
5828
5837
, https://doi.org/10.1002/2015GL064845.
113.
Stevens
,
V.L.
, and
Avouac
,
J.P.
,
2016
,
Millenary Mw > 9.0 earthquakes required by geodetic strain in the Himalaya
:
Geophysical Research Letters
 , v.
43
, p.
1118
1123
, https://doi.org/10.1002/2015GL067336.
114.
Styron
,
R.H.
,
Taylor
,
M.H.
, and
Murphy
,
M.A.
,
2011
,
Oblique convergence, arc-parallel extension, and the role of strike-slip faulting in the High Himalaya
:
Geosphere
 , v.
7
, p.
582
596
, https://doi.org/10.1130/GES00606.1.
115.
Sundell
,
K.E.
,
Taylor
,
M.H.
,
Styron
,
R.H.
,
Stockli
,
D.F.
,
Kapp
,
P.
,
Hager
,
C.
,
Liu
,
D.
, and
Ding
,
L.
,
2013
,
Evidence for constriction and Pliocene acceleration of east-west extension in the North Lunggar rift region of west central Tibet
:
Tectonics
 , v.
32
, p.
1454
1479
, https://doi.org/10.1002/tect.20086.
116.
Taylor
,
M.H.
,
2016
,
Tectonics: Tales of Himalayan topography
:
Nature Geoscience
 , v.
9
, p.
649
651
, https://doi.org/10.1038/ngeo2805.
117.
Taylor
,
M.
, and
Yin
,
A.
,
2009
,
Active structures of the Himalayan-Tibetan orogen and their relationships to earthquake distribution, contemporary strain field, and Cenozoic volcanism
:
Geosphere
 , v.
5
, p.
199
214
, https://doi.org/10.1130/GES00217.1.
118.
Thiede
,
R.C.
,
Arrowsmith
,
J.R.
,
Bookhagen
,
B.
,
McWilliams
,
M.
,
Sobel
,
E.R.
, and
Strecker
,
M.R.
,
2006
,
Dome formation and extension in the Tethyan Himalaya, Leo Pargil, northwest India
:
Geological Society of America Bulletin
 , v.
118
, p.
635
650
, https://doi.org/10.1130/B25872.1.
119.
Tikoff
,
B.
, and
Teyssier
,
C.
,
1994
,
Strain modeling of displacement-field partitioning in transpressional orogens
:
Journal of Structural Geology
 , v.
16
, p.
1575
1588
, https://doi.org/10.1016/0191-8141(94)90034-5.
120.
Tullis
,
J.
, and
Yund
,
R.
,
1992
,
The brittle-ductile transition in feldspar aggregates: An experimental study
:
International Geophysics
 , v.
51
, p.
89
117
, https://doi.org/10.1016/S0074-6142(08)62816-8.
121.
Tyagi
,
A.K.
,
Chaudhary
,
S.
,
Rana
,
N.
,
Sati
,
S.P.
, and
Juyal
,
N.
,
2009
,
Identifying areas of differential uplift using steepness index in the Alaknanda basin, Garhwal Himalaya, Uttarakhand
:
Current Science
 , v.
97
, p.
1473
1477
.
122.
Vannay
,
J.C.
,
Grasemann
,
B.
,
Rahn
,
M.
,
Frank
,
W.
,
Carter
,
A.
,
Baudraz
,
V.
, and
Cosca
,
M.
,
2004
,
Miocene to Holocene exhumation of metamorphic crustal wedges in the NW Himalaya: Evidence for tectonic extrusion coupled to fluvial erosion
:
Tectonics
 , v.
23
, TC1014, https://doi.org/10.1029/2002TC001429.
123.
Velasco
,
A.A.
,
Gee
,
V.L.
,
Rowe
,
C.
,
Grujic
,
D.
,
Hollister
,
L.S.
,
Hernandez
,
D.
,
Miller
,
K.C.
,
Tobgay
,
T.
,
Fort
,
M.
, and
Harder
,
S.
,
2007
,
Using small, temporary seismic networks for investigating tectonic deformation: Brittle deformation and evidence for strike-slip faulting in Bhutan
:
Seismological Research Letters
 , v.
78
, p.
446
453
, https://doi.org/10.1785/gssrl.78.4.446.
124.
Webb
,
A.A.G.
,
2013
,
Preliminary balanced palinspastic reconstruction of Cenozoic deformation across the Himachal Himalaya (northwestern India)
:
Geosphere
 , v.
9
, p.
572
587
, https://doi.org/10.1130/GES00787.1.
125.
Webb
,
A.A.G.
,
Yin
,
A.
,
Harrison
,
T.M.
,
Célérier
,
J.
, and
Burgess
,
W.P.
,
2007
,
The leading edge of the Greater Himalayan Crystalline complex revealed in the NW Indian Himalaya: Implications for the evolution of the Himalayan orogen
:
Geology
 , v.
35
, p.
955
958
, https://doi.org/10.1130/G23931A.1.
126.
Whipp
,
D.M.
, Jr
,
Beaumont
,
C.
, and
Braun
,
J.
,
2014
,
Feeding the “aneurysm”: Orogen-parallel mass transport into Nanga Parbat and the western Himalayan syntaxis: Journal of Geophysical Research
:
Solid Earth
 , v.
119
, p.
5077
5096
, https://doi.org/10.1002/2013JB010929.
127.
Whipple
,
K.X.
,
2004
,
Bedrock rivers and the geomorphology of active orogens
:
Annual Review of Earth and Planetary Sciences
 , v.
32
, p.
151
185
, https://doi.org/10.1146/annurev.earth.32.101802.120356.
128.
Whipple
,
K.X.
,
Shirzaei
,
M.
,
Hodges
,
K.V.
,
Arrowsmith
,
J.R.
,
2016
,
Active shortening within the Himalayan orogenic wedge implied by the 2015 Gorkha earthquake
:
Nature Geoscience
 , v.
9
, p.
711
716
, https://doi.org/10.1038/ngeo2797.
129.
Wobus
,
C.W.
,
Hodges
,
K.V.
, and
Whipple
,
K.X.
,
2003
,
Has focused denudation sustained active thrusting at the Himalayan topographic front?
:
Geology
 , v.
31
, p.
861
864
, https://doi.org/10.1130/G19730.1.
130.
Wobus
,
C.
,
Heimsath
,
A.
,
Whipple
,
K.
, and
Hodges
,
K.
,
2005
,
Active out-of-sequence thrust faulting in the central Nepalese Himalaya
:
Nature
 , v.
434
, p.
1008
1011
, https://doi.org/10.1038/nature03499.
131.
Wobus
,
C.
,
Whipple
,
K.X.
,
Kirby
,
E.
,
Snyder
,
N.
,
Johnson
,
J.
,
Spyropolou
,
K.
,
Crosby
,
B.
, and
Sheehan
,
D.
,
2006
,
Tectonics from topography: Procedures, promise, and pitfalls
, in
Willett
,
S.D.
,
Hovius
,
N.
,
Brandon
,
M.T.
, and
Fisher
,
D.M.
, eds.,
Tectonics, Climate, and Landscape Evolution: Geological Society of America Special Paper 398
 , p.
55
74
, https://doi.org/10.1130/2006.2398(04).
132.
Wu
,
C.
,
Nelson
,
K.D.
,
Wortman
,
G.
,
Samson
,
S.D.
,
Yue
,
Y.
,
Li
,
J.
,
Kidd
,
W.S.F.
, and
Edwards
,
M.A.
,
1998
,
Yadong cross structure and South Tibetan Detachment in the east central Himalaya (89°–90°E)
:
Tectonics
 , v.
17
, p.
28
45
, https://doi.org/10.1029/97TC03386.
133.
Yin
,
A.
,
2000
,
Mode of Cenozoic east-west extension in Tibet suggesting a common origin of rifts in Asia during the Indo-Asian collision
:
Journal of Geophysical Research
 , v.
105
, p.
21,745
21,759
, https://doi.org/10.1029/2000JB900168.
134.
Yin
,
A.
,
2006
,
Cenozoic tectonic evolution of the Himalayan orogen as constrained by along-strike variation of structural geometry, exhumation history, and foreland sedimentation
:
Earth-Science Reviews
 , v.
76
, p.
1
131
, https://doi.org/10.1016/j.earscirev.2005.05.004.
135.
Yu
,
H.
,
Webb
,
A.A.G.
, and
He
,
D.
,
2015
,
Extrusion vs. duplexing models of Himalayan mountain building 1: Discovery of the Pabbar thrust confirms duplex-dominated growth of the northwestern Indian Himalaya since Mid-Miocene
:
Tectonics
 , v.
34
, p.
313
333
, https://doi.org/10.1002/2014TC003589.
Gold Open Access: This paper is published under the terms of the CC-BY-NC license.

Supplementary data