SciELO - Scientific Electronic Library Online

 
 número91Influencia de contraventeos concéntricos en cruz o chevrón en el comportamiento no lineal de marcos dúctiles de concreto reforzadoIdentificación de propiedades dinámicas de un modelo estructural sometido a vibración ambiental y vibración forzada empleando mesa vibradora índice de autoresíndice de materiabúsqueda de artículos
Home Pagelista alfabética de revistas  

Servicios Personalizados

Revista

Articulo

Indicadores

Links relacionados

  • No hay artículos similaresSimilares en SciELO

Compartir


Ingeniería sísmica

versión impresa ISSN 0185-092X

Ing. sísm  no.91 Ciudad de México jul./dic. 2014

 

Articles

Desagregación del peligro sísmico para algunos sitios seleccionados de México

Seismic hazard deaggregation for some selected sites of Mexico

Adrián Pozos-Estrada1 

Roberto Gómez Martínez1 

Hanping Hong2 

1Investigador, Instituto de Ingeniería, Universidad Nacional Autónoma de México, Coyoacán, C.P. 04510, México, D.F. México, APozosE@iingen.unam.mx RGomezM@iingen.unam.mx

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


RESUMEN

Describimos la desagregación del peligro sísmico con base en un modelo aplicable a una región de la costa del Pacífico Mexicano, tomando en cuenta las zonas de fuentes sísmicas, el modelo de recurrencia, y relaciones de atenuación. También desarrollamos relaciones de atenuación aplicables a la ciudad de México (i.e., la estación de registro de CU) y las integramos en el modelo de peligro sísmico adoptado. Proveemos mapas de peligro sísmico para la región de la costa del Pacífico, y espectros de peligro uniforme para algunos sitios seleccionados, incluyendo la estación CU; notamos que las ordenadas del espectro de diseño para la Ciudad de México son mayores a los valores del espectro de peligro uniforme para un periodo de retorno de 125 años. Los escenarios sísmicos identificados muestran que, en ocasiones, la contribución al peligro sísmico para un periodo de retorno de 125 o 475 años es dominada por un solo evento y en otros casos el peligro sísmico es debido a diferentes escenarios sísmicos y tipo de temblores; en particular, para la estación CU, se observa que para un periodo de vibrar menor que aproximadamente 0.7 (s) el peligro sísmico es dominado por sismos intraslab y que los sismos interplaca dominan el peligro sísmico para periodos de vibrar mayores que éste. Para la costa del Pacífico, el peligro sísmico es dominado por sismos interplaca, y por sismos intraslab e interplaca para zonas alejadas de la costa.

Palabras clave: desagregación sísmica; espectros de peligro uniforme; sismos interplaca; sismos intraslab

ABSTRACT

We describe the deaggregation of seismic hazard based on a model applicable to a region of the Pacific coast of Mexico, taking into account the seismic source zones, the recurrence model, and the attenuation relations. We also developed attenuation relations for Mexico City (i.e., the recording station of CU) and integrated them in the adopted seismic hazard model. We provide contour maps of the seismic hazard for the Pacific coast of Mexico, and uniform hazard spectra for some selected sites, including the CU station; we noted that the ordinates of the design spectrum for Mexico City are greater than those for the uniform hazard spectrum for a return period of 125 years. The seismic scenarios identified show that, on occasion, the contribution to the seismic hazard for a return period of 125 or 475 years is dominated by a single event and in other cases the seismic hazard is due to different seismic scenarios and earthquake types; in particular for the CU station, it is observed that for a vibration period less than about 0.7 (s) the seismic hazard is dominated for inslab events and that interplate events dominate the seismic hazard for vibration periods greater than that. For the Pacific coast, the seismic hazard is dominated by interplate earthquakes, and by intraslab and interplate events for zones away from the coast.

Keywords: seismic deaggregation; uniform hazard spectra; interplate earthquakes; inslab earthquakes

Introducción

La desagregación del peligro sísmico es una forma eficiente de identificar combinaciones de valores de parámetros tales como la magnitud y la distancia de la fuente al sitio de interés, que contribuyen a un nivel de peligro sísmico específico (McGuire1995; Chapman 1995; Bazzurro y Cornell 1999; Harmsen et al. 1999). El peligro sísmico puede expresarse en términos del espectro de respuesta. Los eventos identificados que están asociados con altas frecuencias pueden ser empleados como escenarios para evaluar estructuras ingenieriles e infraestructura nueva y existente. Las características de los escenarios también pueden ser considerados para seleccionar registros del movimiento del terreno que se necesitan para analizar estructuras y verificar el diseño sísmico. Una aplicación de esto fue hecha por Reinoso y Jaimes (2009) para obtener registros sísmicos del movimiento del terreno que son compatibles con el espectro de diseño para la ciudad de México.

En lugar de considerar la desagregación basada en el espectro de respuesta, Baker et al. (2005) consideraron la desagregación del peligro sísmico en términos de la distorsión máxima de entrepiso. Hong y Goda (2006) argumentaron que debido a que las estructuras son diseñadas y construidas para proveer servicio en lugar de sólo sostener movimientos sísmicos, la desagregación del riesgo sísmico podría ser igual o más importante que aquella del peligro sísmico. Ellos aplicaron el concepto de desagregación sísmica a algunos sitios de Canadá, donde el riesgo sísmico es expresado como el costo esperado del ciclo de vida, el cual depende de la probabilidad de la fluencia incipiente y el colapso incipiente. Sin embargo, una función del costo del ciclo de vida no está disponible para la construcción mexicana, excepto la dada por Rosenblueth y Jara (1991), la cual fue empleada por otros para seleccionar el diseño óptimo (García-Pérez et al. 2005; Pozos-Estrada et al. 2008). Aunque las reglas de demanda de ductilidad necesarias considerando sismos interplaca e intraslab para la desagregación del riesgo sísmico están disponibles, desafortunadamente, nuestra experiencia con la función de costo dada por Rosenblueth y Jara (1991) indica que esta expresión tiene un comportamiento no realista conforme el estado de daño estructural se acerca al colapso. Por lo tanto, esta función no es considerada y la desagregación del riesgo sísmico no se lleva a cabo en el presente trabajo.

