Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Bulletin Seismological Society 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 181 records · Page 10Linked to original sources
USGS2004

Frequency-dependent Lg Q within the continental United States

Frequency-dependent crustal attenuation (1/ Q ) is determined for seven distinct physiographic/tectonic regions of the continental United States using high-quality Lg waveforms recorded on broadband stations in the frequency band 0.5 to 16 Hz. Lg attenuation is determined from time-domain amplitude measurements in one-octave frequency bands centered on the frequencies 0.75, 1.0, 3.0, 6.0, and 12.0 Hz. Modeling errors are determined using a delete-j jackknife resampling technique. The frequency-dependent quality factor is modeled in the form of Q = Q 0 f η. Regions were initially selected based on tectonic provinces but were eventually limited and adjusted to maximize ray path coverage in each area. Earthquake data was recorded on several different networks and constrained to events occurring within the crust (<40 km depth) and at least m b 3.5 in size. A singular value decomposition inversion technique was applied to the data to simultaneously solve for source and receiver terms along with Q for each region at specific frequencies. The lowest crustal Q was observed in northern and southern California where Q is described by the functions Q = 152(±37) f 0.72(±0.16) and Q = 105(±26) f 0.67(±0.16) , respectively. The Basin and Range Province, Pacific Northwest, and Rocky Mountain states also display lower Q and a strong frequency dependence characterized by the functions Q = 200(±40) f 0.68(±0.12) , Q = 152(±49) f 0.76(±0.18) , and Q = 166(±37) f 0.61(±0.14) , respectively. In contrast, in the central and northeast United States Q functions are Q = 640(±225) f 0.344(±0.22) and Q = 650(±143) f 0.36(±0.14) , respectively, show a high crustal Q and a weaker frequency dependence. These results improve upon previous Lg modeling by subdividing the United States into smaller, distinct tectonic regions and using significantly more data that provide improved constraints on frequency-dependent attenuation and errors. A detailed attenuation map of the continental United States can provide significant input into hazard map mitigation. Both scattering and intrinsic attenuation mechanisms are likely to play a comparable role in the frequency range considered in the study.

Bulletin of the Seismological Society of America↗
USGS2004

Paleoseismology and neotectonics of the Shivwits section of the Hurricane Fault, northwestern Arizona

The Shivwits section of the Hurricane Fault in northwestern Arizona has been largely ignored in evaluating the seismic hazard posed to the rapidly growing populations of southwestern Utah. To assess this hazard, we conducted studies along the Shivwits section using field observations and geomorphic modeling to understand the Quaternary tectonism of this portion of the Hurricane Fault. We have found evidence that it ruptured with up to 2 to 3 m of vertical displacement per event and likely produced ???M 7 earthquakes. Our results suggest that the slip rate along the southern Hurricane Fault has not decreased during the Quaternary. The Moriah Knoll basalt, offset 150 to 200 m along the Hurricane Fault, yielded a maximum long-term slip rate of 0.15 to 0.25 mm/yr, estimated using a new 40Ar/39Ar age of 0.85 ?? 0.06 Ma. The late Quaternary slip rates on alluvial fan surfaces offset 2 to 7.4 m were estimated using pedogenic carbonate rind thickness as a calibrated proxy for age. These observations yielded slip rates of ??? 0.05 to 0.3 mm/yr. Paleoseismic trench investigations showed that two surface-rupturing events occurred in the past 15-78 k.y. A radiocarbon sample from a most-recent event (MRE) fissure yielded a calibrated age of 8900-10,400 years B.P.; the penultimate event was likely ??? 10 k.y. before the most recent event. Slip-rate estimates using the Moriah Knoll basalt (???0. 15-0.24 mm/yr; 850 ka), surface offset (???0.05-0.3 mm/yr; < 100 ka), morphologic modeling (??? 0.06-0.21 mm/yr; <100 ka), and observations from the trench (??? 0.06-0.34 mm/yr; 15-75 ka) suggest that there has been no detectable change in slip rate in the past 1 million years or so. This implies a constant deformation rate for this portion of the Colorado Plateau Margin; therefore, Basin and Range extension is actively encroaching.

Bulletin of the Seismological Society of America↗
USGS2004

Revisiting the 23 February 1892 Laguna Salada earthquake

