Forward modeling reconstructions and data derived from the Himalayan thrust belt and the foreland basin of far western Nepal tie the erosional unroofing and associated deposition to the kinematics and age of fault motion. We reproduce the deformation identified at the surface through a forward-propagating, linked fold-and-thrust belt–foreland basin system. This approach permits estimates of the magnitude of erosion at each time step and the extent, depth, and age of the associated foreland basin. The model reconstructions reveal that the units that supplied the sediment to the foreland basin changed through time: 25–13 Ma, erosion of the Tethyan Himalaya; ca. 12 Ma, first exposure of the Greater Himalaya; ca. 11 Ma, first exposure of the Lesser Himalaya. In our model, exposure of Greater Himalaya and Lesser Himalaya rock is associated with the formation of a thrust ramp that cuts through 7 km of footwall Lesser Himalaya stratigraphy and translates >7 km of Lesser Himalaya rock over the ramp, forming a Lesser Himalaya duplex. An increase in structural relief focuses erosion over the region of the ramp and facilitates exposure of Greater Himalaya and Proterozoic Lesser Himalaya rocks. As the Lesser Himalaya ramp propagates southward, more Lesser Himalaya thrust sheets are incorporated into the Lesser Himalaya duplex. Although uniquely dating thrust events is challenging, these model reconstructions allow us to associate time steps with an age of deposition or exhumation. What emerges is a tempo of deformation that varies with time, marked by periods of rapid shortening during propagation of the Main Central thrust, Ramgarh thrust, and middle stages of the development of the Lesser Himalaya duplex (∼25–30 mm/yr). After emplacement of the Ramgarh thrust, early and late stages of Lesser Himalaya duplex development are marked by periods of slow shortening (∼13–14 mm/yr). Although long-term and modern (geodetic) rates of deformation agree at ∼20 mm/yr, rates of shortening through time have varied from 4 to 33 mm/yr.

Fold-and-thrust belts such as the Himalaya are intrinsically linked to their frontal foreland basins. Understanding the geometry and kinematics of this paired system is a requirement for evaluating the interplay between deformation and climate and their influence on erosion and deposition (Koons, 1990; Willett et al., 1993; Burbank et al., 2003; Thiede et al., 2004; Upton et al., 2009). The Himalaya were formed as the Indian plate moved northward and collided with the Asian plate. After initial collision, deformation propagated southward and built an early fold-and-thrust belt in Neoproterozoic–Paleozoic Tibetan or Tethyan Himalaya rocks (Ratschbacher et al., 1994; Murphy and Yin, 2003) before incorporating the highly metamorphic Neoproterozoic–Ordovician Greater Himalayan rocks and the more proximal to India, Proterozoic Lesser Himalayan rocks into the growing orogen (e.g., Gansser, 1964; LeFort, 1975). This evolving system of folded and faulted rocks shed detritus into the adjacent foreland basin, which was then assimilated into the orogenic system as the Foreland Basin sequence and the Subhimalaya (Burbank et al., 1996; DeCelles et al., 1998b, 2004). The Himalayan Mountains continue to shed sediment into the modern Indo-Gangetic foreland basin.

Because the Himalaya is the largest active continent-continent collision on Earth, it is a testing ground for kinematic and dynamic models focused on understanding collisional orogens, the feedback between erosion and tectonics, and the principal drivers that partition deformation and convergence. A common approach to documenting the structural geometry of the Himalaya is through balanced cross sections (e.g., Long et al., 2011, and references therein). These cross sections illustrate a possible geometry of the thrust belt, while the restored sections are a representation of undeformed stratigraphy prior to the Himalayan collision. In these static cross sections, undeformed and deformed sections are presented without showing how the thrust belt developed kinematically. However, implied in every cross section is a kinematic sequence necessary to reproduce the illustrated deformation, which can be shown through a series of sequential reconstructions that are balanced at every step (e.g., Robinson, 2008). Although viable, these modeling reconstructions do not have absolute time assigned to the frames, and no effort is made to account for erosion of the stratigraphic overburden or flexure of the Indian plate at each step in the reconstruction. Determining the time of fault initiation and duration in order to place an absolute time frame in sequential deformation steps is difficult because the age of faulting can rarely be directly measured. However, geochronologic ages of mineral systems provide records of burial, and thermochronologic ages of mineral systems provide records of exhumation. Bracketing the interval between deposition and exhumation provides an indirect age for the initiation of motion on thrust faults. In addition, if a well-dated foreland basin can be tied to a period of time when the fold-and-thrust belt was deforming and eroding, then the foreland basin succession may be used to place age brackets on the deformation (e.g., Najman et al., 2004, 2009; van der Beek et al., 2006; Szulc et al., 2006; Ravikant et al., 2011).

Fold-and-thrust belts can be characterized to a first order by critical taper wedge theory (Chapple, 1978; Davis et al., 1983; Dahlen, 1990). The metamorphosed and penetratively strained core of the Himalayan orogen is thought by many to preclude its inclusion in a coupled thrust belt–foreland basin system or to be viewed as part of a critically tapered system (e.g., Searle et al., 2008). Although Greater Himalaya rocks have a pronounced flattening fabric (Law et al., 2004; Long and McQuarrie, 2010) and are penetratively strained (e.g., Larson et al., 2010), Greater Himalaya rock was emplaced over 125 km of Lesser Himalaya rocks along a singular structure, the Main Central thrust (Brunel, 1986; Schelling and Arita, 1991; Robinson et al., 2003, 2006). Although we recognize the ductile nature of Greater Himalaya rocks, we do not attempt to include penetrative deformation. The first-order constraint of Main Central thrust and Greater Himalaya evolution is the motion on the Main Central thrust required by the southernmost cutoffs of the Greater Himalaya rocks to emplace 125 km of Greater Himalaya over Lesser Himalaya rock. Precedence for including highly metamorphosed and strained “basement” rocks as part of fold-and-thrust belt systems exists in the Pyrenees (Sinclair et al., 2005), northern Scandinavian Caledonides (Northrup, 1996), and Appalachians (Mitra, 1978; Boyer and Elliott, 1982; Hatcher and Hooper, 1992).

The purpose of this study is to link the deformation, erosion, and depositional history of the Himalayan orogen in far western Nepal to forward models of thrust belt deformation through the incorporation of multiple data sets. We constructed a sequentially restored, structural model based on the kinematic sequence inherent in a balanced cross section using the reconstruction program 2D Move to evaluate the co-evolution of the fold-and-thrust belt and its associated foreland basin. Using 2D Move, we modeled the dimensions of a flexural foreland basin by determining the structural elevation and associated isostatic loading from thrust sheet emplacement, followed by isostatic unloading from erosion. By assuming a topographic wedge that mimics the modern increase in elevation from the Indo-Gangetic foreland to the Tibetan Plateau over ∼150 km (∼2° topographic slope), we can estimate the magnitude of rock eroded in each time step and the expected stratigraphy exposed in the thrust belt. Each increment of thrust emplacement deforms the topographic surface. The difference between the deformed topography and the original topography becomes a gravitational load that flexes the Indian plate, increases the angle of the basal décollement, and creates a foreland basin. The difference between the deformed topography and the new topographic ramp from the deformation front to the hinterland defines the amount of material removed by erosion. Although this approach does not account for topographic details, we show that this first-order approximation allows us to predict depths to thermochronologic and geochronologic samples that can be compared with closure temperatures and ages as well as predicting the age, detrital input, and location of an evolving foreland basin. We compare the burial, exhumational, and depositional history of our sequential reconstruction to published geochronologic and thermochronologic data from the hinterland, as well as detrital zircon U-Pb ages, epsilon Nd values, trace-element data, bulk petrography, sedimentology, and magnetostratigraphic data from the foreland basin. This approach documents the tempo of both deformation and exhumation in the Himalayan thrust belt of far western Nepal from 25 Ma to the present, even within the uncertainty inherent in the evolution of topography with time.

Nepal is located in the central part of the Himalayan arc, which extends ∼2400 km between the eastern to western syntaxes (Fig. 1A). Far western Nepal is situated in the apex of the Himalayan arc and is the location of the largest magnitude of horizontal shortening estimated in the Himalaya (Robinson et al., 2006; see review in Long et al., 2011).

