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 973 records · Page 54Linked to original sources

Spectral structure in the Earth's microseismic background between 20 and 40 seconds

The long-term behavior of a narrow band peak near 26 sec in the microseismic background of the Earth has been studied over many years at several sites. The amplitude of this peak has been determined to be a function of the Earth's weather seasons; it is larger during the southern hemisphere winter. In addition, the microseismic background between 20 and 40 sec has been scanned for peaks in the spectra; three more lower-level peaks have been identified with the possible presence of more.

Bulletin of the Seismological Society of America↗

Coseismic temporal changes of slip direction: The effect of absolute stress on dynamic rupture

We investigate the dynamics of rupture at low-stress level. We show that one main difference between the dynamics of high- and low-stress events is the amount of coseismic temporal rake rotation occurring at given points on the fault. Curved striations on exposed fault surfaces and earthquake dislocation models derived from ground-motion inversion indicate that the slip direction may change with time at a point on the fault during dynamic rupture. We use a 3D boundary integral method to model temporal rake variations during dynamic rupture propagation assuming a slip-weakening friction law and isotropic friction. The points at which the slip rotates most are characterized by an initial shear stress direction substantially different from the average stress direction over the fault plane. We show that for a given value of stress drop, the level of initial shear stress (i.e., the fractional stress drop) determines the amount of rotation in slip direction. We infer that seismic events that show evidence of temporal rake rotations are characterized by a low initial shear-stress level with spatially variable direction on the fault (possibly due to changes in fault surface geometry) and an almost complete stress drop. Our models motivate a new interpretation of curved and cross-cutting striations and put new constraints on their analysis. The initial rake is in general collinear with the initial stress at the hypocentral zone, supporting the assumptions made in stress-tensor inversion from first-motion analysis. At other points on the fault, especially away from the hypocenter, the initial slip rake may not be collinear with the initial shear stress, contradicting a common assumption of structural geology. On the other hand, the later part of slip in our models is systematically more aligned with the average stress direction than the early slip. Our modeling suggests that the length of the straight part of curved striations is usually an upper bound of the slip-weakening distance if this parameter is uniform over the fault plane, and the direction of the late part of slip of curved striations should have more weight in the estimate of initial stress direction.

Bulletin of the Seismological Society of America↗

Bayesian statistics applied to the location of the source of explosions at Stromboli Volcano, Italy

We present a method for determining the location and spatial extent of the source of explosions at Stromboli Volcano, Italy, based on a Bayesian inversion of the slowness vector derived from frequency-slowness analyses of array data. The method searches for source locations that minimize the error between the expected and observed slowness vectors. For a given set of model parameters, the conditional probability density function of slowness vectors is approximated by a Gaussian distribution of expected errors. The method is tested with synthetics using a five-layer velocity model derived for the north flank of Stromboli and a smoothed velocity model derived from a power-law approximation of the layered structure. Application to data from Stromboli allows for a detailed examination of uncertainties in source location due to experimental errors and incomplete knowledge of the Earth model. Although the solutions are not constrained in the radial direction, excellent resolution is achieved in both transverse and depth directions. Under the assumption that the horizontal extent of the source does not exceed the crater dimension, the 90% confidence region in the estimate of the explosive source location corresponds to a small volume extending from a depth of about 100 m to a maximum depth of about 300 m beneath the active vents, with a maximum likelihood source region located in the 120- to 180-m-depth interval.

Stromboli Volcano↗

Seismic-reflection evidence that the hayward fault extends into the lower crust of the San Francisco Bay Area, California

This article presents deep seismic-reflection data from an experiment across San Francisco Peninsula in 1995 using large (125 to 500 kg) explosive sources. Shot gathers show a mostly nonreflective upper crust in both the Franciscan and Salinian terranes (juxtaposed across the San Andreas fault), an onset of weak lower-crustal reflectivity beginning at about 6-sec two-way travel time ( TWTT ) and bright southwest-dipping reflections between 11 and 13 sec TWTT . Previous studies have shown that the Moho in this area is no deeper than 25 km (∼8 to 9 sec TWTT ). Three-dimensional reflection travel-time modeling of the 11 to 13 sec events from the shot gathers indicates that the bright events may be explained by reflectors 15 to 20 km into the upper mantle, northeast of the San Andreas fault. However, upper mantle reflections from these depths were not observed on marine-reflection profiles collected in San Francisco Bay, nor were they reported from a refraction prifile on San Francisco Peninsula. The most consistent interpretation of these events from 2D raytracing and 3D travel-time modeling is that they are out-of-plane reflections from a high-angle (dipping ∼70° to the southwest) impedance contrast in the lower crust that corresponds with the surface trace of the Hayward fault. These results suggest that the Hayward fault truncates the horizontal detachment fault suggested to be active beneath San Francisco Bay.

