skip to main content

Explore Research Products in the NSF-PAR

  • NSF Public Access
  • Search Results
  • Compilation and Forecasting of Paleoliquefaction Evidence for the Strength of Ground Motions in the U.S. Pacific Northwest: A Digital Dataset (Version 2)
Title: Compilation and Forecasting of Paleoliquefaction Evidence for the Strength of Ground Motions in the U.S. Pacific Northwest: A Digital Dataset (Version 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
Award ID(s):
1751216
NSF-PAR ID:
10344902
Author(s) / Creator(s):
; ;
Publisher / Repository:
Designsafe-CI
Date Published:
Subject(s) / Keyword(s):
["Paleoliquefaction","Paleoseismology","Seismic Hazard","Pacific Northwest","Washington","Oregon"]
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ( , The University of Oklahoma Libraries)
    null (Ed.)
    Nonstructural components within mission-critical facilities such as hospitals and telecommunication facilities are vital to a community's resilience when subjected to a seismic event. Building contents like medical and computer equipment are critical for the response and recovery process following an earthquake. A solution to protecting these systems from seismic hazards is base isolation. Base isolation systems are designed to decouple an entire building structure from destructive ground motions. For other buildings not fitted with base isolation, a practical and economical solution to protect vital building contents from earthquake-induced floor motion is to isolate individual equipment using, for example, rolling-type isolation systems (RISs). RISs are a relatively new innovation for protecting equipment. These systems function as a pendulum-like mechanism to convert horizontal motion into vertical motion. An accompanying change in potential energy creates a restoring force related to the slope of the rolling surface. This study seeks to evaluate the seismic hazard mitigation performance of RISs, as well as propose and test a novel double RIS. A physics-based mathematical model was developed for a single RIS via Lagrange's equation adhering to the kinetic constraint of rolling without slipping. The mathematical model for the single RIS was used to predict the response and characteristics of these systems. A physical model was fabricated with additive manufacturing and tested against multiple earthquakes on a shake table. The system featured a single-degree-of-freedom (SDOF) structure to represent a piece of equipment. The results showed that the RIS effectively reduced accelerations felt by the SDOF compared to a fixed-base SDOF system. The single RIS experienced the most substantial accelerations from the Mendocino record, which contains low-frequency content in the range of the RIS's natural period (1-2 seconds). Earthquakes with these long-period components have the potential to cause impacts within the isolation bearing that would degrade its performance. To accommodate large displacements, a double RIS is proposed. The double RIS has twice the displacement capacity of a single RIS without increasing the size of the bearing components. The mathematical model for the single RIS was extended to the double RIS following a similar procedure. Two approaches were used to evaluate the double RIS's performance: stochastic and deterministic. The stochastic response of the double RIS under stationary white noise excitation was evaluated for relevant system parameters, namely mass ratio and tuning frequency. Both broadband and filtered (Kanai-Tajimi) white noise excitation were considered. The response variances of the double RIS were normalized by a baseline single RIS for a comparative study, from which design parameter maps were drawn. A deterministic analysis was conducted to further evaluate the double RIS in the case of nonstationary excitation. The telecommunication equipment qualification waveform, VERTEQ-II, was used for these numerical simulations. Peak transient responses were compared to the single RIS responses, and optimal design regions were determined. General design guidelines based on the stochastic and deterministic analyses are given. The results aim to provide a framework usable in the preliminary design stage of a double RIS to mitigate seismic responses. 
    more » « less
  2. ; ; ; ; ( , Earthquake Engineering & Structural Dynamics)
    Summary

    Studies of recorded ground motions and simulations have shown that deep sedimentary basins can greatly increase the intensity of earthquake ground motions within a period range of approximately 1–4 s, but the economic impacts of basin effects are uncertain. This paper estimates key economic indicators of seismic performance, expressed in terms of earthquake‐induced repair costs, using empirical and simulated seismic hazard characterizations that account for the effects of basins. The methodology used is general, but the estimates are made for a series of eight‐ to 24‐story residential reinforced concrete shear wall archetype buildings in Seattle, WA, whose design neglects basin effects. All buildings are designed to comply with code‐minimum requirements (i.e., reference archetypes), as well as a series of design enhancements, which include (a) increasing design forces, (b) decreasing drift limits, and (c) a combination of these strategies. As an additional reference point, a performance‐based design is also assessed. The performance of the archetype buildings is evaluated for the seismic hazard level in Seattle according to the 2018 National Seismic Hazard Model (2018 NSHM), which explicitly considers basin effects. Inclusion of basin effects results in an average threefold increase in annualized losses for all archetypes. Incorporating physics‐based ground motion simulations to represent the large‐magnitude Cascadia subduction interface earthquake contribution to the hazard results in a further increase of 22% relative to the 2018 NSHM. The most effective of the design strategies considered combines a 25% increase in strength with a reduction in drift limits to 1.5%.

     
    more » « less
  3. ; ( , Earthquake Spectra)

    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. ; ; ; ; ( , Journal of Geophysical Research: Solid Earth)
    Abstract

    The Húsavík‐Flatey Fault Zone (HFFZ) is the largest strike‐slip fault in Iceland and poses a high seismic risk to coastal communities. To investigate physics‐based constraints on earthquake hazards, we construct three fault system models of varying geometric complexity and model 79 3‐D multi‐fault dynamic rupture scenarios in the HFFZ. By assuming a simple regional prestress and varying hypocenter locations, we analyze the rupture dynamics, fault interactions, and the associated ground motions up to 2.5 Hz. All models account for regional seismotectonics, topo‐bathymetry, 3‐D subsurface velocity, viscoelastic attenuation, and off‐fault plasticity, and we explore the effect of fault roughness. The rupture scenarios obey earthquake scaling relations and predict magnitudes comparable to those of historical events. We show how fault system geometry and segmentation, hypocenter location, and prestress can affect the potential for rupture cascading, leading to varying slip distributions across different portions of the fault system. Our earthquake scenarios yield spatially heterogeneous near‐field ground motions modulated by geometric complexities, topography, and rupture directivity, particularly in the near‐field. The average ground motion attenuation characteristics of dynamic rupture scenarios of comparable magnitudes and mean stress drop are independent of variations in source complexity, magnitude‐consistent and in good agreement with the latest regional empirical ground motion models. However, physics‐based ground motion variability changes considerably with fault‐distance and increases for unilateral compared to bilateral ruptures. Systematic variations in physics‐based near‐fault ground motions provide important insights into the mechanics and potential earthquake hazard of large strike‐slip fault systems, such as the HFFZ.

     
    more » « less
  5. ; ; ( , Geophysical Journal International)
    SUMMARY

    We have simulated 0–5 Hz deterministic wave propagation for a suite of 17 models of the 2014 Mw 5.1 La Habra, CA, earthquake with the Southern California Earthquake Center Community Velocity Model Version S4.26-M01 using a finite-fault source. Strong motion data at 259 sites within a 148 km × 140 km area are used to validate our simulations. Our simulations quantify the effects of statistical distributions of small-scale crustal heterogeneities (SSHs), frequency-dependent attenuation Q(f), surface topography and near-surface low-velocity material (via a 1-D approximation) on the resulting ground motion synthetics. The shear wave quality factor QS(f) is parametrized as QS, 0 and QS, 0fγ for frequencies less than and higher than 1 Hz, respectively. We find the most favourable fit to data for models using ratios of QS, 0 to shear wave velocity VS of 0.075–1.0 and γ values less than 0.6, with the best-fitting amplitude drop-off for the higher frequencies obtained for γ values of 0.2–0.4. Models including topography and a realistic near-surface weathering layer tend to increase peak velocities at mountain peaks and ridges, with a corresponding decrease behind the peaks and ridges in the direction of wave propagation. We find a clear negative correlation between the effects on peak ground velocity amplification and duration lengthening, suggesting that topography redistributes seismic energy from the large-amplitude first arrivals to the adjacent coda waves. A weathering layer with realistic near-surface low velocities is found to enhance the amplification at mountain peaks and ridges, and may partly explain the underprediction of the effects of topography on ground motions found in models. Our models including topography tend to improve the fit to data, as compared to models with a flat free surface, while our distributions of SSHs with constraints from borehole data fail to significantly improve the fit. Accuracy of the velocity model, particularly the near-surface low velocities, as well as the source description, controls the resolution with which the anelastic attenuation can be determined. Our results demonstrate that it is feasible to use fully deterministic physics-based simulations to estimate ground motions for seismic hazard analysis up to 5 Hz. Here, the effects of, and trade-offs with, near-surface low-velocity material, topography, SSHs and Q(f) become increasingly important as frequencies increase towards 5 Hz, and should be included in the calculations. Future improvement in community velocity models, wider access to computational resources, more efficient numerical codes and guidance from this study are bound to further constrain the ground motion models, leading to more accurate seismic hazard analysis.

     
    more » « less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Save / Share this Record
  • Send to Email