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 379 records · Page 21Linked to original sources
USGS2024

Uncertainty and spatial correlation in station measurements for mb magnitude estimation

The body‐wave magnitude (⁠⁠) is a long‐standing network‐averaged, amplitude‐based magnitude used to estimate the magnitude of seismic sources from teleseismic observations. The U.S. Geological Survey National Earthquake Information Center (NEIC) relies on in its global real‐time earthquake monitoring mission. Although waveform modeling‐based moment magnitudes are the modern standard to characterize earthquake size, is important because (1) in many cases, waveform modeling is not possible (e.g., low signal‐to‐noise events), (2) is applicable over a broad range of magnitudes, ∼M 4–7, and (3) there is a many decades‐long history of estimating magnitudes. We use the NEIC Preliminary Determination of Epicenters earthquake catalog to investigate the uncertainty in NEIC station measurements. We show that measurements are spatially correlated, which can bias event ⁠, and we describe an empirical relation between this spatial correlation and station‐to‐station distance. We further describe an approach to mitigate bias from the spatial correlation. Accounting for the spatial covariance of observations can change the event from −0.15 to 0.07 units (10th to 90th percentile) for smaller events (⁠⁠). These smaller events have the largest standard deviations ranging from 0.05 to 0.15 units (10th to 90th percentile).

The Seismic Record↗
USGS2024

The 17 January 1994 Northridge, California, earthquake: A retrospective analysis

The 17 January 1994 Northridge, California, earthquake was a watershed event, with far-reaching societal and scientific impacts. The earthquake, which occurred in the early days of both broadband seismic networks and the Internet, spurred advances in seismic monitoring, real-time systems, and development of data products. Motivated by the 30th anniversary of the earthquake, we present a brief retrospective of the earthquake and its impact, and reconsider both ground motions and the aftershock distribution using modern tools and best-available data. With improvements in instrumentation and analysis methodology, recent earthquakes continue to reveal increasing complexity of ground motions, fault systems, and earthquake ruptures. Even in the absence of data from state-of-the art instrumentation, a retrospective consideration of ground motion data from the Northridge earthquake reveals complexities beyond what could be characterized (and modeled) thirty years ago. Aftershock relocations for both the 1971 Sylmar and 1994 Northridge earthquakes also reveal an updated view of fault complexity. Our study does provide a cautionary tale regarding legacy data sets and research results that are not easily accessible, which can result in discrepancies between catalog data and products from best-available science. We also briefly describe outreach products produced as part of the anniversary commemoration.

California↗
USGS2024

Preliminary observations of the April 5th, 2024, Mw4.8 New Jersey earthquake

On 5 April 2024, 10:23 a.m. local time, a moment magnitude 4.8 earthquake struck Tewksbury Township, New Jersey, about 65 km west of New York City. Millions of people from Virginia to Maine and beyond felt the ground shaking, resulting in the largest number (>180,000) of U.S. Geological Survey (USGS) “Did You Feel It?” reports of any earthquake. A team deployed by the Geotechnical Extreme Events Reconnaissance Association and the National Institute of Standards and Technology documented structural and nonstructural damage, including substantial damage to a historic masonry building in Lebanon, New Jersey. The USGS National Earthquake Information Center reported a focal depth of about 5 km, consistent with a lack of signal in Interferometric Synthetic Aperture Radar data. The focal mechanism solution is strike slip with a substantial thrust component. Neither mechanism’s nodal plane is parallel to the primary northeast trend of geologic discontinuities and mapped faults in the region, including the Ramapo fault. However, many of the relocated aftershocks, for which locations were augmented by temporary seismic deployments, form a cluster that parallels the general northeast trend of the faults. The aftershocks lie near the Tewksbury fault, north of the Ramapo fault.

New Jersey↗
USGS2024

The value of hyperparameter optimization in phase-picking neural networks

