nueva página del texto (beta)
Español (pdf)
Artículo en XML
Referencias del artículo
Traducción automática
Enviar artículo por email
Citado por SciELO 
Similares en
SciELO 
Permalink
1. Introducción
El 24 de agosto de 2020 a las 21:51:08 (hora UTC), se presentó un evento sísmico, registrado por la red sísmica del Observatorio Vulcanológico y Sismológico de Costa Rica de la Universidad Nacional (OVSICORI-UNA), 13 km al sur de la ciudad de Jacó, Puntarenas, Pacífico Central de Costa Rica, con magnitud de 6.3 Mwp (Tsuboi et al., 1995), magnitud Richter Ml 6.1 y a una profundidad de 18 km. Con base en los reportes de la población en la página oficial de Facebook del OVSICORI-UNA, se determinó que el sismo fue sentido desde muy fuerte a fuerte en gran parte del territorio nacional. Fue reportado como percibido en forma oscilatoria y mecida, las ondas sísmicas fueron registradas por los instrumentos de baja ganancia instalados en distintos puntos del país y que mostraron la aceleración máxima del suelo en las componentes horizontales; no sobrepasó 15% del valor de la gravedad.
En este trabajo se llevó a cabo la relocalización del sismo principal, así como el cálculo de la inversión de parámetros de la fuente sísmica, como el centroide, ubicación y tiempo de origen, mecanismo focal, plano de falla (principal y auxiliar), distribución de las réplicas, distribución de deslizamientos en la falla, cambio de esfuerzos estáticos de Coulomb y el análisis del contexto sismotectónico donde se localizó el sismo. Además, con base en las aceleraciones pico (PGA) registradas en estaciones sísmicas, se determinaron las intensidades instrumentales y su correspondencia en la escala de Intensidad Mercalli Modificada (IMM), lo que permitió generar un mapa de intensidades para todo el país, el valor máximo registrado en la escala IMM fue de VI. Esto se comparó con los reportes de la población en la red social Facebook del OVSICORI-UNA, en su mayoría el sismo fue reportado desde muy fuerte a fuerte, además se reportó la caída de objetos en áreas pobladas cercanas al epicentro (13 km), lo que coincide con una intensidad máxima IMM de VI.
Aunque la idea principal en este trabajo fue analizar el sismo de Herradura del 2020 y teniendo que en el 2017 se localizó un sismo, con magnitud Mw 6.4 (llamado sismo de Esterillos en este trabajo) y muy cercano al sismo de Herradura es que también procedemos a la caracterización de esa fuente sísmica.
2. Contexto tectónico
La placa de Cocos se subduce perpendicular a la trinchera mesoamericana (TMA) originando un régimen de contracción en la placa superior (Mescua et al., 2017). Actualmente las variaciones morfotectónicas a lo largo del antearco de Costa Rica han sido atribuidas principalmente a cambios composicionales y morfológicos presentes en la placa de Cocos. Estos cambios han llevado a muchos autores a proponer una segmentación de la placa de Cocos, la cual ha sido definida por anomalías magnéticas (Barckhausen et al., 2001), batimetría multibeam (von Huene et al.,1995; von Huene et al., 2000), estudios sísmicos (Walther y Flueh 2002; Sallares y Charvis 2003; Walther, 2003; Wilson et al., 2003; Sallares et al., 2005) y muestreo geológico (Hauff et al., 1997; Werner et al., 1999; Wilson et al., 2003; Hoernle y Hauff, 2007).
Como resultado, en esta área se han definido tres segmentos de litosfera oceánica con edades variables, denominados CNS-1, CNS-2 y CNS-3, originados en el centro de expansión Cocos-Nazca, que limitan al norte con la Dorsal del Pacífico Este y al sur con la Placa Nazca.
El segmento CNS-1, que se subduce bajo el centro y sur de la Península de Nicoya, se caracteriza por una batimetría rugosa asociada a un fragmento de litosfera creada durante el Mioceno inferior (21-22,5 ma) (Barckhausen et al., 2001). El segmento CNS-2 comprende el dominio de los montes submarinos que chocan con la costa del centro y sur de Costa Rica e incluye la litosfera engrosada y rugosa del punto caliente de las Galápagos que contiene la Dorsal de Cocos (Gardner et al., 1992; Fisher et al., 1998). Este segmento exhibe una batimetría rugosa suave con numerosos montes submarinos cónicos y otras características batimétricas tales como el Plateau de Quepos en el centro (von Huene et al., 1995) y un área de corteza oceánica engrosada con un fondo marino elevado y un espesor de sedimentos variable de hasta 1 km, asociado a la cresta de Cocos en el sector sur (Walther et al., 2000; von Huene et al., 2000; Kolarsky et al., 1995).
En cuanto a la sismicidad asociada al segmento CNS-2, Protti et al. (1995a) correlacionan el acoplamiento sísmico con las características del fondo marino subducido y observan que históricamente los montes submarinos aislados pueden reducir el acoplamiento entre las placas y generar sismos con una magnitud máxima de Mw 7.0. Superficialmente, la sismicidad es frecuente y de menor magnitud, lo cual sugiere que los altos batimétricos que ingresan a lo largo del margen central de Costa Rica forman asperezas aisladas que solo pueden generar sismos de tamaño moderado (Mw 6.0-7.3) (Protti et al., 1995b; Bilek et al., 2003; Proyecto Global Centroid Moment Tensor (GCMT), 2020).
Cada uno de estos montes submarinos subducidos deja una depresión a su paso, que puede ser rastreada hasta 55 km hacia el frente de deformación (Sak et al., 2009). Ese rastro es creado por la deformación y el levantamiento del talud continental cuando el monte submarino pasa por debajo de él, seguido por el colapso cuando se retira el soporte del área levantada (Hühnerbach et al., 2005). Las áreas de levantamiento sobre los montes submarinos se caracterizan por patrones complejos de fallas normales y de deslizamiento (Sak et al., 2009; Hühnerbach et al., 2005).
Estos montes submarinos pueden actuar como asperezas o barreras temporales en la zona de acople entre las placas (Watts et al., 2010), lo que dificulta la propagación de ruptura o genera esfuerzos locales que contribuyen a la ocurrencia de sismos fuertes en la zona. Husen et al., (2002) asocian el sismo del 25 de marzo de 1990, Mw 7.3 (GCMT, 2020) dentro del segmento CNS-2 en la entrada del Golfo de Nicoya, ubicado a ~25 km al oeste del sismo del 24 de agosto de 2020, con un monte submarino que podría estar actuando como aspereza.
Entre los sismos más importantes ocurridos recientemente en la zona, se encuentran el del 13 de noviembre de 2017 a las 02:28:20.3 hora UTC, Mw 6.4, que se ubicó a 14 km al sureste del epicentro del sismo de Herradura o el sismo de Cóbano en marzo de 1990, -13:22:55 hora UTC, Mw 7.3, que también se ubicó cercano a este.
3. Metodología
Se realizó una caracterización de los parámetros de la fuente sísmica, para ello se utilizó principalmente el catálogo sísmico y formas de onda de la red sísmica del OVSICORI-UNA. Además, se revisaron los catálogos disponibles en línea del International Seismological Center (ISC) (2020), del Proyecto Global Centroide Momento Tensor (GCMT) (2020), del USGS (2020) y Guendel, (1986). Los sismos se localizaron con el programa GenLoc de Fan et al., (2006); para la inversión, el centroide y tensor de momento (CTM), se utilizó el paquete ISOLA (Sokos y Zahradnik, 2008 y 2013) y FOCMEC (Snoke et al., 1984), esto permitió obtener el mecanismo focal por medio de las polaridades de los primeros arribos de la onda P. Para la distribución del deslizamiento en el plano de falla se implementó el paquete computacional de inversión lineal de desplazamiento (Gallovic y Zahradnik 2011; Gallovic et al., 2014). Para la localización de los sismos, cálculo del CMT, distribución de deslizamiento en la falla y mecanismo focal se usó el modelo de velocidad de Quintero y Kissling (2001) y en el caso del cálculo de cambio de esfuerzos estático de Coulomb (ACFS) se empleó el paquete computacional Coulomb 3.3 (Toda et al., 2005; Lin et al., 2004). Se usó el método de cambio de esfuerzos con el fin de encontrar áreas donde se ve un aumento de ΔCFS y señalar la posible ocurrencia de futuros eventos sísmicos.
Para el catálogo sísmico del OVSICORI-UNA, los sismos entre 1986-1992, fueron localizados con el programa HYPOINVERSE (Klein, 1984). A partir de 1992 hasta 03 de marzo del 2010 se procesó por medio del programa HYPOCENTER (Lienert y Haskov, 1995). Siendo el 04 de marzo del 2010, cuando las asociaciones de sismos y localizaciones fueron hechas con ANTELOPE (BRTT, 2010), el cual hace uso del programa de localización llamado LOCSAT (Bratt y Bache, 1988), basado en el modelo de velocidad IASPEI91 (Kennett y Engdahl, 1991) para sismos regionales y telesismos, y las librerías GENLOC (Generalized Earhquake-Location) de Fan et al. (2006), basado en el modelo de velocidad de Quintero y Kissling (2001) para la sismicidad local. El uso de un modelo 1D para la localización de sismos y la metodología de inversión puede producir que sismos superficiales aparezcan en el aire o en la capa de agua; estos sismos son un poco más profundos; pero en este trabajo son graficados en las localizaciones obtenidas sin hacerles modificación.
En el catálogo OVSICORI-UNA, las magnitudes de los sismos que ocurrieron antes del 10 de marzo del 2010 se obtuvieron de la siguiente forma: la magnitud de los eventos sísmicos fue obtenida del promedio de las magnitudes parciales calculadas para cada estación, según la duración de las trazas en las estaciones de corto período. Estas magnitudes parciales se calcularon según la relación obtenida por el Servicio Geológico de los Estados Unidos (USGS) para sismos en Alaska:
Mc = -1.16+2.01*log10 T+0.0035*Δ+ ΔM
T: duración del sismo
Δ: distancia epicentral
ΔM: corrección para cada estación
A partir de 1992 se implementó un sistema digital en el centro de registro, que además contaba con estaciones de banda ancha, con respuestas instrumentales calibradas.
Estas permitían simular la respuesta de un instrumento Wood-Anderson para el cálculo de la magnitud local Richter. En este cálculo se utiliza la siguiente relación propuesta por Hutton y Boore (1987):
ML = log10(A) + 1.11*log10 (Δ) + 0.00189* Δ - 2.09
Δ: Distancia epicentral.
A: Amplitud pico-pico de la onda.
Después de implementado el sistema de adquisición sismológico ANTELOPE, se dejó de cuantificar la magnitud coda y, en vez de esta y para sismos locales, se cuantifica la magnitud local (Ml) que representa la magnitud Richter, donde el paquete computacional (BRTT, 2010) aplica la metodología y funciones de corrección propuestas por Richter (1935). Es a partir de la implementación de ANTELOPE que se comienza con el cálculo de la magnitud momento de la onda P: Mwp (Tsuboi et al., 1995). Las magnitudes son calculadas automáticamente por el paquete ANTELOPE y se reporta la mediana de los valores de magnitud cuantificados en cada estación. En los años anteriores, se calculaba el valor medio de las magnitudes obtenidas en cada estación. Las magnitudes calculadas según sea la metodología van a variar; para el caso del sismo del 11 de noviembre del 2017, USGS (2020) y GCMT (2020) lo reportan con magnitudes entre 6.1 y 6.6 y el OVSICORI-UNA con Ml 6.4 y en este estudio se encontró Mw 6.4; y la diferencia con otros estudios que hacen cálculo de magnitud momento global es por el ancho de banda usado en la inversión; siendo usadas frecuencias más altas en el caso de usar datos locales. Para el sismo del 24 de agosto de 2020, el GCMT (2020) y USGS (2020) reportaron magnitudes entre 5.5 y 6.1 y el OVSICORI-UNA con Ml 6.1, Mwp 6.3 y en este trabajo se encuentra magnitud momento usando datos locales de Mw 5.9.
El mapa de intensidades instrumentales, para el sismo de Herradura, en la escala de Mercalli Modificada (MMI) se generó mediante una interpolación estadística (Krigging) a través del software Arc Gis 10.7, para ello se utilizaron las aceleraciones pico (PGA) registradas en 21 estaciones distribuidas en el territorio nacional.
La conversión de PGA a MMI se llevó a cabo aplicando la conversión propuesta por Wald et al. (2005) (Tabla 1).
Tabla 1 Intensidades instrumentales en la escala de Mercalli Modificada (MMI).
| Intensidad instrumental (IMM) |
Aceleración (%g) |
Velocidad (cm/s) |
Sacudida percibida |
Potencial de daño |
|---|---|---|---|---|
| I | <0.17 | <0.1 | No sentida | Ninguno |
| II-III | 0.17-1.4 | 0.1-1.1 | Débil | Ninguno |
| IV | 1.4-3.9 | 1.1-3.4 | Ligera | Ninguno |
| V | 3.9-9.2 | 3.4-8.1 | Moderada | Muy ligero |
| VI | 9.2-18 | 8.1-16 | Fuerte | Ligero |
| VII | 18-34 | 16-31 | Muy fuerte | Moderado |
| VIII | 34-65 | 31-60 | Severa | Moderado a fuerte |
| XI | 65-124 | 60-116 | Violenta | Fuerte |
| X + | >124 | >116 | Extrema | Muy fuerte |
Los resultados obtenidos fueron comparados con los reportes de la población en la página oficial del Facebook del OVSICORI-UNA para determinar la relación entre las intensidades instrumentales en la escala MMI y las características de cómo fue sentido el sismo en diferentes ciudades y poblados del país.
4. Resultados
4.1. RÉPLICAS DEL SISMO DE HERRADURA
El sismo de Herradura de las 21:51:08.1, hora UTC, se localizó a una latitud 9.501 y longitud -84.665, a una profundidad de 18 km, con magnitud Ml 6.1 . Las réplicas de este evento (para dos días posteriores) fueron localizadas principalmente frente a la costa del Pacífico Central de Costa Rica, en un área de 400 km2 aproximadamente (Figura 1). En total se registraron 91 réplicas cuyas magnitudes se encontraron entre Ml -0.8 y Ml 3.6. Referente a las profundidades, estas oscilaron entre 3.5 y 32.7 km del total de esas réplicas, 80 se ubicaron a profundidades mayores a 15 km, concentradas en una banda de 20 a 25 km de profundidad y con epicentros desde la costa hacia la TMA.
Figura 1 Sismicidad asociada al sismo de Herradura, el sismo principal se indica por una estrella amarilla y las réplicas por círculos, donde el color indica la profundidad, los anaranjados con profundidades entre 0 y 20 km y los rojos presentan profundidades mayores a 20 km. Las estaciones sísmicas usadas en la localización del sismo principal se indican por triángulos, el color indica el residuo; verde si es menor a 0.5 segundos, amarillo entre 0.5 y 1.5 segundos y rojo para residuos mayores a 1.5 segundos. La estación sísmica más cercana al epicentro del sismo principal está ubicada en Herradura, Jacó. Para el cálculo de las réplicas también se usaron estas estaciones sísmicas.
4.1.1. MECANISMO FOCAL DEL SISMO DE HERRADURA
La solución del mecanismo focal (MF) del evento principal, calculado usando polaridades de los primeros arribos y el programa FOCMEC (Snoke et al., 1984) muestran un plano de falla inclinado al NE con un ángulo de 31°, el cual corresponde a una falla inversa con azimut de 287° y ángulo de deslizamiento de 71° (287°/31°/71°).
La solución del MF obtenida de la inversión de formas de onda y el paquete ISOLA (Sokos y Zahradnik, 2008 y 2013), también da como resultado un mecanismo inverso, con un plano de falla buzando al NE, con ángulo de inclinación de 20°, rumbo de 311° y deslizamiento de 98° (311°/20°/98°) (Figura 2B). El sismo está relacionado con el acople entre la placa de Cocos y la placa Caribe (Co-Ca) en el Pacífico Central de Costa Rica.
Figura 2 A) Correlación de las formas de onda de desplazamiento reales (negro) y sintéticas (rojo) usando un ancho de banda de 0.05 a 0.12 Hz. El gráfico para cada componente de la estación sísmica indica tiempo (segundos) vs desplazamiento (metros), B) Mecanismo focal calculado usando inversión de formas de onda (línea continua en la esfera focal) y polaridades de primeros arribos, indicando los planos de falla con líneas discontinuas. Las estaciones usadas en la inversión de las formas de onda muestran la respectiva polaridad; con rectángulos blancos se indica la dilatación (D) y con rectángulos negros la compresión (C), el eje de Presión (P) está indicado por un círculo negro y el de tensión (T) por un círculo blanco. El eje de presión (215/25) tiene dirección noreste y el de tensión con parámetros 27/65, C) Correlación vs profundidad para el sismo de Herradura. El porcentaje doble par de fuerzas (DP) se indica por una barra de color a la derecha y en este caso es de 81% DP.
La correlación entre los datos sintéticos y observados usando estaciones en el campo cercano entre 18 y 129 km del epicentro es alta (Figura 2A) y el mecanismo focal obtenido está en concordancia con los mecanismos calculados por el USGS (2020). Asimismo, el centroide calculado usando datos de la red del OVSICORI-UNA, se ubicó debajo del epicentro a una profundidad de 25 km y +1.3 segundos después del tiempo de origen y una magnitud Mw 5.9.
Además, la correlación por profundidad también es bastante alta, alcanzando un valor de 0.83 (Figura 2C) para 25 km de profundidad; la cual corresponde a una profundidad similar a los sismos ubicados entre la interfaz CO-CA para esta zona.
Se hizo búsqueda en el plano, pero la mejor correlación se obtuvo para un punto ubicado por debajo del epicentro; lo que indica que los deslizamientos en la falla fueron focalizados y no se produjo rompimiento a lo largo del plano de falla, como lo evidencian las réplicas, que se mantuvieron concentradas y cercanas al hipocentro.
4.1.2. INTENSIDADES INSTRUMENTALES EN LA ESCALA DE MERCALLI MODIFICADA (MMI)
Las máximas PGA se registraron en la estación sísmica JACO, ubicada en Herradura, Garabito, Puntarenas, estas alcanzaron 11.52%g, que corresponde a una intensidad de VI en la escala IMM. Sin embargo, sismos con estas características tienen la particularidad de no producir daños importantes en infraestructuras, aunque la sacudida es fuerte y sentida por la población.
Para establecer la relación planteada por Wald et al., (2005) se utilizaron datos de PGA de 21 estaciones sísmicas, distribuidas por todo el país, y se determinó que las intensidades se corresponden con valores entre II y VI (Tabla 2, Figura 3B).
Tabla 2 PGA en estaciones del OVSICORI-UNA.
| Estación | E (%g) | N(%g) | Z (%g) | Intensidad |
|---|---|---|---|---|
| Valor máximo | Valor máximo | Valor máximo | I-XII | |
| BATAN | 0.74 | 0.63 | 0.18 | III |
| CCOL | 0.80 | 0.56 | 0.04 | III |
| CDITO | 0.18 | 0.24 | 0.14 | III |
| COVE | 0.47 | 0.56 | 0.21 | III |
| CDM | 1.03 | 0.99 | 0.57 | III |
| CPMI | 0.78 | 0.55 | 0.33 | III |
| CTCR | 0.59 | 0.62 | 0.22 | III |
| DUNO | 1.19 | 1.78 | 0.50 | IV |
| HDC3 | 0.34 | 0.32 | 0.29 | III |
| JACO | 11.52 | 11.03 | 3.76 | VI |
| JTS | 1.09 | 0.65 | 0.30 | III |
| OCM | 0.85 | 0.75 | 0.31 | III |
| ORTG | 0.11 | 0.14 | 0.11 | II |
| PEZE | 0.93 | 1.37 | 0.31 | III |
| PNPB | 0.39 | 0.33 | 0.13 | III |
| POTG | 0.58 | 0.60 | 0.16 | III |
| RIFO | 0.53 | 0.40 | 0.21 | III |
| RIMA | 1.29 | 0.99 | 0.35 | III |
| RIOS | 1.97 | 2.10 | 0.77 | IV |
| TIGR | 0.86 | 0.70 | 0.32 | III |
| VMAR | 0.33 | 0.01 | 0.25 | III |
Figura 3 A) Intensidades instrumentales en la escala de Mercalli Modificada (MMI) registradas para el sismo de Herradura del 24 de agosto de 2020, Ml 6,1. Para determinar las intensidades IMM se utilizó la relación de Wald et al., 2005; para conversión de PGA a intensidad. Los puntos blancos indican ciudades y poblados, B) Estaciones sísmicas, señaladas con triángulos rojos, donde se registraron las intensidades instrumentales.
Las mayores intensidades (VI) se concentraron en la ciudad de Jacó y en Herradura; en Cóbano y Paquera, en el extremo sureste de la Península de Nicoya, se registró una intensidad máxima de V, las intensidades menores a III se registraron en poblados como San Isidro del General, Limón, Upala y en ciudades de la Gran Área metropolitana (GAM) destacando centros urbanos como San José (capital del país), Alajuela, Cartago y Heredia (Figura 3A).
Con base en la información recabada en el Facebook del OVSICORI-UNA se identificaron algunas particularidades relevantes, ya que la población reportó el sismo sentido desde muy fuerte a fuerte. Se destaca que en la ciudad de Jacó, ubicada a 13 km del epicentro, se reportó la caída de objetos en supermercados. Con la información facilitada por los usuarios de la red social (Facebook oficial del OVSICORI-UNA) se lograron sistematizar las principales características del sismo en las siete provincias del país (Tabla 3). Se hace énfasis en que la información recabada correspondió a reportes en el Facebook del OVSICORI-UNA y no a encuestas. En total se registraron 2679 reportes, los cuales se clasificaron por provincia y según la descripción macrosísmica.
Tabla 3 Descripción macrosísmica del sismo de Herradura y su distribución por provincias.
| Provincia | Ciudades o poblados | Descripción macrosísmica |
|---|---|---|
| Alajuela | Ciudad de Alajuela | Sumamente fuerte y largo. Movimiento de objetos y caída de ellos en algunas viviendas. |
| Turrucares | Demasiado fuerte y largo. La población reportó que percibieron dos sismos largos de abajo hacia arriba y luego hacia los lados. |
|
| Poás, Naranjo, La Garita, Ciudad Quesada, Sarchí, Río Segundo de Alajuela, Grecia, San Ramón, Palmares y La Fortuna de San Carlos. |
Muy fuerte y de larga duración. | |
| Guácima de Alajuela | Se sintió hacia arriba y abajo y después hacia los lados. Larga duración. |
|
| Cartago | Ciudad de Cartago, Orosí, Caballo Blanco, Tres Ríos de La Unión y Villas de Ayarco. |
Muy fuerte y largo, en la ciudad de Cartago se reportó movimiento de objetos. |
| Juan Viñas | Muy fuerte y duró aproximadamente 13 segundos. |
|
| Heredia | Ciudad de Heredia, Santo Domingo, San Rafael, San Pedro de Barva y San Pablo. |
Muy fuerte y de larga duración, en San Rafael se reportó movimiento de objetos. |
| Varablanca | Muy leve pero de larga duración. | |
| Mercedes Norte y Barva | Se reportó como dos sismos sentidos muy fuertes y de larga duración. |
|
| Guanacaste | Liberia, Sardinal, Nosara | Leve pero de larga duración. |
| Tilarán, Nandayure, Cañas, Guayabo y Santa Cruz. |
Fuerte y de larga duración. | |
| Limón | Ciudad de Limón, Siquirres, Cocles y Pueblo Nuevo de Guácimo. |
Fuerte y de larga duración. |
| Río Jiménez de Guácimo | Muy fuerte. | |
| Guápiles y Pensurth | Leve, en Guápiles provocó el movimiento de objetos. |
|
| Puntarenas | Ciudad de Puntarenas, Barranca | Muy fuerte |
| Monteverde | Fuerte y de larga duración, movimiento de objetos. |
|
| Parrita | Muy fuerte y de larga duración, se reportaron objetos caídos y quebrados. |
|
| San Vito | Leve. | |
| Quepos y Esparza | Fuerte y de larga duración. | |
| Río Claro de Golfito, | Leve y de larga duración. | |
| Orotina, ciudad de Jacó, Herradura y Uvita | Muy fuerte y de larga duración. | |
| San José | Ciudad de San José, Pavas, Moravia Guadalupe, Santa Ana, San Antonio de Coronado, Alajuelita, Tibás, San Pablo León Cortés, San Francisco de Dos Ríos, Sabanilla Montes de Oca, Desamparados, Pérez Zeledón y Acosta. |
Fuerte y de larga duración. En Moravia reportaron sentir dos sismos, una más fuerte que el otro. En Sabanilla Montes de Oca se reportó la caída de objetos. |
| Vueltas del Jorco, Aserrí | Muy fuerte y de larga duración. |
5. Análisis de la sismicidad en el área de estudio
En el área de estudio se han registrado otros sismos importantes, tales como a) el sismo del 25 de marzo de 1990 a las 13:22:55 hora UTC, Mw 7.3 (GCMT, 2020), b) el sismo del 20 de agosto 1999 a las 10:02:20 hora UTC, Mw 6.9 (GCMT, 2020) y c) el sismo del 13 de noviembre 2017, Mw 6.4 (este trabajo) y se presentó a las 02:28:20 hora UTC (Figura 4). Algunos de los sismos registrados en el área presentan dos particularidades importantes: a) microsismicidad abundante y b) que en un período de tiempo corto los eventos pre-monitores y sus réplicas, alcanzan una magnitud significativa, incluso en ocasiones cercanas a la del evento principal. Por ejemplo, el sismo ocurrido el 21 diciembre de 1939 (20:54:53, hora UTC) con magnitud Ms 7.3, y al día siguiente (04:44:03 hora UTC) se presentó una réplica magnitud Ms 6.8; el 19 de julio 1956 (23:26:30 y 23:38:04 hora UTC) se presentaron dos sismos de igual magnitud (Ms 6.2). El 23 agosto de 1978 (00:38:31 y 00:50:27 hora UTC), se ubicaron dos eventos en Sámara, frente a la Península de Nicoya, con magnitud de Ms 7.0 y 6.9, respectivamente; el 25 de marzo de 1990 (13:16:05 hora UTC), ocurrió un evento pre-monitor de Ms 6.4 y un evento principal de Mw 7.3 a las 13:22:55 hora UTC; el 20 de agosto de 1999 (10:02:20 hora UTC), un sismo de magnitud Mw 6.9 y una réplica al día siguiente (10:49:50 hora UTC) de Mw 5.8.
Figura 4 Se muestran con círculos rojos los epicentros para sismos que se han registrado en el Pacífico Central de Costa Rica, desde 1984 hasta agosto de 2020, con magnitudes Ml > 5.9. Se destacan cuatro principales eventos con su respectiva fecha y magnitud asociada a) el sismo del 25 de marzo de 1990 (Mw 7.3), b) el sismo del 05 de septiembre de 2012 (Mw 7.6), conocido como el terremoto de Nicoya, c) el sismo del 13 de noviembre de 2017 (Ml 6.4) y d) el sismo del 01de junio de 2019 (Ml 4.3), e) el sismo del 25 de agosto de 2020 (Ml 6,1). Los MF para el sismo del 05 setiembre 2012 (Quintero et al., 2014) y para el sismo de 25 de marzo 1990 (GCMT, 2020); los otros tres se realizaron en este trabajo. Los MF indican que los sismos ocurren como consecuencia de fallas inversas asociadas a la subducción de la placa de Cocos bajo la placa Caribe, buzando al nor-noreste y en bajo ángulo.
El penúltimo evento fuerte registrado en esta área corresponde al sismo de Esterillos del 13 de noviembre de 2017 (02:28:20 hora UTC), con magnitud Mw 6.4 y con una réplica de Ml 5.2 a las 02:32:31 hora UTC (Güendel, 1986; ISC, 2020; OVSICORI-UNA, 2020; Global Catalogue, 2020; para información sobre sismos históricos en la zona) (Tabla 4).
Tabla 4 Sismos históricos importantes entre 1939 y 2020.
| Sismo | Año | Mes-día | Hora UTC | Latitud | Longitud | Profundidad | Magnitud | Magnitud | Fuente |
|---|---|---|---|---|---|---|---|---|---|
| 1 | 1939 | 12-21 | 20:54:47.20 | 10.140 | -84.630 | 33.0 | 7.3 | 7.6MISC | 1, 2 |
| 2 | 1939 | 12-22 | 04:43:58.10 | 09.800 | -84.550 | 33.0 | 6.8 | 6.8SISC | 1, 2 |
| 3 | 1956 | 07-19 | 23:26:37.00 | 09.550 | -84.140 | 93.0 | 6.2 | 5.8SISC | 1, 2 |
| 4 | 1956 | 07-19 | 23:38:14.00 | 09.520 | -83.960 | 81.0 | 6.2 | 1, 2 | |
| 5 | 1978 | 08-23 | 00:38:27.00 | 09.760 | -85.570 | 26.0 | 7.0SUSGS | 1, 3 | |
| 6 | 1978 | 08-23 | 00:50:24.40 | 09.890 | -85.500 | 42.0 | 6.9SISC | 1, 2 | |
| 7 | 1990 | 03-25 | 13:16:05.30 | 09.585 | -84.948 | 27.9 | 6.4SUSGS | 4, 3 | |
| 8 | 1990 | 03-25 | 13:22:55.20 | 09.628 | -84.923 | 23.3 | 7.3WGCMT | 4, 5 | |
| 9 | 1999 | 08-10 | 14:55:42.80 | 09.263 | -84.120 | 35.0 | 5.7LUNA | 5.5WGCMT | 4, 5 |
| 10 | 1999 | 08-10 | 14:55:42.80 | 09.263 | -84.120 | 35.0 | 5.7LUNA | 5.5WGCMT | 4, 5 |
| 11 | 1999 | 08-20 | 10:02:20.00 | 09.038 | -84.156 | 14.0 | 6.2LUNA | 6.9WGCMT | 4, 5 |
| 12 | 1999 | 08-20 | 10:42:15.60 | 08.884 | -84.107 | 21.6 | 5.3LUNA | 5.2bUSGS | 4, 3 |
| 13 | 1999 | 08-21 | 10:49:50.00 | 08.976 | -84.014 | 45.3 | 5.8WGCMT | 4, 5 | |
| 14 | 2017 | 11-13 | 02:28:20.30 | 09.463 | -84.546 | 27.0 | 6.4LUNA | 4 | |
| 15 | 2017 | 11-13 | 02:32:30.90 | 09.454 | -84.546 | 17.3 | 5.2LUNA | 4 | |
| 16 | 2019 | 06-01 | 06:53:29.20 | 09.542 | -84.573 | 12.3 | 4.3LUNA | 4 | |
| 17 | 2020 | 08-24 | 21:51:08.00 | 09.501 | -84.665 | 18.3 | 6.1LUNA | 4 |
2 International Seismological Centre, (2020)
Referente al sismo de Esterillos, se obtuvo que los parámetros de la fuente sísmica lo ubican en la latitud 9.463 y longitud -84.546, con una profundidad de 27 km. La magnitud es de Ml 6.4, con el centroide a 21 km de profundidad debajo del epicentro y a +4.0 segundos después del tiempo de origen. El plano de falla se orienta en dirección 289°, y buza 17°, su deslizamiento es de 85° (Figura 5). La magnitud momento obtenida usando datos de la red local es dos décimas inferiores a la obtenida usando datos de la red global por el GCMT (2020), que la ubica en Mw 6.6, cuya diferencia en magnitud se debe al ancho de frecuencias usadas en la inversión. Por otra parte, se muestran (Figura 6A) los epicentros de algunos sismos registrados entre el 2010 y agosto de 2020, con su respectivo perfil indicado a lo largo de la línea 1, así como la sismicidad asociada a los sismos de Esterillos y Herradura. La sismicidad se concentró en los primeros 25 km de profundidad (Figura 6B), con epicentros al sur de la ciudad de Jacó. Un perfil de la sismicidad a lo largo de la línea 1 muestra que el sismo se encuentra en la parte superior de la zona sismogénica y en contacto CO-CA y el sismo principal corresponde a una falla inversa. El sismo de Esterillos se caracteriza por tener réplicas hacia arriba del epicentro y el de Herradura por estar alrededor de la zona epicentral.
Figura 5 Inversión por el centroide y tensor de momento del sismo del 13 de noviembre de 2017. A) Gráfico de formas de ondas observadas y sintéticas en estación sísmicas usadas, en color negro y rojo, respectivamente, B) Mecanismo focal y estaciones usadas en la inversión, C) Gráfico que muestra la correlación a diferentes profundidades, con diagrama de barras indicando el porcentaje de doble par de fuerzas.
Figura 6 Epicentros de sismos registrados por la red sísmica del OVSICORI-UNA, entre 2010 hasta septiembre del 2020 representados por círculos negros, verdes y rojos, los dos últimos indican epicentros para los primeros 10 días de réplicas después de ocurrido el sismo de Esterillos y Herradura, respectivamente, B) Proyección vertical de los sismos que se encuentran a una distancia máxima de 30 km (líneas grises gruesas en mapa), a ambos lados de la línea 1 y las proyecciones de sus MF en la línea 1, se muestra el eje P y T, con círculos blancos y negros, respectivamente. Con estrellas azules se indican las localizaciones de los sismos de Herradura y Esterillos. Se presenta una ampliación en la esquina inferior izquierda de los eventos principales.
Es posible que estos sismos estén asociados con la subducción de los altos batimétricos que ingresan a lo largo del margen central de Costa Rica y que forman asperezas aisladas que pueden generar terremotos que podrían alcanzar una magnitud Mw 7.3 (GCMT, 2020).
Después del sismo de Esterillos (Mw 6.4) se ubicó sismicidad cortical superficial, donde ocurrió un sismo magnitud Ml 4.3 el 01 de junio de 2019 a las 06:53:29 hora UTC, cuyo origen es atribuido a una falla normal con componente de desplazamiento y parámetros para un plano nodal azimut=18°, buzamiento= 61° y deslizamiento= -139°, parámetros obtenidos con inversión por el CMT usando datos sísmicos de la red del OVSI-CORI-UNA. En el pasado sismicidad como la del sismo de Cóbano del 25 de marzo 1990, Mw 7.3 disparó sismicidad en la placa continental; que ocurrió en superficie (Protti et al., 1995), después del sismo del 13 de noviembre de 2017 no es muy clara en una relación directa con este; pero existe sismicidad localizada en la placa continental, principalmente en fuentes sísmicas, cercanas al epicentro.
En profundidad y a unos 100 km de distancia de la trinchera, en la zona de contacto de las placas CO-CA se localizan pocos sismos y los sismos como el de Esterillos y el de Herradura ubican el centroide por debajo del epicentro, con réplicas hacia la TMA; por lo que la zona que evidencia poca sismicidad puede ser una barrera que impida rompimiento a lo largo del plano de inmersión de la falla (Figura 6). El cambio de esfuerzos estáticos de Coulomb para los sismos de Esterillos (Mw 6.4) y de Herradura (Mw 5.9) graficados con perfiles que cortan la TMA y en dirección al buzamiento de la placa de Cocos muestran que existe una acumulación de esfuerzos en los extremos superior e inferior de la falla inversa que originó estos sismos (Figura 7). El extremo inferior de esta zona de interfaz sismogénica entre las placas de Cocos y Caribe, ubicada al noroeste de donde ocurrieron los sismos en mención, ha mostrado poca sismicidad en el pasado y donde existe acumulación positiva de esfuerzos puede indicar una zona con potencial para generar futuros eventos sísmicos.
Figura 7 A) Cambios de esfuerzos estáticos de Coulomb debido al sismo del Esterillos (parámetros de falla 289°/17°/85°, Mo=4.35e+18 Nm, Mw 6.4, profundidad 21 km, epicentro indicado por estrella blanca, B) Cambios de esfuerzos estáticos de Coulomb debido al sismo de Herradura (parámetros de falla 11 °/20°/98°, Mo=9.719e+17, Mw 5.9, profundidad 25 km), epicentro indicado por estrella blanca, se asume un coeficiente de fricción de 0.4. La barra de colores indica los cambios de esfuerzos de Coulomb en bares y con círculos negros se muestran algunos de las localizaciones del catálogo sísmico del OVSICORI-UNA. Una línea de color rojo muestra la TMA y en negro y verde se indica la falla en profundidad y la proyección en superficie, un cuadrado en rojo indica la dimensión de falla.
Resultados preliminares de la distribución de deslizamiento final en el plano de falla 289°/17°/85° para el sismo de Esterillos, y para el sismo de Herradura en el plano de falla 311°/20°/98°, muestran una concentración de deslizamientos concéntricos alrededor del centroide con deslizamientos máximos de 0.6 y 0.1 metros, respectivamente y con una ligera concentración hacia la parte inferior de la falla para el sismo de Esterillos y hacia la superficie para el sismo de Herradura (Figura 8).
Figura 8 A) Modelo de deslizamiento obtenido a partir de la inversión de deslizamiento lineal de fallas finitas para el sismo de Esterillos, Mw 6.4, B) Modelo de deslizamiento obtenido a partir de la inversión de deslizamiento lineal de fallas finitas para el sismo de Herradura, Mw 5.9. La escala de la derecha representa la variabilidad de la amplitud de deslizamiento. El asterisco en color gris denota el centroide proyectado en el plano de falla (no se usa en la inversión). La distancia es a lo largo del azimut.
6. Discusión
El sismo del 24 de agosto del 2020, 21:51:08 hora UTC, conocido como el sismo de Herradura, ubicado en el Pacífico Central de Costa Rica, corresponde a una falla inversa, con hipocentro a 18 km de profundidad, su centroide se encuentra debajo del epicentro a 25 km de profundidad y +1.3 segundos después del tiempo de origen. Sismos con este tipo de mecanismo focal están asociados a la interfaz Cocos-Nazca, tal y como se observa en el perfil en profundidad (Figura 6). La ubicación del centroide y de las réplicas muestran que los hipocentros se encuentran concentrados en un área pequeña y que no existe rompimiento a lo largo de la inmersión de la placa de Cocos bajo la placa Caribe.
Lo mismo ocurrió con el sismo de Esterillos con centroide por debajo del epicentro.
Estos sismos inversos de mayor magnitud están asociados con la subducción de los altos batimétricos que ingresan a lo largo del margen central del antearco de Costa Rica y forman asperezas aisladas que tienen potencial para generar sismos con una magnitud de Mw 7.3 (GCMT, 2020). El modelo de deformación de la placa superior propuesto para esta región por Domínguez et al. (1998) indica que se desarrollan fallas inversas arqueadas delante del monte submarino, fracturas subverticales que se extienden desde arriba del monte submarino, algunas con movimiento de desplazamiento de rumbo.
El monte submarino acrecionado proporciona una aspereza local, pero podría transferir la tensión a las fracturas preexistentes en la parte superior, en relación directa con los sismos que abarcan un amplio rango de profundidades.
Estos sismos con mecanismo focal normal estarían asociados al colapso producido por la depresión creada durante la deformación del talud continental cuando el monte submarino pasa por debajo, similar a los modelos propuestos por Dominguez et al. (1998), von Huene et al. (2000), Husen et al. (2002), Hühnerbach et al. (2005) y Sak et al. (2009).
Esta región sobre el monte submarino se caracteriza por patrones complejos de fallas normales y deslizamientos (Sak et al., 2009; Hühnerbach et al., 2005), por lo que se esperarían fallas normales producidas por un monte submarino en subducción, debido al colapso del material de la placa superior (Dominguez et al., 1998).
7. Conclusiones
El sismo de Herradura, 21:51:08 hora UTC, ubicado en el Pacífico Central de Costa Rica, corresponde a una falla inversa, con hipocentro a 18 km de profundidad, su centroide se encuentra debajo del epicentro a 25 km de profundidad y +1,3 segundos después del tiempo de origen, con deslizamiento en forma circular, distribución uniforme y máximos deslizamientos de 0.1 m y magnitud Mw 5.9. Sismos con este tipo de mecanismo focal están asociados a la interfaz CO-CA y asociados con la subducción de la placa de Cocos bajo la placa Caribe. La disminución de la sismicidad que se observa debajo de la interfaz podría estar relacionado a un monte submarino cortado en la placa subducida.
La ubicación del centroide y de las réplicas del sismo de Herradura muestran que los hipocentros se encuentran concentrados en un área pequeña y que no existe rompimiento a lo largo de la inmersión de la placa de Cocos bajo la placa Caribe; lo mismo ocurrió con el sismo de Esterillos, a las 02:28:20.3 hora UTC, con magnitud Mw 6.4 con centroide por debajo del epicentro y máximos deslizamientos de 0.6 m y con concentración de réplicas hacia la superficie y en zonas con cambio de esfuerzos de Coulomb negativos.
La acumulación de esfuerzos estáticos como lo muestra el cambio de esfuerzos de Coulomb se da en las fronteras superior e inferior de la falla en la zona de interfaz sismogénica entre la placa de Cocos y la placa Caribe, zona que se caracteriza por haber subducido montes submarinos que han sido generadores de sismos fuertes en el pasado. La distribución de deslizamientos muestra concentración alrededor del centroide y sin deslizamiento hacia la parte inferior de la interfaz sismogénica; por lo que la distribución de esfuerzos y deslizamiento en la falla producto de estos sismos puede indicar que la parte inferior de la zona sismogénica en esta parte del Pacífico Central de Costa Rica tiene alto potencial sísmico. Es posible que un monte submarino en la zona esté actuando como barrera que impide la propagación de la ruptura a lo largo de la zona de acople CO-CA.
En la proyección de la sismicidad por profundidad se observa que los sismos en la zona de contacto de las cortezas CO-CA se encuentran concentradas en la zona más superficial quedando la zona inferior de la interfaz con capacidad de romper con otro sismo fuerte o por el contrario actuar como barrera que impida el rompimiento hacia debajo de la zona sismogénica.
La magnitud del sismo de Herradura fue fuerte (Ml 6.1; Mw 5.9), sin embargo, al no ser un sismo superficial y al romper en la TMA y concentrar sus réplicas en esta zona, no generó daños en ciudades y poblados cercanos al epicentro, tal como la ciudad de Jacó (ciudad más próxima). Las intensidades instrumentales en la escala de Mercalli Modificada (MMI) máximas de VI y aceleraciones pico de ~12%g en la estación más cercana al epicentro, provocó caída de objetos y pánico en las personas en el área de mayor aceleración del suelo, sin embargo, no se reportaron daños a estructuras o pérdidas de vidas.
