Geology topics

Valerie J. Sahakian

Publications and source records attributed to Valerie J. Sahakian.

13 recordsLinked to original sources

USGS2026

The Cascadia Region Earthquake Science Center (CRESCENT): Advancing understanding of Cascadia's earthquake hazards

The Cascadia Subduction Zone (CSZ) is capable of producing large megathrust earthquakes and tsunamis, posing significant risks to the Pacific Northwest. This paper examines the seismic characteristics of the CSZ, including its offshore coupling, sparse seismicity, and a complex paleoseismic record. While recent advances in geophysical methods and data collection have improved our understanding of subduction zone processes, important uncertainties remain regarding rupture dynamics, ground shaking, and tsunami impacts. The Cascadia Region Earthquake Science Center (CRESCENT) has been established to address these challenges through coordinated research, education, and community engagement. By integrating scientific studies, fostering collaboration, and leveraging new technologies, CRESCENT seeks to improve hazard assessment and resilience in the region. This paper outlines the primary challenges and opportunities in advancing subduction zone science and enhancing societal preparedness for future CSZ events.

California, Oregon, Washington↗

USGS2026

Cascadia Subduction Zone science: Call for the next generation community seismic velocity model

The Cascadia subduction zone (CSZ) hosts major seismic and tsunami hazards, yet key questions persist about the relationship between margin structure, fluid distribution, episodic tremor and slip, shallow megathrust behavior, shaking and tsunamigenesis, and the resulting hazard estimates. Addressing these problems requires an empirically grounded, three‐dimensional seismic velocity model to illuminate subsurface structure and properties and to provide a basis for geophysical studies such as earthquake simulations and ground‐motion estimation. In May 2024, the National Science Foundation‐funded Cascadia Region Earthquake Science Center (CRESCENT) community velocity model (CVM) working group, with U.S. Geological Survey and regional partners, convened a workshop to identify priorities for such a model. Participants emphasized the features necessary for addressing key science questions, including implementing findability, accessibility, interoperability, and reusability (FAIR) access, capturing along‐strike and along‐dip structural heterogeneity, resolving shallow offshore–onshore structure, constraining elastic properties and quantifying their uncertainties for numerical wave propagation simulations, their validation benchmarks, and supporting associated accurate earthquake ground‐motion simulations and hazard assessments. This article describes the priorities defined in the workshop, and a description of how, guided by these needs, CRESCENT plans to develop multiple generations of a CVM to advance CSZ science and improve seismic and tsunami hazard modeling across the Pacific Northwest. The CVM will span the CSZ from the surface to ∼100 km depth, offshore and east of the Cascades into Idaho (∼132°–110° W) and the southern and northern tectonic regime transitions (∼36°–52° N) to capture the entire tectonic system as well as its surroundings.

Cascadia Subduction Zone↗

USGS2025

Slow rupture, long rise times, and multi-fault geometry: The 2020 M6.4 southwestern Puerto Rico mainshock

The M 6.4 mainshock of the southwestern Puerto Rico seismic sequence on 7 January 2020, was one of the most impactful modern earthquakes in the northeastern Caribbean. Due to its offshore location and complex aftershock distribution, its source kinematics remain poorly constrained. This active sequence illuminated a complex set of previously unrecognized structures that indicate multiple causative faults may have slipped during its rupture. Here, we utilize seismic and geodetic observations to enhance model resolution, estimate the finite slip of the mainshock, and test a multi-segment, geologically realistic fault geometry. Our refined model finds a lower rupture velocity and longer rise times than typical for an event of this magnitude. This indicates a slow-evolving rupture process that resembles characteristics of a tsunami earthquake. Although this normal/strike-slip faulting event was not tsunamigenic, these qualities, if pervasive for this region, have important implications for future seismic monitoring and hazards in southwestern Puerto Rico.

Puerto Rico↗

USGS2025

Submarine avalanche deposits hold clues to past earthquakes