California↗

Using safety inspection data to estimate shaking intensity for the 1994 Northridge earthquake

We map the shaking intensity suffered in Los Angeles County during the 17 January 1994, Northridge earthquake using municipal safety inspection data. The intensity is estimated from the number of buildings given red, yellow, or green tags, aggregated by census tract. Census tracts contain from 200 to 4000 residential buildings and have an average area of 6 km 2 but are as small as 2 and 1 km 2 in the most densely populated areas of the San Fernando Valley and downtown Los Angeles, respectively. In comparison, the zip code areas on which standard MMI intensity estimates are based are six times larger, on average, than the census tracts. We group the buildings by age (before and after 1940 and 1976), by number of housing units (one, two to four, and five or more), and by construction type, and we normalize the tags by the total number of similar buildings in each census tract. We analyze the seven most abundant building categories. The fragilities (the fraction of buildings in each category tagged within each intensity level) for these seven building categories are adjusted so that the intensity estimates agree. We calibrate the shaking intensity to correspond with the modified Mercalli intensities ( MMI ) estimated and compiled by Dewey et al. (1995); the shapes of the resulting isoseismals are similar, although we underestimate the extent of the MMI = 6 and 7 areas. The fragility varies significantly between different building categories (by factors of 10 to 20) and building ages (by factors of 2 to 6). The post-1940 wood-frame multi-family (≧5 units) dwellings make up the most fragile building category, and the post-1940 woodframe single-family dwellings make up the most resistant building category.

California↗

Faulting apparently related to the 1994 Northridge, California, earthquake and possible co-seismic origin of surface cracks in Potrero Canyon, Los Angeles County, California

Apparent southward-dipping, reverse-fault zones are imaged to depths of about 1.5 km beneath Potrero Canyon, Los Angeles County, California. Based on their orientation and projection to the surface, we suggest that the imaged fault zones are extensions of the Oak Ridge fault. Geologic mapping by others and correlations with seismicity studies suggest that the Oak Ridge fault is the causative fault of the 17 January 1994 Northridge earthquake (Northridge fault). Our seismically imaged faults may be among several faults that collectively comprise the Northridge thrust fault system. Unusually strong shaking in Potrero Canyon during the Northridge earthquake may have resulted from focusing of seismic energy or co-seismic movement along existing, related shallow-depth faults. The strong shaking produced ground-surface cracks and sand blows distributed along the length of the canyon. Seismic reflection and refraction images show that shallow-depth faults may underlie some of the observed surface cracks. The relationship between observed surface cracks and imaged faults indicates that some of the surface cracks may have developed from nontectonic alluvial movement, but others may be fault related. Immediately beneath the surface cracks, P-wave velocities are unusually low (<400 m/sec), and there are velocity anomalies consistent with a seismic reflection image of shallow faulting to depths of at least 100 m. On the basis of velocity data, we suggest that unconsolidated soils (<800 m/sec) extend to depths of about 15 to 20 m beneath our datum (<25 m below ground surface). The underlying rocks range in velocity from about 1000 to 5000 m/sec in the upper 100 m. This study illustrates the utility of high-resolution seismic imaging in assessing local and regional seismic hazards.

California↗

Variability in nonlinear sediment response during the 1994 Northridge, California, earthquake

Spectral ratios between sediment sites and rock sites for strong and weak ground motion are compared to infer the extent of nonlinear sediment response in the Los Angeles area during the 1994 Northridge earthquake. Northridge mainshock records form the strong-motion database and aftershock records form the weak-motion database. The degree of nonlinearity is evaluated at 25 locations in and around the San Fernando and Los Angeles Basins based on the ratio of weak-to-strong spectral amplification factors (WSAF). Values of WSAF vary greatly from site to site and, for a given site, from one frequency band to another over the range from 0.75 to 20 Hz. This complexity is attributed to differences in several factors: material properties under the sites and the amplitude and frequency content of the input ground motions. Given these observations, it is not possible to assign an acceleration or velocity level at which nonlinearity initiates for all the stiff soil sites analyzed. However, no nonlinearity is seen below a peak acceleration of about 200 to 300 cm/sec 2 or peak velocity of 20 to 30 cm/sec. Nonlinearity is seen at all sites in or near the San Fernando Valley. Nonlinearity in the Los Angeles Basin is limited to a few isolated sites. For the Northridge earthquake, even sites with strong nonlinearity still show amplification levels greater than 1.0 relative to rock sites for frequencies below 10 Hz.

