Quick
Search: 
 
advanced search
 GSW Home    GeoRef Home    My GSW Alerts    Contact GSW    About GSW    Journals List    Help 
Seismological  Research Letters Email Content Delivery
JOURNAL HOME HELP CONTACT PUBLISHER SUBSCRIBE ARCHIVE SEARCH TABLE OF CONTENTS

Seismological Research Letters; November/December 2007; v. 78; no. 6; p. 635-648; DOI: 10.1785/gssrl.78.6.635
© 2007 Seismological Society of America
This Article
Right arrow Extract
Right arrow Figures Only
Right arrow Full Text (PDF)
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Singh, S. K.
Right arrow Articles by Velázquez, M.
Right arrow Search for Related Content
GeoRef
Right arrow GeoRef Citation

A Report on the Atoyac, Mexico, Earthquake of 13 April 2007 (Mw 5.9)

S. K. Singh1,2, M. Ordaz2, J. F. Pacheco1, L. Alcántara2, A. Iglesias1, S. Alcocer2, D. García1, X. Pérez-Campos1, C. Valdes1, and D. Almora2

L. A. Aguilar2, M. Ambriz2, J. G. Anderson3, M. Ayala2, C. Cárdenas1, G. Castro2, J. L. Cruz1, R. Delgado2, L. Domínguez4, J. Estrada1, S. I. Franco1, L. E. Flores4, C. Gutierrez4, M. A. Macías2, I. Molina2, C. Morquecho4, J. Ortíz1, M. A. Pacheco-Martínez4, A. Quezada1, R. Quaas4, L. Quintanar1, C. Pérez2, J. Perez1, A. L. Ruiz2, H. Sandoval2, M. Torres2, R. Vázquez2, J. M. Velasco2, J. Velázquez2, and M. Velázquez2


    INTRODUCTION
 TOP
 INTRODUCTION
 DATA AND ANALYSIS
 SOURCE SPECTRUM AND STRESS...
 RADIATED SEISMIC ENERGY AND...
 PEAK GROUND MOTIONS
 GROUND MOTIONS IN ACAPULCO
 DAMAGE
 PERFORMANCE OF THE SEISMIC...
 NEAR REAL-TIME GROUND-MOTION MAP...
 CONCLUSIONS
 ACKNOWLEDGMENTS
 REFERENCES
 
An inslab earthquake occurred below Atoyac, Guerrero, Mexico, on 13 April 2007 (Mw 5.9) at a depth of 37 km. It was strongly felt in Acapulco ({Delta} ~65 km) and Mexico City ({Delta} ~270 km); the peak ground accelerations (PGAs) at soft sites in these two cities were as high as 200 cm/s2 and 12 cm/s2, respectively. This was the first significant earthquake in Mexico since the recent strengthening of strong-motion and broadband networks in the country. This was also the first moderate event for which ground-motion maps for Mexico City were produced and distributed without human intervention in near real time. In this paper, we report on the source characteristics and the tectonic implications of the earthquake and the ground motions that it produced. We also discuss the performance of the seismic alert system (SAS) and near real-time generation and distribution of expected ground-motion maps for Mexico City during the earthquake.


    DATA AND ANALYSIS
 TOP
 INTRODUCTION
 DATA AND ANALYSIS
 SOURCE SPECTRUM AND STRESS...
 RADIATED SEISMIC ENERGY AND...
 PEAK GROUND MOTIONS
 GROUND MOTIONS IN ACAPULCO
 DAMAGE
 PERFORMANCE OF THE SEISMIC...
 NEAR REAL-TIME GROUND-MOTION MAP...
 CONCLUSIONS
 ACKNOWLEDGMENTS
 REFERENCES
 
Our analysis of the Atoyac earthquake is based on broadband seismograms and accelerograms from stations operated by the National Seismological Service (SSN) of the Instituto de Geofísica (IG) at the Universidad Nacional Autónoma de México (UNAM), and the Instituto de Ingeniería (II), UNAM. We also used data from an accelerograph operated by the Centro Nacional de Prevención de Desastres (CENAPRED) in Acapulco.

About 35 aftershocks with M ≥ 3.4 were reported by SSN in the first four days of the earthquake, the largest being one of magnitude 5.3. The SSN network is sparse. However, the accelerograph network operated by the II is relatively dense in Guerrero (Anderson et al. 1994) where a mature seismic gap has been identified (Singh et al. 1981). Merged data from both networks permit reliable location of earthquakes in the region. As the accelerograms from the II network are not available in real time, two field crews were sent to retrieve them. Much of this report is based on the data collected up to 17 April 2007. We also analyze two events that occurred after 17 April.

The mainshock and the largest aftershock were well-recorded by the networks. The accerelograph at Atoyac (ATYC), which is located in the epicentral region (figure 1), also recorded many of the smaller aftershocks. However, only four aftershocks triggered three or more II accelerographs. Here we analyze the mainshock and these four aftershocks (table 1). We also consider six other events that occurred in the region (table 2), as they have some relevance to this study.


View this table:
[in this window]
[in a new window]

 
TABLE 1 Source parameters of the 13 April, 2007 Atoyac earthquake (Mw 5.9) and four of its aftershocks

 


