USGS Science⌕ Search

Geology topics

Faults & earthquakes

Explore USGS earthquake and fault publications. Search fault names, locations and seismic hazards; follow the original report for maps, methods and publication context.

At least 1,315 records · Page 73Linked to original sources

Distinct yet adjacent earthquake sequences near the Mendocino Triple Junction: 20 December 2021 Mw 6.1 and 6.0 Petrolia, and 20 December 2022 Mw 6.4 Ferndale

Two earthquake sequences occurred a year apart at the Mendocino Triple Junction in northern California: first the 20 December 2021 �w 6.1 and 6.0 Petrolia sequence, then the 20 December 2022 �w 6.4 Ferndale sequence. To delineate active faults and understand the relationship between these sequences, we applied an automated deep‐learning workflow to create enhanced and relocated earthquake catalogs for both the sequences. The enhanced catalog newly identified more than 14,000 M 0–2 earthquakes and also found 852 of 860 already cataloged events. We found that deep‐learning and template‐matching approaches complement each other to improve catalog completeness because deep learning finds more M 0–2 background seismicity, whereas template‐matching finds the smallest M < 0 events near already known events. The enhanced catalog revealed that the 2021 Petrolia and 2022 Ferndale sequences were distinct in space and time, but adjacent in space. Though both the sequences happened in the downgoing Gorda slab, the shallower Ferndale sequence ruptured within the uppermost slab near the subduction interface, while the onshore Petrolia sequence occurred deeper in the mantle. Deep‐learning‐enhanced earthquake catalogs could help monitor evolving earthquake sequences, identify detailed seismogenic fault structures, and understand space–time variations in earthquake rupture and sequence behavior in a complex tectonic setting.

California↗

Rapid characterization of the 2025 Mw 8.8 Kamchatka, Russia earthquake

The 29 July 2025 M w 8.8 Kamchatka, Russia, earthquake was the sixth largest instrumentally recorded earthquake. This event was seismically well observed at regional and teleseismic distances, but publicly available near‐source data were sparse at the time of the event, presenting unique challenges for rapid source and impact characterization. The U.S. Geological Survey (USGS) National Earthquake Information Center provides global real‐time monitoring for earthquakes, including rapid response information products that estimate source characteristics, shaking, and the resulting impacts. We describe the USGS rapid response earthquake information products following the Kamchatka event and discuss their implications for ongoing hazards in the region. We describe potential improvements to our response workflows motivated by this event, including more rapid constraints on source geometries and the automated selection of fault geometries for finite‐fault inversions. The rapid response products together support the interpretation of a unilateral southwestward rupture with significant slip on the southwestern end of the rupture extent. The M w 8.8–9.0 event in 1952, which ruptured a comparable extent of the Kuril–Kamchatka subduction interface, has many similarities to the 2025 rupture. This illustrates that slip deficits may remain following great earthquakes and highlights the usefulness of comparative studies between historic and modern events.

Kuril-Kamchatka subduction zone↗

Criticism of some forecasts of the National Earthquake Prediction Evaluation Council

The Working Group on California Earthquake Probabilities has assigned probabilities for rupture in the interval from 1988 to 2018 to various segments of the San Andreas fault on the basis of the lognormal distribution of recurrence times of characteristic earthquakes postulated by Nishenko and Buland (1987). I question the validity of those probabilities on the basis of three separate arguments: (1) The distributions of recurrence times of the four, best-observed, characteristic-earthquake sequences are each only marginally consistent with the Nishenko - Buland Iognormal distribution. (2) The range of possible 30-year conditional probabilities for many of the fault segments is so great due to uncertainty in the average recurrence time for that segment that the assigned probability is virtually meaningless. (3) The 1988 forecasts not subject to the foregoing objection are those in which there is a low probability of an earthquake in the near future (e.g., only a 5 per cent chance of rupture of the North Coast segment before the year 2049 and of the Carrizo segment before the year 2018). The same reasoning would assign only a 5 per cent chance of rupture before mid-1993 to the southern Santa Cruz Mountains segment, the segment that failed in October 1989. Finally, the forecast of the next Parkfield earthquake (95 per cent probability before 1993.0) by Bakun and Lindh (1985) depends upon an ad hoc explanation of the out-of-sequence 1934 earthquake. A less-contrived forecast would have assigned a conditional probability of about 60 ± 20 per cent to the 1985.0 to 1993.0 interval and 30 ± 15 per cent to the 1990.0 to 1993.0 interval.

