Search USGS⌕ Search

SEARCH · Search USGS

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 289 records · Page 16Linked to original sources

Is there a basis for preferring characteristic earthquakes over a Gutenberg–Richter distribution in probabilistic earthquake forecasting?

The idea that faults rupture in repeated, characteristic earthquakes is central to most probabilistic earthquake forecasts. The concept is elegant in its simplicity, and if the same event has repeated itself multiple times in the past, we might anticipate the next. In practice however, assembling a fault-segmented characteristic earthquake rupture model can grow into a complex task laden with unquantified uncertainty. We weigh the evidence that supports characteristic earthquakes against a potentially simpler model made from extrapolation of a Gutenberg–Richter magnitude-frequency law to individual fault zones. We find that the Gutenberg–Richter model satisfies key data constraints used for earthquake forecasting equally well as a characteristic model. Therefore, judicious use of instrumental and historical earthquake catalogs enables large-earthquake-rate calculations with quantifiable uncertainty that should get at least equal weighting in probabilistic forecasting.

Bulletin of the Seismological Society of America↗

Differential energy radiation from two earthquakes in Japan with identical M w : The Kyushu 1996 and Tottori 2000 earthquakes

We examine two closely located earthquakes in Japan that had identical moment magnitudes M w but significantly different energy magnitudes M e . We use teleseismic data from the Global Seismograph Network and strong-motion data from the National Research Institute for Earth Science and Disaster Prevention's K-Net to analyze the 19 October 1996 Kyushu earthquake (M w 6.7, M e 6.6) and the 6 October 2000 Tottori earthquake (M w 6.7, M e 7.4). To obtain regional estimates of radiated energy E S we apply a spectral technique to regional (<200 km) waveforms that are dominated by S and Lg waves. For the thrust-fault Kyushu earthquake, we estimate an average regional attenuation Q(f) 230f 0:65 . For the strike-slip Tottori earthquake, the average regional attenuation is Q(f) 180f 0:6 . These attenuation functions are similar to those derived from studies of both California and Japan earthquakes. The regional estimate of E S for the Kyushu earthquake, 3:8 ?? 10 14 J, is significantly smaller than that for the Tottori earthquake, E S 1:3 ?? 10 15 J. These estimates correspond well with the teleseismic estimates of 3:9 ?? 10 14 J and 1:8 ?? 10 15 J, respectively. The apparent stress (T a = ??E s /M 0 with ?? equal to rigidity) for the Kyushu earthquake is 4 times smaller than the apparent stress for the Tottori earthquake. In terms of the fault maturity model, the significantly greater release of energy by the strike-slip Tottori earthquake can be related to strong deformation in an immature intraplate setting. The relatively lower energy release of the thrust-fault Kyushu earthquake can be related to rupture on mature faults at a subduction environment. The consistence between teleseismic and regional estimates of E S is particularly significant as teleseismic data for computing E S are routinely available for all large earthquakes whereas often there are no near-field data.

Bulletin of the Seismological Society of America↗

A constant stress-drop model for producing broadband synthetic seismograms: Comparison with the next generation attenuation relations

Broadband (0.1-20 Hz) synthetic seismograms for finite-fault sources were produced for a model where stress drop is constant with seismic moment to see if they can match the magnitude dependence and distance decay of response spectral amplitudes found in the Next Generation Attenuation (NGA) relations recently developed from strong-motion data of crustal earthquakes in tectonically active regions. The broadband synthetics were constructed for earthquakes of M 5.5, 6.5, and 7.5 by combining deterministic synthetics for plane-layered models at low frequencies with stochastic synthetics at high frequencies. The stochastic portion used a source model where the Brune stress drop of 100 bars is constant with seismic moment. The deterministic synthetics were calculated using an average slip velocity, and hence, dynamic stress drop, on the fault that is uniform with magnitude. One novel aspect of this procedure is that the transition frequency between the deterministic and stochastic portions varied with magnitude, so that the transition frequency is inversely related to the rise time of slip on the fault. The spectral accelerations at 0.2, 1.0, and 3.0 sec periods from the synthetics generally agreed with those from the set of NGA relations for M 5.5-7.5 for distances of 2-100 km. At distances of 100-200 km some of the NGA relations for 0.2 sec spectral acceleration were substantially larger than the values of the synthetics for M 7.5 and M 6.5 earthquakes because these relations do not have a term accounting for Q. At 3 and 5 sec periods, the synthetics for M 7.5 earthquakes generally had larger spectral accelerations than the NGA relations, although there was large scatter in the results from the synthetics. The synthetics showed a sag in response spectra at close-in distances for M 5.5 between 0.3 and 0.7 sec that is not predicted from the NGA relations.