The Tibetan part of the Himalayan thrust belt system consists of the high peaks of the Himalaya northward to the Indus suture and contains Tethyan Himalaya rocks (Fig. 1B). Most researchers agree that deformation began at the initiation of the collision of India with Asia at ca. 50–55 Ma (Zhu et al., 2005; Najman, 2006; Green et al., 2008; Dupont-Nivet et al., 2010; Najman et al., 2010; Wang et al., 2011) and continued deforming and shortening Tethyan Himalaya rocks until 30–25 Ma (Fig. 2; e.g., Ratschbacher et al., 1994; Webb et al., 2011). As the Tibetan thrust belt thickened and deformed Tethyan Himalaya rocks, Greater Himalaya rocks underneath the thrust belt were buried, metamorphosed, and transported southward via a shear zone, the Main Central thrust (LeFort, 1986; Vannay and Hodges, 1996; Catlos et al., 2001; Daniel et al., 2003). As this fault advanced above the brittle-ductile transition zone, it evolved into a more brittle fault, as found in the synformal klippen in Nepal (Fig. 1B). In far western Nepal, 40Ar/39Ar cooling ages on muscovite suggest that the Main Central thrust was active at ca. 25–21 Ma (Fig. 2; DeCelles et al., 2001; Robinson et al., 2006). Greater Himalaya rocks consist of a 5–20-km-thick assemblage of Late Proterozoic–early Paleozoic metamorphic rock (LeFort, 1975; Parrish and Hodges, 1996; DeCelles et al., 2000; Gehrels et al., 2003; Martin et al., 2005).

Emplacement of the Greater Himalaya rock on top of the cooler Lesser Himalaya rock along the Main Central thrust initiated burial and metamorphism of Lesser Himalaya strata (e.g., Hodges, 2000). In western and central Nepal, Lesser Himalaya rocks contain three superimposed rock sequences, the Lesser Himalayan sequence, the Gondwana sequence, and the Foreland Basin sequence. Detrital zircon data indicate that the Lesser Himalaya sequence was derived from the Indian craton in Paleoproterozoic to Mesoproterozoic time (e.g., DeCelles et al., 2000; Gehrels et al., 2011, and references therein). This study divides the Lesser Himalaya sequence into the Paleoproterozoic Lesser Himalaya rock, containing Kushma and Ranimata Formations and the ca. 1.73–1.86 Ga Proterozoic augen gneiss (dark gray, Fig. 1B), and the Mesoproterozoic Lesser Himalaya rock, containing the Sangram, Galyang, and Syangia Formations and the Lakharpata Group (light gray, Fig. 1B; Upreti, 1996, 1999; DeCelles et al., 2001; Robinson et al., 2006; Martin et al., 2011). In central Nepal, the Mesoproterozoic Lesser Himalaya units are overlain by the late Paleozoic and Mesozoic Gondwana sequence (Sakai, 1983, 1989). In far western Nepal, a thin unit similar to the Gondwana sequence is poorly exposed and intensely deformed (Robinson et al., 2006). The Foreland Basin sequence is composed of synorogenic sediment shed from the Himalayan thrust belt beginning in mid- to late Eocene time (Sakai, 1983) and consists of the Bhainskati Formation and the Dumri Formation (Sakai, 1983; DeCelles et al., 1998b). In western Nepal, the Dumri and Bhainskati Formations are found north of the Main Boundary thrust and are incorporated into the thrust belt (Robinson et al., 2006). The youngest synorogenic sediment is placed into a different tectonostratigraphic zone, the Subhimalaya, which is a 10–30-km-wide belt of imbricated Neogene Siwalik Group strata (Burbank et al., 1996) located exclusively south of the Main Boundary thrust at the toe of the thrust belt.

Major faults and structures south of the Main Central thrust in far western Nepal include the Ramgarh (Munsiari) thrust, Lesser Himalayan duplex, Dadeldhura thrust, Lesser Himalayan imbricate zone, Main Boundary thrust, and Main Frontal thrust (Figs. 1B, 2, and 3). The Ramgarh thrust is exposed in the footwall of the Main Central thrust, conformably underlies the Dadeldhura klippe in the north and south, and displaces Paleoproterozoic Lesser Himalaya rock >120 km to the south (Robinson et al., 2006). The Dadeldhura klippe (Fig. 1B), underlain by the Dadeldhura thrust, is part of the Lesser Himalaya “crystalline klippe.” These synformal klippen have Greater Himalaya affinity (e.g., Robinson et al., 2001) and are interpreted to have been emplaced via motion on the Main Central thrust (Johnson et al., 2001; Johnson, 2005; Pearson and DeCelles, 2005), which in turn was passively carried by the Ramgarh thrust. The Lesser Himalaya duplex is a hinterland-dipping antiformal duplex between the Main Central thrust and Dadeldhura thrust. These thrust sheets were once continuous but became erosionally separated via the development of the Lesser Himalaya duplex, which folded the overlying Greater Himalaya rock into broad antiformal (over the duplex) and synformal (in front of the duplex) structures (Fig. 3) (Robinson et al., 2001, 2003). Erosion through Greater Himalaya rocks stranded the “crystalline klippe” and cast uncertainty on the relationship between the Dadeldhura thrust and the Main Central thrust to the north. South of the Dadeldhura klippe, there is a sequence of complexly deformed Lesser Himalaya rock in the frontal part of the thrust belt called the Lesser Himalaya imbricate zone, of which the southernmost fault is the Main Boundary thrust, active as late as 4–5 Ma (DeCelles et al., 1998b) (Fig. 3). South of the Main Boundary thrust, synorogenic Siwalik Group strata are deformed by the Subhimalayan thrust system (Fig. 3). These Subhimalaya imbricated panels are separated from the Quaternary to modern sediments of the Indo-Gangetic plain by the active Main Frontal thrust (e.g., Lavé and Avouac, 2000). All faults are assumed to sole into the main décollement at the base of the thrust belt, the Main Himalayan thrust.

Pre–Foreland Basin Rocks

In late Paleozoic to Mesozoic time, the Gondwana sequence was deposited unconformably on top of the Proterozoic Lesser Himalaya sequence. Detrital zircon U-Pb ages from the Charchare conglomerate, a unit in the lower part of the Gondwana sequence, indicate a maximum depositional age of 125 ± 1 Ma (DeCelles et al., 2004) (Table 1). The overlying Amile Formation was deposited between Early Cretaceous to possibly early Eocene time (Sakai, 1983). Table 1 shows detrital zircon U-Pb age populations that indicate the sedimentary source was the Early Proterozoic to Archean rocks of the northern Indian craton located to the south (DeCelles et al., 2000, 2004).

Bhainskati Formation

Deposition of the Bhainskati Formation marks the beginning of the Foreland Basin sequence (Fig. 2) (DeCelles et al., 2004). In far western Nepal, the thickness of the Bhainskati Formation is estimated at ∼100 m. The suggested age for the lower portion of the formation is late Paleocene–early Eocene (Fuchs and Frank, 1970), and the upper portion has a mid–late Eocene age based on foraminifera and marine fossils (Sakai, 1983). Najman et al. (2005) records a 45 Ma zircon population in the upper half of the Bhainskati Formation in central Nepal at Dumri Bridge (Fig. 1B), requiring that the upper half be younger than 45 Ma.

Table 1 shows the detrital zircon U-Pb age populations present in the Bhainskati Formation. These age groups are found in both Tethyan Himalaya and Greater Himalaya rocks and indicate a Himalayan thrust belt source. The average ɛNd(T) value is −10 for the Bhainskati Formation (DeCelles et al., 2004), which is less negative than the average ɛNd(T) value of Greater Himalaya rocks (−16; Robinson et al., 2001), suggesting either a Tethyan Himalaya source, or an additional component of less-evolved detritus. In the correlative Subathu Formation (NW India), the lower half of the formation is more enriched in Ni and Cr, and is interpreted as originating in the mafic arc and ophiolitic rocks of the Indus suture (Najman and Garzanti, 2000; Najman et al., 2005). The upper half of the formation has lower concentrations of Ni and Cr, which may indicate a shift in the source material to sedimentary rocks of the Tethyan Himalaya (Najman and Garzanti, 2000; Najman et al., 2005).