The effectiveness of using neural networks for picking seismic phase arrival times has been demonstrated through several case studies, and seismic monitoring programs are starting to adopt the technology into their workflows. However, published models were designed and trained using rather arbitrary choices of hyperparameters, limiting their performance. In this study, we use phase picks from both routine and template-matching analyses from multiple regions (Ridgecrest, California; Kilauea, Hawaii; Yellowstone, Wyoming-Montana-Idaho) to test a hyperparameter optimization scheme for phase-picking neural networks and to evaluate their performance. We show that a published model, namely PhaseNet (Zhu and Beroza, 2019), can be simplified and improved with reasonable effort and there are preferred choices of hyperparameters that increase the performance. We also show that models optimized based on the arrival times reported in routine event catalogs consistently perform well when picking arrival times of smaller events, which is crucial for certain tasks from microseismicity to explosion monitoring.

California, Hawaii, Idaho, Montana, Wyoming↗
USGS2025

Reply to, “Comment on ‘The 1886 Charleston, South Carolina, earthquake: Relic railroad offset reveals rupture,’ by Roger Bilham and Susan E. Hough”

We welcome this opportunity to respond to Pratt et al. (2024) (hereinafter P24). Bilham and Hough (2023) proposed a “first-cut” elastic deformation model for the 1886 earthquake, a quantitative source model constrained by identified coseismic constraints. A key observation was the measurement of a lateral offset of a railroad line south of Summerville, leading to a model with predominately dextral slip and minor convergence, from which we concluded that active faulting had raised the Penholoway Marine Terrace >6 m since ∼770 ka. P24 questioned these constraints and proposed an alternative rupture model with predominantly reverse slip. This alternative model is neither consistent with coseismic constraints nor with other geophysical data. In a revised model presented here, we recognize that uplift of the Penholoway Terrace is confined to the eastern edge of the terrace, which we conclude results from active folding and tectonic transpression centered on the dextral fault that offset the railroad in 1886.

South Carolina↗
USGS2024

Reduced injection rates and shallower depths mitigated induced seismicity in Oklahoma

The proximity of wastewater disposal to the Precambrian basement is a critical factor influencing induced earthquake rates in the Central United States, but the impact of reducing injection depths has not been widely demonstrated. Beginning in 2015, state regulatory efforts in Oklahoma and Kansas mandated that wells injecting into the lower Arbuckle Group, a basal sedimentary unit, be backfilled with cement (i.e. “plugged back”) so that they inject into shallower formations. This plug back activity gives us a unique opportunity to investigate the relationship between injection depth and induced seismicity rate. To evaluate the impact that decreased injection rates and plug backs had on the seismicity rates, we create a suite of rate-state earthquake models. Observed seismicity rates are best fit when only lower Arbuckle volumes are considered, suggesting the lower Arbuckle injectors were primarily responsible for the seismicity and that plug backs were effective at isolating the injected volumes to shallower formations. Our models demonstrate that if these wells had not been plugged back, seismicity rates would be multiple times larger than they are today. We find that the combination of well plug backs and injection volume decreases can be effective strategies for reducing induced seismicity rates.

Oklahoma↗
USGS2025

Remote single-station seismic monitoring of the July–October 2022 earthquake swarm at Tau volcano, American Samoa

From July to October 2022, a non-eruptive volcanic earthquake swarm occurred within ~15 km of Taʻū Island, located in eastern American Samoa. Felt reports from local residents were the only available information about the swarm when it started, as American Samoa lacked a seismic monitoring network. We developed a consistent single-station catalog for the entire swarm, using seismic data from the nearest station IU.AFI, ~250 km away. We applied the EQTransformer deep-learning model (Mousavi et al., 2020), automatically picking Pn and Sn arrivals on IU.AFI continuous data. We retained only events with Sn-Pn times of 22.5–25 seconds, consistent with the expected locations based on felt reports, then detected smaller swarm events with subsequent template-matching. This single-station catalog characterized the swarm’s onset and escalation to peak activity before a multi-agency field response team installed a local seismic network in mid-August 2022. This permanent seismic network captured the swarm’s decline. EQTransformer identified short S-P times on the first two locally deployed seismometers, both Raspberry Shake sensors, to constrain the swarm’s distance from Taʻū Island. Modern seismological processing methods, combined with basic observations such as felt reports, can quickly contribute useful information during an earthquake response in a poorly monitored region.