California↗

Deaggregation of probabilistic ground motions in the central and eastern United States

Probabilistic seismic hazard analysis (PSHA) is a technique for estimating the annual rate of exceedance of a specified ground motion at a site due to known and suspected earthquake sources. The relative contributions of the various sources to the total seismic hazard are determined as a function of their occurrence rates and their ground-motion potential. The separation of the exceedance contributions into bins whose base dimensions are magnitude and distance is called deaggregation . We have deaggregated the hazard analyses for the new USGS national probabilistic ground-motion hazard maps (Frankel et al. , 1996). For points on a 0.2° grid in the central and eastern United States (CEUS) , we show color maps of the geographical variation of mean and modal magnitudes ( M , Mˇ ) and distances ( D , Dˇ ) for ground motions having a 2% chance of exceedance in 50 years. These maps are displayed for peak horizontal acceleration and for spectral response accelerations of 0.2, 0.3, and 1.0 sec. We tabulate M , D , Mˇ , and Dˇ for 49 CEUS cities for 0.2- and 1.0-sec response. Thus, these maps and tables are PSHA-derived estimates of the potential earthquakes that dominate seismic hazard at short and intermediate periods in the CEUS.

Bulletin of the Seismological Society of America↗

Deformation and quaternary faulting in southeast Missouri across the Commerce geophysical lineament

High-resolution seismic-reflection data acquired at three sites along the surface projection of the Commerce geophysical lineament in southeast Missouri reveal a complex history of post-Cretaceous faulting that has continued into the Quaternary. Near Qulin, Missouri, approximately 20 m of apparent vertical fault displacement has occurred in the Quaternary. Reflection data collected at Idalia Hill, about 45 km to the northeast, reveal a series of reverse and possibly right-lateral strike-slip faults with Quaternary displacement. In the Benton Hills, 45 km northeast of Idalia Hill, seismic data image a complicated series of anticlinal and synclinal fault-bounded blocks immediately north of the Commerce fault. We infer that most of the deformation imaged in the upper 400 m of these three data sets occurred since post-Cretaceous time, and a significant portion of it occurred during Quaternary time. Collectively, these seismic data along with geomorphic and surface-geologic evidence suggest (1) the existence of at least one potential seismogenic structure in southeastern Missouri outside the main zones of New Madrid seismicity, and (2) these structures have been active during the Quaternary. The geographic location of the imaged deformation suggests it is related to structures along with the Commerce geophysical lineament.

Missouri↗

Analysis of earthquake recordings obtained from the Seafloor Earthquake Measurement System (SEMS) instruments deployed off the coast of southern California

For more than 20 years, a program has been underway to obtain records of earthquake shaking on the seafloor at sites offshore of southern California, near oil platforms. The primary goal of the program is to obtain data that can help determine if ground motions at offshore sites are significantly different than those at onshore sites; if so, caution may be necessary in using onshore motions as the basis for the seismic design of oil platforms. We analyze data from eight earthquakes recorded at six offshore sites; these are the most important data recorded on these stations to date. Seven of the earthquakes were recorded at only one offshore station; the eighth event was recorded at two sites. The earthquakes range in magnitude from 4.7 to 6.1. Because of the scarcity of multiple recordings from any one event, most of the analysis is based on the ratio of spectra from vertical and horizontal components of motion. The results clearly show that the offshore motions have very low vertical motions compared to those from an average onshore site, particularly at short periods. Theoretical calculations find that the water layer has little effect on the horizontal components of motion but that it produces a strong spectral null on the vertical component at the resonant frequency of P waves in the water layer. The vertical-to-horizontal ratios for a few selected onshore sites underlain by relatively low shear-wave velocities are similar to the ratios from offshore sites for frequencies less than about one-half the water layer P -wave resonant frequency, suggesting that the shear-wave velocities beneath a site are more important than the water layer in determining the character of the ground motions at lower frequencies.

