Seismic site effects

Fig.1 : Seismic site effects / wave amplification in a horizontal layer (SH-waves): various wavefields.

When propagating, the seismic waves are reflected and refracted at the interface between the various geological layers (Fig.1).

The example of Figure 1 depicts the seismic wave amplification in horizontal geological layers. We consider a homogeneous elastic half-space (in green) over which an elastic alluvial layer of constant thickness ${\displaystyle h}$ is located (in gray). A shear wave (${\displaystyle SH}$) of amplitude ${\displaystyle A_{2}}$ reaches the interface between the half-space and the alluvial layer with an incidence ${\displaystyle \theta _{2}}$. It thus generates:

• a reflected wave in the half-space with amplitude ${\displaystyle A_{2}^{'}}$ and incidence ${\displaystyle \theta _{2}}$
• a refracted wave in the superficial layer with amplitude ${\displaystyle A_{1}}$ and incidence ${\displaystyle \theta _{1}}$

The refracted wave originates a reflected wave when reaching the free surface; its amplitude and incidence are denoted ${\displaystyle A_{1}^{'}}$ and ${\displaystyle \theta _{1}}$ respectively. This latter wave will be reflected and refracted several times at the base and the top of the surficial layer. If the layer is softer than the half-space, the surface motion amplitude can be larger than ${\displaystyle A_{2}}$ thus leading to the amplification of seismic waves or seismic site effects. When the geological interfaces are not horizontal, it also possible to study seismic site effects by considering the basin effects due to the complex geometry of the alluvial filling[2] For small inclinations of the subsurface layers and/or low impedance contrasts, the assumption of horizontal layering (i.e. the 1D assumption) can still be used to predict site response.[3]

In this article, we propose several examples of seismic site effects (observed or simulated during large earthquakes) as well as a theoretical analysis of the amplification phenomenon.

Fig.2 : Site effects in Mexico city: recordings from the 1985 earthquake

Seismic site effects have been first evidenced during the 1985 Mexico City earthquake.[4] The earthquake epicenter was located along the Pacific Coast (several hundreds kilometers from Mexico-City), the seismic shaking was however extremely strong leading to very large damages.

Figure 2 displays the recordings performed at different distances from the epicenter during the earthquake sequence. The acceleration amplitude measured at different distances changes drastically:

• Campos station: this station is located very close to the epicenter and recorded a maximum acceleration of ${\displaystyle 150~cm/s^{2}}$,
• Teacalco station: this station is located at more than 200 km from the epicenter and recorded a much lower acceleration (about ${\displaystyle 18~cm/s^{2}}$). This amplitude decay is due to the wave attenuation during

Figure 2a: The effect of the resonance: the natural frequency of the upper geological bed of the modern area of Mexico City (SCT) is 0.5Hz (period of 2 seconds), that's why the PGA reached its maximum amplitude at the same period. The lower curve is the UNAM spectral accelerations.

the propagation process:[1] geometrical attenuation due to the expansion of the wavefront and material (or intrinsic) attenuation due to the energy dissipation within the medium (e.g. grains friction),

• UNAM station: this station is located at more than 300 km from the epicenter and recorded a maximum acceleration of ${\displaystyle 35~cm/s^{2}}$, larger than that recorded at the Teacalco station,
• SCT station: this station is located in Mexico City at approximately 400 km from the epicenter and recorded a very strong maximum acceleration (about ${\displaystyle 170~cm/s^{2}}$).

We may notice that the acceleration amplitude strongly decreases first and then increases when the seismic waves reach the alluvial deposit on which Mexico City has been founded.

Figure 2a: Shows the effect of resonance: The thickness of the upper geological bed of the modern area of Mexico City is 40 m. The velocity of the shear waves through that layer is 80 m/sec.[5] It means that the natural frequency of that formation is 0.5 Hz (period of 2 seconds).[6] When the shear waves of the same frequency arrived to that area, the resonance was responsible for that huge seismic site effect.

In case of horizontal soil layering (constant thickness, cf Fig.1), we may analyze seismic site effects theoretically. One considers a shear wave (${\displaystyle SH}$) (i.e. polarized perpendicularly to the figure) reflected and refracted wave at the interface between both media and reflected at the free surface.

Considering Fig.1, we may analyze the propagation of the various waves in the sedimentary layer (${\displaystyle i=1}$) and in the half-space (${\displaystyle i=2}$). Assuming both media as linear elastic and writing the continuity conditions at the interface (displacement and traction) as well as the free surface conditions, we may determine the spectral ratio ${\displaystyle {\bar {T}}(\omega )}$ between the surface motion and the motion at the top of the half-space without any sedimentary layer:

${\displaystyle {\bar {T}}(\omega )={\frac {2A_{1}}{2A_{2}}}={\frac {1}{\cos k_{z_{1}}h+i{\bar {\chi }}\sin k_{z_{1}}h}}}$

where ${\displaystyle k_{z_{1}}={\frac {\omega \theta _{i}}{V_{S_{i}}}}}$  ; ${\displaystyle {\bar {\chi }}={\sqrt {\frac {\mu _{1}\rho _{1}}{\mu _{2}\rho _{2}}}}{\frac {\cos \theta _{1}}{\cos \theta _{2}}}}$ and :

• ${\displaystyle h}$ is the layer thickness,
• ${\displaystyle \theta _{i}}$ is the wave incidence in layer ${\displaystyle i}$,
• ${\displaystyle \rho _{i}}$ is the mass density in layer ${\displaystyle i}$,
• ${\displaystyle \mu _{i}}$ is the shear modulus in layer ${\displaystyle i}$,
• ${\displaystyle k_{z_{1}}}$ is the vertical wave number in layer 1,
• ${\displaystyle V_{S_{i}}={\sqrt {\frac {\mu _{i}}{\rho _{i}}}}}$ is the shear wave velocity.

Fig. 3: Seismic site effects in a single sedimentary layer (SH waves): spectral ratio for various layer/bedrock velocity ratios

Fig.3 displays the variations of the spectral ratio ${\displaystyle {\bar {T}}}$ with respect to frequency for different mechanical features of the half-space (with ${\displaystyle V_{S_{1}}=200~m/s}$ for the sedimentary layer). We notice that the motion amplification may be very strong at certain frequencies. The amplification level depends on the velocity contrast ${\displaystyle {\bar {\chi }}}$ and takes the following maximum values:

• ${\displaystyle |{\bar {T}}_{max}|=2}$ for ${\displaystyle V_{S_{2}}=800~m/s}$ (blue curve),
• ${\displaystyle |{\bar {T}}_{max}|\simeq 3.5}$ for ${\displaystyle V_{S_{2}}=2000~m/s}$ (green curve),
• ${\displaystyle |{\bar {T}}_{max}|\simeq 6}$ for ${\displaystyle V_{S_{2}}=5000~m/s}$ (yellow curve).

The red curve corresponds to a large velocity contrast between the layer and the half-space (${\displaystyle {\bar {\chi }}\gg 1}$); the amplification is thus very large. As displayed in Fig.3, the maximum amplification is reached at certain frequencies corresponding to the resonance of the sedimentary layer. The fundamental frequency of the layer (or 1st resonance frequency) may be easily calculated[1] under the form: ${\displaystyle f_{0}={\frac {V_{S_{1}}}{4h}}}$. The fundamental mode thus corresponds to a quarter wavelength resonance. The "quarter wavelength" approach can be used to estimate site amplifications due to the impedance contrast.[7]

When the sedimentary layers are not horizontal (e.g. sedimentary basin), the analysis is more complex since surface waves generated by the lateral heterogeneities (e.g. basin edges) should be accounted for. In such cases, it is possible to perform empirical studies but also theoretical analyses for simple geometries[8] or numerical simulations for more complex cases.[9]

Fig.4: Seismic site effects in Caracas (BEM simulations in the frequency domain).

In sedimentary basins, site effects also lead to the generation of surface waves at the basin edges. This phenomenon may significantly strengthen the amplification of the seismic motion. The aggravation of the amplification level when compared to the case of horizontal layering may be up to a factor of 5 or 10. It depends on the velocity contrast between the layers and the geometry of the basin.[9] Such phenomena are named basin effects and we may consider the analogy with the vibrations in a bowl of jelly.

The theoretical analysis of site effects in canyons or semi-circular sedimentary basins has been performed through semi-analytical methods in the early 80's.[8] Recent numerical simulations [10] allowed the analysis of site effects in ellipsoidal sedimentary basins. Depending on the basin geometry, the aggravation of site effects is different from that of the horizontally layered case.

When the mechanical properties of the sedimentary basin are known, we may simulate site effects numerically. Figure 4 depicts the amplification phenomenon for the city of Caracas.[11][12] The amplification level of a plane wave (${\displaystyle SH}$) is computed by the Boundary Element Method in the frequency domain.[13] Each color map displays the amplification level ${\displaystyle A_{0}}$ at a given frequency ${\displaystyle f_{0}}$ :

• top: ${\displaystyle f_{0}=0.3Hz~;~A_{0}=2.53}$. Site effects due to the topography clearly occur at the top of the hill (right). Nevertheless, site effects due to the sedimentary basin lead to a larger amplification.
• middle: ${\displaystyle f_{0}=0.4Hz~;~A_{0}=8.83}$. Topographical site effects are negligible when compared to that due to the basin (4 times larger than at 0.3 Hz).
• bottom: ${\displaystyle f_{0}=0.6Hz~;~A_{0}=7.11}$. Site effects in the basin are of the same order than at 0.4 Hz but we notice a much shorter wavelength.

Numerous geological sites have been investigated by various researchers for weak earthquakes as well as for strong ones (cf synthesis[1]). In the latter case, it is necessary to account for the nonlinear behavior of the soil under large loading[14] or even the soil liquefaction which may lead to the soil failure.

• Earthquakes
• Earthquake engineering
• Plate tectonics
• Seismic hazard
• Seismic intensity scales
• Seismic magnitude scales
• Seismic risk
• Seismic waves
• Seismology
• Near-surface geophysics