I.    Antecedentes
El 19 de diciembre de 2016, a las 02h11 (TL), se registró un sismo de 5.3 Mw de magnitud (5.8 Ml, recalculada IG-EPN*), a una profundidad de 6.4 km frente a las costas del cantón Atacames (referencia 1). El sismo se produjo en una falla cuyo movimiento es inverso tal como lo muestra el mecanismo de la ruptura en la figura 1 (referencia 2).

Este evento junto a los ocurridos en esta zona desde inicios de diciembre del año anterior, probablemente están asociados a una falla tectónica ubicada en la placa continental, que posiblemente fue activada como resultado de la variación de esfuerzos generados por el sismo de subducción del 16 de abril de 2016 (ver figura 2).

* La nueva localización del IGEPN, obtenida a partir de los datos locales recuperados más tarde en el campo, ubica el epicentro frente a las costas del cantón Atacames con una latitud de 0.9274 N y una longitud de -79.8431 W.

La hipótesis más plausible, considerando los mecanismos focales inversos ubicados en la figura 3 junto con la línea sísmica de gran ángulo procesada por Agudelo et al. (2009, referencia 4 y figura 4), es que una falla tipo splay (SF, es decir una estructura ubicada sobre la corteza continental) es la fuente de la sismicidad registrada.

Según el informe preliminar de evaluación de intensidades (referencia 5), la parroquia Tonsupa perteneciente al cantón Atacames, provincia de Esmeraldas fue la más afectada presentando un valor de 7 EMS, este valor indica que es posible evidenciar daño en elementos estructurales de pocas edificaciones, así como daños superficiales (a nivel de mampostería) en muchas estructuras. Por otro lado, el tipo de suelo juega un papel importante en la respuesta de las edificaciones ante cargas sísmicas. A este fenómeno se le conoce como “Efecto de sitio” y para entender el comportamiento de los diferentes suelos encontrados en Tonsupa y Atacames, se utilizó el método de la razón espectral mHVSR (Nakamura, 1989), según como se describe en los siguientes apartados.

Tomando en cuenta la información antes citada, el IG-EPN decidió realizar un trabajo de inspección en campo para verificar los daños ocurridos en las estructuras por el sismo del 19/12/2016 en las ciudades de Atacames y Tonsupa; este reporte busca exponer los resultados obtenidos en dicha inspección.

II.    Descripción del trabajo realizado
La intensidad sísmica es un índice no instrumental que refleja la fuerza del movimiento del suelo mediante el efecto que este produce en las personas, en los objetos y en las construcciones, en una determinada posición durante un sismo. Este índice se evalúa mediante la aplicación de la Escala Macrosísmica Europea (EMS-98), la cual tiene 12 grados de intensidad e intenta crear un vínculo entre la sismología y la ingeniería (referencia 6), ampliando el número de tipos de edificación, definiendo mejor el grado de daño y sustentándose estadísticamente (ver figura 5).

Para el caso de la ciudad de Atacames, se conformó un equipo de trabajo que visitó la zona afectada durante los días 26, 27 y 28 de diciembre de 2016. En Tonsupa, por ejemplo, se pudo observar varios niveles de daño según el sector analizado. Una primera zona identificada consiste en viviendas residenciales de hasta dos pisos con autoconstrucciones mixtas o de hormigón armado, una segunda zona presenta edificaciones de entre 2 a 7 pisos de hormigón armado sin diseño sismoresistente o con diseño sismoresistente moderado, y una tercera zona con edificaciones de hasta 30 pisos con diseño sismoresistente alto.

La zona I presenta paredes desvinculadas de la estructura, mampostería agrietada y ciertas paredes exteriores caídas (ver figura 6).

La zona II que concentra los mayores daños, presenta pocas edificaciones colapsadas por patologías estructurales como piso blando o con mala configuración estructural, algunas edificaciones con daño estructural substancial a severo es decir falla a compresión del concreto, torcedura de varillas de refuerzo, grandes grietas en paredes, y muchas edificaciones con daño leve a moderado presentando grietas delgadas en enlucidos y grietas en paredes interiores (ver figura 7).

En la zona III, las edificaciones con diseño sismoresistente alto presentan un grado de daño 1, es decir sin daño estructural y daño no estructural leve (ver figura 8).

La ciudad de Atacames presenta tan solo una edificación con daño estructural severo, y pocas edificaciones con daño leve a moderado (ver figura 9).

