Shear surface undulations modulate clayey gouge strength and contribute to divergent landslide acceleration

William H. Schulz; Gonghui Wang; Yao Jiang; Brian D. Collins; Mark E. Reid; Mason Muir Einbund

doi:10.1016/j.enggeo.2025.108353

Home/Archive/USGS

USGS · 70271754

Shear surface undulations modulate clayey gouge strength and contribute to divergent landslide acceleration

William H. Schulz·Gonghui Wang·Yao Jiang·Brian D. Collins·Mark E. Reid·Mason Muir Einbund

View original source ↗Export metadata ↓

Abstract

Landslides display a spectrum of speeds for incompletely known reasons. Sliding occurs along slickensided undulatory shear surfaces within boundary shear gouge. Laboratory tests reveal that gouge shear strength generally decreases with finite cumulative displacement during relatively rapid failure and may increase or decrease with increasing shear rate; these behaviors can result in accelerating or decelerating landslide motion, which significantly affects consequent hazards. However, mechanisms responsible for such behaviors are poorly understood. We performed advanced ring shear strength testing that revealed such variable strength of a landslide near Oso, Washington, USA. We hypothesized that millimeter-scale undulations along shear surfaces caused the strength variability by imparting shear strength but while also modifying stresses that locally increase and decrease the typically considered particle-scale shear strength. We tested our hypotheses in the laboratory and with finite element soil deformation modeling. Lab results suggest that undulations contribute strength that decays with finite cumulative displacement. Modeling similarly reveals this, and that rapid shearing across undulations locally reduces effective normal stress by persistently elevating pore-water pressure and causing dilation. Consequent effects on strength differ by material with granular-rich, high-friction gouge losing strength and clay-rich, low-friction gouge losing little or gaining strength as shear rates increase. Ample testing by others reveals similar patterns. Hence, the propensity for gouge-controlled accelerating or decelerating failure may be estimated from simple index tests. Our findings on the effects of undulations reveal previously unknown mechanisms that may help to explain why some landslides reactivate catastrophically while others do not.

Explore related subjects

Engineering Geology Article

Source-reported bounding extent: 48.26490761362555° to 48.28444110597036° latitude; -121.94130726116578° to -121.90311819756626° longitude. This indicates report coverage, not an exact sampling location. View area on OpenStreetMap.

Keep this discovery

Explore connections, maps & timelines

BibTeX RIS

William H. Schulz, Gonghui Wang, Yao Jiang, Brian D. Collins, Mark E. Reid, Mason Muir Einbund. 2025. Shear surface undulations modulate clayey gouge strength and contribute to divergent landslide acceleration. https://doi.org/10.1016/j.enggeo.2025.108353

Cite the original work for its findings. Save a collection to share your selection of sources.

KEEP EXPLORING

Related USGS reports

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↗

USGS2026

From landslide susceptibility to risk assessment in the conterminous U.S.

Understanding the spatial distribution of landslide prone-areas and what consequences they may have is important for risk management and land-use planning. In the United States, although landslides occur in every state, a comprehensive landslide risk assessment is still missing. Existing efforts, such as the Federal Emergency Management Agency (FEMA)’s National Risk Index, rely on aggregated products and coarse cartographic units, limiting their geomorphological and practical accuracy. In this study, we present a methodological advance for landslide risk assessment across large areas with incomplete and sparse data. We apply our procedures to the conterminous United States by integrating geomorphologically meaningful partitions and spatial and temporal probability data-driven models. Landslide susceptibility is estimated using a Generalized Additive Mixed Model incorporating a bias capture/correction scheme to account for inventory inaccuracies (reference Area Under the Curve = 0.75). The exceedance probabilities of landslide occurrence are defined for three temporal scenarios (2, 5, and 10 year). Then, we explore the associated potential economic consequences for human settlements and agricultural areas. The findings indicate that the spatial variability of risk is primarily controlled by exposure rather than by susceptibility/hazard alone. The mean risk increases by ∼170% from the 2-year to the 10-year scenario. Beyond its quantitative outcomes, this study offers a blueprint for continental or sub-continental scale landslide risk assessments, demonstrating both the opportunities and current limitations.

Engineering Geology↗

USGS2025

Global patterns of coseismic landslide runout mobility differ from aseismic landslide trends

Coseismic landslides significantly contribute to human and economic losses during and immediately following earthquakes, yet very little data on the runout of such landslides exist. While well-established behavior of aseismic (e.g., hydrologically triggered) landslide runout mobility suggests strong correlation between landslide size and mobility, limited studies of coseismic landslide runout find conflicting mobility trends. We present a global dataset of runout lengths produced from a new automated method for estimating landslide runout, developed and validated using 1726 manually mapped landslides from five unique earthquakes. We then apply the automated runout tool to 23 global earthquake-induced landslide inventories, producing a compiled database of 73,665 measured and estimated runout lengths of coseismic landslides to assess mobility trends. We find a significant divergence between well-established aseismic mobility trends and that of coseismic landslides, with far greater scatter and more complex mobility patterns in earthquake-triggered landslides. As a function of landslide size, we observe global coseismic landslide mobility patterns are bilinear, becoming increasingly less mobile with increasing size above some threshold. This discordance between aseismic and coseismic landslide mobility may be a function of landslide type, kinematics, hydrology, and or setting that systematically differ between triggering mechanisms and should be explored in more depth to develop predictive models of these unique runout patterns. These results suggest hazard and risk models for coseismic landslides may significantly under-predict or over-predict impacts, depending on the size of triggered landslides.

Engineering Geology↗