SEARCH · USGS Science

Results for “Seismological Society of America”

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

USGS2009

Scenario liquefaction hazard maps of Santa Clara Valley, Northern California

Maps showing the probability of surface manifestations of liquefaction in the northern Santa Clara Valley were prepared with liquefaction probability curves. These curves were based on complementary cumulative frequency distributions of the liquefaction potential index (LPI) for surficial geologic units in the study area. LPI values were computed with extensive cone penetration test soundings. Maps were developed for three earthquake scenarios, an M 7.8 event on the San Andreas fault comparable to the 1906 event, an M 6.7 event on the Hayward fault comparable to the 1868 event, and an M 6.9 event on the Calaveras fault. Ground motions were estimated with the Boore and Atkinson (2008) attenuation relation. Liquefaction is predicted for all three events in young Holocene levee deposits along the major creeks. Liquefaction probabilities are highest for the M 7.8 earthquake, ranging from 0.33 to 0.37 if a 1.5 m deep water table is assumed, and from 0.10 to 0.14, if a 5 m deep water table is assumed. Liquefaction probabilities of the other surficial geologic units are less than 0.05. Probabilities for the scenario earthquakes are generally consistent with observations during historical earthquakes.

Bulletin of the Seismological Society of America↗

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↗

USGS2009

Seismicity rate changes along the central California coast due to stress changes from the 2003 M 6.5 San Simeon and 2004 M 6.0 Parkfield earthquakes

We investigated the relationship between seismicity rate changes and modeled Coulomb static stress changes from the 2003 M 6.5 San Simeon and the 2004 M 6.0 Parkfield earthquakes in central California. Coulomb stress modeling indicates that the San Simeon mainshock loaded parts of the Rinconada, Hosgri, and San Andreas strike-slip faults, along with the reverse faults of the southern Los Osos domain. All of these loaded faults, except for the San Andreas, experienced a seismicity rate increase at the time of the San Simeon mainshock. The Parkfield earthquake occurred 9 months later on the loaded portion of the San Andreas fault. The Parkfield earthquake unloaded the Hosgri fault and the reverse faults of the southern Los Osos domain, which both experienced seismicity rate decreases at the time of the Parkfield event, although the decreases may be related to the decay of San Simeon-triggered seismicity. Coulomb stress unloading from the Parkfield earthquake appears to have altered the aftershock decay rate of the southern cluster of San Simeon after-shocks, which is deficient compared to the expected number of aftershocks from the Omori decay parameters based on the pre-Parkfield aftershocks. Dynamic stress changes cannot explain the deficiency of aftershocks, providing evidence that static stress changes affect earthquake occurrence. However, a burst of seismicity following the Parkfield earthquake at Ragged Point, where the static stress was decreased, provides evidence for dynamic stress triggering. It therefore appears that both Coulomb static stress changes and dynamic stress changes affect the seismicity rate.

Bulletin of the Seismological Society of America↗

USGS2009

Earthquake stress drops and inferred fault strength on the Hayward Fault, east San Francisco Bay, California

We study variations in earthquake stress drop with respect to depth, faulting regime, creeping versus locked fault behavior, and wall-rock geology. We use the P-wave displacement spectra from borehole seismic recordings of M 1.0-4.2 earthquakes in the east San Francisco Bay to estimate stress drop using a stack-and-invert empirical Green's function method. The median stress drop is 8.7 MPa, and most stress drops are in the range between 0.4 and 130 MPa. An apparent correlation between stress drop and magnitude is entirely an artifact of the limited frequency band of 4-55 Hz. There is a trend of increasing stress drop with depth, with a median stress drop of ~5 MPa for 1-7 km depth, ~10 MPa for 7-13 km depth, and ~50 MPa deeper than 13 km. We use S=P amplitude ratios measured from the borehole records to better constrain the first-motion focal mechanisms. High stress drops are observed for a deep cluster of thrust-faulting earthquakes. The correlation of stress drops with depth and faulting regime implies that stress drop is related to the applied shear stress. We compare the spatial distribution of stress drops on the Hayward fault to a model of creeping versus locked behavior of the fault and find that high stress drops are concentrated around the major locked patch near Oakland. This also suggests a connection between stress drop and applied shear stress, as the locked patch may experience higher applied shear stress as a result of the difference in cumulative slip or the presence of higher-strength material. The stress drops do not directly correlate with the strength of the proposed wall-rock geology at depth, suggesting that the relationship between fault strength and the strength of the wall rock is complex.

Bulletin of the Seismological Society of America↗

USGS2009

Laboratory and field testing of commercial rotational seismometers