American Somoa↗
USGS2025

Long-term trends in microseismicity during operational shut-ins at the Coso Geothermal Field, California

Pausing injection and production can lead to induced seismicity in a variety of settings, with some of the largest events occurring during these so-called shut-ins. In geothermal fields, shut-ins are periodically conducted for maintenance on wells and surface infrastructure, thereby offering recurring means of estimating stress changes in the subsurface that lead to increased seismicity rates. Here, we examine daily production and microseismicity data from the Coso Geothermal Field (CGF) in California between 1996 and 2010 to investigate the repetitive trends of operational shut-in microseismicity. Using 24 local seismic stations, we first analyze spatial and temporal trends of over 60,000 earthquakes with magnitudes between -0.4 to 3.8. We find that the northern region exhibits no significant seismicity changes during shut-ins, whereas the rest of the field experiences induced seismicity during almost every shut-in with an increasing intensity towards the southern and eastern portions of the field, highlighting local differences in stress within the CGF. Additionally, we cluster the seismicity using waveform cross-correlation, revealing several earthquake clusters primarily occurring during shut-in periods. These observations suggest that certain fracture and fault sections respond quicker to changes in pore pressure and poroelastic stresses within the geothermal system, possibly highlighting main fluid pathways.

California↗
USGS2025

Quantitative evaluations of earthquake early warning performance using “Did You Feel It?” and post-alert surveys

We examine responses to the U.S. Geological Survey’s “Did You Feel It?” (DYFI) survey and its companion earthquake early warning (EEW) questionnaire to assess the performance of the U.S. ShakeAlert EEW system directly from the alert recipients’ perspectives. ShakeAlert rapidly detects earthquakes and develops alert information, but as official alert delivery partners issue these alerts, it is thus difficult to determine how many people were alerted and when. We investigate DYFI reports for six California earthquakes that had EEW alerts and substantial responses to the DYFI EEW questionnaire. Comparisons of ShakeAlert predictions to reported intensities demonstrate that magnitude estimation accuracy is not necessarily indicative of ground-motion prediction accuracy. Perceived warning time distributions indicate that estimating maximum-expected warning times using the S-wave arrival is a reasonable assumption when discussing public EEW performance. However, we also find many reports of shorter warning times, late alerts, and missed alerts than expected based on ShakeAlert publication times, suggesting alert delivery latencies are substantial and highly variable. The novelty of our analysis is that we demonstrate that the DYFI EEW survey provides useful EEW efficacy information—independent of the specific alerting pathway—that can be used to inform our choices for conveying EEW performance.

California↗
USGS2025

Neotectonic origins for the Meadow Bank scarp, Wabash Valley seismic zone USA

The Meadow Bank scarp (MBS) in southeastern Illinois is a linear geomorphic expression, ∼10 km long and ∼8 m high above a relatively flat landscape. It parallels an underlying northeast‐oriented Late‐Precambrian–Early‐Cambrian structural fabric, called the Wabash Valley fault zone, and is within an area of modern, historic, and paleo seismicity, called the Wabash Valley seismic zone. In addition, the MBS acts as a boundary of the Wabash River floodplain, as well as Pleistocene glacial outwash channels, which show evidence of frequent outburst flood events. To better understand the MBS’s equivocal origin in this complex geologic environment, we acquired a 917‐m‐long seismic‐reflection survey across its axis to assess the subsurface geologic configuration. The resultant image indicates a complex set of faults that antiformally fold and displace the top of Paleozoic bedrock by ∼12 m across the survey. Moreover, fault and/or fold deformation extends into the shallowest imaged Quaternary strata at ∼6 m below the ground surface. This suggests the MBS origin is related to underlying Quaternary reactivated, positively inverted faults rather than exclusively to glacial outburst flood erosion. These results provide rare paleoearthquake spatial constraints for central U.S. regional seismic hazard consideration.

