Topographic control of asynchronous glacial advances: A case study from Annapurna, Nepal

Differences in the timing of glacial advances, which are commonly attributed to climatic changes, can be due to variations in valley topography. Cosmogenic 10Be dates from 24 glacial moraine boulders in 5 valleys define two age populations, late-glacial and early Holocene. Moraine ages correlate with paleoglacier valley hypsometries. Moraines in valleys with lower maximum altitudes date to the lateglacial, whereas those in valleys with higher maximum altitudes are early Holocene. Two valleys with similar equilibrium-line altitudes (ELAs), but contrasting ages, are \u3c 5 km apart and share the same aspect, such that spatial differences in climate can be excluded. A glacial mass-balance cellular automata model of these two neighboring valleys predicts that change from a cooler-drier to warmer-wetter climate (as at the Holocene onset) would lead to the glacier in the higher altitude catchment advancing, while the lower one retreats or disappears, even though the ELA only shifted by ~120 m

Himalayan landslide-dam lake record

Abstract About 5400 cal yr BP, a large landslide formed a N400-m-tall dam in the upper Marsyandi River, central Nepal. The resulting lacustrine and deltaic deposits stretched N 7 km upstream, reaching a thickness of 120 m. 14 C dating of 7 wood fragments reveals that the aggradation and subsequent incision occurred remarkably quickly (∼ 500 yr). Reconstructed volumes of lacustrine (∼ 0.16 km 3 ) and deltaic (∼ 0.09 km 3 ) deposits indicate a bedload-to-suspended load ratio of 1:2, considerably higher than the ≤1:10 that is commonly assumed. At the downstream end of the landslide dam, the river incised a new channel through ≥ 70 m of Greater Himalayan gneiss, requiring a minimum bedrock incision rate of 13 mm/ yr over last 5400 yr. The majority of incision presumably occurred over a fraction of this time, suggesting much higher rates. The high bedload ratio from such an energetic mountain river is a particularly significant addition to our knowledge of sediment flux in orogenic environments

Pre- and post-seismic deformation related to the 2015, M_w 7.8 Gorkha earthquake, Nepal

We analyze time series from continuously recording GPS stations in Nepal spanning the pre- and post-seismic period associated to the M_w7.8 Gorkha earthquake which ruptured the Main Himalayan Thrust (MHT) fault on April 25th, 2015. The records show strong seasonal variations due to surface hydrology. After corrections for these variations, the time series covering the pre- and post-seismic periods do not show any detectable transient pre-seismic displacement. By contrast, a transient post-seismic signal is clear. The observed signal shows southward displacements consistent with afterslip on the MHT. Using additional data from stations deployed after the mainshock, we invert the time series for the spatio-temporal evolution of slip on the MHT. This modelling indicates afterslip dominantly downdip of the mainshock rupture. Two other regions show significant afterslip: a more minor zone updip of the rupture, and a region between the mainshock and the largest aftershock ruptures. Afterslip in the first ~ 7 months after the mainshock released a moment of [12.8 ± 0.5] × 10^(19) Nm which represents 17.8 ± 0.8% of the co-seismic moment. The moment released by aftershocks over that period of time is estimated to 2.98 × 10^(19) Nm. Geodetically observed post-seismic deformation after co-seismic offset correction was thus 76.7 ± 1.0% aseismic. The logarithmic time evolution of afterslip is consistent with rate-strengthening frictional sliding. According to this theory, and assuming a long-term loading velocity modulated on the basis of the coupling map of the region and the long term slip rate of 20.2 ± 1.1 mm/yr, afterslip should release about 34.0 ± 1.4% of the co-seismic moment after full relaxation of post-seismic deformation. Afterslip contributed to loading the shallower portion of the MHT which did not rupture in 2015 and stayed locked afterwards. The risk for further large earthquakes in Nepal remains high both updip of the rupture area of the Gorkha earthquake and West of Kathmandu where the MHT has remained locked and where no earthquake larger than M_w7.5 has occurred since 1505

Slip pulse and resonance of the Kathmandu basin during the 2015 Gorkha earthquake, Nepal.