El principal objetivo de este estudio es aplicar los métodos de desagregación del peligro sísmico para identificar eventos o escenarios frecuentes que producen un nivel de peligro predefinido para varios sitios seleccionados de México. Los análisis se llevan a cabo considerando diferentes zonas de fuente sísmicas y tipos de sismos. Para los análisis, se considera un modelo de peligro sísmico que es ensamblado por García-Soto et al. (2012) para estimar espectros de peligro uniforme (EPU) y desarrollar mapas de contorno para una región de México. Este modelo de peligro sísmico toma en cuenta: a) las zonas de fuentes sísmicas y el modelo de recurrencia descrito en Ordaz y Reyes (1999) y Jaimes y Reinoso (2006), y b) coeficientes de relaciones de atenuación (o ecuaciones de predicción de movimiento fuerte), las cuales son desarrolladas con base en la media geométrica de las respuestas en dos direcciones ortogonales horizontales dadas por Pozos-Estrada et al. (2008) y Hong et al. (2009b). El uso de la relación de atenuación con base en la media geométrica es justificado ya que el término de error en la relación de atenuación es aplicable para una orientación aleatoria. Como las relaciones de atenuación basadas en la media geométrica y aplicables para la Ciudad de México (i.e., las estaciones de registro de Ciudad Universitaria, las cuales serán referidas como la estación CU por sencillez) no están disponibles para sismos interplaca e intraslab, incluyendo sismos más recientes, éstas son desarrollados e integradas al modelo en este estudio. Los resultados de la desagregación basados en el modelo de peligro sísmico adoptado y sus implicaciones para el diseño sísmico y la planeación para emergencia son discutidos.

Modelo de peligro sísmico y desagregación

Modelo de peligro sísmico

La estimación del peligro sísmico requiere información de sismicidad definida por la caracterización de modelos de ocurrencia de los sismos, las zonas sismogénicas (fuentes sísmicas) y sus relaciones de magnitud-recurrencia correspondientes, y leyes o relaciones de atenuación (o ecuaciones de predicción de movimiento fuerte). Dicha información se puede inferir a través de catálogos y registros de movimiento del terreno históricos.

La ocurrencia de un sismo es comúnmente considerada estadísticamente homogénea en magnitud y en ubicación dentro de una zona de la fuente sísmica. La ocurrencia de un sismo en sitios de México es comúnmente modelada como un proceso Poisoniano (Esteva 1968), aunque procesos no Poisonianos han sido también considerados en la literatura (Singh et al. 1983; Hong y Rosenblueth 1988; Rosenblueth y Jara 1991; Pozos-Estrada et al. 2008). La definición de las zonas sismogénicas para regiones sísmicas activas se puede llevar a cabo basada en la actividad sísmica histórica e incorporando características geológicas y tectónicas. Para cada zona sismogénica, la magnitud de los sismos es tratada como una variable aleatoria. Es común considerar que las magnitudes siguen la relación magnitud-recurrencia de Gutenberg-Richter (i.e., distribuciones exponencial o exponencial truncada). Para algunas zonas sismogénicas en México, otros modelos probabilísticos han sido también considerados (Rosenblueth y Ordaz 1987; Ordaz y Reyes 1999) debido a la ocurrencia de sismos característicos (Singh et al. 1983). El empleo de zonas sismogénicas homogéneas con fronteras definidas es comúnmente empleado y es considerado en el presente estudio también. Esto puede ser justificado al considerar que la mayoría de las zonas sismogénicas consideradas son zonas sismogénicas altamente activas y dominadas por fallas bien caracterizadas en México (García 2006).

El modelo de peligro sísmico adoptado es el mismo modelo considerado por García-Soto et al. (2012), excepto que una relación de atenuación adicional para eventos intraslab es considerada para la Ciudad de México, así como la relación de atenuación para eventos interplaca considerando eventos sísmicos más recientes. El modelo considera las fuentes sísmicas mostradas en la Figura 1, en donde las fuentes símicas son separadas en tres grupos, dependiendo del tipo de sismo y magnitud considerados. Los dos primeros grupos, ubicados cerca de la costa del Pacífico Mexicano, están asociados con sismos interplaca. El primero de ellos considera una fuente sísmica definida por 4 polígonos, los cuales abarcan una zona de subducción acotada aproximadamente por los estados de Jalisco y Chiapas. Las magnitudes de los sismos dentro de esta zona sísmica se consideran inferiores a 7 (Ordaz y Reyes 1999). El segundo grupo considera una fuente sísmica conformada por 14 polígonos con extensión similar a la primera zona definida. De acuerdo con Ordaz y Reyes (1999), las magnitudes de los sismos dentro de esta zona sísmica se consideran mayores a 7. El tercer grupo, que considera los sismos intraslab, tiene una fuente sísmica definida por dos grandes polígonos ubicados dentro de la placa oceánica subducida (Jaimes y Reinoso 2006), (ver también Nishenko y Singh 1987; Reinoso y Jaimes 2009). Parte de esta zona sísmica se ubica por debajo del Distrito Federal.

Figura 1 Zonas sismogénicas: a) sismos interplaca con M W ≤ 7, b) sismos interplaca con M W > 7, c) sismos intraslab, d) Cinturón Volcánico Mexicano (véase Arroyo et al. (2010))  

Además de las fuentes sísmicas anteriormente descritas, también se encuentran las fuentes sísmicas local y continental. De acuerdo con Rosenblueth et al. (1989), los eventos locales afectan sólo una pequeña zona, y los eventos que impactarían cualquier ubicación del Distrito Federal tienen periodos de retorno significativamente mayores a aquellos correspondientes a la región completa. Asimismo, los eventos continentales, están asociados con periodos de retorno superiores a los 1000 años. Como el peligro sísmico asociado con este tipo de eventos se considera bajo, éstos no son considerados en este estudio.

La distribución de probabilidad de la magnitud y la tasa de ocurrencia anual para cada una de las zonas mostradas en la Figura 1 son presentadas en la Tabla 1. Estos valores fueron tomados de Rosenblueth et al. (1989), Ordaz y Reyes (1999), y Jaimes y Reinoso (2006) (las coordenadas para las zonas que son inferidas de las dos últimas referencias están tabuladas en García-Soto et al. 2012).

Tabla 1 Relación magnitud-recurrencia para las zonas sismogénicas mostradas en la Figura 1 (después de Ordaz y Reyes 1999, Jaimes y Reinoso 2006

Notas: El valor en la Tabla 2 en Ordaz y Reyes (1999) es interpretado como λ7; λ(MW ) representa la tasa de ocurrencia de un sismo con magnitud mayor o igual a MW . Los parámetros λ0, λ7 y β se muestran en la Figura 1.

La forma funcional de las relaciones de atenuación adoptadas en el presente estudio para sitios fuera del Valle de México y del Cinturón Volcánico Mexicano (CVM) (véase Figura 1d) (Arroyo et al. 2010) se muestra en la Tabla 2, donde Y denota la aceleración espectral, SA (T n,ξ) en Gal; T n (s) es el periodo fundamental de vibrar; ξ es el porcentaje de amortiguamiento crítico, el cual es igual a 5% en este estudio. Las formas funcionales para logaritmo de SA (T n,ξ) correspondientes a sismos interplaca e intraslab son las sugeridas por García (2006) , y García et al. (2005), respectivamente.

Tabla 2 Coeficientes de relaciones de atenuación 

Notas para la Tabla 2:

1)ε en las ecuaciones es para una orientación aleatoria y es considerado como una variable normal con media cero y desviación estándar ϭ.

