USGS Science⌕ Search

SEARCH · USGS Science

Results for “The Seismic Record”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 811 records · Page 45Linked to original sources
USGS1989

Some significant records from instrumented structures in California - USGS program

The main objective in seismic instrumentation of structures is to facilitate response studies that lead to improved understanding of the dynamic behavior and the potential for damage to structures under seismic loading. The purpose of this paper is: (1) to review the status of the programs for strong-motion instrumentation of structures in the United States and discuss various procedures and instrumentation schemes designed to best acquire response data from buildings and (2) to discuss preliminary results derived from recorded response data obtained from a well-instrumented structure during the recent Whittier Narrows earthquake of October 1, 1987 (M3 = 5.6).

Conference Paper↗
USGS2025

2024 Surprise Inlet landslides: Insights from a prototype landslide‐triggered tsunami monitoring system in Prince William Sound, Alaska

Alaska's coastal communities face growing landslide hazards owing to glacier retreat and extreme weather intensified by the warming climate, yet hazard monitoring remains challenging. As part of ongoing experimental monitoring in Prince William Sound, we detected three large landslides (0.5–2.3 M m 3 ) at Surprise Inlet on 20 September 2024, within the span of an hour. These events were identified in near real-time through seismic data and later confirmed using satellite imagery, tidal records, and infrasound. The landslides generated a modest tsunami, and a 4 cm wave was recorded by a tide gauge 18 km away, marking the first recorded landslide to reach water since monitoring began in this region in 2021. Here, we examine the detection and interpretation of these landslides using multiple data sources and modeling. We demonstrate the effectiveness of this regional seismic monitoring system and show how complementary instrumentation, where available, can enhance detection capabilities.

Alaska↗
USGS1996

Sublake geologic structure from high-resolution seismic-reflection data from four sinkhole lakes in the Lake Wales Ridge, central Florida

Seismic-reflection profiles from Lake Wales, Blue Lake, Lake Letta, and Lake Apthorp located along the Lake Wales Ridge in central Florida provide local detail within the regional hydrogeologic framework as described by litho- and hydrostratigraphic cross sections. Lakes located with the mantled karst region have long been considered to be sinkhole lakes, originating from subsidence activity. High-resolution seismic- reflection data confirm this origin for these four lakes. The geologic framework of the Lake Wales Ridge has proven to be a suitable geologic setting for continuous high-resolution seismic-reflection profiling in lakes; however, the nature of the lake-bottom sediments largely controls the quality of the seismic data. In lakes with significant organic-rich bottom deposits, interpretable record was limited to areas where organic deposits were minimal. In lakes with clean, sandy bottoms, the seismic-reflection methods were highly successful in obtaining data that can be correlated with sublake subsidence features. These techniques are useful in examining sublake geology and providing a better understanding of how confining units are affected by subsidence in a region where their continuity is of significant importance to local lake hydrology. Although local geologic control around each lake generally corresponds to the regional geologic framework, local deviations from regional geologic trends occur in sublake areas affected by subsidence activity. Each of the four lakes examined represents a unique set of geologic controls and provides some degree of structural evidence of subsidence activity. Sublake geologic structures identified include: (1) marginal lake sediments dipping into bathymetric lows, (2) lateral discontinuity of confining units including sags and breaches, (3) the disruption and reworking of overlying unconsolidated siliciclastic sediments as they subside into the underlying irregular limestone surface, and (4) sublake regions where confining units appear to remain intact and unaffected by nearby subsidence activity. Each lake likely is underlain by several piping features rather than one large subsidence feature.

Florida↗
USGS2008

Cenozoic climate history from seismic reflection and drilling studies on the Antarctic continental margin