View this table:
[in this window]
[in a new window]

 
TABLE 2 Source parameters of some earthquakes in the region that are relevant to this study

 
All earthquakes were located using the crustal model of Campillo et al. (1996) and the Seisan program of Lienart and Havskov (1995). For the mainshock and three of the four aftershocks, we inverted the near-source displacement traces (obtained by integration of acceleration or velocity records) at ATYC and Cayaco (CAIG) stations for focal mechanism and seismic moment. The inversion assumes that the events may be approximated by a point-source shear dislocation. Synthetic seismograms include near- and intermediate-field contributions (Singh et al. 2000). The effect of free surface is approximately taken into account by multiplying the infinite-space synthetics by two. This approximation is acceptable if the epicentral distance, {Delta}, is smaller than the depth, h. Figure 2 shows the fit between the observed and synthetic seismograms. The focal mechanisms and the moments are listed in table 1. The source parameters of the mainshock obtained from near-source seismograms and those reported in the Global Centroid Moment Tensor (CMT) catalog are very similar. The fit between observed and synthetic P wave on the Z component for the mainshock at ATYC, however, is poor. ATYC is nearly nodal for the P wave. The observed P pulse is negative for about 1 s and then becomes positive, while the synthetic P wave, as expected, is unipolar (positive). One possible explanation for the change in the polarity of the P wave is a change in the focal mechanism during rupture so that the first motion at ATYC changes from dilatational to compressional.

For the largest aftershock, which occurred on 13 April at 08:43 (event 1, table 1), and all earthquakes listed in table 2, we performed a regional moment tensor analysis following the procedure of Randall et al. (1995). The method uses a time-domain moment-tensor inversion scheme described by Langston (1981). The inversion is performed on bandpass filtered seismograms over a time window that begins with the P wave and includes surface waves. Synthetic seismograms are computed using the reflection-matrix method of Kennett (1983) and Randall (1994). The best centroid depth is obtained through grid search.


Figure 1
View larger version (34K):
[in this window]
[in a new window]

 
Figure 1. (A) Locations and focal mechanisms of Mw ≥ 4.0 earthquakes in the Guerrero region relevant to this study. M denotes Atoyac mainshock and 1, 3, and 4 refer to the aftershocks (table 1). Letters associated with the events are keyed to table 2. Except for events a, b, and f, all others occurred between 31 March and 30 April 2007. Rectangle: accelerographic station, triangle: broadband station. Some critical stations are identified. (B) Section along AA'. Dashed line indicates plate interface. Locations, focal mechanisms and associated stress axes of inslab earthquakes are indicated. Atoyac mainshock and its aftershocks (h ~37 km) represent downdip compression. The earthquakes that occur below this zone of compression exhibit downdip tension. The stress axis of a typical such event is taken from Pacheco and Singh (2007).

 
The locations and the focal mechanisms of the earthquakes are plotted in figure 1(A). Figure 1(B) illustrates a vertical section along the direction of relative convergence of the Cocos and North America (NOAM) plates (N36°E) passing through Atoyac. We note that one of the nodal planes of the Atoyac earthquake sequence is nearly vertical (table 1). The same is true for two other earthquakes that occurred near Atoyac in 2001 (events a and b, table 2; figure 1A), as well as events c and d whose epicenters are close to Acapulco (table 2; figure 1A). All of these events have similar depths (~35-43 km) and occur in the subducted Cocos plate immediately below the plate interface. It is important to note that these events reveal downdip compression. Shallow-angle thrust earthquakes occur on the plate interface at depths of less than about 25 km, updip of this inslab zone of compression (Pacheco and Singh 2007). As the vertical section in figure 1(B) illustrates, just below the zone of downdip compression there occur normal-faulting earthquakes in the slab whose stress axes show downdip tension. The plate interface above these zones of downdip tension and compression may define the downdip edge of the seismogenic, shallow-dipping thrust zone (e.g., Tichelaar and Ruff 1993). The aftershocks of large/great earthquakes may not extend below this edge (although the seismic slip during large/great earthquakes might). If so, then the aftershock area should not extend inland much beyond 10 km from the coast. A bound on downdip extent of rupture during large/great earthquakes in Guerrero is of critical importance in hazard estimation.


Figure 2
View larger version (34K):
[in this window]
[in a new window]

 
Figure 2. Observed (continuous) and synthetic (dashed) displacement seismograms of the Atoyac mainshock and three of its aftershocks at stations ATYC and CAIG (figure 1A). Observed seismograms have been inverted for seismic moments and focal mechanisms (see text).

 
Seismicity in the slab is nonexistent beyond this zone until about 180 km from the coast, where the occurrence of normal-faulting events (indicating downdip tension) resumes. These earthquakes occur at a depth of ~ 50 km over a small lateral extent (<10 km). Further inland, the inslab seismicity completely ceases and the seismic identity of the subducted plate is lost. We note that inslab downdip compression has also been reported in subduction zones of Peru and central Chile (Lemoine et al. 2002).

The seismicity pattern and state of the stress in the region can be explained by the initial bending of the Cocos plate at a shallow angle as its subduction at the middle America trench commences. The plate begins to unbend about 10 km landward from the coast, which gives rise to earthquakes with downdip compression in the upper part of the slab (like those of the Atoyac earthquake sequence) and downdip tension immediately below it. Beyond this point the plate becomes subhorizontal, reaching a depth of 50 km about 180 km from the coast where the plate again bends down. This bending is reflected in normal-faulting earthquakes with downdip tension. Thus, the inslab seismicity is mostly confined to regions of unbending and bending. The state of the stress in the slab is almost neutral in the portion that is subhorizontal. This model is similar to those proposed earlier (e.g., Suarez et al. 1990; Singh and Pardo 1993) but is now much better constrained from improved earthquake locations and focal mechanisms (Pacheco and Singh, in preparation) and from receiver function analysis using data from a linear, portable seismograph array deployed from Acapulco in the south up to 400 km in the north (Pérez-Campos et al. 2007).

Inslab earthquakes with downdip compression also occur ~10 km north of Acapulco (Pacheco and Singh 2007). Indeed, events c and d (table 2) are such earthquakes. The largest such event so far identified is the 10 December 1994 earthquake, which occurred about 41 km north of Zihuatanejo (Cocco et al. 1997).