There are a small number of commercially available sensors to measure rotational motion in the frequency and amplitude ranges appropriate for earthquake motions on the ground and in structures. However, the performance of these rotational seismometers has not been rigorously and independently tested and characterized for earthquake monitoring purposes as is done for translational strong- and weak-motion seismometers. Quantities such as sensitivity, frequency response, resolution, and linearity are needed for the understanding of recorded rotational data. To address this need, we, with assistance from colleagues in the United States and Taiwan, have been developing performance test methodologies and equipment for rotational seismometers. In this article the performance testing methodologies are applied to samples of a commonly used commercial rotational seismometer, the eentec model R-1. Several examples were obtained for various test sequences in 2006, 2007, and 2008. Performance testing of these sensors consisted of measuring: (1) sensitivity and frequency response; (2) clip level; (3) self noise and resolution; and (4) cross-axis sensitivity, both rotational and translational. These sensor-specific results will assist in understanding the performance envelope of the R-1 rotational seismometer, and the test methodologies can be applied to other rotational seismometers.

Bulletin of the Seismological Society of America↗

USGS2009

A Viscoelastic earthquake simulator with application to the San Francisco Bay region

Earthquake simulation on synthetic fault networks carries great potential for characterizing the statistical patterns of earthquake occurrence. I present an earthquake simulator based on elastic dislocation theory. It accounts for the effects of interseismic tectonic loading, static stress steps at the time of earthquakes, and postearthquake stress readjustment through viscoelastic relaxation of the lower crust and mantle. Earthquake rupture initiation and termination are determined with a Coulomb failure stress criterion and the static cascade model. The simulator is applied to interacting multifault systems: one, a synthetic two-fault network, and the other, a fault network representative of the San Francisco Bay region. The faults are discretized both along strike and along dip and can accommodate both strike slip and dip slip. Stress and seismicity functions are evaluated over 30,000 yr trial time periods, resulting in a detailed statistical characterization of the fault systems. Seismicity functions such as the coefficient of variation and a- and b-values exhibit systematic patterns with respect to simple model parameters. This suggests that reliable estimation of the controlling parameters of an earthquake simulator is a prerequisite to the interpretation of its output in terms of seismic hazard.

California↗

USGS2009

On the use of high-resolution topographic data as a proxy for seismic site conditions (V S30 )

An alternative method has recently been proposed for evaluating global seismic site conditions, or the average shear velocity to 30 m depth (V S30 ), from the Shuttle Radar Topography Mission (SRTM) 30 arcsec digital elevation models (DEMs). The basic premise of the method is that the topographic slope can be used as a reliable proxy for V S30 in the absence of geologically and geotechnically based site-condition maps through correlations between VS30 measurements and topographic gradient. Here we evaluate the use of higher-resolution (3 and 9 arcsec) DEMs to examine whether we are able to resolve V S30 in more detail than can be achieved using the lower-resolution SRTM data. High-quality DEMs at resolutions greater than 30 arcsec are not uniformly available at the global scale. However, in many regions where such data exist, they may be employed to resolve finer-scale variations in topographic gradient, and consequently, V S30 . We use the U.S. Geological Survey Earth Resources Observation and Science (EROS) Data Center's National Elevation Dataset (NED) to investigate the use of high-resolution DEMs for estimating V S30 in several regions across the United States, including the San Francisco Bay area in California, Los Angeles, California, and St. Louis, Missouri. We compare these results with an example from Taipei, Taiwan, that uses 9 arcsec SRTM data, which are globally available. The use of higher-resolution NED data recovers finer-scale variations in topographic gradient, which better correlate to geological and geomorphic features, in particular, at the transition between hills and basins, warranting their use over 30 arcsec SRTM data where available. However, statistical analyses indicate little to no improvement over lower-resolution topography when compared to V S30 measurements, suggesting that some topographic smoothing may provide more stable V S30 estimates. Furthermore, we find that elevation variability in canopy-based SRTM measurements at resolutions greater than 30 arcsec are too large to resolve reliable slopes, particularly in low-gradient sedimentary basins.

Bulletin of the Seismological Society of America↗

USGS2009

Major earthquakes recorded by Speleothems in Midwestern U.S. caves

