A site effect study in Acapulco, Guerrero, Mexico: Comparison of results from strong-motion and microtremor data

1992 ◽  
Vol 82 (2) ◽  
pp. 642-659 ◽  
Author(s):  
Carlos Gutierrez ◽  
Shri Krishna Singh

Abstract The city of Acapulco is located near or above the mature seismic gap of Guerrero along the Mexican subduction zone. With the purpose of studying the character of strong ground motion on soft sites, four digital accelerographs have been installed in the city on such sites. These instruments have been in operation since 1988. Two additional instruments, part of the Guerrero Accelerograph Array, are located on hard sites in the area. One of these, VNTA, has been in operation since 1985 and the other, ACAN, since 1989. These stations have recorded several earthquakes. We use data from eight events (4.2 ≤ M ≤ 6.9) to study spectral amplification of seismic waves at the soft sites with respect to VNTA. The S waves are amplified by a factor of 6 to 25 at the soft sites in a fairly broad range of frequencies; both the amplification and the frequency band over which it occurs depend upon the site. Although the largest earthquake in our data set (M = 6.9) gave rise to a peak horizontal acceleration exceeding 0.3 g at one of the soft sites, no clear evidence of nonlinear behavior of the subsoil is found. Spectral amplifications of S-wave coda are very similar to those of S waves. We also measured microtremors at the strong-motion sites. The microtremor spectra were interpreted, using reasonable assumptions, to test the feasibility of this technique in reproducing the spectral amplifications observed during earthquakes. Our results show that only a rough estimate of site response can be obtained from this technique, at least in Acapulco; caution is warranted in its use elsewhere.

1999 ◽  
Vol 36 (2) ◽  
pp. 195-209 ◽  
Author(s):  
John F Cassidy ◽  
Garry C Rogers

Three-component, digital recordings of two recent moderate earthquakes provide valuable new insight into the response to seismic shaking in the greater Vancouver area, particularly on the Fraser River delta. The 1996 M = 5.1 Duvall, Washington, earthquake (180 km southeast of Vancouver) triggered strong-motion seismographs at seven sites and the 1997 M = 4.3 Georgia Strait earthquake (37 km west of Vancouver) triggered instruments at 13 sites in the greater Vancouver area. The latter data set is especially important because it contains the first three-component recordings made on bedrock in greater Vancouver. Both data sets represent weak ground motion, with peak horizontal accelerations of 0.5-1.5% gravity (g) for the Duvall earthquake, and 0.2-2.4% g for the Georgia Strait earthquake. Using the method of spectral ratios, we estimate the site response for each of the strong-motion instrument soil sites. On the Fraser River delta amplification is observed over a relatively narrow frequency range of 1.5-4 Hz (0.25-0.67 s period), with peak amplification of 4-10 (relative to competent bedrock) for the thick soil delta centre sites, and about 7-11 for the delta edge sites. Relative to firm soil, the peak amplification ranges from 2 to 5 for the thick soil delta centre sites, and 2 to 6 for the delta edge sites. At higher frequencies, little or no amplification, and in many cases slight attenuation, is observed.Key words: seismic site response, Fraser delta, earthquakes.


1995 ◽  
Vol 85 (6) ◽  
pp. 1821-1834
Author(s):  
Toshimi Satoh ◽  
Toshiaki Sato ◽  
Hiroshi Kawase

Abstract We evaluate the nonlinear behavior of soil sediments during strong ground shaking based on the identification of their S-wave velocities and damping factors for both the weak and strong motions observed on the surface and in a borehole at Kuno in the Ashigara Valley, Japan. First we calculate spectral ratios between the surface station KS2 and the borehole station KD2 at 97.6 m below the surface for the main part of weak and strong motions. The predominant period for the strong motion is apparently longer than those for the weak motions. This fact suggests the nonlinearity of soil during the strong ground shaking. To quantify the nonlinear behavior of soil sediments, we identify their S-wave velocities and damping factors by minimizing the residual between the observed spectral ratio and the theoretical amplification factor calculated from the one-dimensional wave propagation theory. The S-wave velocity and the damping factor h (≈(2Q)−1) of the surface alluvial layer identified from the main part of the strong motion are about 10% smaller and 50% greater, respectively, than those identified from weak motions. The relationships between the effective shear strain (=65% of the maximum shear strain) calculated from the one-dimensional wave propagation theory and the shear modulus reduction ratios or the damping factors estimated by the identification method agree well with the laboratory test results. We also confirm that the soil model identified from a weak motion overestimates the observed strong motion at KS2, while that identified from the strong motion reproduces the observed. Thus, we conclude that the main part of the strong motion, whose maximum acceleration at KS2 is 220 cm/sec2 and whose duration is 3 sec, has the potential of making the surface soil nonlinear at an effective shear strain on the order of 0.1%. The S-wave velocity in the surface alluvial layer identified from the part just after the main part of the strong motion is close to that identified from weak motions. This result suggests that the shear modulus recovers quickly as the shear strain level decreases.