Table 3 lists the types of seismic sources in the vicinity of Acapulco and other coastal towns in Guerrero. The table also gives observed maximum Mw in Mexico of each type of event. For the purpose of seismic hazard estimation, it will be desirable to split the instrumental seismicity catalog by different types of sources. This, of course, is a very difficult task.


View this table:
[in this window]
[in a new window]

 
TABLE 3 Seismic sources in the vicinity of Acapulco and other towns in Guerrero, and observed maximum Mw of each type of earthquake in Mexico

 

    SOURCE SPECTRUM AND STRESS DROP
 TOP
 INTRODUCTION
 DATA AND ANALYSIS
 SOURCE SPECTRUM AND STRESS...
 RADIATED SEISMIC ENERGY AND...
 PEAK GROUND MOTIONS
 GROUND MOTIONS IN ACAPULCO
 DAMAGE
 PERFORMANCE OF THE SEISMIC...
 NEAR REAL-TIME GROUND-MOTION MAP...
 CONCLUSIONS
 ACKNOWLEDGMENTS
 REFERENCES
 
We estimated source displacement and acceleration spectra, M0(f) and f2 M0(f), of the mainshock and the largest aftershock of the Atoyac earthquake from the analysis of the S-wave group recorded at hard sites within a hypocentral distance of 165 km (19 stations). The Fourier acceleration spectral amplitude of the intense part of the ground motion at a station, under far-field, point-source approximation, may be written as

Formula(1)

where

Formula(2)

In the equations above, M0(f) is the moment rate spectrum so that M0(f) -> M0 as f -> R = hypocentral distance, R{theta}{varphi} = average radiation pattern (0.55), F = free surface amplification (2.0), P takes into account the partitioning of energy in the two horizontal components Formula, β = shear-wave velocity at the source (3.95 km/s), {rho} = density in the focal region (2.9 kg/m3), and Q(f) = quality factor, which includes both anelastic absorption and scattering. The appropriate geometrical spreading term, G(R), for inslab Mexican earthquakes is R-1 and the corresponding Q(f) is 251f0.58 (García et al. 2004). Taking logarithms of equation 1 we obtain

Formula(3)

We solved equation (3) in the least-squares sense to obtain log [f2M0(f)].

The source displacement and acceleration spectra of the mainshock and the largest aftershock are shown in figure 3. We interpret these spectra within the framework of a {omega}2-source model and obtain an estimation of the seismic moment (M0) and corner frequency (fc). The stress drop ({Delta}{sigma}) is computed using the Brune (1970) model. The source spectra of the mainshock and the aftershock in figure 3 can be fit by M0 and {Delta}{sigma} of 8.9 x 1017 N-m and 47 MPa (fc = 0.727 Hz), and 1.2 x 1017 N-m and 29 MPa (fc = 1.19 Hz), respectively. We note that the stress drop of the mainshock is somewhat higher than the median {Delta}{sigma} of 30 MPa reported by García et al. (2004) for inslab normal-faulting Mexican earthquakes.


    RADIATED SEISMIC ENERGY AND APPARENT STRESS
 TOP
 INTRODUCTION
 DATA AND ANALYSIS
 SOURCE SPECTRUM AND STRESS...
 RADIATED SEISMIC ENERGY AND...
 PEAK GROUND MOTIONS
 GROUND MOTIONS IN ACAPULCO
 DAMAGE
 PERFORMANCE OF THE SEISMIC...
 NEAR REAL-TIME GROUND-MOTION MAP...
 CONCLUSIONS
 ACKNOWLEDGMENTS
 REFERENCES
 
We estimated radiated seismic energy, ER, of the mainshock from regional data by a) using aftershocks as empirical Green's functions (for details of the method see, e.g., Venkataraman et al. 2002; Singh et al. 2004), and b) direct integration of squared velocity spectra (Singh and Ordaz 1994). We also estimated ER from teleseismic data following the method of Boatwright and Choy (1986). For brevity, the details are omitted here. Estimates of ER from the three methods are 2.5, 1.6, and 0.87 x 1014 J, respectively. We consider ER = 2.5 x 1014 J, obtained from the empirical Green's function (EGF) technique, to be the best estimate. This yields ER/M0 of 2.5 x 10-4 and, with µ = 4.5 x 104 MPa, an apparent stress {sigma}a = µ ER/M0 of ~11 MPa. If, however, we allow ER to be in the range of 1 to 3 x 1014 J, then we get ER/M0 of 1 to 3 x 10-4 and {sigma}a of 4.5 to 13.5 MPa.

Choy and Kirby (2004) studied the behavior of {sigma}a of normal-faulting subduction earthquakes from a worldwide database. These authors found that the highest {sigma}a for inslab events occurred at 35-70 km depth, in proximity to zones of intense deformation such as sharp bends in the slab. The apparent stresses of these events were up to 5 MPa, significantly higher than for the interplate thrust earthquakes. Inslab earthquakes in Mexico are also significantly more energetic than interplate events (García et al. 2004). The fact that {sigma}a of the Atoyac earthquake is higher than those reported by Choy and Kirby (2004) for such events probably reflects the discrepancy between teleseismic and local estimates of the radiated energy.


Figure 3
View larger version (21K):
[in this window]
[in a new window]

 
Figure 3. Average source displacement, M0(f), and acceleration, f2M0(f), spectra (mean and ±1 standard deviation curves). (Left) mainshock, (right) largest aftershock. Superimposed on the spectra are predicted curves from an {omega}2 source model (smooth curves).

 

    PEAK GROUND MOTIONS
 TOP
 INTRODUCTION
 DATA AND ANALYSIS
 SOURCE SPECTRUM AND STRESS...
 RADIATED SEISMIC ENERGY AND...
 PEAK GROUND MOTIONS
 GROUND MOTIONS IN ACAPULCO
 DAMAGE
 PERFORMANCE OF THE SEISMIC...
 NEAR REAL-TIME GROUND-MOTION MAP...
 CONCLUSIONS
 ACKNOWLEDGMENTS
 REFERENCES
 
