SEARCH · USGS Science

Results for “The Seismological Record”

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

SourceResearch connectionFrom yearThrough year

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

At least 361 records · Page 20Linked to original sources

USGS2014

Rapid earthquake characterization using MEMS accelerometers and volunteer hosts following the M 7.2 Darfield, New Zealand, Earthquake

We test the feasibility of rapidly detecting and characterizing earthquakes with the Quake‐Catcher Network (QCN) that connects low‐cost microelectromechanical systems accelerometers to a network of volunteer‐owned, Internet‐connected computers. Following the 3 September 2010 M 7.2 Darfield, New Zealand, earthquake we installed over 180 QCN sensors in the Christchurch region to record the aftershock sequence. The sensors are monitored continuously by the host computer and send trigger reports to the central server. The central server correlates incoming triggers to detect when an earthquake has occurred. The location and magnitude are then rapidly estimated from a minimal set of received ground‐motion parameters. Full seismic time series are typically not retrieved for tens of minutes or even hours after an event. We benchmark the QCN real‐time detection performance against the GNS Science GeoNet earthquake catalog. Under normal network operations, QCN detects and characterizes earthquakes within 9.1 s of the earthquake rupture and determines the magnitude within 1 magnitude unit of that reported in the GNS catalog for 90% of the detections.

Darfield↗

USGS2019

An investigation of seismicity induced by hydraulic fracturing in the Sichuan basin of China based on data from a temporary seismic network

Hydraulic fracturing has been inferred to trigger the majority of injection‐induced seismicity in the Zhaotong and Changning shale gas field, Sichuan basin of China, in contrast to the Midwestern United States, where massive wastewater disposal has been the dominant triggering mechanism. More than 15,000 earthquakes, with magnitudes ranging up to M w 4.7, were recorded by a temporary network in the Sichuan basin for three years, with a completeness magnitude of M L 1.1. A catalog of earthquakes relocated with code tomoDD, combined with comprehensive injection data during an eight day period, shows that the earthquakes were highly correlated spatiotemporally with hydraulic fracturing activities mostly from a single well pad. Three M L ≥4.0 events occurred during hydraulic fracturing operations from 12 to 19 January 2017, followed by the fourth and largest event, with moment magnitude (⁠M w⁠ ) 4.7, on 28 January. The hypocenters of the four largest events were located in dolomite of Cambrian age, between a shale gas reservoir and the top of the crystalline basement rocks. This was found to be similar to 60% of the smaller earthquakes in this cluster, at depths from 2.5 to 4.0 km.

Sichuan basin↗

USGS2020

Rapid geodetic observations of spatiotemporally varying postseismic deformation following the Ridgecrest earthquake sequence: The U.S. Geological Survey response

The U.S. Geological Survey’s geodetic response to the 4–5 July 2019 (Pacific time) Ridgecrest earthquake sequence comprised primarily the installation and/or reoccupation of Global Navigation Satellite System (GNSS) monumentation. Our response focused primarily on the United States’ Navy’s China Lake Naval Air Weapons Station base (NAWSCL). This focus was because much of the surface rupture occurred on the NAWSCL and because of NAWSCL access restrictions only permitting Federal and State of California personnel. In total, we measured or are still measuring at 24 sites, 14 of which were on the NAWSCL and, as of this writing, operational. The majority of sites were set up as continuous stations logging at either 1 sample per second or 1 sample per 15 s. Two stations were recording a 200 m cross‐rupture aperture starting ∼10 hr after the M 6.4 event, and they recorded the coseismic displacements of the M 7.1. Approximately, 1 hr after the M 7.1 event, two new stations were recording a ∼200 m cross‐rupture aperture of the surface rupture. In the days following, we established the rest of the stations ranging to a distance of ∼15 km from the M 7.1 principal rupture trace. The lack of differential displacement across the M 6.4 rupture during the M 7.1 event suggests that it did not reactivate the M 6.4 plane. The lack of differential cross‐fault displacement for both events suggests that rapid shallow afterslip did not occur at those two locations. The postseismic time series from these stations shows centimeters of horizontal displacement over periods of a few months. They record a mixture of fault‐parallel and fault‐normal displacements that, in conjunction with analysis of more spatially complete Interferometric Synthetic Aperture Radar displacement fields, suggest that both poroelastic and afterslip phenomena occur along the M 6.4 and 7.1 rupture planes. Using preliminary data from these and other regional stations, we also explore the Ridgecrest sequence’s effect on regional GNSS time series and the differentiation of long‐term postseismic motions and secular deformation rates. We find that redefining a common‐mode noise filter using different GNSS stations that are assumed to be unaffected by the earthquakes results in small but systematic differences in the regional velocity field estimate.