California↗

U.S. Geological Survey Earthquake Science Center

The mission of the U.S. Geological Survey (USGS) Earthquake Science Center is to collect a wide range of data on earthquakes, faults, and crustal deformation; conduct research to increase our understanding of earthquake source processes, occurrence, and effects; and synthesize this knowledge into probabilistic seismic hazard assessments, aftershock forecasts, and ground-shaking scenarios for anticipated major earthquakes. We disseminate these data products, hazard assessments, and research discoveries to scientists, engineers, emergency managers, and the public. With university and State partners, we operate the California Integrated Seismic Network and the Pacific Northwest Seismic Network, as well as geodetic networks throughout the western United States. We also lead the National Strong Motion Project and the ShakeAlert earthquake early warning (EEW) system; house renowned rock mechanics laboratories and deep borehole geophysics facilities; and conduct extensive geophysical, geologic, and paleoseismic investigations along active faults. We are funded primarily by the USGS Earthquake Hazards Program, with additional support from the USGS Volcano Hazards and Energy Resources Programs, other Federal and State agencies, private foundations, and public and private utilities and corporations.

California↗

Map and database of Quaternary faults in Venezuela and its offshore regions

As part of the International Lithosphere Program’s “World Map of Major Active Faults,” the U.S. Geological Survey is assisting in the compilation of a series of digital maps of Quaternary faults and folds in Western Hemisphere countries. The maps show the locations, ages, and activity rates of major earthquake-related features such as faults and fault-related folds. They are accompanied by databases that describe these features and document current information on their activity in the Quaternary. The project is a key part of the Global Seismic Hazards Assessment Program (ILP Project II-0) for the International Decade for Natural Hazard Disaster Reduction. The project is sponsored by the International Lithosphere Program and funded by the USGS’s National Earthquake Hazards Reduction Program. The primary elements of the project are general supervision and interpretation of geologic/tectonic information, data compilation and entry for fault catalog, database design and management, and digitization and manipulation of data in †ARCINFO. For the compilation of data, we engaged experts in Quaternary faulting, neotectonics, paleoseismology, and seismology.

Open-File Report↗

Data for Quaternary faults in western Montana

