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 361 records · Page 20Linked to original sources

USGS2013

Ground motion in the presence of complex topography: Earthquake and ambient noise sources

To study the influence of topography on ground motion, eight seismic recorders were deployed for a period of one year over Poverty Ridge on the east side of the San Francisco Bay Area, California. This location is desirable because of its proximity to local earthquake sources and the significant topographic relief of the array (439 m). Topographic amplification is evaluated as a function of frequency using a variety of methods, including reference‐site‐based spectral ratios and single‐station horizontal‐to‐vertical spectral ratios using both shear waves from earthquakes and ambient noise. Field observations are compared with the predicted ground motion from an accurate digital model of the topography and a 3D local velocity model. Amplification factors from the theoretical calculations are consistent with observations. The fundamental resonance of the ridge is prominently observed in the spectra of data and synthetics; however, higher‐frequency peaks are also seen primarily for sources in line with the major axis of the ridge, perhaps indicating higher resonant modes. Excitations of lateral ribs off of the main ridge are also seen at frequencies consistent with their dimensions. The favored directions of resonance are shown to be transverse to the major axes of the topographic features.

California↗

USGS2004

Radiated energy and the rupture process of the Denali fault earthquake sequence of 2002 from broadband teleseismic body waves

Displacement, velocity, and velocity-squared records of P and SH body waves recorded at teleseismic distances are analyzed to determine the rupture characteristics of the Denali fault, Alaska, earthquake of 3 November 2002 ( M W 7.9, M e 8.1 ). Three episodes of rupture can be identified from broadband ( ∼0.1–5.0 Hz ) waveforms. The Denali fault earthquake started as a M W 7.3 thrust event. Subsequent right-lateral strike-slip rupture events with centroid depths of 9 km occurred about 22 and 49 sec later. The teleseismic P waves are dominated by energy at intermediate frequencies (0.1-1 Hz) radiated by the thrust event, while the SH waves are dominated by energy at lower frequencies (0.05-0.2 Hz) radiated by the strike-slip events. The strike-slip events exhibit strong directivity in the teleseismic SH waves. Correcting the recorded P -wave acceleration spectra for the effect of the free surface yields an estimate of 2.8 × 10 15 N m for the energy radiated by the thrust event. Correcting the recorded SH -wave acceleration spectra similarly yields an estimate of 3.3 × 10 16 N m for the energy radiated by the two strike-slip events. The average rupture velocity for the strike-slip rupture process is 1.1 β –1.2 β . The strike-slip events were located 90 and 188 km east of the epicenter. The rupture length over which significant or resolvable energy is radiated is, thus, far shorter than the 340-km fault length over which surface displacements were observed. However, the seismic moment released by these three events, 4 × 10 20 N m, was approximately half the seismic moment determined from very low-frequency analyses of the earthquake. The difference in seismic moment can be reasonably attributed to slip on fault segments that did not radiate significant or coherent seismic energy. These results suggest that very large and great strike-slip earthquakes can generate stress pulses that rapidly produce substantial slip with negligible stress drop and little discernible radiated energy on fault segments distant from the initial point of nucleation. The existence of this energy-deficient rupture mode has important implications for the evaluation of the seismic hazard of very large strike-slip earthquakes.

Alaska↗

USGS1999

Attenuation and source properties at the Coso Geothermal area, California