2) Para la relación de atenuación de sismos interplaca para sitios fuera del Valle de México y del Cinturón Volcánico Mexicano, R (km) es la distancia más cercana a la superficie de falla para eventos con Mw > 6.0, o la distancia hipocentral para el resto; H (km) es la profundidad focal; c 4 = 1.82 - 0.16Mw ; y c 1, c 2, c 3, y c 7 son coeficientes de regresión. Se consideró un total de 413 acelerogramas con dos componentes horizontales obtenidos de 40 eventos sísmicos.

3) Para los sismos intraslab para sitios fuera del Valle de México y del Cinturón Volcánico Mexicano, c i, i = 1,2,3,4,5, representa coeficientes de regresión; c 4 =1; R=(R cld 20)1/2; R cld (km) es la distancia más cercana a la superficie de falla para eventos con Mw > 6.5, o la distancia hipocentral para el resto; H (km) es la profundidad focal; Δ0=0.0075x100.507M w es un término de saturación de la fuente definido por Atkinson y Boore (2003). Se consideraron 275 acelerogramas con dos componentes horizontales obtenidos de 16 eventos sísmicos.

4) Para la relación de atenuación para sismos interplaca para Ciudad Universitaria, R (km) es la distancia más cercana a la superficie de falla; c i, i = 1,...,5 representan coeficientes de regresión. Se consideró un total de 18 acelerogramas con dos componentes horizontales obtenidos de 18 eventos sísmicos. La ecuación T2.c es válida para M W ≤ 8.1; en caso de eventos con M W > 8.1, sus magnitudes son consideradas igual a 8.1 (Reyes, comunicación personal 2007). La misma consideración es adoptada para la ecuación T2.a.

5)Para la relación de atenuación para sismos intraslab para Ciudad Universitaria ver 3) para las definiciones de los parámetros. Se consideraron 18 y12 acelerogramas con dos componentes horizontales obtenidos de 18 y 12 eventos sísmicos, respectivamente. La magnitud máxima para los sismos intraslab se considera menor o igual a 7.4, magnitud máxima observada en la base de datos empleada.

6)Se consideran profundidades deterministas de 22.22 km para sismos interplaca y 54.61 km para sismos intraslab. Estos valores representan la profundidad promedio de los sismos asociados con los registros utilizados en el desarrollo de las leyes de atenuación para la estación CU y son similares a los valores obtenidos con los registros utilizados para el desarrollo de las leyes de atenuación para sitios fuera del CVM y del Valle de México. Lo anterior se realiza ya que no se dispone de modelos probabilísticos asociados con esta variable. Nótese que en García-Soto et al. (2012), el valor de la profundidad promedio para sismos interplaca, 10.45 (km), debería ser leído como 22.33 (km).

Los coeficientes para las Ecuaciones (T2.a) y (T2.b), tomados de Hong et al. (2009b), se ilustran en la Tabla 2 sólo para 3 periodos de vibrar, aunque los coeficientes para otros periodos de vibrar están disponibles en la referencia citada. Estos coeficientes fueron estimados usando el algoritmo propuesto por Joyner y Boore (1993) para la media geométrica. El uso de la media geométrica se justifica, ya que en este caso el término de error ε representa la incertidumbre asociada en la predicción del movimiento del terreno para una orientación aleatoria (Boore et al. 1997). Nótese que en todos los casos, la desviación estándar de ε, σ, para las ecuaciones de regresión desarrolladas con base en la media geométrica (Hong et al. 2009b) es mayor que aquella basada en la media cuadrática dada por García et al. (2005) y García (2006).

Aunque relaciones de atenuación para la estación CU que consideran eventos interplaca están disponibles (Reyes et al. 2002, Pozos-Estrada et al. 2008), éstas no son consideradas debido a que no incluyen eventos más recientes. Sin embargo, la forma funcional de la relación de atenuación dada por Reyes et al. (2002), la cual es mostrada en la Tabla 2, es adoptada. Para sismos intraslab, Montalvo-Arrieta (2002) propuso la forma funcional que es idéntica a la ecuación (T2.b), mostrada en la Tabla 2, excepto que el término asociado a la profundidad focal (i.e., c 5) es considerado igual a cero. Es importante notar que al no tomarse en cuenta la profundidad focal, existen implicaciones al estimar aceleraciones espectrales para distancias pequeñas entre la falla y el sitio de interés. Más específicamente, debido a que parte de la zona sismogénica adoptada para sismos intraslab se ubica por debajo del Distrito Federal, al no considerarse la profundidad focal se sobrestima el peligro sísmico para tales eventos. Para evitar esta sobrestimación, en lugar de adoptar la ecuación recomendada por Montalvo-Arrieta (2002), la ecuación sugerida por García et al. (2005) (i.e., ecuación T2.b, ver Tabla 2) fue considerada.

Para desarrollar los coeficientes de las relaciones de atenuación para la estación CU, se emplearon los eventos sísmicos de las Tablas 3 y 4 correspondientes a sismos interplaca e intraslab, respectivamente. La lista de eventos mostrados en las Tablas fue desarrollada con base en estudios de Reyes et al. (2002), Montalvo-Arrieta et al. (2002), y Singh et al. (2013). Debido a que hay varios registros disponibles de estaciones de registro cercanas entre ellas dentro de Ciudad Universitaria para los sismos ocurridos el 19 y 21 de septiembre de 1985, sólo un registro de la estación denotada como CU01 para cada uno de los eventos mencionados es usado. Para los registros considerados, un análisis de regresión empleando el algoritmo de Joyner y Boore (1993) (Hong et al. 2009a) se lleva a cabo considerando los eventos mostrados en la Tabla 3, y sus correspondientes acelerogramas obtenidos para la estación CU.

Tabla 3 Eventos sísmicos empleados en el desarrollo de la relación de atenuación para Ciudad Universitaria considerando sismos interplaca (h = profundidad) 

Tabla 4 Eventos sísmicos empleados en el desarrollo de la relación de atenuación para Ciudad Universitaria considerando sismos intraslab (h = profundidad) 

Los coeficientes de la regresión considerando los sismos interplaca se presentan en la Figura 2 así como en la Tabla 2 para algunos valores de T n. Es interesante notar que los valores de σ de la relación de atenuación para los acelerogramas obtenidos en la estación CU para los sismos interplaca son de alrededor del 60% de aquellos asociados con la ecuación (T2.a). Esto es consistente con los reportados por Pozo-Estrada et al. (2008), y concuerda con la observación hecha por Atkinson (2006) que indica que para el mismo tipo de sismo σ en una sola estación es de alrededor del 60% que aquel obtenido de los acelerogramas en una región.