Se decidió en campo además, realizar unas medidas de ruido ambiental en la parroquia de Tonsupa, buscando responder de manera alternativa los daños encontrados en las edificaciones dado un sismo de baja magnitud como el suscitado el 19/12/16.

Las medidas de ruido ambiental permiten obtener la frecuencia fundamental de resonancia del suelo (o las frecuencias de resonancia), con la presencia de fuertes contrastes de impedancia. El análisis se hace a través de las relaciones espectrales de las componentes horizontales del sensor sobre la componente vertical, mHVSR (referencia 7).

Es importante especificar que esta metodología permite obtener solamente las frecuencias de resonancia y no la amplitud de las funciones de transferencia de los sitios. La amplitud obtenida con mHVSR es solamente relativa y permite dar una primera idea de los sitios con mayor amplificación (referencia 8).

En total se logró tomar 7 medidas de ruido del ambiente en Tonsupa, valorando sobre todo en los lugares que presentaron edificaciones con mayores daños (ver figuras 10 y 11).

III.    Interpretación de los registros de aceleración
Al analizar los registros de las estaciones acelerográficas ubicadas en la zona epicentral, se observa que las estaciones de Tonsupa (ATON) y Atacames (AATC) registraron las aceleraciones más altas (ver figura 12 y tabla 1).

En la figura 13, se puede apreciar 11 de los acelerogramas con el registro del sismo del 19 de diciembre a diferentes distancias (5 a 175km). Se observa los diferentes tiempos de arribo de las ondas en función de la distancia. Normalmente, la amplitud baja en función de la distancia, excepto en contexto donde existan efectos de sitios que puedan generar amplitudes y duraciones más largas.

Se nota especialmente en el caso de los registros de la estación ASAM que presenta amplitud y duración más pequeñas que las estaciones AMA1 y AES1 que están más lejos del epicentro (ver figuras 13 y 14).

En el contexto de las estaciones de la RENAC, es posible estudiar los efectos de sitio a partir de eventos con la metodología de la razón espectral entre las componentes horizontales y la componente vertical [eHVSR], que permite como en el caso de la metodología Mhvsr únicamente encontrar la frecuencia fundamental de resonancia, en presencia de una diferencia importante de tipo de suelo (e.g., Lermo and Chavez-Garcia, 1993, referencia 10).

Todos los eventos que están reportados en el catálogo del IGEPN para los cuales tenemos un registro y que presentan un cociente entre la señal del evento y del nivel de ruido superior a 3, al menos para una frecuencia han sido usados. Además, para determinar la respuesta lineal se usa solamente los eventos que producen una amplitud máxima (PGA), inferior a 0.1 g (muchos estudios mantienen la hipótesis que este nivel de PGA puede separar los efectos lineales de los efectos no lineales en el suelo).

La figura 15 presenta para 6 estaciones, las curvas eHVSR que corresponden a la respuesta lineal obtenida con varios eventos y  las curvas eHVSR del evento del 19/12/2016.

Primero, con las curvas eHVSR lineales (en negro), es posible observar claramente que la estación ATON es la que presenta las amplitudes más altas sobre un rango de frecuencia amplio entre 1 y 5 Hz aproximadamente. La estación AATC presenta también una amplitud alta (pero menos que ATON) en el rango de frecuencia de 2 a 7 Hz. La estación ASAM presenta una respuesta muy plana que está de acuerdo con su ubicación sobre un suelo tipo roca. Las estaciones AMA1 y AES1 no presentan amplitudes muy altas, pero una ligera amplificación con un máximo alrededor de 2 Hz. La estación AES2 presenta un pico claro a 1 Hz. Además, solamente en el caso de la estación ATON, se nota una diferencia de las curvas eHVSR lineales con la respuesta eHVSR del evento del 19/12/2016. En el caso de ATON, se registró aceleraciones altas de 0.9g. Con este evento se observa picos más marcados, uno primero a 0.3 Hz; después el pico a 1.9 Hz se observa claramente a una frecuencia más baja alrededor de 1.3 Hz. Esta diferencia puede ser característica de un efecto no lineal pero no se observa en el campo evidencias de estos efectos.

IV.    Interpretación de las intensidades
De acuerdo a los daños observados en las diferentes estructuras, la máxima intensidad evaluada es 7 EMS en Tonsupa, lo que corresponde a que algunas edificaciones presentan daño en los elementos estructurales, llegando a fallas por compresión en el hormigón de las columnas e inclinación de las mismas (ver figura 16).

En Atacames, se evaluó una intensidad de 6 EMS, ya que pocas edificaciones sufrieron un daño estructural leve (ver figura 16).  