California↗

Basin waves on a seafloor recording of the 1990 Upland, California, earthquake: Implications for ground motions from a larger earthquake

The velocity and displacement time series from a recording on the seafloor at 74 km from the 1990 Upland earthquake ( M = 5.6) are dominated by late-arriving waves with periods of 6 to 7 sec. These waves are probably surface waves traveling across the Los Angeles basin. Response spectra for the recording are in agreement with predictions from empirical regression equations and theoretical models for periods less than about 1 sec but are significantly larger than those predictions for longer periods. The longer-period spectral amplitudes are controlled by the late-arriving waves, which are not included in the theoretical models and are underrepresented in the data used in the empirical analyses. When the motions are scaled to larger magnitude, the results are in general agreement with simulations of wave propagation in the Los Angeles basin by Graves (1998).

California↗

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↗

Seismic activity of the San Francisco Bay region

Moment magnitude M with objective confidence-level uncertainties are estimated for felt San Francisco Bay region earthquakes using Bakun and Went-worth's (1997) analysis strategy for seismic intensity observations. The frequency-magnitude distribution is well described for M ≧ 5.5 events since 1850 by a Gutenberg-Richter relation with a b -value of 0.90. The seismic moment rate Σ M 0 /yr since 1836 is 2.68 × 10 18 N-m/yr (95% confidence range = 1.29 × 10 18 N-m/yr to 4.07 × 10 18 N-m/yr); the seismic moment rate since 1850 is nearly the same. Σ M 0 /yr in the 56 years before 1906 is about 10 times that in the 70 years after 1906. In contrast, Σ M 0 /yr since 1977 is about equal that in the 56 years before 1906. 80% (1σ = 14%) of the plate-motion moment accumulation rate is available for release in earthquakes. The historical Σ M 0 /yr and the portion of the plate-motion moment accumulation rate available for release in earthquakes are used in a seismic cycle model to estimate the rate of seismic activity in the twenty-first century. High and low rates of future seismic activity are both permissible given the range of possible seismic-cycle recurrence times T and the uncertainties in the historical Σ M 0 and in the percentage of plate motion available for release in earthquakes. If the historical seismic moment rate is not greater than the estimated 2.68 × 10 18 N-m/yr and the percentage of the plate-motion moment accumulation available for release in earthquakes is not less than the estimated 80%, then for all T , the rate of seismic moment release from now until the next 1906-sized shock will be comparable to the rate from 1836 to 1905 when M 6 1/2 shocks occurred every 15 to 20 years.

California↗

Test of two methods for faulting on finite-difference calculations

Tests of two fault boundary conditions show that each converges with second order accuracy as the finite-difference grid is refined. The first method uses split nodes so that there are disjoint grids that interact via surface traction. The 3D version described here is a generalization of a method I have used extensively in 2D; it is as accurate as the 2D version. The second method represents fault slip as inelastic strain in a fault zone. Offset of stress from its elastic value is seismic moment density. Implementation of this method is quite simple in a finite-difference scheme using velocity and stress as dependent variables.

Bulletin of the Seismological Society of America↗

Foreshocks and aftershocks of the great 1857 California earthquake

The San Andreas fault is the longest fault in California and one of the longest strike-slip faults anywhere in the world, yet we know little about many aspects of its behavior before, during, and after large earthquakes. We conducted a study to locate and to estimate magnitudes for the largest foreshocks and aftershocks of the 1857 M 7.9 Fort Tejon earthquake on the central and southern segments of the fault. We began by searching archived first-hand accounts from 1857 through 1862, by grouping felt reports temporally, and by assigning modified Mercalli intensities to each site. We then used a modified form of the grid-search algorithm of Bakun and Wentworth, derived from empirical analysis of modern earthquakes, to find the location and magnitude most consistent with the assigned intensities for each of the largest events. The result confirms a conclusion of Sieh that at least two foreshocks (“dawn” and “sunrise”) located on or near the Parkfield segment of the San Andreas fault preceded the mainshock. We estimate their magnitudes to be M ≈ 6.1 and M ≈ 5.6, respectively. The aftershock rate was below average but within one standard deviation of the number of aftershocks expected based on statistics of modern southern California mainshock-aftershock sequences. The aftershocks included two significant events during the first eight days of the sequence, with magnitudes M ≈ 6.25 and M ≈ 6.7, near the southern half of the rupture; later aftershocks included a M ≈ 6 event near San Bernardino in December 1858 and a M ≈ 6.3 event near the Parkfield segment in April 1860. From earthquake logs at Fort Tejon, we conclude that the aftershock sequence lasted a minimum of 3.75 years.