Figura 2 Coeficientes de regresión para sismos interplaca para la estación CU 

Para el caso de eventos intraslab, los espectros de los acelerogramas considerados se presentan en la Figura 3. En la Figura se incluyen los espectros para la media geométrica (Figura 3a y 3b ) y para cada componente horizontal del terreno (Figuras 3c, 3d , 3e y 3f ). Se observa de la Figura 3 que los espectros presentan picos predominantes por debajo de 1 (s). En general, se observan picos asociados a periodos de vibrar cercanos a 0.2 y 0.5 (s). También se observa que los espectros de algunos sismos presentan un contenido importante de energía para periodos mayores a 1 (s). Lo anterior es consistente con lo observado por Montalvo-Arrieta (2002).

Figura 3 Espectros de respuesta para sismos intraslab para la estación CU (ξ = 5%): a) y b) para la media geométrica; c) y d) primer componente horizontal; e) y f) segunda componente horizontal 

Los resultados del análisis de regresión indican que los coeficientes de la regresión varían significativamente con el periodo natural de vibrar T n como se muestra en la Figura 4 (ver también la Tabla 2). También se muestra en la Figura 4 la desviación estándar ϭ. La Figura indica que ϭ aumenta conforme T n aumenta, sus valores aumentan hasta 0.5 lo cual es significativamente mayor que aquellos para la ecuación (T2.b) que fueron desarrollados con base en eventos intraslab, pero para registros de muchas estaciones. Esto es inesperado ya que ϭ para una sola estación considerando el mismo tipo de sismo debe ser menor que aquel asociado con muchas estaciones (Atkinson 2006). Una inspección de los residuos ε del análisis de regresión indica que éstos se desvían significativamente de una variable normal, lo que indica que el modelo de regresión es inadecuado para los registros considerados.

Figura 4 Coeficientes de regresión para sismos intraslab para la estación CU, considerando todos los acelerogramas 

Para mejorar el ajuste usando el mismo modelo adoptado, los registros cuyo valor máximo de SA (T n,ξ), incluyendo la PGA, es menor que 10 (Gal) fueron eliminados. Esto resultó en sólo 12 registros para la estación CU. En este caso, los coeficientes y ϭ obtenidos se muestran en la Figura 5 así como en la Tabla 2 para algunos valores seleccionados de T n. La comparación de los resultados mostrados en las Figuras 4 y 5 indica que al eliminar los registros cuyo valor máximo de SA (T n,ξ) es menor que 10 (Gal), los valores de ϭ son menores a los obtenidos considerado todos los registros. También se observa que el valor máximo de ϭ para el caso seleccionado es aproximadamente igual a 0.30. Este valor es comparable con aquellos obtenidos fuera del Valle de México y del CVM basados en una amplia red de estaciones para sismos intraslab como se muestra en la Tabla 2.

Figura 5 Coeficientes de regresión para sismos intraslab para la estación CU, considerando acelerogramas cuyo valor máximo de SA (T n,ξ) es mayor que 10 (Gal) 

Aunque se puede desarrollar una relación de atenuación con base en los registros con excitaciones mayores que un límite inferior mayor que 10 (Gal), esto no fue considerado ya que se reduciría el número de acelerogramas que pueden ser usados para desarrollar las relaciones de atenuación. La consideración de un límite inferior de 10 (Gal) en el desarrollo de la relación de atenuación puede ser una suposición conservadora. Debido a que la selección de un límite inferior de 10 (Gal) es de alguna forma arbitraria, con fines de comparación, la evaluación del peligro sísmico para la estación CU es llevada a cabo en las siguientes secciones para los casos CU-1 y CU-2. El caso CU-1 se refiere a la evaluación del peligro sísmico empleando la relación de atenuación desarrollada para registros de eventos intraslab con valores máximos de SA (T n,ξ) mayores que 10 (Gal), mientras que el caso CU-2 se refiere a la evaluación del peligro sísmico considerando la relación de atenuación desarrollada con base en todos los eventos intraslab mostrados en la Tabla 4.

Mapas de contorno, espectros de peligro uniforme y resultados de la desagregación

Procedimiento para estimar y desagregar el peligro sísmico

Una vez que un modelo de peligro sísmico es asignado, el peligro sísmico en el sitio de interés, definido por la distribución de probabilidad de las mediciones del movimiento del terreno, puede ser evaluado, y el análisis de desagregación se puede llevar a cabo. Uno de los procedimientos más populares para la estimación probabilista del peligro sísmico (EPPS) fue desarrollado por Esteva (1968) y Cornell (1968). Este procedimiento se ilustra en la Figura 6. El algoritmo para evaluar el peligro sísmico en un sitio puede ejecutarse con base en integración numérica (Cornell 1968; McGuire 2004) o técnicas de Monte Carlo (Musson 2000; Hong et al. 2006). El uso de técnicas de Monte Carlo ofrece la ventaja de que éste puede ser empleado directamente para tratar con modelos de ocurrencia Poisonianos o no Poisonianos homogéneos o no homogéneos, y puede ser fácilmente adaptado para estimar espectros de peligro uniforme (EPU) y el análisis de desagregación de peligro o riesgo sísmico (Hong et al. 2006; Hong y Goda 2006).

