Large earthquakes initiate chains of surface processes that last much longer than the brief moments of strong shaking. Most moderate‐ and large‐magnitude earthquakes trigger landslides, ranging from small failures in the soil cover to massive, devastating rock avalanches. Some landslides dam rivers and impound lakes, which can collapse days to centuries later, and flood mountain valleys for hundreds of kilometers downstream. Landslide deposits on slopes can remobilize during heavy rainfall and evolve into debris flows. Cracks and fractures can form and widen on mountain crests and flanks, promoting increased frequency of landslides that lasts for decades. More gradual impacts involve the flushing of excess debris downstream by rivers, which can generate bank erosion and floodplain accretion as well as channel avulsions that affect flooding frequency, settlements, ecosystems, and infrastructure. Ultimately, earthquake sequences and their geomorphic consequences alter mountain landscapes over both human and geologic time scales. Two recent events have attracted intense research into earthquake‐induced landslides and their consequences: the magnitude M 7.6 Chi‐Chi, Taiwan earthquake of 1999, and the M 7.9 Wenchuan, China earthquake of 2008. Using data and insights from these and several other earthquakes, we analyze how such events initiate processes that change mountain landscapes, highlight research gaps, and suggest pathways toward a more complete understanding of the seismic effects on the Earth's surface.

No abstract

We propose a physical model for the high‐frequency (>1 Hz) spectral distribution of seismic power generated by debris flows. The modeled debris flow is assumed to have four regions where the impact rate and impulses are controlled by different mechanisms: the flow body, a coarser‐grained snout, a snout lip where particles fall from the snout on the bed, and a dilute front composed of saltating particles. We calculate the seismic power produced by this impact model in two end‐member scenarios, a thin‐flow and thick‐flow limit, which assume that the ratio of grain sizes to flow thicknesses are either near unity or much less than unity. The thin‐flow limit is more appropriate for boulder‐rich flows that are most likely to generate large seismic signals. As a flow passes a seismic station, the rise phase of the seismic amplitude is generated primarily by the snout while the decay phase is generated first by the snout and then the main flow body. The lip and saltating front generate a negligible seismic signal. When ground properties are known, seismic power depends most strongly on both particle diameter and average flow speed cubed, and also depends on length and width of the flow. The effective particle diameter for producing seismic power is substantially higher than the median grain size and close to the 73rd percentile for a realistic grain size distribution. We discuss how the model can be used to estimate effective particle diameter and average flow speed from an integrated measure of seismic power. © 2019 The Authors. Earth Surface Processes and Landforms Published by John Wiley & Sons Ltd. © 2019 The Authors. Earth Surface Processes and Landforms Published by John Wiley & Sons Ltd.

Earthquake‐induced landslide (EQIL) inventories are essential tools to extend our knowledge of the relationship between earthquakes and the landslides they can trigger. Regrettably, such inventories are difficult to generate and therefore scarce, and the available ones differ in terms of their quality and level of completeness. Moreover, access to existing EQIL inventories is currently difficult because there is no centralized database. To address these issues, we compiled EQIL inventories from around the globe based on an extensive literature study. The database contains information on 363 landslide‐triggering earthquakes and includes 66 digital landslide inventories. To make these data openly available, we created a repository to host the digital inventories that we have permission to redistribute through the U.S. Geological Survey ScienceBase platform. It can grow over time as more authors contribute their inventories. We analyze the distribution of EQIL events by time period and location, more specifically breaking down the distribution by continent, country, and mountain region. Additionally, we analyze frequency distributions of EQIL characteristics, such as the approximate area affected by landslides, total number of landslides, maximum distance from fault rupture zone, and distance from epicenter when the fault plane location is unknown. For the available digital EQIL inventories, we examine the underlying characteristics of landslide size, topographic slope, roughness, local relief, distance to streams, peak ground acceleration, peak ground velocity, and Modified Mercalli Intensity. Also, we present an evaluation system to help users assess the suitability of the available inventories for different types of EQIL studies and model development.

