USGS Science⌕ Search

SEARCH · USGS Science

Results for “Bulletin of the Seismological Society of America”

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 883 records · Page 49Linked to original sources
USGS1989

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

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

California↗
USGS1989

On the similarity of theories of anelastic and scattering attenuation

We point out basic parallels between theories of anelastic and scattering attenuation. We consider approximations to scattering effects presented by O'Doherty and Anstey (1971), Sato (1982), and Wu (1982). We use the linear theory of anelasticity. We note that the frequency dependence of Q can be related to a distribution of scales of physical properties of the medium. The frequency dependence of anelastic Q is related to the distribution of relaxation times in exactly the same manner as the frequency dependence of scattering Q is related to the distribution of scatterer sizes. Thus, the well-known difficulty of separating scattering from intrinsic attenuation is seen from this point of view as a consequence of the fact that certain observables can be interpreted by identical equations resulting from either of two credible physical theories describing fundamentally different processes.

Bulletin of the Seismological Society of America↗
USGS1990

The theory and prototype development of a stress-monitoring system

A new approach has been developed by which changes in rock stress can be measured directly. The measurement of stress change depends on the reaction of a slender, compliant cavity to the transient variations of rock stress about the cavity. A liquid-filled pressurized cell, emplaced in the cavity, provides the means by which changes in pressure can be monitored. The change in pressure in the cell is equivalent to the change in rock stress in the direction normal to the cell if the compliance contrast between cavity and rock is large and the aspect ratio of the cavity is small. Earthquake-related variations in the stress field in a highly fractured quartz-monzonite stock have been observed to occur in a characteristic way prior to a microseismic event. The precursory stress change (decompression followed by compression) lasted about 5 hours and led directly to a stress drop of 15 × 10−3 bar. Seismic noise monitored during the same period of time suggests that failure in the fractured rock mass resulted from rapid encroachment of meltwater and subsequent increase in fluid pressure.

Bulletin of the Seismological Society of America↗
USGS1990

The 1989 earthquake swarm beneath Mammoth Mountain, California: An initial look at the 4 May through 30 September activity

Mammoth Mountain is a 50,000- to 200,000-year-old cumulovolcano standing on the southwestern rim of Long Valley in eastern California. On 4 May 1989, two M = 1 earthquakes beneath the south flank of the mountain marked the onset of a swarm that has continued for more than 6 months. In addition to its longevity, noteworthy aspects of this persistent swarm include (1) an exponential-like increase in the rate of activity through the first month; (2) a vertically oriented, planar distribution of hypocenters at depths between 6 and 9 km with a north-northeast strike (roughly perpendicular to the average T-axis orientation for the swarm earthquakes); (3) recurring spasmodic bursts (rapid-fire sequences of similar-sized earthquakes with overlapping coda) and occasional earthquakes with enhanced low-frequency energy; (4) a uniform temporal distribution of the four largest (M ≈ 3) events over the first 4 months of the swarm with a cumulative seismic moment for the entire sequence through 30 September corresponding to a single M ≈ 4 earthquake; (5) a b-value of 1.2; and (6) submicrostrain perturbations on the nearby borehole dilatometer, the first of which led the onset of swarm activity by more than 2 weeks. These aspects of the swarm, together with its location along the southern extension of the youthful Mono-Inyo volcanic chain, which last erupted 500 to 600 years ago, point to a magmatic source for the modest but persistent influx of strain energy into the crust beneath Mammoth Mountain.

Bulletin of the Seismological Society of America↗
USGS1990

The 3 December 1988 Pasadena, California earthquake: Evidence for strike-slip motion on the Raymond Fault

The Pasadena earthquake (ML = 4.9) occurred on 3 December 1988, at a depth of 16 km. The hypocenters of the earthquake and its aftershocks define a east-northeast striking, steeply northwest-dipping surface that projects up to the active surficial trace of the Raymond fault. One of the nodal planes of the focal mechanism of the earthquake parallels the Raymond fault with left-lateral strike-slip movement on that plane, and is consistent with geomorphic and paleoseismic evidence that the Raymond fault is dominantly a left-lateral strike-slip fault. The existence of a component of sinistral slip along the Raymond fault had been suspected prior to the earthquake, but the northward dip of the fault and the prominent scarp along the western portion of its trace had led most workers to conclude that slip along the fault was dominantly reverse.

