Estimación del peligro sísmico debido a sismos interplaca e inslab y sus implicaciones en el diseño sísmico

Adrián David García Soto¹, Adrián Pozos-Estrada¹, Hanping Hong², y Roberto Gómez Martínez¹

¹ Instituto de Ingeniería de la UNAM, Ciudad Universitaria, CP 04510, México DF. AGarciaSo@iingen.unam.mx, APozosE@iingen.unam.mx, RGomezM@iingen.unam.mx

² Departamento de Ingeniería Civil y Ambiental, The University of Western Ontario, London, Ontario, Canadá, N6A, 5B9. hongh@eng.uwo.ca

Resumen

La ocurrencia de grandes sismos causa excitaciones muy intensas del terreno que pueden provocar el colapso de edificios y puentes. Para tomar en cuenta la incertidumbre en esta ocurrencia y en las solicitaciones sísmicas de diseño, la evaluación probabilística del peligro sísmico es empleada para desarrollar espectros de peligro uniforme (EPU) y mapas de contorno de peligro sísmico. Aunque existe mucha información en relación al cálculo del peligro sísmico en México y su inclusión en reglamentos, no hay en la literatura un reporte detallado, que incluya un modelo de peligro sísmico útil para obtener los EPU, para desarrollar mapas de contorno, y para determinar respuestas inelásticas. Además, ya que la demanda de ductilidad de desplazamiento puede diferir para sismos interplaca e inslab, esta diferencia debería incorporarse en la estimación de la confiabilidad estructural. En este estudio se integra un modelo de peligro sísmico útil para estimar los EPU y para desarrollar mapas de contorno para una parte de México; se calcula la contribución al peligro sísmico de cada tipo de sismo; y se estima el nivel de carga sísmica requerido para diseño, considerando las diferencias en la demanda de ductilidad causada por sismos interplaca e inslab.

Palabras Clave: peligro sísmico; mapas de contorno; respuesta inelástica; confiabilidad; diseño sísmico.

Abstract

Large earthquakes cause intensive ground excitations and possibly results in collapse of buildings and bridges. To take into account the uncertainty in earthquake occurrence and seismic load effect in design codes, the probabilistic seismic hazard assessment approach is employed to develop the uniform hazard spectra (UHS) and seismic hazard contour maps. Although there is extensive information concerning Mexican seismic hazard assessment and code making, a detailed report on a hazard model that can be used to evaluate the UHS, to develop the seismic hazard contour maps and to assess inelastic responses is lacking. Furthermore, since the displacement ductility demand could differ for interplate and inslab earthquakes, this difference should be incorporated in estimating the reliability of structures. In the present study a workable seismic hazard model that can be use to estimate UHS and to develop seismic hazard contour maps for part of Mexico is assembled; the contribution of each individual earthquake type to the overall seismic hazard is assessed; and the seismic design load level that is required considering the differences in the ductility demand caused by interplate and inslab earthquakes is estimated.

Key Words: seismic hazard; contour maps; inelastic response; reliability; seismic design.

INTRODUCCIÓN

La ocurrencia de grandes sismos, que son poco frecuentes y aleatorios, causa excitaciones muy intensas del terreno que pueden provocar el colapso de edificios y puentes. La medición del movimiento del terreno está usualmente basada en la aceleración máxima del terreno o en la seudoaceleración (SA). Para tomar en cuenta la incertidumbre en la ocurrencia de un sismo y en la medición del movimiento del terreno, se emplea la evaluación probabilística del peligro sísmico (EPPS) (Cornell 1968, Esteva 1968, McGuire 2004). La EPPS incorpora información de las zonas sísmicas, de los modelos de ocurrencia, de las relaciones magnitud-recurrencia, y de las ecuaciones de predicción del movimiento del terreno (o leyes de atenuación) para estimar el peligro sísmico (e.g., espectros sísmicos de respuesta); también considera la incertidumbre del modelo a través de la actualización Bayesiana. La EPPS es usada para desarrollar espectros de peligro uniforme (EPU) y mapas de contorno de peligro sísmico. Ejemplos de éstos incluyen el cálculo del peligro sísmico en Estados Unidos (Frankel 1995, Frankel et al. 1996) y en Canadá (Adams and Halchuk 2003), donde los modelos de peligro sísmico adoptados (zonas sísmicas, modelos de ocurrencia, relaciones magnitud-recurrencia y leyes de atenuación) se describen detalladamente. Aunque existe una vasta cantidad de publicaciones y actividades de investigación en relación al cálculo del peligro sísmico en México y su inclusión en reglamentos de diseño (Rosenblueth y Gómez 1991, Zúñiga et al. 1997, Ordaz y Reyes 1999, Reyes et al. 2002, Esteva et al. 2002, García 2006, Ordaz et al. 2007, Bojórquez et al. 2008, y Tena-Colunga et al. 2008), no está disponible en la literatura un reporte público y detallado, que contenga la información de un modelo de peligro sísmico en México que pueda ser empleado directamente para obtener los EPU, para desarrollar mapas de contorno de peligro sísmico, y para determinar respuestas inelásticas. Los códigos de diseño de edificios y puentes tienen la intención de alcanzar niveles de confiabilidad consistentes en dichas estructuras o la eficiencia económica de éstas, por lo menos este es actualmente el caso para las mismas categorías estructurales. Con el desarrollo de un código para México, Rosenblueth y Esteva (1972) fueron los pioneros en sentar las bases probabilísticas para seleccionar el nivel de carga sísmica y para la elaboración de reglamentos de diseño. Si bien se dispone de nuevos métodos para evaluar la confiabilidad estructural, así como de nuevas teorías que pueden utilizarse para la toma de decisiones bajo incertidumbre, las bases comúnmente aceptadas y reconocidas para el desarrollo de reglamentos, son las que se sustentan en alcanzar una confiabilidad consistente, y en la teoría del mínimo costo esperado para el ciclo de vida útil (Ellingwood et al. 1980, Madsen et al. 1986, Rackwitz 2000, Ellingwood 2001, Goda y Hong 2006). No obstante, el nivel de confiabilidad implícito para el diseño contra excitaciones sísmicas mediante el uso de un reglamento, no está siempre dilucidado, aún a nivel nacional. Un ejemplo es el actual código de diseño de puentes (SCT 2001), que provee mapas de contorno de la aceleración máxima del terreno y la forma del espectro de diseño, los cuales provienen de los estipulados en el Manual de Obras Civiles (MOC 1993).