Illinois↗
USGS2026

Only very strong shaking can break a tree

A handful of credible accounts, both recent and historical, describe trees being broken during earthquakes. Apart from landslides and tsunamis, such accounts are uncommon, suggesting that trees break only when ground acceleration is extreme or when trees are especially vulnerable to breakage. Here, I consider the question: what ground acceleration is required to break a living tree? I present simplified first‐principles calculations, considering the properties of common trees, focusing on tall, thin trees that can be approximated as uniform vertical cantilevers. The results suggest that, possibly excepting tall trees with heavy canopies, standing trees can only be broken by shaking approaching or exceeding 1 g , as anticipated by the rarity of documented accounts of snapped trees. I briefly consider several well‐documented instances of tree damage during strong earthquakes, drawing inferences about local ground motions.

The Seismic Record↗
USGS2026

Erratum to earthquake recurrence estimates for northern Caribbean faults from combinatorial optimization

In the original article (ten Brink and Geist, 2025), we derived on‐fault magnitude–frequency relationships with uncertainties from a combinatorial optimization method that assigns individual earthquakes from a regional catalog to the different faults. Our model used a magnitude versus area relationship from Shaw (2023). His relationship was published with an error, M=f(1/√A)⁠, instead of M=f(1/A) (Shaw, 2025). Here, we recalculate the magnitude–frequency relationship with uncertainties using the Shaw (2025) corrected formula. The remaining parameters of the model, the model itself, and the number of repeated optimization runs have not changed. Figures 1 and 2 show the corrected results and replace figures 3 and 4 in the original article.

The Seismic Record↗
USGS2026

Investigating high‐resolution spatial variability of site response in the Los Angeles Basin using a dense sensor network

We characterize seismic site response across a densely instrumented region of the Los Angeles Basin by combining data from the Community Seismic Network (CSN) and regional networks in the Advanced National Seismic System (ANSS). Mixed-effects regression is used to partition ground-motion model residuals and compute pseudo-spectral acceleration site residuals (η s ) from the NGA-West2 ground motion database and a CSN dataset. For short periods ( T < 0.3 s), CSN and ANSS site residuals agree closely, while observed differences at longer periods are likely related to signal-to-noise limitations of the CSN sensors rather than sensor accuracy. Correlations between CSN and ANSS site residuals are high, and the dense CSN network reveals localized amplification not resolved by the ANSS network. Along the northeastern basin margin, PSA at 1 second reaches amplitudes up to approximately twice those predicted for reference rock conditions, consistent with basin-edge amplification effects. Systematic variations in site response are observed across major fault zones, particularly between [BSAS1.1]areas north of the Lower Elysian Park fault and between [BSAS2.1]the Lower Elysian Park and Puente Hills faults. These findings demonstrate the value of dense seismic networks, such as the CSN, for resolving localized amplification and improving seismic hazard assessment in urban sedimentary basins.

California↗
USGS1976

Effects of local geological conditions in the San Francisco Bay region on ground motions and the intensities of the 1906 earthquake

