Earthquake Ground Motion and Structural Response in Alluvial Valleys

This paper is concerned with the problem of soil amplification and structural damage due to local site conditions in sedimentary valleys during earthquakes. It focuses on a small valley in Kirovakan, for which 1D wave propagation analyses have failed to provide adequate answers for the large extent and spatial distribution of damage during the 1988 Armenia earthquake. A more realistic 2D finite-element analysis is performed herein in search of an explanation of the observed behavior. Using as input an inferred rock accelerogram, the response of the valley is calculated for a vertically incident SH-wave. Synthetic accelerograms of the surface ground motion are presented for different sites; these accelerograms are then used to determine the amplification ratios of the surface response with respect to that of the free-field motion of the rock outcrop, for different frequencies and for a continuous set of sites. In addition, response spectra are evaluated for simple oscillators representing structures located at various sites. Results of the 2D simulations show striking differences with respect to those from 1D analyses. In particular, (1) while the resonant frequencies exhibited by a 1D model for a given site also appear in the 2D model, the peak ground response and structural response are almost twice as large for the 2D as for the 1D model; (2) the 2D model exhibits, in addition, a new set of resonant frequencies and concomitant “mode shapes” across the valley, which are directly related to its finite width; and (3) due to these additional resonances the ground amplification ratio tends to oscillate very rapidly, both spatially and with frequency, leading to the observation that two identical structures located in the same vicinity or two slightly different structures located essentially on the same site can be subjected to significantly different seismic forces, even if the underlying soils have very similar characteristics. These results provide a meaningful explanation for the observed damage, and thus serve to exemplify a situation in which site effects caused by the finite lateral extent of a valley must be taken into consideration in order to model satisfactorily seismic behavior.

[1]  Paul Christiano,et al.  Implementation of effective seismic input for soil−structure interaction systems , 1988 .

[2]  Z. Alterman,et al.  Propagation of elastic waves in layered media by finite difference methods , 1968 .

[3]  Keiiti Aki,et al.  Local site effects on weak and strong ground motion , 1993 .

[4]  Jonathan Richard Shewchuk,et al.  Triangle: Engineering a 2D Quality Mesh Generator and Delaunay Triangulator , 1996, WACG.

[5]  H. L. Wong,et al.  Surface motion of a semi-elliptical alluvial valley for incident plane SH waves , 1974, Bulletin of the Seismological Society of America.

[6]  Jean-Pierre Bardet,et al.  Engineering observations on ground motion at the Van Norman Complex after the 1994 Northridge earthquake , 1996 .

[7]  V. W. Lee.  Scattering of Plane Sh‐Waves By A Semi‐Parabolic Cylindrical Canyon In an Elastic Half‐Space , 1990 .

[8]  George Gazetas,et al.  1988 ARMENIA EARTHQUAKE. II: DAMAGE STATISTICS VERSUS GEOLOGIC AND SOIL PROFILES. DISCUSSION AND CLOSURE , 1995 .

[9]  Unified Symmetric BEM‐FEM for Site Effects on Ground Motion—SH Waves , 1991 .

[10]  A. Bayliss,et al.  Radiation boundary conditions for wave-like equations , 1980 .

[11]  The spitak (armenia, ussr) earthquake of 7 december 1988: A summary engineering seismology report , 1989 .

[12]  Arthur Frankel,et al.  Three-dimensional simulations of ground motions in the San Bernardino Valley, California, for hypothetical earthquakes on the San Andreas fault , 1993, Bulletin of the Seismological Society of America.

[13]  Arthur Frankel,et al.  Site response for urban Los Angeles using aftershocks of the Northridge earthquake , 1996, Bulletin of the Seismological Society of America.

[14]  M. Cara,et al.  The Spitak (Armenia) earthquake of 7 December 1988: field observations, seismology and tectonics , 1989, Nature.

