Geology topics

Kate E. Allstadt

Publications and source records attributed to Kate E. Allstadt.

At least 19 recordsLinked to original sources

USGS2027

Ground failure resulting from earthquake ground motions

The ground failure resulting from an earthquake, at all scales, is ultimately controlled by the details of near-surface earthquake ground motions, the spatially distributed mechanical characteristics of surface materials and how the two interact. This chapter introduces the aspects of ground-motion seismology that are most relevant to the triggering of different types of ground failure, such as landslides and liquefaction. In particular, the chapter focuses on different sources of local modification (e.g. amplification) of earthquake ground motions in the shallow subsurface, known as site effects. The chapter ends with an overview of the controls on, general spatial patterns of, and approaches to regional modelling of the manifestation of ground failure that may be triggered by an earthquake.

Book chapter↗

USGS2026

Frictional weakening in the highly mobile 2025 Blatten rock and ice avalanche in Switzerland

Cascading slope failures are increasing as glaciers retreat and slopes adjust to a warming climate, yet reconstructing their full dynamics from precursor to deposition remains challenging. Here we show that the 28 May 2025 Blatten rock-ice avalanche, which mobilized about 9.3 million cubic metres, can be reconstructed from seismic, geomorphological, and geotechnical evidence and numerical modeling. The observations reveal the timing, source, volume, and friction of over two weeks of accelerating rockfalls and minor glacier collapses before the main failure, followed by four stages of avalanche motion. Only with drastically reduced friction can both depth-averaged and three-dimensional models reproduce the center-of-mass force history and the deposit distribution. These friction values reflect effective frictional weakening. Geotechnical observations suggest that transient excess pore pressures possibly contributed to this weakening, but are insufficient to rule out additional mechanisms such as ice frictional reduction, segregation, and fragmentation. This work provides a reproducible and robust framework for analyzing catastrophic avalanches under a changing climate.

Birch Glacier↗

USGS2026

Landslide volume estimation from seismic waveform features using a global catalog of seismogenic events

Seismic waves generated by large, rapid landslides encode information about the source and can be analyzed rapidly following an event. Even remote landslides can pose hazards to downstream communities, so rapid detection and characterization using existing seismic monitoring networks could be beneficial. In this study, we expand on past regionally limited work by presenting a globally applicable method for estimating landslide volume from seismic features that could be integrated into future landslide seismic monitoring frameworks. We train the model using multivariable linear regression and five seismic features derived from recordings of 129 landslide events of a range of styles and locations with independently estimated volumes. We present two preferred models, one that combines long‐period (LP) and high‐frequency (HF) features and one for use on smaller landslides without observed long‐period signals. We find that our best‐performing model, applicable to landslides larger than 100,000 m 3 with signals containing observable long‐period energy, requires only two features: LP (20–100 s) absolute maximum amplitude and HF (1–5 Hz) rise time (time between the signal start and the maximum envelope amplitude) and has an R 2 score of 0.79. This model predicted volumes within one order of magnitude for 55 out of 58 events. We find that this combination mitigates the trade‐off between mass and acceleration and other variations of landslide style that limit methods based on amplitude alone.

Seismological Research Letters↗

USGS2026

Rapid seismic and infrasound assessment of large landslides: A case study from Denali National Park and Preserve (Alaska)

Large, rapid landslides are a global hazard that can occur in remote, mountainous areas. Eyewitness reports of landslides and satellite imagery can often be limited or delayed, particularly during inclement weather. However, landslide-generated seismic and infrasound (low-frequency atmospheric sound) waves can be remotely detected in near real-time. This information can significantly expedite characterization and possible landslide response activities. Here, we highlight these capabilities using a > 4 million m³ ice–rock avalanche in Denali National Park and Preserve (Alaska). This event was detected via a landslide-specific seismic location and volume estimation algorithm deployed in Alaska, and — notably — by standard earthquake monitoring systems. Following rapid detection of this event, we combined its seismic and infrasound dataset with optical, synthetic aperture radar, and oblique aerial imagery, multitemporal digital elevation models, and a numerical flow model to reconstruct its failure timeline and dynamics. We apply array processing to infrasound signals traveling > 250 km and find that two precursory events occurred minutes prior to the main failure. We use long-period seismic signals to infer the force exerted by the landslide on the Earth and constrain the rheological parameters of our numerical flow simulation with this result and deposit morphology. The main failure produced a steeply-dipping impulsive initial downward force and reached speeds exceeding 60 m/s. This impulsive force generated relatively strong seismic body waves, which contributed to the earthquake system detection. This large, remote Alaska landslide underscores the key value of seismic and infrasound analysis for rapid landslide assessment and motivates efforts to further operationalize these approaches.