La suposición de que los edificios o puentes pueden ser modelados como sistemas de un solo grado de libertad, sirve de base para establecer la fuerza sísmica de diseño en los reglamentos, ya que el nivel de carga sísmica está basado en la SA. En algunos códigos actuales, la fuerza sísmica está definida por el EPU, mientras que en otros se define como una proporción de un espectro de diseño basado en formas paramétricas predeterminadas, que son función de la aceleración máxima del terreno y otros parámetros. El reglamento canadiense para edificios (NBCC 2005) es un ejemplo donde la primera alternativa mencionada anteriormente es adoptada, mientras que la segunda alternativa es considerada en el código canadiense de diseño de puentes (CSA-S6-06 2006), así como en el código AASHTO LRFD (AASHTO 2007), y en el Manual de Obras Civiles de México (MOC 2008). Debido a que sitios dentro de una región grande podrían ser impactados por diferentes tipos de sismo con características de excitación distintas, el uso de un espectro de diseño basado en una forma paramétrica predeterminada, podría no garantizar un nivel de confiabilidad consistente o eficiencia económica. Aunado a esto, la demanda de ductilidad de desplazamiento puede diferir para diferentes tipos de sismo, como es el caso de sismos interplaca e inslab (Hong et al. 2010), y esta diferencia debería incorporarse en la estimación de la confiabilidad estructural.

Cabe señalar que los sismos aquí indicados como interplaca, corresponden a los sismos de subducción, es decir, a los sismos que presentan mecanismos de falla inversa debido a la liberación de energía acumulada entre las placas tectónicas a lo largo de los límites convergentes de las mismas. El otro tipo de sismos, aquí referidos como sismos inslab, corresponden a los sismos de falla normal, de profundidad intermedia, que se localizan dentro de las placas oceánicas que subducen bajo el continente, y que ocurren cuando la roca se rompe por las grandes deformaciones que producen los esfuerzos de tensión aplicados sobre el plano de la roca.

Los principales objetivos de este estudio son: integrar un modelo de peligro sísmico que facilite la estimación de los EPU y el desarrollo de mapas de contorno de peligro sísmico, para una parte de México con posibilidades de ser significativamente afectada por sismos interplaca e inslab; calcular la contribución al peligro sísmico de cada uno de estos tipos de sismos en diferentes sitios; estimar el nivel de carga sísmica requerido para diseño, considerando las diferencias en la demanda de ductilidad de desplazamiento causada por los tipos de sismos en cuestión. Para los análisis, no son considerados sismos que no correspondan a los de interplaca e inslab, tales como los sismos "crustal" (Zuñiga et al. 1997, Ordaz et al. 2007). Tampoco se ha incorporado la contribución de la actividad sísmica en los Estados Unidos que afecta el peligro sísmico de México. Este estudio se enfoca sólo en los efectos de los sismos interplaca e inslab en sitios de suelo firme (roca). Otros tipos de sismos como los de corteza, intraplaca o "crustal", que pudieran regir el diseño en periodos estructurales cortos, están fuera del alcance de este estudio.

Para alcanzar los objetivos referidos, el cuerpo de este trabajo se conforma de dos partes principales. Por una parte se aborda el cálculo del peligro sísmico describiendo las zonas sismogenéticas, las relaciones de magnitud-recurrencia, y las leyes de atenuación empleadas para tal efecto. Con la información descrita, es posible emplear el procedimiento de la EPPS desarrollada por Cornell (1968) y Esteva (1968), y se indica en esta parte que en lugar de utilizar integración numérica, se adoptó un procedimiento basado en simulación, el cual también se describe. Lo anterior conduce a la obtención de mapas de contorno de peligro sísmico, así como espectros de peligro uniforme considerando los dos tipos de sismos empleados en este trabajo. La segunda parte principal de la que consta el cuerpo del presente artículo, hace uso de los resultados obtenidos en la primera parte, para abordar la respuesta inelástica de estructuras modeladas como un sistema de un grado de libertad y el cálculo de su nivel de confiabilidad, así como su implicación desde el punto de vista reglamentario. En esta segunda parte se utiliza información desarrolladla en trabajos previos, que al complementarse con lo desarrollado en este trabajo, culmina en la consecución de los objetivos planteados.

ANÁLISIS DEL PELIGRO SÍSMICO

Metodología y modelo de peligro sísmico

El proceso de ocurrencia de un sismo es comúnmente modelado como Poisoniano, aunque también se ha considerado como un proceso no Poisoniano. En el primer caso, el tiempo de llegada entre dos temblores sucesivos se asume como una variable con distribución exponencial con un parámetro determinado por la tasa de ocurrencia sísmica, mientras que para el segundo, dicho tiempo entre eventos es usualmente considerado una variable del tipo lognormal, gama, o Gaussiana inversa (Rosenblueth y Jara 1991). Más aún, se asume frecuentemente que la ocurrencia de un sismo es estadísticamente homogénea en magnitud y en ubicación dentro de una zona sísmica. Para definir las características del peligro sísmico en una parte de México, consideramos que son adecuadas las zonas sísmicas proporcionadas por Ordaz y Reyes (1999) cerca de la Costa del Pacífico Mexicano. Estas zonas se separan en dos grupos como se aprecia en las figuras 1a y 1b, y cuyos parámetros se enlistan en la tabla 1. El primer grupo corresponde a sismos pequeños y moderados con magnitudes menores a 7, y el segundo grupo a sismos con magnitudes mayores a 7. Cada zona sísmica está definida por un polígono; sus contornos fueron digitalizados a partir de las figuras graficadas por Ordaz y Reyes (1999). La ocurrencia de sismos interplaca (i.e., eventos en la frontera de la placa) se modela dentro de estas zonas sísmicas. Se considera que los sismos inslab, que ocurren dentro de la placa oceánica subducida debajo de la placa continental, tienen origen en las grandes zonas sísmicas mostradas en la figura 1c y en la tabla 1 (Jaimes y Reinoso 2006).

Para las zonas sísmicas mostradas en la figura 1a con M_W < 7, la relación magnitud-recurrencia, λ(M_W), está definida por,

y para las zonas sísmicas mostradas en la figura 1b con M_W > 7, por,