Bulletin of the Seismological Society of America↗

Fragmentary evidence of great-earthquake subsidence during holocene emergence, Valdivia estuary, South Central Chile

A reconnaissance of Holocene stratigraphy beneath fringing marshes of the Valdivia estuary, where an M 9.5 earthquake caused 1-2 m of regional coseismic subsidence in 1960, shows only fragmentary evidence of prehistoric coseismic subsidence. In most of the 150 hand-driven cores that were examined, a distinct unconformity separates 0.5-1.5 m of late Holocene tidal and floodplain mud, peat, and sand from underlying middle Holocene subtidal mud and sand. At the Las Coloradas site, where stratigraphy is best preserved, two A horizons of marsh and meadow soils abruptly overlain by sand and mud probably record coseismic subsidence shortly followed by tsunamis. The amount of subsidence during the earthquakes proved difficult to reconstruct with a diatom transfer function because of differences between modern and fossil diatom assemblages. Maximum 14C ages on macrofossils from the two A horizons at the Las Coloradas site of 1.7-1.3 ka and 2.7-1.7 ka allow correlation of the younger horizon with either of two of six 14C-dated A horizons buried by tsunami sand or post-tsunami tidal sand 200 km to the south at Maull??n, and with a lake-wide mass wasting event in Lago Puyehue, 100 km to the southeast. Tidal records of prehistoric coseismic subsidence at Valdivia are scarce because of a sea-level fall of 3-8 m over the past 6000 years, erosion of marsh and meadow soils during subsidence-induced flooding of the estuary, and largely complete land-level recovery during cycles of coseismic subsidence and postseismic uplift.

Bulletin of the Seismological Society of America↗

Finite-fault analysis of the 2004 Parkfield, California, earthquake using P nl waveforms

We apply a kinematic finite-fault inversion scheme to P nl displacement waveforms recorded at 14 regional stations (&Delta;<2&deg;) to recover the distribution of coseismic slip for the 2004 Parkfield earthquake using both synthetic Green&rsquo;s functions (SGFs) calculated for one-dimensional (1D) crustal-velocity models and empirical Green&rsquo;s functions (EGFs) based on the recordings of a single M w 5.0 aftershock. Slip is modeled on a rectangular fault subdivided into 2&times;2 km subfaults assuming a constant rupture velocity and a 0.5 sec rise time. A passband filter of 0.1&ndash;0.5 Hz is applied to both data and subfault responses prior to waveform inversion. The SGF inversions are performed such that the final seismic moment is consistent with the known magnitude ( M w 6.0) of the earthquake. For these runs, it is difficult to reproduce the entire P nl waveform due to inaccuracies in the assumed crustal structure. Also, the misfit between observed and predicted vertical waveforms is similar in character for different rupture velocities, indicating that neither the rupture velocity nor the exact position of slip sources along the fault can be uniquely identified. The pattern of coseismic slip, however, compares well with independent source models derived using other data types, indicating that the SGF inversion procedure provides a general first-order estimate of the 2004 Parkfield rupture using the vertical P nl records. The best-constrained slip model is obtained using the single-aftershock EGF approach. In this case, the waveforms are very well reproduced for both vertical and horizontal components, suggesting that the method provides a powerful tool for estimating the distribution of coseismic slip using the regional P nl waveforms. The inferred slip model shows a localized patch of high slip (55 cm peak) near the hypocenter and a larger slip area (~50 cm peak) extending between 6 and 20 km to the northwest.

California↗

San Andreas fault geometry at Desert Hot Springs, California, and its effects on earthquake hazards and groundwater