Historic earthquakes generated by the New Madrid seismic zone represent some of the largest recorded in the United States, yet prehistoric events are recognized only through deformation in late-Wisconsin to Holocene-age, near surface sediments (liquefaction, monoclinal folding, and changes in river meanders). In this article, we show that speleothems in caves of southwestern Illinois and southeastern Missouri may constitute a previously unrecognized recorder of large earthquakes in the U.S. midcontinent region. The timing of the initiation and regrowth of stalagmites in southwestern Illinois and southeastern Missouri caves is consistent with the historic and prehistoric record of several known seismic events in the U.S. midcontinent region. We conclude that dating the initiation of original stalagmite growth and later postearthquake rejuvenation constitutes a new paleoseismic method that has the potential for being applied to any region around the world in the vicinity of major seismic zones where caves exist. Use of this technique could expand the geographical distribution of paleoseimic data, document prehistoric earthquakes, and help improve interpretations of paleoearthquakes.

Bulletin of the Seismological Society of America↗

USGS2009

A case study of two M ~5 mainshocks in Anza, California: Is the footprint of an aftershock sequence larger than we think?

It has been traditionally held that aftershocks occur within one to two fault lengths of the mainshock. Here we demonstrate that this perception has been shaped by the sensitivity of seismic networks. The 31 October 2001 M w 5.0 and 12 June 2005 M w 5.2 Anza mainshocks in southern California occurred in the middle of the densely instrumented ANZA seismic network and thus were unusually well recorded. For the June 2005 event, aftershocks as small as M 0.0 could be observed stretching for at least 50 km along the San Jacinto fault even though the mainshock fault was only &sim;4.5 km long. It was hypothesized that an observed aseismic slipping patch produced a spatially extended aftershock-triggering source, presumably slowing the decay of aftershock density with distance and leading to a broader aftershock zone. We find, however, the decay of aftershock density with distance for both Anza sequences to be similar to that observed elsewhere in California. This indicates there is no need for an additional triggering mechanism and suggests that given widespread dense instrumentation, aftershock sequences would routinely have footprints much larger than currently expected. Despite the large 2005 aftershock zone, we find that the probability that the 2005 Anza mainshock triggered the M 4.9 Yucaipa mainshock, which occurred 4.2 days later and 72 km away, to be only 14%±1%. This probability is a strong function of the time delay; had the earthquakes been separated by only an hour, the probability of triggering would have been 89%.

California↗

USGS2009

Reassessment of probabilistic seismic hazard in the Marmara region

In 1999, the eastern coastline of the Marmara region (Turkey) witnessed increased seismic activity on the North Anatolian fault (NAF) system with two damaging earthquakes (M 7.4 Kocaeli and M 7.2 D??zce) that occurred almost three months apart. These events have reduced stress on the western segment of the NAF where it continues under the Marmara Sea. The undersea fault segments have been recently explored using bathymetric and reflection surveys. These recent findings helped scientists to understand the seismotectonic environment of the Marmara basin, which has remained a perplexing tectonic domain. On the basis of collected new data, seismic hazard of the Marmara region is reassessed using a probabilistic approach. Two different earthquake source models: (1) the smoothed-gridded seismicity model and (2) fault model and alternate magnitude-frequency relations, Gutenberg-Richter and characteristic, were used with local and imported ground-motion-prediction equations. Regional exposure is computed and quantified on a set of hazard maps that provide peak horizontal ground acceleration (PGA) and spectral acceleration at 0.2 and 1.0 sec on uniform firm-rock site condition (760 m=sec average shear wave velocity in the upper 30 m). These acceleration levels were computed for ground motions having 2% and 10% probabilities of exceedance in 50 yr, corresponding to return periods of about 2475 and 475 yr, respectively. The maximum PGA computed (at rock site) is 1.5g along the fault segments of the NAF zone extending into the Marmara Sea. The new maps generally show 10% to 15% increase for PGA, 0.2 and 1.0 sec spectral acceleration values across much of Marmara compared to previous regional hazard maps. Hazard curves and smooth design spectra for three site conditions: rock, soil, and soft-soil are provided for the Istanbul metropolitan area as possible tools in future risk estimates.

Bulletin of the Seismological Society of America↗

USGS2009

The 1911 M ~6.6 Calaveras earthquake: Source parameters and the role of static, viscoelastic, and dynamic coulomb stress changes imparted by the 1906 San Francisco earthquake

