Evaluación del transporte de contaminantes atmosféricos en La Habana mediante el uso del modelo HYSPLIT

Caracterización del área de estudio

La Habana (Figura 1) es la ciudad más poblada de Cuba, así como el mayor centro industrial y portuario del país, contando con muchas industrias fuertemente contaminantes (refinerías, termoeléctricas, fábrica de cemento, farmacéuticas, etc.) en su territorio, cuyas emanaciones a la atmósfera se depositan posteriormente en el suelo. El polvo acumulado en zonas urbanas contiene una diversidad de materiales que van, desde compuestos minerales hasta materiales orgánicos e inorgánicos de origen antropogénico, los que pueden ser depositados en la superficie de las calles y techos. Los metales pesados de la superficie del suelo y el polvo de las calles son buenos indicadores de la acumulación de PM en el entorno urbano (Díaz Rizo et al., 2019Díaz Rizo, O., Olivares Rieumont, S., Rudnikas, A. G., D´Alessandro Rodríguez, K., Lima Cazorla, L., García Céspedes, D., Casanova Díaz, A. O., García Trápaga, C., y Manduca Artiles, M. (2019). Estudios de contaminación ambiental en La Habana mediante técnicas nucleares y conexas. Nucleus, 66, 58-65. http://nucleus.cubaenergia.cu/index.php/nucleus/article/view/693 ).

Modelo HYSPLIT

En el sistema HYSPLIT se puede observar el transporte y dispersión de los contaminantes que se originan en una variedad de fuentes. Este simula la dispersión y la trayectoria de las sustancias a través de nuestra atmósfera, desde escalas locales hasta globales (HYSPLIT, s. f.HYSPLIT. (s.f). Air Resources Laboratory. https://www.arl.noaa.gov/hysplit ).

Es un sistema completo para calcular trayectorias simples de paquetes aéreos, así́ como simulaciones complejas de transporte, dispersión, transformación química y deposición. Una de las aplicaciones más comunes de este modelo es el análisis de trayectoria trasera para determinar el origen de las masas de aire y establecer relaciones fuente- receptor. De igual manera dicho modelo también se ha utilizado en una variedad de simulaciones que describen el transporte atmosférico, la dispersión y la deposición de contaminantes y materiales peligrosos. Algunos ejemplos de las aplicaciones incluyen el seguimiento y pronóstico de la liberación de material radiactivo, humo de incendios forestales, polvo soplado por el viento, contaminante de diversas fuentes de emisión estacionarias y móviles, alérgenos y ceniza volcánica (Stein et al., 2015Stein, A. F., Draxler, R. R., Rolph, G. D., Stunder, B. J. B., Cohen, M. D., y Ngan, F. (2015). NOAA’s HYSPLIT Atmospheric Transport and Dispersion Modeling System. Bulletin of the American Meteorological Society, 96(12), 2059-2077. https://doi.org/10.1175/bams-d-14-00110.1 ).

Gracias a la sencillez, rapidez y efectividad del modelo de trayectoria se pueden establecer las posibles direcciones que llevan las masas de aire y aplicar las decisiones oportunas para el cuidado de la población. Las retrotrayectorias del modelo, combinadas con imágenes de satélite (por ejemplo, de los satélites MODIS de la NASA), pueden proporcionar información sobre si los altos niveles de contaminación del aire son causados por fuentes locales de contaminación del aire o si un problema de contaminación del aire fue arrastrado por el viento.

Casos de estudio

Mediciones in situ

En el año 2019 fue realizada una campaña en la estación de Cubaenergía. Para las mediciones de PM realizadas fue utilizado el muestreador de bajo volumen con cabezales de cortes independientes para el muestreo de las partículas PM2.5 y PM10: TCR TECORA ECHO PM. Dicho equipo cumple con los procedimientos de muestreo regulados por las normas europeas UNE EN 12341: 1999 y UNE EN 14907: 2006, para el monitoreo de partículas PM2.5 y PM10 respectivamente.

En las mediciones fueron utilizados un total de 42 filtros para ambos tipos de partículas con el fin de determinar las concentraciones de estas y apreciar qué días excedían los límites diarios establecidos según las normas cubanas. Para este caso se calcularon las retrotrayectorias y frecuencias de retrotrayectorias, para determinar la procedencia de las partículas colectadas en el sitio de muestreo, principalmente aquellas de origen regional, y de esta manera comprobar la existencia de eventos regionales o locales, además de fuentes fuera de las existentes normalmente (Figura 2), que pudieran incidir en las concentraciones elevadas de PM obtenidas para esos días.

Para este caso se analizaron retrotrayectorias y frecuencia de retrotrayectorias para 72 horas que tuvieron como punto final el sitio Cubaenergía para los días donde las concentraciones de PM excedieron los límites diarios establecidos por las normas cubanas. Y para la meteorología se utilizaron los datos del modelo Sistema global de Asimilación de Datos (GDAS).

Datos satelitales

Fueron descargadas imágenes de la plataforma Copernicus, correspondientes a las columnas de dióxido de azufre (SO₂) (European Union’s Space programme, 2023aEuropean Union’s Space programme(2023). Copernicus Sentinel-5P data products. Sentinel Online. https://sentinels.copernicus.eu/web/sentinel/data-products ) captadas por el sensor TROPOMI del satélite Sentinel 5P (European Union’s Space programme, 2023bEuropean Union’s Space programme(2023). Sentinel-5P. Sentinel Online. https://sentinels.copernicus.eu/web/sentinel/missions/sentinel-5p ).

De igual forma se obtuvieron imágenes para el espesor óptico. Estas fueron captadas por el sensor Modis, del satélite Aqua, descargadas todas de la plataforma Giovanni NASA (NASA, 2023NASA. (2023). Giovanni. https://www.usillife.pe/investigacion-y-desarrollo/web-giovanni-nasa/ ).

Resultados para caso de estudio de erupción del volcán Cumbre Vieja

En la Figura 3 y la Figura 4 se muestran las posibles trayectorias obtenidas por el modelo Hysplit, que seguían las cenizas expulsadas durante la erupción para 24 horas y cada 10 días respectivamente.

Para la totalidad de los días analizados las trayectorias presentan direcciones Este-Oeste aproximadamente, excepto para las semanas del 19 al 28 de septiembre de 2021, del 29 de octubre al 8 de noviembre de 2021 y del 19 de noviembre al 28 de noviembre de 2021, donde las masas de aire siguen direcciones Este-Noreste aproximadamente (Figura 4). Todos estos contaminantes sobrepasaron las alturas iniciales que utilizó el modelo en su cálculo y viajaron principalmente sobre los 1500 m de altura.

En la Figura 5 se obtuvo la posible dispersión que presentaban las cenizas expulsadas en la erupción del volcán Cumbre Vieja, para 24 horas, a partir de los días 19 de septiembre, 1 de octubre, 1 de noviembre y 1 de diciembre.

Fueron estimadas 1000 toneladas por día basado en publicaciones de la prensa internacional. Utilizando como altura sobre el nivel del mar del volcán 6114 pies.

Las mayores concentraciones se estimaron al sureste de la isla con valores de 1 mg/m³ y se dirigen hacia el oeste con respecto a la fuente de emisión.

En las dispersiones analizadas se puede observar la llegada de la pluma de cenizas y demás contaminantes expulsados hasta el norte de África y el Caribe, coincidentemente con los modelos de trayectoria anteriormente presentados para este caso, en direcciones Este-Oeste mayormente.

Resultados para caso de estudio Central Termoeléctrica Ernesto Che Guevara.

Como se puede observar en la Figura 6, las salidas del modelo de trayectoria para el día 15 de octubre de 2022, la pluma de contaminantes expulsados por la CTE Ernesto che Guevara, en todos los casos se mueve en dirección Oeste.

Se obtuvo la posible dispersión que presentaban los contaminantes expulsados por la CTE Ernesto Che Guevara para el día 15 de octubre 2022, para dos horarios de máxima emisión 0500 UTC y 1400 UTC. (Figura 7)

Se estimaron las emisiones de la CTE, considerando factores de emisión obtenidos. Las mayores concentraciones se registraron hacia el oeste de la fuente de emisión, con valores inferiores a 0,01 mg/m³.

Caso de estudio de datos in situ recopilados en la estación de CUBAENERGÍA

Con los datos que fueron recogidos de las mediciones realizadas en la estación de CUBAENERGÍA se realizaron los siguientes gráficos:

Los días que excedieron las normas cubanas establecidas de PM2.5 (50 µg/m³) y PM10 (25 µg/m³), así como los valores de las concentraciones correspondientes se arrojan en las siguientes tablas:

Seguidamente para dichas fechas se analizaron las retrotrayectorias (Figura 8) y sus frecuencias (Figura 9) para analizar la procedencia de estos contaminantes.

En este caso la frecuencia de trayectoria inició una trayectoria desde una sola ubicación (CUBAENERGÍA) y altura cada 6 horas para 6 días, y luego sumó la frecuencia con que la trayectoria pasó por dicha ubicación. En el caso de las retrotrayectorias solamente fueron calculadas para 72 horas, por tal motivo no existe una coincidencia exacta entre las direcciones que muestra el modelo de trayectoria y las frecuencias de trayectoria, específicamente para los días 17 de mayo del 2019, donde las frecuencias de trayectoria indica direcciones del Norte y para el día 26 de junio del 2019, donde indica direcciones del Este.

Resultados de los datos satelitales

Datos in situ de la estación de Cubaenergía

En este caso de estudio solo resultó significativo la influencia de una nube de polvo del Sahara en las mediciones de junio de 2019, un evento muy característico en esta época del año, como se muestra en la Figuras 13, tomadas del satélite Sentinel, en la plataforma Copernicus.

Comparación de las salidas del modelo con las imágenes satelitales

El cálculo de las trayectorias, dispersión y retrotrayectorias en combinación con imágenes satelitales y mediciones in situ ayuda al diagnóstico de la contaminación, aportando una visión global de la zona de estudio, indicando la fuente de la contaminación y las zonas afectadas donde se deben implementar las medidas de prevención y control.

En el caso de la erupción del Volcán Cumbre Vieja gracias a la correspondencia que existe entre las imágenes satelitales y la salida del modelo de trayectoria para 10 días, se puede observar que en la semana correspondiente al 29 de septiembre de 2021 hasta el 8 de octubre de 2021 debe ocurrir la llegada a Cuba, y a territorio habanero, de la pluma de cenizas. Esta conclusión debida a que el Hysplit es un modelo langrangiano que calcula los movimientos de las masas de aire en el tiempo.

Se puede constatar la llegada del gas y la pluma de contaminantes a Cuba sobre los días 4 y 5 de octubre de 2021 trasladada por una nube de polvo del Sahara, donde fue predominante sobre nuestro país y sobre territorio habanero la existencia de partículas.

Teniendo en cuenta las alturas que se utilizaron para el cálculo de las trayectorias, se estimó́ la mayor concentración de contaminantes sobre los 1,5 Km de altura, pero los mismos, debido a la situación sinóptica de los vientos que se encontraban en dirección Noreste-Este, la provincia de la Habana no fue una de las más afectadas por este evento regional.

Para la CTE Ernesto che Guevara, como se observó́ en la Figura 6, todas las trayectorias que se realizaron constatan que las masas de aire se desplazan mayormente en direcciones Oeste-Noroeste, y como se observa en la imagen satelital (Figura 12), no existe gran concentración de partículas sobre la provincia de La Habana para la fecha del 15 de octubre 2022. Con estas comparaciones se observa que para dicha fecha no existe un gran evento de contaminación afectando a La Habana proveniente de esta ubicación. No se descarta, sin embargo, su aporte a la contaminación reinante en la provincia de La Habana.

Con respecto a las mediciones in situ, en este caso se encontró, como fue mencionado anteriormente, solamente la incidencia de un evento de polvo del Sahara en las mediciones de junio 2019, mientras que para el resto de los días que excedieron las normas establecidas para las concentraciones de material particulado se intuye la incidencia de otras fuentes de contaminación, como puede ser: la quema de combustibles fósiles, el alto nivel de tráfico vehicular o la salinidad presente en el ambiente, producto de los mares que rodean a Cuba.