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.

SourceResearch connectionFrom yearThrough year

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

USGS1974

Time-dependent spectral analysis of thirty-one strong-motion earthquake records

Velocity response spectral analysis as a function of time for 5 percent critical damping is presented for the horizontal components of thirty-one earthquake records (61 components). A nomograph that allows conversion of the velocity response envelope to a displacement response envelope and a pseudo-absolute acceleration response envelope spectrum is also presented. Included among the records analyzed are such important accelerograms as those from El Centro, 1940, Olympia, 1949, Taft, 1952, San Francisco, 1957, Park-field, 1966, San Fernando, 1971 and Managua, 1972. Time-duration spectra showing the number of cycles for a particular period of the response envelope above different levels are presented. The effects of peak ground acceleration on the response envelope are examined.

Open-File Report↗

USGS2015

Ground motion simulation for the 23 August 2011, Mineral, Virginia earthquake using physics-based and stochastic broadband methods

Three broadband simulation methods are used to generate synthetic ground motions for the 2011 Mineral, Virginia, earthquake and compare with observed motions. The methods include a physics‐based model by Hartzell et al. (1999, 2005), a stochastic source‐based model by Boore (2009), and a stochastic site‐based model by Rezaeian and Der Kiureghian (2010, 2012). The ground‐motion dataset consists of 40 stations within 600 km of the epicenter. Several metrics are used to validate the simulations: (1) overall bias of response spectra and Fourier spectra (from 0.1 to 10 Hz); (2) spatial distribution of residuals for GMRotI50 peak ground acceleration (PGA), peak ground velocity, and pseudospectral acceleration (PSA) at various periods; (3) comparison with ground‐motion prediction equations (GMPEs) for the eastern United States. Our results show that (1) the physics‐based model provides satisfactory overall bias from 0.1 to 10 Hz and produces more realistic synthetic waveforms; (2) the stochastic site‐based model also yields more realistic synthetic waveforms and performs superiorly for frequencies greater than about 1 Hz; (3) the stochastic source‐based model has larger bias at lower frequencies (<0.5  Hz) and cannot reproduce the varying frequency content in the time domain. The spatial distribution of GMRotI50 residuals shows that there is no obvious pattern with distance in the simulation bias, but there is some azimuthal variability. The comparison between synthetics and GMPEs shows similar fall‐off with distance for all three models, comparable PGA and PSA amplitudes for the physics‐based and stochastic site‐based models, and systematic lower amplitudes for the stochastic source‐based model at lower frequencies (<0.5  Hz).

Bulletin of the Seismological Society of America↗

USGS2021

Regional calibration of hybrid ground‐motion simulations in moderate seismicity areas: Application to the Upper Rhine Graben

This study presents the coupling of the spectral decomposition results for anelastic attenuation, stress drop, and site effects with the Graves‐Pitarka (GP) hybrid ground‐motion simulation methodology, as implemented on the Southern California Earthquake Center (SCEC) broadband platform (BBP). It is targeted to applications in the Upper Rhine graben (URG), which is among the seismically active areas in western Europe, yet a moderate seismicity area. Our development consists of three main steps: (1) calibration of regional high‐frequency (HF) attenuation properties; (2) modification of the hybrid approach to add compressional waves in the HF computation and examine various strategies to evaluate site amplification factors in the Fourier domain (e.g., V S 30 "> V S 30 VS30 ‐based or site‐specific factors); (3) testing of the simulations using earthquake records from the URG ( ⁠ 3.7 < M w < 5 "> 3.7 < M w < 5 3.7<Mw<5 ⁠ ). The validation process of the simulated time histories is performed first on rock sites, and, then subsequently at all stations, whatever their site conditions. The performance of the simulations for rock sites is assessed through the standard validation technique in the BBP (comparison of the waveforms, intensity measures, and estimation of the response spectra model bias). We additionally compare the Fourier amplitude spectrum of the simulations and observations, and compute their corresponding bias. The results show that the simulated ground motions match the general characteristics of the recorded motions, and that the model bias generally fluctuates around zero across the broadband frequency range. Hence, the hybrid ground‐motion methodology implemented in the SCEC BBP can be successfully applied outside high‐seismicity areas and outside those areas for which it had been generally calibrated. Our results also show that HF modification and calibration were necessary to improve the fits with the observation, and demonstrate the potential benefits of using site‐specific amplification factors compared to V S 30 "> V S 30 VS30 ‐based amplification factors.

