Jérôme Lavé

Active folding of fluvial terraces across the Siwaliks Hills, Himalayas of central Nepal

We analyze geomorphic evidence of recent crustal deformation in the sub-Himalaya of central Nepal, south of the Kathmandu Basin. The Main Frontal Thrust fault (MFT), which marks the southern edge of the sub-Himalayan fold belt, is the only active structure in that area. Active fault bend folding at the MFT is quantified from structural geology and fluvial terraces along the Bagmati and Bakeya Rivers. Two major and two minor strath terraces are recognized and dated to be 9.2, 2.2, and 6.2, 3.7 calibrated (cal) kyr old, respectively. Rock uplift of up to 1.5 cm/yr is derived from river incision, accounting for sedimentation in the Gangetic plain and channel geometry changes. Rock uplift profiles are found to correlate with bedding dip angles, as expected in fault bend folding. It implies that thrusting along the MFT has absorbed 21 ± 1.5 mm/yr of N-S shortening on average over the Holocene period. The ±1.5 mm/yr defines the 68% confidence interval and accounts for uncertainties in age, elevation measurements, initial geometry of the deformed terraces, and seismic cycle. At the longitude of Kathmandu, localized thrusting along the Main Frontal Thrust fault must absorb most of the shortening across the Himalaya. By contrast, microseismicity and geodetic monitoring over the last decade suggest that interseismic strain is accumulating beneath the High Himalaya, 50–100 km north of the active fold zone, where the Main Himalayan Thrust (MHT) fault roots into a ductile decollement beneath southern Tibet. In the interseismic period the MHT is locked, and elastic deformation accumulates until being released by large (M_w > 8) earthquakes. These earthquakes break the MHT up to the near surface at the front of the Himalayan foothills and result in incremental activation of the MFT.

Fluvial incision and tectonic uplift across the Himalayas of central Nepal

The pattern of fluvial incision across the Himalayas of central Nepal is estimated from the distribution of Holocene and Pleistocene terraces and from the geometry of modern channels along major rivers draining across the range. The terraces provide good constraints on incision rates across the Himalayan frontal folds (Sub-Himalaya or Siwaliks Hills) where rivers are forced to cut down into rising anticlines and have abandoned numerous strath terraces. Farther north and upstream, in the Lesser Himalaya, prominent fill terraces were deposited, probably during the late Pleistocene, and were subsequently incised. The amount of bedrock incision beneath the fill deposits is generally small, suggesting a slow rate of fluvial incision in the Lesser Himalaya. The terrace record is lost in the high range where the rivers are cutting steep gorges. To complement the terrace study, fluvial incision was also estimated from the modern channel geometries using an estimate of the shear stress exerted by the flowing water at the bottom of the channel as a proxy for river incision rate. This approach allows quantification of the effect of variations in channel slope, width, and discharge on the incision rate of a river; the determination of incision rates requires an additional lithological calibration. The two approaches are shown to yield consistent results when applied to the same reach or if incision profiles along nearby parallel reaches are compared. In the Sub-Himalaya, river incision is rapid, with values up to 10–15 mm/yr. It does not exceed a few millimeters per year in the Lesser Himalaya, and rises abruptly at the front of the high range to reach values of ∼4–8 mm/yr within a 50-km-wide zone that coincides with the position of the highest Himalayan peaks. Sediment yield derived from the measurement of suspended load in Himalayan rivers suggests that fluvial incision drives hillslope denudation of the landscape at the scale of the whole range. The observed pattern of erosion is found to closely mimic uplift as predicted by a mechanical model taking into account erosion and slip along the flat-ramp-flat geometry of the Main Himalayan Thrust fault. The morphology of the range reflects a dynamic equilibrium between present-day tectonics and surface processes. The sharp relief together with the high uplift rates in the Higher Himalaya reflects thrusting over the midcrustal ramp rather than the isostatic response to reincision of the Tibetan Plateau driven by late Cenozoic climate change, or late Miocene reactivation of the Main Central Thrust.

Interseismic strain accumulation on the Himalayan crustal ramp (Nepal)

The Departement of Mines and Geology has been monitoring the seismicity of the Central Himalayas of Nepal since 1985. Intense microseismicity and frequent medium‐size earthquakes (mL<4) tend to cluster beneath the topographic front of the Higher Himalaya. This 10–20km deep seismicity also correlates with a zone of localized uplift that has been evidenced from geodetic data. Both microseismic and geodetic data indicate strain accumulation on a mid‐crustal ramp that had been previously inferred from geological and geophysical evidence. This ramp connects a flat decollement under the Lesser and Sub‐Himalaya with a deeper decollement under the Higher Himalaya, and probably acts as a geometric asperity where strain and stress build up during the interseismic period. The large Himalayan earthquakes could nucleate there and probably activate the whole flat‐and‐ramp system up to the blind thrusts of the Sub‐Himalaya.