Seismological Research Letters↗

USGS2020

Kinematics of fault slip associated with the July 4-6 2019 Ridgecrest, Californai earthquakes sequence

The 2019 Ridgecrest, California, earthquake sequence produced observable crustal deformation over much of central and southern California, as well as surface rupture over several tens of kilometers. To obtain a detailed picture of the fault slip involved in the 4 July M 6.4 foreshock and 6 July M 7.1 mainshock, we combine strong‐motion seismic waveforms with crustal deformation observations to obtain kinematic and static slip models of both events. We sample the regional seismic wavefield for both the foreshock and mainshock with three‐component records from 31 stations of the California Integrated Seismic Network. The deformation observations include Global Positioning System (GPS), Interferometric Synthetic Aperture Radar (InSAR), and borehole strainmeter recordings of the dynamic strain field. These data collectively constrain the kinematic coseismic slip distributions of the events, with measurements variously observing coseismic slip from one event (e.g., seismic waveforms, kinematic solutions from continuous GPS, and strainmeter time series) or coseismic slip from both events combined (InSAR). We find that the foreshock ruptured two separate faults, one with left‐lateral strike slip on a northeast–southwest‐trending fault and the other with right‐lateral strike slip on an orthogonal fault, with unilateral rupture propagation along both. The mainshock ruptured a series of northwest–southeast‐trending faults with right‐lateral strike slip concentrated in the uppermost 6 km with exceptionally low‐rupture velocity averaging 1.0 – 1.5    km / s ">

California↗

USGS2021

Seismic wave propagation and basin amplification in the Wasatch Front, Utah

Ground‐motion analysis of more than 3000 records from 59 earthquakes, including records from the March 2020 M w "> M w 5.7 Magna earthquake sequence, was carried out to investigate site response and basin amplification in the Wasatch Front, Utah. We compare ground motions with the Bayless and Abrahamson (2019 ; hereafter, BA18) ground‐motion model (GMM) for Fourier amplitude spectra, which was developed on crustal earthquake records from California and other tectonically active regions. The Wasatch Front records show a significantly different near‐source rate of distance attenuation than the BA18 model, which we attribute to differences in (apparent) geometric attenuation. Near‐source residuals show a period dependence of this effect, with greater attenuation at shorter periods ( ⁠ T < 0.5    s "> T < 0.5 s ) and a correlation between period and the distance over which the discrepancy manifests ( ⁠ ∼ 20 – 50    km "> ∼ 20 – 50 km ⁠ ). We adjusted the recorded ground motions for these regional path effects and solved for station site terms using linear mixed‐effects regressions, with groupings for events and stations. We analyzed basin amplification by comparing the site terms with the basin geometry and basin depths from two seismic‐velocity models for the region. Sites over the deeper parts of the sedimentary basins are amplified by factors of 3–10, relative to sites with thin sedimentary cover, with greater amplification at longer periods ( ⁠ T ≳ 1    s "> T ≳ 1 s ). Average ground‐motion variability increases with period, and long‐period variability exhibits a slight increase at the basin edges. These results indicate regional seismic wave propagation effects requiring further study, and potentially a regionalized GMM, as well as highlight basin amplification complexities that may be incorporated into seismic hazard assessments.

