USGS Science⌕ Search

SEARCH · USGS Science

Results for “The Seismic Record”

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 613 records · Page 34Linked to original sources

The relocation of microearthquakes in the northern Mississippi Embayment

Three-component seismograms, recorded by a small array of digital instruments in the northern Mississippi embayment, consistently show a high-amplitude phase on the vertical component that arrives approximately 0.8 s before the shear wave. On the basis of its timing and apparent velocity, this phase is identified as an S-P conversion from the boundary between the unconsolidated Cenozoic sediments and the underlying Paleozoic rocks. Synthetic models of ground motion in the Mississippi embayment indicate that vertical displacement amplitudes are higher for S-P conversions than for transmitted shear waves at all angles of incidence. The models and their agreement with observations of the three-component seismograms suggest that true shear waves cannot be reliably identified from vertical component seismograms recorded in this area. The travel times of converted phases are used, together with crustal velocity models derived from a recent seismic refraction survey, to relocate approximately 500 microearthquakes recorded by the Central Mississippi Valley regional seismic network. Since the network data are recorded by vertical component seismographs, we assume that the S phases, observed at the stations sited on unconsolidated sediments, are S-P conversions. This assumption significantly reduces the average rms residual and provides well-constrained hypocentral locations for 350 earthquakes. The relocated microearthquakes cluster tightly along previously identified epicentral trends. Focal depths range from 0.5 to 22 km, but events deeper than 14 km are rare, and events shallower than 3 km are confined, almost exclusively, to the area between Ridgely, Tennessee, and New Madrid, Missouri. Focal mechanisms are generally consistent with the results of previous studies. The axial seismicity trend, which extends from Caruthersville, Missouri, to Marked Tree, Arkansas, is characterized by right-lateral strike-slip motion on a northeast trending, nearly vertical fault plane. The events between Ridgely, Tennessee, and New Madrid, Missouri, commonly show strike-slip focal mechanisms, with minor components of both normal and reverse dip slip.

Northern Mississippi Embayment↗

Earthquake source parameters determined by the SAFOD Pilot Hole seismic array

We estimate the source parameters of #3 microearthquakes by jointly analyzing seismograms recorded by the 32-level, 3-component seismic array installed in the SAFOD Pilot Hole. We applied an inversion procedure to estimate spectral parameters for the omega-square model (spectral level and corner frequency) and Q to displacement amplitude spectra. Because we expect spectral parameters and Q to vary slowly with depth in the well, we impose a smoothness constraint on those parameters as a function of depth using a linear first-differenfee operator. This method correctly resolves corner frequency and Q, which leads to a more accurate estimation of source parameters than can be obtained from single sensors. The stress drop of one example of the SAFOD target repeating earthquake falls in the range of typical tectonic earthquakes. Copyright 2004 by the American Geophysical Union.

Geophysical Research Letters↗

Frequency response of the USGS short period telemetered seismic system and its suitability for network studies of local earthquakes

The USGS telemetered seismic system was intended primarily to record small to moderate earthquakes (magnitude 0 to 4) at distances of a few km to several hundred km. Its frequency response is such that the recorded background noise at a moderately quite Coast Range site has a relatively flat 'record' spectrum from about 1/3 Hz to about 20 Hz. With the system magnification set so that the background noise is clearly recorded (about 1 mm peak-to-peak) one can anticipate that any seismic signal that exceeds background noise appreciably in this spectral region will be large enough to be seen on the seismogram. This response represents the highest sensitivity and broadest bandwidth that we were able to attain with a 1-Hz seismometer, a simple amplifier VCO employing very low-power integrated circuits, and an 8-channel constant-bandwidth FM subcarrier multiplex system for use with commercial voice-grade phone lines.

Open-File Report↗

Frequency-dependent effects of rupture for the 2004 Parkfield mainshock, results from UPSAR