We use a multiple-empirical Green's function method to determine source properties of small ( M −0.4 to 1.3) earthquakes and P - and S -wave attenuation at the Coso Geothermal Field, California. Source properties of a previously identified set of clustered events from the Coso geothermal region are first analyzed using an empirical Green's function ( EGF ) method. Stress-drop values of at least 0.5-1 MPa are inferred for all of the events; in many cases, the corner frequency is outside the usable bandwidth, and the stress drop can only be constrained as being higher than 3 MPa. P - and S -wave stress-drop estimates are identical to the resolution limits of the data. These results are indistinguishable from numerous EGF studies of M 2-5 earthquakes, suggesting a similarity in rupture processes that extends to events that are both tiny and induced, providing further support for Byerlee's Law. Whole-path Q estimates for P and S waves are determined using the multiple-empirical Green's function ( MEGF ) method of Hough (1997), whereby spectra from clusters of colocated events at a given station are inverted for a single attenuation parameter, κ, with source parameters constrained from EGF analysis. The κ estimates, which we infer to be resolved to within 0.01 sec or better, exhibit almost as much scatter as a function of hypocentral distance as do values from previous single-spectrum studies for which much higher uncertainties in individual κ estimates are expected. The variability in κ estimates determined here therefore suggests real lateral variability in Q structure. Although the ray-path coverage is too sparse to yield a complete three-dimensional attenuation tomographic image, we invert the inferred κ value for three-dimensional structure using a damped least-squares method, and the results do reveal significant lateral variability in Q structure. The inferred attenuation variability corresponds to the heat-flow variations within the geothermal region. A central low- Q region corresponds well with the central high-heat flow region; additional detailed structure is also suggested.

California↗

USGS2018

Opportunities to enhance seismic demand parameters for future editions of the AS1170.4

Geoscience Australia has recently released its 2018 National Seismic Hazard Assessment (NSHA18). Results from the NSHA18 indicate significantly lower seismic hazard across almost all Australian localities at the 1/500 annual exceedance probability level relative to the factors adopted for the current Australian Standard AS1170.4–2007 (R2018). These new hazard estimates, coupled with larger probability factors (kp) for long return periods, have challenged notions of seismic hazard in Australia in terms of the recurrence of damaging ground motions. As a consequence, the new hazard estimates have raised questions over the appropriateness of the prescribed National Construction Code probability level as used in the AS1170.4 to determine appropriate seismic demands for the design of ordinary-use structures. Therefore, it is suggested that the ground-motion exceedance probability used in the current AS1170.4 be reviewed in light of the recent hazard assessment and the expected performance of modern buildings for rarer ground motions. Whilst adjusting the AS1170.4 exceedance probability level would be a major departure from previous earthquake loading standards, it would bring it into line with other international building codes in similar tectonic environments. Additionally, it would offer opportunities to further modernise how seismic demands are considered in Australian building design. In particular, the authors highlight the following additional opportunities: 1) the use of uniform hazard spectra to replace and simplify the spectral shape factors, which do not deliver uniform hazard across all natural periods; 2) updated site amplification factors to ensure continuity with modern ground-motion models, and; 3) the potential to define design ground motions in terms of uniform collapse risk rather than uniform hazard. Estimation of seismic hazard at any location is an uncertain science. However, as our knowledge improves, our estimates of the hazard will converge more closely to the actual – but unknowable – (time independent) hazard. It is therefore prudent to regularly update the estimates of the seismic demands in our building codes using the best available evidence-based methods and models.

Conference Paper↗

USGS2009

Along-Arc and Back-Arc Attenuation, Site Response, and Source Spectrum for the Intermediate-Depth 8 January 2006 M 6.7 Kythera, Greece, Earthquake

An M 6.7 intermediate-depth (66 km), in-slab earthquake occurring near the island of Kythera in Greece on 8 January 2006 was well recorded on networks of stations equipped with acceleration sensors and with broadband velocity sensors. All data were recorded digitally using recording instruments with resolutions ranging from almost 11 to 24 bits. We use data from these networks to study the distance dependence of the horizontal-component Fourier acceleration spectra (FAS) and horizontal-component pseudoabsolute response spectral acceleration (PSA). For purposes of simulating motions in the future, we parameterize the distance decay using several forms of the geometrical-spreading function, for each of which we derive Q as a function of frequency. By extrapolating the distance decay back to 1 km, we obtain a reference spectrum that can be used in future simulations. This spectrum requires a more complicated spectral shape than the classic single-corner-frequency model; in particular, there appears to be an enhancement of motion around 0.2-0.3 Hz that may be due to the radiation of a 3-5 sec pulse from the source. We infer a ?? 0 value of about 0.055 sec for rock stations and a stress parameter in the range of 400-600 bars. We also find distinctive differences in the site response of stations on soft soil and soil; both the FAS and the 5% damped PSA amplifications have similar peak amplitudes (about 2 and 4 for soil and soft-soil sites, respectively, relative to the rock sites) at similar frequencies (between about 0.4 and 2.0 Hz, with the soft-soil amplifications peaking at somewhat lower frequencies than the soil amplifications). One of the most distinctive features of the data is the clear difference in the motions for along-arc and back-arc stations, with the former being significantly higher than the latter over a broad range of frequencies at distances beyond about 250 km. The motions from the Kythera earthquake are roughly comparable to those from intermediate-depth earthquakes elsewhere, but they appear to be significantly higher than those from recordings of shallow earthquakes in Greece of comparable magnitude and hypocentral distance.

