Search USGS⌕ Search

Geology topics

Research about Nepal

Source-linked reports with geographic coverage including Nepal.

At least 19 recordsLinked to original sources
USGS2026

Effects of human activity on Khumbu Glacier: Towards a sustainable Everest Base Camp

The Everest Base Camp (EBC), a critical staging area for mountaineers attempting to summit Mount Everest, has seen a significant increase in human activity over the past few decades. The increasing number of climbers, trekkers, and support staff, particularly at the campsite area, has intensified environmental impacts on the Khumbu Glacier. This study investigates the effects of human activity at EBC, focusing on three contributors to glacier melt: global climate change, local fossil fuel consumption, and human urinary discharge. Through field data collected during the spring 2023 climbing season, this study examines the energy released through the use of liquefied petroleum gas (LPG), kerosene, and petrol for cooking, heating, and electricity, as well as the heat generated by human urinary discharge. The total heat energy released from these activities at EBC in the spring of 2023 was 849,174 ± 179,774 MJ, which would be sufficient to melt 2492 ± 528 tons of glacier ice and snow. Trends in land surface temperature at EBC and adjacent environments during the 1991–2023 period were calculated from Landsat 5, 7, and 9 satellite data. The 32-year Landsat record reveals that EBC surface temperature increased by 0.28 °C year −1 , which was roughly twice the rate of warming at the surface of Khumbu Glacier adjacent to EBC and about 15% higher than the temperature increase of the debris-covered terrain immediately to the north of EBC. The findings suggest that addressing these anthropogenic influences could help to preserve the Khumbu Glacier and support the sustainability of mountaineering in the region. This paper also takes a transformative approach and explores stable and safer locations in case there is a necessity to relocate the current EBC. Two sites have been identified southwest of the existing EBC, which, unlike the current supraglacial site, are situated on stable ground.

Everest Base Camp, Khumbu Glacier↗
USGS2026

Leopard occupancy and habitat use in the multi-use Chitwan-Annapurna Landscape, Nepal

We estimated leopard ( Panthera pardus fusca ) occupancy in a multi-use region within Nepal’s Chitwan-Annapurna Landscape to evaluate leopard habitat use and inform conservation planning in areas where most of the species’ habitat occurs outside protected areas. In 2021, sign surveys were conducted along 1277 km of transects distributed among 145 grid cells of 7×7 km within a 7105 km² study area, where 226 leopard signs (pugmarks and scats) were documented. We used an occupancy modeling framework to evaluate the influence of environmental and anthropogenic factors on leopard habitat use. We found that leopard occupancy (ψ = 0.73 ± 0.17 CI) was strongly and positively associated with areas used by wild prey such as red muntjac ( Muntiacus muntjak ), rhesus macaques ( Macaca mulatta ), chital ( Axis axis ), and wild boars ( Sus scrofa ). Our results provide evidence that large carnivores like leopards can persist in human-dominated landscapes when native prey remains abundant, underscoring the need for community-based conservation that sustains both prey and predator populations beyond protected areas. By estimating leopard occupancy outside of protected areas, the research establishes a baseline for developing management strategies to ensure the continued existence of leopards in Nepal's multi-use landscapes.

Chitwan-Annapurna Landscape↗
USGS2025

If you build it, will they come? Assessing the response of tiger populations to elevated conservation efforts in lowland Nepal