donde M_W representa el momento de magnitud, M_L = 4.5, M_U = 7, m_M = 7.5, s_M = 0.3, Φ( ) denota la función de distribución de probabilidad normal estándar, y λ₀, β, y λ(7) son parámetros de modelación. Los valores de los parámetros incluyendo λ₀, β, y λ(7) son tomados directamente del estudio de Ordaz y Reyes (1999) y listados en la tabla 1; λ(7), en lugar de λ₀, se utiliza como la tasa de ocurrencia sísmica para sismos con magnitud mayor a 7.0. Para los sismos inslab dentro de las zonas sísmicas mostradas en la figura 1c, se considera que la ec. 1 también se puede emplear, pero con M_L = 4.5, y M_U = 7.9 para la región denominada Profundidad Intermedia Oeste, y M_U = 7.88 para la región denominada Profundidad Intermedia Centro (Jaimes y Reinoso 2006).

La ec. 1 implica que la función de distribución de probabilidad de la magnitud M_w, F_M(m), se puede expresar como,

mientras que la función de distribución de probabilidad asociada con la ec. 2 está dada por,

De las ecs. 1 y 2 también se infiere que el tiempo de llegada de un sismo está distribuido exponencialmente, su media es igual a 1/λ₀ y 1/λ(7) para las zonas símicas correspondientes, siempre que la ocurrencia esté representada por un proceso Poisoniano. En el presente estudio se adopta el modelo de ocurrencia Poisoniano. Esta suposición ha sido ampliamente adoptada en la literatura (ver las referencias enlistadas en Pozos-Estrada et al. 2008). Lo anterior conlleva a la utilización de caracterizaciones simples del peligro sísmico, independientes del tiempo.

Si bien existe un estudio sobre la traza de la interface de las placas continental y subducida (Pardo y Suárez 1995), debido a que existe incertidumbre acerca de la profundidad de los sismos, y por lo tanto acerca de la superficie que define la profundidad de los mismos, en este estudio se consideraron profundidades deterministas de 10.45 km para sismos interplaca, y de 64.56 km para sismos inslab, las que representan, respectivamente, la profundidad promedio de los sismos asociados con los registros utilizados en el desarrollo de las leyes de atenuación descritas en las líneas subsecuentes. Las formas funcionales adoptadas para las relaciones de atenuación son,

para la estación de CU (localizada en Ciudad Universitaria), desarrollada para sismos interplaca (Reyes et al. 2002),

para sismos interplaca (García 2006), y

En la ec. 5, R (km) es la distancia más cercana a la superficie de falla; y c_i, i = 1,...,5, representa coeficientes de regresión. En la ec. 6, R (km) es la distancia más cercana a la superficie de falla para eventos con M_w > 6.0, o la distancia hipocentral para el resto; H (km) es la profundidad focal; c₄ = 1.82 - 0.16M_w; y c₁, c₂, c₃, y c₇ son coeficientes de regresión. En la ec. 7, c_i, i = 1,2,3,4,5, también representa coeficientes de regresión; c₄ =1; R=(R_cld²+Δ₀)^1/2; R_cld (km) es la distancia más cercana a la superficie de falla para eventos con M_w > 6.5, o la distancia hipocentral para el resto; H (km) es la profundidad focal; y Δ₀=0.0075x10^0.507Mw es un término de saturación de la fuente definido por Atkinson y Boore (2003).

Los coeficientes, c_i, obtenidos por Reyes et al. (2002), García et al. (2005), y García (2006) para las ecs. 5 a 7, están basados en la media cuadrática; nosotros preferimos usar las leyes de atenuación desarrolladas con base en la media geométrica (definida como la raíz cuadrada del producto de las seudoaceleraciones para cada componente), ya que éstas derivan en la desviación estándar en las ecuaciones de predicción de movimiento del terreno para una orientación aleatoria. Los coeficientes que provienen de emplear la media geométrica para las ecs. 5 a 7, se pueden encontrar en las publicaciones de Pozos-Estrada et al. (2008) y Hong et al. (2009) y se muestran en la tabla 2. Estos coeficientes se obtuvieron recurriendo a la misma serie de registros seleccionados por Reyes et al. (2002), García et al. (2005), y García (2006), pero empleando el algoritmo de regresión de Joyner y Boore (1993). Nótese que una corrección de línea base (Boore 2005) fue aplicada a los registros considerados por García et al. (2005) y García (2006); nótese además que menos de 10 registros fueron excluidos, ya que éstos conducían a respuestas poco realistas (Hong et al. 2010). Los valores de Y que se predicen utilizando los coeficientes mostrados en la tabla 2, son comparables con aquellos cuya predicción se obtiene mediante los coeficientes que se basan en la media cuadrática.

Los coeficientes de la ley de atenuación para la estación localizada en Ciudad Universitaria en la Cd. de México (suelo firme) fueron desarrollados para sismos interplaca (subducción) cuyos registros corresponden a momentos de magnitud (M_w) entre 6.1 y 8.1, y distancias epicentrales comprendidas entre 280 y 466 km (Reyes et al. 2002). En lo referente a los coeficientes para la ley de atenuación para sismos interplaca (subducción), estos fueron desarrollados para sitios sobre roca firme comprendidos entre la costa y el Cinturón Volcánico Mexicano; el número de eventos sísmicos es 40, con momentos de magnitud (M_w) entre 5 y 8, y profundidad focal (H) entre 8 y 29 km, y se consideró una distancia más próxima al plano de ruptura de hasta 400 km (García 2006). Finalmente, los coeficientes para la ley de atenuación para sismos inslab (falla normal, profundidad intermedia), fueron desarrollados para sitios en México Central sobre roca firme; el número de eventos sísmicos es 16 (registrados a distancias locales y regionales de hasta 400 km), con momentos de magnitud (M_w) entre 5.2 y 7.4, y profundidad focal (H) entre 35 y 138 km (García et al. 2005). Nótese que a diferencia de la ley de atenuación propuesta por Reyes et al. (2002), las relaciones funcionales propuestas por García (2006) y por García et al. (2005), incluyen la profundidad focal H, debido a que en esos estudios se encontró dependencia entre la aceleración máxima del terreno (o la seudoaceleración) y la profundidad focal.

Debe enfatizarse que las ecs. 6 y 7 no rigen para sitios dentro del Valle de México ni dentro del Cinturón Volcánico Mexicano, debido a las características particulares de propagación a lo largo del Cinturón Volcánico Mexicano y de respuesta de los sitios dentro del Valle de México (García 2006). La ec. 5 es válida para M_W ≤ 8.1; en caso de eventos con M_W > 8.1, sus magnitudes son consideradas iguales a 8.1 (Reyes, comunicación personal 2007). La misma condición se considera para la ec. 6.

