Search

Rapid estimation of earthquake ground shaking and proper accounting of associated uncertainties in such estimates when conditioned on strong-motion station data or macroseismic intensity observations are crucial for downstream applications such as ground failure and loss estimation. The U.S. Geological Survey ShakeMap system is called upon to fulfill this objective in light of increased near-real-time access to strong-motion records from around the world. Although the station data provide a direct constraint on shaking estimates at specific locations, these data also heavily influence the uncertainty quantification at other locations. This investigation demonstrates methods to partition the within- (phi) and between-event (tau) uncertainty estimates under the observational constraints, especially when between-event uncertainties are heteroscedastic. The procedure allows the end users of ShakeMap to create separate between- and within-event realizations of ground-motion fields for downstream loss modeling applications in a manner that preserves the structure of the underlying random spatial processes.

The US Geological Survey’s ShakeMap is used domestically and globally for post-earthquake emergency management and response, engineering analyses, financial instruments, and other decision-making activities. Recent developments in the insurance, reinsurance, and catastrophe bond sectors link payouts of potentially hundreds of millions of dollars to ShakeMap products. Similarly, building codes, post-earthquake building damage forensic evaluations, and geotechnical evaluations often rely on estimated peak response-spectral values for site-specific evaluations that may lead to costly analyses, retrofits, or other expenditures. Given such activities, financial, engineering, and other technical users demand processing specifications and a metadata trail for actuarial, escrow, and forensic purposes for each significant earthquake. Recent inquiries include how and why maps change with time, how to interpret metadata, and how to obtain the creation and update history of various map layers. Similarly, the collection of ShakeMap scenarios and historical ShakeMaps—either created in earlier versions or rerun as part of the latest version of the ShakeMap Atlas—warrant a full explanation of the inputs, processing, and archiving given their contribution to fragility curve development and loss model calibration. For these reasons, in addition to event-specific ShakeMap metadata and a comprehensive online ShakeMap Manual, we have crafted this practice paper to answer several of the most frequently asked technical questions. We also describe an application programming interface (API) for accessing site-specific shaking metrics and their uncertainties for earthquake forensic purposes in a consistent fashion. In all, we describe the advantages of employing ShakeMaps for these critical purposes as well as describe their limitations and uncertainties, offering an extensive set of instructions and disclaimers that can be referenced by ShakeMap users.

Time-averaged shear wave velocity over the upper 30 m of the earth’s surface ( V S30) is a key parameter for estimating ground motion amplification as both a predictive and a diagnostic tool for earthquake hazards. The first-order approximation of V S30 is commonly obtained through a topographic slope–based or terrain proxy due to the widely available nature of digital elevation models. However, better-constrained V S30 maps have been developed in many regions. Such maps preferentially employ various combinations of V S30 measurements, higher-resolution elevation models, lithologic, geologic, geomorphic, and other proxies and often utilize refined interpolation schemes. We develop a new hybrid global V S30 map database that defaults to the global slope-based V S30 map, but smoothly inserts regional V S30 maps where available. In addition, we present comparisons of the default slope-based proxy maps against the new hybrid version in terms of V S30 and amplification ratio maps, and uncertainties in assigned V S30 values.

In the minutes to hours after a major earthquake, such as the recent 2018 Mw 7.1 Anchorage event, the U.S. Geological Survey (USGS) produces a suite of interconnected earthquake products that provides diverse information ranging from basic earthquake source parameters to loss estimates. The 2018 Anchorage earthquake is the first major domestic earthquake to occur since several new USGS products have been developed, thus providing an opportunity to discuss the newly expanded USGS earthquake product suite, its timeliness, performance, and reception. Overall, the products were relatively timely, accurate, well received, and widely used, including by the media, who used information and visualizations from many products to frame their early reporting. One downside of the codependence of multiple products is that reasonable updates to upstream products (e.g., magnitude and source characterization) can result in significant changes to downstream products; this was the case for the Anchorage earthquake. However, the coverage of strong‐motion stations and felt reports was so dense that the ShakeMap and downstream products were relatively insensitive to changes in magnitude or fault‐plane orientation once the ground‐motion data were available. Shaking and loss indicators initially fluctuated in the first hour or two after the earthquake, but they stabilized quickly. To understand how the products are being used and how effectively they are being communicated, we analyze the media coverage of USGS earthquake products. Most references to USGS products occurred within the first 48 hr after the event. The lack of coverage after 48 hr could indicate that longer‐term products addressing what actions the USGS is taking or what early reconnaissance has revealed might be useful for those people wanting additional information about the earthquake.

