USGS Science⌕ Search

SEARCH · USGS Science

Results for “Earthquake Spectra”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 343 records · Page 19Linked to original sources

High-pass corner frequency selection and review tool for use in ground-motion processing

Raw seismological waveform data contain noise from the instrument’s surroundings and the instrument itself that can dominate recordings at low and high frequencies. To use these data in ground‐motion modeling, the effects of noise on the signals must be reduced and the signals’ usable frequency range identified. We present automated procedures to efficiently reduce low‐frequency noise that are implemented in the software package gmprocess. These procedures check for, and as needed remove, low‐frequency artifacts in the displacement record using polynomial fits, which can be used in combination with existing signal‐to‐noise ratio (SNR)‐based corner‐frequency selection procedures. The automated selections are then efficiently verified and refined using a graphical user interface (GUI) that plots relevant ground‐motion time series and spectra and tracks modifications to signal processing parameters. We demonstrate these procedures using recordings from the 2020 M 5.1 Sparta, North Carolina, and the 2013 M 4.7 southern Ontario earthquakes. Data processed with the SNR‐only and polynomial criteria for these events contain displacement artifacts in 37% and 23% of processed traces, respectively. Records with remaining artifacts are corrected manually using the GUI. These processing steps illustrate the workflow for efficient data processing with quality control.

Seismological Research Letters↗

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↗

Complex earthquake rupture and local tsunamis

In contrast to far-field tsunami amplitudes that are fairly well predicted by the seismic moment of subduction zone earthquakes, there exists significant variation in the scaling of local tsunami amplitude with respect to seismic moment. From a global catalog of tsunami runup observations this variability is greatest for the most frequently occurring tsunamigenic subduction zone earthquakes in the magnitude range of 7 < M w < 8.5. Variability in local tsunami runup scaling can be ascribed to tsunami source parameters that are independent of seismic moment: variations in the water depth in the source region, the combination of higher slip and lower shear modulus at shallow depth, and rupture complexity in the form of heterogeneous slip distribution patterns. The focus of this study is on the effect that rupture complexity has on the local tsunami wave field. A wide range of slip distribution patterns are generated using a stochastic, self-affine source model that is consistent with the falloff of far-field seismic displacement spectra at high frequencies. The synthetic slip distributions generated by the stochastic source model are discretized and the vertical displacement fields from point source elastic dislocation expressions are superimposed to compute the coseismic vertical displacement field. For shallow subduction zone earthquakes it is demonstrated that self-affine irregularities of the slip distribution result in significant variations in local tsunami amplitude. The effects of rupture complexity are less pronounced for earthquakes at greater depth or along faults with steep dip angles. For a test region along the Pacific coast of central Mexico, peak nearshore tsunami amplitude is calculated for a large number ( N = 100) of synthetic slip distribution patterns, all with identical seismic moment ( M w = 8.1). Analysis of the results indicates that for earthquakes of a fixed location, geometry, and seismic moment, peak nearshore tsunami amplitude can vary by a factor of 3 or more. These results indicate that there is substantially more variation in the local tsunami wave field derived from the inherent complexity subduction zone earthquakes than predicted by a simple elastic dislocation model. Probabilistic methods that take into account variability in earthquake rupture processes are likely to yield more accurate assessments of tsunami hazards.

Journal of Geophysical Research B: Solid Earth↗

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↗

3D near-surface soil response from H/V ambient-noise ratios

The applicability of the horizontal-to-vertical (H/V) ambient-noise spectral ratio for characterizing earthquake site effects caused by nearsurface topography and velocity structures was evaluated at sites underlain by thick (i.e. >100 m) sediment deposits near the southern-end of the New Madrid seismic zone in the central United States. Three-component ambient-noise and velocity models derived from seismic (shearwave) refraction/reflection surveys showed that a relatively horizontal, sharp shear-wave velocity interface in the soil column resulted in an H/V spectral ratio with a single well-defined peak. Observations at sites with more than one sharp shear-wave velocity contrast and horizontally arranged soil layers resulted in at least two well-defined H/V spectral ratio peaks. Furthermore, at sites where there were sharp shear-wave velocity contrasts in nonhorizontal, near-surface soil layers, the H/V spectra exhibited a broad-bandwidth, relatively low amplitude signal instead of a single well-defined peak.