Thirteen countries within the distributional range of tigers adopted the St. Petersburg Declaration in 2010, committing to double their tiger populations by 2022. As a signatory to this document, Nepal elevated its tiger conservation efforts soon after the declaration was adopted. Using capture-mark-recapture (CMR) analyses of tiger survey data (2013–2022), we assessed the response of tiger populations to enhanced conservation measures in Banke and Bardia national parks (NP) in western lowland of Nepal. In Bardia NP, estimated tiger numbers increased from 49 in 2013 to 122 in 2022, and estimated population density increased from 3.39 in 2013 to 8.47 tigers/100 km 2 in 2022). Female survival rate was consistently higher than male, with 87.06 % during 2013–2018 and 89.08 % during 2018–2022, compared to 72.9 % and 76.6 % for males in the respective periods. Additionally, the female population grew at a faster rate (6–7 % per year vs. ∼ 2 % per year for males). In Banke NP, tiger abundance increased from 3 to 22 between 2013 and 2022, and population density increased from 0.34 to 2.42/100 km 2 over the same period. The remarkable 7-fold increase in tiger abundance within the Banke NP was facilitated by high male and female survival rates (≥ 90 % per year), as well as high recruitment rates (0.27 individual −1 year −1 ) from 2013 to 2018, leading to a robust (18–20 %) annual population growth rate. Tigers were recolonizing Banke NP at the beginning of this study, and the area offered ample space and a recovering prey base. This facilitated the recruitment of dispersers from the larger population in Bardia NP and allowed survival to be high by eliminating mortality due to density-dependent intra-specific aggression. Our results not only suggest that Nepal has met or exceeded its commitment to the St. Petersburg Declaration in Banke and Bardia NPs, but they highlight the ability of certain tiger populations to respond quickly and positively to enhanced conservation measures given suitable habitat and a healthy prey base.

Banke National Park, Bardia National Park↗
USGS2024

Dynamic occupancy modelling of Asian elephants (Elephas maximus) reveals increasing landscape use in Nepal

Large mammals with general habitat needs can persist throughout mixed used landscapes, however, human-wildlife conflict frequently leads to their restriction to protected areas. Conservation efforts, especially for reducing conflicts with humans, can enhance tolerance of humans towards species like Asian elephants ( Elephas maximus ) in human-dominated landscapes. Here, we examine how elephant use in the Chure Terai Madhesh Landscape (CTML) covering the entire elephant range of Nepal changed between 2012 and 2020 in relationship to protection status and environmental conditions. We systematically surveyed ~ 42,000 km 2 of potential habitat, by dividing the study area into 159 grid cells of 15 × 15 km 2 and recorded elephant signs during the cool, dry season in three years (2012, 2018 and 2020). We analyzed the survey data in a single-species, multi-season (dynamic) occupancy modeling framework to test hypotheses regarding the influence of environmental conditions and protected area status on landscape use by elephants over time. The best-supported model included protected area effects on initial use, colonization, and detection probability as well as temporal variation in colonization and detection probability. Initial use and colonization rates were higher in protected areas, however elephants increasingly used cells located both inside and outside the protected areas, and the difference in use between protected areas and outside declined as elephants use became prevalent across most of the landscape. While elephant use was patchily distributed in the first year of surveys consistent with past descriptions of four sub-populations, elephant use consolidated into a western and eastern region in subsequent years with a gap in their distribution occurring between Chitwan and Bardiya National Parks. Our manuscript highlights the increasing landscape use by elephants in both protected areas and areas outside protected areas and suggests that management interventions that focus on reducing conflicts can promote greater use of both protected areas and areas outside of protected areas.

Scientific Reports↗
USGS2023

Pardus in the press: Drivers of leopard (Panthera pardus fusca) attack occurrence on humans in Nepal

The negative impact of large carnivore presence in human-dominated landscapes manifests as livestock depredation and in extreme cases as attacks on humans. In the case of conflict with leopards in Nepal, attacks resulting in human fatality have become more frequent over time, thus creating an urgent socio-ecological and management issue. We estimated the occurrence of leopard attacks in Nepal from human-leopard conflict cases reported in the media. We used occupancy models to analyse data collected from online news reports on incidents of leopard attacks on humans to explore drivers of leopard attacks on a landscape scale. Our results suggest that the probability of occurrence of leopard attack is associated with human population density, terrain ruggedness and livestock density. The human population density effect may be indicative of a density-dependent relationship, where attacks are more likely in areas where an increased abundance of humans increases encounter rates with leopards. The positive effect of livestock density suggests that livestock may be drawing leopards into human settlements, and consequently increasing the likelihood of attacks on humans. Terrain ruggedness might be offering ideal conditions to facilitate attacks on humans, for example remoteness and high amounts of cover to launch ambush attacks. We provide inference and insights into key determinants of leopard attacks on humans on a landscape scale. These insights can be used to guide future research, inform mitigation measures to reduce leopard attacks and foster a better understanding of the interaction between people and leopards. This study demonstrates the applicability and novelty of using a hierarchical modelling framework applied to freely and publicly available media reports to inform the applied management of human-wildlife conflict at a national scale.