Figura 6 Ilustración del método del cálculo del peligro sísmico (Modificado de García-Soto et al. 2012, y Hong 2012

La estimación del EPU consiste en encontrar las respuestas espectrales para diferentes periodos naturales de vibrar con una probabilidad de excedencia específica (i.e., fractiles de SA (Tn ,ξ) denotados por SAf (Tn ,ξ)). La desagregación del peligro sísmico consiste en identificar la combinación de valores de las variables aleatorias que caracterizan a los eventos sísmicos que contribuyen a SAf(Tn ,ξ). A menudo, las variables consideradas para la desagregación son la magnitud, la distancia de la fuente al sitio de interés, y ε en la relación de atenuación. Además de estas tres variables, también incluimos el tipo de sismo y la profundidad del sismo en el análisis de desagregación. Los pasos para el análisis son los siguientes:

  1. Llevar a cabo la simulación para determinar muestras de SA (Tn ,ξ) para un periodo de vibrar y amortiguamiento dados;

  2. Encontrar SAf (Tn ,ξ) mediante un análisis probabilístico de las muestras;

  3. Recolectar todos los eventos cuya respuesta sea igual o cercana al valor de SAf (Tn ,ξ), y llevar a cabo un análisis estadístico con los parámetros que caracterizan a los eventos recolectados.

Los primeros dos pasos para el análisis son iguales a los mostrados en la Figura 6, y son empleados en los análisis numéricos.

Mapas de contorno y espectros de peligro uniforme

Se elaboraron mapas de contorno de SAf (Tn ,ξ) para evaluar el peligro sísmico en algunas regiones de México, incluyendo al Distrito Federal en terreno firme (Ciudad Universitaria). Para el cálculo del peligro símico se empleó el método de simulación descrito en secciones anteriores. Los valores simulados representan 187,500 años de actividad sísmica. Los resultados de SAf (T n,ξ) para un periodo de retorno igual a 125 años, considerando la opción CU-1 para la estación CU, se presentan en la Figura 7. Para evitar confusión en la lectura de los niveles de aceleración presentados en las gráficas de la Figura 7, éstas deben asociarse con la graduación de aceleración correspondiente, cuyo valor máximo difiere para cada caso. La Figura 7 muestra el CVM sombreado en blanco desvanecido (ver Figura 2d). Es importante recalcar que los resultados de la simulación del peligro sísmico, determinados con las relaciones en la Tabla 2, no deben ser extrapolados al CVM ni al Valle de México (excepto para la estación CU), debido a las características particulares de propagación a lo largo de estas regiones (García 2006). Se observa de la Figura 7 que, como se esperaba, el peligro sísmico es mayor para sitios cercanos a la costa. Se observa también que el impacto del peligro sísmico es mayor para estructuras con Tn = 0.2 (s) y que éste disminuye para estructuras con Tn igual a 0.5 y 2 (s), en particular para Tn = 2 (s), se observa una reducción significativa del peligro sísmico.

Figura 7 Mapas de contorno considerando sismos interplaca e intraslab para un periodo de retorno igual a 125 años 

Aunque mapas para periodos de retorno iguales a 475, 1000 y 2475 años también fueron evaluados para diferentes valores de T n, estos no son presentados ya que son idénticos (excepto para la estación en Ciudad Universitaria) a los presentados por García-Soto et al. (2012).

Para apreciar mejor el peligro sísmico para diferentes periodos naturales de vibrar, de forma similar a García-Soto et al. (2012), calculamos el EPU para la estación CU considerando un periodo de retorno igual a 125 años y dos opciones para los sismos de intraslab (CU-1 y CU-2). También calculamos el EPU para otros cuatro lugares que se localizan cerca de Morelia, Acapulco, Oaxaca y Salina Cruz (véase la Figura 1). Estos lugares se refieren como el Sitio I, II, III y IV en lo siguiente. Los resultados obtenidos son mostrados en la Figura 8. Se observa de la Figura que, al igual que en los mapas de contorno, el peligro sísmico es más significativo para los sitios más cercanos a la costa y para Tn cercano a 0.2 (s). También se observa que la forma de los EPU para los Sitios I a IV es similar. El EPU para la estación CU mostrado en la Figura 8 es debido a sismos intraslab e interplaca. Se observa que este EPU tiene una forma muy diferente que los EPU para los Sitios I a IV, principalmente para Tn menores a 1.7 (s), en donde se aprecian periodos de vibrar característicos con un peligro sísmico importante para el valle de México. El EPU para el Sitio CU, usando CU-1, es comparable en magnitud con los EPU para los Sitios I a IV para Tn grandes. Asimismo, la Figura muestra que hay diferencias marcadas entre el EPU obtenido usando la opción CU-1 y CU-2. El EPU empleando CU-2 es mucho más grande que aquel obtenido con CU-1; esto indica que el peligro sísmico en la estación CU es altamente afectado por la relación de atenuación adoptada, incluyendo ε, para sismos intraslab.

Figura 8 Espectros de peligro uniforme para los Sitios I, II, III, IV y para la estación CU (empleando las opciones CU-1 y CU-2) para un periodo de retorno de 125 años 

Para apreciar mejor la influencia de la relación de atenuación y los sismos intraslab e interplaca en el peligro sísmico, calculamos nuevamente el EPU debido sólo a sismos interplaca, sólo a sismos intraslab (con CU-1 o CU-2) y los sismos interplaca e intraslab para periodos de retorno iguales a 125 y 475 años. Los valores obtenidos son mostrados en la Figura 9 y comparados con el espectro de diseño recomendado por las Normas Técnicas Complementarias para el Diseño por Sismo (NTCDS) (Diario Oficial de la Federación, 2004) para la Ciudad de México para terreno firme. Las gráficas mostradas en la Figuras 9a y 9b indican que los sismos interplaca e intraslab contribuyen al peligro sísmico en la estación CU. Asimismo, el uso de CU-2 para sismos intraslab sobrestima el peligro sísmico comparado con aquellos usando CU-1. Esto claramente indica la importancia de usar una relación de atenuación (y ε) apropiada para la evaluación del peligro sísmico; el uso de una relación de atenuación tal como la CU-2 es cuestionable y puede llevar a obtener un peligro sísmico erróneo. Por lo tanto, el uso de CU-2 no es considerado más en lo siguiente.

Figura 9 Comparación de espectros de peligro uniforme con el espectro de diseño para terreno firme de las Normas Técnicas Complementarias 

Nótese que al considerar CU-1 para eventos intraslab, se observa en la Figura 9a que las ordenadas espectrales del EPU para sismos intraslab son similares a las del espectro de diseño hasta un periodo de vibrar cercano a 0.2 (s), después de este periodo las ordenadas espectrales disminuyen hasta Tn cercano a 1 (s), donde se observa nuevamente un leve incremento en las ordenadas espectrales, después de este periodo de vibrar las ordenadas espectrales del EPU disminuyen gradualmente para periodos grandes. Se observa en la misma Figura que el EPU para sismos interplaca presenta ordenadas espectrales máximas para periodos de vibrar cercanos a 1 y a 1.7 (s). Los EPU para sismos intraslab e interplaca se cruzan en Tn =0.7 (s), aproximadamente, indicando que el peligro sísmico para una estructura con un periodo fundamental menor que 0.7 (s) es dominado por sismos intraslab mientras el peligro sísmico para estructuras con un periodo fundamental mayor que 0.7 (s) es dominado por sismos interplaca. Las ordenadas espectrales de los EPU de sismos intraslab o de sismos interplaca son menores a las del espectro de diseño para el periodo de retorno de 125 años.

Similar a la Figura 9a, la Figura 9b presenta la comparación de los EPU pero para un periodo de retorno de 475 años. En general, las tendencias observadas en Figura 9a son igualmente aplicables a Figura 9b, excepto que, y como era de esperarse, las ordenadas del EPU en la Figura 9b son mayores que aquellas de la Figura 9a.

La Figura 9c muestra los EPU considerando la contribución de los sismos intraslab e interplaca para periodos de retorno iguales a 125 y 475 años, empleando sólo CU-1 para sismos intraslab. La Figura indica que el espectro de diseño para terreno firme de NTCDS (Diario Oficial de la Federación, 2004) es mayor que el EPU para a un periodo de retorno de 125 años y, el espectro de diseño es más grande que el EPU para un periodo de retorno de 475 años para Tn entre 0.25 y 1.75 (s). Esta observación difiere de la conclusión que el espectro de diseño corresponde a un EPU para un periodo de retorno de 125 años (Ordaz et al. 2000). La diferencia podía ser explicada al notar que registros adicionales de eventos sísmicos más recientes son considerados en el presente estudio.

Resultados de la desagregación

Como se mencionó previamente, la desagregación sísmica consiste en identificar combinaciones de variables aleatorias que definen eventos sísmicos que contribuyen a un nivel de peligro sísmico predefinido. Estas variables aleatorias incluyen Mw, R y ε, empleados en las relaciones de atenuación.

La desagregación del peligro sísmico se realiza para los sitios identificados como I, II, III y IV, junto con la estación CU, para la aceleración máxima del terreno y la aceleración espectral para Tn = 0.2, 0.5, y 2 (s). Los resultados obtenidos de la desagregación del peligro sísmico para los sitios considerados se presentan en las Figuras 10 a 14 para periodos de retorno iguales a 125 y 475 años. En las Figuras se presenta la contribución de cada tipo de sismo (intraslab o interplaca) al peligro sísmico del sitio en estudio. Las gráficas mostradas en las Figuras consideran un intervalo de 0.5 para Mw y de 25 (km) para R. También en la Figuras se identifican escenarios sísmicos con valores modales de (Mw, R, ε).

Figura 10 Desagregación sísmica para la estación CU para: (a) a (d) periodo de retorno de 125 años; (e) a (h) periodo de retorno de 475 años 

Figura 11 Desagregación sísmica para el Sitio I (cerca de Morelia): (a) a (d) periodo de retorno de 125 años; (e) a (h) periodo de retorno de 475 años 

Figura 12 Desagregación sísmica para el Sitio II (cerca de Acapulco): (a) a (d) periodo de retorno de 125 años; (e) a (h) periodo de retorno de 475 años 

Figura 13 Desagregación sísmica para el Sitio III (cerca de Oaxaca): (a) a (d) periodo de retorno de 125 años; (e) a (h) periodo de retorno de 475 años 

Figura 14 Desagregación sísmica para el Sitio IV (cerca de Salina Cruz): (a) a (d) periodo de retorno de 125 años; (e) a (h) periodo de retorno de 475 años 

La Figura 10 muestra la desagregación sísmica para la estación CU, en donde se observa que la función de masa de probabilidad parece ser unimodal. Los parámetros modales que identifican a los escenarios sísmicos son (7.4, 150, 0.18), (7.3, 147, 0.23), (7.3, 149, 0.21), (7.8, 296, 0.34) para la aceleración máxima del terreno y las aceleraciones espectrales para Tn = 0.2, 0.5, y 2 (s), respectivamente. Estos parámetros están asociados con un periodo de retorno de 125 años. Una inspección más detallada indica que la mayoría de los escenarios sísmicos asociados con la aceleración máxima del terreno y las aceleraciones espectrales para Tn = 0.2 y 0.5 (s) provienen de las zonas sísmicas de eventos intraslab, aunque también se observa contribución de sismos interplaca para la aceleración máxima del terreno y la aceleración espectral para Tn = 0.5 (s). El escenario sísmico asociado con Tn = 2 (s) tiene su contribución más importante de eventos interplaca, más particularmente de las zonas sísmicas de Guerrero (ver Figura 1b). Los escenarios símicos identificados para un periodo de retorno de 475 años son diferentes pero no significativamente a aquellos para un periodo de retorno de 125 años.

En general, para los sitios I a IV (ver Figuras 11 a 14 ), se observa contribución de eventos interplaca e intraslab al peligro sísmico, la contribución más importante de eventos interplaca es para los Sitios II y IV (cerca de Acapulco y Salina Cruz, respectivamente), mientras que los eventos interplaca e intraslab dominan el peligro sísmico para los Sitios I y III (cerca de Morelia y Oaxaca). Los escenarios sísmicos con parámetros modales identificados en la desagregación muestran que en algunos casos la contribución al peligro sísmico es dominada por un solo evento y en otros el peligro sísmico es debido a diferentes escenarios sísmicos y diferentes tipos de temblores. Esto es importante ya que la demanda de ductilidad podría diferir para diferentes tipos de sismos (Hong et al. 2010).

Dependiendo el sitio y el periodo de retorno, el intervalo de Mw de los escenarios sísmicos identificados varía de 6.4 a 7.9, con valores de R entre 18 y 270 km, mientras que el intervalo de valores de ε varía entre 0.10 y 0.86. La variación de ε indica que la desagregación del peligro sísmico es muy sensible al valor de ε que se asocia con las relaciones de atenuación. Los resultados más detallados de la desagregación para el Sitio I se muestran en la Figura 11, en donde se presentan las funciones de masa de probabilidad y los escenarios sísmicos identificados. Se observa en esta Figura que para los periodos de retorno de 125 y 475 años, los eventos intraslab e interplaca contribuyen al peligro sísmico, aunque los eventos intraslab contribuyen con mayor frecuencia al peligro sísmico para la aceleración máxima del terreno y Tn = 0.2 (s), y en general los eventos interplaca contribuyen más al peligro sísmico para Tn = 0.5 y 2 (s). Los parámetros modales de los escenarios sísmicos identificados, para un periodo de retorno de 125 años, son (6.7, 70, 0.53), (7.5, 67, 0.22), (7.7, 185, 0.53), (7.8, 175, 0.48) para la aceleración máxima del terreno y las aceleraciones espectrales para Tn = 0.2, 0.5, y 2 (s), respectivamente. Los escenarios sísmicos identificados provienen de la zona de fuente sísmica de Profundidad intermedia Oeste Nueva (ver Figura 1c) y de las costas de Guerrero y Michoacán (ver Figuras 1b y 1c ).

La Figura 12 presenta los resultados de la desagregación sísmica para el Sitio II. La Figura muestra que el peligro sísmico es dominado por eventos interplaca para los periodos de retorno de 125 y 475 años. Los eventos sísmicos identificados para un periodo de retorno de 125 años indican que existe un sismo importante, con Mw entre a 7.5 y 7.7, R entre 30 y 45 km, y ε aproximadamente igual a 0.45. Los valores específicos de Mw, R y ε dependen del valor de T n. Los eventos sísmicos identificados para un periodo de retorno de 475 años son similares en magnitud a aquéllos para 125 años, excepto que R es menor y ε aumenta. Una inspección de los resultados generados de la simulación indica que los eventos sísmicos identificados provienen de las zonas de fuentes sísmicas de Guerrero (ver Figura 1b), como se esperaba.

La Figura 13 presenta los resultados de la desagregación sísmica del Sitio III, en donde se observa que el peligro sísmico es dominado por eventos intraslab e interplaca para los periodos de retorno de 125 y 475 años. Se observa una mayor contribución de los sismos intraslab al peligro sísmico para la aceleración máxima del terreno y para la aceleración espectral para Tn = 0.2 (s), mientras que para aceleraciones espectrales para Tn = 0.5 y 2 (s), se observa una mayor contribución de los sismos interplaca. Los escenarios sísmicos identificados para un periodo de retorno de 125 años son (6.8, 79, 0.10), (6.7, 58, 0.43), (7.7, 113, 0.74), (7.7, 133, 0.61) para la aceleración máxima del terreno y las aceleraciones espectrales para Tn = 0.2, 0.5, y 2 (s), respectivamente. Los escenarios sísmicos para un periodo de retorno de 475 años corresponden a eventos con Mw entre 6.6 y 7.9, R entre 49 y 131 km y ε entre 0.25 y 0.67. Las zonas de fuentes sísmicas asociadas con los eventos identificados corresponden a la zona de fuente sísmica de Profundidad intermedia Centro (ver Figura 1c) y las zonas de fuentes sísmicas de la costa de Oaxaca (ver Figura 1b).

La Figura 14 muestra los resultados de la desagregación sísmica del Sitio IV, donde el peligro sísmico es dominado por eventos interplaca para los periodos de retorno de 125 y 475 años. Los escenarios sísmicos identificados para un periodo de retorno de 125 años son (6.8, 35, 0.56), (7.5, 38, 0.58), (7.7, 60, 0.63), (7.6, 59, 0.52) para la aceleración máxima del terreno y las aceleraciones espectrales para Tn = 0.2, 0.5, y 2 (s), respectivamente. Los eventos sísmicos para un periodo de retorno de 475 años comprenden valores de Mw entre 6.4 y 7.8, R entre 33 y 100 km y ε entre 0.66 y 0.77. Los escenarios sísmicos identificados se asocian con las zonas de fuentes sísmicas de las costas de Oaxaca (ver Figuras 1a y 1b ).

Conclusiones

Se realizó la desagregación del peligro sísmico empleando un modelo de peligro sísmico aplicable a la región de la costa del Pacífico Mexicano y a la estación Ciudad Universitaria (estación CU) en la Ciudad de México. En el modelo se incluyen nuevas relaciones de atenuación desarrolladas para la estación CU. Con base en el modelo adoptado, se desarrollaron mapas de peligro sísmico para la región de la costa del Pacífico, y espectros de peligro uniforme (EPU) para diferentes periodos de retorno y algunos sitios seleccionados, incluyendo la estación CU. Para esta última se estudió la contribución de sismos interplaca e intraslab al espectro de peligro uniforme.

Las conclusiones principales se presentan a continuación:

  • Se mostró que el uso de la relación de atenuación desarrollada empleando todos los registros obtenidos en la estación CU para los eventos intraslab moderados y grandes conlleva a un nivel de peligro sísmico considerable. Sin embargo, este nivel de peligro sísmico es reducido drásticamente si la relación de atenuación se desarrolla con base en los registros con excitaciones mayores que 10 (Gal). Una inspección adicional de las relaciones de atenuación indica que la basada en todos los registros es inadecuada debido a que provoca que los residuos se desvíen de la normalidad, y que la desviación estándar de los residuos sea mucho mayor que aquella obtenida de la relación de atenuación con base en registros de una red de estaciones más amplia. Lo anterior simplemente enfatiza la importancia de la selección de datos para desarrollar las relaciones de atenuación y la importancia de contar con una mayor cantidad de muestras.

  • Aunque se puede desarrollar una relación de atenuación con base en los registros con excitaciones mayores que un límite inferior mayor que 10 (Gal), esto no fue considerado ya que se reduciría el número de acelerogramas que puede ser usado para desarrollar las relaciones de atenuación. La consideración del límite inferior de 10 (Gal) en el desarrollo de las relaciones de atenuación puede ser una suposición conservadora.

  • Los sismos interplaca e intraslab contribuyen al peligro sísmico en la estación CU, en particular, el peligro sísmico es dominado por eventos intraslab para periodos de vibrar menores a 0.7 (s), mientras que los eventos interplaca dominan el peligro sísmico para periodos de vibrar mayores a 0.7 (s). Esta observación es confirmada por los resultados de desagregación del peligro sísmico.

  • Los valores del espectro de diseño para la Ciudad de México son significativamente mayores a los valores del espectro de peligro uniforme para un periodo de retorno de 125 años. Se observa también que los valores del espectro de diseño para T n entre 0.25 y 1.75 (s) son mayores a los valores del espectro de peligro uniforme para un periodo de retorno de 475 años.

  • Los resultados de la desagregación, basados en el modelo de peligro sísmico adoptado, indican que ésta presenta escenarios sísmicos identificados en donde la contribución al peligro sísmico para un periodo de retorno de 125 o 475 años es dominada por un solo evento y en otros casos el peligro sísmico es debido a diferentes escenarios sísmicos y tipos de temblores. La identificación del tipo de temblor para escenarios sísmicos puede ser importante ya que las estadísticas de demanda de ductilidad pueden diferir para diferentes tipos de temblores. Los escenarios identificados dependen del periodo fundamental de vibrar de la estructura considerada. En general, como es de esperar, para los sitios en la región de la costa del Pacífico Mexicano, el peligro sísmico es dominado por sismos interplaca, mientras que para los sitios alejados de la costa del Pacífico, los sismos intraslab e interplaca contribuyen al peligro sísmico

Reconocimientos

El apoyo financiero recibido del Consejo Nacional de Ciencia y Tecnología (CONACyT) de México, The Natural Science and Engineering Research Council of Canada y del Instituto de Ingeniería de la UNAM es ampliamente agradecido. Agradecemos de manera especial los comentarios constructivos y sugerencias del Dr. J. A. Escobar y de dos revisores anónimos. También agradecemos a la Coordinación de Instrumentación Sísmica del IIUNAM por proveer algunos registros sísmicos empleados en este trabajo.

REFERENCIAS

Arroyo, D, D García, M Ordaz, MA Moray y SK Singh (2010), "Strong ground-motion relations for Mexican interplate earthquakes", Journal of Seismology, vol. 14, pp. 769-785. [ Links ]

Atkinson GM y DM 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 ]

