Abstract
The main-shock (Mw, 6.3) and the aftershocks of the ‘Les Saintes’ earthquake sequence (French Indies) were analyzed to quantify high-frequency directivity effects. A correction method was applied to isolate source spectra within a large frequency range (0.5 to 25 Hz). Most of the aftershocks source spectra are fully consistent with a Brune spectrum point-source shape and do not show any azimuthal dependence. The main-shock (Mw, 6.3) and the two largest aftershocks (Mw, 5.8, 5.3) show, however, a clear azimuthal dependence that indicates significant directivity effect. The discrepancy of the radiated spectral energy and the change in the corner frequencies introduced by directivity effects show that such an effect is significant at high frequency (from 1 to 25 Hz). Our data suggest that the amplitudes in the main-shock Fourier spectrum at directive sites are around a factor of 2.5 higher with respect to anti-directive sites.
Similar content being viewed by others
1 Introduction
Finite-source effects such as low-frequency directivity effects have been well known for decades (Haskell 1964; Boore and Joyner 1989; Somerville et al. 1997; Seekins and Boatwright 2010). However, high-frequency (above 1 to 2 Hz) directivity remains an open debate (Boatwright 2007; Cultrera et al. 2009). The controversy concerning whether directivity has a relevant impact on high frequencies is widely discussed, including the theoretical considerations that support opposite opinions: While some studies claim that steady-state rupture propagation would enhance directivity effects by the assumption of quasi-deterministic rupture behavior (Boore and Joyner 1989; Gallovic and Burjanek 2007; Ruiz et al. 2011), other studies claim that incoherencies in the rupture, because subevents are located randomly on the fault plane with random microscale rupture directivity (stochastic rupture behavior), might strongly attenuate directivity effects (Bernard and Herrero 1994; Somerville et al. 1997; Boatwright et al. 2002; Spudich and Chiou 2008; Ameri et al. 2012; Ruiz et al. 2011). From an observational perspective, Boatwright et al. (2002) evaluated the variation in the corrected velocity spectra with the azimuth, and observed directivity up to 1 Hz. Ameri et al. (2012) suggested directivity effects observed up to 2 Hz at very close stations for strong ground motions recorded during the Mw 6.3 2009 L’Aquila earthquake.
Several fundamental earthquake properties can be measured from the spectral content of seismic-wave arrivals. For example, directivity effects cause variations in the height of the spectra plateau, as the earthquake energy will focus along the rupture direction (Gallovic and Burjanek 2007), as well as a shift in the spectral corner frequency (Boore and Joyner 1989) and the duration of the apparent source-time function is shorter than the non-directive sites. Figure C1 in Ruiz et al. (2011) demonstrates clearly how the acceleration source spectra are expected to change for directive and anti-directive sites under different assumptions at high frequency. Under quasi-deterministic rupture behavior assumption, the expected spectral plateau is the highest at directive site and the lowest at anti-directive site. On the contrary, under stochastic rupture behavior condition, the spectral plateaus are merged together at high frequency from a theoretical perspective. However, because of the variety of masking mechanisms, such as source–receiver path effects (which include the effects of geometric spreading and anelastic attenuation along the ray path), and because of the station site responses, earthquake source spectra are inaccessible to direct observation (Boore and Joyner 1989). The general lack of data covering the whole focal sphere also results in ambiguous observational evidence for the effects of directivity.
On November 21, 2004, at 11:41 UTC, a magnitude Mw 6.3 earthquake (Harvard Global Centroid Moment Tensor Catalog http://www.globalcmt.org/CMTsearch.html) struck offshore (10 km south) of the ‘Les Saintes’ islands in Guadeloupe, French West Indies (Fig. 1). There were more than 30,000 aftershocks recorded over the following years, most of which occurred at shallow depths near the islands of the archipelago. The main-shock and its aftershocks were recorded by the French Accelerometric Network stations, plus the accelerometric stations of the ‘Conseil Général de Martinique’. All the accelerometric data are freely available at http://www-rap.obs.ujf-grenoble.fr/ (Péquegnat et al. 2008).
The resulting high-quality dataset provides a good opportunity to access information on earthquake sources, regional attenuation, and local site effects. The recent study of Drouet et al. (2011) provided updated regional attenuation relationships and amplification factors of the accelerometric stations in this area. Through this new information, we have been able to analyze the corrected source spectra of the Les Saintes earthquake sequence and to evaluate how the source affects control of the observed spectra. Based on the main-shock study of Feuillet et al. (2011), accelerometric stations are either located in directive sites (Guadeloupe) or anti-directive sites (Martinique), and this suitable station coverage allows us to analyze the azimuthal dependency of the radiated energy.
Our research strategy is straightforward. We first computed the corrected source spectra by removing path and site effects. Secondly, as the source spectra can be affected by the combined effects of radiation pattern effects and directivity effects, we removed potential radiation pattern effects. Finally, we compared and analyzed the high-frequency azimuthal dependency of the Fourier spectra of the main-shock and all of the aftershock sequence earthquakes characterized by Mw > 5 and high-quality recordings for both directive and anti-directive sites (Table 1).
2 Analysis of propagation and site-corrected spectra: evidence of a source-directivity effect?
The ability to compute isolated source spectra by removing path and site effects is a key prior need for the detection of high-frequency directivity effects. Given the shape of the Brune spectrum, acceleration spectra are more sensitive to high-frequencies (f > 1 Hz) directivity effects. Thus, we computed and analyzed the average S-wave acceleration Fourier spectra from the Les Saintes earthquake sequence.
Each observed S-wave acceleration Fourier spectrum from source i and receiver j is a product of the source term, the regional propagation path effect, and the local site transfer function, as shown in Eq. (1):
Where, A ijk(r ij,f k) is the acceleration Fourier spectra observed from source i at receiver j, for a given frequency fk. The source term is described using the Brune’s source model, as shown in Eq. (2):
where, \( {M}_{0_i} \) is the seismic moment, and \( {f}_{c_i} \) is the corner frequency of event i.
The attenuation along the ray path rij includes the effects of geometric spreading and anelastic attenuation, as shown in Eq. (3):
where, VS is the average S-wave velocity along the path (assumed to be 3.5 km/s here), Q(fK) is the frequency-dependent quality factor, γ is the geometric spreading exponent, and Sj(fk) is the local site effect at each receiver j.
In the log domain, Eq. (1) can be written as given in Eq. (4):
where, \( {y}_{\mathrm{i}\mathrm{j}\mathrm{k}}\kern-0.3em =\kern-0.3em { \log}_{10}\left[{\mathrm{A}}_{\mathrm{i}\mathrm{j}\mathrm{k}}\left({\mathrm{r}}_{\mathrm{i}\mathrm{j}},{\mathrm{f}}_{\mathrm{k}}\right)\right],{\mathrm{m}}_{0_{\mathrm{i}}}\kern-0.3em =\kern-0.3em { \log}_{10}\left[{\mathrm{M}}_{{\mathrm{o}}_{\mathrm{i}}}\kern-0.3em \times \kern-0.3em \frac{2\left\langle {\mathrm{R}}_{\uptheta \upphi}\right\rangle }{4\pi \rho {v_{\mathrm{S}}}^{\mathrm{S}}}\right]\kern-0.45em \), and 〈Rθϕ〉 refers to the average source radiation pattern (0.55 for S-waves).
To compute the source term, we can rewrite Eq. (4) as Eq. (5):
S-wave Fourier spectra are computed from the time window that starts at the S-wave arrival time and ends where it includes 80 % of the energy computed from the S-wave arrival time. The spectra were smoothed using the Konno and Ohmachi (1998) smoothing procedure, following Drouet et al. (2011). Then, we applied the correction described in Eq. (4), to remove propagation and site terms and to isolate the source spectra. We adopt the attenuation models determined by Drouet et al. (2011) (Table 2) to define the anelastic attenuation (Qo, α) and geometric spreading (γ) parameters. The station site terms are also taken from Drouet et al. (2011). The attenuation parameters differ for travel paths towards Guadeloupe and Martinique because the paths are crossing different regions and different portions of the crust. Indeed, the paths towards Martinique are longer than those towards Guadeloupe and sample a deeper part of the crust which leads to a lower attenuation.
The path- and site-corrected source spectra clearly fit the respective Brune’s source spectral models for small-to-moderate events (Fig. 2a–e). The models developed by Drouet et al. (2011) explain these point-source spectra well, and there is no evident divergence between the various observations from different azimuthal station groups. Thus, for small-to-moderate earthquakes, directivity effects are not observed. This good fit also confirms the accuracy of regional attenuation parameters and the local site amplification derived by Drouet et al. (2011).
The records of the three largest shocks show different behavior (Fig. 2f–h). As the magnitude increases, significant source spectra azimuthal variations are observed. The height of the amplitude plateau is a function of the azimuth, which suggests a directivity effect. The main-shock event shows the clearest azimuthal variations. However, when trying to detect directivity effects, and thus any associated azimuthal dependency in the corrected source spectra, we need to first remove any possible radiation pattern influence.
3 Do source radiation patterns explain the observed azimuthal variations?
The radiation patterns refer to the angular dependency of the wave amplitudes from the seismic source, and they are thus a function of the take-off angles, focal mechanism, and the azimuth angle from the event to the station events.
To estimate the regional take-off angles, we adopt the widely used velocity model for the Les Saintes area (Dorel 1978). The main-shock centroid moment tensor focal mechanism was taken from Courboulex et al. (2010), and the radiation pattern coefficients were calculated using the equations given by Aki and Richards (2002).
The path and site-corrected source spectra are finally corrected from the radiation pattern effect using Eq. (6):
where R(θ,φ) is the radiation pattern coefficient at take-off angle φ and azimuth angle θ.
The corresponding parameters for each station are shown in Table 3. Here, we compute the radiation pattern coefficients for S H and S V waves and then take the geometric mean as the final radiation coefficients (RAD) since in the study, we have analyzed average S-wave acceleration data.
These computed radiation patterns are only valid at low frequencies (Liu and Helmberger 1985). Indeed, Takenaka et al. (2003) suggested that the frequency-dependent distortion of the S-wave radiation pattern might be caused by the mixing and coupling of the horizontal (S H) and vertical (S V) S-waves in the heterogeneous structure near the source region. By analyzing dense KiK-net array observations from the Tottori-Ken Seibu earthquake, Mw 6.6, and its aftershocks, Takemura et al. (2009) demonstrated the collapse of the S-wave front due to seismic-wave scattering in a heterogeneous structure and showed that the radiation pattern is more isotropic at high frequencies (>2 Hz). Castro et al. (2006) suggested that the S H-wave radiation approaches the theoretical radiation only for frequencies below 0.5 Hz. The evaluation of the frequency range at which radiation patterns start to become isotropic is thus still a debated issue.
Consequently, we evaluated the radiation pattern correction of the main-shock for two cases. First, we applied the correction within the whole frequency range (0.1 to 25 Hz), then, only for the low frequency range (0.1 to 1 Hz).
Figure 3 shows the raw source spectra for the main-shock as well as the source spectra corrected according to the radiation pattern. The radiation pattern correction highlights even more the azimuthal variations of the source spectra, especially when the correction is applied to the whole frequency band. These results confirm that directivity effects are the only physical phenomena that can explain the observed main-shock azimuthal variations.
4 Results
To analyze further these directivity effects, we use the slip and kinematic source model developed by Feuillet et al. (2011) (Fig. 1). This slip model is characterized by two main slip zones that are located 5 to 10 km to the southeast and northwest of the hypocenter. Only a small portion of the total seismic moment was released in the hypocentral area. The largest and first rupture patch (southeast of the hypocenter) propagated towards the surface. The second patch rupture propagation to the northwest (towards the western part of Guadeloupe Island) is supported by the spectral amplitudes azimuthal variations shown in this study and by the analysis of Drouet et al. (2011). Indeed, the Drouet et al. (2011) analysis suggests a potential directivity effect to explain the azimuthal variation of the residuals and the large discrepancy between the residuals observed in Guadeloupe and Martinique. Moreover, Boatwright (2007) has shown that the aftershocks pattern is partly controlled by directivity effects. The fact that most of the aftershocks are located north of the main-shock epicenter then also suggests a directivity effect toward the northwest.
As the northwestern patch ruptures towards Guadeloupe (more or less along the strike direction), the stations located in Guadeloupe (Fig. 1, red and green triangles) are directive sites (especially the stations located along the strike). Stations sited in the Martinique area (Fig. 1, blue triangles) can be considered as anti-directive sites for the first patch. Moreover, the southeastern patch ruptures towards the surface and the faults dip to the east, so we also consider the station located on Marie-Galante (Fig. 1, GBGA, yellow triangle) as an anti-directive site (Fig. 4a). To quantify these directivity effects, the mean values of the source spectra of each station group (Fig. 4b) were computed for the main-shock. As stated previously, many studies have shown that the radiation pattern becomes isotropic at high frequencies, and these results were obtained by applying the radiation pattern correction at low frequency only (<1 Hz).
Following the aftershock distribution, the inferred rupture direction is then consistent with the corrected source spectra behavior (Fig. 1): We observe the largest spectral amplitude at the directive sites, and especially at stations sited along the rupture direction. The lowest spectral levels are shown at both of the anti-directive sites, which confirm that the azimuthal dependence of the corrected spectra is caused by directivity effects.
To quantify these directivity effects, the means of the source spectra of each station group (Fig. 4) were computed for the main-shock. The frequency band below 0.5 Hz is dominated by noise for small events, and therefore the propagation and site-effect corrections were not calculated by Drouet et al. (2011). At intermediate frequencies (between 0.5 and 2 Hz), both directive and anti-directive sites have a ‘sag’ in the source spectra (Fig. 4). At high frequencies (>2 Hz), the spectral plateaus obtained for the different station groups show large differences, and there is a ratio of about 2.0 to 2.5 between the directive and non-directive stations. These observations confirm the source models, which predict high directivity effects at high frequencies (Gallovic and Burjanek 2007). We computed the best fit for the source spectrum at the directive and anti-directive sites according to the Brune’s source spectrum with different corner frequencies. We observe that at directive sites, the fitted corner frequency is higher than the corner frequency at anti-directive sites (Fig. 4). We also note that the source spectral shapes obtained at directive and anti-directive sites cannot be explained by a Brune’s source spectral shape anymore, probably because of the source complexity, and so the analysis of these corner frequencies in terms of stress drop or rupture velocity is difficult.
5 Conclusion
Applying appropriate corrections for path and site effects, and for radiation patterns, allows us to detect the impact of directivity effects on S-wave Fourier spectra. The point-source Brune’s model explains the observed source spectra of small earthquakes (Mw < 5.3), while the discrepancies between the observations and the predicted Brune’s models caused by finite-source effects are significant for the main-shock and largest aftershocks (Mw, 6.3, 5.8, 5.3, respectively). Our study shows significant amplification (2.0- to 2.5-fold) of the acceleration spectral plateau due to directivity effects at high frequencies (2 to 25 Hz). This evidence of high-frequency (>1 Hz) directivity effects contributes to the few observational studies that have been suggesting such effects (Boatwright et al. (2002), up to 1 Hz (Ameri et al. (2012)), and up to 2 Hz (Courboulex et al. (2013)).
References
Aki K, Richards PG (2002) Quantitative seismology. Univ Science Books
Ameri G, Gallovič F, Pacor F (2012) Complexity of the Mw 6.3 2009 L’Aquila (Central Italy) earthquake: 2. Broadband strong-motion modeling.
Bernard P, Herrero A (1994) Slip heterogeneity, body-wave spectra, and directivity of earthquake ruptures. Ann Geophys 37(6)
Boatwright J (2007) The persistence of directivity in small earthquakes. Bull Seismol Soc Am 97(6):1850–1861
Boatwright J, Choy GL, Seekins LC (2002) Regional estimates of radiated seismic energy. Bull Seismol Soc Am 92(4):1241–1255
Boore DM, Joyner WB (1989) The effect of directivity on the stress parameter determined from ground motion observations. Bull Seismol Soc Am 79(6):1984–1988
Castro RR, Franceschina G, Pacor F, Bindi D, Luzi L (2006) Analysis of the frequency dependence of the S-wave radiation pattern from local earthquakes in central Italy. Bull Seismol Soc Am 96(2):415–426
Courboulex F, Converset J, Balestra J, Delouis B (2010) Ground-motion simulations of the 2004 mw 6.4 les saintes, guadeloupe, earthquake using ten smaller events. Bull Seismol Soc Am 100(1):116–130
Courboulex F, Dujardin A, Vallée M, Delouis B, Sira C, Deschamps A, Honoré L, Thouvenot F (2013) High‐frequency directivity effect for an Mw 4.1 earthquake, widely felt by the population in southeastern France. Bull Seismol Soc Am 103(6):3347–3353
Cultrera G, Pacor F, Franceschina G, Emolo A, Cocco M (2009) Directivity effects for moderate-magnitude earthquakes (Mw 5.6–6.0) during the Umbria-Marche sequence, central Italy. Tectonophysics 476(1–2):110–120
Dorel J (1978) Sismicité et structure de l’arc des Petites Antilles et du bassin atlantique.
Drouet S, Bouin MP, Cotton F (2011) New moment magnitude scale, evidence of stress drop magnitude scaling and stochastic ground motion model for the French West Indies. Geophys J Int 187:1625–1644. doi:10.1111/j.1365-246X.2011.05219.x
Feuillet N, Beauducel F, Jacques E, Tapponnier P, Delouis B, Bazin S, Vallée M, King G (2011) The Mw = 6.3, November 21, 2004, Les Saintes earthquake (Guadeloupe): Tectonic setting, slip model and static stress changes. Journal of Geophysical Research: Solid Earth (1978–2012) 116 (B10)
Gallovic F, Burjanek J (2007) High frequency directivity in strong ground motion modeling methods. Ann Geophys 50:203–211
Haskell N (1964) Total energy and energy spectral density of elastic wave radiation from propagating faults. Bull Seismol Soc Am 54(6A):1811–1841
Konno K, Ohmachi T (1998) Ground-motion characteristics estimated from spectral ratio between horizontal and vertical components of microtremor. Bull Seismol Soc Am 88(1):228–241
Liu HL, Helmberger DV (1985) The 23: 19 aftershock of the 15 October 1979 Imperial Valley earthquake: more evidence for an asperity. Bull Seismol Soc Am 75(3):689–708
Péquegnat C, Guéguen P, Hatzfeld D, Langlais M (2008) The French accelerometric network (RAP) and national data centre (RAP-NDC). Seismol Res Lett 79(1):79–89
Ruiz J, Baumont D, Bernard P, Berge‐Thierry C (2011) Modelling directivity of strong ground motion with a fractal, k − 2, kinematic source model. Geophys J Int 186(1):226–244
Seekins LC, Boatwright J (2010) Rupture directivity of moderate earthquakes in Northern California. Bull Seismol Soc Am 100(3):1107–1119
Somerville PG, Smith NF, Graves RW, Abrahamson NA (1997) Modification of empirical strong ground motion attenuation relations to include the amplitude and duration effects of rupture directivity. Seismol Res Lett 68(1):199–222
Spudich P, Chiou BSJ (2008) Directivity in NGA earthquake ground motions: analysis using isochrone theory. Earthquake Spectra 24:279
Takemura S, Furumura T, Saito T (2009) Distortion of the apparent S‐wave radiation pattern in the high‐frequency wavefield: Tottori‐Ken Seibu, Japan, earthquake of 2000. Geophys J Int 178(2):950–961
Takenaka H, Mamada Y, Futamure H (2003) Near-source effect on radiation pattern of high-frequency S waves: strong SH–SV mixing observed from aftershocks of the 1997 Northwestern Kagoshima, Japan, earthquakes. Phys Earth Planet Inter 137(1):31–43
Acknowledgments
The comments of the associate editor and two anonymous reviewers significantly improved the submitted manuscript. We thank Ross Stein and Mathieu Causse for their useful suggestions and encouraging comments. This study originates from Yen-Shin’s master’s degree dissertation research project in Earthquake Engineering and Engineering Seismology (MEEES Consortium: www .meees.org). This study would not have been possible without the huge amount of work from the RAP network staff who provides the seismological community with high-quality accelerometric data. We especially thank Marie-Paule Bouin and our colleagues of the Observatoire Volcanologique et Sismologique de Guadeloupe for their support.
Author information
Authors and Affiliations
Corresponding author
Additional information
Data and resources
All accelerometric data are freely available at http://www-rap.obs.ujf-grenoble.fr/.
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution License which permits any use, distribution, and reproduction in any medium, provided the original author(s) and the source are credited.
About this article
Cite this article
Chen, Y., Letort, J., Cotton, F. et al. High-frequency directivity effects: evidence from analysis of the Les Saintes records. J Seismol 18, 457–466 (2014). https://doi.org/10.1007/s10950-014-9419-2
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s10950-014-9419-2