According to some compilations, the Laguna Salada, Baja California, earthquake of 23 February 1892 ranks among the largest earthquakes in California and Baja California in historic times. Although surface rupture was not documented at the time of the earthquake, recent geologic investigations have identified and mapped a rupture on the Laguna Salada fault that can be associated with high probability with the 1892 event (Mueller and Rockwell, 1995). The only intensity-based magnitude estimate for the earthquake, M 7.8, was made by Strand (1980) based on an interpretation of macroseismic effects and a comparison of isoseismal areas with those from instrumentally recorded earthquakes. In this study we reinterpret original accounts of the Laguna Salada earthquake. We assign modified Mercalli intensity (MMI) values in keeping with current practice, focusing on objective descriptions of damage rather than subjective human response and not assigning MMI values to effects that are now known to be poor indicators of shaking level, such as liquefaction and rockfalls. The reinterpreted isoseismal contours and the estimated magnitude are both significantly smaller than those obtained earlier. Using the method of Bakun and Wentworth (1997) we obtain a magnitude estimate of M 7.2 and an optimal epicenter less than 15 km from the center of the mapped Laguna Salada rupture. The isoseismal contours are elongated toward the northwest, which is qualitatively consistent with a directivity effect, assuming that the fault ruptured from southeast to northwest. We suggest that the elongation may also thus reflect wave propagation effects, with more efficient propagation of crustal surface (Lg) waves in the direction of the overall regional tectonic fabric.

Bulletin of the Seismological Society of America↗
USGS2004

Prediction of nonlinear soil effects

Mathematical models of soil nonlinearity in common use and recently developed nonlinear codes compared to investigate the range of their predictions. We consider equivalent linear formulations with and without frequency-dependent moduli and damping ratios and nonlinear formulations for total and effective stress. Average velocity profiles to 150 m depth with midrange National Earthquake Hazards Reduction Program site classifications (B, BC, C, D, and E) in the top 30 m are used to compare the response of a wide range of site conditions from rock to soft soil. Nonlinear soil models are compared using the amplification spectrum, calculated as the ratio of surface ground motion to the input motion at the base of the velocity profile. Peak input motions from 0.1 g to 0.9 g are considered. For site class B, no significant differences exist between the models considered in this article. For site classes BC and C, differences are small at low input motions (0.1 g to 0.2 g ), but become significant at higher input levels. For site classes D and E the overdamping of frequencies above about 4 Hz by the equivalent linear solution with frequency-independent parameters is apparent for the entire range of input motions considered. The equivalent linear formulation with frequency-dependent moduli and damping ratios under damps relative to the nonlinear models considered for site class C with larger input motions and most input levels for site classes D and E. At larger input motions the underdamping for site classes D and E is not as severe as the overdamping with the frequency-independent formulation, but there are still significant differences in the time domain. A nonlinear formulation is recommended for site classes D and E and for site classes BC and C with input motions greater than a few tenths of the acceleration of gravity. The type of nonlinear formulation to use is driven by considerations of the importance of water content and the availability of laboratory soils data. Our average amplification curves from a nonlinear effective stress formulation compare favorably with observed spectral amplification at class D and E sites in the Seattle area for the 2001 Nisqually earthquake.

Bulletin of the Seismological Society of America↗
USGS2005

Unusually low rates of slip on the Santa Rosa Range fault zone, northern Nevada

The Santa Rosa Range fault zone (SRRFZ) is one of the most topographically prominent normal fault systems in the northern Basin and Range province of the western United States. It has been assigned high rates of vertical slip by others and has been identified as a possible site of the future extension of the central Nevada seismic belt (CNSB). We use detailed trench mapping and luminescence dating to estimate displacements and timing of the last several large-magnitude paleoearthquakes on the southern part of the SRRFZ at a trench site near Orovada, Nevada. Coseismic vertical displacements ranged from 1 to 2.8 m for each of the last four events. Luminescence ages provide time limits for the last three events of 125-155 ka, 90-108 ka, and 11-16 ka. These data yield recurrence intervals of 17-65 k.y. and 74-97 k.y. and an elapsed time of 11-16 k.y. since the youngest event. Slip-rate determinations at the Orovada site are complicated by multiple fault strands, but rates calculated from a variety of data are surprisingly low (0.01-0.16 mm/yr), given the topographic prominence of the Santa Rosa Range. A lack of compelling patterns in a comparison of paleoseismic parameters indicate that the SRRFZ is no more likely a location for a large-magnitude earthquake than previously identified seismic gaps or along faults that lie directly north of the CNSB.

Bulletin of the Seismological Society of America↗
USGS2005

Strong ground motion in the Taipei basin from the 1999 Chi-Chi, Taiwan, earthquake