Upper Rhine Graben↗

USGS2014

A generalization of the double-corner-frequency source spectral model and its use in the SCEC BBP validation exercise

The stochastic method of simulating ground motions requires the specification of the shape and scaling with magnitude of the source spectrum. The spectral models commonly used are either single-corner-frequency or double-corner-frequency models, but the latter have no flexibility to vary the high-frequency spectral levels for a specified seismic moment. Two generalized double-corner-frequency ω 2 source spectral models are introduced, one in which two spectra are multiplied together, and another where they are added. Both models have a low-frequency dependence controlled by the seismic moment, and a high-frequency spectral level controlled by the seismic moment and a stress parameter. A wide range of spectral shapes can be obtained from these generalized spectral models, which makes them suitable for inversions of data to obtain spectral models that can be used in ground-motion simulations in situations where adequate data are not available for purely empirical determinations of ground motions, as in stable continental regions. As an example of the use of the generalized source spectral models, data from up to 40 stations from seven events, plus response spectra at two distances and two magnitudes from recent ground-motion prediction equations, were inverted to obtain the parameters controlling the spectral shapes, as well as a finite-fault factor that is used in point-source, stochastic-method simulations of ground motion. The fits to the data are comparable to or even better than those from finite-fault simulations, even for sites close to large earthquakes.

Bulletin of the Seismological Society of America↗

USGS2016

Seismic site characterization of an urban dedimentary basin, Livermore Valley, California: Site tesponse, basin-edge-induced surface waves, and 3D simulations

Thirty‐two accelerometers were deployed in the Livermore Valley, California, for approximately one year to study sedimentary basin effects. Many local and near‐regional earthquakes were recorded, including the 24 August 2014 M w 6.0 Napa, California, earthquake. The resulting ground‐motion data set is used to quantify the seismic response of the Livermore basin, a major structural depression in the California Coast Range Province bounded by active faults. Site response is calculated by two methods: the reference‐site spectral ratio method and a source‐site spectral inversion method. Longer‐period (≥1  s) amplification factors follow the same general pattern as Bouguer gravity anomaly contours. Site response spectra are inverted for shallow shear‐wave velocity profiles, which are consistent with independent information. Frequency–wavenumber analysis is used to analyze plane‐wave propagation across the Livermore Valley and to identify basin‐edge‐induced surface waves with back azimuths different from the source back azimuth. Finite‐element simulations in a 3D velocity model of the region illustrate the generation of basin‐edge‐induced surface waves and point out strips of elevated ground velocities along the margins of the basin.

Bulletin of the Seismological Society of America↗

USGS2016

Adjusting central and eastern North America ground-motion intensity measures between sites with different reference-rock site conditions

Adjustment factors are provided for converting ground‐motion intensity measures between central and eastern North America (CENA) sites with different reference‐rock site conditions ( V S 30 =760, 2000, and 3000 m/s) for moment magnitudes ranging from 2 to 8, rupture distances ranging from 2 to 1200 km, Fourier amplitude spectra (FAS) for frequencies ranging from 0.01 to 100 Hz, response spectra for periods ranging from 0.01 to 10.0 s, peak ground acceleration, and peak ground velocity. The adjustment factors are given for a wide range of the site diminution parameters ( κ 0 ) for sites with V S 30 =760 m/s and for a κ 0 of 0.006 s for two harder rock sites. Fourteen CENA velocity profiles with V S 30 values within a factor of 1.1 of 760 m/s were used to derive average FAS amplification factors as a function of frequency, which were then used in simulations of peak ground‐motion parameters and response spectra to derive the adjustment factors. The amplification function differs from that used in western North America (e.g., Campbell and Boore, 2016 ) in having a peak near 9 Hz, due to the resonance of motions in the relatively thin low‐velocity material over hard rock that characterizes many CENA sites with V S 30 near 760 m/s. We call these B/C sites, because this velocity marks the boundary between National Earthquake Hazards Reduction Program site classes B and C ( Building Seismic Safety Council, 2004 ). The adjustments for short‐period motions are sensitive to the value of κ 0 , but there are very few if any determinations of κ 0 for CENA B/C sites. For this reason, we determined κ 0 from multiple recordings at Pinyon Flat Observatory (PFO), California, which has a velocity‐depth profile similar to those of CENA B/C sites. The PFO and other results from the literature suggest that appropriate values of κ 0 for CENA B/C sites are expected to lie between 0.01 and 0.03 s.