California↗

Nonlinear soil response in the vicinity of the Van Normal Complex following the 1994 Northridge, California, earthquake

Ground-motion recording obtaineds at the Van Norman Complex from the 1994 Northridge, California, mainshock and its aftershocks constitute an excellent data set for the analysis of soil response as a function of ground-motion amplitude. We searched for nonlinear response by comparing the Fourier spectral ratios of two pairs of sites for ground motions of different levels, using data from permanent strong-motion recorders and from specially deployed portable instruments. We also compared the amplitude dependence of the observed ratios with the amplitude dependence of the theoretical ratios obtained from 1-D linear and 1-D equivalent-linear transfer functions, using recently published borehole velocity profiles at the sites to provide the low-strain material properties. One pair of sites was at the Jensen Filtration Plant ( JFP ); the other pair was the Rinaldi Receiving Station ( RIN ) and the Los Angeles Dam ( LAD ). Most of the analysis was concentrated on the motions at the Jensen sites. Portable seismometers were installed at the JFP to see if the motions inside the structures housing the strong-motion recorders differed from nearby free-field motions. We recorded seven small earthquakes and found that the high-frequency, low-amplitude motions in the administration building were about 0.3 of those outside the building. This means that the lack of high frequencies on the strong-motion recordings in the administration building relative to the generator building is not due solely to nonlinear soil effects. After taking into account the effects of the buildings, however, analysis of the suite of strong- and weak-motion recordings indicates that nonlinearity occurred at the JFP . As predicted by equivalent-linear analysis, the largest events (the mainshock and the 20 March 1994 aftershock) show a significant deamplification of the high-frequency motion relative to the weak motions from aftershocks occurring many months after the mainshock. The weak-motion aftershocks recorded within 12 hours of the mainshock, however, show a relative deamplification similar to that in the mainshock. The soil behavior may be a consequence of a pore pressure buildup during large-amplitude motion that was not dissipated until sometime later. The motions at ( RIN ) and ( LAD ) are from free-field sites. The comparison among spectral ratios of the mainshock, weak-motion coda waves of the mainshock, and an aftershock within ten minutes of the mainshock indicate that some nonlinearity occurred, presumably at ( RIN ) because it is the softer site. The spectral ratio for the mainshock is between that calculated for pure linear response and that calculated from the equivalent-linear method, using commonly used modulus reduction and damping ratio curves. In contrast to the Jensen sites, the ratio of motions soon after the high-amplitude portion of the mainshock differs from the ratio of the mainshock motions, indicating the mechanical properties of the soil returned to the low-strain values as the high-amplitude motion ended. This may indicate a type of nonlinear soil response different from that affecting motion at the Jensen administration building.

California↗

Fault-slip distribution of the 1995 Colima-Jalisco, Mexico, earthquake

Broadband teleseismic P waves have been analyzed to recover the rupture history of the large ( M S 7.4) Colima-Jalisco, Mexico, shallow interplate thrust earthquake of 9 October 1995. Ground-displacement records in the period range of 1-60 sec are inverted using a linear, finite-fault waveform inversion procedure that allows a variable dislocation duration on a prescribed fault. The method is applied using both a narrow fault that simulates a line source with a dislocation window of 50 sec and a wide fault with a possible rise time of up to 20 sec that additionally allows slip updip and downdip from the hypocenter. The line-source analysis provides a spatio-temporal image of the slip distribution consisting of several large sources located northwest of the hypocenter and spanning a range of rupture velocities. The two-dimensional finite-fault inversion allows slip over this rupture-velocity range and indicates that the greatest coseismic displacement (3-4 m) is located between 70 and 130 km from the hypocenter at depths shallower than about 15 km. Slip in this shallow region consists of two major sources, one of which is delayed by about 10 sec relative to a coherent propagation of rupture along the plate interface. These two slip sources account for about one-third of the total P -wave seismic moment of 8.3 × 10 27 dyne-cm ( M w 7.9) and may have been responsible for the local tsunami observed along the coast following the earthquake.

Colima-Jalisco↗