Figure 4 shows the quadratic mean horizontal components of peak ground acceleration (PGA) and peak ground velocity (PGV) as a function of R. These plots show data from hard sites. For Acapulco, however, we also include data from soft sites (triangles). The predicted PGA and PGV from a regression analysis of García et al. (2005) are shown in the figures. The observed values, generally, exceed those predicted by the model, which is based on inslab, normal-faulting Mexican earthquakes. This is consistent with the relatively large stress drop of the Atoyac event mentioned before. We recall that higher ground motions are expected in the epicentral area of inslab earthquakes in comparison with interplate events in Mexico (e.g., García et al. 2005).

PGA and, especially, PGV in figure 4 show a somewhat slower decay than that predicted by the model, which cannot be explained as a source effect. This results in PGAs > 1 cm/s2 and PGVs > 0.1 cm/s even at R > 300 km.


    GROUND MOTIONS IN ACAPULCO
 TOP
 INTRODUCTION
 DATA AND ANALYSIS
 SOURCE SPECTRUM AND STRESS...
 RADIATED SEISMIC ENERGY AND...
 PEAK GROUND MOTIONS
 GROUND MOTIONS IN ACAPULCO
 DAMAGE
 PERFORMANCE OF THE SEISMIC...
 NEAR REAL-TIME GROUND-MOTION MAP...
 CONCLUSIONS
 ACKNOWLEDGMENTS
 REFERENCES
 
PGAs recorded in Acapulco during the Atoyac earthquake and its major aftershock ({Delta} ~ 60-75 km) are listed in tables 4 and 5. The locations of the stations are shown in figure 5. Generally, PGAs during the mainshock at soft sites were about 100 cm/s2, although it reached 200 cm/s2 at Diana (ACAD). At hard sites (ACAJ and VNTA), the PGA was about 35 cm/s2.


View this table:
[in this window]
[in a new window]

 
TABLE 4 PGA in Acapulco during Atoyac earthquake of 13 April 2007, 05:42 (MW 5.9)

 


View this table:
[in this window]
[in a new window]

 
TABLE 5 PGA in Acapulco during Atoyac aftershock of 13 April 2007, 08:43 (MW 5.3)

 
As mentioned earlier, an earthquake of characteristics similar to the Atoyac event could occur near Acapulco. The magnitude of such an earthquake could easily reach Mw 6.5. It is of practical importance to estimate ground motions in Acapulco from postulated Mw 6.0 and 6.5 earthquakes of this type. An appropriate ground-motion prediction equation is not available for this type of earthquake. Thus, we use recordings of events c and d (table 1) in Acapulco as empirical Green's functions and apply a technique proposed by Ordaz et al. (1995) to synthesize expected ground motions. The focal mechanisms, depths, and magnitudes of events c and d qualify them as acceptable EGFs. The method assumes an {omega}2-source model. Apart from the seismic moments of the EGF and the target event, the only other needed parameters are their stress drops, which we take as 20 MPa. The results are not very sensitive to reasonable choice of this parameter. Synthesized PGA at ACAJ and ACAD as a function of Mw is shown in figure 6. The figure indicates that an Atoyac-type Mw 6.5 earthquake could cause horizontal PGA of ~ 1/2 g at soft sites such as the Acapulco Diana.


    DAMAGE
 TOP
 INTRODUCTION
 DATA AND ANALYSIS
 SOURCE SPECTRUM AND STRESS...
 RADIATED SEISMIC ENERGY AND...
 PEAK GROUND MOTIONS
 GROUND MOTIONS IN ACAPULCO
 DAMAGE
 PERFORMANCE OF THE SEISMIC...
 NEAR REAL-TIME GROUND-MOTION MAP...
 CONCLUSIONS
 ACKNOWLEDGMENTS
 REFERENCES
 
The earthquake caused damage in some municipalities in the state of Guerrero, mainly in Atoyac de Álvarez, Coyuca de Benítez, and Benito Juárez. Horizontal PGA in the first two municipalities was 215 cm/s2 and 128 cm/s2, respectively. By 22 April, at least some damage had been reported at a total of 750 houses. A quick reconnaissance of the area showed that 97% of damaged houses had been built with adobe or bahareque (rammed earth) while the remainder were built with timber, cardboard, or more durable materials such as brick and block. Damaged adobe houses in the most affected municipalities (Benito Juárez and Atoyac de Álvarez) were between 5% and 6% of the total number of houses.


Figure 4
View larger version (17K):
[in this window]
[in a new window]

 
Figure 4. Horizontal quadratic mean PGA (top) and PGV (bottom) at hard sites as a function of hypocentral distance. Superimposed are predicted curves (mean and ±1 standard deviation) for normal-faulting Mexican inslab earthquakes (García et al. 2005). Triangles: soft sites in Acapulco.

 
In general, damage in dwellings was characterized by vertical cracking at wall intersections, inclined cracking above door and window lintels, and wall-inclined cracking. Few partial collapses of roofs and tile sliding occurred. The type of damage and the known high vulnerability of adobe and rammed-earth houses again calls attention to the need to improve the earthquake performance of these houses through effective, yet inexpensive, rehabilitation techniques.

No damage to roads and bridges was recorded, except for a nine-ton rock that fell on a road in Atoyac de Álvarez.