Utah↗

USGS2016

A rare moderate‐sized (Mw 4.9) earthquake in Kansas: Rupture process of the Milan, Kansas, earthquake of 12 November 2014 and its relationship to fluid injection

The largest recorded earthquake in Kansas occurred northeast of Milan on 12 November 2014 ( M w 4.9) in a region previously devoid of significant seismic activity. Applying multistation processing to data from local stations, we are able to detail the rupture process and rupture geometry of the mainshock, identify the causative fault plane, and delineate the expansion and extent of the subsequent seismic activity. The earthquake followed rapid increases of fluid injection by multiple wastewater injection wells in the vicinity of the fault. The source parameters and behavior of the Milan earthquake and foreshock–aftershock sequence are similar to characteristics of other earthquakes induced by wastewater injection into permeable formations overlying crystalline basement. This earthquake also provides an opportunity to test the empirical relation that uses felt area to estimate moment magnitude for historical earthquakes for Kansas.

Seismological Research Letters↗

USGS1988

Inversion for slip distribution using teleseismic P waveforms: North Palm Springs, Borah Peak, and Michoacan earthquakes

We have inverted the teleseismic P waveforms recorded by stations of the Global Digital Seismograph Network for the 8 July 1986 North Palm Springs, California, the 28 October 1983 Borah Peak, Idaho, and the 19 September 1985 Michoacan, Mexico, earthquakes to recover the distribution of slip on each of the faults using a point-by-point inversion method with smoothing and positivity constraints. In the inversion procedure, a fault plane with fixed strike and dip is placed in the region of the earthquake hypocenter and divided into a large number of subfaults. Rupture is assumed to propagate at a constant velocity away from the hypocenter, and synthetic ground motions for pure strike-slip and dip-slip dislocations are calculated at the teleseismic stations for each subfault. The observed seismograms are then inverted to obtain the distribution of strike-slip and dip-slip displacement for the earthquake.

California↗

USGS1989

On the use of volumetric strain meters to infer additional characteristics of short-period seismic radiation

Volumetric strain meters (Sacks-Evertson design) are installed at 15 sites along the San Andreas fault system, to monitor long-term strain changes for earthquake prediction. Deployment of portable broadband, high-resolution digital recorders (GEOS) at several of the sites extends the detection band for volumetric strain to periods shorter than 5 × 10 −2 sec and permits the simultaneous observation of seismic radiation fields using conventional short-period pendulum seismometers. Simultaneous observations establish that the strain detection bandwidth extends from periods greater than 10 7 seconds to periods near 5 × 10 −2 sec with a dynamic range exceeding 140 dB. Measurements of earth-strain noise for the period band, 10 7 to 10 −2 sec, show that ground noise, not instrument noise, currently limits the measurement of strain over a bandwidth of more than eight orders of magnitude in period. Comparison of the short-period portion of earth-strain, noise spectra (20 to 5 × 10 −2 sec) with average spectra determined from pendulum seismometers, suggest that observed noise is predominantly dilatational energy. Recordings of local and regional earthquakes indicate that dilatometers respond to P energy but not direct shear energy and that straingrams can be used to resolve superimposed reflected P and S waves for inference of wave characteristics not permitted by either sensor alone. Simultaneous measurements of incident P - and S -wave amplitudes are used to introduce a technique for single-station estimates of wave field inhomogeneity, free-surface reflection coefficients and local material P velocity. Estimates of these parameters derived for the North Palm Springs earthquake ( M w 5.9) respectively for an incident P wave of 29° are −85°, 1.71, 2.9 km/sec, and for an incident S wave of 17° are 79°, 0.85, 2.9 km/sec. The empirical estimates of reflection coefficients are consistent with model estimates derived using an anelastic half-space model with incident inhomogeneous wave fields.

California↗

USGS2017

Apparent late Quaternary fault slip rate increase in the southwestern Lower Rhine Graben, central Europe