The Mission Creek and Banning faults are two of the principal strands of the San Andreas fault zone in the northern Coachella Valley of southern California. Structural characteristics of the faults affect both regional earthquake hazards and local groundwater resources. We use seismic, gravity, and geological data to characterize the San Andreas fault zone in the vicinity of Desert Hot Springs. Seismic images of the upper 500 m of the Mission Creek fault at Desert Hot Springs show multiple fault strands distributed over a 500 m wide zone, with concentrated faulting within a central 200 m wide area of the fault zone. High-velocity (up to 5000 m=sec) rocks on the northeast side of the fault are juxtaposed against a low-velocity (<2000 m=sec) basin on the southwest side within the upper few hundred meters. Near-surface strands of the Mission Creek fault dip steeply southwestward and northeastward and merge at depth to form a narrower, steeply southwestward-dipping or near-vertical (80?? to 90??) fault zone. The Banning fault, in contrast, dips northeastward (45?? to 70??) toward the Mission Creek fault, and the two faults likely merge into a single San Andreas fault zone at depth, indicating a transtensional fault system. Mainshock hypocenters for two of the historically largest (M >6.0) earthquakes in the area (in 1948 and 1986) occurred at or near the depths (~10 to 12 km) of the merged (San Andreas) fault. Large-magnitude earthquakes that nucleate at or below the merged fault will likely generate strong shaking from guided waves along both fault zones and from amplified seismic waves in the low-velocity basin between the two fault zones. The Mission Creek fault zone is a groundwater barrier with the top of the water table varying by 60 m in depth and the aquifer varying by about 50 m in thickness across a 200 m wide zone of concentrated faulting.

Bulletin of the Seismological Society of America↗

Holocene vertical displacement on the central segments of the Wasatch fault zone, Utah

Compiled per-event vertical-displacement observations from 17 paleoseismic sites along the six central segments of the Wasatch fault zone (WFZ) highlight possible biases and trends in displacement along the fault. The displacement data are consistent with a model of characteristic-type slip, but anomalous and variable displacements indicate that significant natural variability in displacement occurs. When combined into a composite distribution of displacement, 79% of the data fit within a displacement envelope that shows displacement decreasing in a half-ellipse shape from 1.4-3.5 m near the segment centers to 0.6-2.5 m near the ends. Additionally, displacements normalized by the distance from the segment centers to ends decrease from means of 2.0-3.0 m near the segment centers to 1.3-1.9 m near the ends, consistent with characteristic-type slip termination. Although several paleoseismic sites exhibit repeated, similar displacements, the data are sparse and both low-valued (0.5-0.8 m) and high-valued (4.2-4.7 m) outliers suggest complex strain release, possibly resulting from segment interaction and/or noncharacteristic events. Although a global, normal-fault-type surface-rupture-length (SRL) average-displacement regression underpredicts observed WFZ displacements, the largest displacements per segment correspond well with a SRL maximum-displacement regression. This correlation, as well as moderate variability in SRL- and displacement-based moment magnitude, suggests that the anomalous displacements represent the intrinsic variability in characteristic displacement per segment. Thus, minor variations to the characteristic slip model to account for exceptional upper- and lower-bound displacements, e.g., a hybrid characteristic-variable slip model, may be appropriate for the WFZ. Additional paleoseismic data are necessary to address data gaps and biases, to facilitate more robust tests of earthquake-slip models, and to reduce uncertainty in SRL, displacement, and magnitude.

Bulletin of the Seismological Society of America↗

Comment on "Revisiting the 1872 owens valley, California, earthquake" by Susan E. Hough and Kate Hutton

Bakun (2009) argues that the conclusions of Hough and Hutton (2008) are wrong because the study failed to take into account the Sierra Nevada attenuation model of Bakun (2006). In particular, Bakun (2009) argues that propagation effects can explain the relatively high intensities generated by the 1872 Owens Valley earthquake. Using an intensity attenuation model that attempts to account for attenuation through the Sierra Nevada, Bakun (2006) infers the magnitude estimate (M w 7.4–7.5) that is currently accepted by National Earthquake Information Center (NEIC).

