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.
Large earthquakes can construct mountainous topography by inducing rock uplift but also erode mountains by causing landslides. Observations following the 2008 Wenchuan earthquake show that landslide volumes in some cases match seismically induced uplift, raising questions about how the actions of individual earthquakes accumulate to build topography. Here we model the two‐dimensional surface displacement field generated over a full earthquake cycle accounting for coseismic deformation, postseismic relaxation, landslide erosion, and erosion‐induced isostatic compensation. We explore the related volume balance across different seismotectonic and topographic conditions and revisit the Wenchuan case in this context. The ratio (Ω) between landslide erosion and uplift is most sensitive to parameters determining landslide volumes (particularly earthquake magnitude
- NSF-PAR ID:
- 10371459
- Publisher / Repository:
- DOI PREFIX: 10.1029
- Date Published:
- Journal Name:
- Journal of Geophysical Research: Solid Earth
- Volume:
- 124
- Issue:
- 5
- ISSN:
- 2169-9313
- Page Range / eLocation ID:
- p. 5101-5133
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract -
more » « less
Abstract An important manifestation of the background or coseismic deformation surrounding megathrust earthquakes is the offshore microseismicity, which are difficult to be directly detected by land‐based seismic instruments. Here, we improve the capability of detecting offshore events by combining two popular techniques: backprojection (BP) imaging and matched‐filter (MF) detection. The BP method is effective in retrieving offshore seismicity (
Mw > 4.5Mw Mw > 4b values (≥ 1.2) are found close to large aftershocks, possibly indicating localized pockets of small differential stresses. At several locations close to the trench,p values (0.93–1.11) are higher than those in the inland area (0.64–0.85). This may be due to the larger coseismic slip and hence larger stress drop of the outer‐rise normal‐faulting events compared to the deeper thrust‐faulting events. -
more » « less
Abstract The heterogeneous seafloor topography of the Nazca Plate as it enters the Ecuador subduction zone provides an opportunity to document the influence of seafloor roughness on slip behavior and megathrust rupture. The 2016 Mw7.8 Pedernales Ecuador earthquake was followed by a rich and active postseismic sequence. An internationally coordinated rapid response effort installed a temporary seismic network to densify coastal stations of the permanent Ecuadorian national seismic network. A combination of 82 onshore short and intermediate period and broadband seismic stations and six ocean bottom seismometers recorded the postseismic Pedernales sequence for over a year after the mainshock. A robust earthquake catalog combined with calibrated relocations for a subset of magnitude ≥4 earthquakes shows pronounced spatial and temporal clustering. A range of slip behavior accommodates postseismic deformation including earthquakes, slow slip events, and earthquake swarms. Models of plate coupling and the consistency of earthquake clustering and slip behavior through multiple seismic cycles reveal a segmented subduction zone primarily controlled by subducted seafloor topography, accreted terranes, and inherited structure. The 2016 Pedernales mainshock triggered moderate to strong earthquakes (5 ≤ M ≤ 7) and earthquake swarms north of the mainshock rupture close to the epicenter of the 1906 Mw8.8 earthquake and in the segment of the subduction zone that ruptured in 1958 in a Mw7.7 earthquake.
-
more » « less
Abstract Post‐seismic debris flows are an important hazard following large earthquakes, propagating destruction downstream from hillslopes where coseismic landslides occur and extending damage for years after shaking stops. Data sets of post‐seismic debris flows are necessary to predict initiation and runout characteristics but are presently scarce. We used satellite imagery supplemented by field observations to compile an inventory of >1,000 debris flows associated with the 2015 Gorkha Earthquake in Nepal. We identified two distinct debris flow types: (1) Material from a coseismic landslide was remobilized in a steep channel during a later monsoon; and (2) a new post‐seismic hillslope failure occurred in saturated conditions and became fluidized and channelized. Runout distance was constrained by channel confluences and may be related to confluence geometry. Unstable landslide debris was largely flushed from steep channels during the first monsoon following the earthquake, and the rate of new hillslope failures tailed off over a few years.
-
more » « less
Abstract We analyzed continuous GPS data collected from 2002–2020 to characterize slow slip events (SSEs) in and near the Nicoya Peninsula, Costa Rica. These data are bisected by the 5 September 2012
M w 7.6 earthquake. The displacement time series contain multiple signals, including plate convergence, plate interface locking, coseismic and postseismic deformation, seasonal oscillations, SSEs, and noise. GPS‐measured coseismic and postseismic displacements associated with theM w 7.6 earthquake are modeled and removed by a step function plus multiple timescale relaxation processes with four characteristic times: 11, 94, 470, and 1,865 days. Seasonal oscillations are eliminated using a multichannel singular spectrum analysis (M‐SSA). Ten major SSEs (M w > 6.6) are observed in the remaining time series, with a constant recurrence interval of 21.7 ± 2.6 months. SSEs occur in both shallow (~10 km) and deep (~35 km) portions of the plate interface, but the latter last longer and have larger magnitudes. There is minimum to no slow slip in theM w 7.6 seismic rupture area and a persistent slow slip patch beneath the Nicoya Gulf entrance. Despite strong earthquake‐related stress perturbations, the inter‐SSE locking status on the megathrust is very similar between the late and early stages of the earthquake cycle and includes locked patches that ruptured in the 2012 earthquake or continue to rupture via SSEs. Some locked patches offshore south of the Nicoya Peninsula did not rupture in 2012, do not participate in SSEs, and may be indicative of supercycle behavior, that is, strain accumulation over several seismic cycles. These areas warrant heightened monitoring.