Procedimiento de análisis

donde n_s denota el número total de zonas sísmicas; λ_Mmin,i es la relación de ocurrencia de los eventos sísmicos con M_W ≥ M_min para la i-ésima zona sísmica; f_Mw,R,i(m_W, r) denota la función de densidad de probabilidad conjunta de M_W y R; Ω_i denota el dominio de M_W y R para la i-ésima zona sísmica; F_Y|Mw,R,i(y) es la función de densidad de probabilidad conjunta de Y condicionada a M_W, R, y la i-ésima zona sísmica (la cual puede ser derivada de las características de la ley de atenuación aplicable); X_i representa a las variables aleatorias epistémicas (e.g., incertidumbre en los parámetros de modelación que se emplean para definir los modelos de la ley de atenuación y zonas sísmicas), para la i-ésima zona sísmica cuya función de densidad de probabilidad conjunta está definida por f_Xi(x); y Ω_Xi denota el dominio de X_i para la i-ésima zona sísmica. En lugar de utilizar integración numérica como fue hecho por McGuire (1976), la ec. 8 puede ser evaluada mediante técnicas de simulación para modelos de ocurrencia tanto Poisonianos como no Poisonianos.

En el presente estudio se adopta un procedimiento basado en simulación, siguiendo los pasos descritos a continuación, para la EPPS y considerando un período de observación igual a T_T (Hong et al. 2006):

1) Muestrear las variables aleatorias epistémicas en X_j para la j-ésima zona símica, j = 1,···,n_s;

2) Muestrear el tiempo de llegada entre sismos τ_j de acuerdo al modelo de ocurrencia adoptado y fijar t_j = τ_j, j = 1,···,n_s; encontrar la i-ésima zona sísmica cuyo t_i es igual a min(t₁,···,t_n);

3) Muestrear M_W (ver ecs. 3 y 4) y la localización del sismo, así como el término de error de la ley de atenuación pertinente (ver ecs. 5 a 7); medir el movimiento del terreno en términos de Y, y; muestrear τ_i de acuerdo a su modelo de ocurrencia y calcular t_i = t_i + τ_i;

4) Repetir el paso 3) si min(t₁,···,t_n) es menor que T_T.

Al repetir los pasos 1) a 4) n_E ciclos, se obtienen muestras para n_E×T_T años que se pueden subdividir en varios intervalos de T_d años, los que están igualmente espaciados sin sobreponerse. Los máximos efectos sísmicos que corresponden a cada intervalo de tiempo, pueden utilizarse para formar una secuencia de efectos máximos del movimiento del terreno para un periodo de T_d años. La incertidumbre epistémica no se considera en los análisis numéricos presentados en las secciones subsecuentes, debido a la carencia de diferentes (pero igualmente creíbles) modelos de peligro sísmico.

Peligro sísmico y mapas de contorno

Los resultados obtenidos se presentan en las figuras 3, 4 y 5 para los casos que consideran únicamente sismos interplaca, únicamente sismos inslab, y ambos tipos de sismos, respectivamente. Para ahorrar espacio, sólo se muestran los mapas de contorno de la aceleración máxima del terreno (T_n= 0 s) y la SA para T_n= 0.5 y 2 s, aunque también se calcularon los resultados para otros valores de T_n como ayuda para realizar algunas observaciones que se presentan más adelante. Si bien los resultados para el valle de la Ciudad de México y para el Cinturón Volcánico Mexicano son incluidos por complementariedad, éstos no deben usarse, ya que como se mencionó anteriormente las ecs. 6 y 7 no aplican para estos sitios (García 2006).

La figura 3 muestra mapas de contorno de peligro símico en términos de la media m_Y, cov v_Y, y valores de la aceleración máxima del terreno y la SA, Y, para un período de retorno igual a 2475 años, y₂₇₄₅. Nótese que el mapa con los valores para un periodo de retorno seleccionado, no proporciona una visión completa del peligro sísmico, mientras que los mapas de la media y el coeficiente de variación, podrían emplearse para proporcionar una caracterización probabilística completa del peligro sísmico, siempre que la distribución de probabilidad sea asignada. Esto último es especialmente importante para estimar la confiabilidad de estructuras sometidas a carga sísmica. Como se esperaba, al observar la figura 3, los mapas indican que, m_Y, y y₂₄₇₅ a lo largo y cerca de la costa del Pacífico Mexicano son los más altos para sismos interplaca

Los valores de v_Y varían aproximadamente entre 2 y 3 para la mayoría de los sitios; estos coeficientes de variación son significantemente mayores que aquéllos asociados con la incertidumbre en las propiedades de los materiales y variables geométricas en los sistemas estructurales, como puede consultarse en el estudio de Ellingwood et al. (1980). También se observan valores de v_Y mayores a 3.0, pero sólo para sitios asociados con valores pequeños de m_Y. El rango de valores de v_Y se reduce cuando T_n se incrementa. Debido a que v_Y no es constante para la aceleración máxima del terreno o la SA (para un valor específico de T_n) en la región considerada, las tendencias en y₂₄₇₅ siguen sólo aproximadamente aquéllas de la media.

Para apreciar mejor el peligro sísmico en función de T_n para el caso considerado en la figura 3 (únicamente sismos interplaca), en la figura 6 se muestran las curvas de los EPU para un periodo de retorno de 2475 años para cuatro sitios seleccionados, los cuales se denotan como Sitios I, II, III y IV, con latitudes y longitudes iguales a (19.5^°, -101^°), (17^°, -100^°), (17^°, -96.5^°) y (16.5^°, -95^°), correspondientes aproximadamente a Morelia, Acapulco, Oaxaca y Salina Cruz, respectivamente (ver figura 1a).

Los resultados indican que para sismos interplaca, conforme el sitio se localiza más lejos de la costa, el peligro sísmico disminuye, y que el peligro sísmico para el Sitio II, el cual es el más cercano a la costa, es significativamente más grande que para los otros sitios. Resultados para otros períodos de retorno, no presentados aquí, exhiben tendencias similares.

