USGS Science⌕ Search

SEARCH · USGS Science

Results for “Earthquake Spectra”

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 253 records · Page 14Linked to original sources
USGS2011

Regional spectral analysis of three moderate earthquakes in Northeastern North America

We analyze Fourier spectra obtained from the horizontal components of broadband and accelerogram data from the 1997 Cap-Rouge, the 2002 Ausable Forks, and the 2005 Rivière-du-Loup earthquakes, recorded by Canadian and American stations sited on rock at hypocentral distances from 23 to 602 km. We check the recorded spectra closely for anomalies that might result from site resonance or source effects. We use Beresnev and Atkinson’s (1997) near-surface velocity structures and Boore and Joyner’s (1997) quarter-wave method to estimate site response at hard- and soft-rock sites. We revise the Street et al. (1975) model for geometrical spreading, adopting a crossover distance of r o =50 km instead of 100 km. We obtain an average attenuation of Q =410±25f 0.50±0.03 for S+L g +surface waves with ray paths in the Appalachian and southeastern Grenville Provinces. We correct the recorded spectra for attenuation and site response to estimate source spectral shape and radiated energy for these three earthquakes and the 1988 M 5.8 Saguenay earthquake. The Brune stress drops range from 130 to 419 bars, and the apparent stresses range from 39 to 63 bars. The corrected source spectral shapes of these earthquakes are somewhat variable for frequencies from 0.2 to 2 Hz, falling slightly below the fitted Brune spectra.

Bulletin of the Seismological Society of America↗
USGS2012

Comparison of ground motion from tremors and explosions in deep gold mines

Seismic body waves, from tamped chemical explosions, two with yields of 50 and one of 150 kg, were compared with corresponding data from three mining-induced tremors with a view to testing methods of discriminating between the two types of events. Detonated at depths of about 2 km, all three explosions generated P waves for which the low-frequency spectral asymptotes agree well with corresponding results scaled down from nuclear shots at the Nevada Test Site. For the two smaller explosions, recorded underground at hypocentral distances ranging from 234 to 871 m, the P wave corner frequencies, defined by the intersections of the low- and high-frequency asymptotes, are also in fine agreement with scaled-down results from the Nevada Test Site; for the 150-kg explosion, which was recorded only at the surface, it appears that 2.39 km of upward propagation caused the corner frequency to be at least a factor of 4 lower than anticipated due to attenuation. All three explosions generated S waves that appear to be a consequence of deviatoric stress release in the immediate environs of the explosions. The three tremors analyzed here generated P and S waves whose spectra and source parameters agree well with standard earthquake source models. With regard to discrimination based on P wave spectra, we conclude that for events of fixed low-frequency spectral asymptotes, the explosions typically have higher corner frequencies than tremors or earthquakes, although counterexamples certainly exist. Interestingly, the 150-kg explosion was identified as such on the basis of P and S wave polarities that are incompatible with the normally expected double-couple source model; instead, these initial motions are consistent with an explosion in conjunction with normal faulting. The body wave spectra of this explosion and those of a nearby tremor, however, were indistinguishable.

Journal of Geophysical Research Solid Earth↗
USGS2019

Site response in the Oklahoma region from seismic recordings of the 2011 Mw 5.7 Prague earthquake

We invert the shear-wave displacement spectra obtained from 30 three-component, broadband waveforms recorded within 300 km of the 6 November 2011 Mw 5.7 Prague, Oklahoma earthquake to recover the site-response contribution using an inversion method that simultaneously inverts for source, path, and site effects. Site-response functions identify resonant frequencies within a range of 0.1-10 Hz that generally coincide with spectral peaks in H/V curves derived from the recorded waveforms. S-wave velocity profiles available for several sites were also used to calculate theoretical SH transfer functions that predict the site amplification due to the near-surface soil structure down to depths of 30-50 m. These transfer functions are generally flat below about 8 Hz in the frequency range sampled by the spectral inversion process, indicating that the spectral peaks in the site response obtained from the waveform analysis result from deeper velocity variations. A 0.3-Hz spectral peak observed at several stations, for example, coincides with the strong, surface-wave amplitudes observed at 3s periods for induced M ≥ 3 earthquakes in Oklahoma and Kansas, suggesting that this resonant peak may be due to surface waves trapped in the upper sedimentary layer of the crust. Both shallow and deep contributions to the site response are important for the characterization of ground motion from Central and Eastern North America (CENA) earthquakes. We obtain a seismic moment of 4.32 × 1024 dyne-cm and a corner frequency of 0.229, consistent with the magnitude of the event. A frequency-dependent attenuation relation of Q(f)=1107f 0.398 consistent with prior measurements of path properties in CENA is also derived.