In regions of low strain, long earthquake recurrence intervals (10 4 –10 6 yrs) and erosive processes limit preservation of Quaternary markers suitable for distinguishing whether faults slip at uniform or secularly varying rates. The Lower Rhine graben in the border region of Germany, The Netherlands, and Belgium provides a unique opportunity to explore Quaternary slip‐rate variations in a region of low strain using the basal (2.29±0.29 Ma) and surface (700±80 ka) contacts of the regionally extensive main terrace (“Hauptterrasse”), deposited by the Rhine and Maas Rivers. These surfaces are vertically offset 3–140 m and 0–68 m, respectively, across individual fault strands within a distributed network of northwest‐trending, slow‐slipping (<0.1 mm/yr) normal faults. In this investigation, we construct Quaternary slip histories for the southern Lower Rhine graben faults using new main terrace surface vertical offset measurements made from light detection and ranging (lidar)‐derived bare‐earth digital terrain models, which we synthesize with existing constraints on the offset basal contact of this fluvial deposit ( n =91 collocated sites with displacement constraints). We find that >80% of the sites record an apparent increase in slip rate for the more recent interval from 700 ka to present, which corresponds to a period of increased uplift of the nearby Rhenish Massif and regional volcanism. However, the apparent increase in slip rate could result, in part, from erosion of the footwall surface below the main terrace, leading to an apparent displacement that is smaller than the total vertical offset since the start of the Quaternary. Prior work focused on characterization of these faults as seismic sources in the Lower Rhine graben has preferentially relied on the average fault‐slip rate constrained using the base of the main terrace. We suggest that average fault‐slip rates calculated using the ∼700 ka main terrace surface are subjected to fewer uncertainties and sample a time interval that is more relevant for seismic‐hazard analysis.

Lower Rhine graben↗

USGS2026

Landslide volume estimation from seismic waveform features using a global catalog of seismogenic events

Seismic waves generated by large, rapid landslides encode information about the source and can be analyzed rapidly following an event. Even remote landslides can pose hazards to downstream communities, so rapid detection and characterization using existing seismic monitoring networks could be beneficial. In this study, we expand on past regionally limited work by presenting a globally applicable method for estimating landslide volume from seismic features that could be integrated into future landslide seismic monitoring frameworks. We train the model using multivariable linear regression and five seismic features derived from recordings of 129 landslide events of a range of styles and locations with independently estimated volumes. We present two preferred models, one that combines long‐period (LP) and high‐frequency (HF) features and one for use on smaller landslides without observed long‐period signals. We find that our best‐performing model, applicable to landslides larger than 100,000 m 3 with signals containing observable long‐period energy, requires only two features: LP (20–100 s) absolute maximum amplitude and HF (1–5 Hz) rise time (time between the signal start and the maximum envelope amplitude) and has an R 2 score of 0.79. This model predicted volumes within one order of magnitude for 55 out of 58 events. We find that this combination mitigates the trade‐off between mass and acceleration and other variations of landslide style that limit methods based on amplitude alone.

Seismological Research Letters↗

USGS1997

Variability of site response in the Los Angeles urban area

This article addresses the variability of site response in the Los Angeles area and possible structural causes for the observations. Aftershock records from 231 sites in the San Fernando and Los Angeles basins and the surrounding mountains are used in this study. Spectral ratios, taken with respect to a low-amplitude reference site, are used to document the variation in site amplification in the frequency range 2 to 6 Hz, both spatially and with backazimuth to the source. At higher frequencies (6 to 10 Hz), spectral ratios are shown to have greater spatial variability. Interstation spectral ratios are used to measure the standard deviation among source as a function of station separation. An increase in the variation in ground motion is shown to take place at a station separation of 1 km. Relative site-response estimates between nearby stations are used to demonstrate that preferred directions of motion can exist even in areas with no surface topographic effects.

California↗

USGS2015

Ground motion observations of the 2014 South Napa earthquake