Si bien el análisis de la desagregación basado en la magnitud, la distancia, y el error de modelación (i.e., ε en las ecs. 5 a 7) está fuera de los alcances de este estudio, los resultados mostrados en la figura 6, de hecho, representan la desagregación del peligro sísmico basado en dos diferentes tipos de sismos. Tal desagregación es de interés, ya que los datos estadísticos de la demanda de ductilidad de desplazamiento pueden diferir y afectar la confiabilidad estructural (ver la sección RESPUESTA INELÁSTICA Y CONFIABILIDAD).

Como era de esperarse, para sitios cercanos al centro de las zonas que originan sismos inslab y mostradas en la figura 1c, se tienen los mayores valores de m_Y y y₂₄₇₅, como se aprecia en la figura 4. En contraste con los resultados mostrados en la figura 3, los v_Y para sismos inslab varían de 5 a 15 para sitios con peligro símico significativo, valores que son mucho más grandes que para el caso de sismos interplaca. El valor de y₂₄₇₅ sobre las zonas sísmicas es relativamente consistente y similar a, o mayor que, aquéllos observados para los sismos interplaca si T_n es aproximadamente menor a 0.5 (s). Esto implica que el peligro sísmico debido a sismos inslab puede regir el diseño estructural para estructuras rígidas, una observación que ya había sido advertida por Jaimes y Reinoso (2006). Sin embargo, conforme T_n es aproximadamente mayor que 0.5 (s), los valores de y₂₄₇₅, en sitios sobre las zonas sísmicas consideradas, son comparables a aquéllos debidos a sismos interplaca. Para investigar aún más la contribución relativa de los sismos inslab al peligro sísmico de la región, el EPU para un periodo de retorno igual a 2475 años se grafica en la figura 6 para los sitios mencionados anteriormente, y se compara con aquéllos debidos a sismos interplaca.

Los mapas de contorno mostrados en la figura 5, asociados a sismos inslab e interplaca, exhiben peligros más grandes que aquéllos mostrados en las figuras 3 y 4, lo que también era de esperarse. Los valores de v_Y se reducen significativamente, si se comparan con aquéllos causados por los sismos inslab, y se incrementan, si la comparación se realiza con los coeficientes de variación causados por sismos interplaca. En general, los v_Y para sismos interplaca e inslab varían de 4 a 7 para sitios con peligro sísmico considerable y para T_n menor que aproximadamente 0.5 (s); el rango de valores de v_Y se reduce para T_n mayor que 0.5 (s). En particular, si T_n es mayor que aproximadamente 1.0 (s), v_Y fluctúa entre 2 y 3 para sitios con peligro sísmico importante. Una vez más, los EPU debidos conjuntamente a sismos interplaca e inslab se muestran en la figura 6 con fines comparativos para los sitios seleccionados. Los EPU debidos a ambos tipos de sismo, son una envolvente de los debidos a un particular tipo de sismo (interplaca o inslab); la forma de los EPU depende de si el peligro es dominado por sismos interplaca o inslab.

Con fines comparativos, también se incluye en la figura 6 el EPU para la estación CU (localizada en la UNAM en la Ciudad de México), el cual se calculó utilizando las zonas sísmicas mostradas en la tabla 1 para sismos interplaca y la ec. 5. La figura muestra que la forma del EPU para la estación de CU es significativamente diferente que para los sitios cercanos a la región de la costa del Pacífico Mexicano, y que las magnitudes del EPU son comparables a aquéllas para la zona costera solamente para estructuras con periodos de vibración largos. Cabe mencionar que Singh et al. (1995) señalaron las diferencias de los efectos sísmicos en el valle de México y en la zona de la costa del Pacífico para algunos sismos históricos individuales.

Para cuantificar aún más las contribuciones relativas de los sismos interplaca e inslab al peligro sísmico, se evalúan los cocientes de y2475 de los resultados mostrados en las figuras 3 ó 4 en relación a aquéllos mostrados en la figura 5. Tales cocientes se grafican en la figura 7, en la que se observa que la contribución de los sismos interplaca es más relevante cerca de la costa, como era de suponerse, y que en general, se vuelve más significativa para valores grandes de T_n.

Un aspecto más que merece estudiarse, es la distribución de probabilidad de la SA. Para ello se extrajeron muestras de las simulaciones de la SA máxima anual y se presentan en la figura 8 (como gráficas de probabilidad lognormal), para el caso de eventos interplaca e inslab contribuyendo conjuntamente al peligro sísmico. En todos los casos, las gráficas muestran que la SA puede ser modelada como una variable lognormal, por lo menos en el extremo superior de la distribución de probabilidad. La suposición de que Y es una variable aleatoria lognormal, es común en el análisis de confiabilidad de estructuras sujetas a excitaciones sísmicas.

RESPUESTA INELÁSTICA Y CONFIABILIDAD

Demanda de ductilidad de desplazamiento

Los edificios y puentes se diseñan para incurrir en deformaciones inelásticas debidas a grandes sismos. En el actual reglamento de puentes (SCT 2001) el coeficiente sísmico basal C_s es igual a a dividido por Q', mientras que en el reglamento de la Ciudad de México (NTCS 2004) se considera un término adicional debido a la sobreresistencia R_o; en dichos códigos, a es el coeficiente sísmico elástico de diseño y Q' es un factor de reducción que depende de T_n y del factor de comportamiento sísmico, Q, cuya cuantificación y significado físico no son del todo claros. De hecho, Rosenblueth y Gómez (1991) subrayaron que Q no es la capacidad de ductilidad de desplazamiento, y que por mucho equivaldría a alrededor del 40% de dicha capacidad. Por lo tanto, podría ser instructivo reemplazar Q' por un factor de reducción relacionado con la capacidad de ductilidad de desplazamiento R_μ, tal como fue considerado en el más reciente MOC (2008). No obstante, a diferencia del caso del MOC (2008), en este estudio se considera que el cociente definido por la media de la fuerza cortante basal entre la mínima fuerza cortante basal de diseño requerida, R_o, incluye el efecto del comportamiento del sistema, así como la diferencia entre la media de la resistencia del material y su resistencia de diseño, y la diferencia entre los tamaños de los miembros estructurales diseñados y los tamaños requeridos dictados por los requerimientos mínimos de diseño. Además, se considera que el uso de una forma constante del espectro elástico de diseño para a podría no ser apropiado para una región amplia de México, ya que tal espectro no garantiza una probabilidad de excedencia consistente para sitios localizados dentro de dicha región, la que podría ser afectada por diferentes tipos de sismo con características de excitación distintas como se estudió en la sección anterior.