People and Nature↗
USGS2022

Quaternary fans and terraces in the Khumbu Himal south of Mount Everest: their characteristics, age and formation

Large fans and terraces are frequent in the Khumbu Himal within the high Himalayan valleys south of Mt. Everest. These features are composed of massive matrix- and clast-supported diamicts that were formed from both hyperconcentrated flows and coarse-grained debris flows. Cosmogenic radionuclide (CRN) exposure ages for boulders on fans and terraces indicate that periods of fan and terrace formation occurred at c. 16, c. 12, c. 8, c. 4 and c. 1.5 ka, and are broadly coincident with the timing of glaciation in the region. The dating precision is insufficient to resolve whether the surfaces formed before, during or after the correlated glacial advance. However, the sedimentology, and morphostratigraphic and geomorphological relationships suggest that fan and terrace sedimentation in this part of the Himalaya primarily occurs during glacier retreat and is thus paraglacial in origin. Furthermore, modern glacial-lake outburst floods and their associated deposits are common in the Khumbu Himal as the result of glacial retreat during historical times. We therefore suggest that Late Quaternary and Holocene fan and terrace formation and sediment transfer are probably linked to temporal changes in discharge and sediment load caused by glacier oscillations responding to climate change. The timing of major sedimentation events in this region can be correlated with fans and terraces in other parts of the Himalaya, suggesting that major sedimentation throughout the Himalaya is synchronous and tied to regional climatic oscillations. Bedrock incision rates calculated from strath terrace ages average c. 3.9 mm a −1 , suggesting that the overall rate of incision is set by regional uplift.

Imja Khola Valley, Khumbu Himal, Mount Everest↗
USGS2020

Constraining central Himalayan (Nepal) fault geometry through integrated thermochronology and thermokinematic modeling

Constraining the subsurface structural geometry of the central Himalaya continues to prove difficult, even after the 2015 Gorkha earthquake and the resulting insights into the trajectory of the Main Himalayan thrust (MHT). To this end, we apply a thermokinematic model to evaluate four possible balanced cross section geometries based on three estimates of the MHT in central Nepal. We compare the effect of different décollement and duplex geometries on predicted cooling ages and compare these to new and published ages. We find that the best‐fit geometry able to reproduce the cooling ages at the surface is a hinterland‐dipping duplex, which has been translated over a mid‐crustal ramp located ~110 km north of the Main Frontal thrust. We find that the temporal evolution of the duplex and MHT mid‐crustal ramp both play an integral role in producing the observed cooling ages, implying that the common assumption that the active décollement and ramp geometry solely control the distribution of cooling ages is incorrect. Furthermore, results indicate that the Ramgarh‐Munsiari thrust was emplaced between 17 and ~10 Ma, followed by the Trishuli thrust. Duplex growth occurs between 6.5 Ma and 0.75 Ma, with its constituent thrust sheets moving at variable rates between 10 – 42 mm/yr. Young out‐of‐sequence thrusting (5 km of displacement) in the hinterland produces a slightly improved fit to the cooling ages. Finally, the resulting thermal field modeled from our best‐fit geometry suggests a possible basis for the nucleation and rupture characteristics of the Gorkha earthquake.

Himalayan Fault↗
USGS2018

Habitat overlap between Asiatic black bear Ursus thibetanus and red panda Ailurus fulgens in Himalaya

Studying habitat overlap between sympatric species is one of the best ways to identify interspecies relationships and to direct conservation efforts so that multiple species can benefit. However, studies exploring interspecies relationships are very limited in Nepal, making it difficult for the government of Nepal and conservation partners to manage wildlife in their habitats, especially in Himalayan protected areas. In this study, we identified habitat overlap between Asiatic black bear ( Ursus thibetanus ) and red panda ( Ailurus fulgens ) as well as important habitat types for both species in the Makalu Barun National Park, Nepal using Maximum Entropy (MaxEnt) modeling. GPS points of species occurrence were collected from the field, and environmental variables were extracted from freely available sources. We found that the study area contained 647 km 2 of Asiatic black bear habitat and 443 km 2 of the red panda habitat. 368 km 2 supported both species, which constituted 57% of the Asiatic black bear habitat and 83% of the red panda habitat. We found that conifer forest was the most important habitat type for both species. Because the largest portions of both species’ habitat were located inside the buffer zone, a peripheral zone of national park, conservation efforts for these sympatric species should be focused inside the buffer zone to be most effective.