Ground motions of the South Napa earthquake (24 August 2014; M 6.0) were recorded at 19 stations within 20 km and 292 stations within 100 km of the rupture surface trace, generating peak ground motions in excess of 50% g and 50  cm/s in and near Napa Valley. This large dataset allows us to compare the ground motion from the earthquake to existing ground‐motion prediction equations (GMPEs) in considerable detail. Using the ground‐motion data compiled and reported by ShakeMap ( Wald et al. , 2000 ), we examine the peak ground acceleration (PGA) and peak ground velocity (PGV), as well as the pseudospectral acceleration (PSA) at periods of 0.3, 1.0, and 3.0 s. At the higher frequencies, especially PGA, data recorded at close distances (within &sim;20  km) are very consistent with the GMPEs, implying a stress drop for this event similar to the median for California, that is, 5 MPa ( Baltay and Hanks, 2014 ). At all frequencies, the attenuation with distance is stronger than the GMPEs would predict, which suggests the attenuation in the Napa and San Francisco Bay delta region is stronger than the average attenuation in California. The spatial plot of the ground‐motion residuals is positive to the north, in both Napa and Sonoma Valleys, consistent with increases in amplitude expected from both the directivity and basin effects. More interestingly, perhaps, there is strong ground motion to the south in the along‐strike direction, particularly for PSA at 1.0 s. These strongly positive residuals align with an older, Quaternary fault structure associated with the Franklin or Southampton fault, potentially indicating a fault‐zone‐guided wave.

California↗

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↗

USGS2016

Ground motion in the presence of complex Topography II: Earthquake sources and 3D simulations

Eight seismic stations were placed in a linear array with a topographic relief of 222 m over Mission Peak in the east San Francisco Bay region for a period of one year to study topographic effects. Seventy‐two well‐recorded local earthquakes are used to calculate spectral amplitude ratios relative to a reference site. A well‐defined fundamental resonance peak is observed with individual station amplitudes following the theoretically predicted progression of larger amplitudes in the upslope direction. Favored directions of vibration are also seen that are related to the trapping of shear waves within the primary ridge dimensions. Spectral peaks above the fundamental one are also related to topographic effects but follow a more complex pattern. Theoretical predictions using a 3D velocity model and accurate topography reproduce many of the general frequency and time‐domain features of the data. Shifts in spectral frequencies and amplitude differences, however, are related to deficiencies of the model and point out the importance of contributing factors, including the shear‐wave velocity under the topographic feature, near‐surface velocity gradients, and source parameters.

California↗

USGS2011

Reassessment of stable continental regions of Southeast Asia

Probabilistic seismic-hazard assessments of the central and eastern United States (CEUS) require estimates of the size of the largest possible earthquake (Mmax). In most of the CEUS, sparse historical seismicity does not provide a record of moderate and large earthquakes that is sufficient to constrain Mmax. One remedy for the insufficient catalog is to combine the catalog of moderate to large CEUS earthquakes with catalogs from other regions worldwide that are tectonically analogous to the CEUS (stable continental regions, or SCRs). After the North America SCR, the largest contribution of earthquakes to this global SCR catalog comes from a Southeast Asian SCR that extends from Indochina to southeasternmost Russia. Integration and interpretation of recently published geological and geophysical results show that most of these Southeast Asian earthquakes occurred in areas exposing abundant alkaline igneous rocks and extensional faults, both of Neogene age (last 23 million years). The implied Neogene extension precludes classification of the areas as SCR crust. The extension also reduces the number of moderate and large Southeast Asian historical earthquakes that are available to constrain CEUS Mmax by 86 percent, from 43 to six.

Seismological Research Letters↗

USGS2014

The earthquake cycle in the San Francisco Bay region: A.D. 1600–2012