As has occurred in other earthquakes, the population of Acapulco and Mexico City overreacted to the real intensity of the earthquake. It has been observed that even a moderate earthquake triggers fear and anxiety among people who immediately recall the consequences of the devastating 19 and 21 September 1985 earthquakes.


    PERFORMANCE OF THE SEISMIC ALERT SYSTEM FOR MEXICO CITY
 TOP
 INTRODUCTION
 DATA AND ANALYSIS
 SOURCE SPECTRUM AND STRESS...
 RADIATED SEISMIC ENERGY AND...
 PEAK GROUND MOTIONS
 GROUND MOTIONS IN ACAPULCO
 DAMAGE
 PERFORMANCE OF THE SEISMIC...
 NEAR REAL-TIME GROUND-MOTION MAP...
 CONCLUSIONS
 ACKNOWLEDGMENTS
 REFERENCES
 
The seismic alert system (SAS) for Mexico City takes advantage of the fact that the city is located more than 300 km from the foci of potentially damaging earthquakes. The system consists of 15 accelerometers located along the coast of Guerrero. An algorithm estimates magnitude of an event from the near-source accelerograms and issues public and restricted alerts for M ≥ 6 and 5 ≤ M < 6, respectively (Espinosa-Aranda and Rodríguez 2003). An evaluation of the SAS's performance by Iglesias et al. (2007) for the period 1991-2004 reveals a surprisingly high failure rate. This poor performance results from an inadequate detection algorithm and the limited area covered by the SAS. Iglesias et al. (2007) propose an alternative strategy for detecting potentially damaging earthquakes to Mexico City. One month of seismic activity in the Guerrero region, beginning 31 March, 2007, provides a further test of the present SAS algorithm and the proposed alternative strategy.

A restricted alert was issued by the SAS for the earthquake of 31 March, which was located near Acapulco (Mw 4.7; event c in table 2). This was a false alert since the magnitude was less than 5. The alert caused annoyance to the alert subscribers since the earthquake was felt by very few people in the city. A public alert was issued for the Atoyac mainshock. In a strict sense, this alert was also a failure, since Mw was less than 6. However, since the earthquake was widely felt in the city, it clearly merited a public alert. In this context, we also note that the initial estimate of the magnitude by the SSN was 6.3. The SAS correctly issued a restricted alert for the largest aftershock (Mw 5.3, event 1, table 1). Although the 19 April earthquake near Petatlan qualified for a restricted alert (Mw 5.0, event e, table 2), none was issued by SAS. Since Mw of the earthquake of 28 April near Acapulco was 5.0 (event d, table 2) the SAS should have issued a restricted alert but did not. More important, it is not clear why SAS did not issue a restricted alert for this event while it did for the 19 March earthquake (Mw 4.7). The events were located very near each other, the 28 April earthquake was greater, and it gave rise to larger accelerations at near-source stations, including those in Acapulco.

In the alternative strategy proposed by Iglesias et al. (2007), the alert is based on expected peak acceleration at station CU, Ared, which is computed from root-mean-square acceleration (Arms) at near-source stations using the relation

Formula

where RS and RCU are hypocentral distances to near-source stations and station CU, respectively. CU is located in the hill zone of Mexico City. This station has been in continuous operation since 1964. The equation is based on extensive strong-motion recordings available at CU. In the above equation, Arms is computed over a 10-s window on bandpass-filtered (0.2-1.0 Hz) accelerograms. Ared is the expected maximum acceleration at CU in the same frequency band. Experience shows that if observed peak acceleration in CU, Acu, is ≥ 2 cm/s2, then the earthquake is widely felt and may cause damage in the city. Iglesias et al. (2007) propose a single level of alert and suggest this should occur when Acu ≥ 2 cm/s2. If we require that the probability of missed alert at this threshold be 1% and that of false alert be 30%, then the alert would have to be issued if Ared ≥ 0.8 cm/s2.


Figure 5
View larger version (22K):
[in this window]
[in a new window]

 
Figure 5. Location of strong-motion stations in Acapulco. VNTA and ACAJ are situated on hard sites, while all other stations are located on soft sites.

 


Figure 6
View larger version (11K):
[in this window]
[in a new window]

 
Figure 6. Simulated PGA at ACAJ (hard site) and ACAD (soft site) in Acapulco for postulated 5.5 ≤ Mw ≤ 6.5 events of the same type as the Atoyac earthquake. Solid symbols: event c (table 2) used as EGF, open symbol: event d (table 2) used as EGF. Predicted PGA is greater when event d is used as EGF. Note that data at ACAD for event d was not available.

 


Figure 7
View larger version (11K):
[in this window]
[in a new window]

 
Figure 7. Observed PGA (ACU) versus expected PGA (Ared) at station CU. ACU and Ared are based on bandpass-filtered CU and nearsource recordings, respectively (see text). The plot includes data since the early 1960s and forms the basis of an alternative seismic alert algorithm for Mexico City (Iglesias et al. 2007). In the proposed scheme, an alert will be issued if the expected PGA at CU, estimated from near-source recording, is ≥ 2 cm/s2. If we require that the probability of missed alert at ACU ≥ 2 cm/s2 be 1% and that of false alert be 30%, then an alert must be triggered if Ared ≥ 0.8 cm/s2. The data for 30 March-30 April 2007 sequence (tables 1 and 2) are shown by dots. Only the Atoyac mainshock merits an alert.

 
Figure 7 shows Ared and ACU for the five earthquakes discussed above, along with data from previous earthquakes (Iglesias et al. 2007). During the earthquake sequence under consideration the Atoyac mainshock is the only event for which Ared ≥ 0.8 cm/s2.


    NEAR REAL-TIME GROUND-MOTION MAP FOR MEXICO CITY
 TOP
 INTRODUCTION
 DATA AND ANALYSIS
 SOURCE SPECTRUM AND STRESS...
 RADIATED SEISMIC ENERGY AND...
 PEAK GROUND MOTIONS
 GROUND MOTIONS IN ACAPULCO
 DAMAGE
 PERFORMANCE OF THE SEISMIC...
 NEAR REAL-TIME GROUND-MOTION MAP...
 CONCLUSIONS
 ACKNOWLEDGMENTS
 REFERENCES
 