The frequency-dependent effects of rupture propagation of the Parkfield, California earthquake (Sept. 28, 2004, M6) to the northwest along the San Andreas fault can be seen in acceleration records at UPSAR (USGS Parkfield Seismic Array) in at least two ways. First, we can see the effects of directivity in the acceleration traces at UPSAR, which is about 11.5 km from the epicenter. Directivity or the seismic equivalent of a Doppler shift has been documented in many cases by comparing short duration, high-amplitude pulses (P or S) in the forward direction with longer duration body waves in the backward direction. In this case we detect a change from a relatively large amplitude, coherent, high-frequency signal at the start of rupture to a low-amplitude, low-coherent, low-frequency signal at about the time the rupture front transfers from the forward azimuth to the back azimuth at about 34-36 s (time is UTC and are the seconds after day 272 and 17 hours and 15 minutes. S arrival is just after 30s) for rays leaving the fault and propagating to UPSAR. The frequency change is obvious in the band about 5 to 30 Hz, which is significantly above the corner frequency of the earthquake (about 0.11Hz). From kinematic source models, the duration of faulting is about 9.2 s and the change in frequency is during faulting as the rupture extends to the northwest. Understanding the systematic change in frequency and amplitude of seismic waves in relation to the propagation of the rupture front is important for predicting strong ground motion. Second, we can filter the acceleration records from the array to determine if the low frequency energy emerges from the same part of the fault as the high frequency signal (e.g. has the same back azimuth and apparent velocity at UPSAR) an important clue to the dynamics of rupture. Analysis of sources of strong motion (characterized by relatively high frequencies) compared to kinematic slip models (relatively low frequency) for the March 11, 2011 Tohoku earthquake as well as Maule (Feb. 27, 2010) and Chi-Chi (Sept. 20, 1999) earthquakes show that high- and low-frequency sources do not have the same locations on the fault. In this paper we filter the accelerograms from UPSAR for the 2004 mainshock in various passbands and then re-compute the cross correlations to determine the vector slowness of the incoming waves. At Parkfield, it appears that for seismic waves with frequencies above 1 Hz there is no discernible frequency-dependent difference in source position (up to 8 Hz) based on estimates of back azimuth and apparent velocity. However at lower frequencies, sources appear to be from shallower depths and trail the high frequencies as the rupture proceeds down the fault. This result is greater than one standard deviation of an estimate of error, based on a new method of estimating error that is a measure of how broad the peak in correlation is and an estimate of the variance of the correlation values. These observations can be understood in terms of a rupture front that is more energetic and coherent near the front of rupture (radiating higher frequencies) and less coherent and less energetic (radiating in a lower frequency band) behind the initial rupture front. This result is a qualitative assessment of changes in azimuth and apparent velocity with frequency and time and does not include corrections to find the source location on the fault.

California↗

Seismic evidence for a slab tear at the Puerto Rico Trench

The fore-arc region of the northeast Caribbean plate north of Puerto Rico and the Virgin Islands has been the site of numerous seismic swarms since at least 1976. A 6 month deployment of five ocean bottom seismographs recorded two such tightly clustered swarms, along with additional events. Joint analyses of the ocean bottom seismographs and land-based seismic data reveal that the swarms are located at depths of 50–150 km. Focal mechanism solutions, found by jointly fitting P wave first-motion polarities and S/P amplitude ratios, indicate that the broadly distributed events outside the swarm generally have strike- and dip-slip mechanisms at depths of 50–100 km, while events at depths of 100–150 km have oblique mechanisms. A stress inversion reveals two distinct stress regimes: The slab segment east of 65°W longitude is dominated by trench-normal tensile stresses at shallower depths (50–100 km) and by trench-parallel tensile stresses at deeper depths (100–150 km), whereas the slab segment west of 65°W longitude has tensile stresses that are consistently trench normal throughout the depth range at which events were observed (50–100 km). The simple stress pattern in the western segment implies relatively straightforward subduction of an unimpeded slab, while the stress pattern observed in the eastern segment, shallow trench-normal tension and deeper trench-normal compression, is consistent with flexure of the slab due to rollback. These results support the hypothesis that the subducting North American plate is tearing at or near these swarms. The 35 year record of seismic swarms at this location and the recent increase in seismicity suggest that the tear is still propagating.

Puerto Rico↗

What can we learn from the Wells, NV earthquake sequence about seismic hazard in the intermountain west?