This is the author accepted manuscript. The final version is available from AAAS via http://dx.doi.org/10.1126/science.aac6383Detailed geodetic imaging of earthquake ruptures enhances our understanding of earthquake physics and associated ground shaking. The 25 April 2015 moment magnitude 7.8 earthquake in Gorkha, Nepal was the first large continental megathrust rupture to have occurred beneath a high-rate (5-hertz) Global Positioning System (GPS) network. We used GPS and interferometric synthetic aperture radar data to model the earthquake rupture as a slip pulse ~20 kilometers in width, ~6 seconds in duration, and with a peak sliding velocity of 1.1 meters per second, which propagated toward the Kathmandu basin at ~3.3 kilometers per second over ~140 kilometers. The smooth slip onset, indicating a large (~5-meter) slip-weakening distance, caused moderate ground shaking at high frequencies (>1 hertz; peak ground acceleration, ~16% of Earth's gravity) and minimized damage to vernacular dwellings. Whole-basin resonance at a period of 4 to 5 seconds caused the collapse of tall structures, including cultural artifacts.The Nepal Geodetic Array was funded by internal funding to JPA from Caltech and DASE and by the Gordon and Betty Moore Foundation, through Grant GBMF 423.01 to the Caltech Tectonics Observatory and was maintained thanks to NSF Grant EAR 13-5136. Andrew Miner and the PAcific Northwest Geodetic Array (PANGA) at Central Washington University are thanked for technical assistance with the construction and operation of the Tribhuvan University-CWU network. Additional funding for the TU-CWU network came from United Nations Development Programme and Nepal Academy for Science and Technology. The high rate data were recovered thanks to a rapid intervention funded by NASA (US) and the Department of Foreign International Development (UK). We thank Trimble Navigation Ltd and the Vaidya family for supporting the rapid response as well. The accelerometer record at KATNP was provided by USGS. Research at UC Berkeley was funded by the Gordon and Betty Moore Foundation through grant GBMF 3024. A portion of this work was carried out at the Jet Propulsion Laboratory, California Institute of Technology, under a contract with the National Aeronautics and Space Administration. The GPS data were processed by ARIA (JPL) and the Scripps Orbit and Permanent Array Center. The effort at the Scripps Institution of Oceanography was funded by NASA grants NNX14AQ53G and NNX14AT33G. ALOS-2 data were provided under JAXA (Japan) PI Investigations 1148 and 1413. JPA thanks the Royal Society for support. We thank Susan Hough, Doug Given, Irving Flores and Jim Luetgert for contribution to the installation of this station

Testing fluvial erosion models using the transient response of bedrock rivers to tectonic forcing in the Apennines, Italy

The transient response of bedrock rivers to a drop in base level can be used to discriminate between competing fluvial erosion models. However, some recent studies of bedrock erosion conclude that transient river long profiles can be approximately characterized by a transport‐limited erosion model, while other authors suggest that a detachment‐limited model best explains their field data. The difference is thought to be due to the relative volume of sediment being fluxed through the fluvial system. Using a pragmatic approach, we address this debate by testing the ability of end‐member fluvial erosion models to reproduce the well‐documented evolution of three catchments in the central Apennines (Italy) which have been perturbed to various extents by an independently constrained increase in relative uplift rate. The transport‐limited model is unable to account for the catchments’response to the increase in uplift rate, consistent with the observed low rates of sediment supply to the channels. Instead, a detachment‐limited model with a threshold corresponding to the field‐derived median grain size of the sediment plus a slope‐dependent channel width satisfactorily reproduces the overall convex long profiles along the studied rivers. Importantly, we find that the prefactor in the hydraulic scaling relationship is uplift dependent, leading to landscapes responding faster the higher the uplift rate, consistent with field observations. We conclude that a slope‐ dependent channel width and an entrainment/erosion threshold are necessary ingredients when modeling landscape evolution or mapping the distribution of fluvial erosion rates in areas where the rate of sediment supply to channels is low

Pre- and post-seismic deformation related to the 2015, M_w 7.8 Gorkha earthquake, Nepal

We analyze time series from continuously recording GPS stations in Nepal spanning the pre- and post-seismic period associated to the M_w7.8 Gorkha earthquake which ruptured the Main Himalayan Thrust (MHT) fault on April 25th, 2015. The records show strong seasonal variations due to surface hydrology. After corrections for these variations, the time series covering the pre- and post-seismic periods do not show any detectable transient pre-seismic displacement. By contrast, a transient post-seismic signal is clear. The observed signal shows southward displacements consistent with afterslip on the MHT. Using additional data from stations deployed after the mainshock, we invert the time series for the spatio-temporal evolution of slip on the MHT. This modelling indicates afterslip dominantly downdip of the mainshock rupture. Two other regions show significant afterslip: a more minor zone updip of the rupture, and a region between the mainshock and the largest aftershock ruptures. Afterslip in the first ~ 7 months after the mainshock released a moment of [12.8 ± 0.5] × 10^(19) Nm which represents 17.8 ± 0.8% of the co-seismic moment. The moment released by aftershocks over that period of time is estimated to 2.98 × 10^(19) Nm. Geodetically observed post-seismic deformation after co-seismic offset correction was thus 76.7 ± 1.0% aseismic. The logarithmic time evolution of afterslip is consistent with rate-strengthening frictional sliding. According to this theory, and assuming a long-term loading velocity modulated on the basis of the coupling map of the region and the long term slip rate of 20.2 ± 1.1 mm/yr, afterslip should release about 34.0 ± 1.4% of the co-seismic moment after full relaxation of post-seismic deformation. Afterslip contributed to loading the shallower portion of the MHT which did not rupture in 2015 and stayed locked afterwards. The risk for further large earthquakes in Nepal remains high both updip of the rupture area of the Gorkha earthquake and West of Kathmandu where the MHT has remained locked and where no earthquake larger than M_w7.5 has occurred since 1505