The base of the formation has phosphatic and organic-rich mud rock and limestone that transitions up stratigraphic section into oolitic ironstone, suggesting that paleobathymetry gradually shoaled until the depositional site lay above wave base but far enough offshore to prevent the influx of siliciclastic detritus (DeCelles et al., 1998b). On top, a well-developed unconformity, a paleosol, spans late Eocene–Oligocene time, and has been interpreted to represent subaerial exposure until ca. 20 Ma (DeCelles et al., 1998b). This paleosol is present at Dumri Bridge but absent at Swat Khola (Fig. 1B), which led DeCelles et al. (1998b) to speculate that the paleosol was erosionally removed before deposition of the Dumri Formation. Using the definitions of DeCelles and Giles (1996), the Bhainskati Formation is interpreted as a back-bulge depozone, with the condensed section indicating the southward migration of the forebulge, and deposition of the overlying Dumri Formation marking the beginning of distal foredeep deposition (DeCelles et al., 1998b).

Dumri Formation

At Dumri Bridge (Fig. 1B) in south-central Nepal, the contact is sharp between the Bhainskati paleosol and the red rocks of the 800–1200-km-thick Dumri Formation (DeCelles et al., 1998b). Here, the youngest zircon population shows that the base is younger than 30–32 Ma (Najman et al., 2005). To the north, near the town of Chainpur (Fig. 1B), the lower half of the formation is a green, strongly indurated, micaceous sandstone. Approximately 30 m up section, a pronounced bauxite paleosol occurs with kaolinite and hematite pisolites, above which the formation returns to a micaceous sandstone. The upper half of the formation is dominated by red and green shale. Detrital micas provide a 40Ar/39Ar age from the lower micaceous sandy section of ca. 20 Ma (SR9), the maximum age of the Dumri Formation in this location (DeCelles et al., 2001). Additional Dumri Formation conglomerate crops out north of this sample, suggesting that there, the formation could be older. In the south, near the Main Boundary thrust, Ojha et al. (2009) suggested an age at Swat Khola of 19.9–15.1 Ma using magnetostratigraphy.

Table 1 shows the detrital zircon U-Pb age populations, which are similar to that of the Bhainskati Formation, except more grains in the 470–500 Ma and ca. 1000–1500 Ma age ranges (DeCelles et al., 2004). Because the detrital zircon age spectra for the Greater Himalaya and Tethyan Himalaya rocks are so similar (e.g., Amidon et al., 2005), they cannot be used to identify when Greater Himalaya rocks began to supply detritus. The average ɛNd(T) value is −14.4 (Robinson et al., 2001), indicating that the formation is more isotopically evolved than Bhainskati Formation. This value is similar to both Tethyan Himalaya and Greater Himalaya average values (Robinson et al., 2001).

Correlative strata in NW India show low Cr and Ni, compatible with a Tethyan Himalaya or Greater Himalaya source (Najman and Garzanti, 2000; Najman et al., 2005). Petrographic studies reveal abundant low-grade metasedimentary lithic grains, indicating the first evidence of eroded metamorphic rocks. DeCelles et al. (1998b) documented significant amounts of plagioclase, which they interpreted as originating from Greater Himalaya affinity rocks from the Dadeldhura klippe (Fig. 1B). To the east of the Dadeldhura klippe, in west-central Nepal (Fig. 1B), Sakai et al. (1999) noted the lack of high-grade metamorphic lithic fragments, and suggested that the phyllitic lithic fragments were supplied from lower Lesser Himalaya and/or Tethyan Himalaya rock. Sakai et al. (1999) also reported a 40Ar/39Ar muscovite cooling age of 16–17 Ma, which they interpreted as cooling from emplacement of the Main Central thrust over these foreland basin rocks.

Siwalik Group

The Siwalik Group is a middle Miocene to Pliocene belt of synorogenic foreland basin sediment at the front of the thrust belt informally divided into upper, middle, and lower units (Fig. 2; Quade et al., 1995). This division is predominantly lithostratigraphic, and so it is useful for broad regional comparisons but cannot be used for detailed chronostratigraphy (DeCelles et al., 1998a). The cross section from which the reconstruction is derived for this study is located along the straight black line in Figure 1B. In the Subhimalaya, this line falls in between two measured sections though the Siwalik Group at Khutia Khola and Karnali River (thick gray lines, Fig. 1B).

Lower Siwalik Unit

At Khutia Khola, deposition of the lower Siwalik unit began by ca. 13.3 Ma (Ojha et al., 2000); however, the base is covered by alluvium, so the unit may be as old as 15.1 Ma, the youngest Dumri Formation rock. The thickness is 862 m of interbedded sandstone and shale (DeCelles et al., 1998a). At the Karnali River, deposition began at 16 Ma, and the thickness of the unit is ∼2000 m (Gautam and Fujiwara, 2000). Detrital zircon U-Pb ages (Table 1) show the same zircon ages as found in the Dumri Formation, in addition to an increase in grains from the ∼1000–1500 Ma and 470–500 Ma populations (DeCelles et al., 2004). The average ɛNd(T) value is −16, the same as the average for Greater Himalaya and Tethyan Himalaya rocks (Huyghe et al., 2001; Robinson et al., 2001). The first appearance of high-grade metasedimentary lithic grains, kyanite and sillimanite, occurs in the upper part of the lower Siwalik unit at Khutia Khola (DeCelles et al., 1998b), which is dated at ca. 11 Ma (Ojha et al., 2000). Thus, Greater Himalaya rocks were at the surface contributing sediment at this time. At 10–11 Ma at the Karnali River, Huyghe et al. (2001) noted a coarsening grain size at the microscopic scale, and Nakayama and Ulak (1999) noted a similar coarsening at the macroscopic scale.

Middle Siwalik Unit

As noted already, the reconstruction we present here is in between the physical locations of Khutia Khola, where the lithologic change from lower to middle Siwalik units is at 10.8 Ma (Ojha et al., 2000), and the Karnali River, where the lower to middle transition is at 9.7 Ma (Gautam and Fujiwara, 2000). Between 11 and 9.7 Ma, a distinct normal polarity chron (5n) can be identified on both paleomagnetic data sets and tied to the geomagnetic polarity time scale (Cande and Kent, 1995). In Khutia Khola, this chron is within the middle Siwalik unit; however, in the Karnali River section, this chron is within the lower Siwalik unit. Because the change in lithology that defines the boundary between the lower and middle Siwalik units occurs at different times in these two sections, we tie our observations to both lithology and age. Table 1 combines both sections and groups data older than 10 Ma in the lower Siwalik unit and younger than 10 Ma in the middle Siwalik unit.

The thickness of the middle Siwalik unit at Khutia Khola is 2468 m and is dominated by stacked channel sandstones punctuated by thin shale beds (DeCelles et al., 1998a). Lesser Himalaya–derived carbonate clasts first appear at the base of the section (at 10.8 Ma), and the magnitude increases up section. An upward excursion in the 87Sr/86Sr ratio of pedogenic carbonate nodules in mudstone at ca. 9 Ma is another indication of Lesser Himalaya carbonate erosion (Quade et al., 1997). The carbonate clasts are derived from the Mesoproterozoic Lakharparta Group and not the Paleoproterozoic lower Lesser Himalaya rocks, indicating derivation from rocks in the Lesser Himalaya duplex. In addition, ɛNd(T) values become more negative, from −16 to −18 between 10 and 9 Ma, and remain at these values throughout middle Siwalik unit deposition (Robinson et al., 2001). The transition between the lower to middle units in the Karnali River section is marked by prevalence of stacked fluvial channel sandstones. The section is ∼1500 m thick; detrital staurolite first occurs at 9–10 Ma, and detrital kyanite first occurs at 8 Ma (Szulc et al., 2006). The ɛNd(T) values from the Karnali River section become markedly negative (−18) at 10 Ma, indicating influx of Lesser Himalaya material (Huyghe et al., 2001; Szulc et al., 2006). Most U-Pb detrital zircon studies do not break out the lower, middle, and upper Siwalik units, and instead discuss trends. In the middle and upper Siwalik units, the proportion of Lesser Himalaya rock having detrital zircon with U-Pb ages older than 2 Ga increases up stratigraphic section (DeCelles et al., 1998b), indicating the unroofing of Lesser Himalaya rocks.

Upper Siwalik Unit