Soil Dynamics and Earthquake Engineering↗

Site response for urban Los Angeles using aftershocks of the Northridge earthquake

Ground-motion records from aftershocks of the 1994 Northridge earthquake are used to estimate site response in the urban Los Angeles area. Over 1300 shear-wave records from 61 sources and 90 sites are used in a linear inversion for source and site-response spectra. The methodology makes no assumptions about the shape of the source spectrum. To obtain a stable unique inverse, a Q model and geometrical spreading factor are assumed. In addition, the site response at a hard-rock site is constrained to be approximately 1.0 with a kappa of 0.02. The site-response spectra compare favorably with the results of previous and on-going investigations in Los Angeles. A couple of first-order effects are lower site response in the surrounding mountains, dominated by Mesozoic and Tertiary rocks, and higher values in the San Fernando and Los Angeles Basins, containing surficial Pleistocene and Holocene alluvial deposits. Results show good correlation of high site-response spectral values with localized areas of severe damage (Interstate 10 collapse, Sherman Oaks, Northridge, Interstate 5/14 collapse). However, widespread trends in site response across the sedimentary basins are not obvious. The data suggest that site responses are lower near the southern margin of the San Fernando Valley for sources to the north, due to north-dipping sedimentary structures. But the general pattern of site response is characterized by high variability on length scales less than a kilometer. Variations of a factor of 2 in site response are observed over the length scale of 200 m and for the same surficial geologic unit. For some of the alluvial basin sites, surface-wave generation is a significant contributor to elevated site response at lower frequencies, below 2 Hz. The total damage pattern for the Northridge earthquake is influenced by strong source directivity to the north and strong local site effects. The correlation of weak-motion site-response estimates with areas of significant damage demonstrates the value of these field measurements in future urban planning and in the reduction of seismic risk in urban areas.

California↗

Building damage survey and microtremor measurements for the source region of the 2015 Gorkha, Nepal, earthquake

We performed a damage survey of buildings and carried out microtremor observations in the source region of the 2015 Gorkha earthquake. Our survey area spans the Kathmandu valley and areas to the east and north of the valley. Damage of buildings in the Kathmandu valley was localized, and the percentage of the totally collapsed buildings was less than 5 %. East of the Kathmandu valley, especially in Sindhupalchok district, damage of buildings was more severe. In the center of Chautara and Bahrabise, towns in Sindhupalchok district, the percentage of the totally collapsed houses exceeded 40 %. North of the Kathmandu valley, the damage was moderate, and 20–30 % of the buildings were totally collapsed in Dhunche. Based on the past studies and our microtremor observations near the strong motion station, the H/V spectrum in Kathmandu has a peak at around 0.3 Hz, which reflects the velocity contrast of the deep sedimentary basin. The H/V spectra in Bahrabise, Chautara, and Dhunche do not show clear peaks, which suggests that the sites have stiff soil conditions. Therefore, the more severe damage outside the Kathmandu valley compared with the relatively light damage levels in the valley is probably due to the source characteristics of the earthquake and/or the seismic performance of buildings, rather than the local site conditions.

Earth, Planets and Space↗

VS30 – A site-characterization parameter for use in building Codes, simplified earthquake resistant design, GMPEs, and ShakeMaps

VS30, defined as the average seismic shear-wave velocity from the surface to a depth of 30 meters, has found wide-spread use as a parameter to characterize site response for simplified earthquake resistant design as implemented in building codes worldwide. VS30 , as initially introduced by the author for the US 1994 NEHRP Building Code, provides unambiguous definitions of site classes and site coefficients for site-dependent response spectra based on correlations derived from extensive borehole logging and comparative ground-motion measurement programs in California. Subsequent use of VS30 for development of strong ground motion prediction equations (GMPEs) and measurement of extensive sets of VS borehole data have confirmed the previous empirical correlations and established correlations of SVS30 with VSZ at other depths. These correlations provide closed form expressions to predict S30 V at a large number of additional sites and further justify S30 V as a parameter to characterize site response for simplified building codes, GMPEs, ShakeMap, and seismic hazard mapping.

