In the days to weeks following an earthquake, landslides can display specific post-seismic motions, including delayed initiations and post-seismic relaxations. These motions have an uncertain origin, ...sometimes attributed to specificities of the landslide basal interface or to fluid transports in the landslide basal shear zone. Here we address this question, by documenting the co- and post-seismic motions of slow-moving landslides accelerated by the Gorkha earthquake (Mw 7.8, 25/04/2015, Nepal). We detect 11 slow-moving landslides over an area of 750 kmFormula: see text in the near field of the earthquake, and monitor their motions thanks to a time-series of Pléiades optical satellite images and SAR Sentinel-1 images. The post-seismic landslide motions are much larger than the co-seismic ones, reaching up to Formula: see text m accommodated over 2 months. A delayed initiation of several days (> 4 days) is also measured for at least two of the landslides. We analyze our findings in regards with all the previous observations on slow-moving landslides accelerated by earthquakes, and propose that the post-seismic motions are caused by diffusion of groundwater from co-seismic material contraction up to the landslide basal shear zone or from internal landslide reconfiguration. Our observations strongly suggest the main control of the hydrology in the landslide processes under seismic forcings.
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 D50 and D84 are 27 (±4) and 123 (±9) μm, 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.
Key Points
Chemical heterogeneity of river sediments must be integrated over space and time
Grain size is correlated to chemical composition and can be modeled
Integrated sediment composition is used to solve continental‐scale budgets
The Himalayan range is commonly presented as largely laterally uniform from west to east. However, geological structures, topography, precipitation rate, convergence rates, and low‐temperature ...thermochronological ages all vary significantly along strike. Here, we focus on the interpretation of thermochronological data sets in terms of along‐strike variations in geometry and kinematics of the main crustal detachment underlying the Himalaya: the Main Himalayan Thrust (MHT). We report new apatite fission track (AFT) ages collected along north‐south transects in western and eastern central Nepal (at the latitudes of the Annapurna and Langtang massifs, respectively). AFT ages are consistently young (<3 Ma) along both N‐S transects in the high‐relief zone of the Higher Himalaya and increase (4 to 6 Ma) toward the south in the Lesser Himalaya. We compare our new data to published low‐temperature thermochronological data sets for Nepal and the Bhutan Himalaya. We use the full data set to perform both forward and inverse thermal kinematic modeling with a modified version of the Pecube code in order to constrain potential along‐strike variations in the kinematics of the Himalayan range. Our results show that lateral variations in the geometry of the MHT (in particular the presence or absence of a major crustal‐scale ramp) strongly control the kinematics and exhumation history of the orogen.
Key Points
We report new apatite fission track ages collected in central Nepal
We perform forward and inverse thermal kinematic modeling with PECUBE
Lateral variations in the geometry of the MHT control kinematics and exhumation
We used episodic GNSS measurements to quantify the present‐day velocity field in the northwestern Himalaya from the Himalayan foreland to the Karakoram Range. We report a progressive N‐S ...compressional velocity gradient with two noticeable exceptions: in the Salt Range, where important southward velocities are recorded, and in Nanga Parbat, where an asymmetrical E‐W velocity gradient is recorded. A review of Quaternary slip along active thrusts both in and out of sequence allows us to propose a 14 mm/yr shortening rate. This constraint, together with a geometrical model of the Main Himalayan Thrust (MHT), allows us to propose estimations of the slip distributions along the active faults. The lower flat of the MHT is characterized by ductile slip, whereas the coupling increases along the crustal ramp and along the upper flat of the MHT. The basal thrust of the Potwar Plateau and Salt Range presents weak coupling, which is interpreted as the existence of a massive salt layer forming an excellent décollement. In the central part of the frontal Salt Range, the velocities suggest the existence of a southward horizontal flux in the massive salt layer. The simulations also suggest that the velocities recorded in Nanga Parbat can be explained by active westward thrusting along the fault that borders the massif to the west. Simulations suggest that the slip along this fault evolves with depth from 5 mm/yr ductile slip near the MHT to no slip along the upper part of the fault.
Key Points
Estimation of interseismic deformation using GNSS velocities
Estimation of coupling along the Main Himalayan Thrust in northwestern Himalaya
Quantification of GNSS velocities in northwestern Himalaya and Karakorum ranges
Reconstructing the spatial and temporal response of mountain glaciers to rapid climate change in the past provides access to the effects of current climate change. Yet, the spatial and temporal ...variability of past glacier fluctuations is not fully understood. In this study, we focus on the timing of glacier fluctuations in the European Alps during the Younger Dryas/Early Holocene (YD/EH) transition. In an effort to elucidate whether glacier fluctuations were synchronous during this period, we present a new chronology of the Alpine Talèfre glacier, based on 14 new 10Be ages of moraines and roches moutonnées. The retreat of Talèfre glacier was initiated during the mid‐YD (~12.4 ka), then it experienced a gradual retreat punctuated by at least three oscillations until ∼11 ka before shrinking substantially within its Little Ice Age limits (13th−19th centuries). Comparison of our findings with published glacier chronologies in the Alpine region highlights broadly synchronous behaviour of glaciers across the Alps between 12 and 10 ka. The coeval glacier fluctuations at a regional scale suggest that common regional climate conditions had a major impact on Alpine glacier variations during the YD/EH transition. The similarity of glacier behaviour and independent temperature records in both the Alpine region and the northern high latitudes suggests a teleconnection between these regions, but differences in the amplitude of the mean annual temperature signals relative to summer temperature indicate pronounced changes in seasonality between the YD and the EH.
Investigation of Holocene extents of mountain glaciers along with the related naturally-driven climate conditions helps improve our understanding of glacier sensitivity to ongoing climate change. ...Here, we present the first Holocene glacial chronology in the Mont-Blanc massif (Argentière glacier) in the French Alps, based on 25 in situ-produced cosmogenic 10Be dates of moraines and glacial bedrocks. The obtained ages from mapped sequences of moraines at three locations reveal that the glacier retreated from its Lateglacial extent and oscillated several times between ∼11.7 ka and ∼10.4 ka, i.e. during the Younger Dryas/Early Holocene (YD/EH) transition, before substantially retreating at ∼10.4 ka. Climate conditions corresponding to the past extents of Argentière glacier during the YD/EH transition (∼11 ka) and the Little Ice Age (LIA) were modelled with two different approaches: by determining summer temperature differences from reconstructed ELA-rises and by using a Positive Degree Day (PDD) mass-balance model. The ELA-rise reconstructions yield a possible range of temperatures for the YD/EH transition that were lower by between 3.0 and 4.8 °C compared to the year 2008, depending on the choice of the ELA sensitivity to temperature. The results from the PDD model indicate temperatures lower by ∼3.6–5.5 °C during the YD/EH transition than during the 1979–2002 period. For the LIA, our findings highlight the role of local precipitation changes, superimposed on the dominant temperature signal, in the detailed evolution of the glacier. Overall, this study highlights the challenge that remains in accurately inferring paleoclimate conditions from past glacier extents.
•Argentière glacier oscillated 5 times between 11.7 ka and 10.5 ka.•Modeled Early Holocene temperature range is 3.6–5.5 °C below modern values.•Local precipitation changes explain detailed Little Ice Age glacier fluctuations.
Abstract
The Mont-Blanc massif, being iconic with its large glaciers and peaks of over 4,000 m, will experience a sharp increase in summer temperatures during the twenty-first century. By 2100, the ...impact of climate change on the cryosphere and hydrosphere in the Alps is expected to lead to a decrease in annual river discharge. In this work, we modelled the twenty-first century evolution of runoff in the Arve river, downstream of Mont-Blanc’s French side. For the first time for this region, we have forced a hydrological model with output from an ice-dynamical glacier model and 16 downscaled climate projections, under RCP4.5 and RCP8.5 scenarios. By 2100, under RCP8.5 (high-emission scenario), the winter discharge of the Arve river remains low but is expected to increase by 80% when compared to the beginning of the century. By contrast, the summer season, currently the most important discharge period, will be marked by a runoff decrease of approximately 40%. These changes are almost similar according to a scenario with a lower warming (RCP4.5) and are mostly driven by glacier retreat. These shifts will have significant downstream impacts on water quantity and quality, affecting hydroelectric generation, agriculture, forestry, tourism and aquatic ecosystems.
The morphological boundary between the Himalayas and the foreland plain is well expressed and most often corresponds to the frontal emergence of the Main Himalayan Thrust (MHT). This boundary is ...affected by surface ruptures during very large Himalayan earthquakes (Mw > 8) that regularly induce (with a recurrence of the order of 500 to 1200 years) the uplift of the foothills relative to the plain.
However, a thrust-fold system is hidden beneath the plain and is displayed by the seismic profiles of oil companies in east/central Nepal and by H/V passive geophysical techniques in Darjeeling. Its long-term kinematic evolution is slow, with a tectonic uplift of the hanging wall that is lower than the subsidence rate of the foreland basin, that is, less than approximately half a millimetre per year. During phases of low sedimentation controlled by climatic fluctuations, the morphological surfaces of the piedmont are incised by large rivers for several tens of metres; therefore, structures hidden under the sediments emerge slightly in the plain.
The evolution of the hidden structures corresponds to an embryonic thrust belt mainly affected by a long-term shortening rate of 1.4 +2.5/−1.2 mm·yr−1, that is, 2–20% of the shortening rate of the entire Himalayan thrust system. Nonetheless, the details of the deformation associated with the embryonic thrust belt are still poorly understood. Several deformation components could affect the central Himalayan and Darjeeling piedmonts. i) Any slow steady-state deformation, such as layer parallel shortening (LPS) is not detected by Global Navigation Satellite System (GNSS) data, and such deformation would therefore absorb less than 0.5 mm·yr−1. The geodetic data that suggest the aseismic growth of some of the structures are highly controversial. ii) For the rest of the deformation of the embryonic thrust wedge, it is yet to be proven whether deformation occurs during rare great earthquakes affecting the piedmont during medium earthquakes and/or during post-seismic deformation related to great earthquakes. The amplitude of this long-term low deformation is too limited to significantly reduce the seismic hazard in the seismic gaps of the Himalayan belt. iii) In some portions of the Himalayan front, such as Darjeeling (India), the thrust deformation related to great earthquakes propagates several tens of kilometres south of the morphological front in the zone previously affected by the long-term low deformation. It induces multi-metre surface ruptures in the piedmont and a mean shortening of 8.5 ± 6.2 mm·yr−1. iiii) Pre-existing faults in the bedrock of the Indian craton, often oblique to the Himalayan structures, are locally reactivated beneath the foreland plain with low deformation rates.
•Quantification of interseismic velocities using GNSS and InSAR data.•Coupling estimation along the Main Himalayan Thrust and the Salt Range Thrust.•Deformation induced by the flow towards the ...surface of the detachment level's salt.•Seismic hazard assessment of the Potwar Plateau – Salt Range system.
Using GNSS horizontal surface velocities and Sentinel-1 interferometry line-of-sight velocities, we quantify the current velocity field of the Potwar Plateau Salt Range fold-and-thrust belt. From this velocity field indicating a creep of the Potwar Plateau along the Main Himalayan Thrust, we inferred a weak subhorizontal décollement level formed by a massive Precambrian salt layer. To the south of the Plateau, the Salt Range is uplifted along the Salt Range Thrust, conditioned by the presence of an inherited normal fault. The Kalabagh Fault, which forms the western boundary of the Salt Range and Potwar Plateau, exhibits a creep rate of 3.3 mm/yr. Numerical modelling enabled us to characterise the slip distribution and coupling along the faults, showing the presence of a large asperity along the décollement level beneath the Potwar Plateau and several asperities along the eastern part of this basal thrust. The model also indicated an alternation of coupled and decoupled zones along the Kalabagh Fault, suggesting that this strike-slip fault can be characterised by creep and the occurrence of earthquakes and/or slow slip events. Considering the lack of instrumental and historical large-magnitude earthquakes in this area, the Main Himalayan Thrust and Kalabagh Fault are likely to be affected by earthquakes of magnitude Mw larger than 7.5. A slip rate of 20 mm/yr is modelled along the southern and superficial parts of the Salt Range Thrust, which is larger than the 14 mm/yr slip rate along the Main Himalayan Thrust at depth. This observation suggests the existence of a southward flow of massive salt along the Salt Range Thrust.