Search USGS⌕ Search
Home/Archive/USGS

USGS · 70195695

Spatial and spectral interpolation of ground-motion intensity measure observations

Charles Worden·Eric M. Thompson·Jack W. Baker·Brendon A. Bradley·Nico Luco·David J. Wald
View original source ↗Export metadata ↓

Abstract

Following a significant earthquake, ground‐motion observations are available for a limited set of locations and intensity measures (IMs). Typically, however, it is desirable to know the ground motions for additional IMs and at locations where observations are unavailable. Various interpolation methods are available, but because IMs or their logarithms are normally distributed, spatially correlated, and correlated with each other at a given location, it is possible to apply the conditional multivariate normal (MVN) distribution to the problem of estimating unobserved IMs. In this article, we review the MVN and its application to general estimation problems, and then apply the MVN to the specific problem of ground‐motion IM interpolation. In particular, we present (1) a formulation of the MVN for the simultaneous interpolation of IMs across space and IM type (most commonly, spectral response at different oscillator periods) and (2) the inclusion of uncertain observation data in the MVN formulation. These techniques, in combination with modern empirical ground‐motion models and correlation functions, provide a flexible framework for estimating a variety of IMs at arbitrary locations.

Explore related subjects

Bulletin of the Seismological Society of AmericaArticle

Keep this discovery

Explore connections, maps & timelines

BibTeXRIS

Charles Worden, Eric M. Thompson, Jack W. Baker, Brendon A. Bradley, Nico Luco, David J. Wald. 2018-02-13. Spatial and spectral interpolation of ground-motion intensity measure observations. https://doi.org/10.1785/0120170201

Cite the original work for its findings. Save a collection to share your selection of sources.

KEEP EXPLORING

Related USGS reports

USGS2026

Re-evaluation of shaking, dynamic strain, and the source of the 1886 Charleston, South Carolina, earthquake

As the largest historical event in the eastern United States, the 1886 Charleston, South Carolina, earthquake, plays a prominent role in our nation’s seismological and human history and an outsized role in our probabilistic ground‐motion hazard estimates, affecting building codes, insurance rates, and geophysical monitoring. Thus, constraining the characteristics of the 1886 earthquake source, as well as the nature of its shaking distribution and impacts to the extent possible, is essential. Given its vintage, however, societal, geological, and geophysical datasets are limited compared with those used in modern seismological analyses. Because of these limitations, some characteristics of the rupture model are fundamentally nonunique. Still, by re‐evaluating all available data, I can constrain several key aspects of this earthquake. I show that several interpretations in recent analyses are inconsistent regarding the inferred source; shaking; ground failure; and, particularly, railway impacts. Importantly, I have shown that although dynamic strain caused the ground failure that led to track deformation, without the intermediate liquefaction step, dynamic strain alone was unlikely to be the source of rail deformation, contrary to earlier studies. Had earlier pronouncements that dynamic strain is a real concern for bending railways held true, there would be serious ramifications for railway engineering as well as for other above‐ and below‐ground linear infrastructure.

South Carolina↗
USGS2026

ETAS-positive: Extracting unbiased epidemic-type aftershock models from incomplete earthquake catalogs

Earthquake catalogs are inherently incomplete, and short-term fluctuations in catalog completeness can strongly bias estimates of seismicity parameters. Here, we extend the positive statistics framework recently introduced for estimating Gutenberg–Richter a and b values to the Epidemic-Type Aftershock Sequence (ETAS) model. The approach focuses on complete magnitude–time intervals defined by pairs of earthquakes in which the second event is larger than the first. Comprehensive synthetic tests demonstrate that the proposed ETAS + method retrieves nearly unbiased parameter estimates even from highly incomplete catalogs, whereas parameters obtained using classical estimation approaches are strongly biased. Application to Southern California reveals a similar trend, with the productivity scaling parameter α approaching the b value, consistent with self-similar earthquake triggering within estimation uncertainties. Owing to its robustness against catalog incompleteness and computational efficiency, the ETAS + method provides a powerful new tool for seismicity analysis and earthquake forecasting.

California↗