The occurrence of a right-lateral strike-slip earthquake in 1911 is inconsistent with the calculated 0.2-2.5 bar static stress decrease imparted by the 1906 rupture at that location on the Calaveras fault, and 5 yr of calculated post-1906 viscoelastic rebound does little to reload the fault. We have used all available first-motion, body-wave, and surface-wave data to explore possible focal mechanisms for the 1911 earthquake. We find that the event was most likely a right-lateral strikeslip event on the Calaveras fault, larger than, but otherwise resembling, the 1984 M w 6.1 Morgan Hill earthquake in roughly the same location. Unfortunately, we could recover no unambiguous surface fault offset or geodetic strain data to corroborate the seismic analysis despite an exhaustive archival search. We calculated the static and dynamic Coulomb stress changes for three 1906 source models to understand stress transfer to the 1911 site. In contrast to the static stress shadow, the peak dynamic Coulomb stress imparted by the 1906 rupture promoted failure at the site of the 1911 earthquake by 1.4-5.8 bar. Perhaps because the sample is small and the aftershocks are poorly located, we find no correlation of 1906 aftershock frequency or magnitude with the peak dynamic stress, although all aftershocks sustained a calculated dynamic stress of ???3 bar. Just 20 km to the south of the 1911 epicenter, we find that surface creep of the Calaveras fault at Hollister paused for ~17 yr after 1906, about the expected delay for the calculated static stress drop imparted by the 1906 earthquake when San Andreas fault postseismic creep and viscoelastic relaxation are included. Thus, the 1911 earthquake may have been promoted by the transient dynamic stresses, while Calaveras fault creep 20 km to the south appears to have been inhibited by the static stress changes.

Bulletin of the Seismological Society of America↗

USGS2010

High-resolution seismic reflection imaging of growth folding and shallow faults beneath the Southern Puget Lowland, Washington State

Marine seismic reflection data from southern Puget Sound, Washington, were collected to investigate the nature of shallow structures associated with the Tacoma fault zone and the Olympia structure. Growth folding and probable Holocene surface deformation were imaged within the Tacoma fault zone beneath Case and Carr Inlets. Shallow faults near potential field anomalies associated with the Olympia structure were imaged beneath Budd and Eld Inlets. Beneath Case Inlet, the Tacoma fault zone includes an &sim;350-m wide section of south-dipping strata forming the upper part of a fold (kink band) coincident with the southern edge of an uplifted shoreline terrace. An &sim;2 m change in the depth of the water bottom, onlapping postglacial sediments, and increasing stratal dips with increasing depth are consistent with late Pleistocene to Holocene postglacial growth folding above a blind fault. Geologic data across a topographic lineament on nearby land indicate recent uplift of late Holocene age. Profiles acquired in Carr Inlet 10 km to the east of Case Inlet showed late Pleistocene or Holocene faulting at one location with &sim;3 to 4 m of vertical displacement, south side up. North of this fault the data show several other disruptions and reflector terminations that could mark faults within the broad Tacoma fault zone. Seismic reflection profiles across part of the Olympia structure beneath southern Puget Sound show two apparent faults about 160 m apart having 1 to 2 m of displacement of subhorizontal bedding. Directly beneath one of these faults, a dipping reflector that may mark the base of a glacial channel shows the opposite sense of throw, suggesting strike-slip motion. Deeper seismic reflection profiles show disrupted strata beneath these faults but little apparent vertical offset, consistent with strike-slip faulting. These faults and folds indicate that the Tacoma fault and Olympia structure include active structures with probable postglacial motion.

Washington↗

USGS2009

Broadband records of earthquakes in deep gold mines and a comparison with results from SAFOD, California

For one week during September 2007, we deployed a temporary network of field recorders and accelerometers at four sites within two deep, seismically active mines. The ground-motion data, recorded at 200 samples/sec, are well suited to determining source and ground-motion parameters for the mining-induced earthquakes within and adjacent to our network. Four earthquakes with magnitudes close to 2 were recorded with high signal/noise at all four sites. Analysis of seismic moments and peak velocities, in conjunction with the results of laboratory stick-slip friction experiments, were used to estimate source processes that are key to understanding source physics and to assessing underground seismic hazard. The maximum displacements on the rupture surfaces can be estimated from the parameter , where is the peak ground velocity at a given recording site, and R is the hypocentral distance. For each earthquake, the maximum slip and seismic moment can be combined with results from laboratory friction experiments to estimate the maximum slip rate within the rupture zone. Analysis of the four M 2 earthquakes recorded during our deployment and one of special interest recorded by the in-mine seismic network in 2004 revealed maximum slips ranging from 4 to 27 mm and maximum slip rates from 1.1 to 6.3 m/sec. Applying the same analyses to an M 2.1 earthquake within a cluster of repeating earthquakes near the San Andreas Fault Observatory at Depth site, California, yielded similar results for maximum slip and slip rate, 14 mm and 4.0 m/sec.

Bulletin of the Seismological Society of America↗

USGS2009

A glossary for rotational seismology

