Print version ISSN 0016-7169
Geofís. Intl vol.48 no.2 México Apr./June 2009
Lowcost accelerograph units as earthquake alert devices for Mexico City: how well would they work?
D. García*, S. K. Singh, A. Iglesias, L. Quintanar and C. Valdés
Instituto de Geofísica, Universidad Nacional Autónoma de México, Del. Coyoacán, 04510 Mexico City, México * Corresponding author: email@example.com
Received: October 17, 2008
Accepted: December 5, 2008
Recientemente se ha propuesto la utilización de unidades acelerográficas de bajo coste como dispositivo de alerta sísmica para la Ciudad de México. Este tipo de unidades dispararía la alarma cuando la aceleración del suelo alcanzara un nivel prefijado, en principio 4 gal, lo que para sismos importantes sucedería durante la onda P o en el inicio de su coda. En este trabajo se ha evaluado el desempeño de este tipo de unidades a partir de los acelerogramas registrados en el periodo 1985 2008 en las estaciones SCT y CDAO, ubicadas en la zona de lago de Ciudad de México. Como se preveía, el tiempo de alerta (la diferencia entre el tiempo de disparo del dispositivo y el de inicio de la parte intensa del movimiento), TA, aumenta, para sismos regionales, con la aceleración máxima, PGA. Para los registros analizados (con PGA > 4 gal) TA varía entre 20 y 45 s, correspondiendo el máximo TA al terremoto de Michoacán de 1985. Salvo algunas excepciones TA es superior a 15 s para todos aquellos sismos que produjeron PGA > 20 gal. Este dispositivo podría ser útil para mitigar los daños humanos de un sismo siempre y cuando un tiempo de alerta entre 15 y 45 s fuese suficiente para poner en práctica medidas de seguridad previamente establecidas, lo cual podría ser práctico para escuelas de una a tres plantas. Además el dispositivo podría utilizarse para detener y/o cerrar de forma segura instalaciones críticas antes de la llegada del movimiento sísmico de mayor amplitud.
Palabras clave: Alerta sísmica, tiempo de alerta, aceleración máxima, Ciudad de México.
Recently, lowcost accelerograph units have been proposed as earthquake alert devices for Mexico City. These units would trigger when the acceleration reaches a preestablished level, presumably 4 gal. For significant earthquakes, this would occur during P wave or in the early part of its coda. We test the performance of such a unit on accelerograms recorded in the period 1985 2008 at SCT and CDAO, two lakebed sites in Mexico City. As expected, the alert time (the time of arrival of intense ground motion minus the trigger time), TA, for regional earthquakes is found to increase with the PGA. TA of the recorded accelerograms (with PGA > 4 gal) ranges between about 20 and 45 sec; the largest values of TA correspond to the 1985 Michoacan earthquake. With some exceptions, TA is greater than 15 sec for earthquakes which produced PGA > 20 gal. The device may be useful in mitigating injuries and loss of lives if an alert time of 15 to 45 sec is sufficient to put in to effect preestablished safety measures. This may be practical for one to threestorey public school buildings. The device would be useful for shutting critical facilities before the arrival of largeamplitude ground motion.
Key words: Seismic alert, alert time, peak ground acceleration, Mexico City.
On 19 September, 2008, the 23rd anniversary of the Michoacan earthquake of 1985, the authorities of Mexicos Federal District announced plans to install 10,000 units of an earthquake alert device in Mexico City. The proposed device would trigger when the acceleration at the site of installation reaches a threshold, presumably 4 gal. For large earthquakes one expects this to occur during Pwave or its coda, before the arrival of largeamplitude Swave group, thus providing an alert with some lead time. Here we investigate the performance of such a device using recorded strongmotion data of earthquakes during 1985 2008 at two lakebed sites in Mexico City: Secretaría de Comunicaciones y Transportes (SCT) and Central de Abastos (CDAO). Peak ground acceleration (PGA) during these recordings exceeded 4 gal on, at least, one of the three orthogonal components. We assume that the trigger level of the proposed device is 4 gal. Our goal is to investigate the alert time that such device would provide. The conclusion based on SCT and CDAO data may, generally, apply to other lakebed sites in the Valley of Mexico also.
As it is well known, an early warning system for Mexico City from subduction earthquakes along the Guerrero coast has been in operation since August 1991 (EspinosaAranda et al., 1995). The system, called the Seismic Alert System (SAS), issues public and restricted alerts for M > 6 and 5 < M < 6 earthquakes, respectively. To some extent the proposed devices would compete and their function would overlap with that of the SAS. For this reason, we also discuss the performance of the SAS based on the same strongmotion dataset.
Expected Alert Time
At regional distances, largeamplitude ground motion follows the arrival of S wave. Thus, if a sensor triggers during P wave or its coda, then there is some lead time before the arrival of the intense part of the ground motion. This lead time can be used as an alert time. If an alert is sought only for earthquakes which produce large ground motions and, hence, pose threat to the society, then the sensor should trigger only when the motion exceeds a certain threshold. At a fixed threshold the sensor may trigger: (1) after the arrival of S wave, (2) in the coda of P wave, or (3) at the arrival of P wave. Case (1) provides no alert time. The alert time for case (2) is less than (S P) time. The largest alert time, equal to (S P) time, is obtained in case (3). For earthquakes occurring at the same focus, the strength of the source, as reflected by the PGA at the site, would determine which of the cases would be realized. In other words, the alert time would be directly related to the PGA. Thus, potentially more destructive earthquakes (larger PGA) would have larger alert times. It is clear, however, that the alert time will always be less or equal than (S P) time. It also follows that for local earthquakes the alert time will always be small, irrespective of PGA.
We note that the intense motion at regional distances, often, does not begin at S wave. For this reason, we will measure the arrival time of intense motion as the time when the energy in the signal becomes 5% of the total energy of the record.
SCT and CDAO accelerographic stations have been in operation since 1985. The instrumentation at both sites, including trigger level and preevent memory, has changed over time. The accelerograph at SCT has been operated at trigger level, Ar and preevent memory that have varied between 2 and 7 gal and between 4 and 40 sec, respectively. At CDAO, AT has varied between 1 and 5 gal, while the preevent memory has remained fixed at 4 sec. We note that even when AT was set at 7 gal, the records include that portion of strong motion where the acceleration first reaches 4 gal. This is because of the preevent memory. In particular, records with AT = 7 gal correspond to local events, thus a preevent memory of 4 sec was enough to record completely the intense part of the ground motion. During the great Michoacan earthquake of 1985, AT and preevent memory at both sites were 4 gal and 4 sec, respectively. Clearly, the values of AT and preevent time would affect the record length before acceleration reaches 4 gal.
Table 1 lists the events which produced at least 4 gal at SCT and/or CDAO and Fig. 1 shows the location of these events. PGA values at both stations are given in Table 2. We note that PGA was greater than or equal to 4 gal on, at least, one component during 31 earthquakes at SCT and 32 at CDAO in the period September 1985 August 2008 (23 years). We discarded five of these earthquakes due to poor quality of the records. We, nevertheless, list them in Tables 1 and 2, and plot them in Fig. 1 for the sake of completeness.
Some of the events with PGA > 4 gal were not recorded because the accelerographs malfunctioned or because the memory was full. An unfortunate example is the earthquake of 21 September, 1985 (Mw 7.6), for which there is no recording at SCT. Thus, 31 and 32 as the number of earthquakes during which the ground motion exceeded 4 gal at SCT and CDAO is a lower bound.
To test the performance of the proposed device, we investigate the alert time, TA, that it would provide based on recorded data. We define TA as the difference between the arrival time of the intense part of the ground motion, TI, and the trigger time,TI As mentioned above, it is convenient to measure TI as the time when the energy in the accelerogram, A(t), reaches 5% of the total energy, where total energy is given by
We define TT as the time when A(t) on any component reaches 4 gal. Thus, the alert time, TA, is given by
Generally, TT values for the three components differ. We take the smallest of the three values of TT to compute TA, and we only compute TA for the horizontal components. Fig. 2 illustrates the concepts defined in the previous paragraph. With one exception, TT is larger on the Z component than on the horizontal components (i.e., the 4gal threshold is exceeded later in the vertical component than in the horizontal ones). The exception is the inslab earthquake of Copalillo (21 July, 2000, Mw 5.9), during which the 4gal threshold is first reached in the P wave recorded on the Z component.
As discussed earlier, A(t) = 4 gal may occur at the beginning of the intense motion. In this case TA will be zero. It is also possible that A(t)= 4 gal occurs after the arrival of the intense part of the motion. In this case TA will be negative. For reasons mentioned above, the alert time, TA, will depend on AT and preevent memory. To test the sensitivity of the results we analyzed those accelerograms which were recorded with small AT and large preevent memory. These records were reprocessed at larger values of AT and smaller values of preevent time. Our tests show that the effect of different values of AT and preevent memory on the alert time is small, typically less than 2 sec. Thus, we present results based on the original accelerograms. Fig. 3 illustrates the accelerograms of the Michoacan and Copalillo earthquakes at SCT along with trigger time, TT, and the arrival time of the intense part of the ground motion (5% of the total energy), Tl The figure also shows the time at which the energy is 50% of the total.
In Fig. 4 the alert time, TA , is plotted as a function of PGA for both SCT and CDAO. The two horizontal components are shown separately for each station. In the figure interplate, inslab, continental crustal (near the coast of Guerrero), and local earthquakes in the Valley of Mexico are shown by different symbols. The plot suggests that all data may be grouped together. The figure shows that greater the eventual PGA, larger is the expected TA. This implies that the amplitudes of Swave group and Pwave group scale together. If for an earthquake the amplitude of Pwave group is large, then it is detected early and, hence, TA is large. For the same earthquake the amplitude of Swave group will also be large and, hence, the relationship between TA and PGA seen in the figure.
Damage to buildings in the lakebed zone of the Valley of Mexico from large, coastal, interplate earthquakes, probably, results from the longduration, nearly harmonic ground motion. The collapse occurs because of fatigue from longduration cyclic loading. We define TD as the difference between the time when energy is 50% of the total, T50, and the smallest of the three values of trigger time, TT (Fig. 2). For coastal earthquakes T50 may be close to the time of the PGA (see Fig. 3 for the Michoacan earthquake recording at SCT). In Fig. 5 TD is plotted as a function of PGA. As expected, TD is greater than TA. There is large dispersion in the data due to the variability of the signal during the intense part of the motion, which, in turn, is a consequence of source and path effects. Note in Fig. 3, for example, the difference in the intense part of the ground motion during the interplate, Michoacan earthquake (H = 17 km, Δ = 402 km) and the inslab, Copalillo earthquake (H = 50 km; Δ = 142 km) (events 1 and 31 in Table 1 and Fig. 1).
Discussion and conclusions
To evaluate the performance of the proposed earthquake alert device we have defined the alert time, TA, as the difference between the time of beginning of the intense part of the ground motion (measured as the time when the energy becomes 5% of the total) and the trigger time. We assume that the trigger level is 4 gal. An analysis based on higher threshold levels would have resulted in shorter alert times and vice versa (however, a lower threshold level is not desirable for increasing alert time since it would also increase the number of false alerts, as it is discussed below).
As one would have expected, the larger the PGA, the greater is the alert time, TA. To understand the performance of the device, let us assume that a unit was installed on the ground floor of a one to threestorey school building near SCT before September 1985. (Several hundred one to threestorey school buildings were damaged and many collapsed during the Michoacan earthquake. Luckily, there were few casualties because the earthquake occurred at 7:18 am, when the schools were still not open. The damaged schools have been retrofitted.) Let us also assume that preestablished safety procedures were being followed during each earthquake which triggered the device. With a 4gal threshold, the device near SCT would have triggered at least 31 times due to earthquakes in the period September 1985 August 2008 (23 years). TA would have been negative for 7 out of the 31 events (i.e., the beginning of the intense part of the motion would have already passed at the time of the alert). In these cases, the alert and the associated safety measures would have been in vain (though they may have served as earthquake drills). Even a large, positive alert time for these earthquakes (say, from a system such as the SAS) would have been of little benefit since the PGA would have been less than 10 gal. For 16 earthquakes, TA would have been between 0 and 15 sec. For these earthquakes, the alert device could have been useful in shutting critical facilities, but it would have been difficult to reach a previouslyassigned safe place within the alert time. During 5 events TA would have been greater than 15 sec and the PGA would have exceeded 20 gal (earthquakes of 19/09/1985, Michoacan, Mw 8.0; 14/09/1995, Copala, Mw 7.3; 15/06/1999, Tehuacan,Mww 6.9; 30/09/1999, Oaxaca, Mw 7.4; and 21/07/2000, Copalillo, Mww 5.9). For a small building and persons trained by previous drills, TA > 15 sec probably would have been sufficient to reach a safe area. For the truly destructive 1985 earthquake of Michoacan, the alert time would have been about 40 sec. We recall, however, that TA refers to the beginning of the intense motion. Thus, even if persons had reached a safe place only after TA, they would have avoided being exposed to the complete train of highamplitude strong motion (Fig. 3) in an unsafe place. Similar conclusions are reached from the data recorded at CDAO.
In Table 3 we compile alert time, TA, at other lakebed sites in the Valley of Mexico where the Michoacan earthquakes of 19 and 21 September, 1985 were recorded. TA would have been in the range of 38 44 sec and 30 36 sec during these two earthquakes except at CDAF site. During these earthquakes CDAF had a trigger level of 8 gal, and a preevent memory of only 4 sec. As a consequence, the records from CDAF start after the arrival of the intense part of the groundmotion (4gal threshold is exceeded in the first 1.5 sec of the horizontal components for both events), thus providing underestimated TA values. Thus, Table 3 gives us some assurance that the results from SCT and CDAO may, generally, be valid for all lakebed sites in the Valley of Mexico.
An issue of concern is that the device may trigger due to cultural noise. A check of the maintenance log book of SCT shows eight triggers at a threshold level of 4 gal between 19 April 2006 and 18 August 2007. Of these eight triggers, three were caused by earthquakes and the rest originated from cultural noise. Clearly, it would be desirable to equip the alert devices with filters that discriminate between earthquake ground motion and the motion caused by cultural noise, thus reducing false alerts. A higher threshold level is not desirable for reducing false alerts since it would also decrease the alert time.
The proposed device seems well suited to stop critical services and equipments during an earthquake before the arrival of intense motion. It would, however, be useful in mitigating injuries and loss of lives only to the extent that an alert time of 15 to 45 sec can be fruitfully utilized to implement safety measures. This issue is critical and must be given careful consideration before any deployment. We reiterate that the device will not be useful for local earthquakes, as they will result in negative or small alert time. Even a repeat of 1912, Acambay earthquake (M 7.0, Δ 100 km) will, at most, give an alert time of 12 sec.
It is instructive to note that, had the SAS been in operation since 1985, it would have issued public alert for only two of the seven earthquakes listed in Table 2 which produced PGA > 20 gal at SCT, namely for the earthquakes of San Marcos (25/04/1989) and Copala (14/09/1995). The SAS, in fact, issued a public alert for the Copala earthquake with an alert time of 72 sec (EspinosaAranda et al., 1995) as compared to 20 sec for the alert device under scrutiny. It is possible that the SAS would have missed the Michoacan earthquakes. The reason, of course, is the limited spatial coverage of the SAS, as the sensors cover only the coast of Guerrero.
Table 4 compares the SAS and the earthquake alert device analyzed here. Both systems have merits and deficiencies. A detailed evaluation of the performance of the SAS and a strategy for its improvement is given in Iglesias et al. (2007). A thorough evaluation of any proposed alert system is essential before its implementation. This is because future changes are costly and, socially, confusing. We suggest that the proposed earthquake alert undergo a rigorous analysis before the authorities make a final decision. It should also be remembered that the SAS, in spite of current deficiencies (Iglesias et al., 2007), is already in place. Would public schools receive alerts from both systems? A careful thought should be given about the public confusion which will inevitably result from the operation of two alert systems.
We thank Roberto Meli and Luis Esteva for useful discussions. Ana Laura Ruiz kindly compiled some of the data used in the analysis. Fig. 1 was created with the GMT free software (Wessel and Smith, 1998). This research was partly supported by DGAPA, UNAM project IN110207. D. García was supported by Postdoctoral Fellowships Program from DGAPA, UNAM.
EspinosaAranda, J. M., A. Jiménez, G. Ibarrola, F. Alcantar, A. Aguilar, M. Hinostroza and S. Mal donado, 1995. Mexico City Seismic Alert System. Seism. Res. Lett. 66, 6, 4253. [ Links ]
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. Bull. Seism. Soc. Am. 97, 5, 17181729. [ Links ]
Wessel, P. and W. H. F. Smith, 1998. New, improved version of Generic Mapping Tools released. EOS 79, 47, 579. [ Links ]