Modeling the spatial correlation of ground motion residuals, caused by coherent contributions from source, path, and site, can provide valuable loss and hazard information, as well as a more realistic depiction of ground motion intensities. The U.S. Geological Survey (USGS) software package, ShakeMap, utilizes a deterministic empirical approach to estimate median ground shaking in conjunction with observed seismic data. ShakeMap-based shaking estimates are used in concert with loss estimation algorithms to estimate fatalities and economic losses after significant seismic events around the globe. Incorporating the spatial correlation of ground motion residuals has been shown to improve seismic loss estimates. In particular, Park, Bazzuro, and Baker (Applications of Statistics and Probability in Civil Engineering, 2007) investigated computing spatially correlated random fields of residuals. However, for large scale ShakeMap grids, computational requirements of the method are prohibitive. In this work, a memory efficient algorithm is developed to compute the random fields and implemented using the ShakeMap framework. This new, iterative parallel algorithm is based on decay properties of an associated ground motion correlation function and is shown to significantly reduce computational requirements associated with adding spatial variability to the ShakeMap ground motion estimates. Further, we demonstrate and quantify the impact of adding peak ground motion spatial variability on resulting earthquake loss estimates. HighlightsSuccessive computations increase the efficiency of computing correlated random fields.The use of parallel processes significantly reduces CPU time.Method reduces memory constraints and computes realizations with near-linear speedup.Method is suitable for large ShakeMap grids and can condition upon station data.Realizations are used to compute valuable seismic loss distributions.

We analyzed the regional dependence of ground motion to intensity conversion equations and derived a new global relationship to improve ground motion and intensity estimates for earthquake hazard applications, including those related to the ShakeMap system. For this purpose, we merged several databases collected by other authors in different geographical regions to highlight any systematic regional effects in the relationship between macroseismic intensities and both peak ground velocity and peak ground acceleration. Our database contains macroseismic intensities derived from expert assignments or from the “Did You Feel It?” database, paired with peak ground motions (PGM) from seismic stations. We constrain our intensity–ground‐motion pairs to those with a maximum 2 km separation. For each region, we derived invertible relationships between intensities and ground motion using an orthogonal regression. We also derived a global relationship to quantify the regional differences. We investigated the dependence of intensity on predictor variables such as PGM, magnitude, and hypocentral distance. Our analyses indicate that PGM is the most robust predictor variable of intensity. Within one standard deviation, our regional and global results are in agreement with the relations of Worden et al. (2012) for California, Faenza and Michelini (2010) for Italy, Tselentis and Danciu (2008) for Greece, and Atkinson and Kaka (2007) for central–eastern United States. The earthquakes in the study ranged in magnitude from 2.5 to 7.3, and the distances ranged from less than a kilometer to about 200 km from the epicenter. Online Material: Table summarizing published relationships and figures showing the relationship between peak ground acceleration and intensity.

Equations that predict intensity as a function of magnitude and distance are useful tools for hazard and risk assessment, and in interpretation of both contemporary and historical earthquake information. The intensity prediction equations of Atkinson and Wald (2007; hereafter AW07) have been remarkably successful in describing the level and intensity of motions reported under the “Did You Feel It?” (DYFI) program over the last several years. Examination of the performance of AW07 for North American earthquakes, evaluated using an extensive compiled database of DYFI observations from 2000 to 2013, suggests that there is little statistical basis for revising these equations. However, a problem with the AW07 equations is that they predict unrealistically large median intensities for large events ( M >6) at close distances. In this study, we revise AW07 to improve the intensity scaling at large magnitudes and close distances, by reconciling intensity equations with ground‐motion prediction equations.

We develop globally applicable macroseismic intensity prediction equations (IPEs) for earthquakes of moment magnitude M W 5.0–7.9 and intensities of degree II and greater for distances less than 300 km for active crustal regions. The IPEs are developed for two distance metrics: closest distance to rupture (R rup) and hypocentral distance (R hyp). The key objective for developing the model based on hypocentral distance—in addition to more rigorous and standard measure R rup—is to provide an IPE which can be used in near real-time earthquake response systems for earthquakes anywhere in the world, where information regarding the rupture dimensions of a fault may not be known in the immediate aftermath of the event. We observe that our models, particularly the model for the R rup distance metric, generally have low median residuals with magnitude and distance. In particular, we address whether the direct use of IPEs leads to a reduction in overall uncertainties when compared with methods which use a combination of ground-motion prediction equations and ground motion to intensity conversion equations. Finally, using topographic gradient as a proxy and median model predictions, we derive intensity-based site amplification factors. These factors lead to a small reduction of residuals at shallow gradients at strong shaking levels. However, the overall effect on total median residuals is relatively small. This is in part due to the observation that the median site condition for intensity observations used to develop these IPEs is approximately near the National Earthquake Hazard Reduction Program CD site-class boundary.

Rapid (3-5 minutes) generation of maps of instrumental ground-motion and shaking intensity is accomplished through advances in real-time seismographic data acquisition combined with newly developed relationships between recorded ground-motion parameters and expected shaking intensity values. Estimation of shaking over the entire regional extent of southern California is obtained by the spatial interpolation of the measured ground motions with geologically based frequency and amplitude-dependent site corrections. Production of the maps is automatic, triggered by any significant earthquake in southern California. Maps are now made available within several minutes of the earthquake for public and scientific consumption via the World Wide Web; they will be made available with dedicated communications for emergency response agencies and critical users.