Bulletin of the Seismological Society of America↗

Observing rotational and translational ground motions at the HGSD station in Taiwan from 2007 to 2008

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

Bulletin of the Seismological Society of America↗

Recording rotational and translational ground motions of two TAIGER explosions in northeastern Taiwan on 4 March 2008

Two explosions were set off on 4 March 2008 at the N3 explosion site in northeastern Taiwan. The code name for the first shot with 3000 kg explosives is N3P and that for the second shot with 750 kg explosives is N3. To record these two explosions, 8 triaxial rotational sensors, 13 triaxial accelerometers, and 12 six-channel, 24 bit dataloggers with Global Positioning System receivers were deployed to continuously record several hours before and after the explosions. These instruments were installed at about 250 m (1 station), 500 m (11 stations), and 600 m (1 station) from the explosions. The 11 stations form a center array with station spacing of about 5 m. Except for one rotational sensor, onscale records were obtained. Although the N3P shot used four times larger amounts of explosives than those used for the N3 shot, the peak ground translational acceleration and rotational velocity at the 13 station sites from the N3P shot are only about 1.5 times larger than those for the N3 shot. We also observed large variations (by tens of percent) of translational accelerations and rotational velocities at the center array with station spacing of about 5 m. The largest peak rotational velocity was observed for the x component: 2.74 and 1.75 mrad/sec at a distance of 254 m from the N3P and N3 shots, respectively. The main purpose of this article is to document our recordings of rotational and translation motions from two explosions in Taiwan and to release the data online for open access. The translational acceleration data from this experiment have been analyzed by Langston et al. (2009) , and we plan to submit an article with analysis of the rotational velocity data in the future.

Bulletin of the Seismological Society of America↗

Array deployment to observe rotational and translational ground motions along the Meishan Fault, Taiwan: A progress report

Measurements in the near field of earthquakes in Japan ( Takeo, 1998 ) and in Taiwan ( Huang et al. , 2006 ; Liu et al. , 2009 ) indicate that rotational ground motions are many times larger than expected from the classical elasticity theory. The Central Weather Bureau (CWB) deployed four instrumentation sets on the campus of the National Chung-Cheng University, Chia-yi, Taiwan, in order to study in detail earthquakes in the near field for both seismology and earthquake engineering. Both rotational and translational ground motions are being monitored along the active Meishan fault, where a major earthquake occurred in 1906, more than a century ago. The deployed instruments are: (1) a 32-element seismic array in the free field, (2) a 32-element seismic array in a building, (3) a six-channel unit with a broadband seismometer and an accelerometer, and (4) a six-channel unit with an accelerometer and an external rotational velocity sensor.

Bulletin of the Seismological Society of America↗

Seismic-wave strain, rotation, and gradiometry for the 4 March 2008 TAIGER explosions

Acceleration spatial gradients, horizontal strains, and horizontal rotation were computed using strong-motion array data from the 4 March 2008 TAIGER explosions in northeastern Taiwan and used in conjunction with the original three component acceleration data to perform a gradiometric analysis of the strong ground motion wave train. The analysis yields a complex, frequency-dependent view of the nature of seismic-wave propagation over short propagation distances that imply significant lateral velocity changes in structure. Areal strain and rotation about the vertical axis have equal amplitudes and suggest significant wave scattering within the confines of the river valley where the experiment was performed and/or significant departure from an axisymmetric explosion source. Gradiometry shows that the P wave arrives at the array 35° off-azimuth clockwise from the straight-line path and appears to have been refracted from the northern side of the valley. Large, slowly propagating secondary surface waves initially arrive 45° counterclockwise from the straight-line path but later arrivals are seen to propagate in all directions, including back toward the explosion source. A frequency-dependent radiation pattern for the triple-borehole explosion in comparison to the single-borehole explosion explains the differences in the maximum amplitudes between the sources seen in the acceleration data. The use of seismic strain and rotation with standard particle motion wave fields at a single location allows for a direct view of seismic-wave propagation that illuminates the true nature of the seismogram.

Bulletin of the Seismological Society of America↗

