Análisis de probabilidad espacial de formación o reactivación de centros de emisión en los alrededores del volcán Michinmahuida, provincia de Palena, región de Los Lagos, Chile.

Abstract

El volcán Michinmahuida corresponde a un estratovolcán ubicado en la comuna de Chaitén, provincia de Palena, región de Los Lagos, Chile. Entre las localidades cercanas, la ciudad de Chaitén es la más poblada, y se encuentra a unos 20 km al oeste del volcán. Este volcán se caracteriza por una marcada elongación en dirección NE-SW, una estructura caldérica cubierta por un casquete glacial con elevación máxima de ~2.450 m s.n.m. y su emplazamiento sobre la traza principal del Sistema de Fallas Liquiñe-Ofqui (SFLO). En sus alrededores es posible identificar varios centros de emisión menores y depósitos que dan cuenta de su actividad pasada.

El Michinmahuida ocupa el décimo noveno lugar dentro del ranking de riesgo específico de volcanes activos de Chile, lo que le ha situado en la mira de la Red Nacional de Vigilancia Volcánica (RNVV) de SERNAGEOMIN, cuyo equipo de profesionales busca desarrollar su mapa de peligros. La evaluación de peligros volcánicos y el desarrollo de mapas de peligros es una tarea compleja, pues requiere del análisis conjunto de varios factores. En los últimos años, se han propuesto e implementado metodologías para elaborar mapas que permitan evaluar la probabilidad espacial de apertura de nuevos centros eruptivos, denominada susceptibilidad volcánica. En general, estos métodos utilizan diversas variables volcano-estructurales (que puedan tener algún grado de relación o control sobre el volcanismo) para llevar a cabo un análisis multiparamétrico. Este trabajo se ha desarrollado utilizando MatHaz para Matlab (Bertin y otros, 2019), que es un código de tipo abierto, diseñado para realizar evaluaciones espacio-temporales de peligro volcánico a partir de múltiples conjuntos de datos. Luego de modificar el código para trabajar en el área alrededor del volcán Michinmahuida, se creó una base de datos volcano-estructurales para la zona, donde se recopilaron todos los datos disponibles. De ella se seleccionaron y agruparon datos en seis datasets (centros de emisión menores, epicentros de sismos y estructuras N-S, NE-SE, NW-SE y E-W) para ser utilizados en el análisis, que consiste en el cálculo de una función de densidad de probabilidad (FDP) para cada dataset. Estas FDPs se obtienen mediante una estimación de densidad de kernel (EDK). Para integrar las seis funciones, se realiza una combinación lineal, obteniendo una FDP ponderada, donde cada set de datos contribuye a la ponderación según la importancia previamente asignada por un panel de expertos, obtenida en este mismo estudio mediante una encuesta.

Este estudio tiene como resultado un mapa de probabilidades que muestra las zonas con probabilidades altas, moderadas y bajas de apertura de centros de emisión, en base a los datos disponibles para el área. Este mapa puede ser utilizado para evaluar la susceptibilidad volcánica de la zona, aunque se deben tener en consideración algunas limitaciones metodológicas y prácticas que se enfrentaron durante su desarrollo.