Bulletin of the Seismological Society of America↗

USGS2006

A stochastic estimate of ground motion at Oceano, California, for the M 6.5 22 December 2003 San Simeon earthquake, derived from aftershock recordings

The U.S. Geological Survey deployed a digital seismic station in Oceano, California, in February 2004, to investigate the cause of damage and liquefaction from the 22 December 2003 M 6.5 San Simeon earthquake. This station recorded 11 M > 2.8 aftershocks in almost 8 weeks. We analyze these recordings, together with recordings of the mainshock and the same aftershocks obtained from nearby stations in Park Hill and San Luis Obispo, to estimate the mainshock ground motion in Oceano. We estimate the Fourier amplitude spectrum using generalized spectral ratio analysis. We test a set of aftershocks as Green's functions by comparing simulated and recorded acceleration amplitude spectra for the mainshock at San Luis Obispo and Park Hill. We convolve the aftershock accelerograms with a stochastic operator to simulate the duration and phase of the mainshock accelerograms. This approximation allows us to extend the range of aftershocks that can be used as Green's functions to events nearly three magnitude units smaller than the mainshock. Our realizations for the mainshock accelerogram at Oceano yield peak ground accelerations distributed as 28% ?? 4%g. We interpret these realizations as upper bounds for the actual ground motion, because our analysis assumes a linear response, whereas the presence of liquefaction indicates that the ground behaved nonlinearly in Oceano.

Bulletin of the Seismological Society of America↗

USGS2002

Regional estimates of radiated seismic energy

We revise the spectral technique for estimating radiated energy from recordings of large earthquakes at regional distances (Δ < 200 km) by correcting for geometric spreading and for site amplification as explicit functions of frequency. We analyze 65 recordings of the 1999 Hector Mine earthquake as functions of frequency, distance, and azimuth. For r > 27.5 km from the source, we model the geometrical spreading of the regional wavefield as r – γ where γ = 0.5 for f ≤ 0.2 Hz and γ = 0.7 for f ≥ 0.25 Hz. We fit the spectral falloff with distance using a frequency-dependent attenuation Q = 400( f /1.5) 0.6 , where Q = 400 for f ≤ 1.5 Hz. There is little directivity apparent in the corrected velocity spectra: the velocity spectra observed to the northwest along strike are amplified by a factor of 2.5 from 0.3 to 1.0 Hz and those to the southeast are amplified by a factor of 1.6 from 0.3 to 0.7 Hz. We group the stations in NEHRP site classes, using average 1-D velocity structures to estimate site amplification as a function of frequency and assuming 0.40 ≤ κ ≤ 0.55 sec for the near-surface attenuation. We increase the amplification of the soft-soil sites from 0.1 to 1.0 Hz by a factor that reaches 1.7 at 0.3 Hz because they are more strongly amplified than the NEHRP-D velocity structure predicts. We combine the 65 single-station estimates of radiated energy using an equal-azimuth weighting scheme that compensates for station distribution and incorporates the observed directivity, yielding a regional estimate of E s = 3.4 ± 0.7 × 10 22 dyne cm. This regional estimate of radiated energy corresponds closely to the teleseismic estimate of E s = 3.2 × 10 22 dyne cm.

Bulletin of the Seismological Society of America↗

USGS1990

Stochastic summation of empirical Green's functions