Conference Paper↗

Preliminary report on aftershock sequence for earthquake of January 31, 1986, near Painesville, Ohio (time period: 2/1/86-2/10/86)

A ten-station array of broad-band digital instrumentation (GEOS) was deployed by the U. S. Geological Survey with partial support provided by Electric Power Research Institute to record the aftershock sequence of the moderate (mb ~ 4.9) earthquake that occurred on January 31, 1986 (16:46:43 UTC) near Painesville, Ohio. The occurrence of the event has raised questions concerning possible contributory factors to the occurrence of the event and questions concerning the character of earthquake-induced high-frequency ground motions in the area. To aid in the timely resolution of the implications of some of these questions, this preliminary report provides copies of the ground motion time-histories and corresponding spectra for the six identified aftershocks and two events, thought to be quarry blasts, recorded as of February 10, 1986. Recording station locations and epicenter locations based on two preliminary estimates of local seismic velocity structure are provided.

Ohio↗

Strong-motion data from the two Pingtung, Taiwan, earthquakes of 26 December 2006

1016 strong-motion records at 527 free-field stations and 131 records at 42 strong-motion arrays at buildings and bridges were obtained for the Pingtung earthquake doublet from the Taiwan Central Weather Bureau's dense, digital strong-motion network. We carried out standard processing of these strong-motion records at free-field stations. A data set, including the originally recorded files, processed data files, and supporting software and information, is archived online http:// tecdc.earth.sinica.edu.tw/data/EQ2006Pingtung/. We have not yet completed the processing of the strong-motion array data at buildings and bridges. However, some preliminary results and the strong-motion array data recorded at the second nearest instrumented building to the Pingtung earthquake doublet are shown. This paper is intended to document our data processing procedures and the online archived data files, so that researchers can efficiently use the data. We also include two preliminary analyses: (1) a comparison of ground motions recorded by multiple accelerographs at a common site, the TAP 117 station in Taipei, and (2) attenuation of the horizontal ground motions (peak acceleration and response spectra at periods of 0.2, 1.0, and 3.0 s) with respect to distance. Our comparison study of multiple recordings at TAP 117 indicates that waveform coherence among 20- and 24-bit accelerograph records is much higher as compared to records from 16-bit or 12-bit accelerographs, suggesting that the former are of better quality. For the 20- and 24-bit accelerographs, waveform coherence is nearly 1 over the frequency range 1 to 8 Hz for all components, and is greater than about 0.9 from 8 to 20 Hz for the horizontal component, but only from 8 to 12 Hz for the vertical component. Plots of pseudo-acceleration response spectra (PSA) as a function of distance, however, show no clear indication for a difference related to the performance level of the accelerographs. The ground-motions of the first event (Mw = 7.0) are comparable, or even somewhat lower, than those from the smaller second event (Mw = 6.9), consistent with the relative difference of the local magnitudes (ML = 6.96 and 6.99 for the first and second events, respectively). The ground motions from the first event are generally lower than those predicted from equations based on other in-slab subduction earthquakes, whereas the ground motions from the second event are closer to the predictions. Ground-motions for soil sites are generally larger than those from rock sites.

Terrestrial, Atmospheric and Oceanic Sciences↗

Prediction of regional broadband strong ground motions using a teleseismic source model of the 18 April 2014 Mw 7.3 Papanoa, Mexico, earthquake