V.    Interpretación de las medidas de ruido ambiental
Los resultados observados en los 7 puntos medidos de ruido ambiental, indican que el subsuelo de Tonsupa es heterogéneo (ver figura 17).

De hecho, las curvas de mHVSR obtenidas presentan una variabilidad importante entre ellas. Hay una diferencia clara entre el noreste de la parroquia, con curvas mHVSR que no presentan picos muy claros, una ligera amplificación alrededor de 3 Hz (PR7, PR3, PR1, PR2 y PR6); y el suroeste (PR5, PR4, PR6) para lo cual las curvas presentan picos de amplificación claros: PR5 y PR4 presentan un primer pico claro alrededor de 1 Hz, y PR4 un segundo pico claro entre 4 y 6 Hz. En el caso de PR6, se observa un pico de amplitud alta para un rango de frecuencia largo entre 1.5 y 6 Hz.

Nota: Los resultados de mHVSR son comparables a los resultados de eHVSR con eventos de aceleraciones pequeñas. Como se observa a la estación ATON, una aceleración alta en un sitio puede generar una respuesta del suelo diferente de lo que se observa en un comportamiento lineal.

VI.    Conclusiones y Recomendaciones

•    El evento del 19 de diciembre de 2016 de 5.3 Mw de magnitud, generó un pico de aceleración máxima de 8,74 m/s2 registrado en la estación ATON, produciendo así una intensidad de 7 EMS en la localidad de Tonsupa y de 6 EMS en Atacames. Este valor de aceleración es alto, considerando la magnitud del sismo y podría deberse a un efecto de suelo muy importante.
•    Se puede ver con el análisis de las curvas mHVSR y eHVSR, que la zona sur occidental de Tonsupa (malecón) puede presentar efectos de sitio ante cargas cíclicas como las ocurridas en un sismo.
•    Los daños observados durante esta campaña de evaluación no corresponden a un evento de baja magnitud como el registrado el 19/12/16; esta conclusión puede ser explicada debido a varios factores como: efectos de sitio debido al tipo de suelo de la zona, también a que ciertas edificaciones se encontraban ya con daños debidos al sismo de abril de 2016 y también por problemas constructivos.
•    Se recomienda realizar estudios geotécnicos y geofísicos en detalle sobre el subsuelo de Atacames y Tonsupa, buscando así mejorar el ordenamiento territorial de las zonas propensas a sufrir un evento sísmico en el futuro, sobre todo por la presencia de estos valores de aceleración tan altos y por la variabilidad de frecuencias de vibración.
•    Se recomienda mejorar la evaluación de edificaciones post-desastre.

VII.    Bibliografía

(1)    http://www.igepn.edu.ec/portal/ultimo-sismo/informe-ultimo-sismo.html
(2)    http://www.igepn.edu.ec/informes-sismicos/sismicos-especiales/sism-e-2016/16697-informesi-smicoespecialn-39/file
(3)    http://www.igepn.edu.ec/informes-sismicos/sismicos-especiales/sism-e-2016/16705-informesi-smicoespecialn-40/file
(4)    Agudelo, W., Ribodetti, A., Collot, J.Y. & Operto, S., (2009). Joint inversion of multichannel seismic reflection and wide-angle seismic data; improved imaging and refined velocity model of the crustal structure of the north Ecuador-south Colombia convergent margin, J. geophys. Res., 114(B2). http://dx.doi.org/10.1029/2008JB005690
(5)    http://www.igepn.edu.ec/informes-sismicos/sismicos-especiales/sism-e-2016/16706-informesi-smicoespecialn-41/file
(6)    Grünthal, G. (Edit.). 1998. European Macroseismic Scale 1998, EMS-98. Conseil de L'Europe Cahiers du Centre Européen de Géodynamique et de Séismologie, Vol. 15. Luxembourg. 99p.
(7)    Nakamura Y. (1989). A method for dynamic characteristics estimations of subsurface using microtremors on the ground surface. QR of RTRI. 1;30:25–33.
(8)    Laurendeau, A., Barros, J-G., Alvarado A., y otros, Estimación de las frecuencias de resonancia de los suelos en varias ciudades de la provincia de Manabí con medidas de ruido del ambiente, Noviembre 2016.
(9)    Lermo, J. and F.J. Chavez-Garcia, (1993). Site effect evaluation using spectral ratios with only one station, Bull. Seism. Soc. Am. 83, 1574-1594.

JGBL, AL, CV, JCSA, DA, AA
Instituto geofísico
Escuela Politécnica Nacional