Two simple strategies are presented that use random delay times for repeatedly summing the record of a relatively small earthquake to simulate the effects of a larger earthquake. The simulations do not assume any fault plane geometry or rupture dynamics, but realy only on the ω−2 spectral model of an earthquake source and elementary notions of source complexity. The strategies simulate ground motions for all frequencies within the bandwidth of the record of the event used as a summand. The first strategy, which introduces the basic ideas, is a single-stage procedure that consists of simply adding many small events with random time delays. The probability distribution for delays has the property that its amplitude spectrum is determined by the ratio of ω−2 spectra, and its phase spectrum is identically zero. A simple expression is given for the computation of this zero-phase scaling distribution. The moment rate function resulting from the single-stage simulation is quite simple and hence is probably not realistic for high-frequency (>1 Hz) ground motion of events larger than ML ∼ 4.5 to 5. The second strategy is a two-stage summation that simulates source complexity with a few random subevent delays determined using the zero-phase scaling distribution, and then clusters energy around these delays to get an ω−2 spectrum for the sum. Thus, the two-stage strategy allows simulations of complex events of any size for which the ω−2 spectral model applies. Interestingly, a single-stage simulation with too few ω−2 records to get a good fit to an ω−2 large-event target spectrum yields a record whose spectral asymptotes are consistent with the ω−2 model, but that includes a region in its spectrum between the corner frequencies of the larger and smaller events reasonably approximated by a power law trend. This spectral feature has also been discussed as reflecting the process of partial stress release (Brune, 1970), an asperity failure (Boatwright, 1984), or the breakdown of ω−2 scaling due to rupture significantly longer than the width of the seismogenic zone (Joyner, 1984).

Bulletin of the Seismological Society of America↗

USGS2010

Coherence of Mach fronts during heterogeneous supershear earthquake rupture propagation: Simulations and comparison with observations

We study how heterogeneous rupture propagation affects the coherence of shear and Rayleigh Mach wavefronts radiated by supershear earthquakes. We address this question using numerical simulations of ruptures on a planar, vertical strike-slip fault embedded in a three-dimensional, homogeneous, linear elastic half-space. Ruptures propagate spontaneously in accordance with a linear slip-weakening friction law through both homogeneous and heterogeneous initial shear stress fields. In the 3-D homogeneous case, rupture fronts are curved owing to interactions with the free surface and the finite fault width; however, this curvature does not greatly diminish the coherence of Mach fronts relative to cases in which the rupture front is constrained to be straight, as studied by Dunham and Bhat (2008a). Introducing heterogeneity in the initial shear stress distribution causes ruptures to propagate at speeds that locally fluctuate above and below the shear wave speed. Calculations of the Fourier amplitude spectra (FAS) of ground velocity time histories corroborate the kinematic results of Bizzarri and Spudich (2008a): (1) The ground motion of a supershear rupture is richer in high frequency with respect to a subshear one. (2) When a Mach pulse is present, its high frequency content overwhelms that arising from stress heterogeneity. Present numerical experiments indicate that a Mach pulse causes approximately an ω −1.7 high frequency falloff in the FAS of ground displacement. Moreover, within the context of the employed representation of heterogeneities and over the range of parameter space that is accessible with current computational resources, our simulations suggest that while heterogeneities reduce peak ground velocity and diminish the coherence of the Mach fronts, ground motion at stations experiencing Mach pulses should be richer in high frequencies compared to stations without Mach pulses. In contrast to the foregoing theoretical results, we find no average elevation of 5%-damped absolute response spectral accelerations (SA) in the period band 0.05–0.4 s observed at stations that presumably experienced Mach pulses during the 1979 Imperial Valley, 1999 Kocaeli, and 2002 Denali Fault earthquakes compared to SA observed at non-Mach pulse stations in the same earthquakes. A 20% amplification of short period SA is seen only at a few of the Imperial Valley stations closest to the fault. This lack of elevated SA suggests that either Mach pulses in real earthquakes are even more incoherent that in our simulations or that Mach pulses are vulnerable to attenuation through nonlinear soil response. In any case, this result might imply that current engineering models of high frequency earthquake ground motions do not need to be modified by more than 20% close to the fault to account for Mach pulses, provided that the existing data are adequately representative of ground motions from supershear earthquakes.

Journal of Geophysical Research B: Solid Earth↗

USGS2002

Instrumental intensity distribution for the Hector Mine, California, and the Chi-Chi, Taiwan, earthquakes: Comparison of two methods