The middle to upper Siwalik unit transition at Khutia Khola is estimated at younger than 4.5 Ma (Ojha et al., 2000). The thickness at Khutia Khola is at least 1000 m thick and is dominated by fluvial conglomerate punctuated by fluvial sandstone (DeCelles et al., 1998a). This unit was not measured at Karnali River (Gautam and Fujiwara, 2000), but was determined to be younger than 5 Ma. The trends of the detrital zircon age populations and ɛNd(T) values are similar to those of the middle Siwalik unit (Table 1). After 4–5 Ma in Khutia Khola, a proximal facies is indicated by the appearance of Lesser Himalaya conglomerate clasts (DeCelles et al., 1998b), suggesting a nearby surface exposure of Lesser Himalaya rocks.

Far western Nepal is remote and difficult to access. As such, existing data sets are sparse (Table 2). DeCelles et al. (2001) reported a 40Ar/39Ar muscovite cooling age of ca. 21 Ma in the Budhiganga gneiss (SR40; DeCelles et al., 2001) from the northern limb of Dadeldhura thrust sheet in the Dadeldhura klippe. Robinson et al. (2006) reported a 40Ar/39Ar muscovite cooling age of ca. 25 Ma in the hanging wall of the Main Central thrust sheet (SR124). These two ages are in Greater Himalaya rocks (Table 2) and are located along the reconstruction line. East of the Karnali River (Fig. 1B), Sakai et al. (1999) reported an Ar plateau age of 25.69 ± 0.13 Ma for biotite in gneiss on the southern limb of the Dadeldhura klippe in Greater Himalaya rocks. These ages indicate that in far western Nepal, Greater Himalaya rocks cooled through the blocking temperature for muscovite and biotite from 25 to 21 Ma, and these are interpreted to represent cooling due to erosional exhumation that accompanied initial motion along the Main Central thrust (Fig. 2).

In far western Nepal, the northernmost exposure of Lesser Himalaya rocks in the Lesser Himalaya duplex indicates peak temperatures of 500–550 °C, identified using Raman spectroscopy of carbonaceous material (Bollinger et al., 2004). Beyssac et al. (2004) reported that the Raman spectroscopy of carbonaceous material temperatures in Lesser Himalaya rock along the Nepal-India border decrease gradually from 540 °C in the footwall of the Main Central thrust to 330 °C in the middle of the Lesser Himalaya duplex. This southward decrease in temperature away from the Main Central thrust indicates that the magnitude of overburden that buried these samples also decreases to the south. In addition, the temperatures indicate that the northernmost Lesser Himalaya rocks reached the temperatures necessary to reset 40Ar/39Ar muscovite cooling ages. Robinson et al. (2006) reported a 40Ar/39Ar muscovite cooling age from the Ramgarh thrust sheet (SR123) of 17–7 Ma in far western Nepal (Table 2). A 40Ar/39Ar muscovite cooling and isochron age of 12.5–12 Ma (SR103a; Robinson et al., 2006) was also reported from the northern part of the Lesser Himalaya duplex (Table 2).

Due to the remote location, no paleobarometric data exist for far western Nepal. Thus, the only data available are via extrapolation from central Nepal, where the Ramgarh thrust experienced peak pressures of 10 kbar at a temperature of ∼550 °C, and the Main Central thrust experienced peak pressures of 12 kbar at a temperature of ∼725 °C (Kohn, 2008), corresponding to 35 and 42 km, assuming a buoyant crust. However, the thickness of the Greater Himalaya rocks in central Nepal (16 km) is doubled by the Langtang thrust. The distance from the Main Central thrust to the Langtang thrust is 8 km (Kohn, 2008). Along the reconstruction line in this study, Greater Himalaya thickness is 8 km and restores to a much shallower ∼25 km depth, with the approximate location shown in Figure 4A.

We link the exhumation documented in the hinterland to the depositional history of the foreland basin using a sequentially restored balanced cross section. We model the sequential deformation and associated erosion and deposition, and use the published data to validate the model, which predicts the magnitude of exhumation, depositional basin extent, and the stratigraphic units that were exposed at the surface. This approach allows us to use published ages of exhumation and deposition to place age constraints on fault motion, but it requires several assumptions that are detailed in the following section.


Geometry of the Original Indian Margin

The initial configuration of the units before 25 Ma is portrayed as flat-lying, undeformed strata on the northern Indian margin. This configuration contains simplifications and assumptions.

(1) Pre-Himalayan deformation of Proterozoic rocks. The Lesser Himalaya rocks in far western Nepal are predominantly Proterozoic in age and may have experienced several periods of deformation. Intrusive bodies 1.73–1.86 Ga in age indicate a Himalayan-wide arc event (for more details, see Kohn et al., 2010). A ca. 500 Ma event has been proposed within Greater Himalaya rocks that may have also affected Lesser Himalaya rocks along the northern Indian margin (DeCelles et al., 2000; Gehrels et al., 2003, 2006; Cawood et al., 2007). However, with these caveats, extensive field mapping of Lesser Himalaya rocks throughout far western Nepal indicates that foliation is parallel to original sedimentary bedding, with sedimentary structures emphasizing that bedding has not been overturned (Robinson and Pearson, 2006). In addition, where preserved, the younger Foreland Basin sequence is parallel to underlying Proterozoic formations (Robinson et al., 2006), highlighting that at least in this part of the Himalaya, the Lesser Himalaya section was not significantly deformed prior to the Himalayan event.

(2) Pre-Himalayan and Himalayan internal deformation of Greater Himalaya rocks. Greater Himalaya rocks are metamorphic rocks with internal folding and faulting, and treating the entire unit as a coherent slab is a gross simplification. To avoid complications with differentiating pre-Himalayan and Himalayan deformation, as well as uncertainties in the magnitude of internal strain within the Greater Himalaya, we treat displacement on the Main Central thrust as just the magnitude needed to account for overlap of Greater Himalaya rocks on Lesser Himalaya rocks. The reconstruction assumes that Greater Himalaya rock north of the Main Central thrust is the same rock contained in the Dadeldhura klippe; thus, the Dadeldhura thrust is also the Main Central thrust (Fig. 3) (Johnson et al., 2001; Johnson, 2005; Pearson and DeCelles, 2005). This assumption places a minimum value on shortening and is the simplest solution. We restore Greater Himalaya rocks at the northern limit of restored Lesser Himalaya rocks (e.g., Robinson et al., 2006), as the original relationship between Lesser Himalaya and Greater Himalaya rocks is uncertain and not preserved. Greater Himalaya rocks have been interpreted as a terrane that accreted onto India during the ∼500 m.y. orogenic event (e.g., DeCelles et al., 2000) or as part of a passive margin sequence that was only deformed during the Cenozoic collision (e.g., Myrow et al., 2003, 2010). For the purpose of this reconstruction, the origin of Greater Himalaya rocks is not a critical factor, but rather the location in the crust from which they were translated. The maximum southward extent of the Greater Himalaya rocks must be north of restored Lesser Himalaya rocks on the Greater Indian margin.

(3) Pre-Himalayan geometry of Tethyan Himalaya rocks. The sparsely exposed Gondwana sequence through western Nepal is assumed to thicken to the north and become part of the Cambrian–Ordovician through Cretaceous Tethyan Himalaya rocks, which overlay Greater Himalaya rocks precollision. In far western Nepal, the age of the Tethyan Himalaya rocks north of the South Tibetan Detachment system (STDS) is not clear. To the east of this reconstruction in the Simikot region, Murphy and Copeland (2005) identified basal Tethyan Himalaya rocks as Cambrian–Ordovician in age. Another unknown factor is the full thickness of the Tethyan Himalaya stratigraphic section. Available data suggest a thickness of 17–20 km (Ratschbacher et al., 1994; Zhang and Guo, 2007) north of the STDS. Thus, in our reconstruction, Tethyan Himalaya rocks are assumed to reach thicknesses of 17 km at their northern extent, and directly overlie Greater Himalaya rocks. For our original configuration, we show a gently north-dipping (1.5°) Indian passive margin, where Cretaceous rocks are at approximate sea level from the modern foreland to the northernmost extent of the restoration.

Early Himalayan Deformation

