Introduction
A detailed analysis of the M w 4.5 earthquake of 6 September 1997, which occurred near Coatzacoalcos-Minatitlán, Veracruz, is of interest for three reasons. First, the surface and the near-surface information from volcanic alignments, borehole elongations, and unpublished PEMEX seismic sections points to active normal-faulting in the region (Suter, 1991). The focal mechanism of an earthquake which occurred nearby on 26 August 1959 (M w 6.4), however, shows thrust faulting at a depth of about 26 km (Wickens and Hodgson, 1967; Suárez, 2000; see Figure 1). It is, therefore, of interest to know whether the 1997 earthquake confirms such a change in the stress regime (permutation of the least and maximum principal stresses) with depth in the region. In most regions, the stress regime at relatively shallow depth agrees with that at mid-crustal depth. There are some exceptions (see, e.g., Zoback and Zoback, 1991) and one such exception appears to be the coastal plain of the Gulf of Mexico (Frohlich, 1982; Zoback and Zoback, 1991; Suter, 1991). Another classical example is the decoupling of the stress field across the basal detachment of the Jura fold-thrust belt; the near-surface stress field is different from that in the basement (Becker et al ., 1987).
Second, the earthquake of 26 August 1959 caused serious damage to the towns of Jáltipan, Coatzacoalcos, and Minatitlán (Figueroa, 1964; Rosenblueth, 1964; Reséndiz, 1964). The latter two towns have beco-me important industrial centers related to the intense activity of PEMEX, the national petroleum company, and population has grown by 22% in the last 10 years to reach more than half million people (INEGI, 2010). For this reason, it is important to estimate ground motions that may be expected in these towns if an earthquake, such as that of 1959, were to recur in the region. We may use the records of the 1997 earthquakes obtained at the near-source broadband station of TUIG as empirical Green's function to simulate the corresponding motions from an M w 6.4 event. Although the station is about 25 km SE of Coatzacoalcos-Minatitlán, the geology of these sites is roughly similar; to a first approximation, the results for the TUIG site may be valid for the entire region in case of an earthquake at about the same focal distance from TUIG as the event of 1997.
Finally, a study of the 1997 earthquake (and other events in and along the Gulf of Mexico) has an important bearing on the seismic safety of the Laguna Verde nuclear power plant (Figure 1) as well as the hydrocarbon exploration and production facilities in this region.
Data and analysis
The 1997 earthquake was recorded by the local broadband station TUIG (S-P time = 5 s) and by seven other broadband stations of the National Seismological Service (SSN), which were located at epicentral distances greater than 260 km. The reported coda-wave magnitude, M c , is 4.3. The epicenter of the event, given by the SSN, is 18.146 °N and 94.499 °W. However, its focal depth could not be constrained. The epicenter and focal depth reported by the National Earthquake Information Center (NEIC), U.S. Geological Survey, are 18.017 °N and 94.396 °W and 33 km, respectively. This depth was fixed by NEIC in the location of the earthquake.
For moderate Mexican earthquakes it is now possible to obtain a regional centroid moment tensor (CMT) solution using relatively long-period regional waveforms (see, e.g., Pacheco and Singh, 1998). Unfortunately, the seismograms of the 1997 earthquake show little energy at periods greater than 10 s because of its relatively small magnitude. At shorter periods, a detailed three-dimensional crustal structure is needed to model the observed seismograms but is currently lacking for the region. For these reasons it was not possible to obtain a CMT solution by inverting the regional waveforms.
A single, near-source, three-component broadband recording can be used to find a reliable location, origin time, and focal mechanism of an earthquake, provided that a rough initial guess of the mechanism is available from other data (see, e.g., Kanamori et al ., 1990; Singh et al ., 1997, 2000a). The calculation of the location and origin time requires clear first arrivals on each of the three components of the ground motion as well as the knowledge of the local crustal structure. Figure 2 shows acceleration, velocity and displacement traces at TUIG during the 6 September 1997 earthquake. Figure 3 illustrates the initial part of the three components of displacement at TUIG. Note that the horizontal components have been amplified by a factor of 20. The incidence of initial P-wave at the station is nearly vertical. Nevertheless, it can be seen that the station is located in the NW quadrant with respect to the source. The amplitude towards west at time 0.172 s is 3.5 times that towards north. It follows that the station azimuth, ϕs, is 286°. Before proceeding further with the analysis of the source parameters of this earthquake, we summarize our knowledge of the local crustal structure.
Crustal structure of the region
The P-wave speed, α, in the shallow crust south and near Coatzacoalcos is available from exploration seismology carried out by PEMEX (A. Camargo, personal communication, 1997). The crust consists of a 1.8 km thick layer with a = 2.5 km/s, overlying a layer of a = 4.25 km/s. The thickness of the second layer exceeds 3.4 km, the maximum depth reached by boreholes. Based on receiver function analysis, Cruz-Atienza (2000) reports sediment thickness of 16 km below TUIG. N. Shapiro (unpublished report) inverted group velocity dispersion curve corresponding to the region between the City of Oaxaca and TUIG. In the inversion, Shapiro fixed the thickness of the first layer and a values of the first and the second layers to the values given by PEMEX. The shear-wave speeds, β, of the first and second layers were taken as 1.4 km/s and 2.4 km/s, respectively. The crustal model adopted from the results of Cruz-Atienza (2000) and N. Shapiro, and used by us in generating synthetic seismograms, is given in Table 1. In this table, the densities and the quality factors, Q , of the layers have been taken arbitrarily; the results are not very sensitive to their choices.
Source parameters of the earthquake
Based on the crustal model in Table 1 and the (S-P) time of 5 s at station TUIG, the maximum depth of the earthquake, Hmax, assuming the station to be located directly above the focus, is 30.9 km. If the thickness of the second layer in Table 1 is taken as 10 km, then Hmax becomes 34.8 km. Near-vertical incidence at TUIG may be a consequence of both small epicentral distance as compared to the source depth, and refraction of waves caused by progressively lower seismic speeds near the surface.
Figure 4, top, shows P-wave first-motion polarities at those Mexican and teleseismic stations where they could be read unequivocally. These data suggest a thrust-faulting earthquake with possible strike-slip component. The first motion polarities provide some constrains on the azimuth (140° ≤ ϕ≤ 190°) and the dip (45° ≤ δ ≤ 85°) of one of the nodal planes but the rake, l, of this plane can vary between 35° and 120°. To determine the focal mechanism, we performed a waveform inversion of the displacement traces recorded at TUIG. For the inversion the event was approximated by a point-source shear dislocation. Synthetic seismograms include near- and intermediate-field contributions (Singh et al ., 2000a). The effect of the free surface is approximately taken into account by multiplying the infinite-space synthetics by two. This approximation is acceptable if the epicentral distance, Δ, is smaller than the depth, H . We took the sta-tion azimuth, ϕ s , as 286°, the take-off angle from the source, i h , as 170°, and the angle of incidence at the surface, i O , as 5° from the vertical. Based on the observed P-pulse on the Z-component (Figure 2c), we chose a triangular source with duration, τ, of 0.38 s. We varied the azimuth, dip, and rake within the ranges mentioned above. No contradiction of the first motion polarities was allowed. The two nodal planes obtained from the inversion are: ϕ = 150°; δ = 70°; λ= 90° and ϕ = 330°; δ = 20°; λ = 90°. The observed and synthetic seismograms are shown in the bottom left of Figure 4.
We generated synthetic seismograms corresponding to this focal mechanism and the crustal model given in Table 1. Bouchon's (1982) discrete wave number algorithm was used in the computation. We again took ϕs = 286° and τ = 0.38 s. A good fit between observed and synthetic seismograms at TUIG was found for an epicentral distance of 7 km, a depth of 30 km, and M 0 = 6.0x1015 N-m (M w 4.5) (Figure 4, bottom right). Table 2 summarizes the relevant source parameters of this earthquake.
No. | Date | Latitude °N | Longitude °W | Depth km | M w | Focal Mechanism ϕ'° | δ'° | λ'° |
---|---|---|---|---|---|---|---|---|
1a | 26 Aug 1959 | 18.26 | 94.43 | 21 | 6.4 | 309 | 32 | 102 |
2b | 11 Mar 1967 | 19.23 | 95.74 | 26 | 5.7 | 250 | 39 | 20 |
3c | 24 Jul 1978 | 26.49 | 88.79 | 15 | 5.0 | 225 | 49 | 111 |
240 | 63 | 52 | ||||||
4d | 06 Sep 1997 | 18.08 | 94.47 | 30 | 4.5 | 330 | 20 | 90 |
5e | 10 Sep 2006 | 26.32 | 86.84 | 30 | 5.9 | 324 | 28 | 117 |
6e | 23 May 2007 | 21.98 | 96.31 | 24 | 5.6 | 102 | 80 | -1 |
(22.02 | 96.27 | 11 | 5.6 | 95 | 71 | -16)f | ||
(21.98 | 96.14 | 44 | 5.5 | 106 | 83 | 8)g | ||
7e | 29 Oct 2009 | 19.14 | 95.58 | 17 | 5.7 | 310 | 25 | 59 |
(18.95 | 95.69 | 16 | 5.4 | 288 | 26 | 4)g |
a Location from International Seismological Summary, ISS; depth, focal mechanism, and Mw from Suárez (2000). b Location from International Seismological Centre, ISC; depth, focal mechanism, and Mw from Suárez (2000). c Location, depth, and focal mechanism from Frohlich (1982). The two mechanisms are extreme types consistent with data first-motion data. d This study. e Global Centroid Moment Tensor (CMT) catalog. f Source parameters listed in http://www.eas.slu.edu/eqc/eqc_mt/MECH.NA/20070523190916/index.htmlg Franco et al. (2013)
Stress regime of the Gulf coast region of Coatzacoalcos-Minatitlán
In spite of some uncertainty in its focal mechanism, there is no doubt that the 1997 earthquake was a thrust event with, probably, some strike-slip component. Its source was below the Comalcalco basin at a depth of about 30 km. As mentioned earlier, the surface and the near-surface information from volcanic alignments, borehole elongations, and unpublished PEMEX seismic sections in the Gulf coast basin region of Coatzacoalcos-Minatitlán suggest active normal-faulting in the upper few kilometers (Suter, 1991). Thus, there is a change of stress regime from extension in the sediments of the upper crust to shortening in the mid- and lower crust indicated by focal mechanisms (Figure 1). A similar change of stress regime with depth is reported in the central Gulf of Mexico by Frohlich (1982) from an analysis of an earthquake in 1978 at the edge of the Mississippi Fan at a depth of 15 km (Figure 1). Other cases of stress change below sedimentary basins are discussed in Zoback and Zoback (1991). Extension in upper part of crust is probably controlled by gravitational load and major topographic gradient, at least in the adjacent Veracruz basin.
Table 2 lists seven intraplate earthquakes located in and near the Gulf of Mexico with known focal mechanisms. It includes the earthquakes of 1978 and 1997. The locations of these events and their focal mechanisms are illustrated in Figure 1. For the earthquakes of 2007 and 2009, more than one solution is available (Table 2). For these two earthquakes, the Global CMT location and focal mechanism is shown. As most intraplate regions, the Gulf is characterized by thrust-type stress regime (Zoback et al ., 1989). The orientations of the P-axes of the 1959, 1967, 2006, 2007 and 2009 earthquakes range between N30°E and N65°E. These are consistent with reported offshore stress orientations near the same region, inferred from breakouts (Zoback et al ., 1990; Suter, 1991; World Stress Map, http://dc-app3-14.gfz-potsdam.de/, 12 May 2015). These orientations fall between the directions of: (1) the relative convergence of the Cocos and North American plates and (2) the absolute motion of the North American plate. Perhaps both the relative convergence and the absolute motion are responsible for the observed P-axis, with the former playing a more dominant role for events in and near the SSW Gulf due to their relative proximity to the middle America subduction zone (earthquakes of 1959, 1967, 1997, and 2009) and the latter being dominant for the events of 2006 and 2007 which are relatively far from the plate boundary. For the 2006 earthquake, sliding of Sigsbee salt and landslide have also been suggested as possible causes (Nettles, 2007; Franco et al ., 2013).
Dewey and Suárez (1991) and Suárez (2000) suggest that the intraplate, mid- and lower-crust compression below the Comalcalco basin, as revealed by the 1959 and 1967 earthquakes, may be a consequence of strong coupling along the plate interface where the Tehuantepec ridge subducts below Mexico. This may also be true for the 1997 and 2009 earthquakes. If so, then, unlike subduction of many other ridges (Kelleher and McCann, 1976), the Tehuantepec ridge does not subduct aseismically. Since there is no clear evidence for major/great earthquakes in the Tehuantepec region during the last two centuries (see, e.g., Singh et al ., 1981), it could mean that the recurrence period of such events in this segment is much greater than in other segments along the Mexican subduction zones where it is ∼ 30 to 75 years (Singh et al ., 1981; Astiz and Kanamori, 1984).
As for the 1978 Mississippi Fan earthquake (Figure 1), the P-axis does not follow the trend of the other events. The compressional nature of this earthquake was interpreted by Frohlich (1982) as a consequence of downwarping of the lithosphere due to sediment accumulation.
Expected ground motions in the Coatzacoalcos-Minatitlán region from a postulated M w 6.4 earthquake
The earthquake of 1959 destroyed a majority of the dwellings in the town of Jáltipan (Rosenblueth, 1964). Many buildings suffered structural or foundation failures in Coatzacoalcos and Minatitlán (Marsal, 1961; Reséndiz, 1964). The land near the port of Coatzacoalcos subsided. Some of the effects of the earthquake were attributed to partial liquefaction of sand and silt (Marsal, 1961). Modified Mercalli (MM) intensities in these towns during this earthquake were VIII (Figueroa, 1964). The 1959 earthquake was not an isolated event. The epicenter of the earthquake of 11 January, 1946 was apparently close to that of 1959 (Figueroa, 1964). The earthquake of 1946 was assigned a magnitude of 6.0 (Figueroa, 1970) and a MM intensity of VII in Coatzacoalcos (Figueroa, 1964). The towns of Coatzacoalcos and Minatitlán are now important centers of national petroleum activity. Thus, it is of significant earthquake engineering interest to estimate the ground motions in these towns during a future local M w 6.4 earthquake.
To estimate the ground motions from an M w 6.4 earthquake, we used the recording of 1997 earthquake as an empirical Green's function (EGF) and a method proposed by Ordaz et al . (1995) which is based on adding N scaled EGF records, each differed in time by a random delay. The probability distribution of the delays is such that, on average, the simulations follow an ω2-spectral scaling at all frequencies. The method requires specification of the seismic moment, M 0 , and the stress drop, Ds, of both the EGF and the target earthquake. In our case, M 0 of the EGF is 6x1015 N-m and that of the target event is 5x1018 N-m (M w 6.4). A rough estimate of the static stress drop of the EGF can be obtained from the following considerations. For a circular rupture, the radius, α , of the fault can be estimated by (Boatwright, 1980):
where τ1/2 is the rise time of the far-field pulse, v is the rupture speed, c is the wave speed, and θ is the take-off angle measured from fault normal. For this event, τ1/2 is about 0.19 s (Table 2, Figure 2c). For S wave, c = β, which we take as 3.75 km/s (Table 1). Assuming f = 330°; δ = 20°; λ = 90° as the fault plane and i h =170°, θ is ∼10°. For v = 0.9β, we obtain a fault radius Δσ = 0.76 km (equation 1). For a circular fault, Ds is related to M 0 and a by (Keilis-Borok, 1959):
which gives Δσ of ∼ 6 MPa for the EGF. In the simulations, we take Δσ of the EGF and the target event to be either 6 MPa or 12 MPa. Typical simulated acceleration, velocity, and displacement traces, corresponding to Δσ = 6 MPa for both events, are illustrated in Figure 5. In the figure, we compare deterministic, synthetic displacement seismograms with ones obtained by random summation of EGF. The synthetics were generated at TUIG for an event with M 0 = 5x1018 N-m (M w 6.4), located at the same focus as the 1997 earthquake, and having the same focal mechanism. The duration of the source time function, τ, of Mexican earthquakes is related to their seismic moment by
(Singh et al ., 2000b). Thus, the estimated t for the target event is 4.2 s. A point source with a triangular source-time function of 4.2 s duration and the crustal structure given in Table 1 were taken for the computation. The PGD values are within a factor of two of each other. Both calculations show important near-field contribution (the ramp-like wave between P and S wave).
Results of simulation for various combinations of Δσ are summarized in Table 3. The expected horizontal PGA , PGV , and PGD range between 120 and 260 gal, 12 and 28 cm/s, and 6 and 11 cm, respectively.
Stres Drop, MPa EFG/Target | PGA, N | cm/s2 E | Z | PGV, N | cm/s E | Z | PGD, N | cm E | Z |
---|---|---|---|---|---|---|---|---|---|
6/6 | 146 | 173 | 103 | 13.2 | 18.7 | 6.8 | 8.8 | 6.5 | 6.1 |
6/12 | 213 | 262 | 161 | 20.0 | 27.6 | 9.2 | 10.7 | 7.2 | 6.4 |
12/12 | 171 | 198 | 124 | 17.9 | 25.8 | 8.7 | 10.6 | 7.1 | 6.3 |
12/6 | 121 | 127 | 88 | 11.8 | 17.5 | 6.0 | 8.5 | 6.4 | 6.0 |
Although the ground motions estimated by random summation of the EGF are reasonable, there are several factors which introduce uncertainties in these results. The stress drop of the EGF event is uncertain and that of the target earthquake is assumed. The directivity of the source may give rise to greater or smaller ground motions than those computed by our method. Finally, a single EGF may not be adequate to sample the entire fault plane of an M w 6.4.
We emphasize that this synthesis of the ground motion is based on the assumption of linear elastic response of the sediments. Almost certainly the shallow sediments of the Gulf basin will behave nonlinearly under such excitation, as was the case during the 1959 earthquake. In view of the uncertainties, the observed ground motions may easily differ by a factor of 2 or 3 from the expected ones. This however is usually the case in ground motion predictions.
Conclusions
Our analysis shows that the 6 September 1997 earthquake (H = 30 km; M w 4.5), like the nearby earthquake of 26 August 1957 (H = 21 km; M w 6.4), was a thrust event. The event confirms a previous conclusion that while the upper sediments of the Gulf coast basins are in an extensional (normal fault-type) stress regime, the mid- and lower crust is in a shortening (thrust fault-type) stress regime (Dewey and Suárez, 1991; Suter, 1991), which implies a permutation between the vertical and maximum horizontal principal stresses. Upper crust is under extension probably because of gravitational loads and high topographic gradient; the stress field of middle to lower crust may be controlled by far-field loads at convergent plate margin and/or absolute plate motion.
Our estimation of ground motion in the epicentral region of Comalcalco basin due to a postulated M w 6.4 earthquake indicates that peak acceleration, velocity, and displacement (assuming linear behavior of the sediments) may be in the range of 120-260 gal, 12-28 cm/s, and 6-11 cm, respectively. These, then, are also our estimations of the ground motions during the 1959, M w 6.4, earthquake. These estimations were obtained using 6 September 1997 event as an empirical Green's function under various simplifying, though reasonable, assumptions. They are valid for sites with local geology similar to that of station TUIG. The estimation of ground motions at other nearby sites, such as Minatitlán, may be calculated if their site responses are known. Subjected to such ground motions the sediments of the Gulf Coast are likely to behave nonlinearly and may liquefy, as was the case during the 1959 earthquake.