Implications of the intensity centers location of the caracas 1967 and the caracas 1812 subevent

Choy, José E; Palme, Christl; Guada, Carlos; Klarica, Stéphanie; Morandi, María

Servicios Personalizados

Revista

Articulo

Indicadores

Citado por SciELO
Accesos

Links relacionados

Similares en SciELO

Otros
Otros

Permalink

Revista de la Facultad de Ingeniería Universidad Central de Venezuela

versión impresa ISSN 0798-4065

Rev. Fac. Ing. UCV v.26 n.2 Caracas jun. 2011

Implicaciones de la localización de los centros de intensidad del subevento de caracas del terremoto de 1812 y del sismo de caracas de 1967

José E. Choy, Christl Palme, Carlos Guada, Stéphanie Klarica, María Morandi.

Universidad de los Andes. Facultad de Ciencias. Departamento de Física. Laboratorio de Geofísica. Mérida. Venezuela. e-mail: choy@ula.ve; palme@ula.ve; carlosg@ula.ve; klarica@ula.ve; maria@ula.ve

RESUMEN

En este trabajo se relocalizan los centros de intensidades del sismo de Caracas de 1967 y del subevento de Caracas de 1812 utilizando el método de Bakun & Wentworth, incorporando información sobre el rumbo del fallamiento e incertidumbres en los valores de intensidades. La relación entre la ubicación de ambos centros, arroja información importante para entender la relación entre ambas rupturas. Los parámetros del centro de intensidades del evento de 1967, son consistentes con la información instrumental. La literaratura especializada mostró que la ruptura del sismo de Caracas de 1967 tuvo al menos cuatro subeventos. Posteriormente otros autores refinan la localización del epicentro (10,558ºN y 67,31ºO) y lo colocan virtualmente sobre la falla de San Sebastián. El centro de intensidades determinado en este trabajo (10,60ºN y 67,20ºO) se ubica entre los dos primeros subventos, lo cual tiene sentido porque ellos liberaron el 72% del momento sísmico total. Los eventos tercero y cuarto no sólo fueron de menor tamaño sino que además ocurrieron más alejados mar adentro y, por ello, tuvieron menor influencia en la distribución de las intensidades. Esto también explica por qué la magnitud momento obtenida a partir de las intensidades (MWI = 6,4 ± 0,3) es algo menor que la magnitud momento obtenida por los especialistas (MW = 6,6). El centro de intensidades para el subevento de Caracas de 1812, se ubica en 10,6ºN y 67,1ºO sobre la falla de San Sebastián y entre los dos primeros subeventos de 1967. La magnitud momento MWI es igual a 7,1 ± 0,33. Esto indica que el subevento de Caracas de 1812 rompió la falla de San Sebastián sin descartar que, al igual que en el sismo de 1967, hayan ocurrido rupturas en otras fallas. Por otra parte, es probable que en 1812 se haya roto el segmento, o partes de él, entre los dos primeros subeventos de 1967 y otros segmentos al oeste de los mismos.

Palabras clave: Caracas, 1967, 1812, Terremoto, Intensidades, Sismicidad, Histórica, Epicentro.

Implications of the intensity centers location of the caracas 1967 and the caracas 1812 subevent

ABSTRACT

In this work, intensity centers of Caracas 1967 and Caracas 1812 sub-event are relocated using the Bakun & Wentworth method with two modifications. The first one introduces intensity intervals to account for uncertainties of the historical information, and the second one uses a priori information of the fault strike. The relative locations between those events give important information to understand the relationship between both ruptures. The intensity center parameters of the 1967 event are consistent with instrumental information. The specialized literature showed that the rupture had at least four subevents. Later, other workers relocate the epicenter (10.558ºN y 67.31ºW) very close to the San Sebastián fault. The intensity center determined in this work (10.60ºN y 67.20ºW) lies between the first two subevents. This make sense as these two subevents released approximately 72% of the total moment produced by the earthquake and they are the ones closer to the damage areas; the other two subevents are not only smaller, but also farther away from the coast and therefore they had less influence on the intensity distribution. This also explains why the moment magnitude obtained from intensities (MWI = 6.4 ± 0.3) is smaller than the one obtained by other specialists (MW = 6.6). The Intensity center obtained for the Caracas 1812 subevent is: 10.6ºN and 67.1º W with magnitude MWI = 7.1 ± 0.33; very close to the San Sebastián fault and in between the first two subevents of 1967. This shows that the Caracas 1812 subevent ruptured the San Sebastián fault and perhaps, as in the case of the 1967 event, also ruptured other faults. On the other hand is very likely that broke the segment located between the first two subevents of 1967 and other segments to the west of the first subevent.

Keywords: Caracas, 1967, 1812, Earthquake, Intensities, Seismicity, Historical, Epicenter

Recibido: octubre de 2009 Recibido en forma final revisado: julio de 2011

INTRODUCCIÓN