Seismic stratigraphic studies and scientific drilling of the Antarctic continental margin have yielded clues to the evolution of Cenozoic climates, depositional paleoenvironments and paleoceanographic conditions. This paper draws on studies of the former Antarctic Offshore Stratigraphy Project and others to review the geomorphic and lithostratigraphic offshore features that give insights into the long-duration (m.y.) and short-term (k.y.) changes that document the great variability of Cenozoic Antarctic paleoenvironments. The lithologic drilling record documents non-glacial (pre-early Eocene) to full-glacial (late Pliocene to Holocene) times, and documents times of cyclic ice-sheet fluctuations at k.y. scales (early Miocene to Pliocene and Holocene). Times of significant change in types and/or amounts of glaciation are also seen in the offshore lithologic record (early Oligocene, mid-Miocene, early Pliocene). Seismic data illustrate large-scale geomorphic features that point to massive sediment erosion and dispersal by ice sheets and paleoceanographic processes (e.g. cross-shelf troughs, slope-fans, rise-drifts). The commonality of these features to East and West Antarctica since late Eocene time points to a continent that has been intermittently covered, partially to completely, by glaciers and ice sheets. The greatest advances in our understanding of paleoenvironments and the processes that control them have been achieved from scientific drilling, and future progress depends on a continuation of such drilling.

Book chapter↗
USGS2004

Precise relocation of earthquakes following the 15 June 1991 eruption of Mount Pinatubo (Philippines)

The 15 June 1991 climactic eruption of Mount Pinatubo (Philippines) was followed by intense seismicity that remained at a high level for several months. We located 10,839 events recorded between 1 July and mid-December 1991. In contrast to the preeruptive seismicity which was focused in two groups below the summit area, posteruptive events were widely distributed below and around the volcano. The classification of the events indicates the presence of several large multiplets, and the application of relative relocation techniques to the similar events by calculating high-precision delays between traces outlines a number of clear seismogenic structures. We used different methods to confirm the validity of our results; these tests indicate that reliable features can be detected with a small monitoring network. While the main cluster of activity can be attributed to an intrusive process starting from below the 15 June crater, the volcanic origin of the seismic activity in the other areas is more difficult to establish. Away from the summit, relocations define streaks or planes which are oriented predominantly southwest-northeast, with in several cases the presence of northwest-southeast conjugate structures. Most of the composite focal mechanisms that we could determine indicate predominantly strike-slip, right-lateral faulting. Our results indicate that most of the seismicity that occurred after the 15 June eruption is related to the east-west regional compressional stress field related to the subduction. We suggest that the regional stress field induces seismicity along new or preexisting faults in the medium surrounding the volcano where the stress field was locally disturbed by the volcanic eruption.

Mount Pinatubo↗
USGS2018

Integrate urban‐scale seismic hazard analyses with the U.S. National Seismic Hazard Model