California↗
USGS1990

Source inversion of the 1988 Upland, California, earthquake: Determination of a fault plane for a small event

We examined short-period P waves to investigate if waveform data could be used to determine which of two nodal planes was the actual fault plane for a small (ML 4.6) earthquake near Upland, California. We removed path and site complications by choosing a small aftershock (ML 2.7) as an empirical Green function. The main shock P waves were deconvolved by using the empirical Green function to produce simple far-field displacement pulses. We used a least-squares method to invert these pulses for the slip distribution on a finite fault. Both nodal planes (strike 125°, dip 85° and strike 221°, dip 40°) of the first-motion focal mechanism were tested at various rupture velocities. The southwest trending fault plane consistently gave better fitting solutions than the southeast-trending plane. We determined a moment of 4.2 × 1022 dyne-cm. The rupture velocity, and thus the source area could not be well resolved, but if we assume a reasonable rupture velocity of 0.87 times the shear wave velocity, we obtain a source area of 0.97 km2 and a stress drop of 38 bars. Choice of a southwest-trending fault plane is consistent with the trend of the nearby portion of the Transverse Ranges frontal fault zone and indicates left-lateral motion. This method provides a way to determine the fault plane for small earthquakes that have no surface rupture and no obvious trend in aftershock locations.

California↗
USGS1990

Near-surface velocities and attenuation at two boreholes near Anza, California, from logging data

To investigate near-surface site effects in granite rock, we drilled 300-m-deep boreholes at two sites which are collocated with stations from the digital array at Anza, California. The first borehole was sited at station KNW (Keenwild fire station), which is located along a ridge line about 8.7 km east of the San Jacinto Fault zone. Station PFO (Piñon Flat Observatory), chosen for the second site, is another 6 km further to the east of station KNW and is located on a gently sloping hillside. We logged each borehole for P- and S-wave velocities, as well as for crack density and orientation.

California↗
USGS1990

Calculation of zero-offset vertical seismic profiles generated by a horizontal point force acting on the surface of an elastic half-space

Impulse responses including near-field terms have been obtained in closed form for the zero-offset vertical seismic profiles generated by a horizontal point force acting on the surface of an anelastic half-space. The method is based on the correspondence principle. Through transformation of variables, the Fourier transform of the elastic impulse response is put in a form such that the Fourier transform of the corresponding anelastic impulse response can be expressed as elementary functions and their definite integrals involving distance, angular frequency, phase velocities, and attenuation factors. These results are used for accurate calculation of shear-wave arrival rise times of synthetic seismograms needed for data interpretation of anelastic-attenuation measurements in near-surface sediment.

Bulletin of the Seismological Society of America↗
USGS1990

Source parameters and effects of bandwidth and local geology on high- frequency ground motions observed for aftershocks of the northeastern Ohio earthquake of 31 January 1986

A 10-station array (GEOS) yielded recordings of exceptional bandwidth (400 sps) and resolution (up to 96 dB) for the aftershocks of the moderate ( m b ≈ 4.9) earthquake that occurred on 31 January 1986 near Painesville, Ohio. Nine aftershocks were recorded with seismic moments ranging between 9 × 10 16 and 3 × 10 19 dyne-cm ( M w : 0.6 to 2.3). The two largest aftershocks (depth 5.3, 5.6 km; oblique right slip, rake ≈30°, strike ≈N25°E) yielded seismic signals above background noise at frequencies as high as 130 Hz at epicentral distances up to 17 km. The aftershock recordings at a site underlain by ≈8 m of lakeshore sediments show significant levels of high-frequency soil amplification of vertical motion at frequencies near 8, 20, and 70 Hz. Viscoelastic models for P and SV waves incident at the base of the sediments yield estimates of vertical P -wave response consistent with the observed high-frequency site resonances, but suggest additional detailed shear-wave logs are needed to account for observed S -wave response. Peak acceleration values obtained from the broadband recordings are about two and four times as large as those that would be recorded on strong-motion recorders or short-period networks with upper bandwidth limits of 30 and 15 Hz, respectively. Attenuation-corrected acceleration spectra are used to reduce the influence of high-frequency (up to 100 Hz) local site effects on corner frequency estimates. The moment versus source radius trend inferred for events with moments as small as 9 × 10 16 dyne-cm, based on the Brune source model, extends previous relations inferred for the central United States, shows little evidence for a minimum source radius, and suggests that stress drops for the smaller events ( M 0 < 10 19 dyne-cm) decrease with decreasing moment.