Some common terms for rotational seismology are compiled from contributions by many authors, including R. DeSalvo, J. R. Evans, E. F. Grekova, C. R. Hutt, H. Igel, C. A. Langston, B. Lantz, E. Majewski, R. Nigbor, J. Pujol, P. Spudich, R. Teisseyre, M. D. Trifunac, R. J. Twiss, and Z. Zembaty. In addition, some glossary terms about earthquakes are excerpted from Aki and Lee (2003), Lee and Wu (2009), and others for the benefit of readers who are not seismologists. It is a first attempt to compile a glossary for rotational seismology (including some terms for continuum mechanics), and it is by no means complete. Several iterations will be required to reach a satisfactory glossary, and this compilation is just a start.

Bulletin of the Seismological Society of America↗

USGS2009

Introduction to the special issue on rotational seismology and engineering applications

Because of a lack of suitable instruments, rotational ground motions have not been observed until the last decade. Rotational measurements in the near field of earthquakes in Japan ( Takeo, 1998 ) indicate that rotational ground motions are many times larger than expected from the classical elasticity theory. After failing to obtain useful rotational ground motions (using similar rotational sensors as Takeo did), we deployed a far more sensitive rotational velocity sensor (R-1) at the HGSD station in eastern Taiwan. From 7 December 2004 to 12 November 2006, several hundreds of earthquakes were recorded during our Phase 1 operation. This was mostly a learning exercise to solve field operation problems; Phase 1 operations ended when our two R-1 sensors ceased to operate. A K2+R1 instrument was deployed in the spring of 2007 to start our Phase 2 operation. From 8 May 2007 to 17 February 2008, we observed 52 local earthquakes with good rotational velocity signals (with signal-to-noise ratio >∼5), together with excellent translational acceleration signals (with signal-to-noise ratio >∼10). Unfortunately, field operation was interrupted due to flooding of the HGSD station site in mid-February 2008; we just resumed normal operation in June 2008.

Bulletin of the Seismological Society of America↗

USGS2009

The puzzle of the 1996 Bárdarbunga, Iceland, earthquake: no volumetric component in the source mechanism

A volcanic earthquake with M w 5.6 occurred beneath the Bárdarbunga caldera in Iceland on 29 September 1996. This earthquake is one of a decade-long sequence of events at Bárdarbunga with non-double-couple mechanisms in the Global Centroid Moment Tensor catalog. Fortunately, it was recorded well by the regional-scale Iceland Hotspot Project seismic experiment. We investigated the event with a complete moment tensor inversion method using regional long-period seismic waveforms and a composite structural model. The moment tensor inversion using data from stations of the Iceland Hotspot Project yields a non-double-couple solution with a 67% vertically oriented compensated linear vector dipole component, a 32% double-couple component, and a statistically insignificant (2%) volumetric (isotropic) contraction. This indicates the absence of a net volumetric component, which is puzzling in the case of a large volcanic earthquake that apparently is not explained by shear slip on a planar fault. A possible volcanic mechanism that can produce an earthquake without a volumetric component involves two offset sources with similar but opposite volume changes. We show that although such a model cannot be ruled out, the circumstances under which it could happen are rare.

Bardarbunga↗

USGS2010

Detailed p- and s-wave velocity models along the LARSE II transect, Southern California

Structural details of the crust determined from P-wave velocity models can be improved with S-wave velocity models, and S-wave velocities are needed for model-based predictions of strong ground motion in southern California. We picked P- and S-wave travel times for refracted phases from explosive-source shots of the Los Angeles Region Seismic Experiment, Phase II (LARSE II); we developed refraction velocity models from these picks using two different inversion algorithms. For each inversion technique, we calculated ratios of P- to S-wave velocities (VP/VS) where there is coincident P- and S-wave ray coverage.We compare the two VP inverse velocity models to each other and to results from forward modeling, and we compare the VS inverse models. The VS and VP/VS models differ in structural details from the VP models. In particular, dipping, tabular zones of low VS, or high VP/VS, appear to define two fault zones in the central Transverse Ranges that could be parts of a positive flower structure to the San Andreas fault. These two zones are marginally resolved, but their presence in two independent models lends them some credibility. A plot of VS versus VP differs from recently published plots that are based on direct laboratory or down-hole sonic measurements. The difference in plots is most prominent in the range of VP = 3 to 5 km=s (or VS ~ 1:25 to 2:9 km/s), where our refraction VS is lower by a few tenths of a kilometer per second from VS based on direct measurements. Our new VS - VP curve may be useful for modeling the lower limit of VS from a VP model in calculating strong motions from scenario earthquakes.

Bulletin of the Seismological Society of America↗