Most researchers agree that the Himalayan collision began at ca. 55 Ma (see review in Najman, 2006), with many studies preferring an age of 52–50 Ma (e.g., Zhu et al., 2005; Green et al., 2008; Dupont-Nivet et al., 2010; Najman et al., 2010; Wang et al., 2011). Prior to 25 Ma, an entire thrust belt was built from Tethyan Himalaya strata, the Tibetan thrust belt, as the Indian margin shortened and thickened from ca. 55–25 Ma (Ratschbacher et al., 1994; Murphy and Yin, 2003; Leech et al., 2005; Ding et al., 2005; Webb et al., 2011). Presumably, southward propagation of these deformed Tethyan Himalaya rocks facilitated the burial of Greater Himalaya rocks to their peak temperature and pressure conditions (Patel et al., 1993; Robinson et al., 2003; Murphy and Yin, 2003; Ding et al., 2005; Aikman et al., 2008; Webb et al., 2011). However, because of the remote nature of the Tethyan Himalaya rocks in far western Nepal and Tibet, studies of the Tibetan fold-and-thrust belt are lacking. Thus, the burial depth of Greater Himalaya rocks and the magnitude of shortening accommodated within the Tethyan Himalaya system are unknown.


Estimating topography through time is fraught with uncertainties; thus, we use a very simple approach to estimating the evolution of topography. To be able to estimate both the topographic load as well as the amount of material removed though erosion, we assume a constant topographic slope from 0 to 4.5 km over a 150 km width, modeled after the modern topographic slope of the Himalaya (e.g., Bookhagen and Burbank, 2006). The transition from no topography to a simple topographic gradient begins at the frontal thrust for a specific time period. The assumption that the Tibetan Plateau had reached its modern height in the part of Tibet north of far western Nepal by 25 Ma is supported by the Nima basin in central Tibet, which indicates an arid climate and high paleoelevations of 4.5–5 km by 26 Ma (DeCelles et al., 2007), as well as negligible exhumation of the region from ca. 45 Ma to present (Hetzel et al., 2011; Rohrmann et al., 2012). An early Himalayan fold-and-thrust belt with a northward-increasing critically tapered slope overlying the Greater Himalaya rocks at the start of the deformation is supported by prograde monazite ages within Greater Himalaya rocks, which indicate burial from 30 to 20 Ma in the Langtang region of central Nepal (Kohn et al., 2004). To test the effect of our topographic assumptions on both the flexural evolution of the Himalayan foreland basin and estimates of erosion, we produced a series of models that are discussed with our simplified topography model in the methods section.

Model Parameters

We bracketed the time steps based on the available synorogenic and hinterland data. Although the data are sparse, seven specific well-dated windows of time bracket magnitudes of deformation documented in sequential reconstructions (Table 3).

A paleosol developed on the exposed Bhainskati Formation at the southern limit of the foreland basin (an exposed forebulge; Fig. 4B) from 40 Ma to ca. 20 Ma (DeCelles et al., 1998b) sets the initial width of the foreland basin after motion on the Main Central thrust ca. 25–21 Ma. These data are from Dumri Bridge in south-central Nepal; although the paleosol is not present in Swat Khola, we use the description at Dumri Bridge because these observations are along strike with the Swat Khola Dumri Formation section. The restored distance between the front of the Main Central thrust after 125 km of motion and the documented paleosol in the Bhainskati Formation is ∼260 km. One major factor that controls the width of the foreland basin is the flexural rigidity of the Indian plate (e.g., Jordan and Watts, 2005). Thus, we would predict a match between the estimated width of the foreland basin based on geological constraints, from both foreland basin sediments (thickness, sedimentology, and age) and structure (restored distance between the thrust front and foreland basin deposits), and the estimated width of the basin from flexural modeling of the thrust load. Flexural rigidity and effective elastic thickness (Te) are two parameters that provide a measure of the long-term strength of the lithosphere (Jordan and Watts, 2005). Te is related to the flexural rigidity through estimates for Young’s modulus and Poisson’s ratio (Turcotte and Schubert, 1982).

The reconstruction program, 2D Move, allows the user to estimate the flexural isostatic deflection of the lithosphere based on the load generated in a sequential reconstruction. The parameters used in the flexural modeling include: effective elastic thickness (Te), mantle density, load density, and Young’s modulus. Although any of these parameters may be changed in successive time steps, our goal was to choose parameters supported by published gravity models (Lyon-Caen and Molnar, 1985; Jordan and Watts, 2005) and find a self-consistent model that accurately reproduces the basin history. Throughout the model we present here, to model the Eocene–early Miocene flexural rigidity, we used a Te of 70 km (Jordan and Watts, 2005), which is approximately equal to a flexural rigidity (D) of 2 × 1024 and is between the values 7 × 1023 and 3 × 1024 used by DeCelles et al. (1998b). Mantle density was set at 3300 kg/m3, while the load density was estimated at 2700 kg/m3. Young’s modulus was set at 0.7 × 105 MPa.

To create the time steps shown in Figure 4, the displacement along the faults was moved sequentially, with both the overlying strata and topography deformed from motion on underlying faults. The deformed topographic surface (area above sea level, gray dot-dash lines on Fig. 4) was modeled as the top of the load, and the base of the load was sea level. For each step, we calculated the isostatic difference from our initial condition with the top of the Cretaceous at sea level. After the program calculated flexural compensation of the load, a new topographic profile was drawn starting at sea level from the southernmost extent of the load and reaching 4.5 km, 150 km to the north (black dashed lines on Fig. 4). Any material above this topographic surface was assumed to have been removed by erosion (area not filled in with a pattern on Fig. 4). Rock that passed through this erosion surface is not included as part of the load put on the plate in the subsequent steps. In each time step, this scenario was repeated, with the new topographic surface being deformed by the southward-propagating deformation, the deformed topographic surface acting as the load, isostatic compensation, and a new topographic surface being generated, above which material was removed by erosion. The eroded material is the area between the deformed topographic surface (the load) and the new topographic slope. The flexure of the plate produced an accommodation space found at the front of the thrust belt into which the synorogenic sediments were shed. Any region below sea level was assumed to be collecting synorogenic sediment. We then compared the thickness of synorogenic sediment at locations where geologic constraints from the foreland basin exist. The critical components that we needed the model to match were the location and depth of foreland basin data, the dip of the décollement, and the exposure of rocks mapped at the surface today. Detailed data from the Dumri Formation (ca. 20–15 Ma) and the Siwalik Group sediments (15 Ma to present), presented previously and summarized in Table 1, were used as constraints for checking both the southern extent of the load and the flexural isostatic parameters as the model developed through time. The final model needed to have a décollement dip consistent with that identified for the Main Himalayan thrust (Ni and Barazangi, 1984; Pandey et al., 1999).

To test the effect of our simplified topographic assumption, we repeated the method steps varying topography, Te, and the density of the load to find the solutions that matched the location and depth of foreland basin data, the dip of the décollement, and the exposure of rocks mapped at the surface today. Because increasing topography and/or density increases the load while an increase in Te supports that load, there are several variations of topography (1°–3° uniform slopes as well as nonuniform slopes), density (2600–2800 kg/m3), and values of Te (45–90) that matched the available constraints. Thus, the absolute value of these variables is not critical or unique for the model. However, within these variations, every scenario that matched our constraints of location and depth of foreland basin data, the dip of the décollement, and the exposure of rocks mapped at the surface today produced a similar history for both the magnitude of exhumation and depth and migration of the foreland basin. Using a best-fit “no topography” model decreased foreland basin thickness (by 1.5–2 km), shallowed the position of the décollement by 1°–1.5°, and increased the magnitude of material removed (2 km). Thus, the rock that was at the surface did not match the data from the hinterland and foreland and highlighted that there are limits to permissible variations in topography, Te, and the density. These sensitivity tests illustrated that while our models are nonunique, the first-order patterns of the location and timing of exhumation and deposition of sediments are controlled by the evolving structural load (Fig. 4). The similarity in the exhumation estimates argues that the largest control on location, time, and amount of exhumation is the vertical component of deformation (Whipp et al., 2007).

We used thermochronologic data from the hinterland to compare the magnitude of exhumation proposed in our model to that measured with available thermochronometers. The closure temperature for any mineral system depends on the exhumation rate (e.g., Brandon et al., 1998), while the depth of that closure temperature is a function of both exhumation rate and the original geothermal gradient. The thermochronologic data for far western Nepal are too meager to warrant detailed thermodynamic modeling; thus, to compare the data to our model, we obtained an average closure temperature from the literature and then used the predicted depth of the sample to calculate the necessary geothermal gradient required to match the data to the model. This approach is not intended to solve for the gradient or make a robust calculation of exhumation, but simply to provide enough information to evaluate whether the match is viable.