The Taipei basin, located in northwest Taiwan about 160 km from the epicenter of the Chi-Chi earthquake, is a shallow, triangular-shaped basin filled with low-velocity fluvial deposits. There is a strong velocity contrast across the basement interface of about 600 m/sec at a depth of about 600-700 m in the deeper section of the basin, suggesting that ground motion should be amplified at sites in the basin. In this article, the ground-motion recordings are analyzed to determine the effect of the basin both in terms of amplifications expected from a 1D model of the sediments in the basin and in terms of the 3D structure of the basin. Residuals determined for peak acceleration from attenuation curves are more positive (amplified) in the basin (average of 5.3 cm/ sec2 compared to - 24.2 cm/sec2 for those stations outside the basin and between 75 and 110 km from the surface projection of the faulted area, a 40% increase in peak ground acceleration). Residuals for peak velocity are also significantly more positive at stations in the basin (31.8 cm/sec compared to 20.0 cm/sec out). The correlation of peak motion with depth to basement, while minor in peak acceleration, is stronger in the peak velocities. Record sections of ground motion from stations in and around the Taipei basin show that the largest long-period arrival, which is coherent across the region, is strongest on the vertical component and has a period of about 10-12 sec. This phase appears to be a Rayleigh wave, probably associated with rupture at the north end of the Chelungpu fault. Records of strong motion from stations in and near the basin have an additional, higher frequency signal: nearest the deepest point in the basin, the signal is characterized by frequencies of about 0.3 - 0.4 Hz. These frequencies are close to simple predictions using horizontal layers and the velocity structure of the basin. Polarizations of the S wave are mostly coherent across the array, although there are significant differences along the northwest edge that may indicate large strains across that edge of the basin. The length of each record after the main S wave are all longer at basin stations compared to those outside. This increase in duration of ground shaking is probably caused by amplification of ground motion at basin stations, although coda Q (0.67 - 1.30 Hz) is slightly larger inside the basin compared to those at local stations outside the basin. Durations correlate with depth to basement. These motions are in the range that can induce damage in buildings and may have contributed to the structural collapse of multistory buildings in the Taipei basin.

Bulletin of the Seismological Society of America↗
USGS2005

Coal-mining seismicity and ground-shaking hazard: A case study in the Trail Mountain area, Emery County, Utah

We describe a multipart study to quantify the potential ground-shaking hazard to Joes Valley Dam, a 58-m-high earthfill dam, posed by mining-induced seismicity (MIS) from future underground coal mining, which could approach as close as ???1 km to the dam. To characterize future MIS close to the dam, we studied MIS located ???3-7 km from the dam at the Trail Mountain coal mine. A 12-station local seismic network (11 stations above ground, one below, combining eight triaxial accelerometers and varied velocity sensors) was operated in the Trail Mountain area from late 2000 through mid-2001 for the dual purpose of (1) continuously monitoring and locating MIS associated with longwall mining at a depth of 0.5-0.6 km and (2) recording high-quality data to develop ground-motion prediction equations for the shallow MIS. (Ground-motion attenuation relationships and moment-tensor results are reported in companion articles.) Utilizing a data set of 1913 earthquakes (M ??? 2.2), we describe space-time-magnitude distributions of the observed MIS and source-mechanism information. The MIS was highly correlated with mining activity both in space and time. Most of the better-located events have depths constrained within ??0.6 km of mine level. For the preponderance (98%) of the 1913 located events, only dilatational P-wave first motions were observed, consistent with other evidence for implosive or collapse-type mechanisms associated with coal mining in this region. We assess a probable maximum magnitude of M 3.9 (84th percentile of a cumulative distribution) for potential MIS close to Joes Valley Dam based on both the worldwide and regional record of coal-mining-related MIS and the local geology and future mining scenarios.

Bulletin of the Seismological Society of America↗
USGS2005

Development of ground-motion prediction equations relevant to shallow-mining-induced seismicity in the Trial Mountain area, Emery County, Utah