PLoS ONE↗
USGS2018

The size, distribution, and mobility of landslides caused by the 2015 Mw7.8 Gorkha earthquake, Nepal

Coseismic landslides pose immediate and prolonged hazards to mountainous communities, and provide a rare opportunity to study the effect of large earthquakes on erosion and sediment budgets. By mapping landslides using high-resolution satellite imagery, we find that the 25 April 2015 M w 7.8 Gorkha earthquake and aftershock sequence produced at least 25,000 landslides throughout the steep Himalayan Mountains in central Nepal. Despite early reports claiming lower than expected landslide activity, our results show that the total number, area, and volume of landslides associated with the Gorkha event are consistent with expectations, when compared to prior landslide-triggering earthquakes around the world. The extent of landsliding mimics the extent of fault rupture along the east-west trace of the Main Himalayan Thrust and increases eastward following the progression of rupture. In this event, maximum modeled Peak Ground Acceleration (PGA) and the steepest topographic slopes of the High Himalaya are not spatially coincident, so it is not surprising that landslide density correlates neither with PGA nor steepest slopes on their own. Instead, we find that the highest landslide density is located at the confluence of steep slopes, high mean annual precipitation, and proximity to the deepest part of the fault rupture from which 0.5–2 Hz seismic energy originated. We suggest that landslide density was determined by a combination of earthquake source characteristics, slope distributions, and the influence of precipitation on rock strength via weathering and changes in vegetation cover. Determining the relative contribution of each factor will require further modeling and better constrained seismic parameters, both of which are likely to be developed in the coming few years as post-event studies evolve. Landslide mobility, in terms of the ratio of runout distance to fall height, is comparable to small volume landslides in other settings, and landslide volume-runout scaling is consistent with compilations of data on larger slope failures. In general, the size ratios of landslide source area to full landslide area are smaller than global averages, and hillslope length seems to largely control runout distance, which we propose reflects a topographic control on landslide mobility in this setting. We find that landslide size dictates runout distance and that more than half of the landslide debris was deposited in direct connection with stream channels. Connectivity, which is defined as the spatial proximity of landslides to fluvial channels, is greatest for larger landslides in the high-relief part of the High Himalaya. Although these failures are less abundant than those at lower elevations, they may have a disproportionate impact on sediment dynamics and cascading hazards, such as landslide reactivation by monsoon rainfall and landslide dams that lead to outburst floods. The overall high fluvial connectivity of coseismic landsliding in the Gorkha event suggests coupling between the earthquake cycle and sediment/geochemical budgets of fluvial systems in the Himalaya.

Geomorphology↗
USGS2018

The 2015 Gorkha (Nepal) Earthquake sequence: I. Source modeling and deterministic 3D ground shaking

To better quantify the relatively long period (< 0.3 Hz) shaking experienced during the 2015 Gorkha (Nepal) earthquake sequence, we study the finite rupture processes and the associated 3D ground motion of the Mw7.8 mainshock and the Mw7.2 aftershock . The 3D synthetics are then used in the broadband ground shaking in Kathmandu with a hybrid approach, summarized in a companion paper (Chen and Wei, 2017, submitted together). We determined the coseismic rupture process of the mainshock by joint inversion of InSAR/SAR, GPS (static and high-rate), strong motion and teleseismic waveforms. Our inversion for the mainshock indicates unilateral rupture towards the ESE, with an average rupture speed of 3.0 km/s and a total duration of ~ 60 s. Additionally, we find that the beginning part of the rupture (5–18 s) has about 40% longer rise time than the rest of the rupture, as well as slower rupture velocity. Our model shows two strong asperities occurring ~ 24 s and ~ 36 s after the origin and located ~ 30 km to the northwest and northeast of the Kathmandu valley, respectively. In contrast, the Mw7.2 aftershock is more compact both in time and space, as revealed by joint inversion of teleseismic body waves and InSAR data. The different rupture features between the mainshock and the aftershock could be related to difference in fault zone structure. The mainshock and aftershock ground motions in the Kathmandu valley, recorded by both strong motion and high-rate GPS stations, exhibited strong amplification around 0.2 Hz. A simplified 3D basin model, calibrated by an Mw5.2 aftershock, can match the observed waveforms reasonably well at 0.3 Hz and lower frequency. The 3D simulations indicate that the basin structure trapped the wavefield and produced an extensive ground vibration. Our study suggests that the combination of rupture characteristics and propagational complexity are required to understand the ground shaking produced by hazardous earthquakes such as the Gorkha event.