The initial configuration of the sequential model is illustrated in Figure 4A. Due to uncertainties in the geometry, thickness, and extent of Tethyan Himalaya shortening, we only show the undeformed Tethyan Himalaya section. As a result of an undeformed Tethyan Himalaya section, as well as uncertainty in the initial pressure and temperature conditions in far western Nepal, the base of Greater Himalaya rocks is initially at 25 km (17 km Tethyan Himalaya + 8 km Greater Himalaya). Initial topography is displayed with elevation gradually increasing toward the north to 4.5 km. The following time steps were chosen based on the criteria outlined in Table 3.

From 25 to 20 Ma (Fig. 4B) in the foreland basin, the Bhainskati Formation paleosol continued to form at the front of the thrust belt. In Figure 4B, the arrow representing the paleosol points to the top of the Lesser Himalaya sequence because the Bhainskati Formation is only 100 m thick and too thin to be represented on the reconstruction. The Main Central thrust carried Greater Himalaya rock 125 km over the Lesser Himalaya in the Ramgarh thrust sheet. During this time, the northern Lesser Himalaya duplex and Ramgarh thrust sheet were buried, and Lesser Himalaya rocks were metamorphosed (dark-gray color in Fig. 4B). Post erosion, sample SR124, immediately adjacent to the Main Central thrust (25 Ma cooling age), and sample SR40, within the Dadeldhura klippe (21 Ma cooling age), are 18 km and 15 km below the surface, respectively. This suggests a geothermal gradient of 22–30 °C/km to reconcile cooling ages with the erosional history implied by the reconstruction, assuming a 400–450 °C closure temperature for muscovite (Hames and Bowring, 1994). The foreland basin extends southward 90 km (Fig. 4B), ∼25 km north of the Bhainskati Formation paleosol described by DeCelles et al. (1998b) at Dumri Bridge. A sample from the Dumri Formation near the village of Chainpur, described earlier, has a detrital mica age of ca. 20 Ma; however, the exact location within the section is unknown. Thus, north of this sample, the Dumri Formation is likely older, 25–20 Ma. The oldest available age on the Dumri Formation is <30 Ma from samples in NW India (Najman et al., 2005).

From 19.99 to 16 Ma (Fig. 4C) in the foreland basin, deposition of the subaerial, fluvial Dumri Formation began at 19.9 Ma at Swat Khola (Ojha et al., 2000, 2009). The Ramgarh thrust was active in several thrust sheets, totaling ∼109 km (Robinson et al., 2006). Post motion on the Ramgarh thrust and after erosion, sample SR123 (17–7 Ma cooling age) restores to a depth of 17 km, but remains between 17 and 13 km depth until ca. 12 Ma. The length of time at this depth (∼7 m.y.) may explain the complicated age spectra. Using a closure temperature range of 400–450 °C, a closure depth of 17–13 km would suggest geothermal gradients of 25–35 °C/km. Alternatively, the muscovite could represent the age that muscovite grew during dynamic recrystallization while the Ramgarh thrust was active, possibly at temperatures as low as <300 °C (Passchier and Trouw, 1996), as suggested by Robinson et al. (2006). However, maximum temperatures for these rocks measured from Raman spectroscopy of carbonaceous material range from 400 °C to 500 °C (Bollinger et al., 2004; Beyssac et al., 2004). Further support for 400–500 °C maximum temperatures comes from the cooling history of sample SR103a, immediately north of sample SR123. Sample SR103a is from the first Lesser Himalaya thrust sheet under the Main Central thrust, a Ramgarh thrust imbricate, and it cooled rapidly from 12.5 to 12 Ma (Fig. 4E), precluding significantly lower temperatures for sample SR123. The foreland basin created by displacement on the Ramgarh thrust extends to 26 km, and suggests ∼0.5 km of Dumri deposition by 16 Ma at Swat Khola. At the same time as thrust faulting was occurring in the thrust belt, normal faulting on the STD occurred in the Tethyan Himalaya. Detrital and geochronological data suggest slip on the STD faults from 20 Ma to 9 Ma (Murphy and Copeland, 2005) in the Simikot region, NW of the village of Chainpur. Because of the uncertainty surrounding the timing and geometry of the STD through time, the fault system is not added to Figure 4. In Figures 4F–4H, the STD at the surface would be at the contact between Greater Himalaya and Tethyan Himalaya rocks.

From 15.99 to 13 Ma (Fig. 4D) in Swat Khola, deposition of the Dumri Formation ended at ca. 15 Ma, and deposition of the lower Siwalik unit occurred from ca. 15 Ma to 10.8 Ma (Ojha et al., 2000, 2009). North of the Dadeldhura klippe, where the Dumri Formation is sandstone, the Ramgarh thrust sheet cuts the Dumri Formation at Arnakoli Khola. Based on the minimum 15 Ma age of the Dumri Formation at Swat Khola, DeCelles et al. (2001) argued that the age of slip on the Ramgarh thrust must be younger than 15 Ma, because that fault cuts the Dumri Formation. Arnakoli Khola is ∼130 km north of Swat Khola (see annotation on Fig. 4D). Our model suggests that initial deposition of the Dumri Formation at Arnakoli Khola was as early as 25 Ma, but it becomes cut by the Ramgarh thrust between 16 and 13 Ma. In Khutia Khola, deposition of the lower Siwalik unit began by 13.3 Ma (Ojha et al., 2000); however, in the Karnali River, Gautam and Fujiwara (2000) documented deposition initiating at 16 Ma. Our reconstruction shows deposition of the lower Siwalik unit from 15.99 to 13 Ma in the location of the frontal two thrusts of the Subhimalayan thrust system. Over this time window, the northern Lesser Himalaya duplex continued to build via the emplacement of lower Lesser Himalaya thrust sheets.

From 12.99 to 11 Ma (Fig. 4E), the Lesser Himalaya duplex continued to build by incorporating more of the Paleoproterozoic Lesser Himalaya sequence. The increase in the thickness of material incorporated into the Lesser Himalaya duplex creates a more pronounced structural high and facilitates erosion to expose the highest-grade Greater Himalaya rocks in the hanging wall of the Main Central thrust at the surface. Proterozoic Lesser Himalaya rocks rapidly exhume from ∼12 km depth to ∼5 km. At the transition between the time steps shown in Figures 4D and 4E, Lesser Himalaya sample SR103a records rapid cooling at ca. 12 Ma. At Khutia Khola, kyanite and sillimanite metasedimentary lithic fragments first appear in the upper portion of the lower Siwalik unit at ca. 11 Ma (DeCelles et al., 1998a). This influx of high-grade metamorphic minerals corresponds with the emplacement of thicker Lesser Himalaya thrust sheets in the Lesser Himalaya duplex.

From 10.99 to 10 Ma (Fig. 4F), the transition from the lower to middle Siwalik unit occurred in Khutia Khola at 10.8 Ma (Ojha et al., 2000). In the hinterland, the first thrust sheet containing the entire Mesoproterozoic Lesser Himalaya sequence is emplaced in the Lesser Himalaya duplex, creating a 7 km ramp in the Main Himalayan thrust. As illustrated in Figure 4F, this structural high has a pronounced effect on the magnitude of erosion immediately above it. This focused erosion exposes Lesser Himalaya rock, creating a new sediment source for the foreland basin. At this time, input from Lesser Himalaya rocks is documented by the change to a more negative ɛNd value (Huyghe et al., 2001; Robinson et al., 2001; Szulc et al., 2006) and the first appearance of Lesser Himalaya–derived carbonate clasts (DeCelles et al., 1998a).

From 9.99 to 5 Ma (Fig. 4G), the middle Siwalik unit continued to be deposited at both Khutia Khola and Karnali River. In the thrust belt, the last three Lesser Himalaya thrust sheets of the Lesser Himalaya duplex were consecutively emplaced, and the locus of thrusting was transferred south of the Dadeldhura klippe, forming the Main Boundary thrust. Emplacement of the Main Boundary thrust, the last Lesser Himalaya thrust sheet, over Subhimalaya rock caused the southern limb of the Dadeldhura klippe to uplift and rotate to the north, giving the Dadeldhura klippe its synclinal shape. DeCelles et al. (1998a) argued for emplacement of the Main Boundary thrust by 4–5 Ma. Map patterns indicate that the Main Boundary thrust cuts middle Siwalik unit rock (Robinson et al., 2006); thus, some component of motion on the Main Boundary thrust occurred after middle Siwalik unit deposition (ca. 5 Ma). Although no cooling ages exist for the Lesser Himalaya duplex in far western Nepal, in central Nepal and NW India, pre-Cenozoic 40Ar/39Ar cooling ages are prevalent in the Lesser Himalaya duplex and in rocks carried by later thrusts to the south (Catlos et al., 2001; Wobus et al., 2003; Bollinger et al., 2004; Kohn et al., 2004; Célérier et al., 2009), suggesting that samples were not buried deeply enough during the Cenozoic orogeny to completely reset muscovite ages.