The February 21, 2008 Wells, NV earthquake (M 6) was felt throughout eastern Nevada, southern Idaho, and western Utah. The town of Wells sustained significant damage to unreinforced masonry buildings. The earthquake occurred in a region of low seismic hazard with little seismicity, low geodetic strain rates, and few mapped faults. The peak horizontal ground acceleration predicted by the USGS National Seismic Hazard Maps is about 0.2 g at 2% probability of exceedance in 50 years, with the contributions coming mostly from the Ruby Mountain fault and background seismicity (M5-7.0). The hazard model predicts that the probability of occurrence of an M>6 event within 50 km of Wells is about 15% in 100 years. Although the earthquake was inside the USArray Transportable Array network, the nearest on-scale recordings of ground motions from the mainshock were too distant to estimate accelerations in town. The University of Nevada Reno, the University of Utah, and the U.S. Geological Survey deployed portable instruments to capture the ground motions from aftershocks of this rare normal-faulting event. Shaking from a M 4.7 aftershock recorded on portable instruments at distances less than 10 km exceeded 0.3 g, and sustained accelerations above 0.1 g lasted for about 5 seconds. For a magnitude 5 earthquake at 10 km distance the NGA equations predict median peak ground accelerations about 0.1 g. Ground motions from normal faulting earthquakes are poorly represented in the ground motion prediction equations. We compare portable and Transportable Array ground-motion recordings with prediction equations. Advanced National Seismic System stations in Utah recorded ground motions 250 km from the mainshock of about 2% g. The maximum ground motion recorded in Salt Lake City was in the center of the basin. We analyze the spatial variability of ground motions (rock vs. soil) and the influence of the Salt Lake Basin in modifying the ground motions. We then compare this data with the September 28, 2004 Parkfield aftershocks to contrast the differences between strike-slip and normal ground motions.

Conference Paper↗

Testing small-aperture array analysis on well-located earthquakes, and application to the location of deep tremor

We have here analyzed local and regional earthquakes using array techniques with the double aim of quantifying the errors associated with the estimation of propagation parameters of seismic signals and testing the suitability of a probabilistic location method for the analysis of nonimpulsive signals. We have applied the zero-lag cross-correlation method to earthquakes recorded by three dense arrays in Puget Sound and Vancouver Island to estimate the slowness and back azimuth of direct P waves and S waves. The results are compared with the slowness and back azimuth computed from the source location obtained by the analysis of data recorded by the Pacific Northwest seismic network (PNSN). This comparison has allowed a quantification of the errors associated with the estimation of slowness and back azimuth obtained through the analysis of array data. The statistical analysis gives ??BP = 10?? and ??BS = 8?? as standard deviations for the back azimuth and ??SP = 0.021 sec/km and ??SS = 0.033 sec /km for the slowness results of the P and S phases, respectively. These values are consistent with the theoretical relationship between slowness and back azimuth and their uncertainties. We have tested a probabilistic source location method on the local earthquakes based on the use of the slowness estimated for two or three arrays without taking into account travel-time information. Then we applied the probabilistic method to the deep, nonvolcanic tremor recorded by the arrays during July 2004. The results of the tremor location using the probabilistic method are in good agreement with those obtained by other techniques. The wide depth range, of between 10 and 70 km, and the source migration with time are evident in our results. The method is useful for locating the source of signals characterized by the absence of pickable seismic phases.

Bulletin of the Seismological Society of America↗

Kinematics of fault slip associated with the July 4-6 2019 Ridgecrest, Californai earthquakes sequence

The 2019 Ridgecrest, California, earthquake sequence produced observable crustal deformation over much of central and southern California, as well as surface rupture over several tens of kilometers. To obtain a detailed picture of the fault slip involved in the 4 July M 6.4 foreshock and 6 July M 7.1 mainshock, we combine strong‐motion seismic waveforms with crustal deformation observations to obtain kinematic and static slip models of both events. We sample the regional seismic wavefield for both the foreshock and mainshock with three‐component records from 31 stations of the California Integrated Seismic Network. The deformation observations include Global Positioning System (GPS), Interferometric Synthetic Aperture Radar (InSAR), and borehole strainmeter recordings of the dynamic strain field. These data collectively constrain the kinematic coseismic slip distributions of the events, with measurements variously observing coseismic slip from one event (e.g., seismic waveforms, kinematic solutions from continuous GPS, and strainmeter time series) or coseismic slip from both events combined (InSAR). We find that the foreshock ruptured two separate faults, one with left‐lateral strike slip on a northeast–southwest‐trending fault and the other with right‐lateral strike slip on an orthogonal fault, with unilateral rupture propagation along both. The mainshock ruptured a series of northwest–southeast‐trending faults with right‐lateral strike slip concentrated in the uppermost 6 km with exceptionally low‐rupture velocity averaging 1.0 – 1.5    km / s ">