[1] Seismic methods can substantially improve the characterization of the dynamics of large and rapid landslides. Such landslides often generate strong long-period seismic waves due to the large-scale acceleration of the entire landslide mass, which, according to theory, can be approximated as a single-force mechanism at long wavelengths. I apply this theory and invert the long-period seismic waves generated by the 48.5 Mm 3 August 2010 Mount Meager rockslide-debris flow in British Columbia. Using data from five broadband seismic stations 70 to 276 km from the source, I obtain a time series of forces the landslide exerted on the Earth, with peak forces of 1.0 × 10 11 N. The direction and amplitude of the forces can be used to determine the timing and occurrence of events and subevents. Using this result, in combination with other field and geospatial evidence, I calculate an average horizontal acceleration of the rockslide of 0.39 m/s 2 and an average apparent coefficient of basal friction of 0.38 ± 0.02, which suggests elevated basal fluid pressures. The direction and timing of the strongest forces are consistent with the centripetal acceleration of the debris flow around corners in its path. I use this correlation to estimate speeds, which peak at 92 m/s. This study demonstrates that the time series recording of forces exerted by a large and rapid landslide derived remotely from seismic records can be used to tie post-slide evidence to what actually occurred during the event and can serve to validate numerical models and theoretical methods.

Earthquake‐triggered landslides are a significant hazard in seismically active regions, but our ability to assess the hazard they pose in near‐real‐time is limited. In this study, we present a new globally applicable model for seismically induced landslides based on the most comprehensive global data set available; we use 23 landslide inventories that span a range of earthquake magnitudes and climatic and tectonic settings. We use logistic regression to relate the presence and distribution of earthquake‐triggered landslides with spatially distributed estimates of ground shaking, topographic slope, lithology, land cover type, and a topographic index designed to estimate variability in soil wetness to provide an empirical model of landslide distribution. We tested over 100 combinations of independent predictor variables to find the best fitting model, using a diverse set of statistical tests. Blind validation tests show that the model accurately estimates the distribution of available landslide inventories. The results indicate that the model is reliable and stable, with high balanced accuracy (correctly versus incorrectly classified pixels) for the majority of test events. A cross‐validation analysis shows high balanced accuracy for a majority of events as well. By combining near‐real‐time estimates of ground shaking with globally available landslide susceptibility data, this model provides a tool to estimate the distribution of coseismic landslide hazard within minutes of the occurrence of any earthquake worldwide for which a U.S. Geological Survey ShakeMap is available.

We focus on the 6 August 2010 Mount Meager landslide that occurred in Southwest British Columbia, Canada. This 48.5 Mm3 rockslide that rapidly changed into a debris flow was recorded by over 25 broadband seismic stations. We showed that the waveform inversion of the seismic signal making it possible to calculate the time history of the force applied by the landslide to the ground is very robust and stable, even when using only data from a single station. By comparing this force with the force calculated through numerical modeling of the landslide, we are able to support the interpretation of seismic data made using a simple block model. However, our study gives different values of the friction coefficients involved and more details about the volumes and orientation of the subevents and the flow trajectory and velocity. Our sensitivity analysis shows that the characteristics of the released mass and the friction coefficients all contribute to the amplitude and the phase of the force. Despite this complexity, our study makes it possible to discriminate the best values of all these parameters. Our results suggest that comparing simulated and inverted forces helps to identify appropriate rheological laws for natural flows. We also show that except for the initial collapse, peaks in the low‐frequency force related to bends and runup over topography changes are associated with high‐frequency generation, possibly due to an increased agitation of the granular material involved.

Shortly before the beginning of the 2017–2018 winter rainy season, one of the largest fires in California (USA) history (Thomas fire) substantially increased the susceptibility of steep slopes in Santa Barbara and Ventura Counties to debris flows. On 9 January 2018, before the fire was fully contained, an intense burst of rain fell on the portion of the burn area above Montecito, California. The rainfall and associated runoff triggered a series of debris flows that mobilized ∼680,000 m3 of sediment (including boulders >6 m in diameter) at velocities up to 4 m/s down coalescing urbanized alluvial fans. The resulting destruction (including 23 fatalities, at least 167 injuries, and 408 damaged homes) underscores the need for improved understanding of debris-flow runout in the built environment, and the need for a comprehensive framework to assess the potential loss from debris flows following wildfire. We present observations of the inundation, debris-flow dynamics, and damage from the event. The data include field measurements of flow depth and deposit characteristics made within the first 12 days after the event (before ephemeral features of the deposits were lost to recovery operations); an inventory of building damage; estimates of flow velocity; information on flow timing; soil-hydrologic properties; and post-event imagery and lidar. Together, these data provide rare spatial and dynamic constraints for testing debris-flow runout models, which are needed for advancing post-fire debris-flow hazard assessments. Our analysis also outlines a framework for translating the results of these models into estimates of economic loss based on an adaptation of the U.S. Federal Emergency Management Agency’s Hazus model for tsunamis.