skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Constraining Cascadia Subduction Zone Ground Motions via Paleoliquefaction Evidence: A Case Study from Kellogg Island, Washington, with Regional Implications
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
Award ID(s):
2103713 1751216
PAR ID:
10521782
Author(s) / Creator(s):
; ;
Publisher / Repository:
American Society of Civil Engineers
Date Published:
ISBN:
9780784485316
Page Range / eLocation ID:
141 to 151
Format(s):
Medium: X
Location:
Vancouver, British Columbia, Canada
Sponsoring Org:
National Science Foundation
More Like this
  1. Given the likelihood of future M9 Cascadia Subduction Zone (CSZ) earthquakes, various estimates of the resulting, regional ground motions have been made, including a suite of 30 physics-based simulations that reflect key modeling uncertainties. However, because the last CSZ interface rupture occurred in 1700 CE, the shaking expected in such an event is especially uncertain, as are the impacts to the built and living environments. Like other coseismic impacts, soil liquefaction poses a significant threat and must be considered by any scenario study used to inform planning and response, or to focus mitigation resources. Liquefaction is also notable for its potential to “ground truth” ground-motion estimates, given that its presence or absence in the geologic record can provide constraint on the intensities of shaking in past events. It is thus an important phenomenon looking both forward and backward. Accordingly, using recent physics-based simulations, this study (1) predicts liquefaction in M9 CSZ ruptures at 400 locations in Oregon, Washington, and British Columbia (BC) using an array of cone-penetration-test based models and (2) uses paleoliquefaction evidence at ten sites spanning from Southern Oregon to Vancouver, BC to constrain possible ground-motion intensities experienced in the 1700 CE earthquake. The forward predictions indicate that liquefaction in M9 events could be pervasive in the region and affect numerous population hubs, with the potential for damage across hundreds of square kilometers. The backward analyses suggest that 1700 CE ground-motion intensities may have been less than expected from M9 simulations in some northern portions of the CSZ (e.g. Seattle), given the paucity of 1700 CE liquefaction evidence in these areas. Ultimately, further discovery and analysis of CSZ paleoliquefaction, or lack thereof, will confirm or modify this possibility and the conclusions drawn herein. 
    more » « less
  2. In the U.S. Pacific Northwest (PNW), the historic earthquake record is often insufficient to provide inputs to seismic-hazard analyses or to inform ground-motion predictions for certain seismic sources (e.g., the Cascadia Subduction Zone, CSZ). As a result, paleoseismic studies are commonly used to infer information about the seismic hazard. However, among the many forms of coseismic evidence, soil liquefaction provides the best, if not only, evidence from which the intensities of previous ground motions may be constrained. Accordingly, the overarching goal of this research is to use paleoliquefaction to elucidate previous ground motions in the PNW – both for CSZ events and others – and to further constrain the locations, magnitudes, and recurrence rates of such ruptures. Towards that goal, this paper: (i) reviews current paleoliquefaction inverse-analysis methods and their limited, prior applications in the PNW; (ii) compiles all PNW paleoliquefaction evidence from the literature into a GIS database, resulting in data from 185 study sites (e.g., feature locations, types, sizes, and ages); and (iii) develops maps – specific to the CSZ – that forecast paleoliquefaction for 30 different simulations of a CSZ event. These maps can be used to guide field explorations for new evidence, such that they are conducted efficiently and strategically, considering the apparent utility of evidence toward constraint of CSZ ground-motion models. Of additional utility, this process provides regional ground-motion predictions for physics-based simulations of an M9 event, to include expected site effects. Collectively, the maps of expected shaking intensity and liquefaction may be useful in downstream hazard modelling, regional loss estimation, policy development, and science communication. Ultimately, as more paleoliquefaction evidence is identified and studied, better constraint of regional ground-motion hazards will result. Version 2 (this posting) supersedes Version 1 (10.17603/ds2-jm19-2w09). Updates include GIS rasters that provide regional ground-motion intensity predictions (PGA, PGV) for 30 physics-based simulations of an M9 event, to include expected site effects 
    more » « less
  3. In regions of infrequent moderate-to-large earthquakes, historic earthquake catalogs are often insufficient to provide inputs to seismic-hazard analyses (i.e. fault locations and magnitude–frequency relations) or to inform ground-motion predictions for certain seismic sources. In these regions, analysis of relic coseismic evidence, such as paleoliquefaction, is commonly used to infer information about the seismic hazard. However, while paleoliquefaction studies have been performed widely, all existing analysis techniques require a priori assumptions about the causative earthquake’s location (i.e. rupture magnitude and ground motions can otherwise not be estimated). This may lead to inaccurate assumptions in some settings, and by corollary, erroneous results. Accordingly, this article proposes an inversion framework to probabilistically constrain seismic-source parameters from paleoliquefaction. Analyzing evidence at regional scale leads to (a) a geospatial likelihood surface that constrains the rupture location and (b) a probability distribution of the rupture magnitude, wherein source-location uncertainty is explicitly considered. Simulated paleoliquefaction studies are performed on earthquakes with known parameters. These examples demonstrate the framework’s potential, even in cases of limited field evidence, as well as important caveats and lessons for forward use. The proposed framework has the potential to provide new insights in enigmatic seismic zones worldwide. 
    more » « less
  4. Abstract From California to British Columbia, the Pacific Northwest coast bears an omnipresent earthquake and tsunami hazard from the Cascadia subduction zone. Multiple lines of evidence suggests that magnitude eight and greater megathrust earthquakes have occurred ‐ the most recent being 321 years ago (i.e., 1700 A.D.). Outstanding questions for the next great megathrust event include where it will initiate, what conditions are favorable for rupture to span the convergent margin, and how much slip may be expected. We develop the first 3‐D fully dynamic rupture simulations for the Cascadia subduction zone that are driven by fault stress, strength and friction to address these questions. The initial dynamic stress drop distribution in our simulations is constrained by geodetic coupling models, with segment locations taken from geologic analyses. We document the sensitivity of nucleation location and stress drop to the final seismic moment and coseismic subsidence amplitudes. We find that the final earthquake size strongly depends on the amount of slip deficit in the central Cascadia region, which is inferred to be creeping interseismically, for a given initiation location in southern or northern Cascadia. Several simulations are also presented here that can closely approximate recorded coastal subsidence from the 1700 A.D. event without invoking localized high‐stress asperities along the down‐dip locked region of the megathrust. These results can be used to inform earthquake and tsunami hazards for not only Cascadia, but other subduction zones that have limited seismic observations but a wealth of geodetic inference. 
    more » « less
  5. Landslides are a significant hazard and dominant feature throughout the landscape of the Pacific Northwest. However, the hazard and risk posed by coseismic landslides triggered by great Cascadia Subduction Zone (CSZ) earthquakes is highly uncertain due to a lack of local and global data. Despite a wealth of other geologic evidence for past earthquakes on the Cascadia Subduction Zone, no landslides have been definitively linked to such earthquakes, even in areas otherwise highly susceptible to failure. While shallow landslides may not leave a lasting topographical signature in the landscape, there are thousands of deep-seated landslides in Cascadia, and these deposits often persist for hundreds of years and multiple earthquake cycles. Synthesizing newly developed inventories of dated large deep-seated landslides in the Oregon Coast Range, we use statistical methods to estimate the proportion of these types of landslides that could have been triggered during past great Cascadia Subduction Zone earthquakes. Statistical analysis of high-precision dendrochronology ages of landslide-dammed lakes and surface roughness-dated bedrock landslides reveal Cascadia Subduction Zone earthquakes may have triggered 0–15 % of large deep-seated landslides in the Oregon Coast Range over multiple earthquake cycles. Our results refine estimates from previous studies and further suggest that coseismic triggering accounts for a small fraction of the total deep-seated bedrock landslides mapped in coastal Cascadia. However, if the real rate of coseismic landslide triggering during CSZ earthquakes is near our estimated upper bound for the 1700 CSZ earthquake, we estimate up to 2400 coseismic large deep-seated landslides could occur in the Oregon Coast Range in a single earthquake. These findings suggest Cascadia is consistent with global observations from other subduction zones and that coseismic landslides may still represent a serious geohazard in the region. 
    more » « less