1. Introducción
La medición de las emisiones constantes de dióxido de azufre (SO2) es una piedra angular del monitoreo de los volcanes activos. Los volcanes se consideran una de las mayores fuentes naturales de SO2 a nivel mundial, emitiendo a la atmósfera 23 ± 2 millones de toneladas métricas anualmente (Carn et al., 2017). Por otro lado, las bajas concentraciones de fondo del SO2 en la atmósfera (~1 ppbv en el aire claro; Breeding et al., 1973) y sus firmas espectrales distintivas tanto en la región del ultravioleta (UV) como en la región del infrarrojo térmico (TIR, por sus siglas en inglés), permiten una detección mucho más eficiente al utilizar métodos de percepción remota, a diferencia de otros volátiles volcánicos más abundantes en la pluma como el agua (H2O) o el dióxido de carbono (CO2), debido a que estos elementos se encuentran en mayor concentración en la atmósfera (Edmonds, 2003; Oppenheimer et al., 2011). Las mediciones continuas de SO2 permiten comprender el almacenamiento y la liberación de volátiles a partir de magmas y, por lo tanto, los controles dinámicos y químicos sobre las erupciones volcánicas. Esto hace que el monitoreo del SO2 sea un componente clave en la mitigación del riesgo volcánico y el pronóstico potencial de una erupción (Edmonds, 2003; Shinohara, 2008; Galle et al., 2010; Jimenez-Escalona et al., 2011; Oppenheimer et al., 2011).
La magnitud física más utilizada para el monitoreo volcánico es generalmente el flujo de SO2 (expresado en T/día o kg/s), el cual, combinado con otros sistemas de monitoreo y datos (p. ej. sismicidad), proporciona una amplia visión de la dinámica del comportamiento interno de un volcán. Los cambios en la emisión del flujo de SO2 se usan frecuentemente como un precursor de una erupción (Jiménez-Escalona et al., 2011; Werner et al., 2012). La interpretación de las variaciones del flujo de SO2 puede ayudar a identificar procesos volcánicos importantes, en particular: (a) reabastecimiento del sistema magmático con material joven (Caltabiano et al., 1994), caracterizado por un aumento gradual del flujo de SO2 medido durante un período de tiempo relativamente largo; (b) agotamiento de los elementos volátiles de un cuerpo de magma, estos procesos están marcados por una disminución gradual del flujo de SO2 (Kazahaya et al., 2004) y, debido a que los gases constituyen la fuerza motriz de las erupciones volcánicas, generalmente este evento precede al final de un período eruptivo; y (c) la obstrucción del sistema de conducto superior por cristalización de magma (Fischer et al., 2002). Este proceso promueve la acumulación de gas bajo el tapón, hasta su destrucción o fallo cuando se supera su resistencia, lo que se cree que alimenta la actividad denominada “vulcaniana” y “estromboliana” (Iguchi et al., 2008).
Las mediciones del flujo de SO2 volcánico, utilizando sensores de percepción remota, inició a mediados de los años setenta con el uso del Espectrómetro de Correlación (COSPEC, por sus siglas en Inglés) (Stoiber et al., 1983). El COSPEC es un sensor que utiliza la radiación ultravioleta entre 300 y 330 nm dispersada por la atmósfera como fuente de radiación. Puede ser utilizado de forma fija para medir los gases de salida del cráter del volcán o se puede montar en un vehículo terrestre o aéreo para la realización de medición en transectos de la pluma. Stoiber et al. (1983) calcularon el error relativo de las mediciones de COSPEC entre valores de ±13% y ±42%, los cuales incluyen errores tanto instrumentales como metodológicos. En este estudio, se obtuvieron mediciones de COSPEC desde un vehículo terrestre el cual estima las concentraciones de SO2 tomando mediciones verticales en transectos que cruzan la pluma volcánica como se observa en la Figura 1. La concentración media de la sección tomada de la pluma es multiplicada por la velocidad de la pluma para determinar el flujo de SO2 de la emisión.
Por otro lado, a principios de los años ochenta se descubrió la capacidad de los sensores satelitales para detectar y estimar las concentraciones de SO2 volcánico emitidas a la atmósfera (Krueger, 1983); desde entonces, se han desarrollado diferentes metodologías tanto para sensores que detectan en el UV (p. ej. TOMS, OMI, GOME-2), como para sensores que permiten la detección en el TIR (p. ej. ASTER, IASI, MODIS).
Al tratarse de técnicas de percepción remota, tanto las que utilizan sensores midiendo desde superficie como los sensores satelitales, presentan valores de incertidumbre en las estimaciones que se realizan, dependiendo fuertemente de la distancia de la pluma al sensor y la interacción que presenta la longitud de onda utilizada para la estimación del SO2 con los elementos de la atmósfera. Por otro lado, un factor muy importante en la calidad de los datos que se obtienen con estas técnicas es la complejidad de la logística que se debe seguir para la adquisición de la información, así como el contar con el personal capacitado para el análisis de ésta. Por esta razón, lo recomendable no es que una técnica sustituya a la otra, sino se debe buscar la manera de complementar la información de un sensor con la que se obtiene con otros.
El objetivo de este trabajo es proponer una metodología para ajustar dos técnicas de percepción remota utilizadas para monitorear el flujo de SO2 en el volcán Popocatépetl, con el propósito de contar con la mayor cantidad de información para la construcción de series de tiempo con promedios diarios. Las mediciones en superficie con el sensor COSPEC y la estimación de SO2 por medio de imágenes satelitales del sensor MODIS (Moderate Resolution Imaging Spectroradiometer) son comparadas. Debido a que las mediciones del sensor COSPEC presenta las ventajas de llevar a cabo una medición a menor distancia de la pluma volcánica (~6 a 8 km) y es menos afectada por las condiciones ambientales presentes, se toma como referencia de calibración para corregir el sesgo que se presenta en las estimaciones derivadas de las imágenes MODIS. Por otro lado, las imágenes MODIS detectan el SO2 en la atmósfera en el rango del infrarrojo térmico permitiendo la obtención de información tanto en el día como en la noche. Este rango de detección permite contar con una resolución temporal de hasta 4 imágenes diarias de la zona del volcán Popocatépetl. Para este ajuste con los valores estimados con COSPEC, se debe tomar en cuenta la baja resolución espacial de las imágenes MODIS (1 km), así como la distancia que debe recorrer la señal que es utilizada para obtener la información (radiación infrarroja de la superficie). Se propone tomar en consideración las ventajas y desventajas de cada sensor con el propósito de buscar su complementariedad sinérgica que ayude a mejorar las series de tiempo que muestran el comportamiento del volcán.
2. Descripción de los sensores
2.1. Sensor COSPEC
El sensor COSPEC fue desarrollado en Canadá a principios de los sesenta con el propósito de monitorear los NOx y el SO2 de la industria (Horton et al., 2006) y posteriormente también se aplicó para medir SO2 en plumas volcánicas (Moffat y Millan, 1971; Stoiber y Jepsen, 1973; Stoiber et al., 1983; Delgado-Granados et al., 2001). Millan (2008) describe de manera clara los principios de COSPEC, que usualmente opera en modo pasivo en estudios de monitoreo volcánico. Este sensor es utilizado para estimar la concentración de SO2 en una pluma volcánica comparando la extinción de radiación electromagnética en el rango de longitud de onda de 300 ≤ λ ≤ 330 nm. COSPEC requiere de la observación de la radiación UV difusa en la atmósfera, la cual puede ser recolectada en un punto estacionario o una plataforma móvil, tanto aérea como por tierra. Cuando el instrumento es transportado por un vehículo, se hacen mediciones en transectos de la pluma perpendiculares al desplazamiento de esta, estimando la abundancia de SO2 en la columna en el campo de vista del instrumento. El sensor que se coloca en un vehículo terrestre apunta de forma vertical a la pluma de SO2, como se muestra en la Figura 1, para estimar la concentración de SO2 en la vertical. El éxito de esta metodología depende de la intensidad de radiación solar en la atmósfera y de la disponibilidad de carreteras o caminos alrededor del volcán que permitan la medición de los transectos de la pluma volcánica. Debido a que el instrumento sólo ve una pequeña porción de la pluma para la estimación del flujo diario de SO2, se toma el promedio de varios transectos medidos durante las horas del día mientras la luz del sol es óptima para el uso del instrumento. Además, para la estimación del flujo de SO2 es muy importante considerar la componente perpendicular del vector de viento a la altitud de la pluma de SO2 en cada uno de los transectos. Este dato podría obtenerse de estaciones meteorológicas o aeropuertos locales relativamente cercanos a la zona de monitoreo. Sin embargo, un método más preciso y confiable implica determinar las características del viento en la zona donde se mide el transecto, ya sea por observación o con apoyo de instrumentos de medición en el caso de ser posible (Williams-Jones et al., 2008).
2.2. Sensor MODIS
El sensor MODIS se encuentra a bordo de los satélites NASA Terra-I y Aqua-I los cuales recorren trayectorias en órbita polar a una altitud de 705 km sobre el nivel del mar (http://modis.gsfc.nasa.gov). El sensor tiene un ángulo de visión de ± 55° que permite cubrir zonas con un ancho de banda de 2330 km lo que permite realizar una cobertura global completa cada dos días. MODIS es un sensor multiespectral que adquiere información en 36 bandas de detección desde la región del visible hasta el infrarrojo térmico (0.405 μm a 14.385 μm) y su resolución espacial varía entre 250, 500 y 1000 m. Las bandas del sensor MODIS utilizadas para las observaciones de productos volcánicos son B28 (7.3 μm), B29 (8.5 μm), B31 (11.0 μm) y B32 (12.0 μm) con una resolución espacial de 1 km. Su amplia cobertura espectral permite que MODIS sea de los pocos sensores que en la actualidad pueden detectar tanto SO2 como ceniza volcánica (Watson et al., 2004; Matiella-Novak et al., 2008; Jiménez-Escalona et al., 2011).
El algoritmo de detección de SO2 en imágenes multiespectrales en el infrarrojo térmico fue descrito por primera vez por Realmuto et al. (1994) para el espectrómetro aerotransportado TIMS (Thermal Infrared Multispectral Scanner). Posteriormente, esta técnica ha sido utilizada por diversos autores para aplicarla a sensores satelitales como es el caso de MODIS (Watson et al., 2004). Las ventajas que presenta la detección en el infrarrojo térmico son, que debido a la naturaleza de la emisión de esta radiación electromagnética, la detección de productos volcánicos se puede hacer con imágenes que son tomadas tanto de día como de noche. Para el caso de la zona de estudio, los sensores MODIS embarcados en Terra-I y en Aqua-I permiten la detección de imágenes durante las horas de día (entre las 11:00 y las 15:00 hora local) y durante las horas de noche (entre las 23:00 y las 4:00 hora local), por lo que se puede tener un máximo de 4 imágenes diariamente. Por otro lado, la desventaja que presenta este tipo de algoritmos para imágenes en el infrarrojo térmico, es la fuerte interacción con los elementos contenidos en la atmósfera, principalmente con el vapor de agua en la zona de detección de 7.3 μm. Por esta razón, en la zona de la troposfera se recomienda el uso de la zona de detección de 8.6 μm (Realmuto y Worden, 2000; Watson et al., 2004).
3. Actividad del volcán Popocatépetl
El volcán Popocatépetl (19.02ºN, 98.62ºW, 5425 msnm) se encuentra a ~70 km al sureste de la Ciudad de México y ~40 km al oeste de la ciudad de Puebla, y se encuentra rodeado de cientos de medianos y pequeños poblados en un radio de 60 km alrededor del volcán (De la Cruz-Reyna y Siebe, 1997; Sheridan et al., 2001). Estudios geológicos muestran que el volcán Popocatépetl ha presentado grandes erupciones en los últimos 15000 años, los cuales son descritos en Siebe et al. (1996). Los eventos eruptivos de mayor intensidad en los últimos 800 años ocurrieron en 1510, 1519, 1540, 1580, 1664 y 2001, ocasionando caídas generalizadas de cenizas que afectaron tanto a las zonas pobladas como rurales (Martin-Del Pozzo et al., 2016). Por otro lado, estos estudios históricos han mostrado que la duración de los episodios eruptivos han presentado periodos de tiempo de entre menos de un año y más de 30 años, separados por periodos de reposo entre 7 y más de 100 años.
El último período eruptivo del volcán Popocatépetl ocurrió entre 1920 ‒ 1927 y después de este, tuvo un lapso de reposo de alrededor de 70 años (Delgado-Granados et al., 2001). El volcán presentó un aumento notorio de su actividad fumarólica a partir de 1990. Durante el año de 1993, se observó un aumento de la actividad interna del volcán, representado por el aumento de frecuencia y magnitud de micro-sismos, así como por las emisiones de gases a la atmósfera. Mediciones hechas por Goff et al. (1998) y Delgado-Granados et al. (2001) muestran que las tasas de emisión de SO2 aumentaron de ~2000 ton/día en 1994 a ~50000 ton/día en 1997 hasta llegar a valores de 100000 ton/día a principios de 2001. A partir de los episodios de enero de 2001, la actividad eruptiva del Popocatépetl comenzó a descender gradualmente hasta 2007. El Centro de Aviso de Ceniza Volcánica de Washington (VAAC, por sus siglas en inglés), reportó en 2001, 49 erupciones y en 2007 sólo se reportaron 8, número que se mantuvo hasta 2011 en donde se reportaron más de 14 erupciones explosivas que depositaron ceniza en la atmósfera.
4. Metodología
4.1. Estimación del flujo de SO2 con COSPEC
Debido al principio de detección del sensor COSPEC, es importante considerar las condiciones de luz solar para determinar el mejor momento para llevar a cabo las mediciones. Para el caso del volcán Popocatépetl, el mejor momento para realizar las mediciones con COSPEC se presenta entre las 10:00 y 14:00 hora local (Delgado-Granados, 2008). Es durante este tiempo que los transectos de la pluma volcánica son obtenidos al pasar al sensor por debajo de ésta. Para lograrlo, es necesario seguir varios caminos que rodean el volcán dado que la dirección de la pluma puede presentarse en varias direcciones a lo largo del año. Como se muestra en la Figura 2, el sistema de carreteras que rodean al volcán permite llevar acabo mediciones en cualquier dirección que presente la pluma.
La información de la medición del transecto de la pluma se registra, así como el tiempo de inicio y final de cada uno de ellos y la distancia de la medición al cráter del volcán. Un punto importante que se debe considerar para la obtención del flujo diario de SO2 a partir de la medición de transectos de una pluma volcánica, es que la emisión de gases en el cráter no mantiene un flujo continuo en el tiempo. Es muy común la presencia de intervalos de aumento y descenso de la cantidad de gases que son emitidos en función de los procesos de la actividad interna del sistema. Estos cambios de la concentración de gases pueden ser identificados a lo largo de la pluma en varias decenas de kilómetros. Por esta razón, para calcular la tasa de emisión promedio diaria, durante el tiempo que se llevan a cabo las mediciones, se realizan el mayor número de transectos posibles de la pluma volcánica. Debido a lo ancho de la pluma y la velocidad que se puede desplazar el vehículo, generalmente se llevan a cabo de 2 a 4 mediciones de transectos.
Para la obtención de los datos de viento, se toma en cuenta que la pluma de gases que emite el volcán Popocatépetl se transporta en altitudes arriba de los 5000 m. Al encontrarse en la atmósfera libre se considera que no hay una influencia importante de la capa límite planetaria y por consiguiente las características del viento no cambian significativamente a nivel de mesoescala. Por lo que, para las estimaciones en la zona del volcán Popocatépetl se utilizan los datos de los sondeos atmosféricos generados por el Servicio Meteorológico Nacional (SMN), en la estación de la Ciudad de México (19.40°N, 99.20°O), reforzando esta información con estimaciones visuales en campo (Delgado-Granados et al., 2001).
El cálculo del flujo de SO2 en cada transecto se obtiene estimando la concentración promedio de este, multiplicado por la componente perpendicular del viento y línea de muestreo. Este procedimiento se describe con la Ecuación 1 propuesta por Okita (1971).
Donde F es el flujo de SO2, α es el ángulo entre la dirección del viento y la línea de muestreo del transecto, C es la abundancia de SO2 en la columna, dS es un pequeño segmento de la línea de medición y th es el umbral de abundancia de la columna de SO2.
4.2. Estimaciones del flujo de SO2 con MODIS
Las imágenes de MODIS utilizadas en este trabajo fueron descargadas de las bases de datos de NASA (https://ladsweb.modaps.eosdis.nasa.gov/). Este sitio web permite al usuario seleccionar el tipo de datos de satélite, la ubicación de la imagen deseada basada en la latitud y longitud, y el tiempo de la imagen deseada. MODIS captura 288 imágenes cada día alrededor de la Tierra. Bajo condiciones ideales, sería posible obtener 4 imágenes al día, aunque en la mayoría de los casos, no se pueden obtener más de 2 imágenes. Sin embargo, las condiciones de baja o nula nubosidad en la zona de estudio generalmente no se satisfacen durante la estación de lluvia (junio a agosto).
Para recuperar la masa total de SO2 y el área de la pluma volcánica de las imágenes MODIS, se empleó el programa MAP_SO2 (Realmuto y Worden, 2000). Este software fue diseñado como una herramienta de análisis de datos de imágenes en el infrarrojo térmico, tanto de sensores aerotransportados como adquiridos desde satélites. El procedimiento de recuperación de la información de la imagen se basan en una versión personalizada del modelo de transmisión atmosférica MODTRAN® (MODerate resolution atmospheric TRANsmission) (Berk et al., 2004; Berk, 2013), el cual es una aplicación estándar reconocida para el modelado de transferencia radiativa.
El programa detecta la cantidad de radiancia que es absorbida por el SO2 presente en la atmósfera en una longitud de onda de ~8.6 μm (Realmuto et al., 1994; Realmuto y Worden, 2000). El éxito de esta técnica depende de una eficiente caracterización de: (i) la superficie debajo de la pluma volcánica, (ii) la atmósfera en la trayectoria de la superficie a sensor, y (iii) la geometría de la pluma (altura y espesor). El espectro de radiación térmica infrarroja observado por el sensor Ls (λ,Ts ), puede ser representado por el modelo simplificado mostrado en la ecuación 2 (Realmuto et al., 1994):
Donde, T s y ε g (λ) representan la temperatura y la emisividad superficial y B(λ,T s ) representa la Ecuación de Planck para la radiancia del cuerpo negro. El producto de ε g (λ) y B(λ,T s ) describe la radiancia superficial, que es atenuada por la transmisión atmosférica τ a (λ), en el camino recorrido al sensor. La presencia de una pluma volcánica se expresa principalmente a través de la absorción de esta radiación de la superficie. Sin embargo, las contribuciones del flujo de radiación dispersa que llega a la superficie L d (λ) y es reflejada, así como la que viaja hacia el sensor L u (λ) también pueden ser significativas en la estimación de la radiancia en el sensor L s (λ,T s ).
Para la estimación del flujo de SO2, se utilizó un derivado del método transversal según la metodología de medición de plumas volcánicas con COSPEC desde una plataforma móvil (Stoiber et al., 1983). Esta técnica se ha aplicado en imágenes de alta y moderada resolución espacial (p. ej. ASTER y MODIS), ya que permite estimar el flujo de SO2 de pequeñas plumas que luego pueden compararse fácilmente con datos obtenidos de mediciones en tierra (Theys et al., 2013). Los datos del viento en la zona de la pluma se obtuvieron de los datos de reanálisis preparados por el Centro Nacional de Predicción Ambiental y el Centro Nacional de Investigación Atmosférica (NCEP/NCAR, por sus siglas en inglés), de la NOAA (National Oceanic and Atmospheric Administration).
Para aplicar esta técnica, se definen secciones transversales al eje de desplazamiento de la pluma volcánica en planos perpendiculares a la superficie. La formulación general del flujo de masa de SO2 a través de una superficie S de la pluma está dada por la Ecuación 3 (Theys et al., 2013):
Donde c es la concentración de masa de SO2 (kg/m3), v es el vector de viento y n es un vector unitario normal a la superficie. Integrando sobre la vertical y asumiendo un campo de viento constante, el flujo de SO2 a través de un transecto puede ser calculado por la Ecuación 4 (Theys et al., 2013):
Donde VC i es la concentración en la columna, l i es la longitud horizontal del i-ésimo píxel del perfil, v la velocida del viento y θ es el ángulo entre la dirección del perfil y la dirección del transporte. La dirección de la pluma puede ser detectada automáticamente buscando los píxeles alrededor del volcán que tienen la cantidad máxima de columna.
Se debe tomar en cuenta que la metodología de COSPEC para la estimación del flujo diario se basa en el promedio de los transectos tomados en la pluma durante el día. Para el caso de una imagen satelital, se podría suponer un número n de transectos que permitan considerar detalles que se presentan en un flujo heterogéneo. Una modificación en las ecuaciones 2 y 3 podrían expresarse respectivamente como (Ecuación 5 y Ecuación 6):
En donde la Ecuación 4 indica la estimación del promedio del flujo lo largo L de la pluma y la Ecuación 5, de forma análoga, indica la estimación del flujo promedio en m transectos dentro de la pluma volcánica.
5. Resultados y discusión
5.1. Comparación de datos COSPEC y MODIS
Los sensores COSPEC y MODIS proporcionan información complementaria que puede usarse como una herramienta de monitoreo más robusta. El sensor COSPEC presenta menos interferencia con los elementos de la atmósfera, y una mayor precisión en los datos debido a la corta distancia entre la pluma volcánica y el instrumento. Por otro lado, las imágenes satelitales del sensor MODIS muestran una cobertura espacial mucho más amplia la cual permite observar plumas volcánicas de varias decenas de kilómetros. También permite obtener información durante el día y la noche con un bajo costo de operación y una alta resolución temporal (hasta 4 imágenes de la zona del volcán diariamente).
Es importante notar que hay varias dificultades para lograr una comparación objetiva de los métodos COSPEC y MODIS. Por un lado, estamos tratando de comparar datos de concentración de SO2 obtenidos en dos ventanas de detección del espectro electromagnético en regiones diferentes (UV y TIR). Además, estos sensores se encuentran ubicados en diferentes tipos de plataformas (tierra vs. satélite) y por consecuencia la interacción de la radiación con los elementos de la atmósfera son diferentes. Trabajos previos de comparación de los dos métodos (Matiella-Novak et al., 2008), usan transectos obtenidos en la imagen MODIS en las zonas donde se llevan a cabo los pases del vehículo que transporta al COSPEC, sin tomar en cuenta el comportamiento irregular de la emisión del flujo de SO2 en el cráter del volcán. También se han desarrollado trabajos de comparación de concentración de SO2 emitida a la atmósfera, con datos obtenidos por redes de monitoreo fijas en tierra utilizando sensores DOAS (Differential Optical Absorption Spectroscopy), y concentraciones estimadas con el uso de imágenes satelitales (Spinetti et al., 2014). En el trabajo desarrollado por Spinetti et al. (2014) reportan que encontraron una diferencia del 10% al comparar los datos obtenidos en tierra con la red de sensores DOAS y los datos obtenidos con imágenes MODIS. Por otro lado, al correlacionar los datos del flujo estimado con ambas técnicas, obtuvieron un coeficiente de 0.84.
En este trabajo, se propone una metodología basada en la reconstrucción del comportamiento histórico del flujo de SO2 extraído de una imagen de satélite mostrado por Merucci et al. (2011), con el propósito de identificar las secciones de la pluma de SO2 en la imagen MODIS que fueron medidas con COSPEC. Para lograr este acoplamiento, se toma en cuenta la velocidad del viento y la diferencia de tiempo de adquisición entre las mediciones COSPEC y las imágenes MODIS. El propósito de identificar este desfase entre el tiempo de adquisición de los datos de cada sensor es para tratar de asegurar que la zona de la pluma volcánica que midió COSPEC sea identificada en la imagen MODIS. Para identificar los transectos de COSPEC dentro de la imagen MODIS es importante considerar el tiempo de inicio y fin de los transectos, así como las mediciones de la velocidad del viento. Esto permite una ubicación aproximada del área de la pluma medida con COSPEC (medición dinámica) en la imagen MODIS (medición instantánea de todo el evento). La pluma se mueve constantemente por efecto del viento predominante a la altura donde se ubica, por lo que al considerar este efecto sobre la imagen satelital se deberá considerar una relativa inclinación del transecto como se muestra en la Figura 3. El ángulo de inclinación (θ) dependerá de la velocidad del viento así como del tiempo que se haya tardado COSPEC en llevar a cabo el recorrido del transecto (Δt = t2 - t1). La resolución de cada píxel de MODIS es de un kilómetro mientras que para COSPEC la resolución es de unos pocos metros, por lo que es posible dibujar los transectos medidos con COSPEC después de viajar a través de la pluma dentro de la imagen MODIS como se muestra en las figuras 3 y 4.
La Figura 4 muestra los transectos de las mediciones hechas con COSPEC estimados en la imagen MODIS adquirida el 7 de noviembre de 2006. Los datos obtenidos en cada uno de los transectos individuales se muestran en la Tabla 1. Para el caso analizado anteriormente, se puede observar que uno de los transectos medido por COSPEC se realizó a las 18:39 UTC (tiempo medio de la medición) a una distancia de 34 km del cráter. La imagen MODIS se obtuvo a las 19:45 UTC lo que representa un Δt de 66 minutos. Debido a que la medición COSPEC se obtuvo antes de la imagen MODIS, la zona medida por COSPEC sobre la pluma se encontrara desplazada por el viento alejándola del cráter durante este Δt con una velocidad reportada de 2 m/s. Como se muestra en la Figura 4, el transecto se estimó en la imagen MODIS a 42 km del cráter del volcán. El segundo transecto tomado por COSPEC se obtuvo después de la imagen MODIS, por lo que en esta ocasión se considera que esta sección de la pluma se encontrara más cerca del cráter en la dirección del desplazamiento del viento. El transecto medido por COSPEC se localizó a 29 km del cráter del volcán en este caso se estimó una velocidad del viento superior de 4 m/s, lo que originó que el transecto de COSPEC sobre la imagen MODIS se ubicara a 14 km del cráter del volcán.
COSPEC | MODIS | ||||
---|---|---|---|---|---|
Fecha | Hora de la Medición (UTC) | SO2 (flujo ton/día) | Hora de la Imagen MODIS (UTC) | Distancia al cráter del volcán | SO2 (flujo ton/día) |
24/01/2006 | 17:14 hrs | 2487 | 17:20 hrs | 16 km | 2854 |
19:03 hrs | 1740 | 22 km | 2518 | ||
09/02/2006 | 17:48 hrs | 4656 | 20:25 hrs | 7 km | 9894 |
20:38 hrs | 4079 | 11 km | 4935 | ||
21:06 hrs | 2247 | 22 km | 2312 | ||
21/02/2006 | 17:50 hrs | 2820 | 17:45 hrs | 30 km | 3855 |
19:51 hrs | 2775 | 32 km | 3133 | ||
21:22 hrs | 3387 | 36 km | 4490 | ||
07/11/2006 | 18:39 hrs | 5335 | 19:45 hrs | 42 km | 6081 |
20:47 hrs | 7632 | 14 km | 11476 | ||
21/11/2006 | 18:01 hrs | 5576 | 19:55 hrs | 47 km | 9834 |
22:16 hrs | 2824 | 13 km | 5852 | ||
05/12/2006 | 19:06 hrs | 3625 | 20:10 hrs | 60 km | 5544 |
20:57 hrs | 3709 | 14 km | 3937 | ||
13/02/2007 | 19:13 hrs | 5449 | 19:35 hrs | 35 km | 4467 |
20:26 hrs | 6562 | 2 km | 11757 | ||
27/02/2007 | 20:12 hrs | 3006 | 19:45 hrs | 14 km | 5350 |
20/03/2007 | 18:19 hrs | 4328 | 16:55 hrs | 34 km | 5833 |
20:33 hrs | 3589 | 38 km | 5004 |
La Tabla 1 presenta el flujo de SO2 (ton/día) obtenidos para cada transecto de COSPEC y el transecto equivalente identificado en la imagen MODIS. Similar a Matiella-Novak et al. (2008), los resultados muestran que los datos obtenidos de COSPEC y los recuperados de imágenes MODIS se encuentran en el mismo orden de magnitud. Y de igual manera, el flujo de SO2 estimado en las imágenes MODIS presenta valores mayores que los estimados con los datos obtenidos con COSPEC. Esta sobre estimación en la técnica basada en la detección en el TIR se relaciona con la presencia de especies radiativamente activas, como ceniza y aerosoles, por errores en la geometría de la pluma y las condiciones atmosféricas y la disparidad entre la temperatura del suelo verdadera y aparente (Watson et al., 2004; Matiella-Novak et al., 2008).
La Figura 5a muestra los valores promedios del flujo de SO2 en cada transecto tomado en los días que se llevaron a cabo campañas de medición con COSPEC. Los valores de cada transecto de la pluma medidos en el mismo día son promediados para obtener el valor del flujo de SO2 diario de la medición. Estos valores de flujo diario promedio y su desviación estándar son mostrados en la Figura 5b, la cual describe el comportamiento del flujo de SO2 del volcán Popocatépetl durante el periodo de estudio.
La Figura 5b muestra que los valores del flujo de SO2, estimados por medio del sensor COSPEC, durante el periodo de estudio se encontraron entre un mínimo de 2300 ton/día y un máximo de 7600 ton/día. Durante este periodo de noviembre de 2006 a febrero de 2007, se llevaron a cabo 7 días de medición, lo que representa un ~6% de tiempo cubierto con mediciones de flujo de SO2 con el sensor COSPEC. Esto es debido a las dificultades logísticas para llevar a cabo las campañas de medición con COSPEC, por lo que en la actualidad se realizan una o dos campañas de medición al mes. Además, como ya se había mencionado, es posible obtener información sólo con luz de día, restringiendo aún más los intervalos de tiempo para vigilar el volcán por este método.
La Figura 6 muestra los valores del flujo de SO2 emitido por el Popocatépetl entre noviembre de 2006 a febrero de 2007. Similar a los datos de COSPEC, en esta serie de tiempo contiene patrones claramente cíclicos del comportamiento del volcán. Una gran ventaja de las imágenes MODIS es de poder construir una base de datos con una alta resolución temporal que permitió obtener datos del flujo de SO2 en el 88% de los días del periodo. Esto permite detectar cambios en el comportamiento de las emisiones del volcán que podrían ser asociados a una posible erupción explosiva (Jimenez-Escalona et al., 2011). Como se puede observar en la Figura 6, un aumento en los valores del flujo de SO2, en algunos casos, anteceden a un evento de erupción explosiva.
Se considera que los valores obtenidos de los transectos de la pluma tomados por COSPEC son comparables a los valores de estos transectos obtenidos en las imágenes MODIS. Por lo que, si ahora se considera el cálculo del flujo promedio en la pluma completa tomada en la imagen MODIS tendremos una aproximación más real con las variaciones del flujo en el volcán. Para estimar el flujo diario de SO2, se promedian los valores de flujo total obtenido en cada imagen MODIS de ese día. La Tabla 2 muestra la comparación entre los valores promedio diario de los flujos obtenidos con COSPEC y las imágenes MODIS para los días en que la información está disponible para ambas metodologías (enero, febrero, noviembre y diciembre de 2006, enero, febrero y marzo de 2007). Es importante mencionar que mientras los datos de COSPEC en cada medida mostrada en la Tabla 2 se obtuvieron en un periodo de 3-5 horas, la instantánea de una imagen MODIS puede representar el gas emitido durante 5-10 horas con una visión general del evento, teniendo en cuenta variaciones en el flujo que no pueden ser detectadas con COSPEC.
Fecha | COSPEC SO2 (flujo ton/día) | COSPEC (desv std) | MODIS SO2 (flujo ton/día) | MODIS (desv std) |
---|---|---|---|---|
24/01/2006 | 2114 | 528 | 5798 | 703 |
09/02/2006 | 3478 | 1090 | 3394 | 336 |
21/02/2006 | 2994 | 341 | 8481 | 297 |
07/11/2006 | 4846 | 2331 | 11376 | 962 |
21/11/2006 | 3630 | 1693 | 18102 | 331 |
05/12/2006 | 3667 | 59 | 2963 | 365 |
13/02/2007 | 6255 | 705 | 14194 | 403 |
27/02/2007 | 2282 | 967 | 6840 | 316 |
20/03/2007 | 3958 | 523 | 5185 | 493 |
Al comparar los valores de los flujos estimados a partir de las mediciones de ambos sensores, se identificó un sesgo sistemático entre los valores que dependen de la magnitud del flujo de SO2. La Figura 7a muestra que la diferencia entre los valores de flujo de SO2 obtenidos por ambos sensores presenta una tendencia positiva con respecto al valor del flujo estimado de las imágenes MODIS. Como se mencionó anteriormente, ambas técnicas de percepción remota (COSPEC y MODIS) tienen sus propias fuentes de error, relacionadas al tipo de radiación utilizada para detectar la presencia del SO2. Sin embargo, ambas técnicas se derivan del análisis de la absorción de radiación que cruza la nube de SO2, la cual es comparada con un valor de calibración en condiciones de ausencia del gas. Esta información transportada por la onda de radiación, es contaminada en el desplazamiento por la atmósfera hasta que ésta llega al sensor. Por lo que, al procesar la información es necesario conocer los elementos que interfieren con la señal o tratar de identificar los efectos que producen en la señal con el fin de proponer un método para eliminar o mitigar dichos elementos.
La diferencia de los valores de flujo entre los sensores varía entre 65 y 5240 (ton/día) durante el intervalo de estudio. A partir de la base de datos mostrada en la Tabla 1, se seleccionaron aleatoriamente 11 puntos para calibrar la ecuación y compensar la diferencia entre las estimaciones del flujo de SO2 con los sensores COSPEC y MODIS (ver Figura 7b). Los restantes 8 puntos se utilizaron para validar la ecuación de corrección de sesgo propuesta (véase la Figura 7c). Como se muestra en la Figura 7a, cuanto más altas son las estimaciones del flujo de SO2 de MODIS, mayor es la diferencia. Debido a la pendiente positiva que presenta la recta de ajuste de los datos, esta tendencia de sesgo se puede representar utilizando, ya sea una función exponencial o una serie de potencias (Gentle, 2009). Ambas técnicas de ajuste de curva fueron probadas para la base de datos de calibración y se observó que el error cuadrático medio (RMSE) y la desviación estándar eran menores cuando se representaba la tendencia usando series de potencia. Para obtener los coeficientes y el número de términos de la curva, se utilizaron las herramientas de MATLAB® versión 2013A. La mejor función de ajuste para representar la tendencia de sesgo es:
donde x representa el flujo de SO2 estimado utilizando imágenes MODIS en toneladas por día. Tanto para los subconjuntos de calibración como para los de validación, la aplicación de la ecuación (6) para corregir las estimaciones MODIS reduce el RMSE y la desviación estándar (ver figuras 7b y 7c). En ambos casos, el RMSE se reduce en un 33% y la desviación estándar mejora en un 40% en comparación con las estimaciones MODIS originales. Observamos que esta ecuación puede aplicarse para corregir las estimaciones de MODIS con el fin de acercarlas lo más posible a las observaciones COSPEC cuando la diferencia entre los sensores oscila entre 65 y 5240 ton/día. Para las diferencias superiores a 5240 ton/día, la ecuación necesita ser revisada con una base de datos más amplia y una nueva ecuación podría ser requerida. La Figura 7d compara las estimaciones de COSPEC, MODIS y MODIS después de eliminar el sesgo mediante el uso de la ecuación (6) para la base de datos completa. Se observa que las estimaciones del flujo de SO2 de las imágenes MODIS corregidas de sesgo son más cercanos a los valores del flujo de SO2 estimados de las mediciones de COSPEC y ambos sensores representan una tendencia similar en el tiempo. Particularmente, los picos obtenidos por MODIS redujeron su amplitud para estar en el orden de magnitud de COSPEC.
5.2. Monitoreo de las emisiones de SO2 del volcán Popocatépetl
La Figura 8 muestra las series temporales de las estimaciones MODIS de la emisión de SO2 (después de aplicar la Ecuación 6 al conjunto de datos) y las mediciones COSPEC, cuando están disponibles. En esta gráfica, los valores del flujo de SO2 estimado de las mediciones de MODIS se observan dentro de las barras de variación de los datos de COSPEC durante todo el período estudiado. La base de datos COSPEC de 2006 a 2007 muestra valores de emisión de SO2 por debajo de 10000 ton/día, mientras que las estimaciones MODIS alcanzan valores de alrededor de 40000 ton /día durante los períodos de ausencia de las mediciones COSPEC. La ventaja de contar con una herramienta que puede proporcionar datos para fechas en las que no hay datos de COSPEC es de suma importancia para contar con un seguimiento del comportamiento del volcán. Por otro lado, esta metodología presenta la ventaja de ser de muy bajo costo debido a que las imágenes MODIS son proporcionadas de forma gratuita y el procesamiento para la obtención de la información no requiere de equipo de cómputo de muy alto desempeño pudiéndose procesar en computadores de escritorio.
En los datos presentados en este estudio, se puede observar que desde comienzos de diciembre y hasta finales de febrero, no se habían llevado a cabo mediciones del flujo de SO2 para el monitoreo de la actividad del volcán Popocatépetl debido a la falta de presupuesto al final del año y la llegada tardía del presupuesto del siguiente año. Esta es una de las razones por la cual se necesita contar con herramientas alternativas que permitan el uso de información de bases de datos de imágenes satelitales. En el presente trabajo, las estimaciones MODIS de las tasas de emisión de SO2 muestran altos valores a finales de 2006 y principios de 2007. A principios de diciembre de 2006, los registros del Centro Nacional de Prevención de Desastres (CENAPRED), muestran que hubo un terremoto volcanotectónico seguido de eventos de largo período de amplitudes variables acompañados de temblores volcánicos en las siguientes fechas. El 6 de enero de 2007, el temblor armónico se observó antes de una explosión que arrojó fragmentos balísticos a los alrededores del volcán. Se produjo un temblor más armónico en los días siguientes y cesaron el 14 de enero. Sin embargo, el 25 de enero hubo otra explosión similar a la del 6 de enero. Esta última explosión ocurrió en el pico de emisión de gases medida por la imagen MODIS en un período en el que no se realizaron mediciones COSPEC entre el 22 y el 27 de enero. La metodología presentada en este trabajo para explotar las imágenes MODIS fue crucial para entender estos eventos. La información sísmica no está incluida porque está fuera de los objetivos de este documento.
6. Conclusiones
Este artículo presenta el análisis comparativo de las mediciones del flujo de SO2 estimado de las mediciones in-situ hechas con el sensor COSPEC y el flujo del SO2 derivado de las imágenes satelitales del sensor MODIS en el periodo de noviembre de 2006 a febrero de 2007 para el volcán Popocatépetl. En este estudio, se detectó una tendencia sistemática entre los datos del flujo de SO2 obtenido de las imágenes MODIS y la diferencia entre los valores de flujo obtenidos por ambos sensores. Se observó que cuando el valor del flujo estimado con los datos de MODIS incrementaba su valor, también aumentaba la diferencia con las mediciones de COSPEC, por lo que se propuso una función de corrección de este sesgo con el propósito de minimizar las diferencias entre estas dos estimaciones derivadas de sensores diferentes en el mismo evento volcánico. La función de corrección de sesgo se basa en la suposición de que los datos de COSPEC presentan una mayor precisión que las estimaciones MODIS ya que el sensor COSPEC se encuentra más cerca de la pluma que el satélite y tiene una mayor resolución espacial.
La metodología descrita en este artículo tiene en cuenta la velocidad del viento y la diferencia en el tiempo de recuperación exacta entre las dos técnicas. Las diferencias entre los datos COSPEC y las estimaciones MODIS fueron modeladas con una función de series de potencia, y después de aplicar esta ecuación, el RMSE se redujo en un 33% entre los dos conjuntos de valores. El RMSE obtenido entre los dos conjuntos fue inferior a 571 ton/día.
Al acoplar los datos del flujo de los sensores COSPEC y las estimaciones de MODIS corregidas, es posible observar que el volcán Popocatépetl emitió 5000 ton/día durante la desgasificación volcánica pasiva. Además, la serie temporal de las estimaciones MODIS permitió identificar un período de actividad volcánica (emisión frecuente de SO2), lo que sugiere la intrusión de material nuevo en el sistema interno del volcán Popocatépetl. Las estimaciones MODIS alcanzan valores de hasta 40000 ton/día durante los períodos de ausencia de las mediciones COSPEC. La ventaja de contar con una herramienta que puede proporcionar datos para fechas en las que no hay datos de COSPEC es crucial para construir una base histórica del comportamiento del volcán los más detallada posible con el propósito de comprender su propio comportamiento. La metodología presentada en este artículo muestra el potencial de acoplamiento sinérgico de dos técnicas para crear una herramienta robusta que permita un monitoreo continuo de la actividad volcánica y su implementación se puede extender fácilmente a otros volcanes.