Earthquakes and other natural events sometimes shake the seafloor near coastlines severely enough to cause underwater avalanches that rush down steep slopes, scouring the seabed and carrying sediment to greater depths. These fast-moving sediment-laden flows, called turbidity currents , have at times damaged underwater infrastructure like pipelines and communications cables, as they did, for example, in snapping transatlantic cables off the coast of Newfoundland after the 1929 Grand Banks earthquake.

EOS Transactions↗

USGS2022

Estimates of k0 and effects on ground motions in the San Francisco Bay area

Ground‐motion studies are a key component of seismic hazard analyses and often rely on information of the source, path, and site. Extensive research has been done on each of these parameters; however, site‐specific studies are of particular interest to seismic hazard studies, especially in the field of earthquake engineering, as near‐site conditions can have a significant impact on the resulting ground motion at a site. There has been much focus on the constraint of site parameters and their application to seismic hazard studies, especially in the development of ground‐motion models (GMMs). Kappa is an observational parameter describing the high‐frequency attenuation of spectra, and its site contribution ( ⁠ κ 0 "> κ 0 ) has shown to be a good predictor of high‐frequency ground motions; however, measurements are often limited. In this study, we develop a κ 0 "> κ 0 dataset for the San Francisco Bay area (SFBA) by estimating κ 0 "> κ 0 for 228 stations, and we produce a continuous regional map of κ 0 "> κ 0 ⁠ . We find κ 0 "> κ 0 to range between 0.003 and 0.072 s, with larger values concentrating on the east, north, and south sides of the bay, and lower values concentrating on the west side. We also evaluate the robustness of κ 0 "> κ 0 as a site parameter and find it to correlate with peak ground acceleration. These estimates of κ 0 "> κ 0 can add predictive power to GMMs, thus increasing the accuracy of predicted ground motion and improving the robustness of ground‐motion studies in the SFBA.

California↗

USGS2022

The occurrence and hazards of great subduction zone earthquakes

Subduction zone earthquakes result in some of the most devastating natural hazards on Earth. Knowledge of where great (moment magnitude M ≥ 8) subduction zone earthquakes can occur and how they rupture is critical to constraining future seismic and tsunami hazards. Since the occurrence of well-instrumented great earthquakes, such as the 2004 M 9.1 Sumatra–Andaman and 2011 M 9.1 Tohoku earthquakes, the hypotheses that plate age and convergence rate influence the ability of subduction zones to host large earthquakes have been dispelled. In this Review, we highlight how certain subduction zone properties might influence the location and characteristics of great earthquake rupture and impact seismic and tsunami hazard. The rupture characteristics of great earthquakes that most heavily impact earthquake hazards include the rupture extent (seaward and landward), location of strong motion-generating areas and earthquake recurrence. By contrast, large slip or displacement at the seafloor is one of the major controls of tsunami hazard. Future improvements in addressing hazards posed by subduction zones depend heavily on sustained geophysical monitoring in subduction zone systems (both onshore and offshore), expanded development of palaeoseismic data sets and improved integration of observations and models across disciplines and timescales.

Nature Reviews Earth & Environment↗

USGS2022

Zoning Verification in Mexico City Using Strong Motions of the M7.1 Puebla-Morelos Earthquake of September 19, 2017

Mexico City suffers extensive damage from large earthquakes that originate at far distances due mainly to densely built areas on a filled lakebed. Seismic design codes in Mexico recognize the site-specific issues in Mexico City by acknowledging zones that represent lakebed as being riskier compared to other Mexico City areas. At the time of the 19 September 1985 M8.1 Michoacán Earthquake, (a) the seismic zoning map comprised only three zones (Hills [now Zone I], Transition [now Zone II] and Lake Zone [now Zone III] and (b) a limited number of permanent strong motion stations (~6) existed. Since, 1985 the seismic zonation maps have evolved: (a) the Lake Zone (now Zone III) into 4 sub-zones (a, b, c and d) and (b) the number of strong motion stations has increased significantly (~25) and recorded the 19 September 2017 M7.1 Puebla-Morelos earthquake. To identify the site frequency [f] (period [T]) at the specific location of a strong motion station within a particular zone or sub-zone (Zone I, II and III a,b,c,d), two well-known methods (H/H and H/V) for determining spectral ratios are employed using these records. The predominant site periods identified by both methods are in good agreement with those interpolated from the zoning-site period maps of Mexico City in the 2004 Seismic Design Code. Such verification is a testament to the benefits of improved zoning and greater distribution of strong motion network stations, and sets a good example for applications in other seismic regions.

Mexico City↗

USGS2022

Estimates of kappa in the San Francisco Bay area

Site characterization is a critical component of seismic hazards studies, especially in the development and use of ground motion models (GMMs). One such parameter, kappa ( Κ 0 ), represents local site attenuation and effectively describes regional variations in ground motion [1]. However, estimates of Κ 0 are limited. We estimate the site parameter Κ 0 for 296 broadband and accelerometer stations in the San Francisco Bay area using 225 local events with M>3.5. We first invert for kappa on the high-frequency slope of every record at a station, which represents full attenuation effects [2]. We then perform weighted regression of kappa versus distance to obtain expected kappa at the site ( Κ 0 ). We find that Κ 0 varies extensively in the San Francisco Bay area, with values ranging from 0.004-0.113 s, thus representing the heterogeneous geologic conditions in the region. Κ 0 estimates from this study will be useful to consider in future seismic hazard analyses and ground motion studies in the Bay Area.

California↗

USGS2021

A ground motion model for GNSS peak ground displacement

We present an updated ground‐motion model (GMM) for M w "> M w Mw 6–9 earthquakes using Global Navigation Satellite Systems (GNSS) observations of the peak ground displacement (PGD). Earthquake GMMs inform a range of Earth science and engineering applications, including source characterization, seismic hazard evaluations, loss estimates, and seismic design standards. A typical GMM is characterized by simplified metrics describing the earthquake source (magnitude), observation distance, and site terms. Most often, GMMs are derived from broadband seismometer and accelerometer observations, yet during strong shaking, these traditional seismic instruments are affected by baseline offsets, leading to inaccurate recordings of low‐frequency ground motions such as displacement. The incorporation of geodetic data sources, particularly for characterizing the unsaturated ground displacement of large‐magnitude events, has proven valuable as a complement to traditional seismic approaches and led to the development of an initial point‐source GMM based on PGD estimated from high‐rate GNSS data. Here, we improve the existing GMM to more effectively account for fault finiteness, slip heterogeneity, and observation distance. We evaluate the limitations of the currently available GNSS earthquake data set to calibrate the GMM. In particular, the observed earthquake data set is lacking in observations within 100 km of large‐magnitude events ( ⁠ M w > 8 "> M w > 8 Mw>8 ⁠ ), inhibiting evaluation of fault dimensions for earthquakes too large to be represented as point sources in the near field. To that end, we separately consider previously validated synthetic GNSS waveforms within 10–1000 km of M w "> M w Mw 7.8–9.3 Cascadia subduction zone scenario ruptures. The synthetic data highlight the importance of fault distance rather than point‐source metrics and improve our preparedness for large‐magnitude earthquakes with spatiotemporal qualities unlike those in our existing data set.

Bulletin of the Seismological Society of America↗

USGS2021

Toward an integrative geological and geophysical view of Cascadia subduction zone earthquakes

The Cascadia subduction zone (CSZ) is an exceptional geologic environment for recording evidence of land level changes, tsunamis, and ground motion that reveals at least 19 great megathrust earthquakes over the past 10 kyr. Such earthquakes are among the most impactful natural hazards on Earth, transcend national boundaries, and can have global impact. Reducing the societal impacts of future events in the U.S. Pacific Northwest and coastal British Columbia, Canada requires improved scientific understanding of megathrust earthquake rupture, recurrence, and corresponding hazards. Despite substantial knowledge gained from decades of research, large uncertainties remain about the characteristics and frequencies of past CSZ earthquakes. In this review, we summarize geological, geophysical, and instrumental evidence relevant to understanding megathrust earthquakes along the CSZ and associated uncertainties. We discuss how the evidence constrains various models of great megathrust earthquake recurrence in Cascadia and identify potential paths forward for the earthquake science community.

British Columbia, Oregon, Washington↗

USGS2019

κ0 and broadband site spectra in Southern California from source model-constrained inversion

Ground-motion modeling requires accurate representation of the earthquake source, path, and site. Site amplification is often modeled by VS30, the time-averaged shear-wave velocity of the top 30 meters of the Earth’s surface, though recent studies find that its ability to accurately predict site effects varies. Another measure of the site is κ0, the attenuation of high frequency energy near the site (Anderson & Hough, 1984). We develop a novel application of the Andrews (1986) method to simultaneously invert the spectra of 3,357 earthquakes in Southern California into source and site components. These earthquakes have magnitudes 2.5 to 5.72 and were recorded on 16 stations for a total of 52,297 records. We constrain the inversion with an individual earthquake demonstrating the most Brune-like shape to preserve the site spectra. We then solve for κ0 site amplification at each station in three frequency bands: 1-6 Hz, 6-14 Hz, and 14-35 Hz. The resulting values of κ0 range from 0.017 seconds at ANZA station PFO to 0.059 seconds at ANZA station SND. We compare our results with values of site κ0 from other studies as well as site residuals from GMPEs. We find good agreement between our site κ0 and previous studies in the region. We find that κ0 and high frequency site amplification (14-35 Hz band) correlates well with independent site residuals, making it a good first-order approximation for the effects of site attenuation or amplification on ground motion.

California↗

USGS2019

Ground-motion residuals, path effects, and crustal properties: A pilot study in southern California

To improve models of ground motion estimation and probabilistic seismic hazard analyses, the engineering seismology field is moving toward developing fully nonergodic ground motion models, models specific for individual source‐to‐site paths. Previous work on this topic has examined systematic variations in ground‐motion along particular paths (from either recorded or simulated earthquake data) and has not included physical properties of the path. We present here a framework to include physical path properties, by seeking correlations between ground motion amplitudes along specific paths and crustal properties, specifically seismic velocity and anelastic attenuation, along that path. Using a large data set of small‐magnitude earthquakes recorded in Southern California, we find a correlation between the gradient of seismic S wave velocity and the path term residual, after accounting for an average geometric spreading and anelastic attenuation, indicating that heterogeneity in crustal velocity primarily controls the path‐specific attenuation. Even in aseismic regions, details of path‐specific ground motion prediction equations can be developed from crustal structure and property data.

California↗

USGS2018

Ground motions from the 7 and 19 September, 2017 Tehuantepec and Puebla-Morelos, Mexico, earthquakes

The 2017 M "> M 8.2 Tehuantepec and M "> M 7.1 Puebla‐Morelos earthquakes were deep inslab normal‐faulting events that caused significant damage to several central‐to‐southern regions of Mexico. Inslab earthquakes are an important component of seismicity and seismic hazard in Mexico. Ground‐motion prediction equations (GMPEs) are an integral part of seismic hazard assessment as well as risk and rapid‐response products. This work examines the observed ground motions from these two events in comparison to the predicted median ground motions from four GMPEs. The residuals between the observed and modeled ground motions allow us to study regional differences in shaking, the effects of each earthquake, and basin effects in Mexico City, Puebla, and Oaxaca. We find that the ground motions from these two earthquakes are generally well modeled by the GMPEs. However, the Tehuantepec event shows larger than expected ground motions at greater distances and longer periods, which suggests a waveguide effect from the subduction zone geometry. Finally, Mexico City and the cities of Puebla and Oaxaca exhibit very large ground motions, indicative of well‐known site and basin effects that are much stronger than the basin terms included in some of the GMPEs. Simple and rapid ground‐motion parameter estimates that include site effects are key for hazard and real‐time risk assessments in regions such as Mexico, where the vast majority of the population lives in areas where the aforementioned effects are relevant. However, GMPEs based on site correction terms dependent on topographic slope proxies underestimate, at least in the three cities tackled in this work, the observed amplification. Therefore, there is a need to improve models of seismic amplification in basins that could be included in GMPEs.

Bulletin of the Seismological Society of America↗