Tectonophysics↗
USGS2017

Using interviews and biological sign surveys to infer seasonal use of forested and agricultural portions of a human-dominated landscape by Asian elephants in Nepal

Understanding how wide-ranging animals use landscapes in which human use is highly heterogeneous is important for determining patterns of human–wildlife conflict and designing mitigation strategies. Here, we show how biological sign surveys in forested components of a human-dominated landscape can be combined with human interviews in agricultural portions of a landscape to provide a full picture of seasonal use of different landscape components by wide-ranging animals and resulting human–wildlife conflict. We selected Asian elephants ( Elephas maximus ) in Nepal to illustrate this approach. Asian elephants are threatened throughout their geographic range, and there are large gaps in our understanding of their landscape-scale habitat use. We identified all potential elephant habitat in Nepal and divided the potential habitat into sampling units based on a 10 km by 10 km grid. Forested areas within grids were surveyed for signs of elephant use, and local villagers were interviewed regarding elephant use of agricultural areas and instances of conflict. Data were analyzed using single-season and multi-season (dynamic) occupancy models. A single-season occupancy model applied to data from 139 partially or wholly forested grid cells estimated that 0.57 of grid cells were used by elephants. Dynamic occupancy models fit to data from interviews across 158 grid cells estimated that monthly use of non-forested, human-dominated areas over the preceding year varied between 0.43 and 0.82 with a minimum in February and maximum in October. Seasonal patterns of crop raiding by elephants coincided with monthly elephant use of human-dominated areas, and serious instances of human–wildlife conflict were common. Efforts to mitigate human–elephant conflict in Nepal are likely to be most effective if they are concentrated during August through December when elephant use of human-dominated landscapes and human–elephant conflict are most common.

Ethology Ecology and Evolution↗
USGS2017

Toward a unified near-field intensity map of the 2015 Mw 7.8 Gorkha, Nepal, Earthquake

We develop a unified near-field shaking intensity map for the 25 April 2015 M w 7.8 Gorkha, Nepal, earthquake by synthesizing intensities derived from macroseismic effects that were determined by independent groups using a variety of approaches. Independent assessments by different groups are generally consistent, with minor differences that are likely due in large part to differences in spatial sampling. Throughout most of the near-field region, European Macroseismic Scale (EMS-98) intensities were generally close to 7 EMS. In the Kathmandu Valley, intensities were somewhat higher (6.5–7.5) along the periphery of the valley and in the adjacent foothills than in the central valley, where they were ≈6. The results are consistent with instrumental intensity values estimated from available data using a published relationship between peak ground acceleration (PGA) and intensity. Using this relationship to convert intensities to PGA, we estimate strong-motion PGA de-amplification factors of ≈0.7 in the central Kathmandu Valley, with amplification of ≈1.6 in adjacent foothills. The results support the conclusion that the Kathmandu Valley experienced a pervasively nonlinear response during the Gorkha main shock.

Earthquake Spectra↗
USGS2017

Using structural damage statistics to derive macroseismic intensity within the Kathmandu valley for the 2015 M7.8 Gorkha, Nepal earthquake