Investigation of the relationships between basin morphology, tectonic uplift, and denudation from the study of an active fold belt in the Siwalik Hills, central Nepal

The present study investigates correlations between an extensive range of geomorphic properties that can be estimated from a digital elevation model and the uplift rate on geological timescales. The analysis focuses on an area in the Siwalik Hills (central Nepal), where lithology and climate can be considered as uniform. This area undergoes rapid tectonic uplift at rates of up to 15 mm yr^(−1), which are derived from the geometric pattern of a fault-bend model of fold growth. The selected geomorphic properties can be divided in two categories, depending on whether or not the vertical dimension is taken into account. None of the planar properties are significantly correlated to uplift rate, unlike those that include the vertical dimension, such as the mean elevation of basins, hypsometric curve, and hypsometric integral, and relief defined by the amplitude factor of length scaling analysis. Correlation between relief and uplift rate is observed for all length scales of topography shorter than 600 m, which suggests that all orders of the streams are able to adjust to the tectonic signal. Simple mass balance considerations imply that the average elevation is only 10% of surface uplift, suggesting that a dynamic equilibrium has been reached quite rapidly. Using a simple two-process model for erosion, we find that fairly high diffusion coefficients (order of 10 m^2 yr^(−1)) and efficient transport of the material by rivers are required. This unusually high value for mass diffusivity at small length scales may be obtained by either a very efficient linear diffusion or by landsliding. Actually, both processes may be active, which appears likely given the nature of the unconsolidated substratum and the favorable climatic conditions. Local relief in the study area may therefore be used to predict either uplift or denudation, but the prediction is calibrated only for that specific climatic and lithologic conditions and cannot be systematically applied to other contexts.

Seismic anisotropy beneath Tibet: evidence for eastward extrusion of the Tibetan lithosphere?

Strong seismic anisotropy beneath Tibet has recently been reported from the study of SKS shear wave splitting. The fast split waves are generally polarized in an easterly direction, close to the present day direction of motion of the Tibetan crust relative to stable Eurasia, as deduced from Holocene slip rates on the major active faults in and around Tibet. This correlation may be taken to suggest that the whole Tibetan lithosphere is being extruded in front of indenting India and that the anisotropic layer is the deforming asthenosphere, that accommodates the motion of the Tibetan lithosphere relative to the fixed mantle at depth. Uncertainties about this motion are at present too large to bring unambiguous support to that view. Assuming that this view is correct however, a simple forward model is used to compute theoretical delay times as a function of the thickness of the anisotropic layer. The observed delay times would require a 50–100 km thick anisotropic layer beneath south-central Tibet and an over 200 km thick layer beneath north-central Tibet, where particularly hot asthenosphere has been inferred. This study suggests that the asthenospheric anisotropy due to present absolute block motion might be dominant under actively deforming continents.

A Rouse‐based method to integrate the chemical composition of river sediments: Application to the Ganga basin

[1] The Ganga River is one of the main conveyors of sediments produced by Himalayan erosion. Determining the flux of elements transported through the system is essential to understand the dynamics of the basin. This is hampered by the chemical heterogeneity of sediments observed both in the water column and under variable hydrodynamic conditions. Using Acoustic Doppler Current Profiler (ADCP) acquisitions with sediment depth profile sampling of the Ganga in Bangladesh we build a simple model to derive the annual flux and grain size distributions of the sediments. The model shows that ca. 390 (±30) Mt of sediments are transported on average each year through the Ganga at Haring Bridge (Bangladesh). Modeled average sediment grain size parameters D 50 and D 84 are 27 (±4) and 123 (±9) mm, respectively. Grain size parameters are used to infer average chemical compositions of the sediments owing to a strong grain size chemical composition relation. The integrated sediment flux is characterized by low Al/Si and Fe/Si ratios that are close to those inferred for the Himalayan crust. This implies that only limited sequestration occurs in the Gangetic floodplain. The stored sediment flux is estimated to c.a. 10% of the initial Himalayan sediment flux by geochemical mass balance. The associated, globally averaged sedimentation rates in the floodplain are found to be ca. 0.08 mm/yr and yield average Himalayan erosion rate of ca. 0.9 mm/yr. This study stresses the need to carefully address the average composition of river sediments before solving large‐scale geochemical budgets. Citation: Lupker, M., C. France‐Lanord, J. Lave, J. Bouchez, V. Galy, F. Metivier, J. Gaillardet, B. Lartiges, and J.-L. Mugnier (2011), A Rouse‐based method to integrate the chemical composition of river sediments: Application to the Ganga basin,

The kinematics of the Zagros Mountains (Iran)