Atkinson GM (2006), "Single-station sigma", Bulletin of the Seismic Society of America, vol. 96, pp. 446-455. [ Links ]

Baker JW, CA Cornell y P Tothong (2005), "Deaggregation of seismic drift hazard", in Proc.9th International Conference on Structural Safety and Reliability (ICOSSAR05), Rome, Italy. [ Links ]

Bazzurro P y CA Cornell (1999), "Disaggregation of seismic hazard", Bulletin of the Seismic Society of America, vol. 89, pp. 501-520. [ Links ]

Boore DM, WB Joyner y TE Fumal (1997), "Equations for estimating horizontal response spectra and peak acceleration from western North America", Seismological Research Letters, vol.68, pp.128-153. (for Erratum: see 2005, Seismological Research Letters, vol.76, pp. 368-369). [ Links ]

Chapman MC (1995), "A probabilistic approach to ground-motion selection for engineering design", Bulletin of the Seismic Society of America, vol. 85, pp. 937-942. [ Links ]

Cornell CA (1968), "Engineering seismic risk analysis", Bulletin of the Seismic Society of America, vol. 58, pp. 1583-1606. [ 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, México. [ 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, SK Singh, M Herraiz, M Ordaz, y JP 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-Pérez J, F Castellanos, O Díaz (2005), "Occupancy importance factor in earthquake engineering", Engineering Structures, vol. 27, pp. 1625-1632. [ Links ]

García-Soto, AD, A Pozos-Estrada, HP Hong y R Gómez-Martínez (2012), "Estimación del peligro sísmico debido a sismos interplaca e inslab y sus implicaciones en el diseño símico", Revista de Ingeniería Sísmica, No. 86, pp. 27 - 54. [ Links ]

Harmsen S, D Perkins y A Frankel (1999), "Deaggregation of probabilistic ground motions in the central and eastern United States", Bulletin of the Seismic Society of America, vol. 89, pp. 1-13. [ Links ]

Hong HP(2012), "Peligro, riesgo, diseño sísmico óptimo", The Mexican Academy of Engineering, Mexico. [ Links ]

Hong HP y K Goda (2006), "A comparison of seismic hazard and risk deaggregation", Bulletin of Seismological Society of America, vol. 96(6), pp. 2021-2039. [ Links ]

Hong HP, AD Garcia-Soto y R Gomez (2010), "Impact of different earthquake types on the statistics of ductility demand", Journal of Structural Engineering, ASCE, vol. 136(7), pp. 770-780. [ Links ]

Hong HP, Y Zhang y K Goda (2009a), "Effect of the spatial correlation on the estimated ground motion prediction equations", Bulletin of the Seismological Society of America, vol. 99(2A), pp. 928-934. [ Links ]

Hong HP, A Pozos-Estraday R Gomez (2009b), "Orientation effect on ground motion measure for Mexican subduction earthquakes", Earthquake Engineering and Engineering Vibration, vol. 8(1), pp. 1-16. [ Links ]

Hong HP, K Goday AG Davenport (2006), "Seismic hazard analysis: a comparative study", Canadian Journal of Civil Engineering, vol. 33(9), pp. 1156-1171. [ Links ]

Hong HP y E Rosenblueth (1988), "The Mexico Earthquake of September 19, 1985 - Model for Generation of Subduction Earthquakes", Earthquake Spectra, vol. 4(3), pp. 481 - 498. [ Links ]

Jaimes MA, y E Reinoso (2006), "Comparación de 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 WB y DM Boore (1993), "Methods for regression analysis of strong-motion data", Bulletin of the Seismic Society of America, vol. 83, pp. 469-487. [ Links ]

McGuire RK (1995), "Probabilistic seismic hazard analysis and design earthquakes: closing the loop", Bulletin of the Seismic Society of America, vol. 85, pp. 1275-1284. [ Links ]

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

Montalvo-Arrieta JC (2002), "La respuesta sísmica del terreno firme en la ciudad de México. Observaciones y Modelos", Tesis Doctoral, Universidad Nacional Autónoma de México, Instituto de Geofísica, Posgrado en Ciencias de la Tierra. [ Links ]

Musson RMW (2000), "The use of Monte Carlo simulations for seismic hazard assessment in the UK". Annali di Geofisica, vol. 43, pp. 1-9. [ Links ]

Nishenko SP y SK Singh (1987), "The Acapulco-Ometepec, México earthquake of 1907-1982: Evidence for a variable recurrence history", Bulletin of the Seismic Society of America, vol. 77, pp. 1359-1367. [ Links ]

NTCDS (2004), "Normas Técnicas Complementarias para Diseño por Sismo", , Gaceta Oficial del Departamento del Distrito Federal, México, 6 de octubre de 2004. [ Links ]

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

Ordaz M, E Miranda y J Avilés (2000), "Propuesta de espectros de diseño por sismo para el DF", Memorias del XII Congreso Nacional de Ingeniería Estructural, León, CD, noviembre. [ Links ]

Pozos-Estrada A, R Gomezy HP Hong (2008), "Seismic design: benefit/cost for overall service time versus per unit service time", 14 Word Conference on Earthquake Engineering, Beijing, China. [ Links ]

Reinoso E y MA Jaimes (2009), "Criterios para obtener acelerogramas de diseño en sitios afectados por varias fuentes sísmicas usando como ejemplo el caso de terreno firme de la ciudad de México", Revista de Ingeniería Sísmica, No. 81, pp. 1 - 18. [ Links ]

Reyes C, E Miranda, M Ordazy 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 M Ordaz (1987), "Use of seismic data from similar regions", Earthquake Engineering &. Structural Dynamics, vol. 15, pp. 619-634. [ Links ]

Rosenblueth E y JM 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 ]

Rosenblueth E, M Ordaz, JF Sánchez-Sesma y SK Singh (1989), "The Mexico Earthquake of September 19, 1985 - Design Spectra for Mexico's Federal District", Earthquake Spectra, vol. 5(1), pp. 273 - 291. [ Links ]

Singh SK, M Rodriguez y L Esteva (1983), "Statistics of small earthquakes and frequency of occurrence of large earthquakes along the Mexican subduction zone", Bulletin of the Seismic Society of America, vol. 73, pp. 1779-1796. [ Links ]

Singh SK, M Ordaz, X Pérez-Campos y A Iglesias (2013), "Intraslab versus Interplate Earthquakes as Recorded in Mexico City: Implications for Seismic Hazard", Earthquake Spectra, doi: http://dx.doi.org/10.1193/110612EQS324M. [ Links ]

Recibido: 15 de Abril de 2014; Aprobado: 27 de Octubre de 2014

Creative Commons License This is an open-access article distributed under the terms of the Creative Commons Attribution License