Oklahoma↗
USGS2011

Moments, magnitudes, and radiated energies of non-volcanic tremor near Cholame, CA, from ground motion spectra at UPSAR

By averaging the spectra of events within two episodes of tremor (on Jan. 21 and 24, 2005) across the 12 stations of UPSAR, we improved the S/N sufficiently to define source spectra. Analysis of eleven impulsive events revealed attenuation-corrected spectra of displacement similar to those of earthquakes, with a low-frequency plateau, a corner frequency, and a high frequency decay proportional to f −2 . Seismic moments, M 0 , estimated from these spectra range from about 3 to 10 × 10 11 N-m or moment magnitudes in the range 1.6 to 1.9. The corner frequencies range from 2.6 to 7.2 Hz and, if interpreted in the same way as for earthquakes, indicate low stress drops that vary from 0.001 to 0.04 MPa. Seismic energies, estimated from the ground motion spectra, vary from 0.2 × 10 5 to 4.4 × 10 5 J, or apparent stresses in the range 0.002 to 0.02 MPa. The low stress parameters are consistent with a weak fault zone in the lower crust at the depth of tremor. In contrast, the same analysis on a micro-earthquake, located near Cholame (depth = 10.3 km), revealed a stress drop of 0.5 MPa and an apparent stress of 0.02 MPa. Residual spectra from ω −2 model fits to the displacement spectra of the non-volcanic tremor events show peaks near 4 Hz that are not apparent in the spectra for the microearthquake nor for the spectrum of earth noise. These spectral peaks may indicate that tremor entails more than shear failure reminiscent of mechanisms, possibly entailing fluid flow, associated with volcanic tremor or deep volcanic earthquakes.

Geophysical Research Letters↗
USGS1994

Estimates of site-dependent response spectra for design (methodology and justification)

Recent borehole-geotechnical data and strong-motion measurements constitute a new empirical basis to account for local geological conditions in earthquake-resistant design and site-dependent, building-code provisions. They provide new unambiguous definitions of site classes and rigorous empirical estimates of site-dependent amplification factors in terms of mean shear-wave velocity. A simple four-step methodology for estimating site-dependent response spectra is specified herein. Alternative techniques and commentary are presented for each step to facilitate application of the methodology for different purposes. Justification for the methodology is provided in terms of definitions for the new site classes and derivations of simple empirical equations for amplification as a function of mean shear-wave velocity and input ground-motion level. These new results provide a rigorous framework for improving estimates of site-dependent response spectra for design, site-dependent building-code provisions, and predictive maps of strong ground shaking for purposes of earthquake hazard mitigation.

California↗
USGS2003

Effect of causal and acausal filters on elastic and inelastic response spectra

With increasing interest in displacement spectra and long-period motions, it is important to check the sensitivity of both elastic and inelastic response spectra to the filtering that is often necessary to remove long period artifacts, even from many modern digital recordings. Using two records of very different character from the M=7.1, 1999 Hector Mine, California, earthquake, we find that the response spectra can be sensitive to the corner periods used in causal filtering, even for oscillator periods much less than the filter corner periods. The effect is most pronounced for inelastic response spectra, where the ratio of response spectra computed from accelerations filtered at 25 and 200 sec can be close to a factor of 2 for oscillator periods less than 5 sec. Published in 2003 by John Wiley and Sons, Ltd.