Con base en lo arriba expuesto, y considerando que la estructura diseñada puede ser modelada como un sistema histerético bilineal de un grado de libertad cuyo diseño está regido por el criterio de resistencia y sujeto a la demanda sísmica S(T_n,ξ) (g), la resistencia de fluencia normalizada ϕ puede ser expresada como,

donde

y donde R_n denota el cociente de la capacidad de fluencia de la estructura diseñada entre la fuerza cortante basal mínima de diseño, y R_μ es un factor de reducción relacionado con la ductilidad. El factor de irregularidad y el de importancia de la estructura se ignoran en las ecs. 9.

Para un valor dado o estimado de ϕ, se dispone de los datos estadísticos de la demanda de ductilidad para los grupos de registros sísmicos aquí considerados y para distintos valores del cociente de la rigidez de postfluencia entre la rigidez inicial, γ (Pozos-Estrada et al. 2008, Hong et al. 2010), por lo que la media de la demanda de ductilidad de desplazamiento μ, denotada por m_μ, puede ser aproximada con,

donde los parámetros de modelación α₁ y b, también están disponibles en los estudios mencionados y como complemento se reproducen en la tabla 3. El cov de μ, v_μ, es una función decreciente de ϕ, excepción hecha para estructuras muy rígidas; en la mayoría de los casos v_μ varía aproximadamente entre 0.4 y 1.2 para ϕ entre 0.1 y 0.5. Gráficas ilustrando la variación de m_μ y v_μ se muestran en las figuras 9 y 10, respectivamente.

De acuerdo a los estudios referidos en las líneas anteriores, en términos generales, el modelo de probabilidad de Frechet puede ser adecuado para μ, excepto para T_n < 0.3, en cuyo caso se prefiere el modelo de probabilidad lognormal. Una de las observaciones principales de estos estudios, es que en general los datos estadísticos de la demanda de ductilidad de desplazamiento difieren para los sismos Mexicanos inslab e interplaca. Tal diferencia requiere ser incorporada en el análisis de confiabilidad y en la calibración de reglamentos de diseño.

Cálculo de confiabilidad e implicación en el reglamento de diseño

El análisis de confiabilidad se simplifica si se aproxima el comportamiento de las estructuras como sistemas bilineales histeréticos de un grado de libertad. La función de estado límite para el colapso incipiente, g_c, puede expresarse como,

donde μ_R denota la capacidad de ductilidad del sistema, y μ, la demanda de ductilidad de desplazamiento, se remplaza por μ(A) para resaltar que ésta depende o es función de A, A = [ϕ, γ, T_n, ξ]. La combinación de valores de ϕ menores que 1.0 y g_c menores que cero implica colapso incipiente.

Antes de llevar a cabo el estudio paramétrico de la probabilidad anual de que g_c ≤0, denotada por P_c, señalamos que R_n/R_o (que representa la capacidad (real) de fluencia entre la capacidad media del sistema) debe ser, en promedio, mayor que 1.0 para el sistema, y que la incertidumbre en R_n/R_o puede ignorarse dado que su impacto en la P_c estimada es usualmente despreciable. Siguiendo la práctica actual en otras jurisdicciones en el mundo y sin pérdida de generalidad, por el momento S_E(T_n,ξ) se iguala con el valor correspondiente a un periodo de retorno de 2475 años. μ_R puede modelarse como una variable aleatoria lognormal (Díaz-López y Esteva 1991) con un valor de cov, v_μR , variando de 0.5 a 1.0. Lo anterior es orientado por el cov de la ductilidad de miembros de acero que varía de 0.5 a 1.0 (Nakashima 1997).

Se asume que el factor de reducción relacionado con la capacidad de ductilidad de desplazamiento R_μ, es igual a la media de μ_R, m_μR . Ya que μ(A) difiere para sismos interplaca e inslab, P_c requiere estimarse sumando la probabilidad de colapso incipiente debido a sismos interplaca, más la probabilidad de colapso incipiente debido a sismos inslab.

Considérese que S_i(T_n, ξ) denota la SA originada por el i-ésimo tipo de sismo, i = 1 y 2 representa respectivamente los sismos interplaca e inslab. Para estimar P_c, se emplea la técnica de simulación (Hong y Hong 2007), donde P_c representa la suma de la probabilidad de colapso incipiente debido a cada tipo de sismo. Básicamente, se simula un valor de S_i(T_n, ξ) y se le utiliza en la ec. 9 para evaluar ϕ; se generan muestras de μ_R y μ(A) y se usan en la ec. 11 para evaluar g cuando ϕ es menor que 1.0. Nótese que los valores estadísticos de S_i(T_n, ξ) para un sitio específico, pueden obtenerse de resultados similares a los mostrados en las figuras 3 a 5 (4) (ver también la tabla 4 para unos cuantos sitios seleccionados).

Para los sitios considerados, mostrados en la figura 1a, incluyendo la estación CU en la Ciudad de México, y cuyos modelos probabilísticos de la SA se enlistan en la tabla 4, los valores estimados de P_c se muestran en la figura 11 para R_μ = m_μR , y combinaciones de R_n/R_o = 1.15, m_μR = 2 y 4, v_μR = 0.5 y 1.0, y T_n = 0.1, 0.2, 0.5, 1 y 2. Los resultados mostrados en la figura y otros no incluidos aquí para ahorrar espacio, sugieren que:

1) Para los Sitios I a IV, el uso consistente de valores de S_E(T_n,ξ) para un período de retorno dado, conlleva a su vez a valores consistentes de P_c para cocientes fijos de R_n/R_o. P_c depende de R_n/R_o como se esperaba, por lo tanto, si se desea aun más consistencia en P_c, deben hacerse recomendaciones de tal manera que R_n/R_o sea uniforme para un rango de casos de diseño.

2) P_c para la estación CU es marcadamente diferente que aquéllas para los Sitios I a IV. Esto puede explicarse al notar que la media de la demanda de ductilidad para la estación de CU difiere significativamente que aquélla para otros Sitios (ver figura 9).

3) Para la estación CU, P_c se incrementa a medida que la media o el cov de μ_R aumenta.

