Observaciones del sismo del 16 de abril de 2016 de magnitud mw 7.8. Intensidades y aceleraciones.
INTRODUCCIÓN
El sismo registrado el sábado 16 de abril a las 18h58 (tiempo local), de magnitud 7.8 (Mw magnitud momento), cuyo hipocentro se ubicó frente a Pedernales (Manabí), a 20 km de profundidad, fue resultado del desplazamiento entre dos placas tectónicas: la placa de Nazca (placa oceánica) que se sumerge bajo la Sudamericana (placa continental). A este proceso se le conoce como subducción, y es el mismo fenómeno que originó los sismos del 31 de enero 1906 (Mw 8.8), que es el más grande registrado en Ecuador y el sexto más grande a escala mundial; el del 14 de mayo 1942 (Mw 7.8); 19 de enero de 1958 (Mw 7.8) y del 12 de diciembre de 1979 (Mw 8.1).
En este reporte se resumen las observaciones de daños en las edificaciones y de los datos recuperados de las estaciones de aceleración.
RESUMEN DE EVALUACIÓN DE INTENSIDADES
La intensidad es un indicador de la fuerza del evento, medido en personas, objetos y edificaciones. Para determinar los niveles de intensidad se utilizó la Escala Macrosísmica Europea (EMS98), que es una actualización de la Escala de Mercalli. Esta escala cuenta con 12 grados.
El resumen que se presenta en esta sección y el mapa de intensidades es una actualización del informe especial N°12, publicado en la página web del IG-EPN el 17 de abril (http://www.igepn.edu.ec/servicios/noticias/1316-informe-sismico-especial-n-12-2016). Como se menciona en dicho informe, los trabajos realizados en el campo permitieron tener una mejor definición de las zonas con mayores daños en las estructuras.
Durante los días posteriores al evento principal, se organizó la campaña de campo para el reconocimiento de daños en las edificaciones con el propósito de definir las áreas más afectadas en términos de la intensidad sísmica. Los grupos de trabajo estuvieron conformados por técnicos del Instituto Geofísico, Docentes de la Facultad de Ingeniería Civil y Ambiental de la Escuela Politécnica Nacional y estudiantes de Ingeniería Civil. En total, se establecieron 6 equipos que recorrieron las zonas que reportaron daños durante el evento principal. Es así que la información levantada en campo sirvió para determinar una intensidad máxima de 9 EMS en las zonas de Pedernales y San José de Chamanga. Este valor refleja el colapso observado en muchas edificaciones vulnerables y los daños en edificaciones construidas bajo las recomendaciones de las normas técnicas.
En ciudades y poblaciones costeras como Bahía de Caráquez, Jama y Canoa y las zonas centrales de Manta y Portoviejo, los daños en las edificaciones muestran una intensidad de 8 EMS. En estas dos últimas ciudades la distribución de daños es muy heterogénea y concentrada en ciertas áreas de la urbe.
En la figura 1, conocida como mapa de isosistas, se observa que los daños están distribuidos en la provincia de Manabí. En las provincias de Esmeraldas, Santa Elena, Guayas, Los Ríos y parte de Santo Domingo de los Tsáchilas, se observaron fisuras en las paredes de las edificaciones, que no comprometen a la estabilidad de la estructura. En general, en las provincias de la Sierra, la intensidad máxima es 4 EMS que indica que el sismo fue sentido ampliamente por la población.
En un informe posterior se presentará con mayor detalle el proceso de evaluación de intensidades.
DATOS DE ACELERACIÓN
La aceleración es un parámetro muy importante en el diseño de edificaciones porque es directamente proporcional a la fuerza que debe soportar una estructura durante un sismo. Esta información es utilizada en los códigos de construcción para generar las recomendaciones que rigen el diseño.
Hasta el momento se ha recuperado la información de varias estaciones de aceleración ubicadas en todo el país, que se detallan con círculos en el mapa de la figura 2. Estas estaciones son parte de la Red Nacional de Acelerógrafos (RENAC) del Instituto Geofísico (http://www.igepn.edu.ec/red-nacional-de-acelerografos). Además, se dispone de los valores de la máxima amplitud (PGA) de 7 sitios de la red de la empresa Oleoducto de Crudos Pesados (OCP). Se cuenta con información de las estaciones cercanas a la ruptura; sin embargo, se está trabajando en la recuperación de los datos de aquellas estaciones que no cuentan con transmisión en tiempo real.
Los acelerogramas de la figura 3 están pre-procesados y presentados en ventanas de 250 segundos; esto garantiza la visualización de las ondas principales del sismo. En algunos casos en que el tiempo GPS no está sincronizado, la selección de la ventana se realizó manualmente. En cada señal, se removió el promedio y la tendencia lineal. También se aplicó una función de encerado en los bordes de cada registro (5% cosine taper en inglés). En la tabla 1, los valores de la amplitud máxima obtenidos con este procedimiento se presentan para cada componente.
Tabla 1: Valores de la máxima amplitud (m/s2) para cada componente para todos los sitios con un registro disponible.
La Figura 2 muestra la distribución espacial de los valores de la máxima amplitud (PGA) de las tres componentes para los registros disponibles. Las mayores aceleraciones están registradas en Pedernales y en el sur del epicentro (estrella negra). Los valores en Manta, Portoviejo, y Chone son más grandes que las observaciones en Esmeraldas que está más cerca del epicentro y de la falla.
Figura 2: Distribución espacial de las estaciones de la Red Nacional de Acelerógrafos y de la red de OCP cuyos valores de aceleraciones sísmicas se muestran en la tabla 1. Los valores presentados corresponden a los máximos de las tres componentes (PGA: Peak Ground Accelerations) en m/s2.
La Figura 2 muestra la aceleración máxima obtenida en los sitios de las estaciones para la componente con la aceleración máxima. En Pedernales (estación APED) se registró un valor de 13.803 m/s2 (1.407 g). En la figura 3 se observa una diferencia importante entre las formas de onda de las estaciones al sur (AMNT, ACHN, APO2) y las estaciones ubicadas en el norte (LGCB, AES2). En el norte, los valores de aceleración pico son más bajos pero con una duración mayor, mientras que en el sur, se observa mayores aceleraciones pico pero con una duración más corta.
Figura 3: Acelerogramas del terremoto de las 18h58 (TL) del 16 de abril del 2016 ordenados con respecto a la distancia epicentral. Se han considerado las componentes con la máxima aceleración (PGA), la cual está indicada a la derecha de la señal. El tiempo de origen (0 s) corresponde al momento de ocurrencia del evento menos 20 s. Se utiliza en todos los casos la misma escala vertical.
En las figuras 4 y 5 se presentan los registros de dos estaciones en la provincia de Manabí. Las señales con los espectros de respuesta de las componentes horizontales están comparadas con los espectros de la Norma Ecuatoriana de Construcción (NEC15). Para las estaciones de Manta y Portoviejo, los espectros de respuesta son ligeramente menores a los de la NEC15, para un tipo de suelo D. Se ha seleccionado este tipo de suelo en base a la aplicación de la metodología de Zhao et al, 2006.
Figura 4: Señales de aceleración y espectros de respuesta para la estación Portoviejo.
Figura 5: Señales de aceleración y espectros de respuesta para la estación Manta.
COMPARACIÓN CON ECUACIONES DE PREDICCIÓN DE MOVIMIENTO DEL TERRENO (GMPEs)
Las ecuaciones de predicción de movimiento del terreno (GMPE por sus siglas en inglés) son funciones que nos permiten determinar valores de aceleración en función de parámetros sismológicos como magnitudes y ubicaciones del sismo.
En la Figura 6 se observan las aceleraciones máximas (PGA) de las señales previamente presentadas y los valores de PGA disponibles por cortesía de la empresa Oleoducto de Crudos Pesados (OCP). Las aceleraciones máximas [PGA] están representadas en función de la distancia a la falla, comparadas con la ecuación de predicción de movimiento (GMPE) de Abrahamson et al. (2015). Éste es el modelo más reciente para la interface de subducción que incorpora los datos de 43 terremotos de todo el mundo. Estudios anteriores han demostrado la buena concordancia de esta GMPE con los datos registrados de otros eventos en el país, pero con magnitudes más bajas (ver por ejemplo Beauval et al., 2012, comunicación personal de Stéphane Drouet, coordinador del grupo GMPEs del proyecto SARA América del Sur). La distancia a la falla se calculó tomando en cuenta el modelo de fuente definido por el programa de procesamiento SWIFT (Kumagai, comunicación personal). Este tipo de modelo arroja una distribución de probabilidades para obtener un valor de aceleración cuando se dispone de una magnitud, una distancia y un tipo de suelo.
La Figura 6 presenta la media geométrica de las dos componentes horizontales en función de la distancia mínima a la falla (RRUP). Se observa un buen ajuste entre los datos registrados y la ecuación mencionada. Los datos que tienen los valores más altos pueden explicarse con un efecto del suelo particular.
Figura 6: Comparación entre las aceleraciones máximas [PGA] con la ecuación de predicción de movimiento [GMPE] de Abrahamson et al. (2015), definida en el contexto de un sismo de interface en función de la distancia a la falla (RRUP, en km). Las aceleraciones máximas corresponden a la media geométrica de las componentes horizontales expresadas en g (aceleración de la gravedad). La GMPE está definida para un sismo de MW 7.8, una profundidad de 17 km, un VS30 de 760 m/s y un evento en la interface entre las dos placas. La línea roja continua representa el promedio de este modelo y las líneas punteadas son el promedio con ± 1 o 2σ (desviación estándar). (Figura preparada con la colaboración de Céline Beauval y Judith Marinière).
Observaciones finales
- Este sismo es el evento más grande en términos de aceleración registrado por la Red Nacional de Acelerógrafos (RENAC).
- La intensidad máxima evaluada es 9 EMS en zonas específicas de la provincia de Manabí (Pedernales y Chamanga) en donde se evidenció daños en la mayoría de las edificaciones. En zonas muy delimitadas de Portoviejo y Manta se alcanzó una intensidad de 8 EMS y en Guayaquil se estimó una intensidad de 6 EMS. Estos valores están relacionados con el tipo de suelo de las ciudades. En general, los daños observados se concentran en la provincia de Manabí.
- La forma de las líneas isosistas indica una atenuación muy marcada en la dirección oeste-este, si se compara con la dirección norte-sur.
- La aceleración máxima del suelo registrada, sobrepasa el valor de la aceleración de la gravedad en la estación APED de Pedernales para la componente este-oeste. En las estaciones de Manta y Portoviejo, estos valores sobrepasan el 30% de la aceleración de la gravedad.
AGRADECIMIENTOS
El trabajo de evaluación de intensidades se realizó con la colaboración de los profesores y estudiantes de la Facultad de Ingeniería Civil y Ambiental tanto de pregrado como de postgrado. Los autores agradecen a los coordinadores de los grupos: Ing. Juan Gabriel Barros, Dr. Matthieu Perrault, Ing. Diego Sosa, Ing. Diego Arévalo, Ing. Francisco Ortiz y por su intermedio a las personas que formaron parte de cada equipo de trabajo.
JCS, AL, CV, MR
Área de Sismología – Instituto Geofísico
Escuela Politécnica Nacional