Alaska↗

USGS2026

Quantifying seismic properties of a river channel at Mount Rainier for use in debris flow monitoring and analysis

Theoretical models that relate debris flow properties to their seismic signature suggest seismic methods may be used to remotely characterize properties of these events. However, the complexity of debris flow sources and poorly constrained material properties in near‐surface environments limit our ability to determine important attributes of debris flows, such as volume and particle size distribution, from seismic records alone. In this study, we explore the sensitivity of debris flow seismic signals to subsurface seismic characteristics using an established debris flow seismicity model and subsurface properties derived from active‐source measurements in the Tahoma Creek stream channel at Mount Rainier, United States. Using refraction and multichannel analysis of surface waves, we estimate 1D primary and secondary wave velocity profiles to a depth of 11 m and calculate the frequency‐varying phase velocity ( ⁠⁠ v c ) and Rayleigh‐wave quality factor (Q R ⁠⁠ ) for frequencies between 9 and 50 Hz. We find that the Tahoma Creek stream channel has low v c ⁠ , varying with frequency between 226 and 434 m/s, and is highly attenuating, with Q R estimated below 13.3 at all analyzed frequencies. We model the seismic signal of a hypothetical debris flow in Tahoma Creek using our measured values and compare the results against the same model but varying v c and Q R over a range of frequency‐independent, single values commonly used in the literature when measurements are not available ( ⁠⁠ v c = 250-750 m/s, Q R = 3-33 ⁠ ). We find that the power spectral densities of the modeled debris flows vary by orders of magnitude within the subset of our test values, highlighting the benefits of measuring material properties when using modeled debris flow seismic signals for quantitative monitoring.

Washington↗

USGS2026

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↗

USGS2026

Sackung at Bald Eagle ridge, central Colorado: An updated interpretation of ridge-spreading movement, structures, and mechanisms from 50 years (1975–2025) of U.S. Geological Survey research

Slow gravitational failures of mountain peaks and ridges are poorly understood. Herein, we report on 50 years of studies at a slowly spreading castellate ridge in the Sawatch Range in central Colorado. The orientations of geomorphic-structural features indicate that the fractured Precambrian granitic rock underlying the ridge has extended and spread northwestward toward the formerly glacier-covered Busk Creek valley. Results from surveying, field-based geomorphic-structural mapping using lidar, rock mass quality measurements, a passive seismic survey, and satellite radar provide a major update to research started by U.S. Geological Survey researchers in the 1970s and 1980s. New insights include a recognition that the entire ridge has slowly moved by concurrent sliding along an inferred northwest dipping, compound basal-slip surface (or zone), and through the formation of multiple grabens by normal faulting and flexural toppling along sets of pre-existing fractures that dip moderately (∼45°) to the southeast. We were unable to distinguish the presence of a sudden and strong contrast in seismic velocity across the inferred slip surface. Movement during the 50-year study period has been episodic and gradually decreasing, in correspondence with decreasing cumulative annual precipitation and increasing mean annual air temperatures. The fastest moving area, just upslope from the glacier trimline, had an average horizontal velocity of 3 – 4 mm/yr. Evidence suggests that movement started as a paraglacial response mechanism, but because of the site’s proximity to the Rio Grande Rift, we cannot exclude earthquake shaking as a mechanism for initiation or enhancement of slope movement. An estimate of longer-term horizontal movement from the exposed basal-slip surface at the uphill side of the ridgetop graben is ∼1.1 mm/yr for the 13 – 14 ky post-glacial period. Broad implications of our work are that: (1) long-term measurements (decades or longer) of slope movement can add insights into how sackungen form and evolve through time; (2) the identification of thrust faults and toes in zones of compression near valley bottoms can be crucial for interpreting sackung failure mechanisms, and (3) the use of passive-seismic techniques to identify the depth to a slip surface may not be successful in granitic terrain dominated by planar fractures and subtle changes in rock-mass characteristics.

Colorado↗