Esta tendencia en el incremento de P_c para la estación CU, en directa proporción a m_μR, se explica si notamos que, a medida que m_μR aumenta, la ϕ debida a los eventos sísmicos probablemente será reducida (ver ec. 9), y un significativo cov de la demanda de ductilidad se asocia con valores más pequeños de ϕ. Sin embargo, esta clara tendencia no se presenta para los Sitios I a IV, excepto para valores pequeños de T_n y diferentes m_μR. Esto se atribuye a las diferencias en las curvas de cov mostradas en las figuras 10b y 10c a aquéllas mostradas en la figura 10a. Por lo tanto, para aumentar todavía más la consistencia en la confiabilidad, un valor de R_μ dependiente del período (similar al caso para el Q implementado en el reglamento de la Ciudad de México (NTCS 2004) pero con diferentes parámetros) debería incorporarse. Sin embargo, este aspecto está fuera de los alcances de este estudio.

RESUMEN Y CONCLUSIONES

Se integró un modelo de riesgo sísmico y fue empleado para evaluar mapas de contorno de peligro uniforme para regiones cercanas a la costa del Pacífico Mexicano que son significativamente afectadas por sismos interplaca e inslab. De los resultados de los mapas de contorno se desprende que, dependiendo del sitio considerado, el peligro sísmico predominante para estructuras rígidas podría provenir principalmente de sismos interplaca, o de sismos inslab. El coeficiente de variación (cov) de la SA varía significativamente con el periodo de vibración, y sus valores están entre aproximadamente 2 a 10 si se consideran ambos tipos de sismos, interplaca e inslab. Este rango de valores de cov es mucho más grande que aquél para la estación CU, el cual es de alrededor de 1.2. Una comparación de los valores estadísticos de la demanda de ductilidad también indica que la media y el cov de la demanda de ductilidad para sitios cerca del área costera son más pequeños que aquéllos para la estación CU. Por lo tanto, los valores estadísticos de peligro sísmico de esta última, no deben aplicarse a la primera (o viceversa).

Al adoptar un formato de diseño que es similar al recomendado en el MOC (2008) y el NBCC (2005), se puede obtener un nivel de confiabilidad relativamente consistente si los factores de reducción son apropiadamente seleccionados y la SA (nominal) es seleccionada con base en un periodo de retorno de 2475 años.

AGRADECIMIENTOS

Se reconoce de manera especial el apoyo financiero del Consejo Nacional de Ciencia y Tecnología (CONACYT) de México, del Consejo de Investigación para las Ciencias Naturales e Ingeniería de Canadá (NSERC por sus iniciales en inglés), y del Instituto de Ingeniería de la UNAM. Agradecemos a K. Goda, Daniel García, Carlos Reyes, y Miguel Jaimes por la información proporcionada, y a J. A. Escobar, Orlando Díaz, y dos revisores anónimos por sus cuantiosos comentarios constructivos, sugerencias, y críticas.

REFERENCIAS

Adams, J, y S Halchuk (2003), "Fourth generation seismic hazard maps of Canada: values for over 650 Canadian localities intended for the 2005 National Building Code of Canada", Open-File 4459, Geological Survey of Canada, Ottawa, Ontario, Canada.         [ Links ]

AASHTO (2007), LRFD bridge design specifications, 4th edition, American Association of State Transportation Officials, Washington D. C.         [ Links ]

Atkinson, G M y D M Boore, (2003), "Empirical ground-motion relations for subduction-zone earthquakes and their application to cascadia and other regions", Bulletin of the Seismic Society of America, Vol. 93, pp. 1703-1729.         [ Links ]

Bojórquez, E, S E Ruiz y A Terán-Gilmore (2008), "Reliability-based evaluation of steel structures using energy concepts", Engineering Structures, Vol. 30, No. 6, pp. 1745-1759.         [ Links ]

Boore, D M, W B Joyner and T E Fumal (1997), "Equations for estimating horizontal response spectra and peak acceleration from western North America", Seismic Research Letters, No. 68, pp. 128-153. (for Erratum: see 2005, Seism Res. Lett., No. 76, pp. 368-369).         [ Links ]

Boore, D M (2005), "On pads and filters: processing strong-ground motion data", Bulletin of the Seismological Society of America, Vol. 95, No. 2, pp. 745-750.         [ Links ]

CAN/CSA-S6 (2006), Canadian Standards Association, Canadian Highway Bridge Design Code, Toronto.         [ Links ]

Cornell, C A (1968), "Engineering seismic risk analysis", Bulletin of the Seismic Society of America, Vol. 58, pp. 1583-1606.         [ Links ]

Díaz-López, O, y L Esteva (1991), "Calibrating simplified models for computation of seismic reliability of multistory frames", in Proceeding of ICASP6 (eds. Esteva, L, and S Ruiz.), Institute of Engineering, Mexico, D.F., Mexico, pp. 328-337.         [ Links ]

Ellingwood, B R (2001), "Earthquake risk assessment of building structures", Reliab. Engng. Syst. Safety, Vol. 74, pp. 251-262.         [ Links ]

Ellingwood, B R, T V Galambos, J G MacGregor y C A Cornell (1980), "Development of a probability based load criterion for American National Standard A58", Spec. Pub. No. 577, National Bureau of Standards, Washington, D. C.         [ Links ]

Esteva, L (1968), "Bases para la formulación de decisiones de diseño sísmico", Tesis doctoral, Universidad Nacional Autónoma de México, Ciudad de México, México.         [ Links ]

Esteva, L, O Díaz-López, J García-Pérez, G Sierra y E Ismael (2002), "Life-cycle optimization in the establishment of performance-acceptance parameters for seismic design", Structural Safety, Vol. 24, pp. 187-204.         [ Links ]

Frankel, A (1995), "Mapping seismic hazard in the central and eastern United States", Seismic Research Letters., Vol. 66, pp. 8-21.         [ Links ]

Frankel, A, C Mueller, T Barnhard, D Perkins, E V Leyendecker, N Dickman, S Hanson, y M Hopper (1996), "National seismic hazard maps", Open-File 96-532, U.S. Department of the Interior, U.S. Geological Survey, Denver, CO.         [ Links ]

García D, (2006), "Estimación de parámetros del movimiento fuerte del suelo para terremotos interplaca e intraslab en México Central", Tesis doctoral, Universidad Complutense de Madrid, Madrid, España.         [ Links ]

García, D, S K Singh, M Herraiz, M Ordaz y J P Pacheco (2005), "Inslab earthquakes of Central Mexico: Peak ground-motion parameters and response spectra", Bulletin of the Seismic Society of America, Vol. 95, pp. 2272-2282.         [ Links ]