1994 ◽  
Vol 84 (1) ◽  
pp. 1-15 ◽  
Author(s):  
John Boatwright

Abstract The vertical components of the S wave trains recorded on the Eastern Canadian Telemetered Network (ECTN) from 1980 through 1990 have been spectrally analyzed for source, site, and propagation characteristics. The data set comprises some 1033 recordings of 97 earthquakes whose magnitudes range from M ≈ 3 to 6. The epicentral distances range from 15 to 1000 km, with most of the data set recorded at distances from 200 to 800 km. The recorded S wave trains contain the phases S, SmS, Sn, and Lg and are sampled using windows that increase with distance; the acceleration spectra were analyzed from 1.0 to 10 Hz. To separate the source, site, and propagation characteristics, an inversion for the earthquake corner frequencies, low-frequency levels, and average attenuation parameters is alternated with a regression of residuals onto the set of stations and a grid of 14 distances ranging from 25 to 1000 km. The iteration between these two parts of the inversion converges in about 60 steps. The average attenuation parameters obtained from the inversion were Q = 1997 ± 10 and γ = 0.998 ± 0.003. The most pronounced variation from this average attenuation is a marked deamplification of more than a factor of 2 at 63 km and 2 Hz, which shallows with increasing frequency and increasing distance out to 200 km. The site-response spectra obtained for the ECTN stations are generally flat. The source spectral shape assumed in this inversion provides an adequate spectral model for the smaller events (Mo < 3 × 1021 dyne-cm) in the data set, whose Brune stress drops range from 5 to 150 bars. For the five events in the data set with Mo ≧ 1023 dyne-cm, however, the source spectra obtained by regressing the residuals suggest that an ω2 spectrum is an inadequate model for the spectral shape. In particular, the corner frequencies for most of these large events appear to be split, so that the spectra exhibit an intermediate behavior (where |ü(ω)| is roughly proportional to ω).


Geophysics ◽  
2001 ◽  
Vol 66 (5) ◽  
pp. 1519-1527 ◽  
Author(s):  
Robert Sun ◽  
George A. McMechan

Reflected P‐to‐P and P‐to‐S converted seismic waves in a two‐component elastic common‐source gather generated with a P‐wave source in a two‐dimensional model can be imaged by two independent scalar reverse‐time depth migrations. The inputs to migration are pure P‐ and S‐waves that are extracted by divergence and curl calculations during (shallow) extrapolation of the elastic data recorded at the earth’s surface. For both P‐to‐P and P‐to‐S converted reflected waves, the imaging time at each point is the P‐wave traveltime from the source to that point. The extracted P‐wave is reverse‐time extrapolated and imaged with a P‐velocity model, using a finite difference solution of the scalar wave equation. The extracted S‐wave is reverse‐time extrapolated and imaged similarly, but with an S‐velocity model. Converted S‐wave data requires a polarity correction prior to migration to ensure constructive interference between data from adjacent sources. Synthetic examples show that the algorithm gives satisfactory results for laterally inhomogeneous models.


2020 ◽  
Vol 221 (2) ◽  
pp. 1029-1042 ◽  
Author(s):  
Hiroo Kanamori ◽  
Zachary E Ross ◽  
Luis Rivera

