We studied the formation of the Himalayan mountain range and the Tibetan Plateau by investigating their lithospheric structure. Using an 800-kilometer-long, densely spaced seismic array, we have constructed an image of the crust and upper mantle beneath the Himalayas and the southern Tibetan Plateau. The image reveals in a continuous fashion the Main Himalayan thrust fault as it extends from a shallow depth under Nepal to the mid-crust under southern Tibet. Indian crust can be traced to 31 degrees N. The crust/mantle interface beneath Tibet is anisotropic, indicating shearing during its formation. The dipping mantle fabric suggests that the Indian mantle is subducting in a diffuse fashion along several evolving subparallel structures.

[1] We document geodetic strain across the Nepal Himalaya using GPS times series from 30 stations in Nepal and southern Tibet, in addition to previously published campaign GPS points and leveling data and determine the pattern of interseismic coupling on the Main Himalayan Thrust fault (MHT). The noise on the daily GPS positions is modeled as a combination of white and colored noise, in order to infer secular velocities at the stations with consistent uncertainties. We then locate the pole of rotation of the Indian plate in the ITRF 2005 reference frame at longitude = À 1.34 AE 3.31 , latitude = 51.4 AE 0.3 with an angular velocity of W = 0.5029 AE 0.0072 /Myr. The pattern of coupling on the MHT is computed on a fault dipping 10 to the north and whose strike roughly follows the arcuate shape of the Himalaya. The model indicates that the MHT is locked from the surface to a distance of approximately 100 km down dip, corresponding to a depth of 15 to 20 km. In map view, the transition zone between the locked portion of the MHT and the portion which is creeping at the long term slip rate seems to be at the most a few tens of kilometers wide and coincides with the belt of midcrustal microseismicity underneath the Himalaya. According to a previous study based on thermokinematic modeling of thermochronological and thermobarometric data, this transition seems to happen in a zone where the temperature reaches 350 C. The convergence between India and South Tibet proceeds at a rate of 17.8 AE 0.5 mm/yr in central and eastern Nepal and 20.5 AE 1 mm/yr in western Nepal. The moment deficit due to locking of the MHT in the interseismic period accrues at a rate of 6.6 AE 0.4 Â 10 19 Nm/yr on the MHT underneath Nepal. For comparison, the moment released by the seismicity over the past 500 years, including 14 M W ≥ 7 earthquakes with moment magnitudes up to 8.5, amounts to only 0.9 Â 10 19 Nm/yr, indicating a large deficit of seismic slip over that period or very infrequent large slow slip events. No large slow slip event has been observed however over the 20 years covered by geodetic measurements in the Nepal Himalaya. We discuss the magnitude and return period of M > 8 earthquakes required to balance the long term slip budget on the MHT.

The rocks of the Indian subcontinent are last seen south of the Ganges before they plunge beneath the Himalaya and the Tibetan plateau. They are next glimpsed in seismic reflection profiles deep beneath southern Tibet, yet the surface seen there has been modified by processes within the Himalaya that have consumed parts of the upper Indian crust and converted them into Himalayan rocks. The geometry of the partly dismantled Indian plate as it passes through the Himalayan process zone has hitherto eluded imaging. Here we report seismic images both of the decollement at the base of the Himalaya and of the Moho (the boundary between crust and mantle) at the base of the Indian crust. A significant finding is that strong seismic anisotropy develops above the decollement in response to shear processes that are taken up as slip in great earthquakes at shallower depths. North of the Himalaya, the lower Indian crust is characterized by a high-velocity region consistent with the formation of eclogite, a high-density material whose presence affects the dynamics of the Tibetan plateau.

No abstract

Movie S1Correction: In table S1, the displacement at station SNDL was reported erroneously. The correct displacement is: east, 0.047 ±0.002 m; north, -0.223 ±0.003 m; vertical, 0.003 ±0.003 m. The PDF has been corrected.