To provide a basis for assessing the seismic hazard to the Joes Valley Dam due to future coal mining in the nearby Cottonwood Tract, central Utah, we developed ground-motion prediction relations using data recorded by a seismic network, established and operated by the University of Utah Seismograph Stations. The network was centered on the Trail Mountain coal mine, located adjacent to the Cottonwood Tract. From late 2000 until early 2001, this network recorded numerous mining-induced events with magnitudes as large as 2.17. The ground motion from these events, recorded at hypocentral distances ranging from about 500 m to approximately 10 km, were well suited to developing new ground-motion prediction relations, especially when augmented by data from a M 4.2 earthquake in the Willow Creek mine, about 50 km north of Trail Mountain. Using a two-stage regression analysis, we determined prediction relations for peak acceleration, peak velocity, and pseudovelocity response spectra, at 5% damping, for periods of 0.1, 0.2, 0.5, 1.0, and 2.0 s. To illustrate the potential seismic hazard at the Joes Valley dam, we used these ground-motion relations to predict a peak velocity of 6.8 cm/s due to an earthquake with the probable maximum magnitude of 3.9, at a hypocentral distance of 1 km, recorded at a rock site typical for this region. This result does not take into account the site response at the dam.

Bulletin of the Seismological Society of America↗
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↗
USGS2005

Upper-crustal structure of the inner Continental Borderland near Long Beach, California

A new P-wave velocity/structural model for the inner Continental Borderland (ICB) region was developed for the area near Long Beach, California. It combines controlled-source seismic reflection and refraction data collected during the 1994 Los Angeles Region Seismic Experiment (LARSE), multichannel seismic reflection data collected by the U.S. Geological Survey (1998-2000), and nearshore borehole stratigraphy. Based on lateral velocity contrasts and stratigraphic variation determined from borehole data, we are able to locate major faults such as the Cabrillo, Palos Verdes, THUMS-Huntington Beach, and Newport Inglewood fault zones, along with minor faults such as the slope fault, Avalon knoll, and several other yet unnamed faults. Catalog seismicity (1975-2002) plotted on our preferred velocity/structural model shows recent seismicity is located on 16 out of our 24 faults, providing evidence for continuing concern with respect to the existing seismic-hazard estimates. Forward modeling of P-wave arrival times on the LARSE line 1 resulted in a four-layer model that better resolves the stratigraphy and geologic structures of the ICB and also provides tighter constraints on the upper-crustal velocity structure than previous modeling of the LARSE data. There is a correlation between the structural horizons identified in the reflection data with the velocity interfaces determined from forward modeling of refraction data. The strongest correlation is between the base of velocity layer 1 of the refraction model and the base of the planar sediment beneath the shelf and slope determined by the reflection model. Layers 2 and 3 of the velocity model loosely correlate with the diffractive crust layer, locally interpreted as Catalina Schist.

Bulletin of the Seismological Society of America↗
USGS2005

Evidence for a global seismic-moment release sequence

Temporal clustering of the larger earthquakes (foreshock-mainshock-aftershock) followed by relative quiescence (stress shadow) are characteristic of seismic cycles along plate boundaries. A global seismic-moment release history, based on a little more than 100 years of instrumental earthquake data in an extended version of the catalog of Pacheco and Sykes (1992), illustrates similar behavior for Earth as a whole. Although the largest earthquakes have occurred in the circum-Pacific region, an analysis of moment release in the hemisphere antipodal to the Pacific plate shows a very similar pattern. Monte Carlo simulations confirm that the global temporal clustering of great shallow earthquakes during 1952-1964 at M ??? 9.0 is highly significant (4% random probability) as is the clustering of the events of M ??? 8.6 (0.2% random probability) during 1950-1965. We have extended the Pacheco and Sykes (1992) catalog from 1989 through 2001 using Harvard moment centroid data. Immediately after the 1950-1965 cluster, significant quiescence at and above M 8.4 begins and continues until 2001 (0.5% random probability). In alternative catalogs derived by correcting for possible random errors in magnitude estimates in the extended Pacheco-Sykes catalog, the clustering of M ??? 9 persists at a significant level. These observations indicate that, for great earthquakes, Earth behaves as a coherent seismotectonic system. A very-large-scale mechanism for global earthquake triggering and/or stress transfer is implied. There are several candidates, but so far only viscoelastic relaxation has been modeled on a global scale.

Bulletin of the Seismological Society of America↗
USGS2005

Neotectonics of the offshore Oak Ridge fault near Ventura, southern California

The Oak Ridge fault is a large-offset, south-dipping reverse fault that forms the south boundary of the Ventura Basin in southern California. Previous research indicates that the Oak Ridge fault south of the town of Ventura has been inactive since 200-400 ka ago and that the fault tip is buried by ??? 1 km of Quaternary sediment. However, very high-resolution and medium-resolution seismic reflection data presented here show a south-dipping fault, on strike with the Oak Ridge fault, that is truncated at 80 m depth by an unconformity that is probably at the base of late Pleistocene and Holocene sediment. Furthermore, if vertically aligned features in seismic reflection data are eroded remnants of fault scarps, then a subsidiary fault within the Oak Ridge system deforms the shallowest imaged sediment layers. We propose that this subsidiary fault has mainly left-slip offset. These observations of Holocene slip on the Oak Ridge fault system suggest that revision of the earthquake hazard for the densely populated Santa Clara River valley and the Oxnard coastal plain may be needed.