USGS2025

Insights into widespread landsliding in southern Appalachia from Hurricane Helene

Between 23 and 27 September 2024, antecedent rain followed by Hurricane Helene produced one of the most damaging weather events in southern Appalachia history. The back-to-back storm events resulted in a maximum cumulative rainfall of 848 mm and hurricane-force wind gusts over 170 km/h in western North Carolina, eastern Tennessee, and southwestern Virginia. The resulting regional flooding, landslides, and tree blowdown caused over 100 fatalities, damaged or destroyed critical infrastructure and thousands of structures, and severed connectivity across the region. Over the next several weeks, a multi-agency landslide response produced a rapid hazard assessment and mapped 2217 landslides, 55% of which damaged infrastructure or property. Orographic uplift enhanced rainfall, resulting in concentrated landsliding along the ~250 km swath of the Blue Ridge escarpment in western North Carolina. Landslides initiated predominantly on windward-facing (southeast-facing) slopes, and localized clustering of initiation points indicated a strong influence of hillslope-scale meteorological and geomorphic factors. Many shallow landslides mobilized into larger, highly mobile, and damaging debris flows that graded into floods. Here, we put our preliminary observations in the context of historical storm-driven landslide events and open new avenues for investigating the nature and extent of landslides and their effects in southern Appalachia and similar environments.

Georgia, North Carolina, South Carolina, Tennessee↗

USGS2025

Updating regional‐scale geospatial liquefaction models with locally available geotechnical data

We present a method to update the geospatial liquefaction model used by the U.S. Geological Survey’s near‐real‐time ground failure product with subsurface geotechnical data. The geospatial model estimates liquefaction probability from peak ground velocity (via ShakeMap) and geospatial susceptibility proxies. In many regions, additional information relevant to constraining liquefaction likelihood is also available, including surface geology maps and subsurface geotechnical measurements. There is currently no mechanism to use these data in the ground failure product liquefaction model, even though these data could provide more precise constraints on spatial variations in the lithologic character of the soil (surface geology) and direct measurements of the subsurface mechanical properties that affect liquefaction occurrence and severity (geotechnical measurements). In this study, we develop a method to integrate these data with the geospatial model and assess how these data can improve regional‐scale predictions. We develop a Bayesian updating framework and apply it to the 1989 magnitude 6.9 Loma Prieta, California, earthquake, for which mapped observations are available to evaluate performance. We constrain the Bayesian framework with 373 Northern California cone penetration tests and liquefaction susceptibility classes based on the mapped surface geology. This Bayesian model incorporates geotechnical information into the geospatial model and more accurately predicts liquefaction occurrences than the geospatial model, while sacrificing less accuracy in terms of predicting the absence of liquefaction than the geotechnical model. In future applications, this approach could be adapted to update other geospatial models using locally available subsurface data.

California↗

USGS2025

Preliminary field report of landslide hazards following Hurricane Helene

Executive Summary This report reflects our knowledge regarding the widespread landslide activity associated with Hurricane Helene observed during the U.S. Geological Survey’s (USGS) mission assignment to North Carolina in October 2024. The material in this report was originally prepared for the Federal Emergency Management Agency under mission assignment DR-4827-NC. The data and commentary in this report are reflective of a report provided to the Federal Emergency Management Agency (FEMA) on October 18, 2024, as well as information provided in briefings at the Buncombe County Emergency Operations Center. The report has been modified for public dissemination. This assessment was based on systematic visual examination and mapping of landslide locations from aerial and satellite imagery, visual and photographic observations from low-level helicopter overflights and conversations with local landslide experts from the North Carolina Geological Survey and Appalachian Landslide Consultants PLLC, and more than 50 years of combined landslide hazard professional experience of the mission-assigned field team. No systematic field investigations were done by the USGS. While responding to the event, the USGS did not identify any landslides that posed an immediate major threat to recovery personnel in parts of nine counties in North Carolina (Avery, Buncombe, Henderson, McDowell, Mitchell, Polk, Rutherford, Watauga, and Yancey); however, threats from renewed landslide activity may remain heightened in localized areas for months or even years. Known areas of the most abundant landslide occurrence include Bat Cave, Lake Lure, Chimney Rock, Swannanoa, Black Mountain, Fairview, steep areas in Asheville, and the Blue Ridge Parkway. The USGS shared detailed locations of known landslides with the Emergency Operations Centers. The thousands of landslide scars on hillsides and landslide deposits on flatter ground may present some threat to recovery activities. Soil and rocks will continue to erode from newly exposed landslide scars and may pose a threat to people and infrastructure who are immediately nearby. In general, the steeper and taller the landslide scar, the greater the potential threat. This threat is heightened during periods of rainfall and increases with the duration and intensity of rainstorms. Very heavy rainfall, or repeated rainfall events during short periods, could also initiate new landslides on steep slopes. Excavation of landslide deposits, particularly excavation of those deposits directly adjacent to steep slopes, may also pose a threat to nearby people and equipment. An interagency collaborative mapping effort led by the USGS that informed this assessment identified 1,155 landslide locations by the October 2024 briefings, but that number increased to 2,217 in a final reviewed version of the locations published in January 2025. Locations were mapped from satellite imagery, fixed-wing and helicopter surveys, media and social media, and field reports in the 3 weeks following the passage of the remnants of Hurricane Helene. USGS products outlined in this report are publicly available and include geotagged photographs from aerial reconnaissance, hazard models, an interactive view of mapped landslide locations, and landslide safety and education resources.