For more than 20 yrs, damage patterns and instrumental recordings have highlighted the influence of the local 3D geologic structure on earthquake ground motions (e.g., M "> M M 6.7 Northridge, California, Gao et al. , 1996 ; M "> M M 6.9 Kobe, Japan, Kawase, 1996 ; M "> M M 6.8 Nisqually, Washington, Frankel, Carver, and Williams, 2002 ). Although this and other local‐scale features are critical to improving seismic hazard forecasts, historically they have not been explicitly incorporated into the U.S. National Seismic Hazard Model (NSHM, national model and maps), primarily because the necessary basin maps and methodologies were not available at the national scale. Instead, the U.S. Geological Survey (USGS), its partners, and external groups developed urban seismic hazard maps (urban models and maps) that consider detailed site effects in local areas (e.g., Wong et al. , 2002 ; Cramer et al. , 2006 ; Frankel et al. , 2007 ; Graves et al. , 2011 ). The disconnect between the urban and national hazard models, however, means that the national models, which underlie U.S. building codes and other applications, do not make use of all of the scientific results informing earthquake ground‐shaking hazards. We recommend that future U.S. national seismic hazard assessment make use of all available regional information, including that in urban models. In this column, we describe the roles of and differences between the urban and national models, and discuss the obstacles to and benefits of integrating the urban models with the national model. Future progress on this issue will require further coordination and implementation efforts between the USGS and external partners.

Seismological Research Letters↗
USGS2009

The source of infrasound associated with long-period events at mount St. Helens

During the early stages of the 2004-2008 Mount St. Helens eruption, the source process that produced a sustained sequence of repetitive long-period (LP) seismic events also produced impulsive broadband infrasonic signals in the atmosphere. To assess whether the signals could be generated simply by seismic-acoustic coupling from the shallow LP events, we perform finite difference simulation of the seismo-acoustic wavefield using a single numerical scheme for the elastic ground and atmosphere. The effects of topography, velocity structure, wind, and source configuration are considered. The simulations show that a shallow source buried in a homogeneous elastic solid produces a complex wave train in the atmosphere consisting of P/SV and Rayleigh wave energy converted locally along the propagation path, and acoustic energy originating from , the source epicenter. Although the horizontal acoustic velocity of the latter is consistent with our data, the modeled amplitude ratios of pressure to vertical seismic velocity are too low in comparison with observations, and the characteristic differences in seismic and acoustic waveforms and spectra cannot be reproduced from a common point source. The observations therefore require a more complex source process in which the infrasonic signals are a record of only the broadband pressure excitation mechanism of the seismic LP events. The observations and numerical results can be explained by a model involving the repeated rapid pressure loss from a hydrothermal crack by venting into a shallow layer of loosely consolidated, highly permeable material. Heating by magmatic activity causes pressure to rise, periodically reaching the pressure threshold for rupture of the "valve" sealing the crack. Sudden opening of the valve generates the broadband infrasonic signal and simultaneously triggers the collapse of the crack, initiating resonance of the remaining fluid. Subtle waveform and amplitude variability of the infrasonic signals as recorded at an array 13.4 km to the NW of the volcano are attributed primarily to atmospheric boundary layer propagation effects, superimposed upon amplitude changes at the source.

Washington↗
USGS1994

Applications of shallow high-resolution seismic reflection to various environmental problems

Shallow seismic reflection has been successfully applied to environmental problems in a variety of geologic settings. Increased dynamic range of recording equipment and decreased cost of processing hardware and software have made seismic reflection a cost-effective means of imaging shallow geologic targets. Seismic data possess sufficient resolution in many areas to detect faulting with displacement of less than 3 m and beds as thin as 1 m. We have detected reflections from depths as shallow as 2 m. Subsurface voids associated with abandoned coal mines at depths of less than 20 m can be detected and mapped. Seismic reflection has been successful in mapping disturbed subsurface associated with dissolution mining of salt. A graben detected and traced by seismic reflection was shown to be a preferential pathway for leachate leaking from a chemical storage pond. As shown by these case histories, shallow high-resolution seismic reflection has the potential to significantly enhance the economics and efficiency of preventing and/or solving many environmental problems.

Journal of Applied Geophysics↗
USGS2020

Combined seismic and geodetic analysis before, during and after the 2018 Mt. Etna eruption

In December 2018, Etna volcano experienced one of the largest episodes of unrest since the installation of geophysical monitoring networks in 1970. The unrest culminated in a short eruption with a small volume of lava erupted, a significant seismic crisis and deformation of the entire volcanic edifice of magnitude never recorded before at Mount Etna. Here we describe the evolution of the 2018 eruptive cycle from the analysis of seismic and geodetic data collected in the months preceding, during, and following the intrusion. We model the space-time evolution of high-rate deformation data starting from the active source previously identified from deformation data and the propagation of seismicity in a 3-D velocity model. The intrusion model suggests emplacement of two dikes: a smaller dike located beneath the eruptive fissure and a second, deeper dike between 1 and 5 km below sea level that opened ~2 m. The rise and eruption of magma from the shallower dike did not interrupt the pressurization of a long-lasting deeper reservoir (~6 km) that induced continuous inflation and intense deformation of the eastern flank. Shortly after the intrusion, on 26 December 2018, a M L 4.8 earthquake occurred near Pisano, destroying buildings and roads in two villages. We propose a time-dependent intrusion model that supports the hypothesis of the inflation inducing flank deformation and that this process has been active since September 2018.

Mt. Etna↗
USGS2005

Moment tensor inversion of ground motion from mining-induced earthquakes, Trail Mountain, Utah

A seismic network was operated in the vicinity of the Trail Mountain mine, central Utah, from the summer of 2000 to the spring of 2001 to investigate the seismic hazard to a local dam from mining-induced events that we expect to be triggered by future coal mining in this area. In support of efforts to develop groundmotion prediction relations for this situation, we inverted ground-motion recordings for six mining-induced events to determine seismic moment tensors and then to estimate moment magnitudes M for comparison with the network coda magnitudes Mc. Six components of the tensor were determined, for an assumed point source, following the inversion method of McGarr (1992a), which uses key measurements of amplitude from obvious features of the displacement waveforms. When the resulting moment tensors were decomposed into implosive and deviatoric components, we found that four of the six events showed a substantial volume reduction, presumably due to coseismic closure of the adjacent mine openings. For these four events, the volume reduction ranges from 27% to 55% of the shear component (fault area times average slip). Radiated seismic energy, computed from attenuation-corrected body-wave spectra, ranged from 2.4 ?? 105 to 2.4 ?? 106 J for events with M from 1.3 to 1.8, yielding apparent stresses from 0.02 to 0.06 MPa. The energy released for each event, approximated as the product of volume reduction and overburden stress, when compared with the corresponding seismic energies, revealed seismic efficiencies ranging from 0.5% to 7%. The low apparent stresses are consistent with the shallow focal depths of 0.2 to 0.6 km and rupture in a low stress/low strength regime compared with typical earthquake source regions at midcrustal depths.

Bulletin of the Seismological Society of America↗
USGS1979

Microearthquake networks and earthquake prediction

A microearthquake network is a group of highly sensitive seismographic stations designed primarily to record local earthquakes of magnitudes less than 3. Depending on the application, a microearthquake network will consist of several stations or as many as a few hundred . They are usually classified as either permanent or temporary. In a permanent network, the seismic signal from each is telemetered to a central recording site to cut down on the operating costs and to allow more efficient and up-to-date processing of the data. However, telemetering can restrict the location sites because of the line-of-site requirement for radio transmission or the need for telephone lines. Temporary networks are designed to be extremely portable and completely self-contained so that they can be very quickly deployed. They are most valuable for recording aftershocks of a major earthquake or for studies in remote areas.

Earthquake Information Bulletin (USGS)↗
USGS1981

The geology of Hess Rise, central north Pacific Ocean

The geology of Hess Rise is interpreted from studies of morphology and structure, igneous petrology, sediment lithofacies, seismic stratigraphy, and paleoenvironments. Hess Rise probably formed at a triple junction, between 116 and 95 m.y. ago, in the equatorial zone of the southern hemisphere. A large archipelago, formed by volcanism, subsequently moved northward as part of the Pacific Plate and subsided. A Late Cretaceous-early Tertiary tectonic and(or) volcanic event created oceanic islands and probably influenced the normal subsidence rate. The lithofacies reflect both the northward movement under various surface water masses of differing productivity and the subsidence history, which influenced calcareous-sediment accumulation and dissolution. Hess Rise is informally divided into three physiographic or morphologic provinces: northern Hess Rise, the central platform, and southern Hess Rise. Southern Hess Rise is further subdivided into the northwestern ridge, Mellish Bank, and the eastern ridge. The structure of Hess Rise is dominated by normal faults which bound grabens and horsts. The west and south boundaries of the rise are fault scarps. Most faults trend about 345 ° on northern Hess Rise, 039° to 060° on southern Hess Rise. Igneous rocks are tholeiitic basalt on northern Hess Rise at Site 464, and trachyte at Site 465 on southern Hess Rise. The tholeiite probably is not a mid-ocean-ridge basalt; more likely it is a transitional basalt erupted during off-ridge volcanism. The trachyte at Site 465 was erupted either in shallow water or subaerially; it represents a late-stage differentiate of alkali-basalt magma. Geochemically, the trachyte is similar to late-stage eruptives on many oceanic islands. The alkali basalt clasts in Upper Cretaceous calcareous ooze at Site 466 probably record a Late Cretaceous-early Tertiary tectonic event that brought old crust above sea level, where it was eroded; alternatively, volcanism may have built the islands. Cores recovered from Holes 464, 310, 310A, 465,465A, and 466 (north to south) show three major lithofacies: limestone of mid-Cretaceous (Albian-Cenomanian) age, calcareous ooze and chalk of Late Cretaceous to Quaternary age, and pelagic clay that is Late Cretaceous to Quaternary in age at Site 464 and Eocene in age at Site 310. Minor lithofacies in the carbonates and pelagic clay include siliceous microfossils and chert. The limestone contains as much as 8.6% by weight of organic carbon. The organic-carbon-rich sediments apparently were deposited along the sides of Hess Rise, above the CCD and within a mid-water oxygen minimum, as the rise crossed the highly productive equatorial divergence. Surface sediments are influenced by present-day near-surface water-mass productivity and the depth of deposition. Above about 3000 meters calcareous ooze is deposited, and below about 4000 meters pelagic clay is dominant. Between 3000 meters and 4000 meters, greater clay contents correlate with increased water depth. The seismic stratigraphy interpreted from single-channel seismic-reflection records can be correlated with the lithofacies. Three major acoustic units apparently vary in age, composition, thickness, and geographic distribution. The lowest acoustic unit generally is acoustically opaque; it corresponds to volcanic basement. In places, however, the top parts of this unit may be well-lithified limestone or chert. The middle acoustic unit, in general strongly layered acoustically, is Albian and Cenomanian limestone, chalk, and chert on southern Hess Rise and Albian and Cenomanian chert, chalk, limestone, and clay stone on northern Hess Rise. The middle unit probably is Albian to at least lower Tertiary chalk, chert, and nannofossil ooze on the central platform. The upper acoustic unit is weakly stratified to acoustically transparent and corresponds to Upper Cretaceous to Recent nannofossil ooze on southern Hess Rise and to pelagic clay of the same age on northern Hess Rise. On the central platform, it probably corresponds to Eocene to Recent nannofossil ooze, marl, and pelagic clay

Initial Reports of the D.S.D.P.↗
USGS1984

Common-depth-point seismic-reflection survey on the Mississippi River in the vicinity of Alton, Illinois

A high-resolution seismic-reflection survey was run on the Mississippi River using marine common-depth-point (CDP), digital, and multichannel techniques in an attempt to locate a possible fault that had been interpreted by others from magnetic and land seismic data to be in the vicinity of Lock and Dam No. 26. The geologic feature is commonly known in the literature as the Cap au Ores Fault or Cap au Ores Uplift. A second objective was to demonstrate that the marine seismic system could be used successfully for making relatively deep penetration seismic-reflection surveys on inland waterways with shallow water depths (less than 10 meters) and organic sediments on the waterway bottom. The survey was made by agreement between the U.S. Army Corps of Engineers (COE), St. Louis District, and the U.S. Geological Survey (USGS). The USGS used a specially equipped, trailerable, 38-ft aluminum boat for the work. The quality of the survey records can be classified as from fair to good. Seismic-reflection events could be interpreted in the upper 200 ms of data. Five reflecting horizons with little slope were recorded between Mississippi River Miles 195 and 210 without any evidence of faulting of the magnitude (50 to 100 m) suggested by the Shannon and Wilson report (1980). Minor faulting and buried channeling could be interpreted. This survey demonstrated the value and economy of obtaining deeper penetration seismic data on inland water using an offshore marine, CDP, digital, multichannel seismic-reflection system. On-going modifications in field procedure, equipment configuration, equipment, and data processing will further improve the results obtained with such equipment. Future marine-style surveys conducted with the newly available hydraulic waterguns (80-, 200-, and 400-in 3 ) would probably produce data that are a significant improvement over data obtained with the 15-in 3 pneumatic watergun and 40-in 3 airgun used during the 1980 survey.

Illinois, Missouri↗
USGS2022

Earthquake scenario selection for portfolio holders in CEUS: A case study with Oklahoma DOT

Portfolio managers of spatially distributed assets in the central and eastern United States (CEUS) and other low- to moderate seismic hazard regions require scenario-based seismic risk assessment for the purpose of emergency management and planning. Uncertainties regarding the long-term seismicity of the region, unknown faults, and limited historical records complicate the selection of an earthquake scenario. Through a case study with the Oklahoma Department of Transportation (ODOT) and their portfolio of bridges, we look at one such exercise, which consists of two scenario earthquakes: one scenario earthquake selected from the U.S. Geological Survey Building Seismic Safety Commission (BSSC) scenario catalog, a magnitude (M) 7.2 event on the Meers fault, and a second aftershock selected by the consequence-driven earthquake scenario selection (Co-DESS) method. The latter is driven by ODOT’s desired service actions to be included during the earthquake drill; in this case, we identify an earthquake that is likely to trigger inspections for bridges across multiple districts, thereby testing not only inspection protocols but also coordination efforts between district groups. We find that the Co-DESS selected event is smaller in magnitude and offers different geographical options than selection through conventional selection methods, while still meeting necessary consequences for an effective earthquake exercise.

Conference Paper↗
USGS2009

Multi-segment earthquakes and tsunami potential of the Aleutian megathrust

Large to great earthquakes and related tsunamis generated on the Aleutian megathrust produce major hazards for both the area of rupture and heavily populated coastlines around much of the Pacific Ocean. Here we use paleoseismic records preserved in coastal sediments to investigate whether segment boundaries control the largest ruptures or whether in some seismic cycles segments combine to produce earthquakes greater than any observed since instrumented records began. Virtually the entire megathrust has ruptured since AD1900, with four different segments generating earthquakes >M8.0. The largest was the M9.2 great Alaska earthquake of March 1964 that ruptured ???800 km of the eastern segment of the megathrust. The tsunami generated caused fatalities in Alaska and along the coast as far south as California. East of the 1964 zone of deformation, the Yakutat microplate experienced two >M8.0 earthquakes, separated by a week, in September 1899. For the first time, we present evidence that earthquakes ???900 and ???1500 years ago simultaneously ruptured adjacent segments of the Aleutian megathrust and the Yakutat microplate, with a combined area ???15% greater than 1964, giving an earthquake of greater magnitude and increased tsunamigenic potential. ?? 2008 Elsevier Ltd. All rights reserved.

Quaternary Science Reviews↗
USGS2022

Empirical map-based nonergodic models of site response in the greater Los Angeles area

We develop empirical estimates of site response at seismic stations in the Los Angeles area using recorded ground motions from 414 M 3–7.3 earthquakes in southern California. The data are from a combination of the Next Generation Attenuation‐West2 project, the 2019 Ridgecrest earthquakes, and about 10,000 newly processed records. We estimate site response using an iterative mixed‐effects residuals partitioning approach, accounting for azimuthal variations in anelastic attenuation and potential bias due to spatial clusters of colocated earthquakes. This process yields site response for peak ground acceleration, peak ground velocity, and pseudospectral acceleration relative to a 760 m/s shear‐wave velocity ( ⁠ V S "> V V S "> s ) reference condition. We employ regression kriging to generate a spatially continuous site response model, using the linear site and basin terms from Boore et al. (2014) as the background model, which depend on V S "> V V S "> s V S "> 30 ⁠ and depth to the 1 km/s V S "> V V S "> s isosurface. This is different from past approaches to nonergodic models, in which spatially varying coefficients are regressed. We validate the model using stations in the Community Seismic Network (CSN) that are in the middle of our model spatial domain but were not considered in model development, finding strong agreement between the interpolated model and CSN data for long periods. Our model could be implemented in regional seismic hazard analyses, which would lead to improvements especially at long return periods. Our site response model also has potential to improve both ground‐motion accuracy and warning times for the U.S. Geological Survey ShakeAlert earthquake early warning (EEW) system. For a point‐source EEW simulation of the 1994 M 6.7 Northridge earthquake, our model produces ground motions more consistent with the ground‐truth ShakeMap and would alert areas with high population density such as downtown Los Angeles at lower estimated magnitudes (i.e., sooner) than an ergodic model for a modified Mercalli intensity 4.5 alerting threshold.

California↗
USGS2019

Seismic design and hazard maps: Before and after

The 1994 Northridge earthquake generated world-record ground motions. At the time, the horizontal peak ground acceleration of 1.8 g measured by a seismometer in Tarzana was the largest ever. The same is true of the peak ground velocity of 148 cm/s measured in Granada Hills. Both measurements were within approximately 15 km of the source of the earthquake; they were also near most of the damage described in other articles of this series. Consequently, the near-source design forces from the seismic zone maps in the Uniform Building Code (UBC) were increased. From the 1994 to 1997 editions, acceleration- and velocity-related near-source factors were introduced. The factors increased the design forces in Zone 4, already the highest seismic zone, by a multiplier as large as 2.0. More enduringly, generational changes were made to the seismic design maps in the NEHRP Recommended Seismic Provisions for New Buildings and Other Structures . The NEHRP maps were – and continue to be – adopted into the International Building Code (IBC), which supplanted the UBC and other model building codes. As described below, the changes to the NEHRP maps took advantage of another post-Northridge change: the modern generation of U.S. Geological Survey (USGS) National Seismic Hazard Maps.

Structure↗
USGS2018

Decomposing leftovers: Event, path, and site residuals for a small magnitude ANZA region GMPE

Ground‐motion prediction equations (GMPEs) are critical elements of probabilistic seismic hazard analysis (PSHA), as well as for other applications of ground motions. To isolate the path component for the purpose of building nonergodic GMPEs, we compute a regional GMPE using a large dataset of peak ground accelerations (PGAs) from small‐magnitude earthquakes ( ⁠0.5 ≤ M ≤ 4.5 with > 10,000 events, yielding ∼120,000 recordings) that occurred in 2013 centered around the ANZA seismic network (hypocentral distances ≤180 km⁠ ) in southern California. We examine two separate methods of obtaining residuals from the observed and predicted ground motions: a pooled ordinary least‐squares model and a mixed‐effects maximum‐likelihood model. Whereas the former is often used by the broader seismological community, the latter is widely used by the ground‐motion and engineering seismology community. We confirm that mixed‐effects models are the preferred and most statistically robust method to obtain event, path, and site residuals and discuss the reasoning behind this. Our results show that these methods yield different consequences for the uncertainty of the residuals, particularly for the event residuals. Finally, our results show no correlation (correlation coefficient [CC] <0.03⁠ ) between site residuals and the classic site‐characterization term V S30 ⁠ , the time‐averaged shear‐wave velocity in the top 30 m at a site. We propose that this is due to the relative homogeneity of the site response in the region and perhaps due to shortcomings in the formulation of V S30 ⁠ and suggest applying the provided PGA site correction terms to future ground‐motion studies for increased accuracy.

California↗