California↗

Largest recent impact craters on Mars: Orbital imaging and surface seismic co-investigation

Two 130+ meter diameter impact craters formed on Mars during the later half of 2021. These are the two largest fresh impact craters discovered by the Mars Reconnaissance Orbiter since operations started 16 years ago. The impacts created two of the largest seismic events (magnitudes greater than 4) recorded by InSight during its three year mission. The combination of orbital imagery and seismic ground motion enables the investigation of subsurface and atmospheric energy partitioning of the impact process on a planet with a thin atmosphere and the first direct test of Martian deep-interior seismic models with known event distances.The impact at 35°N excavated blocks of water ice, which is the lowest latitude ice has been directly observed on Mars.

Science↗

Tracking secondary lahar flow paths and characterizing pulses and surges using infrasound array networks at Volcán de Fuego, Guatemala

Lahars are one of the greatest hazards at many volcanoes, including Volcán de Fuego (Guatemala). On 1 December 2018 at 8:00pm local Guatemala time (2:00:00 UTC), an hour-long lahar event was detected at Volcán de Fuego by two permanent seismo-acoustic stations along the Las Lajas channel on the southeast side. To establish the timing, duration, and speed of the lahar, infrasound array records were examined to identify both the source direction(s) and the correlated energy fluctuations at the two stations. Co-located seismic and acoustic signals were also examined, which indicated at least 5 distinct energy pulses within the lahar record. We infer that varying sediment load and/or changes in flow velocity is shown by clear fluctuations in the acoustic and seismic power recorded at one of the stations. This particular event studied with infrasound provides insight into how lahars occur around Volcán de Fuego.

Volcán de Fuego↗

Travel times of P and S from the global digital seismic networks: Implications for the relative variation of P and S velocity in the mantle

We present new data sets of P and S arrival times which have been handpicked from long-period vertical and transverse component recordings of the various global seismic networks. Using events which occurred from 1976 to 1994 results in ∼38,000 globally well-distributed measurements of teleseismic P and ∼41,000 measurements of S . These data are particularly useful for looking at the relative variation of S and P velocities in the lower mantle. We describe both the measurement techniques and the gross characteristics of the data sets. The size of our data sets allows us to exploit the internal consistency of the data to identify outliers using a summary ray analysis. Since the polarity of each arrival is also known, we can construct fault plane solutions and/or compare with polarities predicted by the Harvard centroid moment tensor solutions to further diagnose phase misidentification. This analysis results in ∼5% of the data being identified as outliers. An analysis of variance indicates that the S residual travel times are dominated by the effects of three-dimensional structure but the P data have comparable contributions from noise and source mislocation effects. The summary ray analysis reveals the basic character of lower mantle structure, and there are large-scale patterns in both the S and P data sets that correlate quite well with each other. This analysis suggests that on average, d ln v s / d ln v p is an increasing function of depth in the mantle going from a value of ∼1.7 at the top of the lower mantle to an apparent value of 4 near the base of the mantle. This latter extreme value of R seems to result mainly from data which sample one region in the lowermost mantle under the central Pacific, where large positive S residuals are associated with very small P residuals. Such an anomaly cannot be thermal in origin.

Journal of Geophysical Research B: Solid Earth↗

Forecasting probabilistic seismic shaking for greater Tokyo from 400 years of intensity observations

The long recorded history of earthquakes in Japan affords an opportunity to forecast seismic shaking exclusively from past shaking. We calculate the time-averaged (Poisson) probability of severe shaking by using more than 10,000 intensity observations recorded since AD 1600 in a 350 km-wide box centered on Tokyo. Unlike other hazard-assessment methods, source and site effects are included without modeling, and we do not need to know the size or location of any earthquake nor the location and slip rate of any fault. The two key assumptions are that the slope of the observed frequency-intensity relation at every site is the same, and that the 400-year record is long enough to encompass the full range of seismic behavior. Tests we conduct here suggest that both assumptions are sound. The resulting 30-year probability of IJMA ??? 6 shaking (??? PGA ??? 0.4 g or MMI ??? IX) is 30%-40% in Tokyo, Kawasaki, and Yokohama, and 10% 15% in Chiba and Tsukuba. This result means that there is a 30% chance that 4 million people will be subjected to IJMA ??? 6 shaking during an average 30-year period. We also produce exceedance maps of PGA for building-code regulations, and calculate short-term hazard associated with a hypothetical catastrophe bond. Our results resemble an independent assessment developed from conventional seismic hazard analysis for greater Tokyo. ?? 2007, Earthquake Engineering Research Institute.