SUMMARY We use KiK-net (NIED) downhole records to estimate the radiated energy, ER, of 29 Japanese inland earthquakes with a magnitude range from Mw = 5.6 to 7.0. The method is based on the work of Gutenberg and Richter in which the time integral of S-wave ground-motion velocity-squared is measured as a basic metric of the radiated energy. Only stations within a distance of 100 km are used to minimize complex path and attenuation effects. Unlike the teleseismic method that uses mainly P waves, the use of S waves which carry more than 95 per cent of the radiated energy allows us to obtain robust results. We calibrate the method using synthetic seismograms to modernize and improve the Gutenberg–Richter method. We compute synthetic seismograms for a source model of each event with a given source function (i.e. known ER), the actual mechanism and the source-station geometry. Then, we compare the given ER with the computed energy metric to correct for the unknown effect of wave propagation and the mechanism. The use of downhole records minimizes the uncertainty resulting from the site response. Our results suggest that the currently available estimates of ER from teleseismic data are probably within a factor of 3, on average, of the absolute value. The scaled energy eR ( = ER/M0) is nearly constant at about 3 × 10−5 over a magnitude range from Mw = 5.6 to 7.0 with a slight increasing trend with Mw. We found no significant difference in eR between dip-slip and strike-slip events.


Author(s):  
Maryam Safarshahi ◽  
Igor B. Morozov

ABSTRACT Empirical models of geometrical-, Q-, t-star, and kappa-type attenuation of seismic waves and ground-motion prediction equations (GMPEs) are viewed as cases of a common empirical standard model describing variation of wave amplitudes with time and frequency. Compared with existing parametric and nonparametric approaches, several new features are included in this model: (1) flexible empirical parameterization with possible nonmonotonous time or distance dependencies; (2) joint inversion for time or distance and frequency dependencies, source spectra, site responses, kappas, and Q; (3) additional constraints removing spurious correlations of model parameters and data residuals with source–receiver distances and frequencies; (4) possible kappa terms for sources as well as for receivers; (5) orientation-independent horizontal- and three-component amplitudes; and (6) adaptive filtering to reduce noise effects. The approach is applied to local and regional S-wave amplitudes in southeastern Iran. Comparisons with previous studies show that conventional attenuation models often contain method-specific biases caused by limited parameterizations of frequency-independent amplitude decays and assumptions about the models, such as smoothness of amplitude variations. Without such assumptions, the frequency-independent spreading of S waves is much faster than inferred by conventional modeling. For example, transverse-component amplitudes decrease with travel time t as about t−1.8 at distances closer than 90 km and as t−2.5 beyond 115 km. The rapid amplitude decay at larger distances could be caused by scattering within the near surface. From about 90 to 115 km distances, the amplitude increases by a factor of about 3, which could be due to reflections from the Moho and within the crust. With more accurate geometrical-spreading and kappa models, the Q factor for the study area is frequency independent and exceeds 2000. The frequency-independent and Q-type attenuation for vertical-component and multicomponent amplitudes is somewhat weaker than for the horizontal components. These observations appear to be general and likely apply to other areas.


2020 ◽  
Author(s):  
Louise Watremez ◽  
Sylvie Leroy ◽  
Elia d'Acremont ◽  
Stéphane Rouzo

<p>The Gulf of Aden is a young and active oceanic basin, which separates the south-eastern margin of the Arabian Plate from the Somali Plate. The rifting leading to the formation of the north-eastern Gulf of Aden passive margin started ca. 34 Ma ago when the oceanic spreading in this area initiated at least 17.6 Ma ago. The opening direction (N26°E) is oblique to the mean orientation of the Gulf (N75°E), leading to a strong structural segmentation.</p><p>The Encens cruise (2006) allowed for the acquisition of a large seismic refraction dataset with profiles across (6 lines) and along (3 lines) the margin, between the Alula-Fartak and Socotra-Hadbeen fracture zones, which define a first order segment of the Gulf. P-wave velocity modelling already allowed us to image the crustal thinning and the structures, from continental to oceanic domains, along some of the profiles. A lower crustal intermediate body is observed in the Ashawq-Salalah segment, at the base of the transitional and oceanic crusts. The nature of this intermediate body is most probably mafic, linked to a post-rift thermal anomaly. The thin (1-2 km) sediment layer in the study area allows for a clear conversion of P-waves to S-waves at the top basement. Thus, most seismic refraction records show very clear S-wave arrivals.</p><p>In this study, we use both P-wave and S-wave arrivals to delineate the crustal structures and segmentation along and across the margin and add insight into the nature of the rocks below the acoustic basement. P-wave velocity modelling allows for the delineation of the structure variations across and along the margin. The velocity models are used as a base for the S-wave modelling, through the definition of Poisson’s ratios in the different areas of the models. Picking and modelling of S-wave arrivals allow us to identify two families of converted waves: (1) seismic waves converted at the basement interface on the way up, just before arriving to the OBS and (2) seismic waves converted at the basement on the way down, which travelled into the deep structures as S-waves. The first set of arrivals allows for the estimation the S-wave velocities (Poisson’s ratio) in the sediments, showing that the sediments in this area are unconsolidated and water saturated. The second set of arrivals gives us constraints on the S-wave velocities below the acoustic basement. This allows for an improved mapping of the transitional and oceanic domains and the confirmation of the mafic nature of the lower crustal intermediate body.</p>