Bulletin of the Seismological Society of America↗

USGS2004

Acceleration and volumetric strain generated by the Parkfield 2004 earthquake on the GEOS strong-motion array near Parkfield, CA

An integrated array of 11 General Earthquake Observation System (GEOS) stations installed near Parkfield, CA provided on scale broad-band, wide-dynamic measurements of acceleration and volumetric strain of the Parkfield earthquake (M 6.0) of September 28, 2004. Three component measurements of acceleration were obtained at each of the stations. Measurements of collocated acceleration and volumetric strain were obtained at four of the stations. Measurements of velocity at most sites were on scale only for the initial P-wave arrival. When considered in the context of the extensive set of strong-motion recordings obtained on more than 40 analog stations by the California Strong-Motion Instrumentation Program (Shakal, et al., 2004 http://www.quake.ca.gov/cisn-edc) and those on the dense array of Spudich, et al, (1988), these recordings provide an unprecedented document of the nature of the near source strong motion generated by a M 6.0 earthquake. The data set reported herein provides the most extensive set of near field broad band wide dynamic range measurements of acceleration and volumetric strain for an earthquake as large as M 6 of which the authors are aware. As a result considerable interest has been expressed in these data. This report is intended to describe the data and facilitate its use to resolve a number of scientific and engineering questions concerning earthquake rupture processes and resultant near field motions and strains. This report provides a description of the array, its scientific objectives and the strong-motion recordings obtained of the main shock. The report provides copies of the uncorrected and corrected data. Copies of the inferred velocities, displacements, and Psuedo velocity response spectra are provided. Digital versions of these recordings are accessible with information available through the internet at several locations: the National Strong-Motion Program web site (http://agram.wr.usgs.gov/), the COSMOS Virtual Data Center Web site (http://www.cosmos-eq.org), and the CISN Engineering and Berkeley data centers (http://www.quake.ca.gov/cisn-edc). They are also accessible together with recordings on the GEOS Strong-motion Array near Parkfield, CA since its installation in 1987 through the USGS GEOS web site ( http://nsmp.wr.usgs.gov/GEOS).

California↗

USGS1995

Relative performance of several inexpensive accelerometers

We examined the performance of several low-cost accelerometers for highly cost-driven applications in recording earthquake strong motion. We anticipate applications for such sensors in providing the lifeline and emergency-response communities with an immediate, comprehensive picture of the extent and characteristics of likely damage. We also foresee their use as 'filler' instruments sited between research-grade instruments to provide spatially detailed and near-field records of large earthquakes (on the order of 1000 stations at 600-m intervals in San Fernando Valley, population 1.2 million, for example). The latter applications would provide greatly improved attenuation relationships for building codes and design, the first examples of mainshock information (that is, potentially nonlinear regime) for microzonation, and a suite of records for structural engineers. We also foresee possible applications in monitoring structural inter-story drift during earthquakes, possibly leading to local and remote alarm functions as well as design criteria. This effort appears to be the first of its type at the USGS. It is spurred by rapid advances in sensor technology and the recognition of potential non-classical applications. In this report, we estimate sensor noise spectra, relative transfer functions and cross-axis sensitivity of six inexpensive sensors. We tested three micromachined ('silicon-chip') sensors in addition to classical force-balance and piezoelectric examples. This sample of devices is meant to be representative, not comprehensive. Sensor noise spectra were estimated by recording system output with the sensor mounted on a pneumatically supported 545-kg optical-bench isolation table. This isolation table appears to limit ground motion to below our system noise level. These noise estimates include noise introduced by signal-conditioning circuitry, the analog-to-digital converter (ADC), and noise induced in connecting wiring by ambient electromagnetic fields in our suburban laboratory. These latter sources are believed to dominate sensor noise in the quieter sensors we tested. Transfer functions were obtained relative to a research grade force-balance accelerometer (a Kinemetrics TM FBA-11) by shaking the sensors simultaneously on the same shake table and taking spectral ratios with the output of the FBA- 11. This reference sensor is said to have 120 db dynamic range (-+20 bits, though we only digitized it at 16 bits resolution and drove it with relatively small signals). We did not test temperature sensitivity, which is thought to be a significant issue at least for the silicon devices. Though these tests were not designed to be definitive (our anticipated applications do not demand research-grade precision), our tests do appear to have been successful in estimating relative transfer functions from about 0.3 to 50 Hz. Most sensors performed adequately in this range, with essentially fiat relative transfer functions. Noise tests appear to measure sensor noise well for the noisier (generally less expensive) instruments from about 0.1 to 50 Hz.

Open-File Report↗

USGS2005

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

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

Bulletin of the Seismological Society of America↗

USGS2003

Estimation of ground motion for Bhuj (26 January 2001; Mw 7.6) and for future earthquakes in India

Only five moderate and large earthquakes ( M w ≥5.7) in India—three in the Indian shield region and two in the Himalayan arc region—have given rise to multiple strong ground-motion recordings. Near-source data are available for only two of these events. The Bhuj earthquake ( M w 7.6), which occurred in the shield region, gave rise to useful recordings at distances exceeding 550 km. Because of the scarcity of the data, we use the stochastic method to estimate ground motions. We assume that (1) S waves dominate at R < 100 km and Lg waves at R ≥ 100 km, (2) Q = 508 f 0.48 is valid for the Indian shield as well as the Himalayan arc region, (3) the effective duration is given by fc -1 + 0.05R, where fc is the corner frequency, and R is the hypocentral distance in kilometer, and (4) the acceleration spectra are sharply cut off beyond 35 Hz. We use two finite-source stochastic models. One is an approximate model that reduces to the ω 2 -source model at distances greater that about twice the source dimension. This model has the advantage that the ground motion is controlled by the familiar stress parameter, Δ σ . In the other finite-source model, which is more reliable for near-source ground-motion estimation, the high-frequency radiation is controlled by the strength factor, sfact , a quantity that is physically related to the maximum slip rate on the fault. We estimate Δ σ needed to fit the observed Amax and Vmax data of each earthquake (which are mostly in the far field). The corresponding sfact is obtained by requiring that the predicted curves from the two models match each other in the far field up to a distance of about 500 km. The results show: (1) The Δ σ that explains Amax data for shield events may be a function of depth, increasing from ∼50 bars at 10 km to ∼400 bars at 36 km. The corresponding sfact values range from 1.0-2.0. The Δ σ values for the two Himalayan arc events are 75 and 150 bars ( sfact = 1.0 and 1.4). (2) The Δ σ required to explain Vmax data is, roughly, half the corresponding value for Amax, while the same sfact explains both sets of data. (3) The available far-field Amax and Vmax data for the Bhuj mainshock are well explained by Δ σ = 200 and 100 bars, respectively, or, equivalently, by sfact = 1.4. The predicted Amax and Vmax in the epicentral region of this earthquake are 0.80 to 0.95 g and 40 to 55 cm/sec, respectively.

Bulletin of the Seismological Society of America↗

USGS2009

Inducing in situ, nonlinear soil response applying an active source

[1] It is well known that soil sites have a profound effect on ground motion during large earthquakes. The complex structure of soil deposits and the highly nonlinear constitutive behavior of soils largely control nonlinear site response at soil sites. Measurements of nonlinear soil response under natural conditions are critical to advancing our understanding of soil behavior during earthquakes. Many factors limit the use of earthquake observations to estimate nonlinear site response such that quantitative characterization of nonlinear behavior relies almost exclusively on laboratory experiments and modeling of wave propagation. Here we introduce a new method for in situ characterization of the nonlinear behavior of a natural soil formation using measurements obtained immediately adjacent to a large vibrator source. To our knowledge, we are the first group to propose and test such an approach. Employing a large, surface vibrator as a source, we measure the nonlinear behavior of the soil by incrementally increasing the source amplitude over a range of frequencies and monitoring changes in the output spectra. We apply a homodyne algorithm for measuring spectral amplitudes, which provides robust signal-to-noise ratios at the frequencies of interest. Spectral ratios are computed between the receivers and the source as well as receiver pairs located in an array adjacent to the source, providing the means to separate source and near-source nonlinearity from pervasive nonlinearity in the soil column. We find clear evidence of nonlinearity in significant decreases in the frequency of peak spectral ratios, corresponding to material softening with amplitude, observed across the array as the source amplitude is increased. The observed peak shifts are consistent with laboratory measurements of soil nonlinearity. Our results provide constraints for future numerical modeling studies of strong ground motion during earthquakes.

Journal of Geophysical Research B: Solid Earth↗

USGS2022

Estimates of k0 and effects on ground motions in the San Francisco Bay area

Ground‐motion studies are a key component of seismic hazard analyses and often rely on information of the source, path, and site. Extensive research has been done on each of these parameters; however, site‐specific studies are of particular interest to seismic hazard studies, especially in the field of earthquake engineering, as near‐site conditions can have a significant impact on the resulting ground motion at a site. There has been much focus on the constraint of site parameters and their application to seismic hazard studies, especially in the development of ground‐motion models (GMMs). Kappa is an observational parameter describing the high‐frequency attenuation of spectra, and its site contribution ( ⁠ κ 0 "> κ 0 ) has shown to be a good predictor of high‐frequency ground motions; however, measurements are often limited. In this study, we develop a κ 0 "> κ 0 dataset for the San Francisco Bay area (SFBA) by estimating κ 0 "> κ 0 for 228 stations, and we produce a continuous regional map of κ 0 "> κ 0 ⁠ . We find κ 0 "> κ 0 to range between 0.003 and 0.072 s, with larger values concentrating on the east, north, and south sides of the bay, and lower values concentrating on the west side. We also evaluate the robustness of κ 0 "> κ 0 as a site parameter and find it to correlate with peak ground acceleration. These estimates of κ 0 "> κ 0 can add predictive power to GMMs, thus increasing the accuracy of predicted ground motion and improving the robustness of ground‐motion studies in the SFBA.

California↗

USGS1985

SENSITIVITY OF STRUCTURAL RESPONSE TO GROUND MOTION SOURCE AND SITE PARAMETERS.

Designing structures to withstand earthquakes requires an accurate estimation of the expected ground motion. While engineers use the peak ground acceleration (PGA) to model the strong ground motion, seismologists use physical characteristics of the source and the rupture mechanism, such as fault length, stress drop, shear wave velocity, seismic moment, distance, and attenuation. This study presents a method for calculating response spectra from seismological models using random vibration theory. It then investigates the effect of various source and site parameters on peak response. Calculations are based on a nonstationary stochastic ground motion model, which can incorporate all the parameters both in frequency and time domains. The estimation of the peak response accounts for the effects of the non-stationarity, bandwidth and peak correlations of the response.

Conference Paper↗

USGS1986

SITE AMPLIFICATION OF EARTHQUAKE GROUND MOTION.

When analyzing the patterns of damage in an earthquake, physical parameters of the total earthquake-site-structure system are correlated with the damage. Soil-structure interaction, the cause of damage in many earthquakes, involves the frequency-dependent response of both the soil-rock column and the structure. The response of the soil-rock column (called site amplification) is controversial because soil has strain-dependent properties that affect the way the soil column filters the input body and surface seismic waves, modifying the amplitude and phase spectra and the duration of the surface ground motion.

Conference Paper↗

USGS2017

Development and utilization of USGS ShakeCast for rapid post-earthquake assessment of critical facilities and infrastructure

The ShakeCast system is an openly available, near real-time post-earthquake information management system. ShakeCast is widely used by public and private emergency planners and responders, lifeline utility operators and transportation engineers to automatically receive and process ShakeMap products for situational awareness, inspection priority, or damage assessment of their own infrastructure or building portfolios. The success of ShakeCast to date and its broad, critical-user base mandates improved software usability and functionality, including improved engineering-based damage and loss functions. In order to make the software more accessible to novice users—while still utilizing advanced users’ technical and engineering background—we have developed a “ShakeCast Workbook”, a well documented, Excel spreadsheet-based user interface that allows users to input notification and inventory data and export XML files requisite for operating the ShakeCast system. Users will be able to select structure based on a minimum set of user-specified facility (building location, size, height, use, construction age, etc.). “Expert” users will be able to import user-modified structural response properties into facility inventory associated with the HAZUS Advanced Engineering Building Modules (AEBM). The goal of the ShakeCast system is to provide simplified real-time potential impact and inspection metrics (i.e., green, yellow, orange and red priority ratings) to allow users to institute customized earthquake response protocols. Previously, fragilities were approximated using individual ShakeMap intensity measures (IMs, specifically PGA and 0.3 and 1s spectral accelerations) for each facility but we are now performing capacity-spectrum damage state calculations using a more robust characterization of spectral deamnd.We are also developing methods for the direct import of ShakeMap’s multi-period spectra in lieu of the assumed three-domain design spectrum (at 0.3s for constant acceleration; 1s or 3s for constant velocity and constant displacement at very long response periods). As part of ongoing ShakeCast research and development, we will also explore the use of ShakeMap IM uncertainty estimates and evaluate the assumption of employing multiple response spectral damping values rather than the single 5%-damped value currently employed. Developing and incorporating advanced fragility assignments into the ShakeCast Workbook requires related software modifications and database improvements; these enhancements are part of an extensive rewrite of the ShakeCast application.

Conference Paper↗

USGS2006

Site response and attenuation in the Puget Lowland, Washington State

Simple spectral ratio (SSR) and horizontal-to-vertical (HN) site-response estimates at 47 sites in the Puget Lowland of Washington State document significant attenuation of 1.5- to 20-Hz shear waves within sedimentary basins there. Amplitudes of the horizontal components of shear-wave arrivals from three local earthquakes were used to compute SSRs with respect to the average of two bedrock sites and H/V spectral ratios with respect to the vertical component of the shear-wave arrivals at each site. SSR site-response curves at thick basin sites show peak amplifications of 2 to 6 at frequencies of 3 to 6 Hz, and decreasing spectra amplification with increasing frequency above 6 Hz. SSRs at nonbasin sites show a variety of shapes and larger resonance peaks. We attribute the spectral decay at frequencies above the amplification peak at basin sites to attenuation within the basin strata. Computing the frequency-independent, depth-dependent attenuation factor (Qs,int) from the SSR spectral decay between 2 and 20 Hz gives values of 5 to 40 for shallow sedimentary deposits and about 250 for the deepest sedimentary strata (7 km depth). H/V site responses show less spectral decay than the SSR responses but contain many of the same resonance peaks. We hypothesize that the H/V method yields a flatter response across the frequency spectrum than SSRs because the H/V reference signal (vertical component of the shear-wave arrivals) has undergone a degree of attenuation similar to the horizontal component recordings. Correcting the SSR site responses for attenuation within the basins by removing the spectral decay improves agreement between SSR and H/V estimates.

Bulletin of the Seismological Society of America↗

USGS2000

SMSIM--Fortran programs for simulating ground motions from earthquakes: Version 2.0.--a revision of OFR 96-80-A

A simple and powerful method for simulating ground motions is based on the assumption that the amplitude of ground motion at a site can be specified in a deterministic way, with a random phase spectrum modified such that the motion is distributed over a duration related to the earthquake magnitude and to distance from the source. This method of simulating ground motions often goes by the name "the stochastic method." It is particularly useful for simulating the higher-frequency ground motions of most interest to engineers, and it is widely used to predict ground motions for regions of the world in which recordings of motion from damaging earthquakes are not available. This simple method has been successful in matching a variety of ground-motion measures for earthquakes with seismic moments spanning more than 12 orders of magnitude. One of the essential characteristics of the method is that it distills what is known about the various factors affecting ground motions (source, path, and site) into simple functional forms that can be used to predict ground motions. SMSIM is a set of programs for simulating ground motions based on the stochastic method. This Open-File Report is a revision of an earlier report (Boore, 1996) describing a set of programs for simulating ground motions from earthquakes. The programs are based on modifications I have made to the stochastic method first introduced by Hanks and McGuire (1981). The report contains source codes, written in Fortran, and executables that can be used on a PC. Programs are included both for time-domain and for random vibration simulations. In addition, programs are included to produce Fourier amplitude spectra for the models used in the simulations and to convert shear velocity vs. depth into frequency-dependent amplification. The revision to the previous report is needed because the input and output files have changed significantly, and a number of new programs have been included in the set.

Open-File Report↗

USGS1999

Site response for Seattle and source parameters of earthquakes in the Puget Sound Region

We analyzed seismograms from 21 earthquakes ( M L 2.0-4.9) recorded by digital seismographs we deployed in urban Seattle to determine site response and earthquake stress drops. The seismometers were situated on a wide variety of geologic units, including artificial fill (e.g., Kingdome, Harbor Island), Pleistocene age soils (glacial till and outwash deposits of Seattle's hills), modified land (downtown Seattle, Space Needle), and Tertiary sedimentary rock. Two mainshock-aftershock sequences were recorded: the June 1997 Bremerton sequence (mainshock M L 4.9) and the February 1997 South Seattle sequence (mainshock M L 3.5), along with other events in the Puget Sound region. We developed a new inversion procedure to estimate site response, source corner frequencies, and seismic moments from the S -wave spectra. This inversion uses corner frequencies determined from spectral ratios of mainshock-aftershock pairs as constraints. The site responses found from the inversion are not relative to the rock site but are relative to an idealized site with a flat frequency response. The response of the rock site is also found from the inversion. The inversion results show high response for the sites on artificial fill, more moderate amplification for most sites on stiff Pleistocene soils or modified land, and low response for the rock site. Some sites display resonances, such as a strong 2-Hz resonance at our site near the Kingdome, which is caused by the surficial layers of fill and younger alluvium. The sites in West Seattle exhibit high amplification, even though they are on relatively stiff soils of glacial outwash. This may be partly caused by basin surface waves produced by conversion of incident S waves. This high response in West Seattle is consistent with damage reports from the 1949 ( m b 7.1) and 1965 ( m b 6.5) earthquakes. Stress-drop estimates for the events we recorded were generally low, between 0.4 and 25 bars, although some of the events may have had higher stress drops that could not be resolved because of the limited passband. We calculated a stress drop of 24 bars for the Bremerton mainshock and 10 bars for the South Seattle mainshock.

Washington↗

USGS2001

Observations of earthquake source parameters at 2 km depth in the Long Valley Caldera, eastern California

To investigate seismic source parameter scaling and seismic efficiency in the Long Valley caldera, California, we measured source parameters for 41 earthquakes (M 0.5 to M 5) recorded at 2 km depth in the Long Valley Exploratory Well. Borehole recordings provide a wide frequency bandwidth, typically 1 to 200–300 Hz, and greatly reduce seismic noise and path effects compared to surface recordings. We calculated source parameters in both the time and frequency domains for P and S waves. At frequencies above the corner frequency, spectra decay faster than ω 3 , indicating that attenuation plays an important role in shaping the spectra (path averaged Q p = 100–400, Q s = 200–800). Source parameters are corrected for attenuation and radiation pattern. Both static stress drops and apparent stresses range from approximately 0.01 to 30 MPa. Although static stress drops do not vary with seismic moment for these data, our analyses are consistent with apparent stress increasing with increasing moment. To estimate tectonic driving stress and seismic efficiencies in the region, we combined source parameter measurements with knowledge of the stress field and a Coulomb failure criterion to infer a driving stress of 40–70 MPa. Subsequent seismic efficiencies are consistent with McGarr's (1999) hypothesis of a maximum seismic efficiency of 6%.

California↗