Stress changes produced by the 1906 San Francisco earthquake had a profound effect on the seismicity of the San Francisco Bay region (SFBR), dramatically reducing it in the twentieth century. Whether the SFBR is still within or has emerged from this seismic quiescence is an issue of debate with implications for earthquake mechanics and seismic hazards. Historically, the SFBR has not experienced one complete earthquake cycle (i.e., the accumulation of stress, its release primarily as coseismic slip during surface‐faulting earthquakes, its re‐accumulation in the interval following, and its subsequent rerelease). The historical record of earthquake occurrence in the SFBR appears to be complete at about M 5.5 back to 1850 ( Bakun, 1999 ). For large events, the record may be complete back to 1776, which represents about half a cycle. Paleoseismic data provide a more complete view of the most recent pre‐1906 SFBR earthquake cycle, extending it back to about 1600. Using these, we have developed estimates of magnitude and seismic moment for alternative sequences of surface‐faulting paleoearthquakes occurring between 1600 and 1776 on the region’s major faults. From these we calculate seismic moment and moment release rates for different time intervals between 1600 and 2012. These show the variability in moment release and suggest that, in the SFBR regional plate boundary, stress can be released on a single fault in great earthquakes such as that in 1906 and in multiple ruptures distributed on the regional plate boundary fault system on a decadal time scale.

California↗

USGS2022

Quantifying the sensitivity of microearthquake slip inversions to station distribution using a dense nodal array

To investigate the sensitivity of slip inversions to station distribution and choice of empirical Green’s function (EGF), we examine three microearthquakes that occurred within the high‐density LArge‐n Seismic Survey in Oklahoma (LASSO) nodal seismic array. The LASSO array’s dense distribution of 1825 geophones provides an exceptional level of spatial and azimuthal coverage, allowing for more accurate inversions of slip than are possible with typical station distributions. The highly accurate slip inversions, in turn, allow for the exploration of the sensitivity of slip inversions to station distribution and parameter choices. We examine the effects of these choices using three well‐recorded strike‐slip microearthquakes ( M L 1.7, 2.3, and 2.7) using an EGF method. From this analysis and the systematic testing of varied network arrangements, we find that station distributions that have uniform coverage of azimuth and distance can retrieve the overall pattern of slip, but the estimated amplitude of slip can vary by 30% for high‐slip regions due to small variations in station location. In addition, we find that the distance range that accurately resolves the overall pattern of slip is the one that contains the takeoff angles of 45°–65°. Concerning azimuthal coverage, a network with >270° performs similarly to having complete coverage. The choice of EGF can shift the location of resolved areas of slip and their amplitude, depending on its similarity in location and radiation pattern to the target earthquake.

Bulletin of the Seismological Society of America↗

USGS2024

Site response in the Walnut Creek–Concord region of the San Francisco Bay, California: Ground motion amplification in a fault-bounded basin

Thirty‐seven portable accelerometers were deployed in the eastern San Francisco Bay communities of Walnut Creek and Concord to study site response in a fault‐bounded, urban, sedimentary basin. Local earthquakes were recorded for a period of two years from 2017 to 2019 resulting in 101 well‐recorded events. Site response is estimated by two methods: the reference site spectral ratio method and a source‐site spectral inversion method. The reference site spectral ratio method allows investigation of the variability of site amplification with source azimuth and frequency. The source‐site spectral inversion method yields the best least‐squares fit to site response for a database of ground‐motion records. Both methods show substantial amplification in the Walnut Creek–Concord basin below 2 Hz indicating strong surface‐wave development. Greater amplification is seen for sources aligned along the long axis of the basin. Inversion using close‐in sources at short distances yields lower amplification at longer periods than the entire data set due to reduced surface‐wave generation for steeper angles of incidence. Inversion of site response spectra for shallow shear‐wave velocity using a global search algorithm yields V S30 values consistent with generalized mapping results based on geology and topography but with greater variability due to local site variations. 3D finite‐element modeling shows greater amplification in the Walnut Creek–Concord basin with a basin‐edge effect likely contributing to higher ground motions. Topography is also seen to lead to increased scattering and shadowing effects.

California↗