Abstract We present a synthesis of recently conducted tectonic, global positioning system (GPS), geomorphological and seismic studies to describe the kinematics of the Zagros mountain belt, with a special focus on the transverse right-lateral strike-slip Kazerun Fault System (KFS). Both the seismicity and present-day deformation (as observed from tectonics, geomorphology and GPS) appear to concentrate near the 1000 m elevation contour, suggesting that basement and shallow deformation are related. This observation supports a thick-skinned model of southwestward propagation of deformation, starting from the Main Zagros Reverse Fault. The KFS distributes right-lateral strike-slip motion of the Main Recent Fault onto several segments located in an en echelon system to the east. We observe a marked difference in the kinematics of the Zagros across the Kazerun Fault System. To the NW, in the North Zagros, present-day deformation is partitioned between localized strike-slip motion on the Main Recent Fault and shortening located on the deformation front. To the SE, in the Central Zagros, strike-slip motion is distributed on several branches of the KFS. The decoupling of the Hormuz Salt layer, restricted to the east of the KFS and favouring the spreading of the sedimentary cover, cannot be the only cause of this distributed mechanism because seismicity (and therefore basement deformation) is associated with all active strike-slip faults, including those to the east of the Kazerun Fault System.

Persistence of full glacial conditions in the central Pacific until 15,000 years ago

The magnitude of atmospheric cooling during the Last Glacial Maximum and the timing of the transition into the current interglacial period remain poorly constrained in tropical regions, partly because of a lack of suitable climate records. Glacial moraines provide a method of reconstructing past temperatures, but they are relatively rare in the tropics. Here we present a reconstruction of atmospheric temperatures in the central Pacific during the last deglaciation on the basis of cosmogenic 3He ages of moraines and numerical modelling of the ice cap on Mauna Kea volcano, Hawaii—the only highland in the central Pacific on which moraines that formed during the last glacial period are preserved. Our reconstruction indicates that the Last Glacial Maximum occurred between 19,000 and 16,000 years ago in this region and that temperatures at high elevations were about 7 °C lower than today during this interval. Glacial retreat began about 16,000 years ago, but temperatures were still about 6.5 °C lower than today until 15,000 years ago. When combined with estimates of sea surface temperatures in the central Pacific Ocean, our reconstruction indicates that the lapse rate during the Last Glacial Maximum was higher than at present, which is consistent with the proposal that the atmosphere was drier at that time. Furthermore, the persistence of full glacial conditions until 15,000 years ago is consistent with the relatively late and abrupt transition to warmer temperatures in Greenland, indicating that there may have been an atmospheric teleconnection between the central Pacific and North Atlantic regions during the last deglaciation.

Le cycle sismique en Himalaya

We discuss the seismic cycle in the Himalayas and its relation to mountain building on the basis of geodetic, seismological and geological data collected in the Himalaya of Nepal. On average over several seismic cycles, localized slip on a major thrust fault, the Main Himalayan Thrust fault, MHT, accommodates the ∼21 mm·yr^(−1) convergence rate between southern Tibet and India. The geodetic data show that the MHT is presently locked from the sub-Himalayas to beneath the front of the high range where it roots into a sub-horizontal ductile shear zone under southern Tibet. Aseismic slip during the interseismic period induces stress accumulation at the southern edge of this shear zone triggering intense microseismic activity and elastic straining of the upper crust at the front of the high range. This deformation is released, on the long term, by major earthquakes on the MHT. Such an event is the M_w 8.4-1934-earthquake that ruptured a 250–300-km long segment. The major seismic events along the Himalayas since the 19th century have released more than 70% of the crustal strain accumulated over that period, suggesting that, if any, aseismic slip on the MHT cannot account for more than 30% of the total slip.

Numerical modeling of mountain building: Interplay between erosion law and crustal rheology

Coupling between erosion and tectonics is thought to play a determinant role in orogenic evolution. Here, we investigate the interplay in this coupling between the assumed erosion law and the crustal rheology at the margin of a collisional plateau, like the Himalaya of Central Nepal. Lithospheric deformation is calculated over a time scale of 100 kyr by a 2D finite element model that incorporates the rheological layering of the crust and the main features of the convergence across the range. For the upper boundary condition, two surface processes were tested: a linear diffusion model and a 1D1/2 integrative model including fluvial incision along the fluvial network and hillslope erosion by landsliding. Model results and their sensitivity to the chosen combinations of erosion law and crustal properties are discussed in light of the constraining geologic and geomorphologic observations. In contrast with the conclusions of Cattin and Avouac [2000] , where a compliant quartz-rich crustal rheology and diffusion law were required, we combine a composite quartz-diabase rheology for the crust with fluvial incision erosion law to account for erosion and elevation profiles across the Himalaya of Central Nepal. More generally, it is proposed that, because of the interplay between the dominant denudation conditions and the rheology of the crust, both well documented erosion rates and processes can provide significant constraints on crustal properties within an active orogen.