Diversos autores (Fiedler, 1961; Centeno Graü, 1969; Fiedler, 1972; Audemard, 2002; Altez, 2005 y 2006) han considerado al sismo del 26 de marzo de 1812 como un sismo multifoco. Centeno Graü (1969) considera que tuvo tres focos: uno en el mar Caribe, entre los Roques y la costa de La Guaira; otro frente a San Felipe, Barquisimeto y El Tocuyo; y un tercero, al sur del Lago de Maracaibo. Fiedler (1961, 1972) también apoya la hipótesis de que tuvo tres focos: uno al norte de la costa de la Guaira, un segundo foco en la falla de Boconó, entre San Felipe y Barquisimeto, y un tercero al suroeste de Mérida. Altez (1998, 2005, 2006) sostiene que en realidad tuvo dos focos, porque considera que en Mérida, el movimiento sísmico ocurrió aproximadamente una hora después de los de Caracas y San Felipe.

Fiedler (1961) asigna una magnitud de ondas superficiales Ms = 6,3 al sismo de Caracas; Ms = 6,2 al de San Felipe y de 7,0 al de Mérida. Posteriormente (Fiedler, 1972) realiza una reevaluación y asigna una magnitud de ondas de cuerpo mC = 7,1 a cada uno de los tres eventos, que corresponde a una magnitud Ms = 7,3 (Richter, 1958). Grases (1990) estima en 7,0; 7,2 y 6,3 las magnitudes para las rupturas en Mérida, Barquisimeto-San Felipe y Caracas, respectivamente. Grases&Rodríguez (2001) dan valores de 6,5 - 6,7; 6,9 - 7,2 y 7,1 para los eventos en Mérida, San Felipe y Caracas, respectivamente. Audemard (2002) quien estudia las rupturas de los grandes sismos históricos a partir del análisis de la sismicidad contemporánea, propone para el subevento de Caracas una ruptura de la falla de San Sebastián con un foco que se encuentra entre Puerto Cabello y La Guaira, al oeste del propuesto por Fiedler, sin hacer inferencias sobre su magnitud. Para el segundo foco propone como fuente la terminación norte de la Falla de Boconó a lo largo del valle de Yaracuy que se extiende desde la cuenca de tracción de Cabudare hacia el nordeste. La magnitud de este segundo subevento la estima en aproximadamente 7,0 a partir de la extensión de la isosista VIII. Laffaille&Ferrer (2003) consideran que la ruptura de Mérida tuvo una magnitud 5,1 – 5,4 Ml y que la falla del Albarregas (Laffaille&Ferrer, 2005) podría haber estado asociada con el evento sísmico.