Sedimentary basin effects in Seattle, Washington: Ground-motion observations and 3D simulations

Seismograms of local earthquakes recorded in Seattle exhibit surface waves in the Seattle basin and basin-edge focusing of S waves. Spectral ratios of Swaves and later arrivals at 1 Hz for stiff-soil sites in the Seattle basin show a dependence on the direction to the earthquake, with earthquakes to the south and southwest producing higher average amplification. Earthquakes to the southwest typically produce larger basin surface waves relative to S waves than earthquakes to the north and northwest, probably because of the velocity contrast across the Seattle fault along the southern margin of the Seattle basin. S to P conversions are observed for some events and are likely converted at the bottom of the Seattle basin. We model five earthquakes, including the M 6.8 Nisqually earthquake, using 3D finite-difference simulations accurate up to 1 Hz. The simulations reproduce the observed dependence of amplification on the direction to the earthquake. The simulations generally match the timing and character of basin surface waves observed for many events. The 3D simulation for the Nisqually earth-quake produces focusing of S waves along the southern margin of the Seattle basin near the area in west Seattle that experienced increased chimney damage from the earthquake, similar to the results of the higher-frequency 2D simulation reported by Stephenson et al. (2006). Waveforms from the 3D simulations show reasonable agreement with the data at low frequencies (0.2-0.4 Hz) for the Nisqually earthquake and an M 4.8 deep earthquake west of Seattle.

Bulletin of the Seismological Society of America↗

Review: Progress in rotational ground-motion observations from explosions and local earthquakes in Taiwan

Rotational motions generated by large earthquakes in the far field have been successfully measured, and observations agree well with the classical elasticity theory. However, recent rotational measurements in the near field of earthquakes in Japan and in Taiwan indicate that rotational ground motions are 10 to 100 times larger than expected from the classical elasticity theory. The near-field strong-motion records of the 1999 M w 7.6 Chi-Chi, Taiwan, earthquake suggest that the ground motions along the 100 km rupture are complex. Some rather arbitrary baseline corrections are necessary in order to obtain reasonable displacement values from double integration of the acceleration data. Because rotational motions can contaminate acceleration observations due to the induced perturbation of the Earth’s gravitational field, we started a modest program to observe rotational ground motions in Taiwan. Three papers have reported the rotational observations in Taiwan: (1) at the HGSD station ( Liu et al. , 2009 ), (2) at the N3 site from two TAiwan Integrated GEodynamics Research (TAIGER) explosions ( Lin et al. , 2009 ), and (3) at the Taiwan campus of the National Chung-Cheng University (NCCU) ( Wu et al. , 2009 ). In addition, Langston et al. (2009) reported the results of analyzing the TAIGER explosion data. As noted by several authors before, we found a linear relationship between peak rotational rate (PRR in mrad/sec) and peak ground acceleration (PGA in m/sec 2 ) from local earthquakes in Taiwan, PRR=0.002+1.301 PGA, with a correlation coefficient of 0.988.

Bulletin of the Seismological Society of America↗

On baseline corrections and uncertainty in response spectrafor baseline variations commonly encountered in digital accelerograph records

Most digital accelerograph recordings are plagued by long-period drifts, best seen in the velocity and displacement time series obtained from integration of the acceleration time series. These drifts often result in velocity values that are nonzero near the end of the record. This is clearly unphysical and can lead to inaccurate estimates of peak ground displacement and long-period spectral response. The source of the long-period noise seems to be variations in the acceleration baseline in many cases. These variations could be due to true ground motion (tilting and rotation, as well as local permanent ground deformation), instrumental effects, or analog-to-digital conversion. Very often the trends in velocity are well approximated by a linear trend after the strong shaking subsides. The linearity of the trend in velocity implies that no variations in the baseline could have occurred after the onset of linearity in the velocity time series. This observation, combined with the lack of any trends in the pre-event motion, allows us to compute the time interval in which any baseline variations could occur. We then use several models of the variations in a Monte Carlo procedure to derive a suite of baseline-corrected accelerations for each noise model using records from the 1999 Chi-Chi earthquake and several earthquakes in Turkey. Comparisons of the mean values of the peak ground displacements, spectral displacements, and residual displacements computed from these corrected accelerations for the different noise models can be used as a guide to the accuracy of the baseline corrections. For many of the records considered here the mean values are similar for each noise model, giving confidence in the estimation of the mean values. The dispersion of the ground-motion measures increases with period and is noise-model dependent. The dispersion of inelastic spectra is greater than the elastic spectra at short periods but approaches that of the elastic spectra at longer periods. The elastic spectra from the most basic processing, in which only the pre-event mean is removed from the acceleration time series, do not diverge from the baseline-corrected spectra until periods of 10-20 sec or more for the records studied here, implying that for many engineering purposes elastic spectra can be used from records with no baseline correction or filtering.