Earthquake Spectra↗

Seismicity and eruptive activity at Fuego Volcano, Guatemala: February 1975 -January 1977

We examine seismic and eruptive activity at Fuego Volcano (14??29???N, 90?? 53???W), a 3800-m-high stratovolcano located in the active volcanic arc of Guatemala. Eruptions at Fuego are typically short-lived vulcanian eruptions producing ash falls and ash flows of high-alumina basalt. From February 1975 to December 1976, five weak ash eruptions occurred, accompanied by small earthquake swarms. Between 0 and 140 (average ??? 10) A-type or high-frequency seismic events per day with M > 0.5 were recorded during this period. Estimated thermal energies for each eruption are greater by a factor of 106 than cumulative seismic energies, a larger ratio than that reported for other volcanoes. Over 4000 A-type events were recorded January 3-7, 1977 (cumulative seismic energy ??? 109 joules), yet no eruption occurred. Five 2-hour-long pulses of intense seismicity separated by 6-hour intervals of quiescence accounted for the majority of events. Maximum likelihood estimates of b-values range from 0.7 ?? 0.2 to 2.1 ?? 0.4 with systematically lower values corresponding to the five intense pulses. The low values suggest higher stress conditions. During the 1977 swarm, a tiltmeter located 6 km southeast of Fuego recorded a 14 ?? 3 microradian tilt event (down to SW). This value is too large to represent a simple change in the elastic strain field due to the earthquake swarm. We speculate that the earthquake swarm and tilt are indicative of subsurface magma movement. ?? 1984.

Journal of Volcanology and Geothermal Research↗

Strong SH-to-Love wave scattering off the Southern California Continental Borderland

Seismic scattering is commonly observed and results from wave propagation in heterogeneous medium. Yet, deterministic characterization of scatterers associated with lateral heterogeneities remains challenging. In this study, we analyze broadband waveforms recorded by the Southern California Seismic Network and observe strongly scattered Love waves following the arrival of teleseismic SH wave. These scattered Love waves travel approximately in the same (azimuthal) direction as the incident SH wave at a dominant period of ~10 s but at an apparent velocity of ~3.6 km/s as compared to the ~11 km/s for the SH wave. Back-projection suggests that this strong scattering is associated with pronounced bathymetric relief in the Southern California Continental Borderland, in particular the Patton Escarpment. Finite-difference simulations using a simplified 2-D bathymetric and crustal model are able to predict the arrival times and amplitudes of major scatterers. The modeling suggests a relatively low shear wave velocity in the Continental Borderland.

California↗

Quantifying permanent uplift due to lithosphere-hotspot interaction

Vertical motions that accompany the passage of the lithosphere over a mantle hotspot can shed light on the nature of the hotspot and its effect on the lithosphere. However, quantifying the temporal vertical and spatial extent, is challenging due to the paucity of evidence in the geological record. Here, we utilize dense seismic and well data covering the intersection of the Great Meteor Hotspot (GMH) track with the U.S. Atlantic continental margin to constrain the surface expression of the hotspot passage under the lithosphere. The continuous sedimentary record of the eastern North American margin during its passage over the hotspot allows determination of the timing, magnitude, width and rate of denudation. We find that a ∼300 km wide region was denuded by up to 850 m between ∼97 and 86 Ma, ∼10 m.y. after the passage of the GMH. Stratigraphic relationships suggest a decaying rock uplift rate with time and no subsequent sagging. The broad, long-lasting, and delayed uplift was modeled as a surface manifestation of either sub-lithospheric mantle depletion, permanently eroded base of the continental lithosphere, or intrusions of depleted magma. We consider sub-lithospheric depletion to be the most likely cause, based on seismic imaging results.

Tectonics↗

Use of a hybrid computer in engineering-seismology research