Earthquake Engineering and Structural Dynamics↗
USGS2003

Strong ground-motion prediction from Stochastic-dynamic source models

In the absence of sufficient data in the very near source, predictions of the intensity and variability of ground motions from future large earthquakes depend strongly on our ability to develop realistic models of the earthquake source. In this article we simulate near-fault strong ground motion using dynamic source models. We use a boundary integral method to simulate dynamic rupture of earthquakes by specifying dynamic source parameters (fracture energy and stress drop) as spatial random fields. We choose these quantities such that they are consistent with the statistical properties of slip heterogeneity found in finite-source models of past earthquakes. From these rupture models we compute theoretical strong-motion seismograms up to a frequency of 2 Hz for several realizations of a scenario strike-slip Mw 7.0 earthquake and compare empirical response spectra, spectra obtained from our dynamic models, and spectra determined from corresponding kinematic simulations. We find that spatial and temporal variations in slip, slip rise time, and rupture propagation consistent with dynamic rupture models exert a strong influence on near-source ground motion. Our results lead to a feasible approach to specify the variability in the rupture time distribution in kinematic models through a generalization of Andrews' (1976) result relating rupture speed to apparent fracture energy, stress drop, and crack length to 3D dynamic models. This suggests that a simplified representation of dynamic rupture may be obtained to approximate the effects of dynamic rupture without having to do full dynamic simulations.

Bulletin of the Seismological Society of America↗
USGS1994

The attenuation of seismic shear waves in quaternary alluvium in Santa Clara Valley, California

We used shear waves, generated by an air-powered source at the ground surface and recorded in a borehole, to estimate the shear-wave quality factor at strong-motion station Gilroy no. 2. We find similar values of Q using both the decay of the spectra with depth and the slope of the spectral ratio at two depths; we find no evidence of a frequency dependence of Q . The mean value of Q over the depth range 10 to 115 m is close to 10. The use of this value over the depth of the borehole and the observed travel time of 0.358 sec gives a cumulative attenuation factor t * of 0.036 sec for the upper 180 m of the Quaternary alluvium. This is comparable to the differential decay between Gilroy no. 2 and a rock site 1.9 km away (Gilroy no. 1), as measured from the decay of the high-frequency spectra of accelerograms from large earthquakes, plotted on a log-linear scale: t * = 0.05, 0.04, and 0.03 sec for the 1979 Coyote Lake, 1984 Morgan Hill, and 1989 Loma Prieta earthquakes, respectively. The similarity between the attenuations measured from the low-strain surface source and those from the larger amplitude earthquake sources suggests that increases of damping due to nonlinear wave propagation effects are limited.

California↗
USGS1993

Seismic response of eccentrically braced tall buildings

Spectral analysis and system identification techniques are used to analyze a set of acceleration reponse records obtained during the Loma Prieta earthquake from the 47-story, moment-resisting framed and eccentrically braced Embarcadero Building (EMB). The EMB was constructed in 1979 based on the 1976 Uniform Building Code requirements and a design response spectra defined by two levels of earthquake performances. The EMB is in the lower market area of San Francisco, which is of great interest to the engineering community because of the area's soft soil characteristics that amplify ground motions originating at long distances, and because the Embarcadero freeway suffered extensive damage during the earthquake and was razed in 1991. The first modal frequencies of the building at approximately 0.19 Hz (north-south) and 0.16 Hz (east-west) are identified. The torsional response and rocking motions of the building are insignificant. Discontinuity of stiffness and mass at the 40th floor level causes significant response issues above that floor such as excessive drift ratios.

Journal of Structural Engineering↗
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 &lt; 0.5 &#x2009;&#x2009; s "> T < 0.5 s ) and a correlation between period and the distance over which the discrepancy manifests ( ⁠ &#x223C; 20 &#x2013; 50 &#x2009;&#x2009; 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 &#x2273; 1 &#x2009;&#x2009; 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↗
USGS2005

An alternative approach to characterize nonlinear site effects