Measurements of ground motion generated by nuclear explosions in Nevada have been completed for 99 locations in the San Francisco Bay region, California. The recordings show marked amplitude variations in the frequency band 0.25 to 3.0 Hz that are consistently related to the local geological conditions of the recording site. The average spectral amplifications observed for vertical and horizontal ground motions are, respectively: (1, 1) for granite, (1.5, 1.6) for the Franciscan Formation, (3.0, 2.7) for the Santa Clara Formation, (3.3, 4.4) for alluvium, and (3.7, 11.3) for bay mud. Spectral amplification curves define predominant ground frequencies in the band 0.25 to 3.0 E for bay mud sites and for some alluvial sites. Amplitude spectra computed from recordings of seismic background noise at 50 sites do not generally define predominant ground frequencies. The intensities ascribed to various sites in the San Francisco Bay region for the California earthquake of April 18, 1906, are strongly dependent on distance from the zone of surface faulting and the geological character of the ground. Considering only those sites (approximately one square city block in size) for which there is good evidence for the degree of ascribed intensity, the intensities for 917 sites on Franciscan rocks generally decrease with the logarithm of distance as Intensity = 2.69 - 1.90 log (Distance in kilometers). (1) For sites on other geological units, intensity increments, derived from this empirical relation, correlate strongly with the Average Horizontal Spectral Amplifications (AHSA) according to the empirical relation Intensity Increment = 0.27 + 2.70 log (AHSA). (2) Average intensity increments predicted for the various geological units are −0.3 for granite, 0.2 for the Franciscan Formation, 0.6 for the Great Valley sequence, 0.8 for the Santa Clara Formation, 1.3 for alluvium, and 2.4 for bay mud. The maximum intensity map predicted on the basis of these data delineates areas in the San Francisco Bay region of potentially high intensity for large earthquakes on either the San Andreas fault or the Hayward fault. The map provides a crude form of seismic zonation for the region and may be useful for certain general types of land-use zonation.

California↗
USGS1978

Stamukhi shoals of the Arctic - some observations from the Beaufort Sea

A number of linear shoals, representing pronounced topographic anomalies on the surface of the Arctic shelf, have been studied in the Prudhoe Bay area. These shoals have been referred to in several previous studies. Based on seismic reflection records, Reimnitz et al., (1972a), stated that the shoals are constructional features younger than the post-Wisconsin transgression. Large chunks of grounded ice have frequently been seen on the linear shoals of the inner shelf. These ice chunks form barriers parallel to the shore (Reimnitz, et al. 1972b). Reimnitz and Barnes (1974) considered the lack of gravel concentrations on the shoals to be evidence that the stream of pack ice drifting past northern Alaska carries very little gravel today. Lewellen (1977) referred to the linear shoals as submerged barrier islands, while Reimnitz, et al., (1977b) pointed out that, although similar in shape to barrier islands, the linear shoals are very different in composition and do not appear to represent drowned harrier islands. They show that the shoals localize the formation of major shear and pressure events in the ice, which in turn cause the formation of linear belts of deformed and grounded ice. Today the shoals appear to be migrating under the influence of ice-bottom interaction, and indeed may have formed in response to ice-bottom interaction within the "stamukhi zone". The shoals migrate rather slowly and retain their shapes over periods of 25 years, yet control the location and stabilize the outer edge of the floating fast ice zone, and provide shelter for the inner shelf and coast. Reimnitz et al. (1977b) surmised that similar artificial structures might be used to modify the ice environment on the arctic shelf.

Alaska↗
USGS1981

Description of seafloor sediment and preliminary geo-environmental report, Shelikof Strait, Alaska

Shelikof Strait, situated between the Kodiak Island group and the Alaska Peninsula (Fig. 1), is included in OCS oil and gas lease area 60. Environmental geologic studies are being conducted by the U.S. Geological Survey prior to the scheduled September 1981 sale date. Seismic-reflection records, collected with 40- to 95-cm3 airgun, 800-joule minisparker, 800-joule boomer, and 3.5-and 12-kilohertz systems and covering over 6,400 km of trackline have been examined to identify geologic conditions at or below the seafloor that might affect petroleum operations. Of the total trackline distance, 865 km were collected in June 1980 aboard the Geological Survey's ship R/V S.P. LEE. The remainder were collected by Nekton Inc. in 1979, under contract to the USGS (Fig. 2). Sediment samples were collected at 42 stations on the June 1980 cruise for geological and geotechnical analysis (Fig. 2). The purpose of this report is to present physical and chemical measurements that have been made on sediment samples and to give a preliminary geo-environmental assessment of Shelikof Strait, based on presently completed analyses of the geophysical records and sediment samples.

Alaska↗