To estimate predicted ground motion from a teleseismic slip model, we use a low‐ and high‐frequency hybrid method to simulate the regional, strong ground motions observed following the 18 April 2014 moment magnitude ( ⁠M w ⁠ ) 7.3 Papanoa, Mexico, earthquake. To generate the regional ground motion at low frequencies (<1 Hz), a teleseismically derived, finite‐fault, kinematic model is used to define the earthquake source, taking into account slip‐model variations identified with a parameter sampling approach that considers possible errors in the fault geometry, the hypocenter depth, and the rupture velocity. A 3D crustal model is used to calculate the low‐frequency ground motions using a finite‐element calculation that includes topography and considers variations in the source model to estimate the uncertainty in the calculations. High frequencies (>1 Hz) are added using a 1D full‐wave propagation code that estimates uncertainties by considering multiple random distributions of slip with different spatial correlation lengths. The synthetic, broadband (0.05–10.0 Hz) ground motions are obtained by combining the low‐ and high‐frequency portions match filtered at 1 Hz. These synthetic ground motions are compared with the regional observations using velocity records, peak ground acceleration, and medians of the orientation‐independent response spectra of the horizontal components (RotD50) calculated at periods of 0.2, 0.3, 0.5, 1.0, 2.0, 3.0, 5.0, 7.5, and 10.0 s. The results indicate that ground motions estimated at these periods using our hybrid approach based primarily on a teleseismically derived source model are comparable to the values observed for the 2014 Papanoa earthquake at regional distances. The approach could be used to estimate strong‐motion spectral levels expected for regions with limited local and regional recordings and could also fill in magnitude or distance gaps in ground‐motion prediction relations utilized in the assessment of seismic hazard.

Papanoa↗

The global aftershock zone

The aftershock zone of each large (M &ge; 7) earthquake extends throughout the shallows of planet Earth. Most aftershocks cluster near the mainshock rupture, but earthquakes send out shivers in the form of seismic waves, and these temporary distortions are large enough to trigger other earthquakes at global range. The aftershocks that happen at great distance from their mainshock are often superposed onto already seismically active regions, making them difficult to detect and understand. From a hazard perspective we are concerned that this dynamic process might encourage other high magnitude earthquakes, and wonder if a global alarm state is warranted after every large mainshock. From an earthquake process perspective we are curious about the physics of earthquake triggering across the magnitude spectrum. In this review we build upon past studies that examined the combined global response to mainshocks. Such compilations demonstrate significant rate increases during, and immediately after (~ 45 min) M > 7.0 mainshocks in all tectonic settings and ranges. However, it is difficult to find strong evidence for M > 5 rate increases during the passage of surface waves in combined global catalogs. On the other hand, recently published studies of individual large mainshocks associate M > 5 triggering at global range that is delayed by hours to days after surface wave arrivals. The longer the delay between mainshock and global aftershock, the more difficult it is to establish causation. To address these questions, we review the response to 260 M &ge; 7.0 shallow (Z &le; 50 km) mainshocks in 21 global regions with local seismograph networks. In this way we can examine the detailed temporal and spatial response, or lack thereof, during passing seismic waves, and over the 24 h period after their passing. We see an array of responses that can involve immediate and widespread seismicity outbreaks, delayed and localized earthquake clusters, to no response at all. About 50% of the catalogs that we studied showed possible (localized delayed) remote triggering, and ~ 20% showed probable (instantaneous broadly distributed) remote triggering. However, in any given region, at most only about 2&ndash;3% of global mainshocks caused significant local earthquake rate increases. These rate increases are mostly composed of small magnitude events, and we do not find significant evidence of dynamically triggered M > 5 earthquakes. If we assume that the few observed M > 5 events are triggered, we find that they are not directly associated with surface wave passage, with first incidences being 9&ndash;10 h later. We note that mainshock magnitude, relative proximity, amplitude spectra, peak ground motion, and mainshock focal mechanisms are not reliable determining factors as to whether a mainshock will cause remote triggering. By elimination, azimuth, and polarization of surface waves with respect to receiver faults may be more important factors.

Tectonophysics↗

Observations of volcanic tremor during January-February 2005 eruption of Mt. Veniaminof, Alaska