We compare two methods of seismic-intensity estimation from ground-motion records for the two recent strong earthquakes: the 1999 (M 7.1) Hector Mine, California, and the 1999 (M 7.6) Chi-Chi, Taiwan. The first technique utilizes the peak ground acceleration (PGA) and velocity (PGV), and it is used for rapid generation of the instrumental intensity map in California. The other method is based on the revised relationships between intensity and Fourier amplitude spectrum (FAS). The results of using the methods are compared with independently observed data and between the estimations from the records. For the case of the Hector Mine earthquake, the calculated intensities in general agree with the observed values. For the case of the Chi-Chi earthquake, the areas of maximum calculated intensity correspond to the areas of the greatest damage and highest number of fatalities. However, the FAS method producees higher-intensity values than those of the peak amplitude method. The specific features of ground-motion excitation during the large, shallow, thrust earthquake may be considered a reason for the discrepancy. The use of PGA and PGV is simple; however, the use of FAS provides a natural consideration of site amplification by means of generalized or site-specific spectral ratios. Because the calculation of seismic-intensity maps requires rapid processing of data from a large network, it is very practical to generate a "first-order" map from the recorded peak motions. Then, a "second-order" map may be compiled using an amplitude-spectra method on the basis of available records and numerical modeling of the site-dependent spectra for the regions of sparse station spacing.

California↗

USGS2009

On baseline corrections and uncertainty in response spectrafor baseline variations commonly encountered in digital accelerograph records

Most digital accelerograph recordings are plagued by long-period drifts, best seen in the velocity and displacement time series obtained from integration of the acceleration time series. These drifts often result in velocity values that are nonzero near the end of the record. This is clearly unphysical and can lead to inaccurate estimates of peak ground displacement and long-period spectral response. The source of the long-period noise seems to be variations in the acceleration baseline in many cases. These variations could be due to true ground motion (tilting and rotation, as well as local permanent ground deformation), instrumental effects, or analog-to-digital conversion. Very often the trends in velocity are well approximated by a linear trend after the strong shaking subsides. The linearity of the trend in velocity implies that no variations in the baseline could have occurred after the onset of linearity in the velocity time series. This observation, combined with the lack of any trends in the pre-event motion, allows us to compute the time interval in which any baseline variations could occur. We then use several models of the variations in a Monte Carlo procedure to derive a suite of baseline-corrected accelerations for each noise model using records from the 1999 Chi-Chi earthquake and several earthquakes in Turkey. Comparisons of the mean values of the peak ground displacements, spectral displacements, and residual displacements computed from these corrected accelerations for the different noise models can be used as a guide to the accuracy of the baseline corrections. For many of the records considered here the mean values are similar for each noise model, giving confidence in the estimation of the mean values. The dispersion of the ground-motion measures increases with period and is noise-model dependent. The dispersion of inelastic spectra is greater than the elastic spectra at short periods but approaches that of the elastic spectra at longer periods. The elastic spectra from the most basic processing, in which only the pre-event mean is removed from the acceleration time series, do not diverge from the baseline-corrected spectra until periods of 10-20 sec or more for the records studied here, implying that for many engineering purposes elastic spectra can be used from records with no baseline correction or filtering.

Bulletin of the Seismological Society of America↗

USGS2007

Multi-component ground motion response spectra for coupled horizontal, vertical, angular accelerations, and tilt

Rotational and vertical components of ground motion are almost always ignored in design or in the assessment of structures despite the fact that vertical motion can be twice as much as the horizontal motion and may exceed 2g level, and rotational excitation may reach few degrees in the proximity of fault rupture. Coupling of different components of ground excitation may significantly amplify the seismic demand by introducing additional lateral forces and enhanced P-?? effects. In this paper, a governing equation of motion is postulated to compute the response of a SDOF oscillator under a multi-component excitation. The expanded equation includes secondary P-?? components associated with the combined impacts of tilt and vertical excitations in addition to the inertial forcing terms due to the angular and translational accelerations. The elastic and inelastic spectral ordinates traditionally generated considering the uniaxial input motion are compared at the end with the multi-component response spectra of coupled horizontal, vertical and tilting motions. The proposed multi-component response spectrum reflects kinematic characteristics of the ground motion that are not identifiable by the conventional spectrum itself, at least for the near-fault region where high intensity vertical shaking and rotational excitation are likely to occur.

ISET Journal of Earthquake Technology↗

USGS1996

Directional topographic site response at Tarzana observed in aftershocks of the 1994 Northridge, California, earthquake: Implications for mainshock motions