A hybrid computer is an important tool in the seismological research conducted by the U.S. Geological Survey in support of the Energy Research and Development Administration nuclear explosion testing program at the Nevada Test Site and the U.S. Geological Survey Earthquake Hazard Reduction Program. The hybrid computer system, which employs both digital and analog computational techniques, facilitates efficient seismic data processing. Standard data processing operations include: (1) preview of dubbed magnetic tapes of data; (2) correction of data for instrument response; (3) derivation of displacement and acceleration time histories from velocity recordings; (4) extraction of peak-amplitude data; (5) digitization of time histories; (6) rotation of instrumental axes; (7) derivation of response spectra; and (8) derivation of relative transfer functions between recording sites. Catalog of time histories and response spectra of ground motion from nuclear explosions and earthquakes that have been processed by the hybrid computer are used in the Earthquake Hazard Research Program to evaluate the effects of source, propagation path, and site effects on recorded ground motion; to assess seismic risk; to predict system response; and to solve system design problems.

Journal of Research of the U.S. Geological Survey↗

Digital recordings of aftershocks of the 17 October 1989 Loma Prieta, California, earthquake

After the 17 October 1989 Loma Prieta, central California, earthquake (291 00:04:15.25 UTC, M s =7.1), the U. S. Geological Survey recorded aftershocks at sites of seismologic and engineering interest. This report describes a seismic-waveform dataset collected from 17 October 1989 (291 UTC) to 14 March 1990 (073 UTC) with GEOS digital seismographs deployed from USGS offices in Menlo Park, California. It is a summary of field and data-playback information that is intended to facilitate use of the waveform data in seismologic and engineering studies. It includes station locations, instrumentation histories (trigger parameters, sensor parameters, clock corrections, etc.), listings of waveform records, preliminary seismicity listings, and information about data availability. Our principal research goal is to understand the influence of seismic-source and wave-propagation phenomena on damaging ground motions, and to use this knowledge to predict hazards from future earthquakes. Source processes encompass the spatial and temporal variations in the excitation of seismic waves at the fault. Wave-propagation processes include amplification and deamplification as seismic energy radiates away from the fault through rocks and soils of variable seismic properties. Since shallow geology at a site (upper tens to hundreds of meters) can exert a relatively great influence on damaging ground motions, it is useful to separate wave-propagation phenomena into whole-path and site parts. In studying these phenomena, aftershock waveform recordings can play an important role. Compared to main shocks, aftershocks are simple earthquakes and aftershock records tend to be relatively dominated by wave-propagation effects. Geometrically, since aftershocks are distributed over the fault zone, aftershock records should resemble part of the main-shock ground motion contributed by localized rupture. This observation forms the basis of the empirical-Green's-function method wherein main-shock ground motion is modeled by summing aftershock records. The Loma Prieta earthquake provided a rare opportunity that motivated the intensive aftershock-recording effort described in this report. In a seismically noisy urban environment like the San Francisco Bay area, such a large set of high quality seismic data would otherwise be unobtainable.

California↗

Amplification of seismic waves by the Seattle basin, Washington state

Recordings of the 1999 M w 7.6 Chi-Chi (Taiwan) earthquake, two local earthquakes, and five blasts show seismic-wave amplification over a large sedimentary basin in the U.S. Pacific Northwest. For weak ground motions from the Chi-Chi earthquake, the Seattle basin amplified 0.2- to 0.8-Hz waves by factors of 8 to 16 relative to bedrock sites west of the basin. The amplification and peak frequency change during the Chi-Chi coda: the initial S -wave arrivals (0–30 sec) had maximum amplifications of 12 at 0.5–0.8 Hz, whereas later arrivals (35–65 sec) reached amplifications of 16 at 0.3–0.5 Hz. Analysis of local events in the 1.0- to 10.0-Hz frequency range show fourfold amplifications for 1.0-Hz weak ground motion over the Seattle basin. Amplifications decrease as frequencies increase above 1.0 Hz, with frequencies above 7 Hz showing lower amplitudes over the basin than at bedrock sites. Modeling shows that resonance in low-impedance deposits forming the upper 550 m of the basin beneath our profile could cause most of the observed amplification, and the larger amplification at later arrival times suggests surface waves also play a substantial role. These results emphasize the importance of shallow deposits in determining ground motions over large basins.

Washington↗