Bulletin of the Seismological Society of America↗
USGS2005

Calculation of broadband time histories of ground motion, Part II: Kinematic and dynamic modeling using theoretical Green's functions and comparison with the 1994 northridge earthquake

In the evolution of methods for calculating synthetic time histories of ground motion for postulated earthquakes, kinematic source models have dominated to date because of their ease of application. Dynamic models, however, which incorporate a physical relationship between important faulting parameters of stress drop, slip, rupture velocity, and rise time, are becoming more accessible. This article compares a class of kinematic models based on the summation of a fractal distribution of subevent sizes with a dynamic model based on the slip-weakening friction law. Kinematic modeling is done for the frequency band 0.2 to 10.0. Hz, dynamic models are calculated from 0.2 to 2.0. Hz. The strong motion data set for the 1994 Northridge earthquake is used to evaluate and compare the synthetic time histories. Source models are propagated to the far field by convolution with 1D and 3D theoretical Green&rsquo;s functions. In addition, the kinematic model is used to evaluate the importance of propagation path effects: velocity structure, scattering, and nonlinearity. At present, the kinematic model gives a better broadband fit to the Northridge ground motion than the simple slip-weakening dynamic model. In general, the dynamic model overpredicts rise times and produces insufficient shorter-period energy. Within the context of the slip-weakening model, the Northridge ground motion requires a short slip-weakening distance, on the order of 0.15 m or less. A more complex dynamic model including rate weakening or one that allows shorter rise times near the hypocenter may fit the data better.

Bulletin of the Seismological Society of America↗
USGS2005

On pads and filters: Processing strong-motion data

Processing of strong-motion data in many cases can be as straightforward as filtering the acceleration time series and integrating to obtain velocity and displacement. To avoid the introduction of spurious low-frequency noise in quantities derived from the filtered accelerations, however, care must be taken to append zero pads of adequate length to the beginning and end of the segment of recorded data. These padded sections of the filtered acceleration need to be retained when deriving velocities, displacements, Fourier spectra, and response spectra. In addition, these padded and filtered sections should also be included in the time series used in the dynamic analysis of structures and soils to ensure compatibility with the filtered accelerations.

Bulletin of the Seismological Society of America↗
USGS2006

Seismic attenuation structure of the Seattle Basin, Washington State from explosive-source refraction data

We used waveform data from the 1999 SHIPS (Seismic Hazard Investigation of Puget Sound) seismic refraction experiment to constrain the attenuation structure of the Seattle basin, Washington State. We inverted the spectral amplitudes of compressional- and shear-wave arrivals for source spectra, site responses, and one- and two-dimensional Q-1 models at frequencies between 1 and 40 Hz for P waves and 1 and 10 Hz for S waves. We also obtained Q-1 models from t* values calculated from the spectral slopes of P waves between 10 and 40 Hz. One-dimensional inversions show that Qp at the surface is 22 at 1 Hz, 130 at 5 Hz, and 390 at 20 Hz. The corresponding values at 18 km depth are 100, 440, and 1900. Qs at the surface is 16 and 160 at 1 Hz and 8 Hz, respectively, increasing to 80 and 500 at 18 km depth. The t* inversion yields a Qp model that is consistent with the amplitude inversions at 20 and 30 Hz. The basin geometry is clearly resolved in the t* inversion, but the amplitude inversions only imaged the basin structure after removing anomalously high-amplitude shots near Seattle. When these shots are removed, we infer that Q-1 values may be ???30% higher in the center of the basin than the one-dimensional models predict. We infer that seismic attenuation in the Seattle basin will significantly reduce ground motions at frequencies at and above 1 Hz, partially countering amplification effects within the basin.

Bulletin of the Seismological Society of America↗
USGS2005

Erratum: Empirical evidence for acceleration-dependent amplification factors

Incorrect versions of Figures 5 and 6 containing normalization errors were accidentally published by Borcherdt (2002) . They should be replaced with the figures shown here. The text and tabulated regression values published in Borcherdt (2002) are correct and refer to the figures shown here.

Bulletin of the Seismological Society of America↗
USGS2006

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

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

Bulletin of the Seismological Society of America↗
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↗