Mt. Veniaminof, Alaska Peninsula, is a stratovolcano with a summit ice-filled caldera containing a small intracaldera cone and active vent. From January 2 to February 21, 2005, Mt. Veniaminof erupted. The eruption was characterized by numerous small ash emissions (VEI 0 to 1) and accompanied by low-frequency earthquake activity and volcanic tremor. We have performed spectral analyses of the seismic signals in order to characterize them and to constrain their source. Continuous tremor has durations of minutes to hours with dominant energy in the band 0.5&ndash; 4.0 Hz, and spectra characterized by narrow peaks either irregularly (non-harmonic tremor) or regularly spaced (harmonic tremor). The spectra of non-harmonic tremor resemble those of low-frequency events recorded simultaneously with surface ash explosions, suggesting that the source mechanisms might be similar or related. We propose that non-harmonic tremor at Mt. Veniaminof results from the coalescence of gas bubbles while low-frequency events are related to the disruption of large gas pockets within the conduit. Harmonic tremor, characterized by regular and quasisinusoidal waveforms, has duration of hours. Spectra containing up to five harmonics suggest the presence of a resonating source volume that vibrates in a longitudinal acoustic mode. An interesting feature of harmonic tremor is that frequency is observed to change over time; spectral lines move towards higher or lower values while the harmonic nature of the spectra is maintained. Factors controlling the variable characteristics of harmonic tremor include changes in acoustic velocity at the source and variations of the effective size of the resonator.

Alaska↗

Earthquake site response in Santa Cruz, California

Aftershocks of the 1989 Loma Prieta, California, earthquake are used to estimate site response in a 12-km 2 area centered on downtown Santa Cruz. A total of 258 S -wave records from 36 aftershocks recorded at 33 sites are used in a linear inversion for site-response spectra. The inversion scheme takes advantage of the redundancy of the large data set for which several aftershocks are recorded at each site. The scheme decomposes the observed spectra into source, path, and site terms. The path term is specified before the inversion. The undetermined degree of freedom in the decomposition into source and site spectra is removed by specifying the site-response factor to be approximately 1.0 at two sites on crystalline bedrock. The S -wave site responses correlate well with the surficial geology and observed damage pattern of the mainshock. The site-response spectra of the floodplain sites, which include the heavily damaged downtown area, exhibit significant peaks. The largest peaks are between 1 and 4 Hz. Five floodplain sites have amplification factors of 10 or greater. Most of the floodplain site-response spectra also have a smaller secondary peak between 6 and 8 Hz. Residential areas built on marine terraces above the floodplain experienced much less severe damage. Site-response spectra for these areas also have their largest peaks between 1 and 4 Hz, but the amplification is generally below 6. Several of these sites also have a secondary peak between 6 and 8 Hz. The response peaks seen at nearly all sites between 1 and 4 Hz are probably caused by the natural resonance of the sedimentary rock column. The higher amplifications at floodplain sites may be caused by surface waves generated at the basin margins. The secondary peak between 6 and 8 Hz at many sites may be a harmonic of the 1- to 4-Hz peaks.

California↗

Comparisons of ground motions from the 1999 Chi-Chi, earthquake with empirical predictions largely based on data from California

This article has the modest goal of comparing the ground motions recorded during the 1999 Chi-Chi, Taiwan, mainshock with predictions from four empirical-based equations commonly used for western North America; these empirical predictions are largely based on data from California. Comparisons are made for peak acceleration and 5%-damped response spectra at periods between 0.1 and 4 sec. The general finding is that the Chi-Chi ground motions are smaller than those predicted from the empirically based equations for periods less than about 1 sec by factors averaging about 0.4 but as small as 0.26 (depending on period, on which equation is used, and on whether the sites are assumed to be rock or soil). There is a trend for the observed motions to approach or even exceed the predicted motions for longer periods. Motions at similar distances (30-60 km) to the east and to the west of the fault differ dramatically at periods between about 2 and 20 sec: Long-duration wave trains are present on the motions to the west, and when normalized to similar amplitudes at short periods, the response spectra of the motions at the western stations are as much as five times larger than those of motions from eastern stations. The explanation for the difference is probably related to site and propagation effects; the western stations are on the Coastal Plain, whereas the eastern stations are at the foot of young and steep mountains, either in the relatively narrow Longitudinal Valley or along the eastern coast-the sediments underlying the eastern stations are probably shallower and have higher velocity than those under the western stations.