Since mid-2006, a new system has been put in operation to produce, in automatic fashion, the expected ground-motion map for Mexico City during earthquakes.

The system, presently in its first stage, relies on recordings of a single station. It is triggered by the detection of an earthquake at the station CU. As mentioned in the previous section, this station, located in the hill zone of Mexico City, has been in operation since the 1960s and has been used as a reference site to compute spectral amplifications of all the other instrumented sites throughout the city (see figure 8).

Once the end of the motion is detected, the system computes the response spectra (pseudoaccelerations, 5% damping) of the horizontal components and decides, based on their absolute size and the ratio between the spectral ordinates at 1 s and the corresponding PGA, whether it is a significant event or not.

If the system decides that it is dealing with a significant event, a process starts to compute estimated average horizontal response spectra at the nodes of a grid of 1,600 points with separation of approximately 500 m, using response spectral transfer functions that are precomputed and stored in the system. These transfer functions have been obtained using those computed from previous earthquakes at the approximately 100 strong-motion stations of the Mexico City Accelerographic Network (figure 8) and a Bayesian interpolation scheme that uses the transfer functions at instrumented sites, plus information from about 2,500 points in which the predominant site period has been measured. Previous studies have shown that the response spectra thus estimated at sites in the Valley of Mexico are very close to the exact response spectra (e.g., Ordaz et al. 1988).