North Carolina, South Carolina, Tennessee, Virgini↗

USGS2025

Utility of a swath laser rangefinder for characterizing mass movement flow depth and landslide initiation

Mass movements such as debris flows and landslides are some of the deadliest and most destructive natural hazards occurring mostly in alpine and volcanic settings. With ever-growing populations located downslope from known debris flow channels, early warning systems can help prevent loss of life. Geophysical and technological advances have improved monitoring and detection capabilities in recent years; however, they can often be cost prohibitive and resource intensive, making them less accessible to disadvantaged populations. We tested and validated a readily available and cost-effective two-dimensional swath laser rangefinder in a controlled experimental setting against two independent flow-depth lasers. The swath laser successfully recorded cross-sectional changes in flow depth from four debris flows and a water-only flood, in addition to geomorphic changes associated with landslide initiation. The results suggest that a swath laser could be integrated into systems for debris flow detection and characterization of mass movements in natural settings, thus improving the ability to monitor these hazards.

Landslides↗

USGS2024

Using the horizontal-to-vertical spectral ratio method to estimate thickness of the Barry Arm landslide, Prince William Sound, Alaska

Conducting detailed investigations of large landslides is difficult, especially in the subsurface, largely due to environmental factors such as steep slopes, difficult access, and numerous objective hazards. These factors have made it challenging to accurately estimate the depth to the failure surface of the Barry Arm landslide, a large (roughly 10 8 cubic meters), deep-seated bedrock landslide in Prince William Sound, Alaska, recognized in 2019. The landslide has exhibited accelerated movement in recent years and poses a potential tsunamigenic hazard if rapid failure occurs. Failure surface depth, equivalent to landslide thickness, is a necessary metric for landslide-volume calculations and associated tsunami wave models. In this report, we used seismic noise recorded by a seismometer located on the Barry Arm landslide in Alaska to calculate the horizontal-to-vertical spectral ratio (HVSR) to investigate the site fundamental frequency ( f 0 ) and depth of the failure surface. To ensure that observed peak frequencies in the spectral ratio were related to the underlying stratigraphy (and not caused by other noise sources like nearby glaciers, topographic resonance, weather, or human activities), we also calculated HVSRs using earthquake signals, HVSRs at other seismic stations within a 2.5-kilometer radius, and a standard spectral ratio between the landslide station and other sites. We observed multiple peaks in the landslide HVSR curves at 1.5 hertz (Hz), 4–5 Hz, and 7–11 Hz. The frequencies of these peaks were consistent at the landslide site through time and across methods and were dissimilar to those identified at other seismic stations in the area, making it unlikely the peaks were caused by local noise. Directional HVSRs calculated at 15-degree intervals showed amplification of the higher frequency peaks in the direction parallel to slip, indicating two-dimensional site effects. We used the distinct frequency peaks in the seismic record to develop a 4-layer conceptual model of the landslide wherein the top of the deepest layer represents the primary failure surface, or the boundary between damaged (mobile) and undamaged material. We inverted Rayleigh wave ellipticity curves within this 4-layer configuration with constraints on S-wave velocity and layer thickness based on analogous material properties identified in the literature. This was necessary absent any site-specific subsurface S-wave velocity data. The best-fitting models indicate a mean slope-normal depth to the failure surface of 188 (±9) meters (m), with additional stratigraphic boundaries at 4 and 20 m below ground surface, potentially representing layered motion. These results agree with and improve upon ranges estimated by previous studies and can support future modeling and assessment efforts at Barry Arm.