The Northridge earthquake caused 1.78 g acceleration in the east-west direction at a site in Tarzana, California, located about 6 km south of the mainshock epicenter. The accelerograph was located atop a hill about 15-m high, 500-m long, and 130-m wide, striking about N78°E. During the aftershock sequence, a temporary array of 21 three-component geophones was deployed in six radial lines centered on the accelerograph, with an average sensor spacing of 35 m. Station C00 was located about 2 m from the accelerograph. We inverted aftershock spectra to obtain average relative site response at each station as a function of direction of ground motion. We identified a 3.2-Hz resonance that is a transverse oscillation of the hill (a directional topographic effect). The top/base amplification ratio at 3.2 Hz is about 4.5 for horizontal ground motions oriented approximately perpendicular to the long axis of the hill and about 2 for motions parallel to the hill. This resonance is seen most strongly within 50 m of C00. Other resonant frequencies were also observed. A strong lateral variation in attenuation, probably associated with a fault, caused substantially lower motion at frequencies above 6 Hz at the east end of the hill. There may be some additional scattered waves associated with the fault zone and seen at both the base and top of the hill, causing particle motions (not spectral ratios) at the top of the hill to be rotated about 20° away from the direction transverse to the hill. The resonant frequency, but not the amplitude, of our observed topographic resonance agrees well with theory, even for such a low hill. Comparisons of our observations with theoretical results indicate that the 3D shape of the hill and its internal structure are important factors affecting its response. The strong transverse resonance of the hill does not account for the large east-west mainshock motions. Assuming linear soil response, mainshock east-west motions at the Tarzana accelerograph were amplified by a factor of about 2 or less compared with sites at the base of the hill. Probable variations in surficial shear-wave velocity do not account for the observed differences among mainshock acceleration observed at Tarzana and at two different sites within 2 km of Tarzana.

California↗

USGS2024

Framework for implementing damping scaling factors in U.S. Geological Survey National Seismic Hazard Models

Traditionally, probabilistic seismic hazard analysis (PSHA) has focused on calculating ground motion hazard curves for elastic, 5%-damped pseudo spectral accelerations, Sa(T,5%), which are used as the basis for engineering design parameters and targets for ground motion selection and modification. However, structures and geotechnical systems can exhibit a wide range of damping ratios both above and below the 5% level, depending on the construction material, structural system, nonstructural elements, or subsurface soil properties. When spectral parameters at such damping levels are required for certain applications, 5%-damped accelerations have traditionally been extracted from PSHA-based hazard curves and adjusted outside of the hazard integral using damping scaling factors (DSF) such as those from Newmark & Hall (1982). Recent advances in the development of more rigorous and comprehensive damping scaling models (e.g., Rezaeian et al., 2014; Rezaeian et al., 2021) have allowed for the modeling of means and standard deviations of DSFs as functions of earthquake source and path properties for crustal, intraslab, and subduction interface tectonic environments. These DSF models can be applied to ground motion model (GMM) estimates of Sa(T,5%) for a given earthquake rupture scenario to produce a corresponding mean and standard deviation Sa at a specified damping ratio β, Sa(T,β). In this study, the DSF models of Rezaeian et al. (2014) and Rezaeian et al. (2021) are implemented within the U.S. Geological Survey National Seismic Hazard Model (NSHM) PSHA framework to calculate probabilistic hazard curves for spectral accelerations at damping ratios from 0.5% to 30%. The DSF models are applied directly to the mean and standard deviation of Sa(T,5%) predictions from each GMM in the NSHM logic tree. Resulting hazard curves and uniform hazard and risk spectra for Sa(T,β) are presented for several geographic locations and compared with corresponding spectra estimated using current design practices by applying the same DSFs outside of the PSHA calculation. Key differences between the two methods for estimating Sa(T,β) are discussed, and potential strategies are presented for the implementation and usage of the hazard-consistent Sa(T,β) in building codes. Comparing the results to those from DSFs used in current design practices that are mainly based on Newmark & Hall (1982) is not explored in this study.

Conference Paper↗

USGS2013

The effect of sampling rate and anti-aliasing filters on high-frequency response spectra