From 4.99 to 0 Ma (Fig. 4H) in Khutia Khola, the upper Siwalik unit was deposited at a projected time of 4.5 Ma (Ojha et al., 2000). The upper Siwalik unit is a proximal facies that contains large Lesser Himalaya conglomerate clasts, suggestive of an immediate source to the north (DeCelles et al., 1998b). South of the Main Boundary thrust, the Subhimalayan thrust system forms exclusively within Subhimalaya rocks. Slip on the Main Frontal thrust may have begun at ca. 2 Ma and continues in Karnali River (van der Beek et al., 2006). After the end of the major deformation sequence in the thrust belt, eight small-scale thrust and normal faults with less than 5 km of displacement were added to the reconstruction program based on field relationships (see Robinson, 2008). Because the sequence and timing are difficult to determine and the offsets are less than 5 km, only the end result is shown.

The gray dot-dash lines on Figure 4 represent the deformed topographic surface. The rock removed by erosion in each time step is estimated as the difference between the deformed topographic surface and the new topographic profile drawn after flexural compensation of the load (black dashed lines on Fig. 4) (see area, Table 4). To produce a denudation rate estimate, both the duration of exhumation and the length scale over which that denudation was applied are needed (Barnes and Pelletier, 2006; McQuarrie et al., 2008). Duration is the length of each frame in Figure 4. The length-scale is the distance along the cross section where rock was “removed” via erosion in the model. The thickness of material removed (exhumation amount) is calculated by dividing the area by the length. The exhumation rate is calculated by dividing the exhumation amount by the duration. Exhumation is a result of either erosion (rivers, glaciers, or landslides) or normal faulting. Only minor normal faulting has affected the thrust belt in far western Nepal (Robinson et al., 2006; Robinson, 2008). Figure 4H illustrates two normal faults that were added in at the end of the reconstruction and that match faults observed in the field. These faults have less than 5 km of offset and are located between the Main Boundary thrust and southern Dadeldhura thrust in the Lesser Himalaya imbricate zone. The ages of these faults are unknown. Because the offsets are small and do not appear in consecutive cross sections along strike (Robinson et al., 2006), we attribute the rate of removal to erosion.

The rate of erosion we calculate varies from 0.7 to 4.2 mm/yr, with most estimates between 0.7 and 1.4 mm/yr (Table 4), comparable to exhumation estimates across the Himalaya and particularly close in magnitude to estimates obtained in western Nepal (Table 5). Figure 4F, from 10.99 to 10 Ma, has an anomalously high rate of erosion (4.2 mm/yr) when compared to periods of time before and after, and it also has the shortest duration. This erosion rate increase occurs when ∼20 km of overburdened is thrust over a 7 km ramp in the Main Himalayan thrust and exposes Lesser Himalaya rock at the surface. The time window over which that happens is not exact; however, sedimentological constraints from Khutia Khola require a Lesser Himalaya source at 10.8 Ma. Thus, at least half of the erosion must be before 10 Ma. By extending the time window 1 m.y., from 10.99 to 9 Ma (see asterisk in Table 4), the exhumation rate drops to 2.1 mm/yr, and the erosion rate in the subsequent time step increases to 1.1 mm/yr. In general, the rate of erosion is fairly constant through the emplacement of major thrust sheets, and more minor faults in far western Nepal. The notable increase in erosion is linked to the creation and southward propagation of 7-km-thick ramps that built the southern half of the Lesser Himalaya duplex.

We can compare our estimates of magnitude and rates of erosion to published estimates through the region (Tables 4 and 5). Bernet et al. (2006) used detrital zircon fission track on samples from Karnali River to determine rates of exhumation, and they reported evidence for widespread cooling at ca. 16.0 ± 1.4 Ma, with continuous exhumation at a rate of ∼1.4 ± 0.2 km/m.y. Bernet et al. (2006) reported a 4 m.y. lag time between the youngest statistic peaks in the samples from the age of the sediment. The youngest statistic peak from 14 Ma sediments is 18.7 ± 3.2 Ma, and from 12 Ma sediments, it is 16.5 ± 2.2 Ma. The zircon fission-track system closure temperature is ∼240 °C, yielding orogenic cooling rates of 15 °C/m.y. (Hurford, 1986; Brandon et al., 1998). Using our model, we can estimate the magnitude of vertical exhumation between these two time windows. From 19 to 14 Ma, an 8–11 km thickness of rock was removed, while from 16 to 12 Ma, 7.5–8.5 km of rock were removed. Assuming a 240 °C closure temperature, our model would suggest geothermal gradients of 22–32°/km to account for the exhumation and then deposition of these minerals. In addition, we can compare times of pronounced exhumation. As discussed previously, from 10 to 11 Ma exhumation dramatically increased due to the development of a significant ramp and the associated vertical component of deformation. Using the lag time of detrital muscovite, Szulc et al. (2006) showed the highest exhumation rates are between 13 and 10 Ma, with rates between 2 and 2.5 mm/yr. Our rates are also the highest over this window of time, with rates that range from 1.4 to 2.1 mm/yr and possibly as high as 4.2 mm/yr.

Table 5 compares the available exhumation rates from across the Himalaya. What becomes quickly visible is that exhumation rates from 2 Ma to the present are notably higher in central Nepal and NW India than the youngest rate (5–0 Ma) that we propose here. Measurements of Quaternary features in central Nepal indicate rates of exhumation as high as 2–5 mm/yr, similar to modern rates (2–3 mm/yr) inferred from river geochemistry (Galy and France-Lanord, 2001). Thiede et al. (2009) also suggested faster exhumation rates (2–3 mm/yr) from 3 to 0 Ma. Two possible explanations exist for this discrepancy. One is that our broad-brush approach for determining magnitudes of exhumation is not sensitive to a recent (2–0 Ma) increase in exhumation that would be best illustrated with a lower-temperature thermochronometer system like apatite fission track. Our approach averages out exhumation both temporally (5 Ma to present) as well as spatially (Fig. 4). More discrete regions on Figure 4 may have lost up to 7 km of vertical overburden within the past 5 m.y., indicating potential exhumation rates of 1.4 mm/yr. If most of that erosion happened in the past 2–3 m.y., rates could be as high as 2–3.5 mm/yr. The second possibility is that deformation and associated exhumation vary spatially with time along the Himalayan arc. Thiede et al. (2009) showed how exhumation rates vary with time in the Sutlej River valley in NW India, with times of fast exhumation (2–3 mm/yr) that fall between 23 and 19 Ma and 11–0 Ma. The window of time in our study that highlights the fastest exhumation (ca. 13–10 Ma) is a period of slower exhumation in Thiede et al. (2009). Their highest exhumation rates start at 11 Ma but are more focused around 6–4 Ma. Presently, we just do not have the data necessary to differentiate between these two hypotheses.

Table 6 shows the amount of shortening in each time period as highlighted in Figure 4. The shortening amount is derived from both balanced cross sections (Robinson et al., 2006) and 2D Move reconstructions detailed in Robinson (2008). The rate is calculated by dividing the shortening amount by the duration. The timing of the late faults (shown in Fig. 4H) is unknown, although they must postdate the faults they cut. The 10.6 km total offset for these faults listed here is the magnitude achieved by adding and subtracting the offset on the thrust and normal faults, respectively. Adding the total magnitude of the late fault displacement to the smallest rates (Figs. 4D and 4H) only increases the rate by 2–4 mm/yr. Thus, with that logic, the error on the rates is ∼3 mm/yr.