Bulletin of the Seismological Society of America↗

Digital seismograms of aftershocks of the Imperial Valley, California, earthquake of October 15, 1979

This report presents the digital seismograms collected for the aftershocksof the Imperial Valley, California earthquake (M L = 6.6) that occurred onL October 15, 1979. A short summary follows that describes the field procedures, instrumentation, and routine computer analysis that was used to produce the ground motion time histories. Boore and Fletcher (1981) may be referenced for selected fault plane solutions, spectra and moments. Requests for digital time histories on magnetic tape may be made to the Branch of Ground Motion and Faulting, Menlo Park, California.

California↗

A partially nonergodic ground-motion model for Fourier amplitude spectra for the San Francisco Bay area, California, USA

We develop a partially nonergodic ground-motion model (GMM) for Fourier amplitude spectra for the San Francisco Bay Area, California, USA, using the Bayless and Abrahamson (2019) GMM as a reference ergodic GMM and developing location-dependent adjustments to the predicted median and variance. We compile regional ground-motion data from moment magnitude (𝑀 w ) >3 earthquakes occurring during 2000–2022 for which magnitude information is available in the U.S. Geological Survey Comprehensive Catalog (Guy et al., 2015). The data set predominantly consists of records from 𝑀 w 3.5–4.5 earthquakes but includes three well-recorded 𝑀 w > 5 events. Ground-motion residuals are evaluated using the time-averaged shear-wave velocity in the top 30 m (𝑉 S30 ) from the California-specific map of Thompson et al. (2018) and basin-depth site parameters from the seismic velocity model of Aagaard and Hirakawa (2021). The 𝑉 S30 dependence and basin-depth scaling of the reference ergodic GMM of Bayless and Abrahamson (2019) are evaluated and modified with the updated data set. We compute maps of site adjustments using a varying-coefficient model that considers the spatial correlation structure and uncertainties at each observation location. The spatial covariance model is developed using ground-motion residuals that are standardized by the uncertainty model, which allows for consideration of the aleatory variability in developing the site adjustments. The covariance model is fit considering the means and standard deviations of the site terms at all locations. The use of partially nonergodic median adjustments results in modified variance components of the within-event variability. Due to the low number of large-magnitude earthquakes that control seismic hazard in the data set, we do not modify between-event variance; however, we present adjustments to site-to-site variability for use in partially nonergodic hazard assessments.

California↗

Fault-zone guided waves from explosions in the San Andreas fault at Parkfield and Cienega Valley, California

Fault-zone guided waves were successfully excited by near-surface explosions in the San Andreas fault zone both at Parkfield and Cienega Valley, central California. The guided waves were observed on linear, three-component seismic arrays deployed across the fault trace. These waves were not excited by explosions located outside the fault zone. The amplitude spectra of guided waves show a maximum peak at 2 Hz at Parkfield and 3 Hz at Cienega Valley. The guided wave amplitude decays sharply with observation distance from the fault trace. The explosion-excited fault-zone guided waves are similar to those generated by earthquakes at Parkfield but have lower frequencies and travel more slowly. These observations suggest that the fault-zone wave guide has lower seismic velocities as it approaches the surface at Parkfield. We have modeled the waveforms as S waves trapped in a low-velocity wave guide sandwiched between high-velocity wall rocks, resulting in Love-type fault-zone guided waves. While the results are nonunique, the Parkfield data are adequately fit by a shallow wave guide 170 m wide with an S velocity 0.85 km/sec and an apparent Q ∼ 30 to 40. At Cienega Valley, the fault-zone wave guide appears to be about 120 m wide with an S velocity 0.7 km/sec and a Q ∼ 30.

California↗