[15]  George Gazetas,et al.  Seismological, soil and valley effects in Kirovakan, 1988 Armenia earthquake , 1994 .

[16]  Donald V. Helmberger,et al.  Elastic finite-difference modeling of the 1971 San Fernando , 1988 .

[17]  M. D. Trifunac,et al.  Erratum to surface motion of a semi-cylindrical alluvial valley for incident plane SH waves , 1972, Bulletin of the Seismological Society of America.

[18]  David M. Boore,et al.  Comparison of two independent methods for the solution of wave-scattering problems: Response of a sedimentary basin to vertically incident SH waves , 1971 .

[19]  D. Boore,et al.  Finite Difference Methods for Seismic Wave Propagation in Heterogeneous Materials , 1972 .

[20]  John Lysmer,et al.  The propagation of Love waves across nonhorizontally layered structures , 1971, Bulletin of the Seismological Society of America.

[21]  1988 ARMENIA EARTHQUAKE. I: SEISMOLOGICAL, GEOTECHNICAL, AND STRUCTURAL OVERVIEW , 1994 .

[22]  Paul Christiano,et al.  On the effective seismic input for non-linear soil-structure interaction systems , 1984 .

[23]  Tatsuo Ohmachi,et al.  Love-wave propagation in a three-dimensional sedimentary basin , 1992, Bulletin of the Seismological Society of America.

[24]  K. Foss.  COORDINATES WHICH UNCOUPLE THE EQUATIONS OF MOTION OF DAMPED LINEAR DYNAMIC SYSTEMS , 1956 .

[25]  Pierre-Yves Bard,et al.  The seismic response of sediment-filled valleys. Part 2. The case of incident P and SV waves , 1980 .

[26]  M. Dravinski,et al.  Application of a hybrid method for scattering of P, SV, and Rayleigh waves by near-surface irregularities , 1987 .

[27]  Ground motion on alluvial valleys under incident plane SH waves , 1991, The Journal of the Australian Mathematical Society. Series B. Applied Mathematics.

[28]  A. H. Hadjian.  The Spitak, Armenia earthquake of 7 December 1988 — why so much destruction , 1993 .

[29]  Keiiti Aki,et al.  Surface motion of a layered medium having an irregular interface due to incident plane SH waves , 1970 .

[30]  Paul C. Jennings.  Earthquake response of tall regular buildings , 1997 .

[31]  P. Bard,et al.  The two-dimensional resonance of sediment-filled valleys , 1985 .

[32]  Jacobo Bielak,et al.  Symmetric Local Absorbing Boundaries in Time and Space , 1991 .

[33]  J. Lawrence Von Thun,et al.  Earthquake engineering and soil dynamics II : recent advances in ground-motion evaluation : proceedings of the specialty conference , 1988 .

[34]  David R. O'Hallaron,et al.  Earthquake ground motion modeling on parallel computers , 1996, Supercomputing '96.

[35]  Francisco J. Sánchez-Sesma,et al.  An indirect boundary element method applied to simulate the seismic response of alluvial valleys for incident P, S and Rayleigh waves , 1993 .

[36]  G. Gazetas,et al.  Ground‐Motion and Soil‐Response Analyses for Leninakan, 1988 Armenia Earthquake , 1994 .

[37]  Hiroshi Kawase,et al.  Time-domain response of a semi-circular canyon for incident SV, P, and Rayleigh waves calculated by the discrete wavenumber boundary element method , 1988 .

[38]  Francisco J. Sánchez-Sesma,et al.  Ground motion on alluvial valleys under incident plane SH waves , 1979, Bulletin of the Seismological Society of America.

[39]  L. Pérez-Rocha,et al.  Diffraction of elastic waves by three-dimensional surface irregularities. Part II , 1989 .

[40]  Masahiro Iida,et al.  The seismic coda, site effects, and scattering in alluvial basins studied using aftershocks of the 1986 North Palm Springs, California, earthquake as source arrays , 1993 .