The "World Map of Major Active Faults" Task Group is compiling published fault data, developing a digital database of the fault data, and preparing a series of maps for the United States and other countries in the western Hemisphere. The data is intended to portray the locations, ages, and activity rates of major earthquake-related features such as faults, folds, and liquefaction features that have geologic evidence of Quaternary (1.6 Ma) deformation. The Western Hemisphere effort is sponsored by International Lithosphere Program (ILP) Task Group II-2; the data compilation, database, and map for the United States is funded largely by the National Earthquake Hazard Reduction Program (NEHRP) through the U.S. Geological Survey. The ILP effort in the Western Hemisphere is coordinated by Michael N. Machette, the digital database is designed and managed Kathleen M. Haller, and map data are digitized and manipulated by Richard L. Dart. In addition to meeting the goals of the Task Group II-2, this effort represents a key contribution to the new Global Seismic Hazards Assessment Program (ILP Task Group II-0) for the International Decade for Natural Disaster Reduction. This compilation, which documents the published data on Quaternary surface faulting in western Montana, is one of many similar state or regional compilations that are planned for the project. Compilations for Arizona (Pearthree, 1998 #2945), Colorado (Widmann and others, 1998 #3441), New Mexico (Machette and others, 1998), and West Texas (Collins and others, 1996 #993) are currently available and the compilation for features east of the Rocky Mountain front will be available in early 2000 (Crone and Wheeler, in press). All are primarily a catalog of data that includes a variety of geographic, geologic, and paleoseismologic parameters for known or assumed Quaternary faults. These data compilations, the digital database, and the companion maps summarize the published information on known tectonic features and present the information in an internally consistent format. The compilations will be available in digital database format on the WorldWide Web in the near future, which will greatly improve their utility. Release of data for individual states and regions within the United States in this text-based format was necessary because of the time required to develop the national database.

Montana↗

Principal facts for gravity data along the Hayward fault and vicinity, San Francisco Bay area, northern California

The U.S. Geological Survey (USGS) established over 940 gravity stations along the Hayward fault and vicinity. The Hayward fault, regarded as one of the most hazardous faults in northern California (Working Group on California Earthquake Probabilities, 1999), extends for about 90 km from Fremont in the southeast to San Pablo Bay in the northwest. The Hayward fault is predominantly a right-lateral strike-slip fault that forms the western boundary of the East Bay Hills. These data and associated physical property measurement were collected as part of on-going studies to help determine the earthquake hazard potential of major faults within the San Francisco Bay region. Gravity data were collected between latitude 37°30' and 38°15' N and longitude 121°45' and 122°30' W. Gravity stations were located on the following 7.5 minute quadrangles: Newark, Niles, San Leandro, Hayward, Dublin, Oakland West, Oakland East, Las Trampas Ridge, Diablo, Richmond, Briones Valley, Walnut Creek, and Clayton. All data were ultimately tied to primary gravity base station Menlo Park A, located on the campus of the U.S. Geological Survey in Menlo Park, Calif. (latitude 37°27.34' N, longitude 122°10.18' W, observed gravity value 979944.27 mGal).

California↗

A preliminary study of the Coyote Lake earthquake of August 6, 1979, and its major aftershocks

The August 6, 1979 earthquake (magnitude 5.7 + 0.2) was located about 1 km east of the Calaveras fault trace near Coyote Lake, at a depth of about 10 km. The spatial distribution of the aftershocks in the first 15 days and the focal mechanism of the main shock and selected aftershocks suggest that faulting took place primarily to the southeast of the main shock along a 25-km segment of the Calaveras fault zone at a depth of 4 to 12 km. The motion was right-lateral strike-slip along a nearly vertical fault plane. Discontinuous surface rupture was observed along a 14.4-km length of the recently active trace of the fault southeast of Coyote Lake, but no surface faulting was seen in the epicentral area. A maximum 5-mm fault offset was recorded near San Felipe Lake, about 16 km south of the main shock. There was no prominent foreshock activity. The cumulative number of aftershocks (N) versus magnitude (M) could be described by the equation log N= 3.33 - 0.75 M, in the magnitude range between 1 and 3.5. The number and size of aftershocks appeared to be much less than that expected for a main shock of magnitude 5.7.

California↗

Energy budgets of mining-induced earthquakes and their interactions with nearby stopes

In the early 1960's, N.G.W. Cook, using an underground network of geophones, demonstrated that most Witwatersrand tremors are closely associated with deep level gold mining operations. He also showed that the energy released by the closure of the tabular stopes at depths of the order of 2 km was more than sufficient to account for the mining-induced earthquakes. I report here updated versions of these two results based on more modern underground data acquired in the Witwatersrand gold fields. Firstly, an extensive suite of in situ stress data indicate that the ambient state of crustal stress here is close to the failure state in the absence of mining even though the tectonic setting is thoroughly stable. Mining initially stabilizes the rock mass by reducing the pore fluid pressure from its initial hydrostatic state to nearly zero. The extensive mine excavations, as Cook showed, concentrate the deviatoric stresses, in localized regions of the abutments, back into a failure state resulting in seismicity. Secondly, there appears to be two distinct types of mining-induced earthquakes: the first is strongly coupled to the mining and involves shear failure plus a coseismic volume reduction; the second type is not evidently coupled to any particular mine face, shows purely deviatoric failure, and is presumably caused by more regional changes in the state of stress due to mining. Thirdly, energy budgets for mining induced earthquakes of both types indicate that, of the available released energy, only a few per cent is radiated by the seismic waves with the majority being consumed in overcoming fault friction. Published by Elsevier Science Ltd.In the early 1960's, N.G.W. Cook, using an underground network of geophones, demonstrated that most Witwatersrand tremors are closely associated with deep level gold mining operations. He also showed that the energy released by the closure of the tabular stopes at depths of the order of 2 km was more than sufficient to account for the mining-induced earthquakes. I report here updated versions of these two results based on more modern underground data acquired in the Witwatersrand gold fields. Firstly, an extensive suite of in situ stress data indicate that the ambient state of crustal stress here is close to the failure state in the absence of mining even though the tectonic setting is thoroughly stable. Mining initially stabilizes the rock mass by reducing the pore fluid pressure from its initial hydrostatic state to nearly zero. The extensive mine excavations, as Cook showed, concentrate the deviatoric stresses, in localized regions of the abutments, back into a failure state resulting in seismicity. Secondly, there appears to be two distinct types of mining-induced earthquakes: the first is strongly coupled to the mining and involves shear failure plus a coseismic volume reduction; the second type is not evidently coupled to any particular mine face, shows purely deviatoric failure, and is presumably caused by more regional changes in the state of stress due to mining. Thirdly, energy budgets for mining induced earthquakes of both types indicate that, of the available released energy, only a few per cent is radiated by the seismic waves with the majority being consumed in overcoming fault friction.

Conference Paper↗

Effects of stochastically-simulated near-fault ground motions on soil liquefaction

The scarcity of historically recorded near-fault ground motions poses a challenge to systematically understanding the influence of near-fault effects on various types of seismic demands for engineering purposes. In particular, the current state of knowledge of the influence of ground-shaking intensity on soil liquefaction and its consequences does not specifically account for the effects of near-fault ground motion characteristics. In this study, the influence of near-fault ground motions on liquefaction triggering and lateral spreading are investigated using non-linear modeling of a hypothetical liquefiable soil column in the finite-element computational platform OpenSees subjected to simulated ground motion time series that represent strong earthquake shaking in the near field. The simulated ground motion time series and resulting datasets are based on a parametric stochastic model and are developed for a range of source and path parameters to represent a realistic variability of ground motion characteristics. Dependencies between ground motion intensity measures (IMs) and liquefaction demand parameters are investigated for near-fault pulse and nonpulse-like ground motion sets. Evolutionary IMs, such as cumulative absolute velocity (CAV) and the time-varying magnitude-adjusted peak ground acceleration (PGAM), are considered in developing liquefaction triggering probability density functions. Post-liquefaction triggering responses such as lateral spreading displacements are examined in relation to PGA M and CAV. The ground motion simulations are validated by comparing their liquefaction-capacity PGA M fragilities and post-triggering CAV vulnerability relationships to historical records from the 1994 Northridge earthquake in California, USA. Finally, a path forward for future studies that includes finding systematic differences in the IM-liquefaction demand relationships between near-fault and far-field stochastic ground motion sets is outlined.

Conference Paper↗

An evaluation of the strong ground motion recorded during the May 1, 2003 Bingol Turkey, earthquake

An important record of ground motion from a M 6.4 earthquake occurring on May 1, 2003, at epicentral and fault distances of about 12 and 9 km, respectively, was obtained at a station near the city of Bing&ouml;l, Turkey. The maximum peak ground values of 0.55g and 36 cm/s are among the largest ground-motion amplitudes recorded in Turkey. From simulations and comparisons with ground motions from other earthquakes of comparable magnitude, we conclude that the ground motion over a range of frequencies is unusually high. Site response may be responsible for the elevated ground motion, as suggested from analysis of numerous aftershock recordings from the same station. The mainshock motions have some interesting seismological features, including ramps between the P- and S-wave that are probably due to near- and intermediate-field elastic motions and strong polarisation oriented at about 39 degrees to the fault (and therefore not in the fault-normal direction). Simulations of motions from an extended rupture explain these features. The N10E component shows a high-amplitude spectral acceleration at a period of 0.15 seconds resulting in a site specific design spectrum that significantly overestimates the actual strength and displacement demands of the record. The pulse signal in the N10E component affects the inelastic spectral displacement and increases the inelastic displacement demand with respect to elastic demand for very long periods.

Journal of Earthquake Engineering↗

The 2011 Mw 7.1 Van (Eastern Turkey) earthquake

We use interferometric synthetic aperture radar (InSAR), body wave seismology, satellite imagery, and field observations to constrain the fault parameters of the M w 7.1 2011 Van (Eastern Turkey) reverse-slip earthquake, in the Turkish-Iranian plateau. Distributed slip models from elastic dislocation modeling of the InSAR surface displacements from ENVISAT and COSMO-SkyMed interferograms indicate up to 9 m of reverse and oblique slip on a pair of en echelon NW 40 &deg;–54 &deg; dipping fault planes which have surface extensions projecting to just 10 km north of the city of Van. The slip remained buried and is relatively deep, with a centroid depth of 14 km, and the rupture reaching only within 8–9 km of the surface, consistent with the lack of significant ground rupture. The up-dip extension of this modeled WSW striking fault plane coincides with field observations of weak ground deformation seen on the western of the two fault segments and has a dip consistent with that seen at the surface in fault gouge exposed in Quaternary sediments. No significant coseismic slip is found in the upper 8 km of the crust above the main slip patches, except for a small region on the eastern segment potentially resulting from the M w 5.9 aftershock on the same day. We perform extensive resolution tests on the data to confirm the robustness of the observed slip deficit in the shallow crust. We resolve a steep gradient in displacement at the point where the planes of the two fault segments ends are inferred to abut at depth, possibly exerting some structural control on rupture extent.

Journal of Geophysical Research↗

Triggered earthquakes and deep well activities

Earthquakes can be triggered by any significant perturbation of the hydrologic regime. In areas where potentially active faults are already close to failure, the increased pore pressure resulting from fluid injection, or, alternatively, the massive extraction of fluid or gas, can induce sufficient stress and/or strain changes that, with time, can lead to sudden catastrophic failure in a major earthquake. Injection-induced earthquakes typically result from the reduction in frictional strength along preexisting, nearby faults caused by the increased formation fluid pressure. Earthquakes associated with production appear to respond to more complex mechanisms of subsidence, crustal unloading, and poroelastic changes in response to applied strains induced by the massive withdrawal of subsurface material. As each of these different types of triggered events can occur up to several years after well activities have begun (or even several years after all well activities have stopped), this suggests that the actual triggering process may be a very complex combination of effects, particularly if both fluid extraction and injection have taken place locally. To date, more than thirty cases of earthquakes triggered by well activities can be documented throughout the United States and Canada. Based on these case histories, it is evident that, owing to preexisting stress conditions in the upper crust, certain areas tend to have higher probabilities of exhibiting such induced seismicity. ?? 1992 Birkha??user Verlag.

Pure and Applied Geophysics PAGEOPH↗

Earthquake stress drops, ambient tectonic stresses and stresses that drive plate motions

A variety of geophysical observations suggests that the upper portion of the lithosphere, herein referred to as the elastic plate, has long-term material properties and frictional strength significantly greater than the lower lithosphere. If the average frictional stress along the non-ridge margin of the elastic plate is of the order of a kilobar, as suggested by the many observations of the frictional strength of rocks at mid-crustal conditions of pressure and temperature, the only viable mechanism for driving the motion of the elastic plate is a basal shear stress of several tens of bars. Kilobars of tectonic stress are then an ambient, steady condition of the earth's crust and uppermost mantle. The approximate equality of the basal shear stress and the average crustal earthquake stress drop, the localization of strain release for major plate margin earthquakes, and the rough equivalence of plate margin slip rates and gross plate motion rates suggest that the stress drops of major plate margin earthquakes are controlled by the elastic release of the basal shear stress in the vicinity of the plate margin, despite the existence of kilobars of tectonic stress existing across vertical planes parallel to the plate margin. If the stress differences available to be released at the time of faulting are distributed in a random, white fasbion with a mean-square value determined by the average earthquake stress drop, the frequency of occurrence of constant stress drop earthquakes will be proportional to reciprocal faulting area, in accordance with empirically known frequency of occurrence statistics. ?? 1977 Birkha??user Verlag.

Pure and Applied Geophysics PAGEOPH↗

Quantitative study of tectonic geomorphology along Haiyuan fault based on airborne LiDAR

High-precision and high-resolution topography are the fundamental data for active fault research. Light detection and ranging (LiDAR) presents a new approach to build detailed digital elevation models effectively. We take the Haiyuan fault in Gansu Province as an example of how LiDAR data may be used to improve the study of active faults and the risk assessment of related hazards. In the eastern segment of the Haiyuan fault, the Shaomayin site has been comprehensively investigated in previous research because of its exemplary tectonic topographic features. Based on unprecedented LiDAR data, the horizontal and vertical coseismic offsets at the Shaomayin site are described. The measured horizontal value is about 8.6 m, and the vertical value is about 0.8 m. Using prior dating ages sampled from the same location, we estimate the horizontal slip rate as 4.0 &plusmn; 1.0 mm/a with high confidence and define that the lower bound of the vertical slip rate is 0.4 &plusmn; 0.1 mm/a since the Holocene. LiDAR data can repeat the measurements of field work on quantifying offsets of tectonic landform features quite well. The offset landforms are visualized on an office computer workstation easily, and specialized software may be used to obtain displacement quantitatively. By combining precious chronological results, the fundamental link between fault activity and large earthquakes is better recognized, as well as the potential risk for future earthquake hazards.

Gansu↗

Gravitational potential as a source of earthquake energy

Some degree of tectonic stress within the earth originates from gravity acting upon density structures. The work performed by this "gravitational tectonics stress" must have formerly existed as gravitational potential energy contained in the stress-causing density structure. According to the elastic rebound theory (Reid, 1910), the energy of earthquakes comes from an elastic strain field built up by fairly continuous elastic deformation in the period between events. For earthquakes resulting from gravitational tectonic stress, the elastic rebound theory requires the transfer of energy from the gravitational potential of the density structures into an elastic strain field prior to the event. An alternate theory involves partial gravitational collapse of the stress-causing density structures. The earthquake energy comes directly from a net decrease in gravitational potential energy. The gravitational potential energy released at the time of the earthquake is split between the energy released by the earthquake, including work done in the fault zone and an increase in stored elastic strain energy. The stress associated with this elastic strain field should oppose further fault slip.

Alaska↗

Ring fault creep drives volcano-tectonic seismicity during caldera collapse of Kīlauea in 2018

Basaltic caldera collapses are episodic, producing very-long-period (VLP) earthquakes up to M w 5.4, with prolific inter-collapse (between collapses) volcano-tectonic (VT) seismicity. During the 2018 caldera collapse of Kīlauea Volcano, VT seismicity ceased following each collapse, and then accelerated to a quasi-steady rate prior to the next collapse, marking a temporal pattern distinct from typical foreshock/aftershock sequences. There is currently no consensus on the mechanism(s) that generates the VT seismicity. Here we demonstrate that inter-collapse ring fault creep, induced by chamber depressurization, was the main driver of VT seismicity at Kīlauea in 2018. This is evidenced by: 1) the correlation between cumulative number of VT events and GNSS-derived ring fault creep; 2) agreement between repeating earthquake and GNSS derived creep rates; and 3) consistency between the time dependence of mechanically modeled, creep-driven seismicity and observations. We further show that, ring fault creep can be explained by velocity strengthening friction alone or in conjunction with viscous shear zone rheology. The simultaneous occurrence of creep and seismicity highlights the spatially heterogeneous velocity weakening/strengthening friction on the ring fault. If the VT seismicity-creep correlation can be replicated at other basaltic volcanoes, it would demonstrate that VT seismicity can be used as a proxy for ring fault creep in the absence of GNSS measurements on subsiding caldera block(s).

Hawaii↗

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↗