Alaska↗

USGS2024

Earthquake-triggered ground-failure inventory associated with the M7.1 2018 Southcentral Alaska earthquake

The 30 November 2018, magnitude (Mw) 7.1 earthquake in Southcentral Alaska triggered substantial landslides, liquefaction, and ground cracking throughout the region, resulting in widespread geotechnical damage to buildings and infrastructure. Despite a challenging reconnaissance and remote-sensing environment, we constructed a detailed digital inventory of ground failure associated with the event from several sources. Sources included information derived from remotely sensed data, and data compiled from literature, social media postings, and earthquake damage information compiled by local, state, and federal agencies. Each instance of ground failure within the inventory contains information on the location and type of observed ground failure, and the methods and data used to document the occurrence. Where high-quality data, such as LIDAR or satellite imagery, were available and showed the ground-failure instance clearly, the extent is mapped as a polygon or polyline. All other locations are mapped as points. There are a total of 886 ground-failure instances documented within the inventory (400 landslides, 286 liquefaction features, and 200 features unattributed to specific processes). A semi-quantitative confidence scheme is used to describe mapping certainty associated with each ground-failure feature. This inventory represents a relatively moderate ground-failure-triggering event that occurred in a subarctic environment. This data paper describes the content within the inventory, the inventory data collection procedures, and limitations of the data. Events of this type are not often documented in detail; thus, adding the inventory data to the US Geological Survey Open Repository of Earthquake-Triggered Ground-Failure Inventories further diversifies the datasets available to the scientific community to be used to better understand and model earthquake-triggered ground failure.

Alaska↗

USGS2024

The influence of anthropogenic regulation and evaporite dissolution on earthquake-triggered ground failure

Remote sensing observations of Searles Lake following the 2019 moment magnitude 7.1 Ridgecrest, California, earthquake reveal an area where surface ejecta is arranged in a repeating hexagonal pattern that is collocated with a solution-mining operation. By analyzing geologic and geotechnical data, here we show that the hexagonal surface ejecta is likely not a result of liquefaction. Instead, we propose dissolution cavity collapse (DCC) as an alternative driving mechanism. We support this theory with pre-event Interferometric Synthetic Aperture Radar data, which reveals differential subsidence patterns and the creation of subsurface void space. We also find that DCC is likely triggered at a lower shaking threshold than classical liquefaction. This and other unknown mechanisms can masquerade as liquefaction, introducing bias into liquefaction prediction models that rely on liquefaction inventories. This paper also highlights the opportunities and drawbacks of using remote sensing data to disentangle the complex factors that influence earthquake-triggered ground failure.

California↗

USGS2024

Slope Unit Maker (SUMak): An efficient and parameter-free algorithm for delineating slope units to improve landslide modeling

Slope units are terrain partitions bounded by drainage and divide lines. In landslide modeling, including susceptibility modeling and event-specific modeling of landslide occurrence, slope units provide several advantages over gridded units, such as better capturing terrain geometry, improved incorporation of geospatial landslide-occurrence data in different formats (e.g., point and polygon), and better accommodating the varying data accuracy and precision in landslide inventories. However, the use of slope units in regional ( > 100 km 2 ) landslide studies remains limited due, in part, to the large computational costs and/or poor reproducibility with current delineation methods. We introduce a computationally efficient algorithm for the parameter-free delineation of slope units that leverages tools from within TauDEM and GRASS, using an R interface. The algorithm uses geomorphic laws to define the appropriate scaling of the slope units representative of hillslope processes, avoiding the often ambiguous determination of slope unit size. We then demonstrate how slope units enable more robust regional-scale landslide susceptibility and event-specific landslide occurrence maps.

Natural Hazards and Earth Systems Sciences (NHESS)↗

USGS2023

The 2022 Chaos Canyon landslide in Colorado: Insights revealed by seismic analysis, field investigations, and remote sensing