El sismo de Caracas en 1967 fue registrado instrumentalmente por las estaciones de la Red Mundial y existen varios estudios de sus procesos focales. Fiedler (1968) da como coordenadas del epicentro 11,0°N y 66,25°O con una magnitud de 6,3 mC, y opina que el epicentro dado por el USGS (United State Geological Survey) (10,56°N y 67,26°O) no es acertado porque la precisión del USGS en ese tiempo no era suficiente a escala local. Además, de acuerdo con Fiedler, un epicentro como el del USGS no coincide con la distribución de los daños observados. Sin embargo, el Servicio Geológico de Moscú (MOS) da la misma localización, y la del ISC también está sólo ligeramente al nor-oeste (10,68°N y 67,40°O, mb 5,5; todos estos datos pueden ser consultados en el portal del International Seismological Centre (ISC). La única localización que coincide con la de Fiedler es la de la Universidad West Indies (TRN: 11,3°N y 67,0°O). Para explicar la dispersión de las localizaciones de los epicentros, Fiedler supone que el evento consistió en dos sismos que ocurrieron casi simultáneamente, el primero a 70 km del Observatorio Cajigal de magnitud mC 5,7 y el segundo a 30 km del observatorio de magnitud mC 6,3. En 1990 Suárez&Nabélêk (1990) utilizan un método de inversión de los sismogramas de ondas corpóreas para obtener un modelo de ruptura que consiste en cuatro focos, la mayor parte de la energía sería liberada por los dos primeros que se ubican en una línea paralela a la costa del Litoral Central con un mecanismo focal rumbo-deslizante dextral si se considera como plano de falla aquel plano nodal que está paralelo a la misma línea que une los dos subeventos. La distancia entre los dos subeventos sería aproximadamente 50 km. La magnitud de momento que obtienen es Mw 6,6. Estos autores no presentan una localización absoluta del evento, solamente dan una localización de los cuatro subeventos relacionada con el epicentro publicado por el ISC. Engdahl et al. (1998) publicaron la relocalización de los hipocentros de casi 100000 eventos sísmicos de todo el mundo, entre ellos el sismo de Caracas de 1967 (estos datos están accesibles en el boletín en línea del ISC) obteniendo las coordenadas 10,55°N, 67,31°O y 22 km profundidad, un resultado muy cercano al que fue publicado originalmente por el USGS.

Para las evaluaciones del terremoto de 1812 todos los autores, salvo Audemard (2002), han utilizado isosistas, lo cual implica inferencias visuales, que potencialmente pueden introducir sesgos y hacen difícil realizar verificaciones estadísticas (Gasperini et al. 1999). Para tratar de obtener inferencias más objetivas a partir de intensidades de sismos históricos. Bakun&Wenworth (1997; B&W en adelante) desarrollaron un método estadístico de inversión que utiliza directamente las observaciones individuales de intensidades en la determinación de epicentros y magnitudes. El método requiere de un proceso previo de calibración entre intensidades, magnitudes y distancias epicentrales, lo cual se logra a través de sismos modernos registrados instrumentalmente. Presenta ciertas ventajas con respecto al tradicionalmente utilizado método de las isosistas: (1) Es objetivo en el análisis de las intensidades, sobre todo cuando la distribución de las mismas no es uniforme o su número es escaso. (2) Todos los datos de intensidades contribuyen con el mismo peso al cálculo de la magnitud y del epicentro, y no solamente los cercanos al foco. (3) El tratamiento estadístico de los datos permite establecer intervalos de confianza.

Palme et al. (2005a) llevaron a cabo el proceso de calibración para el Centro Occidente de Venezuela. Las relaciones obtenidas fueron utilizadas por Palme et al. (2005b) para analizar la información de intensidades asociada con el evento de Mérida, y obtuvieron un centro de intensidades entre las ciudades de Mérida y Tabay, y una magnitud de momento MWI = (6,5 ± 0,5). Este valor de magnitud solapa con el intervalo obtenido por Laffaille&Ferrer (2003) porque convirtiendo sus valores a magnitud de momento, se obtiene MWI = (5,6 – 6,0), aunque el intervalo de valores probables obtenido por Palme et al. (2005b) también acepta magnitudes algo mayores.

En años recientes se ha incrementado notablemente la disponibilidad de información sobre los daños causados por los subeventos de Caracas y de San Felipe en 1812, debido a las publicaciones de Altez (2006), Grases et al. (1999) y Palme et al. (2005b), lo cual permitió incrementar el número de localidades en las que se les puede asignar intensidades y, además, establecer rangos de incertidumbres para las intensidades previamente reportadas.

En este trabajo se utiliza el método de B&W para procesar las intensidades y sus rangos de incertidumbre con el objeto de obtener epicentros y magnitudes del sismo de 1967 y del subevento de Caracas de 1812. Los resultados obtenidos son bastante robustos, porque para diferentes escenarios (diferentes intensidades asignadas por diversos autores a algunas localidades, o la inclusión o no de algunas localidades en el análisis) los resultados no cambian de manera significativa.

MÉTODO

Palme et al. (2005a) utilizaron eventos del centro-occidente de Venezuela para calibrar la relación entre magnitud momento, MWI, y la intensidad MSK, IMSK, aplicando el método propuesto por B&W. El resultado de dicha calibración es:

I_MSK = -2.2237+1.6684*M_WI–0.04121*(med Δ) (1)

En donde: M_WI es la magnitud momento calculada a partir de intensidades y calibrada con magnitudes momento determinadas a partir de sismogramas, y (med Δ) es la mediana de las distancias entre el centro de intensidades y los lugares en donde se ha producido una intensidad I_MSK dada. Como se demostró en Palme et al. (2005a), esta relación produce buenos resultados para distancias de hasta aproximadamente 150 km. Para distancias mayores, la relación distanciaintensidad deja de ser lineal para la región centro-occidental de Venezuela. Por ello, en el proceso de inversión sólo se toman en cuenta intensidades de lugares que se encuentren a distancias de hasta aproximadamente 150 km. Por otra parte, no se llevaron a cabo correcciones por efectos de sitio, ya que el número de lugares con valores de intensidad es muy pequeño como para una evaluación estadística de tales efectos y los estudios sistemáticos en Venezuela de la influencia de los suelos sobre las intensidades apenas están en sus inicios (Rocabado et al. 2005; Schmitz et al. 2005; Altez & Laffaille, 2006). Para evitar sesgos en la ecuación (1) se excluyeron de la calibración aquellos lugares con notorios efectos de sitio, como por ejemplo, Los Palos Grandes, durante el sismo de Caracas de 1967. Es interesante notar que este lugar no estaba habitado en 1812.

En este trabajo se utilizan intensidades EMS-98 (European Macroseismic Scale 1998; Grünthal, 1998). La determinación de los centros de intensidades y las magnitudes momento- intensidad, M_WI, se lleva a cabo siguiendo el procedimiento dado en B &W y en Hinzen & Oemisch (2001). Primero se construye una malla de epicentros tentativos (centros de intensidades) alrededor de la zona epicentral. Denotemos por j a un punto arbitrario de la malla y denotemos por k a una localidad para la cual se conoce la intensidad. Utilizando la ecuación (1), realizamos el siguiente cálculo para cada punto de la malla:

M_WIjk = [I_EMSk+2.2237+0.04121*(Δ_jk)]/1.6684 (2)

En donde: I_EMSk es la intensidad en el k-ésimo lugar y Δ_jk su distancia hasta el j-ésimo punto de la malla. Luego se calcula el valor medio de los M_WIjk (que llamaremos M_WIj):

M_WIj = media (M_WIjk) (3) y

rms(M_WIj) rms(M_WIj-M_WIjk)–rms0(M_WIj–M_WIjk) (4)

donde: rms0(M_WIj–M_WIjk) es el valor mínimo de rms(M_WIj- M_WIjk). El centro de intensidades corresponde al punto de la malla donde rms(M_WIj) = 0. La magnitud momento-intensidad, M_WI, corresponde al M_WIj del centro de intensidad. Los niveles de confianza para el rms y para los valores de magnitud se tomaron de B &W (1999).

CORRECCIÓN POR EL RUMBO DE LA FALLA PARA SISMOS GRANDES

Las ecuaciones del método de B&W (1997) sólo dependen de las distancias entre los lugares con intensidades y el centro de intensidades, y son independientes del rumbo de la ruptura. Sin embargo, sismos de magnitud 7 pueden producir rupturas sub-superficiales del orden de los 50 km (Wells & Coppersmith, 1994). En consecuencia, las distribuciones de intensidades no se pueden aproximar con isosistas circulares. Cuando los puntos de intensidades están distribuidos isotrópicamente con respecto a la falla, el método de B&W funciona bien, a pesar del mencionado efecto del rumbo de la falla, porque promedia sobre todas las direcciones. Sin embargo, en los casos que presentan un intervalo grande de acimuts (~180º) alrededor del epicentro sin ningún o muy pocos lugares que tengan intensidades asignadas, el método tiende a alejar el epicentro hacia la región donde faltan puntos de intensidades (Choy et al. 2007), lo cual podría pasar con la mayoría de los sismos históricos que han ocurrido cerca de Caracas, porque la principal estructura sismogénica, la falla de San Sebastián, se ubica esencialmente en la línea de la costa y hay muy pocos datos de intensidad de las islas cercanas.

Para tratar de corregir este efecto de la asimetría en la distribución de las intensidades, Choy et al. (2007) proponen el uso de una función de ponderación de las intensidades, que depende del rumbo de la falla pero no de su ubicación. Esto es: se considera un conjunto tentativo de fallas paralelas que cubran el área epicentral. Cada una de esas fallas se asocia con un epicentro tentativo. La función error se calcula al igual que en B&W (1997) pero la función adicional de ponderación depende de la distancia desde cada centro tentativo de intensidades hasta la traza de la falla. Se escogió una forma exponencial para la función peso, dada por:

pjk = wjk*exp(-c*djk) (5)

Donde: w_jk es la función de ponderación dada en B&W (1997), c es un parámetro determinado empíricamente en Choy et al. (2007) y d_jk es la distancia desde la falla que pasa por el j-ésimo centro tentativo de intensidades hasta la k-ésima localidad con valor conocido de intensidad. Por consiguiente, la función error se escribe como en B&W:

rms(M_WIj –M_Wijk)={å_k[(M_WIj–M_Wijk)*p_jk]²/Pj}^1/2 (6)

con:

P = å_k (p_jk)² (7)

Para sismos de magnitud ≈ 7, un valor de c entre 0,03 km^-1 y 0,035 km^-1 produce buenos resultados.

UN CASO DE PRUEBA: EL SISMO DE LOMA PRIETA

Para validar este enfoque, utilizamos el sismo de Loma Prieta. B&W (1999) aplicaron su método a este evento, cuyo epicentro instrumental se ubica cerca de la costa de California, y la mayoría de los 560 puntos de intensidad se encuentran al este de la falla de San Andrés. Las intensidades se obtuvieron del sitio web de Bakun (http://quake.wr.usgs.gov/~bakun/, última consulta en 24-05-2008). Para realizar los cálculos, se utilizaron las ecuaciones publicadas en BW (1997), válidas para California. En una primera fase se corrió el programa sin la corrección de rumbo de falla, y se obtuvo 37ºN y 122,1ºO para el centro de intensidad, y M_WI = 6,92; lo que coincide con los resultados de B&W (1999). Luego se corrió el programa incluyendo la corrección de rumbo, suponiendo un rumbo de falla de 130º en sentido horario desde el norte. De esta manera se obtuvo un centro de intensidad en 37,05ºN y 121,95ºO, el cual es muy cercano al epicentro determinado instrumentalmente: 37,04ºN y 121,88ºO (Figura 1), y una magnitud 6,8 M_WI, la cual es consistente con la determinada instrumentalmente, de 6,9 M_W.

PROCEDIMIENTO PARA INCORPORAR INTERVALOS DE INTENSIDADES EN EL MÉTODO DE BAKUN Y WENTWORTH

En primer lugar se establece una tabla con las coordenadas de los lugares y los valores máximos y mínimos de intensidades que fueron reportados en la literatura, pero si en un lugar se tiene un solo valor, éste se utiliza como máximo y mínimo. Con esta información se procede de la siguiente manera: se debe suponer que para un determinado sismo se dispone de intensidades correspondientes a N sitios de los cuales a M < N sitios se le han asignado dos valores de intensidad por dos o más autores distintos y en N – M los diferentes autores coinciden en un único valor. Se realiza una permutación aleatoria de los M sitios y a los primeros P se les asignan los valores máximos de intensidades reportados, en donde: P es un número aleatorio que cumple con la condición 0 ≤ P ≤ M. A los restantes M – P se les asignan los valores mínimos. Se aplica de manera usual el método a las intensidades resultantes para los N sitios. El proceso se repite un número suficientemente de veces como para evaluar el efecto de las incertidumbres en las intensidades sobre la incertidumbre en la determinación del centro de intensidades y de la magnitud. En cada repetición se realiza una nueva permutación de los M sitios y se selecciona un nuevo valor aleatorio de P. Los valores finales se obtienen tomando el promedio rms de todas las repeticiones.

DATOS

Subevento de Caracas de 1812

Para el sismo de 1812 se realizó un análisis de los datos de intensidad que han sido publicados. Estos datos fueron tomados del STSHV (Sistema de Teleinformación de la Sismología Histórica de Venezuela, 2005). En algunos casos encontramos diferencias entre los valores que han sido asignados por los diferentes autores, entonces volvimos a las descripciones históricas de Altez (1998, 2005, 2006) e hicimos nuestro propio análisis. Además, haciendo uso del STSHV pudimos estimar intensidades en algunos lugares adicionales. Muchas de las intensidades publicadas se refieren a la escala MMI (Modified Mercalli Intensity) mientras los datos más recientes y los estimados por nosotros se refieren a la escala EMS-98 (Grünthal, 1998). Los grados de ambas escalas son equivalentes, la diferencia está en las descripciones cualitativas de los grados particulares. En algunos casos complementamos para las nuevas asignaciones las descripciones de la escala EMS-98 con las de la escala MMI y con los efectos ambientales descritos en la nueva escala ESI-07 (INQUA Environmental Seismic Intensity Scale, 2007).

En este proceso confrontamos cuatro problemas que son comunes en el análisis sismológico de los datos históricos.

Primero, en las fuentes históricas se tiende a sobreestimar los daños ocurridos en un terremoto, muchas veces por el interés de los afectados de obtener ayuda económica para la reconstrucción de los edificios destruidos o también por interés político. Este fenómeno fue analizado con detalle por Laffaille et al. (2007). Para comprender el contexto histórico de Venezuela en el cual se desarrolló el terremoto de 1812 fue de suma utilidad el libro publicado por Altez (2006). En caso que todavía tuviéramos dudas sobre un valor en particular optamos por definir un intervalo de intensidades, porque el método utilizado permite procesarlo.

El segundo problema surge por las características de la escala EMS-98, que requiere de una clasificación de la vulnerabilidad de los edificios. Las construcciones comunes en la época de 1812 eran de adobe, bahareque o tapias. Centeno Graü (1900) comenta acerca del comportamiento de estas construcciones durante el terremoto del 29 de octubre de 1900: ... las de adobe y tapial son los que peores se comportan; las de bahareque aún las muy viejas resisten bien. La escala EMS-98 clasifica las construcciones de adobe como A, la clase de mayor vulnerabilidad. Sobre las construcciones de tapias, dice Castillo (1983) que presentan muy poca resistencias a fuerzas de cizallas y que, gracias a sus pesados techos de tejas, son muy propicias a colapsar. Así consideramos tanto las construcciones de adobe como de tapia como vulnerables de acuerdo a la clase A de la escala EMS-98. Si bien las casas de bahareque en algunos casos parecen resistir mejor, esto no puede ser generalizado. De acuerdo con el informe de FUNVISIS (1997) acerca del terremoto de Cariaco, la mayoría de las casas colapsadas y aquellas que fueron fuertemente dañadas eran de bahareque. Sin embargo, hay que tomar en cuenta también que difícilmente las casas de bahareque causan víctimas mortales. Como no conocemos los detalles de la calidad de construcción de las casas de bahareque, optamos por clasificarlas igualmente como A para evitar una sobreestimación de las intensidades. Muchos de los datos de intensidad del sismo de 1812 se refieren a edificios monumentales (iglesias y edificios públicos). Probablemente, la construcción de estos edificios fue de mejor calidad porque ciertos elementos de reforzamiento fueron utilizados, sin embargo hay que considerar de acuerdo con Grünthal (sección 3.5, 1998) ...the structural and non-structural complexity of such buildings is such that they may be more likely to be damaged than ordinary buildings, even though they may be better built. También el comentario de Alejandro Ibarra (citado por Altez, 2006, p.325) ...la calidad de aquellas construcciones fue altamente responsable de los daños producidos por el terremoto. Si suponemos que la calidad de las construcciones en los pueblos no fue mejor que en Caracas podemos llegar a la conclusión de que la vulnerabilidad de los edificios monumentales era también muy alta. Yamazaki et al. (2005) proponen un modelo para Caracas de 1812 con un 75% de los edificios de vulnerabilidad clase A y 25% de clase B. Nosotros suponemos que en el interior del país el porcentaje de edificios de alta vulnerabilidad era aún mayor.

El tercer problema que se presenta es el hecho de que en muchos sitios disponemos solamente de descripciones que se refieren a los daños de uno o pocos edificios monumentales, más que todo en iglesias. En estos casos asignamos intervalos de intensidad.

El cuarto problema se relaciona con los efectos de sitios. A pesar de que en los últimos años se está realizando un gran esfuerzo en Venezuela por realizar estudios sistemáticos de las respuestas de los suelos, estos estudios se están haciendo solamente en las grandes ciudades y no tenemos información en las poblaciones pequeñas. Son famosos los efectos de sitio observados en Caracas en el terremoto de 1967 (Los Palos Grandes, San Bernardino); sin embargo, estos lugares no estaban poblados en 1812. Yamazaki et al. (2005) publican un mapa de distribución de intensidades en el sismo de 1812 en Caracas que muestra valores entre VII+ y IX MMI. Estas variaciones se tomaron en cuenta como se observa en la Tabla 1. En otros lugares, a falta de más información no se tomaron en cuenta efectos de sitio para nuestro estudio y en su lugar se utilizaron intervalos de incertidumbre para las intensidades. Es importante mencionar que los autores del presente trabajo modificaron el método de B&W para tomar en cuenta este tipo de situaciones, al considerar que el uso de intervalos de intensidades en lugar de valores únicos disminuye la influencia de posibles efectos de sitio sobre los resultados de la inversión. De todos modos B&W (1997) comparan los resultados de datos que fueron procesados con correcciones de sitio con aquellos que fueron procesados sin tomar en cuenta estas correcciones. Ellos llegaron a la conclusión de que en caso de no considerar efectos de sitio, el área limitada por un determinado nivel de confianza aumenta significativamente, pero los resultados siguen siendo buenos.

En el portal de Guada (2009) se han publicado los datos de las intensidades utilizadas en este trabajo para el subevento de Caracas de 1812 conjuntamente con los valores propuestos por otros autores. También se publicaron en este sitio las citas más importantes y sus referencias que justifican la asignación de los valores de intensidad. Los valores utilizados por nosotros se pueden ver en la Tabla 1 y en la Figura 2. Para representar los valores de intensidad y los intervalos en las figuras se utiliza un símbolo de reloj que puede indicar un solo valor o un intervalo como se muestra en la Figura 3.

Sismo de Caracas de 1967

Para las intensidades del terremoto de 1967 en Caracas se utilizaron dos fuentes independientes, la publicación de Suárez&Nábelek (1990) y los datos publicados en la base de datos de CERESIS (Centro Regional de Sismología de América del Sur). Combinamos las dos listas utilizando los valores mínimos y máximos en cada sitio. El resultado de esta fusión se puede ver en la Tabla 2 y la Figura 4.

RESULTADOS

Tanto para el sismo de 1812 como para el de 1967, se calculó la función dada en la ecuación (5) suponiendo que las rupturas fueron paralelas a la falla de San Sebastián (Audemard, 2002). El centro de intensidades hallado para el evento de 1812 fue 10,60°N y 67,10°O y la magnitud obtenida fue 7,1 ± 0,33 MWI. En la Figura 5 se muestra este resultado junto con el de Fiedler (1961). Para el evento de 1967, se obtuvo un centro de intensidades en: 10,60°N y 67,20°O y una magnitud de 6,4 ± 0,3 MWI. El resultado de nuestra localización se puede ver en la Figura 6a. En la misma figura se muestra el resultado de la relocalización hecha por Engdahl et al. (1998); las coordenadas de su epicentro son 10,558°N y 67,31°O, alrededor de trece kilómetros al oeste de nuestro centro de intensidades.

DISCUSIÓN

Con esta relocalización del sismo de 1967 los primeros dos subeventos determinados por Suárez&Nábêlek (1990) se ubican virtualmente sobre la traza de la falla de San Sebastián y nuestro centro de intensidades entre estos dos subeventos (Figura 6b). Esto demuestra que nuestros resultados son consistentes con los de Engdahl et al. (1998) y los de Suárez&Nábêlek. Es importante notar que el método de B&W (1997) ubica el centro de intensidades, el cual para sismos pequeños prácticamente coincide con el epicentro. Sin embargo, para sismos grandes puede haber una separación entre ambos, dependiendo de la longitud de la ruptura y de las características de la misma (unilateral, bilateral o multifoco). Por ello el centro de intensidades no coincide con el epicentro sino que se ubica entre los dos primeros subeventos de Suárez&Nábêlek porque ellos liberaron el 72% del total de la energía sísmica y fueron los más cercanos a las áreas dañadas. Los otros subeventos no sólo fueron menores sino que además ocurrieron costa afuera y por consiguiente tuvieron una influencia mucho menor en las intensidades. Esto también es consistente con una magnitud MWI ligeramente menor que la magnitud MW, aunque es de notar que ambos valores son compatibles si se toma en cuenta los intervalos de incertidumbre.

Los eventos de Loma Prieta 1989, Caracas 1967 y Caracas 1812, tienen en común una asimetría en la distribución de las intensidades con respecto a la ruptura, debido a que sus epicentros se ubican cerca de la costa. Hemos verificado en el presente trabajo que el método de B&W (1997) con las modificaciones de Choy et al. (2007) produjeron excelentes resultados en los dos primeros casos. Por ello creemos que los resultados aquí obtenidos para 1812 también deben ser confiables. La magnitud MWI de 7,1 ± 0,33 es consistente con los valores determinados y/o utilizados por otros autores (Audemard, 2002; Yamazaki et al. 2005). El centro de intensidades se ubica casi al norte de Carayaca y muy cerca del centro de intensidades del evento de 1967 y por lo tanto es consistente con una ruptura en la falla de San Sebastián, como lo propuso Audemard. Durante el sismo de 1967, no hubo ruptura en el segmento de falla ubicado entre los dos primeros subeventos determinados por Suárez & Nábêlek (1990). La ubicación del centro de intensidades del subevento de Caracas de 1812 sugiere que esto fue debido a que este segmento ya había sido roto en 1812. Por otra parte, es probable que al igual que en 1967, la ruptura haya sido discontinua y que la ruptura haya incluido un segmento de la falla al oeste del epicentro de 1967.

REFERENCIAS

1. Altez, R. (1998). Cronometrización extemporánea: Los sismos del 26 de Marzo de 1812 en Caracas y en Mérida. Revista Geográfica Venezolana. 39, pp. 297 – 326. [ Links ]

2. Altez, R. (2005). Todo lo que se movió en 1812 en la placa del Caribe: sismos, volcanes y transmisión de energía. Revista Geográfica Venezolana. Número especial, pp. 143-170. [ Links ]

3. Altez, R. (2006). El desastre de 1812 en Venezuela: sismos, vulnerabilidades y una patria no tan boba. Universidad Católica Andrés Bello y Fundación Polar, Caracas, Venezuela. 522 pp. [ Links ]

4. Altez, Rogelio & Laffaille, Jaime. (2006). La microzonificación sismo-histórica como complemento fundamental de la evaluación de la amenaza sísmica. Revista de la Facultad de Ingeniería. UCV, 21, pp.117-127. [ Links ]

5. Audemard, F. (2002). Ruptura de los grandes sismos históricos de los siglos XIX y XX revelados por la sismicidad instrumental contemporánea. Memorias XI Congreso Venezolano de Geofísica, 17-19 de noviembre. Caracas, Venezuela. [ Links ]

6. Bakun, W. H. & Wentworth, C. M. (1997). Estimating earthquake locations and magnitudes from seismic intensity data. Bull. Seismol. Soc. Am. 87, pp. 1502- 1521. [ Links ]

7. Bakun, W. H. & Wentworth, C. M. (1999). Erratum to estimating earthquake locations and magnitudes from seismic intensity data. Bull. Seismol. Soc. Am. 89 p. 557. [ Links ]

8. Castillo, W. (1983). El reforzamiento estructural de edificaciones de tapia, un caso en los Andes venezolanos. Boletín IMME. 72-73, pp. 53-82. [ Links ]

9. Centeno GraÜ, M. (1900). Interesantísimo estudio. La Linterna Mágica, Caracas, N° 250, p. 1. [ Links ]

10. Choy, J. E., Palme, C., Guada, C., Morandi, M. Klarica, S. (2007). Incorporación de intervalos de intensidades y de corrección por rumbo de fallas al método de Bakun & Wentworth para la evaluación de sismos históricos. IX Congreso Geológico Venezolano, Caracas, Venezuela. [ Links ]

11. Engdahl, E., Vanderhilst, R., Buland, R. (1998). Global teleseismic earthquake relocation with improved travel times and procedures. Bull. Seis. Soc. Amer., 88, pp. 722-743. [ Links ]

12. Fiedler, G. (1961). Áreas afectadas por terremotos en Venezuela. Memoria del III Congreso Geológico Venezolano, 3,pp. 1791-1810. [ Links ]

13. Fiedler, G. (1968). Estudio sismológico de la región de Caracas con relación al terremoto del 29 de Julio de 1967 (terremoto cuatricentenario): reporte y evaluaciones, Boletín Técnico IMME, 6, pp. 23-24. [ Links ]

14. Fiedler, G. (1972). La liberación de energía sísmica en Venezuela, volúmenes sísmicos y mapas de Isosistas. Memoria del IV Congreso Geológico Venezolano, pp. 2441-2462. [ Links ]

15. Funvisis (1997). Evaluación preliminar del sismo de Cariaco del 9 de julio de 1997, estado Sucre, Venezuela. Versión revisada. INFORME, CARACAS, VENEZUELA, 122 PP. [ Links ]

16. Gasperini, P., Bernardini, F., Valensise, G., Boschi, E. (1999). Defining seismogenic sources from historical earthquakes felt reports, Bull. Seism. Soc. Am. 89, pp. 94-110. [ Links ]

17. Grases, J. (1990). Terremotos destructores del Caribe. 1502-1990. 1º ed., Orcyt-Unesco. Montevideo, Uruguay. 132 pp. [ Links ]

18. Grases, J., Altez, R., Lugo, M. (1999). Catálogo de Sismos Sentidos o Destructores Venezuela 1530/1998. Academia de Ciencias Física, Matemáticas y Naturales Facultad de Ingeniería, UCV, Caracas, Venezuela. 650 pp. [ Links ]

19. Grases, J. & Rodríguez, J.A. (2001). Estimaciones de sismos venezolanos a partir de mapas de isosistas. Memorias del 2do Seminario Iberoamericano de Ingeniería Sísmica. Madrid, España. [ Links ]

20. GrÜnthal, G. (Ed.) (1998). European Macroseismic Scale 1998, Cahiers du Centre Europèen de Gèodynamique et de Seismologie, 15, Conseil de lΈurope, Luxembourg, 99 pp. [ Links ]

21. Hinzen, K. & Oemisch, M. (2001). Locations and magnitude from seismic intensity data of recent and historic earthquake data in the northern Rhine area, central Europe. Bull. Seismol. Soc. Am. 91, pp. 40-56. [ Links ]

22. Laffaille, J. & Ferrer, C. (2003). El terremoto del jueves santo en Mérida: año 1812. Revista Geográfica Venezolana, 44, pp. 107-123. [ Links ]

23. Laffaille, J. & Ferrer, C. (2005). El terremoto de Mérida de 1812: escrudiñando entre las páginas de una novela inconclusa en busca de información acerca de una historia real. Revista Geográfica Venezolana, número especial, pp. 217-232. [ Links ]

24. Laffaille, J., Klarica, S., Alvarado, M., Guerrero, O., Laffaille, M. (2007). La reevaluación de los terremotos históricos de Venezuela como una propuesta multidisciplinaria. Boletín de Historia de las Geociencias en Venezuela. 98, pp. 53-60. [ Links ]

25. Palme, C., Morandi, M., Choy, J. (2005a). Determinación de una relación lineal entre intensidad, magnitud y distancia epicentral para el occidente de Venezuela. INTERCIENCIA, 30, pp. 195-204. [ Links ]

26. Palme, C., Morandi, M., Choy, J. (2005b). Re-evaluación de las intensidades de los grandes sismos históricos de la región de la cordillera de Mérida utilizando el método de Bakun & Wentworth. Revista Geográfica Venezolana, número especial, pp. 233-253. [ Links ]

27. Rocabado, V., Sánchez, J., Schmitz, M., De Marco, R., Alvarado, L., Andrade, L. (2005). Caracterización Geofísica de la terraza de Barquisimeto y de la cuenca de tracción de Cabudare. Memorias del IV Coloquio Sobre Microzonificación Sísmica, Serie Técnica de FUNVISIS Nº 1, pp 224-229. [ Links ]

28. Richter, C. (1958). Elementary Seismology. Freeman, San Francisco. 768 pp. [ Links ]

29. Schmitz, M., Hernández, J.J., Audemard, F.A., Malavé, G., Andrade, L. (2005). Proyecto de Microzonificación sísmica en las ciudades de Caracas y de Barquisimeto. Memorias del IV Coloquio Sobre Microzonificación Sísmica, Serie Técnica de FUNVISIS Nº 1, pp 260- 263. [ Links ]

30. Suárez, G. & NÁbÊlek, J. (1990). The 1967 Caracas Earthquake: Fault Geometry, Direction of Rupture Propagation and Seismotectonic Implications. J. Geophys. Res. 95, pp. 17.459-17.474. [ Links ]

31. Wells, D. L. & Coppersmith, K.J. (1994). New Empirical Relationships among Magnitude, Rupture Length, RuptureWidth, Rupture Area, and Surface Displacement, Bull. Seismol. Soc. Am. 84, pp. 974–1002. [ Links ]

32. Yamazaki, Y., Audemard, F., Altez, R., Hernández, J., Orihuela, N., Safina, S., Schmitz, M., Tanaka, I., Kagawa, H. (2005). Estimation of the seismic intensity in Caracas during 1812 earthquake using seismic microzonation methodology. Revista Geográfica Venezolana, Número Especial, pp. 199-216. [ Links ]

REFERENCIAS ELECTRÓNICAS

Centro Regional de Sismología para América del Sur Ceresis. HTTP://WWW.CERESIS.ORG. ÚLTIMA CONSULTA: MAYO 2011.

GrÜnthal, G. (1998). EUROPEAN MACROSEISMIC SCALE EMS-98. HTTP://WWW.GFZ-POTSDAM.DE/PORTAL/-?$PART=CMSPART&DOCID=2051128. ÚLTIMA CONSULTA: MAYO 2011.

Guada, C. (2009). PORTAL DEL PROFESOR. HTTP://WEBDELPROFESOR.ULA.VE/CIENCIAS/CARLOSG/HTML/1812/TABLES/ANEXOSTABLE_3B.HTML, ÚLTIMA CONSULTA: MAYO 2011.

Inqua Environmental Seismic Intensity Scale 2007. ESI 2007. HTTP://WWW.ISPRAMBIENTE.GOV.IT/.../INQUA/ESI_2007_INQUA_SCALE_SPAGNOLO_2010_1.PDF. ÚLTIMA CONSULTA: MAYO 2011.

Sistema de Teleinformación de la Sismología Histórica de Venezuela, STSHV (2005). HTTP://WWW.SISMICIDAD.CIENS.ULA.VE. ÚLTIMA CONSULTA: MAYO 2011.

United States Geological Survey (USGS). HTTP://EARTHQUAKE.USGS.GOV/REGIONAL/QFAULTS/. ÚLTIMA CONSULTA: NOVIEMBRE 2008.