2021 ◽  
Author(s):  
Fatma Sevil Malcıoğlu ◽  
Hakan Süleyman ◽  
Eser Çaktı

Abstract An MW 4.5 earthquake took place on September 24, 2019 in the Marmara Sea. Two days after, on September 26, 2019, Marmara region was rattled by an MW5.7 earthquake. With the intention of compiling an ample strong ground motion data set of recordings, we have utilized the stations of Istanbul Earthquake Rapid Response and Early Warning System operated by the Department of Earthquake Engineering of Boğaziçi University and of the National Strong Motion Network operated by AFAD. All together 438 individual records are used to calculate the source parameters of events; namely, corner frequency, radius, rupture area, average source dislocation, source duration and stress drop. Some of these parameters are compared with empirical relationships and discussed extensively. Duration characteristics are analyzed in two steps; first, by making use of the time difference between P-wave and S-wave onsets and then, by considering S-wave durations and significant durations. It is observed that they yield similar trends with global models. PGA, PGV and SA values are compared with three commonly used ground motion prediction models. At distances closer than about 60 km observed intensity measures mostly conform with the GMPE predictions. Beyond 60 km their attenuation is clearly faster than those of GMPEs. Frequency-dependent Q models are developed for both events. Their consistency with existing regional models are confirmed.


Author(s):  
Rong Xu ◽  
Lanmin Wang

AbstractThe horizontal-to-vertical spectral ratio (HVSR) has been extensively used in site characterization utilizing recordings from microtremor and earthquake in recent years. This method is proposed based on ground pulsation, and then it has been applied to both S-wave and ambient noise, accordingly, in practical application also different. The main applications of HVSR are site classification, site effect study, mineral exploration, and acquisition of underground average shear-wave velocity structure. In site response estimates, the use of microtremors has been introduced long ago in Japan, while it has long been very controversial in this research area, as there are several studies reporting difficulties in recognizing the source effects from the pure site effects in noise recordings, as well as discrepancies between noise and earthquake recordings. In practice, the most reliable way is the borehole data, and the theoretical site response results were compared with the HVSR using shear wave to describe site response. This paper summarizes the applications of the HVSR method and draws conclusions that HVSR has been well applied in many fields at present, and it is expected to have a wider application in more fields according to its advantages.


2020 ◽  
Author(s):  
Eser Çakti ◽  
Fatma Sevil Malcioğlu ◽  
Hakan Süleyman

<p>On 24<sup>th</sup> and 26<sup>th</sup>  September 2019, two earthquakes of M<sub>w</sub>=4.5 and M<sub>w</sub>=5.6 respectively took place in the Marmara Sea. They were associated with the Central Marmara segment of the North Anatolian Fault Zone, which is pinpointed by several investigators as the most likely segment to rupture in the near future giving way to an earthquake larger than M7.0. Both events were felt widely in the region. The M<sub>w</sub>=5.6 event, in particular, led to a number of building damages in Istanbul, which were larger than expected in number and severity. There are several strong motion networks in operation in and around Istanbul. We have compiled a data set of recordings obtained at the stations of the Istanbul Earthquake Rapid Response and Early Warning operated by the Department of Earthquake Engineering of Bogazici University and of the National Strong Motion Network operated by AFAD. It consists of 148 three component recordings, in total.  444 records in the data set, after correction, were analyzed to estimate the source parameters of these events, such as corner frequency, source duration, radius and rupture area, average source dislocation and stress drop. Duration characteristics of two earthquakes were analyzed first by considering P-wave and S-wave onsets and then, focusing on S-wave and significant durations. PGAs, PGVs and SAs were calculated and compared with three commonly used ground motion prediction models (i.e  Boore et al., 2014; Akkar et al., 2014 and Kale et al., 2015). Finally frequency-dependent Q models were estimated using the data set and their validity was dicussed by comparing with previously developed models.</p>


Sign in / Sign up

Export Citation Format

Share Document