An unusual, high-alpine, rapid debris slide originating in ice-rich debris occurred on June 28, 2022, at 16:33:16 MDT at the head of Chaos Canyon, a formerly glacier-covered valley in Rocky Mountain National Park, CO, USA. In this study, we integrate eyewitness videos and seismic records of the event with meteorological data, field observations, pre- and post-event satellite imagery, and uncrewed aircraft vehicle imagery to characterize the event and future hazards it may pose. Deformation of the eventual slide mass preceded rapid failure by decades, starting in the early to mid-2000s, accelerating in 2018 (the warmest year on record), and reaching ~ 20 m/year in 2021. The main event, which was preceded by smaller sliding episodes earlier that day, had a volume of ~ 2.1 million m 3 , reached peak velocities of about 5 m/s, slid on a surface up to 80 m deep, and moved up to ~ 245 m downslope in < 2 min. We observed blocks of frozen debris (permafrost) in the landslide deposits. Within ~ 2 weeks, these blocks had melted and became dry, conical debris mounds (molards). We hypothesize that the rapid slide was induced by gradually increasing long-term air temperatures that thawed ice and increased pore pressures. The presence and suspected influence of permafrost on the occurrence of this landslide indicate other slopes in the park, and other moderate-to-low latitude alpine regions may experience similar slope stability issues as temperatures continue to warm.

Colorado↗

USGS2023

Mapping landslide susceptibility over large regions with limited data

Landslide susceptibility maps indicate the spatial distribution of landslide likelihood. Modeling susceptibility over large or diverse terrains remains a challenge due to the sparsity of landslide data (mapped extent of known landslides) and the variability in triggering conditions. Several different data sampling strategies of landslide locations used to train a susceptibility model are used to mitigate this challenge. However, to our knowledge, no study has systematically evaluated how different sampling strategies alter a model's predictor effects (i.e., how a predictor value influences the susceptibility output) critical to explaining differences in model outputs. Here, we introduce a statistical framework that examines the variation in predictor effects and the model accuracy (measured using receiver operator characteristics) to highlight why certain sampling strategies are more effective than others. Specifically, we apply our framework to an array of logistic regression models trained on landslide inventories collected at sub-regional scales over four terrains across the United States. Results show significant variations in predictor effects depending on the inventory used to train the models. The inconsistent predictor effects cause low accuracies when testing models on inventories outside the domain of the training data. Grouping test and training sets according to physiographic and ecological characteristics, which are thought to share similar triggering mechanisms, does not improve model accuracy. We also show that using limited landslide data distributed uniformly over the entire modeling domain is better than using dense but spatially isolated data to train a model for applications over large regions.

California, Kentucky, North Carolina, West Virgini↗

USGS2023

Preliminary documentation of coseismic ground failure triggered by the February 6, 2023 Türkiye earthquake sequence

The devastating Kahramanmaraş earthquake sequence occurred on February 6, 2023. Two main events, M w 7.8 and M w 7.5 occurred 9 hours apart, affected 11 cities in Turkey, and subjected an area of ∼90,000 km 2 to shaking levels known to trigger landslides (peak ground acceleration > 0.08 g). Extensive landsliding was expected given the hilly terrain affected by this significant ground shaking—about 15% of the topography is steeper than 20°—but was not initially apparent in early satellite imagery, mostly because of obscuring snow that fell just after the earthquakes. However, after a more detailed investigation using high-resolution satellite images, aerial photos, and a field survey, we confirmed that this earthquake sequence did, indeed, trigger numerous landslides. In this study, we present those findings and provide a preliminary characterization of the spatial distribution, general characteristics, and dominant types of landslides and hillslope deformation triggered by the earthquake sequence. We mapped 3673 coseismic landslides, mostly concentrated in the northern half of the impacted area. Rock falls are the most abundant landslide type, but bedrock rotational landslides, translational slides and lateral spreads are also numerous. Surface rupture through mountainous terrain caused several large, and in some cases fatal, landslides. Incipient landslides and ground cracks are also widespread, especially in the north. Lithology, spatial variability of ground shaking, and topographic relief appear to be the main variables controlling the spatial distribution of coseismic landslides. There are few detailed studies of earthquake-triggered landslides in arid and semi-arid regions such as this one, nor for such complex earthquake sequences. Therefore, this contribution provides valuable information for future hazard and modeling efforts in arid and semi-arid regions.

Engineering Geology↗