Ohio↗
USGS1990

Coda duration magnitudes in central California: An empirical approach

A new empirical coda magnitude M D is presented for the central California seismic network (CALNET), in which station corrections are introduced to reduce the influence of hypocenter and station distributions. The new magnitude scale is obtained by relating the lapse-time τ^ to local magnitudes M L using data from 55 earthquakes with 1.1 ≦ M L ≤ 5.6 from June 1977 through 1981. Magnitude estimates are significantly improved by introducing corrections for the seismograph instrument attenuation setting α ij and a site correction δ j . The magnitude of the i th event 〈 M D 〉, is the median value of individual M D ij = -1.03 + 2.10log 10 τˆ ij + 0.00268 τˆ ij + α ij + δ j , where τ^ ij is the lapse-time measured at the j th station. Local site effects influence the duration of the coda and cause stations to either overestimate or underestimate 〈 M D 〉, typically up to 0.5 magnitude units. The distribution of the site corrections is spatially correlated, changing coherently across major tectonic and geologic features such as the San Andreas fault. It appears that the site corrections are influenced by the physical attenuation properties near the station, seemingly contradicting the concept that the coda is related to back-scattering and attenuation in a volume with dimensions on the order of the hypocenter to receiver distance.

California↗
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↗
USGS1990

Seismic site effects and the spatial interpolation of earthquake seismograms: Results using aftershocks of the 1986 North Palm Springs, California, earthquake

We address the following two questions. Can a microearthquake's ground motions be modeled by incident P and S waves that excite a site transfer-function that is a smooth function of incidence angle? Given recorded ground motions from a set of earthquakes having known locations and mechanisms, can we derive such a site transfer-function and use it to obtain the ground motions that would result from an earthquake source occurring somewhere in the same volume but having a location and mechanism that are different from the recorded events? Although many factors will cause two distinct microearthquake sources to have different seismograms at a common station, in this paper we concentrate only upon the differences caused by source mechanisms, P- and S-wave travel-time variations and by variations in the site transfer-function. We specifically exclude the effects of waves scattered from heterogeneities in the geologic structure away from the seismic site. We express the site transfer-function as a sum of several terms having simple dependences upon incidence angle and azimuth. Each term is an independent function of time. Given a set of seismograms observed at the site, we solve a linear system of equations for the time dependences of each term. These time series may be used to calculate the seismograms that would have resulted from an earthquake having arbitrary mechanism and location. This step is an interpolation. We have applied this technique to seismograms after aftershocks of the 1986 North Palm Springs earthquake. Our interpolation technique works fairly well within the volume occupied by the recorded events, but the method is not very successful at providing accurate seismograms for sources located outside the aftershock volume. The primary causes of the inaccuracy are the inadequacy of our chosen angular functions to model the site response fully and the likely scattering of seismic waves by geological heterogeneities (in this case, the Banning and Mission Creek faults) near the seismic stations. Our methods could be used to determine the effects of single scattering from lateral heterogeneities in geologic structure.

California↗
USGS1990

Relationship between strong-motion array parameters and the accuracy of source inversion and physical waves

We develop general rules for a strong-motion array layout on the basis of our method of applying a prediction analysis to a source inversion scheme. A systematic analysis is done to obtain a relationship between fault-array parameters and the accuracy of a source inversion. An overdetermined least-squares inversion scheme is used, where the spatial resolution is predetermined by the subfault size and the model variance is used as the accuracy of the source inversion. The accuracy of source inversion is efficiently estimated by using Wolberg's prediction analysis.