Figure 8
View larger version (58K):
[in this window]
[in a new window]

 
Figure 8. Map of Mexico City showing its geotechnical subdivision in hill, transition, and lake-bed zones. Rectangles, dots, and triangles are accelerographic stations operated by three different institutions. Relative spectral amplifications of these sites with respect to CU (large dot) are known. Interpolated relative spectral amplifications are available at 2,500 grid points. Real-time accelerograms at CU are used to generate expected ground-motion maps of the city via random vibration theory.

 
Once the average horizontal response spectra are computed at the 1,600 points, maps of spectral accelerations at four periods (0 s, 0.2 s, 1 s, and 2 s) are produced. These maps, along with key parameters of the ground motion, are: 1) uploaded to an Internet site (http://aplicaciones.iingen.unam.mx/webSAPS/Defaultsn.aspx); 2) sent by e-mail to a list of city officials and interested persons; 3) sent by cellular phone to the same list of persons; and 4) sent by radio links to the Civil Protection City Authority (still in a test phase).

As an example, figure 9 shows two spectral acceleration maps of the 19 September 1985 earthquake. These maps, of course, have been recreated, a posteriori, using the recordings obtained at station CU at that time.

Although the system had been triggered by several minor events, the earthquake of 13 April 2007 was the first moderate event for which ground-motion maps were produced without human intervention. The total time employed by the system was seven minutes after it first detected the motion. In figure 10 we present the maps corresponding to PGA and the spectral acceleration response at 2 s. Note that the general color of the maps is blue, indicating moderate accelerations throughout the city. Compare the colors of these maps with those of the 1985 maps in figure 9.


Figure 9
View larger version (50K):
[in this window]
[in a new window]

 
Figure 9. (A) Distribution of estimated PGA (in cm/s2) in Mexico City for the earthquake of 19 September 1985. (B) Distribution of estimated spectral acceleration at 2 s (in cm/s2) in Mexico City for the earthquake of 19 September 1985.

 


Figure 10
View larger version (25K):
[in this window]
[in a new window]

 
Figure 10. (A) Distribution of estimated PGA (in cm/s2) in Mexico City for the Atoyac earthquake of 13 April 2007, computed and distributed in near real time. (B) Distribution of estimated spectral acceleration at 2 s (in cm/s2) in Mexico City for the Atoyac earthquake of 13 April 2007, computed and distributed in near real time.

 

    CONCLUSIONS
 TOP
 INTRODUCTION
 DATA AND ANALYSIS
 SOURCE SPECTRUM AND STRESS...
 RADIATED SEISMIC ENERGY AND...
 PEAK GROUND MOTIONS
 GROUND MOTIONS IN ACAPULCO
 DAMAGE
 PERFORMANCE OF THE SEISMIC...
 NEAR REAL-TIME GROUND-MOTION MAP...
 CONCLUSIONS
 ACKNOWLEDGMENTS
 REFERENCES
 
The Atoyac earthquake (Mw 5.9), located about 10 km inland from the coast at a depth of 37 km, occurred in the subducted Cocos plate. The focal mechanisms of the event and its aftershocks correspond to downdip compression. A detailed study of the seismicity of the region shows that smaller events of similar focal mechanism occur frequently at about the same distance from the coast in the 35-40 km depth range (Pacheco and Singh 2007). The occurrence of this type of earthquake may be taken as evidence of unbending of the subducting Cocos plate after an initial bending, which begins at the trench. The unbending results in along-slab compression in the upper part of the subducted slab and gives rise to both steeply dipping thrust and slab-push normal-faulting earthquakes.

Shallow-dipping thrust events on the plate interface cease to occur at a depth of ~25 km, just above the zone of Atoyactype earthquakes. This may define the downdip edge of the seismogenic, shallow-dipping thrust zone. It has been suggested that aftershocks of large/great earthquakes may not extend any farther below this edge. If so, then the aftershock area should not extend inland much beyond 10 km from the coast. A bound on downdip extent of rupture during large/great earthquakes is of critical importance to seismic hazard estimation in Guerrero.

The source spectrum of the earthquake is reasonably well-fit by the {omega}2-source model along with a stress drop of 47 MPa. This stress drop is somewhat larger than the median value of 30 MPa estimated for inslab normal-faulting Mexican earthquakes. This relatively large stress drop also explains larger than predicted PGA and PGV recorded during the earthquake.

Radiated seismic energy, ER, during the earthquake was in the range of 1 to 3 x 1014 J. This yields ER/M0 of 1 to 3 x 10-4 and apparent stress of 4.5 to 13.5 MPa. These values are higher than those reported for inslab earthquakes in subduction zones using teleseismic data ({sigma}a < 5 MPa). The difference probably reflects systematic error in teleseismic estimation of ER.

The damage in the epicentral region was minor and limited to adobe houses.

Synthesis of ground motion in Acapulco from a postulated Atoyac-type Mw 6.5 earthquake suggests that such an event might give rise to a horizontal PGA of ~1/2 g at soft sites such as the Acapulco Diana.

The performance of the seismic alert system (SAS) during a one-month period beginning 31 March 2007 was mixed. Public and restricted alerts for the Atoyac mainshock (Mw 5.9) and its largest aftershock (Mw 5.3), respectively, may be considered successes. The system, however, issued a false restricted alert for the Acapulco earthquake of 31 March (Mw 4.7) and missed restricted alerts for the earthquakes of 19 and 28 April (Mw 5.0). We tested a recently proposed alternative strategy for SAS (Iglesias et al, 2007) on the same earthquake sequence and found that it performed well. We think that this alternative strategy merits serious consideration.

Estimated ground-motion maps for Mexico City were automatically generated and distributed in seven minutes after an accelerograph at CU, a station in the hill zone of the city, first detected the earthquake. Presently in the first stage of its development, this project will provide quick and reliable assessment of possible damage in the city.

The Atoyac earthquake demonstrates that permanent seismic instrumentation in some parts of Mexico, especially Guerrero, is reasonably dense for fruitful research in seismology and earthquake engineering. The near real-time availability of strong-motion and other geophysical data (e.g., GPS and tide-gauge), however, still remains elusive and a goal to strive for.


    ACKNOWLEDGMENTS
 TOP
 INTRODUCTION
 DATA AND ANALYSIS
 SOURCE SPECTRUM AND STRESS...
 RADIATED SEISMIC ENERGY AND...
 PEAK GROUND MOTIONS
 GROUND MOTIONS IN ACAPULCO
 DAMAGE
 PERFORMANCE OF THE SEISMIC...
 NEAR REAL-TIME GROUND-MOTION MAP...
 CONCLUSIONS
 ACKNOWLEDGMENTS
 REFERENCES
 
The broadband network of SSN and the accelerographic network of II are partially funded by Secretaria de Gobernacion, Mexico. The II network is also partially funded by the National Science Foundation. We thank CENAPRED for providing us with data from its Acapulco station, ACAJ. The research was partly supported by DGAPA (project IN110207) and CONACyT (project 42671-T).

1 Instituto de Geofísica, Universidad Nacional Autónoma de México, C.U., Coyoacán 04510, Mexico City, Mexico Back

2 Instituto de Ingeniería, Universidad Nacional Autónoma de México, C.U., Coyoacán 04510, Mexico City, Mexico Back

3 Seismological Laboratory, Mackay School of Mines, University of Nevada, Reno, Nevada 89557 Back

4 Centro Nacional de Prevención de Desastres, Delfin Madrigal 665, Coyoacán 04360, Mexico City, Mexico Back


    REFERENCES
 TOP
 INTRODUCTION
 DATA AND ANALYSIS
 SOURCE SPECTRUM AND STRESS...
 RADIATED SEISMIC ENERGY AND...
 PEAK GROUND MOTIONS
 GROUND MOTIONS IN ACAPULCO
 DAMAGE
 PERFORMANCE OF THE SEISMIC...
 NEAR REAL-TIME GROUND-MOTION MAP...
 CONCLUSIONS
 ACKNOWLEDGMENTS
 REFERENCES
 

Anderson, J. G., J. N. Brune, J. Prince, R. Quaas, S. K. Singh, D. Almora, P. Bodin, M. Oñate, R. Vásquez, and J. M. Velasco (1994). The Guerrero accelerograph network. Geofisica International 33,341 -372.

Boatwright, J., and G. L. Choy (1986). Teleseismic estimates of the energy radiated by shallow earthquakes. Journal of Geophysical Research 91,2,095 -2,112.[CrossRef]

Brune, J. N. (1970). Tectonic stress and the spectra of seismic shear waves from earthquakes. Journal of Geophysical Research 75,4,997 -5,009.[CrossRef]

Campillo, M., S. K. Singh, N. Shapiro, J. Pacheco, and R. B. Hermann (1996). Crustal structure of the Mexican volcanic belt, based on group velocity dispersion. Geophysical Journal International 35,361 -370.

Choy, G. L., and S. H. Kirby (2004). Apparent stress, fault maturity and seismic hazard for normal-fault earthquakes at subduction zones. Geophysical Journal International 159,991 -1012.[CrossRef][Web of Science][GeoRef]

Cocco, M., J. F. Pacheco, S. K. Singh, and F. Courboulex (1997). The Zihuatanejo, Mexico earthquake ofDecember 10, 1994 (M=6.6): Source characteristics and tectonic implications. Geophysical Journal International 131,135 -145.[CrossRef][Web of Science][GeoRef]

Espinosa-Aranda, J. M., and F. H. Rodríguez (2003). The seismic alert system of Mexico City. In International Handbook of Earthquake and Engineering Seismology, ed. W. H. K. Lee,1,253 -1,259 (chapter 76). Amsterdam and Boston: Academic Press.

García, D., S. K. Singh, M. Herraiz, J. F. Pacheco, and M. Ordaz (2004). Inslab earthquakes of Central Mexico: Q, source spectra and stress drop. Bulletin of the Seismological Society of America 94,789 -802.[Abstract/Free Full Text][CrossRef][Web of Science][GeoRef]

García, D., S. K. Singh, M. Herraiz, M. Ordaz, and J. F. Pacheco (2005). Inslab earthquakes of Central Mexico: peak ground-motion parameters and response spectra. Bulletin of the Seismological Society of America 95,2,272 -2,282.[Abstract/Free Full Text][CrossRef][Web of Science][GeoRef]

Iglesias, A., S. K. Singh, M. Ordaz, M. A. Santoyo, and J. F. Pacheco (2007). The seismic alert system for Mexico City: An evaluation of its performance and a strategy for its improvement. Bulletin of the Seismological Society of America 97(5),1,718 -1,729.[Abstract/Free Full Text][CrossRef][Web of Science][GeoRef]

Kennett, B. L. N. (1983). Seismic Wave Propagation in Stratified Media. Cambridge: Cambridge University Press, 342 pps.

Langston, C. (1981). Source inversion of seismic waveforms: The Koyna, India, earthquake of 13September , 1967. Bulletin of the Seismological Society of America 71,1 -24.[Abstract/Free Full Text][Web of Science][GeoRef]

Lemoine, A., R. Madariaga, and J. Campos (2002). Slab-pull and slabpush earthquakes in the Mexican, Chilean and Peruvian subduction zones. Physics of the Earth and Planetary Interiors 132,157 -175.[CrossRef][Web of Science][GeoRef]

Lienart, B. R., and J. Haskov (1995). A computer program for locating earthquakes both locally and globally. Seismological Research Letters 66, 26-36.[GeoRef]

Ordaz, M., J. Arboleda, and S. K. Singh (1995). A scheme of random summation of an empirical Green's function to estimate ground motions from future large earthquakes. Bulletin of the Seismological Society of America 85,1,635 -1,647.[Abstract/Free Full Text][Web of Science][GeoRef]

Ordaz, M., S. K. Singh, E. Reinoso, J. Lermo, J. M. Espinosa, and T. Domínguez (1988). Estimation of response spectra in the lake bed zone of the Valley of Mexico. Earthquake Spectra 4,815 -834.[CrossRef][GeoRef]

Pacheco, J. F., and S. K. Singh (2007). Seismicity and state of stress in the subduction zone of Guerrero, Mexico. Jt. Assem. Suppl., Abstract: S43B-01, EOS—Transactions of the American Geophysical Union 88,23 .

Pérez-Campos, X., R. W. Clayton, A. Iglesias, A. Husker, and C. M. Valdés-González (2006). MASE: A seismological perspective of the sub-horizontal subduction of the Cocos Plate under North America, Jt. Assem. Suppl., Abstract: T13F-0,. EOS—Transactions of the American Geophysical Union 87,52 .

Randall, G. E. (1994). Efficient calculation of composite differential seismograms for laterally homogeneous Earth models. Geophysical Journal International 118,255 -259.[Web of Science][GeoRef]

Randall, G. E., C. J. Ammon, and T. J. Owens (1995). Moment tensor estimation using regional seismograms from a Tibetan plateau portable network deployment. Geophysical Research Letters 22,1,665 -1,668.[CrossRef][Web of Science][GeoRef]

Singh, S. K., J. Havskov, and L. Astiz (1981). Seismic gaps and recurrence periods of large earthquakes along the Mexican subduction zone. Bulletin of the Seismological Society of America 71,827 -843.[Abstract/Free Full Text][Web of Science][GeoRef]

Singh, S. K., and M. Ordaz (1994). Seismic energy release in Mexican subduction zone earthquakes. Bulletin of the Seismological Society of America 84,1,533 -1,550.[Abstract/Free Full Text][Web of Science][GeoRef]

Singh, S. K., M. Ordaz, J. F. Pacheco, and F. Courboulex (2000). A simple source inversion scheme for displacement seismograms recorded at short distances. Journal of Seismology 4,267 -284.[CrossRef][Web of Science][GeoRef]

Singh, S. K., J. F. Pacheco, B. K. Bansal, X. Pérez-Campos, R. S. Dattatrayam, and G. Suresh (2004). A source study of the Bhuj, India, earthquake of 26 January 2001 (Mw 7.6). Bulletin of the Seismological Society of America 94,1,195 -1,206.[Abstract/Free Full Text][CrossRef][Web of Science][GeoRef]

Singh, S. K., and M. Pardo (1993). Geometry of the Benioff zone and state of stress in the overriding plate in central Mexico. Geophysical Research Letters 20,1,483 -1,486.[Web of Science][GeoRef]

Suárez, G., T. Monfret, G. Wittlinger, and C. David (1990). Geometry of subduction and depth of the seismogenic zone in the Guerrero gap, Mexico. Nature 345,336 -338.[CrossRef][GeoRef]

Tichelaar, B. and L. J. Ruff (1993). Depth of seismic coupling along subduction zones. Journal of Geophysical Research 98,2,017 -2,037.[CrossRef]

Venkataraman, A., L. Rivera, and H. Kanamori (2002). Radiated energy from 16 October 1999 Hector Mine earthquake; regional and teleseismic estimates. Bulletin of the Seismological Society of America 92,1,256 -1,265.[Abstract/Free Full Text][CrossRef][Web of Science][GeoRef]

Instituto de Geofísica
Universidad Nacional Autónoma de México
C.U. 04510. México D.F.
krishna{at}ollin.igeofcu.unam.mx
(S. K. S.)





This Article
Right arrow Extract
Right arrow Figures Only
Right arrow Full Text (PDF)
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Singh, S. K.
Right arrow Articles by Velázquez, M.
Right arrow Search for Related Content
GeoRef
Right arrow GeoRef Citation


JOURNAL HOME HELP CONTACT PUBLISHER SUBSCRIBE ARCHIVE SEARCH TABLE OF CONTENTS
Copyright © 2009 by Seismological Society of America