Between 25 and 16 Ma and 13 and 10 Ma, the shortening rate averages ∼27 mm/yr, with shortening rates as high as 33 mm/yr from 13 to 11 Ma (Table 6). In contrast, from 16 to 13 Ma and 10 to 5 Ma, shortening rates are 13–14 mm/yr. From 5 to 0 Ma, shortening rates are 4 mm/yr. Even summing all of the shortening accommodated by the last-stage faults, shortening rates during the last time period (5–0 Ma) are still low (6.4 mm/yr). The slowing of shortening post–5 Ma is a function of the age of displacement on the Main Boundary thrust. The Main Boundary thrust cuts the middle Siwalik unit, indicating that some component of displacement occurred post–5 Ma. An alternative interpretation of the Main Boundary thrust and recent rates is to combine displacement from 10 Ma to present, for an 8.5 mm/yr shortening rate, possibly as high as 9.6 mm/yr, if all of the late-stage faults developed during this time.

The average rate of shortening from 25 Ma to the present is calculated by dividing the total shortening by the total duration, providing a long-term average of 20 mm/yr. This rate agrees with active global positioning system (GPS) shortening rates across the Himalaya of 15–20 mm/yr (e.g., Bettinelli et al., 2006), as well as the Holocene rates of 21.5 ± 2 mm/yr slip rate on the Main Frontal thrust determined from deformed Holocene terraces (Lavé and Avouac, 2000). The matches among the long-term rate, Holocene rates, and GPS rates led to the prevailing assumption that the shortening rate along the Main Himalayan thrust has been constant through time (e.g., Herman et al., 2010). However, Table 6 shows that the rate of shortening varies significantly through time, leading to periods of India-Asia convergence where more or less deformation is being accommodated by the Indian plate, as well as pulsed propagation of the thrust front. The amount of convergence accommodated by India is greater during emplacement of the Main Central thrust, Ramgarh thrust, and during the middle stage of Lesser Himalaya duplex development. The rate is slower during the initial building of the northern Lesser Himalaya duplex, when the first few thrust sheets are being emplaced over smaller ramps in the Main Himalayan thrust and as the toe of the thrust belt propagates to the south, forming the Main Boundary thrust and Subhimalayan thrust system.

Previous to this study, balanced cross sections through the Himalaya were used to provide a long-term average of shortening solely from the initiation of motion on the Main Central thrust to present. By linking both hinterland and foreland data to the emplacement of major structures, we present a method with which to extract a more detailed view of shortening rates and how they vary through time. What emerges is a tempo of deformation that varies, with periods of rapid shortening during the propagation of the Main Central thrust and development of the Lesser Himalaya duplex (∼30 mm/yr). Both early (after the emplacement of the Ramgarh thrust) and late stages of Lesser Himalaya duplex development are marked by windows of slow shortening (∼13–14 mm/yr). As noted by Rahl et al. (2011) in the Pyrenees, fault activity occurs in pulses, with slip occurring contemporaneously on multiple faults throughout the wedge. This reconstruction of fault activity and erosion suggests that orogenic wedges may not evolve in a steady fashion but instead exhibit significant changes in rates of deformation.

An important implication of varying shortening rates is that the partitioning of deformation within the India-Asia collision zone may have changed through time. If the result of 30+ mm/yr of motion of individual faults is robust, then that means that certain faults (Main Central thrust, Ramgarh thrust, and middle stage of Lesser Himalaya duplex development) accommodated a significant component of the 55–40 mm/yr India-Asia convergence rates (e.g., Copley et al., 2010). The approach we use here can be also applied to other orogens with robust data sets in the foreland and hinterland (e.g., Appalachians, Andes, and western United States) to evaluate the presence and possible mechanisms for variations in shortening rates. While proposing variable shortening through time is a new view of Himalayan deformation, it has been documented in other orogens. Linking detailed cross sections with a well-dated foreland basin showed shortening variations between 4 and 13 mm/yr over the past 8 m.y. in the Subandean fold-and-thrust belt of Argentina (Echavarria et al., 2003). Over longer times scales (45 Ma to present), shortening rates though the central Andes have also varied between 4 and 13 mm/yr (McQuarrie et al., 2005; Oncken et al., 2006), although integrating detailed basin studies, thermochronology, and structure with a linked structural and flexural model may facilitate extracting a more detailed history of shortening and exhumation for this region as well.

Using balanced cross sections to determine exhumation has also been limited to long time scales averaged over the duration of deformation due to the uncertainty in partitioning the details of exhumation magnitude with time (Barnes and Pelletier, 2006; McQuarrie et al., 2008). However, by accounting for flexure and making reasonable assumptions about the topographic slope in the past, we can calculate first-order denudation magnitudes through time that can be compared to both thermochronometer systems as well as the associated basin history. This approach is particularly powerful in the Himalaya because the Himalayan foreland basin has preserved minerals and isotopic signatures that can be linked to specific tectonostratigraphic packages.

Exhumation is most strongly linked to the vertical component of deformation (e.g., Whipp et al., 2007), and sequentially restored balanced cross sections are the most accurate models for the locations of active ramps in the past. By linking the sequentially restored cross sections to spatially constrained areas of rapid uplift, we can test the accuracy of the geometry displayed in the balanced cross sections. To a first order, the presence and location of the Lesser Himalaya ramps are supported by the exhumation history. An obvious benefit to linking sequentially restored cross sections to thermochronometry and basin development is that these two data sets provide temporal constraints on the deformation. The only way to determine rates is to link the shortening calculated in the fold-and-thrust belt to the time over which that shortening happens. In this study, exhumation appears consistent except for a period from 10 to 13 Ma, which exhibits an increase in exhumation. This coincides with the development of a 7 km ramp in the Main Himalayan thrust that lifted the overburden to the erosional surface. Thus, development of structures with a significant vertical component affects the exhumation rate.

This study is the first to fully integrate deformation and erosion of the fold-and-thrust belt with flexure of the foreland basin to produce a three-dimensional model, length, height, and time, of the unroofing of the central Himalaya. The result is the ability to link magnitudes of shortening and exhumation to the concurrent deposition within the linked foreland basin–thrust belt system.

A shift in provenance from only Greater Himalaya–Tethyan Himalaya contribution to sediment that contains Lesser Himalaya input is widely recognized in the central Himalaya at ca. 10 Ma. We illustrate that during this time period, Lesser Himalaya and Greater Himalaya rocks were lifted over a large (7 km) Lesser Himalaya ramp. The ramp in the Main Himalayan thrust caused Lesser Himalaya rocks to be elevated, focused erosion in this region, and allowed eroded Lesser Himalaya sediment to be deposited in the Siwalik Group.

In previous studies, it was suggested that Greater Himalaya rocks were exposed by 20 Ma, providing sediment into the foreland basin. By analyzing the data without a bias and by using the hinterland and foreland basin data in combination with the reconstruction, we found that Greater Himalaya rocks were not exposed until ca. 12 Ma along the reconstruction line.

The average shortening rate is ∼20 mm/yr, similar to other studies. However, by identifying critical windows of time that can be accurately dated, we bracket the magnitude of shortening that occurs in these windows of time. What is illuminated is a shortening rate that has varied significantly through time. The highest rates from 25 to 33 mm/yr were during emplacement of the Main Central thrust, Ramgarh thrust, and the middle part of the Lesser Himalaya duplex. After emplacement of the Ramgarh thrust, the early and late stages of Lesser Himalaya duplex development are marked by windows of slow shortening (∼13–14 mm yr) as well as formation of the Main Boundary thrust and Subhimalayan thrust system.

By accounting for both flexure and an evolving topography with time, we can estimate the magnitude of denudation and obtain exhumation rates. The exhumation history we elucidate through our model is consistent with bedrock and detrital thermochronology obtained through the region. Our average exhumation rates are within the values obtained by other studies in the Himalaya, validating the techniques and assumptions used in the model. Even though our approach allows us to get the most information out of existing data, our understanding of the deformation and exhumational history of far western Nepal is still limited by a lack of data. Future work combining more detailed thermochronology and reconstructions that extend along the Himalayan arc can illuminate the four-dimensional mountain building in the Himalaya, providing a more detailed database for dynamic models.

Informal discussions regarding the data presented were held with P. DeCelles, A. Martin, and O. Pearson. Reviews of an early version of this manuscript were provided by D. Burbank, P. van der Beek, and an anonymous reviewer. Reviews of this manuscript were provided by An Yin, an anonymous reviewer, and Editor John Goodge. Robinson acknowledges support from Cooperative Institute for Research and Science (CIRES) at the University of Colorado at Boulder. McQuarrie received support from the Alexander von Humboldt Foundation. Robinson and McQuarrie thank Midland Valley for the use of the program 2D Move.