Bulletin of the Seismological Society of America↗
USGS1990

The effect of S-wave arrival times on the accuracy of hypocenter estimation

Well-constrained hypocenters (latitude, longitude, depth, and origin time) are required for nearly all studies that use earthquake data. We have examined the theoretical basis behind some of the widely accepted “rules of thumb” for obtaining accurate hypocenter estimates that pertain to the use of S phases and illustrate, in a variety of ways, why and when these “rules” are applicable. Results of experiments done for this study show that epicentral estimates (latitude and longitude) are typically far more robust with respect to data inadequacies; therefore, only examples illustrating the relationship between S phase arrival time data and focal depth and origin time estimates are presented. Most methods used to determine earthquake hypocenters are based on iterative, linearized, least-squares algorithms. Standard errors associated with hypocenter parameters are calculated assuming the data errors may be correctly described by a Gaussian distribution. We examine the influence of S-phase arrival time data on such algorithms by using the program HYPOINVERSE with synthetic datasets. Least-squares hypocenter determination algorithms have several shortcomings: solutions may be highly dependent on starting hypocenters, linearization and the assumption that data errors follow a Gaussian distribution may not be appropriate, and depth/origin time trade-offs are not readily apparent. These shortcomings can lead to biased hypocenter estimates and standard errors that do not always represent the true error. To illustrate the constraint provided by S-phase data on hypocenters determined without some of these potential problems, we also show examples of hypocenter estimates derived using a probabilistic approach that does not require linearization. We conclude that a correctly timed S phase recorded within about 1.4 focal depth's distance from the epicenter can be a powerful constraint on focal depth. Furthermore, we demonstrate that even a single incorrectly timed S phase can result in depth estimates and associated measures of uncertainty that are significantly incorrect.

Bulletin of the Seismological Society of America↗
USGS1991

Application of an iterative least-squares waveform inversion of strong-motion and teleseismic records to the 1978 Tabas, Iran, earthquake

An iterative least-squares technique is used to simultaneously invert the strong-motion records and teleseismic P waveforms for the 1978 Tabas, Iran, earthquake to deduce the rupture history. The effects of using different data sets and different parametrizations of the problem (linear versus nonlinear) are considered. A consensus of all the inversion runs indicates a complex, multiple source for the Tabas earthquake, with four main source regions over a fault length of 90 km and an average rupture velocity of 2.5 km/sec.

Bulletin of the Seismological Society of America↗
USGS1991

The southeastern Illinois earthquake of 10 June 1987: The later aftershocks

The 10 June 1987 southeastern Illinois earthquake (mbLg = 5.2) was the seventh in a series of moderate magnitude (≧ 4.5) MMI VII shocks to occur in the Wabash Valley seismic zone of southeastern Illinois and southwestern Indiana since 1891. Located about 200 km east of St. Louis, Missouri, the shock caused minor damage in the epicentral area, had a contiguous felt area of about 433,000 km2, and had a total felt area over 1 million km2. Within 47 hours after the main shock, a cooperative effort between the U.S. Geological Survey and Memphis State University installed a 15-station aftershock monitoring network (later expanded to 21 instruments) that recorded more than 100 aftershocks in the following 4-day period. Results from the 56 aftershocks that were well located indicate a compact, cylindrically shaped aftershock volume about 1.7 km long, 0.8 km wide, and with a vertical distribution between about 9 and 12 km in depth. Composite focal mechanism solutions of the aftershocks suggest that the predominant mode of faulting is reverse slip, but some strike-slip type motion occurred similar to the mechanism for the main shock as determined from teleseismic data. The maximum principal compressive stress (P axes) is oriented easterly to east-southeasterly and is subhorizontal in plunge. These results, combined with historical observations and the recent discovery of earthquake-induced liquefaction features in the area, imply that the seismic potential of the zone may be greater than generally perceived.

Illinois↗