Bulletin of the Seismological Society of America↗

Impediments to predicting site response: Seismic property estimation and modeling simplifications

We compare estimates of the empirical transfer function (ETF) to the plane SH-wave theoretical transfer function (TTF) within a laterally constant medium for invasive and noninvasive estimates of the seismic shear-wave slownesses at 13 Kiban-Kyoshin network stations throughout Japan. The difference between the ETF and either of the TTFs is substantially larger than the difference between the two TTFs computed from different estimates of the seismic properties. We show that the plane SH-wave TTF through a laterally homogeneous medium at vertical incidence inadequately models observed amplifications at most sites for both slowness estimates, obtained via downhole measurements and the spectral analysis of surface waves. Strategies to improve the predictions can be separated into two broad categories: improving the measurement of soil properties and improving the theory that maps the 1D soil profile onto spectral amplification. Using an example site where the 1D plane SH-wave formulation poorly predicts the ETF, we find a more satisfactory fit to the ETF by modeling the full wavefield and incorporating spatially correlated variability of the seismic properties. We conclude that our ability to model the observed site response transfer function is limited largely by the assumptions of the theoretical formulation rather than the uncertainty of the soil property estimates.

Bulletin of the Seismological Society of America↗

Scenario liquefaction hazard maps of Santa Clara Valley, Northern California

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

Bulletin of the Seismological Society of America↗

Along-Arc and Back-Arc Attenuation, Site Response, and Source Spectrum for the Intermediate-Depth 8 January 2006 M 6.7 Kythera, Greece, Earthquake

An M 6.7 intermediate-depth (66 km), in-slab earthquake occurring near the island of Kythera in Greece on 8 January 2006 was well recorded on networks of stations equipped with acceleration sensors and with broadband velocity sensors. All data were recorded digitally using recording instruments with resolutions ranging from almost 11 to 24 bits. We use data from these networks to study the distance dependence of the horizontal-component Fourier acceleration spectra (FAS) and horizontal-component pseudoabsolute response spectral acceleration (PSA). For purposes of simulating motions in the future, we parameterize the distance decay using several forms of the geometrical-spreading function, for each of which we derive Q as a function of frequency. By extrapolating the distance decay back to 1 km, we obtain a reference spectrum that can be used in future simulations. This spectrum requires a more complicated spectral shape than the classic single-corner-frequency model; in particular, there appears to be an enhancement of motion around 0.2-0.3 Hz that may be due to the radiation of a 3-5 sec pulse from the source. We infer a ?? 0 value of about 0.055 sec for rock stations and a stress parameter in the range of 400-600 bars. We also find distinctive differences in the site response of stations on soft soil and soil; both the FAS and the 5% damped PSA amplifications have similar peak amplitudes (about 2 and 4 for soil and soft-soil sites, respectively, relative to the rock sites) at similar frequencies (between about 0.4 and 2.0 Hz, with the soft-soil amplifications peaking at somewhat lower frequencies than the soil amplifications). One of the most distinctive features of the data is the clear difference in the motions for along-arc and back-arc stations, with the former being significantly higher than the latter over a broad range of frequencies at distances beyond about 250 km. The motions from the Kythera earthquake are roughly comparable to those from intermediate-depth earthquakes elsewhere, but they appear to be significantly higher than those from recordings of shallow earthquakes in Greece of comparable magnitude and hypocentral distance.

Bulletin of the Seismological Society of America↗