We make and analyze structural damage observations from within the Kathmandu valley following the 2015 M7.8 Gorkha, Nepal earthquake to derive macroseismic intensities at several locations including some located near ground motion recording sites. The macroseismic intensity estimates supplement the limited strong ground motion data in order to characterize the damage statistics. This augmentation allows for direct comparisons between ground motion amplitudes and structural damage characteristics and ultimately produces a more constrained ground shaking hazard map for the Gorkha earthquake. For systematic assessments, we focused on damage to three specific building categories: (a) low/mid-rise reinforced concrete frames with infill brick walls, (b) unreinforced brick masonry bearing walls with reinforced concrete slabs, and (c) unreinforced brick masonry bearing walls with partial timber framing. Evaluating dozens of photos of each construction type, assigning each building in the study sample to a European Macroseismic Scale (EMS)-98 Vulnerability Class based upon its structural characteristics, and then individually assigning an EMS-98 Damage Grade to each building allows a statistically derived estimate of macroseismic intensity for each of nine study areas in and around the Kathmandu valley. This analysis concludes that EMS-98 macroseismic intensities for the study areas from the Gorkha mainshock typically were in the VII–IX range. The intensity assignment process described is more rigorous than the informal approach of assigning intensities based upon anecdotal media or first-person accounts of felt-reports, shaking, and their interpretation of damage. Detailed EMS-98 macroseismic assessments in urban areas are critical for quantifying relations between shaking and damage as well as for calibrating loss estimates. We show that the macroseismic assignments made herein result in fatality estimates consistent with the overall and district-wide reported values.

Kathmandu Valley↗
USGS2017

Source modeling of the 2015 Mw 7.8 Nepal (Gorkha) earthquake sequence: Implications for geodynamics and earthquake hazards

The Gorkha earthquake on April 25th, 2015 was a long anticipated, low-angle thrust-faulting event on the shallow décollement between the India and Eurasia plates. We present a detailed multiple-event hypocenter relocation analysis of the Mw 7.8 Gorkha Nepal earthquake sequence, constrained by local seismic stations, and a geodetic rupture model based on InSAR and GPS data. We integrate these observations to place the Gorkha earthquake sequence into a seismotectonic context and evaluate potential earthquake hazard. Major results from this study include (1) a comprehensive catalog of calibrated hypocenters for the Gorkha earthquake sequence; (2) the Gorkha earthquake ruptured a ~ 150 × 60 km patch of the Main Himalayan Thrust (MHT), the décollement defining the plate boundary at depth, over an area surrounding but predominantly north of the capital city of Kathmandu (3) the distribution of aftershock seismicity surrounds the mainshock maximum slip patch; (4) aftershocks occur at or below the mainshock rupture plane with depths generally increasing to the north beneath the higher Himalaya, possibly outlining a 10–15 km thick subduction channel between the overriding Eurasian and subducting Indian plates; (5) the largest Mw 7.3 aftershock and the highest concentration of aftershocks occurred to the southeast the mainshock rupture, on a segment of the MHT décollement that was positively stressed towards failure; (6) the near surface portion of the MHT south of Kathmandu shows no aftershocks or slip during the mainshock. Results from this study characterize the details of the Gorkha earthquake sequence and provide constraints on where earthquake hazard remains high, and thus where future, damaging earthquakes may occur in this densely populated region. Up-dip segments of the MHT should be considered to be high hazard for future damaging earthquakes.

Gorkha↗
USGS2017

Characterizing the Kathmandu Valley sediment response through strong motion recordings of the 2015 Gorkha earthquake sequence

We analyze strong motion records and high-rate GPS measurements of the M 7.8 Gorkha mainshock, M 7.3 Dolakha, and two moderate aftershock events recorded at four stations on the Kathmandu basin sediments, and one on rock-outcrop. Recordings on soil from all four events show systematic amplification relative to the rock site at multiple frequencies in the 0.1–2.5 Hz frequency range, and de-amplification of higher frequencies ( >2.5–10 Hz). The soil-to-rock amplification ratios for the M 7.8 and M 7.3 events have lower amplitude and frequency peaks relative to the ratios of the two moderate events, effects that could be suggestive of nonlinear site response. Further, comparisons to ground motion prediction equations show that 1) both soil and rock mainshock recordings were severely depleted of high frequencies, and 2) the depletion at high frequencies is not present in the aftershocks. These observations indicate that the high frequency deamplification is additionally related to characteristics of the source that are not captured by simplified ground motion prediction equations, and allude to seismic hazard analysis models being revised – possibly by treating isolated high frequency radiation sources separately from long period components to capture large magnitude near-source events such as the 2015 Gorkha mainshock.

