An official website of the United States government
ABSTRACT We explore the response of ground motions to topography during large crustal fault earthquakes by simulating several magnitude 6.5–7.0 rupture scenarios on the Seattle fault, Washington State. Kinematic simulations are run using a 3D spectral element code and a detailed seismic velocity model for the Puget Sound region. This model includes realistic surface topography and a near-surface low-velocity layer; a mesh spacing of ∼30 m at the surface allows modeling of ground motions up to 3 Hz. We simulate 20 earthquake scenarios using different slip distributions and hypocenter locations on a planar fault surface. Results indicate that average ground motions in simulations with and without topography are similar. However, shaking amplification is common at topographic highs, and more than a quarter of all sites experience short-period (≤2 s) ground-motion amplification greater than 25%–35%, compared with models without topography. Comparisons of peak ground velocity at the top and bottom of topographic features demonstrate that amplification is sensitive to period, with the greatest amplifications typically manifesting near a topographic feature’s estimated resonance frequency and along azimuths perpendicular to its primary axis of elongation. However, interevent variability in topographic response can be significant, particularly at shorter periods (<1 s). We do not observe a clear relationship between source centroid-to-site azimuths and the strength of topographic amplification. Overall, our results suggest that although topographic resonance does influence the average ground motions, other processes (e.g., localized focusing and scattering) also play a significant role in determining topographic response. However, the amount of consistent, significant amplification due to topography suggests that topographic effects should likely be considered in some capacity during seismic hazard studies.
more »
« less
Evidence of Seattle Fault Earthquakes from Patterns in Deep-Seated Landslides
ABSTRACT Earthquake-induced landslides can record information about the seismic shaking that generated them. In this study, we present new mapping, Light Detection and Ranging-derived roughness dating, and analysis of over 1000 deep-seated landslides from the Puget Lowlands of Washington, U.S.A., to probe the landscape for past Seattle fault earthquake information. With this new landslide inventory, we observe spatial and temporal evidence of landsliding related to the last major earthquake on the Seattle fault ∼1100 yr before present. We find spatial clusters of landslides that correlate with ground motions from recent 3D kinematic models of Seattle fault earthquakes. We also find temporal patterns in the landslide inventory that suggest earthquake-driven increases in landsliding. We compare the spatial and temporal landslide data with scenario-based ground motion models and find stronger evidence of the last major Seattle fault earthquake from this combined analysis than from spatial or temporal patterns alone. We also compare the landslide inventory with ground motions from different Seattle fault earthquake scenarios to determine the ground motion distributions that are most consistent with the landslide record. We find that earthquake scenarios that best match the clustering of ∼1100-year-old landslides produce the strongest shaking within a band that stretches from west to east across central Seattle as well as along the bluffs bordering the broader Puget Sound. Finally, we identify other landslide clusters (at 4.6–4.2 ka, 4.0–3.8 ka, 2.8–2.6 ka, and 2.2–2.0 ka) in the inventory which let us infer potential ground motions that may correspond to older Seattle fault earthquakes. Our method, which combines hindcasting of the surface response to the last major Seattle fault earthquake, using a roughness-aged landslide inventory with forecasts of modeled ground shaking from 3D seismic scenarios, showcases a powerful new approach to gleaning paleoseismic information from landscapes.
more »
« less
- PAR ID:
- 10521749
- Publisher / Repository:
- Seismological Society of America
- Date Published:
- Journal Name:
- Bulletin of the Seismological Society of America
- Volume:
- 114
- Issue:
- 2
- ISSN:
- 0037-1106
- Page Range / eLocation ID:
- 1084 to 1102
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
ABSTRACT We simulate shaking in Tacoma, Washington, and surrounding areas from Mw 6.5 and 7.0 earthquakes on the Tacoma fault. Ground motions are directly modeled up to 2.5 Hz using kinematic, finite-fault sources; a 3D seismic velocity model considering regional geology; and a model mesh with 30 m sampling at the ground surface. In addition, we explore how adjustments to the seismic velocity model affect predicted shaking over a range of periods. These adjustments include the addition of a region-specific geotechnical gradient, surface topography, and a fault damage zone. We find that the simulated shaking tends to be near estimates from empirical ground-motion models (GMMs). However, long-period (T = 5.0 s) shaking within the Tacoma basin is typically underpredicted by the GMMs. The fit between simulated and GMM-derived short-period (T = 0.5 s) shaking is significantly improved with the addition of the geotechnical gradient. From comparing different Mw 6.5 earthquake scenarios, we also find that the response of the Tacoma basin is sensitive to the azimuth of incoming seismic waves. In adding surface topography to the simulation, we find that average ground motion is similar to that produced from the nontopography model. However, shaking is often amplified at topographic highs and deamplified at topographic lows, and the wavefield undergoes extensive scattering. Adding a fault damage zone has the effect of amplifying short-period shaking adjacent to the fault, while reducing far-field shaking. Intermediate-period shaking is amplified within the Tacoma basin, likely due to enhanced surface-wave generation attributable to the fault damage zone waveguide. When applied in the same model, the topography and fault damage zone adjustments often enhance or reduce the effects of one another, adding further complexity to the wavefield. These results emphasize the importance of improving near-surface velocity model resolution as waveform simulations progress toward higher frequencies.more » « less
-
Tens of thousands of landslides were generated over 10,000 km2 of North Canterbury and Marlborough as a consequence of the 14 November 2016, Mw7.8 Kaikōura Earthquake. The most intense landslide damage was concentrated in 3500 km2 around the areas of fault rupture. Given the sparsely populated area affected by landslides, only a few homes were impacted and there were no recorded deaths due to landslides. Landslides caused major disruption with all road and rail links with Kaikōura being severed. The landslides affecting State Highway 1 (the main road link in the South Island of New Zealand) and the South Island main trunk railway extended from Ward in Marlborough all the way to the south of Oaro in North Canterbury. The majority of landslides occurred in two geological and geotechnically distinct materials reflective of the dominant rock types in the affected area. In the Neogene sedimentary rocks (sandstones, limestones and siltstones) of the Hurunui District, North Canterbury and around Cape Campbell in Marlborough, first-time and reactivated rock-slides and rock-block slides were the dominant landslide type. These rocks also tend to have rock material strength values in the range of 5-20 MPa. In the Torlesse ‘basement’ rocks (greywacke sandstones and argillite) of the Kaikōura Ranges, first-time rock and debris avalanches were the dominant landslide type. These rocks tend to have material strength values in the range of 20-50 MPa. A feature of this earthquake is the large number (more than 200) of valley blocking landslides it generated. This was partly due to the steep and confined slopes in the area and the widely distributed strong ground shaking. The largest landslide dam has an approximate volume of 12(±2) M m3 and the debris from this travelled about 2.7 km2 downslope where it formed a dam blocking the Hapuku River. The long-term stability of cracked slopes and landslide dams from future strong earthquakes and large rainstorms are an ongoing concern to central and local government agencies responsible for rebuilding homes and infrastructure. A particular concern is the potential for debris floods to affect downstream assets and infrastructure should some of the landslide dams breach catastrophically. At least twenty-one faults ruptured to the ground surface or sea floor, with these surface ruptures extending from the Emu Plain in North Canterbury to offshore of Cape Campbell in Marlborough. The mapped landslide distribution reflects the complexity of the earthquake rupture. Landslides are distributed across a broad area of intense ground shaking reflective of the elongate area affected by fault rupture, and are not clustered around the earthquake epicentre. The largest landslides triggered by the earthquake are located either on or adjacent to faults that ruptured to the ground surface. Surface faults may provide a plane of weakness or hydrological discontinuity and adversely oriented surface faults may be indicative of the location of future large landslides. Their location appears to have a strong structural geological control. Initial results from our landslide investigations suggest predictive models relying only on ground-shaking estimates underestimate the number and size of the largest landslides that occurred.more » « less
-
Physics-based ground motion simulations are a valuable tool for studying seismic sources with missing historical records, such as Cascadia Subduction Zone (CSZ) interface earthquakes. The last such event occurred in 1700 CE and is believed to be an M8-M9 rupture. The United States Geological Survey recently developed 30 physics-based simulations of a CSZ rupture to predict ground motions across the Pacific Northwest. Consideration of key modeling uncertainties across these simulations leads to estimates of ground motion intensity that vary by ~100% in some areas (e.g., Seattle). Paleoliquefaction, or soil liquefaction from past earthquakes, provides the best geologic evidence for constraining or "ground truthing" the intensity of past shaking, yet while paleoliquefaction has been documented throughout Cascadia, limited analyses have been performed to exploit this evidence. This study focuses on Kellogg Island, 2 mi south of Seattle, where liquefaction has been documented from several earthquakes, but not from the 1700 CE event. Therefore, using the CSZ simulations and in situ cone penetration test data, this study predicts the probability of surficial liquefaction manifestation at Kellogg Island during an M9 CSZ event. As part of this effort, velocity profiles are developed from multichannel analysis of surface waves, and non-linear site response analyses are used to propagate simulated motions to the surface. Results show a high probability of liquefaction near Kellogg Island for most simulations, whereas to date no evidence of 1700 CE liquefaction has been discovered at Kellogg Island, nor at any other location in the Puget Sound. The discrepancy between predictions and observations might indicate that the 1700 CE ground motions were less intense in Seattle than most predictions of M9 earthquakes indicate. Toward the goal of elucidating the expected impacts of future CSZ earthquakes, similar analyses are ongoing at additional sites across the region.more » « less
-
Abstract 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 magnitudeMw, seismic energy source depth, and failure susceptibility, as well as the seismological factor responsible for triggering landslides), and is moderately sensitive to the effective elastic thickness of lithosphere,Te. For a specified magnitude, more erosive events (higher Ω) tend to occur at shallower depth, in thicker‐Telithosphere, and in steeper, more landslide‐prone landscapes. For given landscape and seismotectonic conditions, the volumes of both landslides and uplift to first order positively scale withMwand seismic momentMo. However, higherMwearthquakes generate lower landslide and uplift volumes per unitMo, suggesting lower efficiency in the use of seismic energy to drive topographic change. With our model, we calculate the long‐term average seismic volume balance for the eastern Tibetan region and find that the net topographic effect of earthquakes in this region tends to be constructive rather than erosive. Overall, destructive events are rare when considering processes over the full earthquake cycle, although they are more likely if only considering the coseismic volume budget (as was the case for the 2008 Wenchuan earthquake where landsliding substantially offset coseismic uplift). Irrespective of the net budget, our results suggest that the erosive power of earthquakes plays an important role in mountain belt evolution, including by influencing structures and spatial patterns of deformation, for example affecting the wavelength of topography.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1785/0120230079