García-Soto, A D, H P Hong y R Gómez (2012), "Effect of the orientation of records on displacement ductility demand", Canadian Journal of Civil Engineering, 39(4), 362-373.         [ Links ]

Goda, K, y H P Hong (2006), "Optimum seismic design considering risk attitude, societal tolerable risk level and life quality criterion", Journal of Structural Engineering, ASCE, Vol. 132, No. 12, pp. 2027-2035.         [ Links ]

Hong, H P, K Goda y A G Davenport (2006), "Seismic hazard analysis: a comparative study", Canadian Journal of Civil Engineering, Vol. 33 No. 9, pp. 1156-1171.         [ Links ]

Hong, H P, y Hong P. (2007), "Assessment of ductility demand and reliability of bilinear single-degree-of-freedom-systems under earthquake loading", Canadian Journal of Civil Engineering, 34(12), 1606-1615.         [ Links ]

Hong, H P, A D García-Soto, y R Gomez (2010), "Impact of different earthquake types on the statistics of ductility demand", Journal of Structural Engineering, ASCE, Vol. 136 No. 7, pp. 770-780.         [ Links ]

Hong, H P, A Pozos-Estrada y R Gomez (2009), "Orientation effect on ground motion measure for Mexican subduction earthquakes", Earthquake Engineering and Engineering Vibration, Vol. 8, No. 1, pp. 1-16.         [ Links ]

Jaimes M A y E Reinoso (2006), "Comparación del comportamiento de edificios en el valle de México ante sismos de subducción y de falla normal", Revista de Ingeniería Sísmica, No. 75, pp. 1-22.         [ Links ]

Joyner, W B y D M Boore (1993), "Methods for regression analysis of strong-motion data", Bulletin of the Seismic Society of America, Vol. 83, pp. 469-487.         [ Links ]

Madsen, H O, S Krenk, y N C Lind (1986), Methods of Structural Safety, Prentice-Hall, Inc., Englewood Cliffs, New Jersey.         [ Links ]

McGuire, R K (1976), "FORTRAN Computer program for seismic risk analysis", Open-File 76-67, U.S. Department of the Interior, U.S. Geological Survey, Denver, CO.         [ Links ]

McGuire, R K (2004), "Seismic hazard and risk analysis", Oakland, Calif.: Earthquake Engineering Research Institute.         [ Links ]

MOC (1993), Manual de diseño de obras civiles, Diseño por sismo, Instituto de Investigaciones Eléctricas, Comisión Federal de Electricidad.         [ Links ]

MOC (2008), Manual de diseño de obras civiles, Diseño por sismo. Recomendaciones y Comentarios, Instituto de Investigaciones Eléctricas, Comisión Federal de Electricidad.         [ Links ]

Nakashima M (1997), Uncertainties associated with ductility performance of steel building structures, Seismic design methodology for the next generation of codes (eds. Fajfar P, Krawinkler H.), Balkema, Rotterdam.         [ Links ]

NBCC (2005), National Building Code of Canada, Canadian Commission for Building and Fire Codes, National Research Council Canada, Ottawa.         [ Links ]

NTCS (2004), Normas técnicas complementarias para diseño por sismo, Reglamento de Construcciones para el Distrito Federal, Gaceta Oficial del Departamento del Distrito Federal.         [ Links ]

Ordaz, M y C Reyes (1999), "Earthquake hazard in Mexico City: Observations versus computations", Bulletin of the Seismic Society of America, Vol. 89, No. 5, pp. 1379-1383.         [ Links ]

Ordaz, M, L E Pérez-Rocha y J Alemán Velásquez (2007), "Definición de sismos máximos para el diseño sísmico de estructuras", XVI Congreso Nacional de Ingeniería Sísmica, Ixtapa-Zihuatanejo, Guerrero, México.         [ Links ]

Pardo, M y G Suárez (1995), "Shape of the subducted Rivera and Cocos plates in southern Mexico: Seismic and tectonic implications", Journal of Geophysical Research, Vol. 100, No. B7, pp. 12,357-12373.         [ Links ]

Pozos-Estrada, A, R Gómez y H P Hong (2008), "Seismic design: benefit/cost for overall service time versus per unit service time", Proceedings of the 14th World Conference on Earthquake Engineering, Beijing.         [ Links ]

Rackwitz, R (2000), "Optimization-the basis of code-making and reliability verification", Structural Safety, Vol. 22, pp. 27-60.         [ Links ]

Reyes, C, E Miranda, M Ordaz y R Meli (2002), "Estimación de espectros de aceleraciones correspondientes a diferentes periodos de retorno para las distintas zonas sísmicas de la ciudad de México", Revista de Ingeniería Sísmica, No. 66, pp. 95-121.         [ Links ]

Rosenblueth, E y L Esteva (1972), "Reliability Basis for Some Mexican Codes Publication", ACI Special Publication, Vol. 31, pp. 1-42.         [ Links ]

Rosenblueth, E y R Gómez (1991), "Comentarios y ejemplos a las normas técnicas complementarias para diseño por sismo", DDF, Series de Instituto de Ingeniería, ES-7, UNAM, México.         [ Links ]

Rosenblueth, E y J M Jara (1991), "Constant versus time dependent seismic design coefficients", Proceedings of the 3rd IFIP WG 7.5 Conference on Reliability and Optimization of Structural Systems '90, Springer-Verlag, Berlin, Germany, Vol. 1, pp. 315-327.         [ Links ]

SCT (2001), Nueva normatividad para diseño de puentes, página electrónica del IMT y la SCT, http://normas.imt.mx/default.htm, N-PRY-CAR-6-01-001/01 a N-PRY-CAR-6-01-006/01, Secretaria de Comunicaciones y transportes.         [ Links ]

Singh, S K, R Quaas, M Ordaz, F Mooser, D Almora, M Torres y R Vásquez (1995), "Is there truly a hard-rock site in the valley of Mexico", Geophysical Research Letters, Vol. 22, No. 4, pp. 481-484.         [ Links ]

Tena-Colunga, A, U Mena, L E Pérez-Rocha, M Ordaz, J Avilés y J I Vilar (2008), "Updated guidelines for the seismic design of buildings in the MOC-2008 code of Mexico", Proceedings of the 14th World Conference on Earthquake Engineering, Beijing, China.         [ Links ]

Zúñiga, R, G Suárez, M Ordaz y V García-Acosta (1997), "Peligro sísmico en Latinoamérica y el Caribe", Reporte Final IPGH.         [ Links ]