The most commonly used intensity measure in ground-motion prediction equations is the pseudo-absolute response spectral acceleration (PSA), for response periods from 0.01 to 10 s (or frequencies from 0.1 to 100 Hz). PSAs are often derived from recorded ground motions, and these motions are usually filtered to remove high and low frequencies before the PSAs are computed. In this article we are only concerned with the removal of high frequencies. In modern digital recordings, this filtering corresponds at least to an anti-aliasing filter applied before conversion to digital values. Additional high-cut filtering is sometimes applied both to digital and to analog records to reduce high-frequency noise. Potential errors on the short-period (high-frequency) response spectral values are expected if the true ground motion has significant energy at frequencies above that of the anti-aliasing filter. This is especially important for areas where the instrumental sample rate and the associated anti-aliasing filter corner frequency (above which significant energy in the time series is removed) are low relative to the frequencies contained in the true ground motions. A ground-motion simulation study was conducted to investigate these effects and to develop guidance for defining the usable bandwidth for high-frequency PSA. The primary conclusion is that if the ratio of the maximum Fourier acceleration spectrum (FAS) to the FAS at a frequency f saa corresponding to the start of the anti-aliasing filter is more than about 10, then PSA for frequencies above f saa should be little affected by the recording process, because the ground-motion frequencies that control the response spectra will be less than f saa . A second topic of this article concerns the resampling of the digital acceleration time series to a higher sample rate often used in the computation of short-period PSA. We confirm previous findings that sinc-function interpolation is preferred to the standard practice of using linear time interpolation for the resamplin

Bulletin of Earthquake Engineering↗

USGS2011

Model and parametric uncertainty in source-based kinematic models of earthquake ground motion

Four independent ground-motion simulation codes are used to model the strong ground motion for three earthquakes: 1994 M w 6.7 Northridge, 1989 M w 6.9 Loma Prieta, and 1999 M w 7.5 Izmit. These 12 sets of synthetics are used to make estimates of the variability in ground-motion predictions. In addition, ground-motion predictions over a grid of sites are used to estimate parametric uncertainty for changes in rupture velocity. We find that the combined model uncertainty and random variability of the simulations is in the same range as the variability of regional empirical ground-motion data sets. The majority of the standard deviations lie between 0.5 and 0.7 natural-log units for response spectra and 0.5 and 0.8 for Fourier spectra. The estimate of model epistemic uncertainty, based on the different model predictions, lies between 0.2 and 0.4, which is about one-half of the estimates for the standard deviation of the combined model uncertainty and random variability. Parametric uncertainty, based on variation of just the average rupture velocity, is shown to be consistent in amplitude with previous estimates, showing percentage changes in ground motion from 50% to 300% when rupture velocity changes from 2.5 to 2.9 km/s. In addition, there is some evidence that mean biases can be reduced by averaging ground-motion estimates from different methods.

Bulletin of the Seismological Society of America↗

USGS2009

Composition and variation of noise recorded at the Yellowknife Seismic Array, 1991-2007

We analyze seismic noise recorded on the 18 short-period, vertical component seismometers of the Yellowknife Seismic Array (YKA). YKA has an aperture of 23 km and is sited on cratonic lithosphere in an area with low cultural noise. These properties make it ideal for studying natural seismic noise at periods of 1-3 s. We calculated frequency-wave number spectra in this band for over 6,000 time windows that were extracted once per day for 17 years (1991-2007). Slowness analysis reveals a rich variety of seismic phases originating from distinct source regions: R g waves from the Great Slave Lake; L g waves from the Atlantic, Pacific, and Arctic Oceans; and teleseismic P waves from the north Pacific and equatorial mid-Atlantic regions. The surface wave energy is generated along coastlines, while the body wave energy is generated at least in part in deep-water, pelagic regions. Surface waves tend to dominate at the longer periods and, just as in earthquake seismograms, L g is the most prominent arrival. Although the periods we study are slightly shorter than the classic double-frequency microseismic band of 4-10 s, the noise at YKA has clear seasonal behavior that is consistent with the ocean wave climate in the Northern Hemisphere. The temporal variation of most of the noise sources can be well fit using just two Fourier components: yearly and biyearly terms that combine to give a fast rise in microseismic power from mid-June through mid-October, followed by a gradual decline. The exception is the R g energy from the Great Slave Lake, which shows a sharp drop in noise power over a 2-week period in November as the lake freezes. The L g noise from the east has a small but statistically significant positive slope, perhaps implying increased ocean wave activity in the North Atlantic over the last 17 years. Copyright 2009 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