This paper examines the rationale of a method of nonstationary processing and analysis, referred to as the Hilbert-Huang transform (HHT), for its application to a recording-based approach in quantifying influences of soil nonlinearity in site response. In particular, this paper first summarizes symptoms of soil nonlinearity shown in earthquake recordings, reviews the Fourier-based approach to characterizing nonlinearity, and offers justifications for the HHT in addressing nonlinearity issues. This study then uses the HHT method to analyze synthetic data and recordings from the 1964 Niigata and 2001 Nisqually earthquakes. In doing so, the HHT-based site response is defined as the ratio of marginal Hilbert amplitude spectra, alternative to the Fourier-based response that is the ratio of Fourier amplitude spectra. With the Fourier-based approach in studies of site response as a reference, this study shows that the alternative HHT-based approach is effective in characterizing soil nonlinearity and nonlinear site response.

Earthquake Spectra↗
USGS2008

Ground-motion prediction equations for the average horizontal component of PGA, PGV, and 5%-damped PSA at spectral periods between 0.01 s and 10.0 s

This paper contains ground-motion prediction equations (GMPEs) for average horizontal-component ground motions as a function of earthquake magnitude, distance from source to site, local average shear-wave velocity, and fault type. Our equations are for peak ground acceleration (PGA), peak ground velocity (PGV), and 5%-damped pseudo-absolute-acceleration spectra (PSA) at periods between 0.01 s and 10 s. They were derived by empirical regression of an extensive strong-motion database compiled by the 'PEER NGA' (Pacific Earthquake Engineering Research Center's Next Generation Attenuation) project. For periods less than 1 s, the analysis used 1,574 records from 58 mainshocks in the distance range from 0 km to 400 km (the number of available data decreased as period increased). The primary predictor variables are moment magnitude (M), closest horizontal distance to the surface projection of the fault plane (RJB), and the time-averaged shear-wave velocity from the surface to 30 m (VS30). The equations are applicable for M=5-8, RJB<200 km, and VS30= 180-1300 m/s. ?? 2008, Earthquake Engineering Research Institute.

Earthquake Spectra↗
USGS2014

NGA-West2 Research Project

The NGA-West2 project is a large multidisciplinary, multi-year research program on the Next Generation Attenuation (NGA) models for shallow crustal earthquakes in active tectonic regions. The research project has been coordinated by the Pacific Earthquake Engineering Research Center (PEER), with extensive technical interactions among many individuals and organizations. NGA-West2 addresses several key issues in ground-motion seismic hazard, including updating the NGA database for a magnitude range of 3.0&ndash;7.9; updating NGA ground-motion prediction equations (GMPEs) for the &ldquo;average&rdquo; horizontal component; scaling response spectra for damping values other than 5%; quantifying the effects of directivity and directionality for horizontal ground motion; resolving discrepancies between the NGA and the National Earthquake Hazards Reduction Program (NEHRP) site amplification factors; analysis of epistemic uncertainty for NGA GMPEs; and developing GMPEs for vertical ground motion. This paper presents an overview of the NGA-West2 research program and its subprojects.

Earthquake Spectra↗
USGS2020

Peak ground motions and site response at Anza and Imperial Valley, California

Power spectra of shear-waves for eighteen earthquakes from the Anza-Imperial Valley region were inverted for source, mid-path Q, site attenuation and site response. The motivation was whether differences in site attenuation (parameterized as t*, r/cQ, where r is distance along ray path near the site, c is shear velocity and Q is the quality factor that parameterizes attenuation) and site response could be correlated with residuals in peak values of velocity or acceleration after removing the affect of distance-dependent attenuation. We decomposed spectra of S-waves from horizontal components of 18 earthquakes from 2010 to 2018 into a common source for each event with ω −2 spectral fall-off at high frequencies and then projected the residuals onto path and site terms following the methodology of Boatwright et al. (Bull Seismol Soc Am 81:1754–1782, 1991). The site terms were constrained to have an amplification at a particular frequency governed by V S30 at two of the sites which had downhole shear-wave logs. The 18 events, 3 < M < 4, had moments between approximately 10 20 and 10 22 dyne-cm, and stress drops between 1 and 100 bars. Average mid-crust attenuation had a Q of 844 reflecting the average path through the crystalline rock of the San Jacinto Mountains. t* for each station corresponded to the geologic environment such that stations on hard rock had low t* (e.g. stations KNW, PFO and RDM) a station in the San Jacinto fault zone (station SND) had a moderate t* of 0.035 s and stations in the Imperial Valley usually had higher t*s . Generally t* correlated with average amplification suggesting that sites characterized by low surface velocities and higher attenuation also have more amplification in the 1–6 Hz band. Residuals of peak values were determined by subtracting the prediction of Boore and Atkinson ( 2008 ). There is a correlation between average amplification and peak velocity, but not peak acceleration. Interestingly, there is less scatter at high values of amplification although there is also less data. Scatter in values of peak velocity and peak acceleration are higher at shorter compared to longer durations. When using a frequency-dependent form for Q , variances are higher, sometimes much higher; the dataset does not support frequency-dependent Q , which is not similar to results from the Imperial Valley and northeastern North America.

California↗
USGS2005

Earthquake and ambient vibration monitoring of the steel-frame UCLA factor building

Dynamic property measurements of the moment-resisting steel-frame University of California, Los Angeles, Factor building are being made to assess how forces are distributed over the building. Fourier amplitude spectra have been calculated from several intervals of ambient vibrations, a 24-hour period of strong winds, and from the 28 March 2003 Encino, California (ML = 2.9), the 3 September 2002 Yorba Linda, California (ML = 4.7), and the 3 November 2002 Central Alaska (Mw = 7.9) earthquakes. Measurements made from the ambient vibration records show that the first-mode frequency of horizontal vibration is between 0.55 and 0.6 Hz. The second horizontal mode has a frequency between 1.6 and 1.9 Hz. In contrast, the first-mode frequencies measured from earthquake data are about 0.05 to 0.1 Hz lower than those corresponding to ambient vibration recordings indicating softening of the soil-structure system as amplitudes become larger. The frequencies revert to pre-earthquake levels within five minutes of the Yorba Linda earthquake. Shaking due to strong winds that occurred during the Encino earthquake dominates the frequency decrease, which correlates in time with the duration of the strong winds. The first shear wave recorded from the Encino and Yorba Linda earthquakes takes about 0.4 sec to travel up the 17-story building. ?? 2005, Earthquake Engineering Research Institute.

Earthquake Spectra↗
USGS2026

An examination of the reasons for, and variability of, the differences between response spectra and pseudo response spectra

For many decades, it has been common practice to overlook the differences between response spectra and pseudo response spectra. This practice stems from analyses indicating that the average difference between response spectra and pseudo response spectra is negligible for a viscous damping ratio of 5%, which is near-universally adopted. As others have shown, the differences depend strongly on the damping ratio, earthquake magnitude, and oscillator period, with larger differences for larger damping, smaller magnitudes, and longer periods. To quantify the full statistical distribution of the ratio of the absolute acceleration response spectrum (SA) to the pseudo absolute acceleration response spectrum (PSA), we fit empirical ratios with the generalized gamma distribution. Data for magnitudes from 4.0 to 5.5 give 95% prediction intervals for the ratio of 1.00–1.02 at 1 s and 1.01–2.24 at 10 s for 5% damping. For 30% damping, these intervals increase to 1.14–2.01 at 1 s and 1.62–11.8 at 10 s, respectively. Further, we demonstrate how the shape of the ground-motion Fourier amplitude spectrum interacts with the natural period of the oscillator to control the differences in the response spectra type. Routinely treating the two measures of acceleration response spectra as interchangeable may lead to unacceptable inconsistencies when the response spectra of individual records are important or when larger damping ratios are used, such as in direct displacement-based design. If one must convert between them, the uncertainty of the conversion should be accounted for during analysis.

Earthquake Spectra↗