skip to main content


This content will become publicly available on November 13, 2024

Title: Coevolution of Rock Slope Instability Damage and Resonance Frequencies From Distinct‐Element Modeling
Abstract

Ambient vibration measurements can detect resonance frequency changes related to rock slope instability damage or boundary condition changes during progressive failure. However, the impact of slope kinematics on resonance changes and the expected form and sensitivity of frequency evolution during destabilization require clarification to improve the implementation of this technique across diverse settings. Since instrumented rock slope failures are rare, numerical modeling is needed to study the anticipated spectral response from in situ monitoring. We used 2D distinct‐element modeling to evaluate the sensitivity and evolution of rock slope resonance behavior for slab toppling, flexural toppling, and planar sliding instabilities during progressive failure. Model simulations revealed that fundamental resonance frequency decreases between 20% and 60% with changes correlated with increasing length of open joints. Changes to higher‐order frequencies associated with landslide sub‐volumes were also detectable for cases with multiple fracture networks. Resonance behavior was most pronounced for failures dominated by steeply dipping open tension cracks, that is, flexural and slab toppling. Additionally, amplification patterns across the slope varied for the flexural toppling and sliding cases, providing potential new information with which to characterize landslide failure mechanisms using ambient vibration array measurements. Our results demonstrate landslide characteristics well‐suited for in situ ambient resonance monitoring and provide new data describing the anticipated changes in resonance frequencies during progressive rock slope failure.

 
more » « less
Award ID(s):
2150896
NSF-PAR ID:
10477501
Author(s) / Creator(s):
 ;  
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Journal of Geophysical Research: Earth Surface
Volume:
128
Issue:
11
ISSN:
2169-9003
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Seismic resonance and surface displacement measurements can be implemented in tandem to improve landslide characterization and progressive failure monitoring. Crack aperture data are frequently used in rock slope stability monitoring and often exhibit recognizable trends prior to failure, such as accelerated crack opening. Alternatively, ambient resonance data offer multiple parameters including modal frequencies, damping, and polarization that can be monitored alongside crack aperture and may respond differently to environmental forcings and complex failure evolution. We analyzed data from continuous ambient vibration monitoring and concomitant crack aperture measurements at the Courthouse Mesa instability, a large toppling sandstone slab in Utah, USA. Three years of data revealed crack aperture increases of 2–4 mm/year with no clearly detectable irreversible changes in modal parameters, including frequency. Annually, frequency and displacement varied by 29% and 19% of the mean, respectively, with average and maximum daily frequency fluctuations of 6.5% and 16%, respectively. These reversible cyclic changes were primarily temperature-driven, but annually, frequency was in-phase with temperature whereas crack aperture lagged temperature changes by ∼37 days. Polarization and damping also varied seasonally but were less strongly correlated with temperature. Conceptual 3D finite element modeling demonstrated consistent frequency decreases associated with crack propagation but variable changes in crack aperture measured at a single point; i.e., crack propagation did not always result in increased crack opening but always generated a resonance frequency decrease. Taken together, our data suggest a possible thermal wedging-ratcheting mechanism at the Courthouse Mesa instability, where annual thermoelastic crack closure is impeded by debris infill but the absence of downward crack propagation during the monitoring period is evidenced by no permanent resonance frequency changes. Our study demonstrates that combined seismic resonance and crack aperture data provide an improved description of rock slope instability behavior, supporting refined characterization and monitoring of changes accompanying progressive failure.

     
    more » « less
  2. Abstract

    Accurate assessments of the internal structure and boundary conditions of unstable rock slopes are imperative for evaluating landslide hazard scenarios. However, instability characterization at depth remains challenging and is often limited by costly or invasive subsurface investigations. Here, we develop a new approach coupling array‐based ambient vibration modal analysis and numerical modeling to improve structural characterization of rock slope instabilities at depth. We used ambient noise cross‐correlation on 4 hr of seismic data recorded by an array of 30 nodal geophones at a 500‐m‐long toppling rock slab in Utah, USA to identify modal frequencies between 0.8 and 3.5 Hz and derive modal displacements. We show that transverse and longitudinal bending modes span the length of the instability, indicating an interconnected slab. Statistical comparison of field results with outputs from >1,000 finite element models with varying boundary conditions showed that the instability depth varies between 40–70 and 10–20 m in the middle and lateral regions, respectively. Our approach yields new information on the structural conditions of rock cliff and column instabilities at depth, which is not easily obtained by other means but is imperative for change detection monitoring and improved hazard assessments.

     
    more » « less
  3. Abstract

    Continuous ambient seismic monitoring is becoming common for structural health assessment of geologic features. However, the ability to detect or predict permanent mechanical change associated with damage requires detailed understanding of reversible drifts in resonance attributes associated with changing meteorological conditions. Here, we analyze the response of 17 sandstone rock arches to changing meteorological conditions during extended vibration‐based monitoring, with a focus on 1 arch in Utah which was continuously monitored for 15 months. Our results show that variations in resonance are correlated with temperature on daily and yearly time scales, but that the temperature sensitivity of frequency changes are variable at different sites and resonance modes, generally ranging from 0.5% to 6% per 10°C. Numerical modeling suggests the primary mechanism governing these frequency drifts is stress stiffening, where confined thermal dilation induces bulk stress changes in the low‐stress, nonlinear elastic regime via fracture closure. Secondary mechanisms affecting resonance attributes are driven by moisture, including formation of shallow pore ice, which can generate sharp frequency changes of up to 17% per −10°C, and moisture‐induced softening. Daily lag times of several hours between temperature and frequency extrema provide constraints on the rock volumes affected by these mechanisms, indicating modal attributes are sensitive to thermally driven changes occurring in the outermost centimeters of the structure. Through improved understanding of the reversible variations of rock mass resonance, our results aid future assessment of irreversible frequency changes associated with damage, and thus structural health monitoring of fragile geologic features.

     
    more » « less
  4. Abstract We deployed a network of 68 three-component geophones on the slow-moving Two Towers earthflow in northern California. We compute horizontal-to-vertical spectral ratios (HVSRs) from the ambient seismic field. The HVSRs have two prominent peaks, one near 1.23 Hz and another between 4 and 8 Hz at most stations. The 1.23 Hz resonance is a property of the background noise field and may be due to a velocity contrast at a few hundred meters depth. We interpret the higher frequency peaks as being related to slide deposits and invert the spectral ratios for shallow velocity structure using in situ thickness measurements as a priori constraints on the inversion. The thickness of the shallowest, low-velocity layer is systematically larger than landslide thicknesses inferred from inclinometer data acquired since 2013. Given constraints from field observations and boreholes, the inversion may reflect the thickness of deposits of an older slide that is larger in spatial extent and depth than the currently active slide. Because the HVSR peaks measured at Two Towers are caused by shallow slide deposits and represent frequencies that will experience amplification during earthquakes, the depth of the actively sliding mass may be less relevant for assessing potential slide volume and associated hazard than the thicknesses determined by our inversions. More generally, our results underscore the utility of combining both geotechnical measurements and subsurface imaging for landslide characterization and hazard assessment. 
    more » « less
  5. Abstract

    The effects of strong ground shaking on hillslope stability can persist for many years after a large earthquake, leading to an increase in the rates of post earthquake land sliding. The factors that control the rate of post‐earthquake land sliding are poorly constrained, hindering our ability to reliably forecast how landscapes and landslide hazards and risk evolve. To address this, we use a unique data set comprising high‐resolution terrestrial laser scans and airborne lidar captured during and after the 2010–2011 Canterbury Earthquake Sequence, New Zealand. This earthquake sequence triggered thousands of rock falls, and rock and debris avalanches (collectively referred to as “rockfall”), resulting in loss‐of‐life and damage to residential dwellings, commercial buildings and other infrastructure in the Port Hills of Christchurch, New Zealand. This unique data set spans 5 years and includes five significant earthquakes. We used these data to (a) quantify the regional‐scale “rockfall” rates in response to these earthquakes and the postearthquake decay in rockfall rates with time; and (b) investigate the site‐specific factors controlling the location of seismic and nonseismic rockfalls using frequency ratios and logistic regression techniques. We found that rockfall rates increased significantly in response to the initial earthquake that generated the strongest shaking in the sequence—The MW6.2 22 February 2011 event—Compared to the long‐term background rates derived from the dating of pre‐2010 talus piles at the toe of the slopes. Non seismic rockfall rates also increased immediately after the 22 February 2011 earthquake and decayed with time following a power‐law trend. About 50% of the decay back to the pre‐earthquake rockfall rates occurred within 1–5 years after the 22 February 2011 earthquake. Our results show that the short‐term decay in rockfall rates over time, after the initial earthquake, was attributed to the subsequent erosion of seismically damaged rock mass materials caused by environmental processes such as rain. For earthquake‐induced rockfall at the regional‐scale, the peak ground accelerations is the most significant variable in forecasting rockfall volume, followed by the relative height above the base of the slope. For both earthquake and non‐seismic conditions at the site‐specific scale, the probability of rockfall increases when the adjacent areas have failed previously, indicating that accrued damage preconditions localized areas of the slope for subsequent failure. Such preconditioning is a crucial factor driving subsequent rockfalls; that is, future rockfalls are likely to cluster near areas that failed in the past.

     
    more » « less