Kathmandu Valley↗
USGS2017

Application and evaluation of a rapid response earthquake-triggered landslide model to the 25 April 2015 Mw 7.8 Gorkha earthquake, Nepal

The 25 April 2015 M w 7.8 Gorkha earthquake produced strong ground motions across an approximately 250 km by 100 km swath in central Nepal. To assist disaster response activities, we modified an existing earthquake-triggered landslide model based on a Newmark sliding block analysis to estimate the extent and intensity of landsliding and landslide dam hazard. Landslide hazard maps were produced using Shuttle Radar Topography Mission (SRTM) digital topography, peak ground acceleration (PGA) information from the U.S. Geological Survey (USGS) ShakeMap program, and assumptions about the regional rock strength based on end-member values from previous studies. The instrumental record of seismicity in Nepal is poor, so PGA estimates were based on empirical Ground Motion Prediction Equations (GMPEs) constrained by teleseismic data and felt reports. We demonstrate a non-linear dependence of modeled landsliding on aggregate rock strength, where the number of landslides decreases exponentially with increasing rock strength. Model estimates are less sensitive to PGA at steep slopes (> 60°) compared to moderate slopes (30–60°). We compare forward model results to an inventory of landslides triggered by the Gorkha earthquake. We show that moderate rock strength inputs over estimate landsliding in regions beyond the main slip patch, which may in part be related to poorly constrained PGA estimates for this event at far distances from the source area. Directly above the main slip patch, however, the moderate strength model accurately estimates the total number of landslides within the resolution of the model (landslides ≥ 0.0162 km 2 ; observed n = 2214, modeled n = 2987), but the pattern of landsliding differs from observations. This discrepancy is likely due to the unaccounted for effects of variable material strength and local topographic amplification of strong ground motion, as well as other simplifying assumptions about source characteristics and their relationship to landsliding.

Tectonophysics↗
USGS2016

A comparison of observed and predicted ground motions from the 2015 M W 7.8 Gorkha, Nepal, earthquake

We use 21 strong motion recordings from Nepal and India for the 25 April 2015 moment magnitude (M W ) 7.8 Gorkha, Nepal, earthquake together with the extensive macroseismic intensity data set presented by Martin et al. (Seism Res Lett 87:957–962, 2015 ) to analyse the distribution of ground motions at near-field and regional distances. We show that the data are consistent with the instrumental peak ground acceleration (PGA) versus macroseismic intensity relationship developed by Worden et al. (Bull Seism Soc Am 102:204–221, 2012 ), and use this relationship to estimate peak ground acceleration from intensities (PGA EMS ). For nearest-fault distances (R RUP < 200 km), PGA EMS is consistent with the Atkinson and Boore (Bull Seism Soc Am 93:1703–1729, 2003 ) subduction zone ground motion prediction equation (GMPE). At greater distances (R RUP > 200 km), instrumental PGA values are consistent with this GMPE, while PGA EMS is systematically higher. We suggest the latter reflects a duration effect whereby effects of weak shaking are enhanced by long-duration and/or long-period ground motions from a large event at regional distances. We use PGA EMS values within 200 km to investigate the variability of high-frequency ground motions using the Atkinson and Boore (Bull Seism Soc Am 93:1703–1729, 2003 ) GMPE as a baseline. Across the near-field region, PGA EMS is higher by a factor of 2.0–2.5 towards the northern, down-dip edge of the rupture compared to the near-field region nearer to the southern, up-dip edge of the rupture. Inferred deamplification in the deepest part of the Kathmandu valley supports the conclusion that former lake-bed sediments experienced a pervasive nonlinear response during the mainshock (Dixit et al. in Seismol Res Lett 86(6):1533–1539, 2015 ; Rajaure et al. in Tectonophysics, 2016 . Ground motions were significantly amplified in the southern Gangetic basin, but were relatively low in the northern basin. The overall distribution of ground motions and damage during the Gorkha earthquake thus reflects a combination of complex source, path, and site effects. We also present a macroseismic intensity data set and analysis of ground motions for the M W 7.3 Dolakha aftershock on 12 May 2015, which we compare to the Gorkha mainshock and conclude was likely a high stress-drop event.

Gorkha District↗