INTRODUCCIÓN

Los incendios forestales constituyen actualmente un problema ambiental importante que con frecuencia llega a ocasionar incluso la pérdida de vidas humanas (^{PROFONAMPE 2000}), además de graves daños económicos (^{Llerena 1991}) y la reducción de áreas boscosas (^{Reátegui 1996}). Desde el punto de vista global son muy importantes las emisiones de gases y partículas a la atmósfera, entre ellos los conocidos gases de efecto invernadero (^{Blake et al. 1996}, ^{Andreae y Merlet 2001}). No obstante, actualmente las emisiones de dióxido de carbono provenientes de los incendios y quemas forestales son responsables de menos del 25 % de las emisiones totales de este gas en el mundo (^{Castillo et al. 2003}). En la región del Amazonas son muy frecuentes también los incendios producidos por la quema de vegetación durante prácticas de agricultura y desforestación (^{Alencar et al. 2006}, ^{Montenegro et al. 2011}). La quema de biomasa, además, es la mayor fuente de emisiones de partículas carbonáceas finas en la troposfera global (^{Akagi et al. 2011}).

A pesar de los daños que causa la quema de biomasa, sólo recientemente se han comenzado a analizar sus impactos atmosféricos y el transporte de contaminantes en Sudamérica (^{Mielnicki et al. 2005a} , ^{Mielnicki et al. 2005b}, ^{Boian 2005}, ^{Edwards et al. 2006}).

Tradicionalmente la modelación de los procesos químicos en los sistemas de calidad del aire se realiza independientemente de la modelación meteorológica, mediante una integración de los mecanismos químicos sobre una malla externa que contiene información sobre el transporte de masa y los campos meteorológicos (^{Fernández et al. 2010}). La ventaja de utilizar el modelo acoplado WRF-CHEM radica precisamente en que éste posee una estructura modular que permite acoplar la modelación meteorológica con los procesos químicos y considerar la interacción que existe entre ambos debido a que en la atmósfera real también los procesos químicos pueden afectar la meteorología. Por ejemplo, a través de sus efectos sobre el balance de radiación, así como la interacción de los aerosoles con las nubes a través de la nucleación de gotas (^{Grell et al. 2005}, ^{Fast et al. 2006}, ^{Chapman et al. 2009}).

El presente trabajo tiene como objetivo fundamental simular el transporte de material particulado de 10 µm de diámetro (PM₁₀), emitido por la quema de biomasa, producida por los incendios forestales que tienen lugar en Perú y en países vecinos, y que son transportadas a largas distancias en dependencia de la dirección y de la fuerza del viento. Con la simulación se persigue contar con una estimación de las zonas del Perú que son afectadas por estas partículas con el fin de implementar una red de monitoreo de las mismas. Para lograr el objetivo se utilizó el modelo Weather Research and Forecasting acoplado con el modelo de química atmosférica (WRF-CHEM).

MATERIALES Y MÉTODOS

Para la implementación del modelo se emplearon las condiciones iniciales y de frontera del Sistema Global de Predicción (GFS, por sus siglas en inglés), con una resolución horizontal de 0.5º y limitada entre los 30º y los 120º de longitud oeste y entre los 60º de latitud sur y los 30º de latitud norte. Se utilizaron los datos pronóstico de las 00 del Tiempo Universal Coordinado (UTC, por sus siglas en inglés) hasta 72 h.

Se emplearon 24 niveles de salida del GFS en el plano vertical (nivel medio del mar (nmm), 1000, 975, 950, 925, 900, 850, 800, 750, 700, 650, 600, 550, 500, 450, 400, 350, 300, 250, 200, 150, 100, 70, y 50 hPa). Las variables tomadas fueron las siguientes:

Para la construcción de los mapas promedio de altura de geopotencial se emplearon los datos del reanálisis desarrollado por el Centro Nacional de Predicciones Ambientales y por el Centro Nacional de Investigaciones Atmosféricas de los Estados Unidos de Norteamérica (NCEP-NCAR, por sus siglas en inglés) para el periodo 1979-2013, correspondiente a los meses de junio, julio y agosto (^{Kalnay et al. 1996}).

Los datos de emisión fueron tomados del sitio Detección Remota y Predicción de la Composición Atmosférica (ACRESP, por sus siglas en inglés) del Inventario de incendios (FINN, por sus siglas en inglés) del NCAR, con información derivada del Espectroradiómetro de Imágenes de Resolución Moderada (MODIS, por sus siglas en inglés), cuyas emisiones son calculadas casi en tiempo real. Se tomaron solamente los datos de emisiones diarias de material particulado del tipo PM₁₀, mostrados de forma explícita en los ficheros.

Se empleó también información del espesor óptico por aerosoles (EOA) derivada a partir de las mediciones que se realizan con un fotómetro solar del tipo CIMEL en el Observatorio "Huancayo" (12.04020 S y 75.32090 W, 3352 msmm), perteneciente al Instituto Geofísico del Perú (IGP). Los datos se tomaron para el periodo junio-agosto de 2015 con el objetivo de ser comparados con las concentraciones de aerosoles obtenidas con el modelo.

Configuración del modelo WRF-CHEM

Las simulaciones se realizaron para un dominio, cuyas características se especifican en el cuadro I. Para el experimento se seleccionó el periodo junio-agosto de 2015, durante el cual tradicionalmente se manifiesta un incremento progresivo del número de focos de incendios en la región de estudio como fue mencionado anteriormente.

CUADRO I Principales características del dominio y datos iniciales y de frontera 

*Sistema Global de Predicción

Como se indica en el cuadro I, las condiciones iniciales y de contorno del dominio externo son tomadas del pronóstico del modelo global GFS con una resolución de 0.5º, inicializado en el mismo horario para el cual se realizan las simulaciones de WRF-CHEM. En este cuadro la característica Frecuencia de CC (Frecuencia de Condiciones de Contorno) se refiere a la frecuencia con que las condiciones de contorno son capturadas de la salida del modelo base, en este caso, el GFS.

Los esquemas de parametrización utilizados se muestran en el cuadro II. Para la convección se seleccionó el esquema de Grell-Freitas ^{(Grell y Freitas 2014}). Otros esquemas de parametrizaciones que se describen en la literatura internacional son los de ^{Kain-Fritsch (1990} y ²⁰⁰⁴), de Betts-Miller (^{ECMWF 1984}), y de ^{Grell y Devenyi (2002)}, entre otros.

CUADRO II Esquemas de parametrización empleados 

*Modelo 2 de asimetría convectiva

** Modelo rápido de transferencia radiativa - global

*** Modelo de superficie de suelo Noah

Para la capa límite de la atmósfera se seleccionó el Modelo 2 de Asimetría Convectiva (ACM2, por sus siglas en inglés), descrito en los trabajos de ^{Pleim (2007}).

El modelo de radiación seleccionado fue el "modelo rápido de transferencia radiativa-global" (RRTMG, por sus siglas en inglés), usado anteriormente para radiación de onda corta y onda larga en el modelo del Centro Europeo de Pronóstico a Mediano Plazo (^{Morcrette et al. 2008}), en los modelos GFS y WRF y en el Modelo Climático ECHAM5 (^{Wild y Roeckner 2006}). RRTMG es una versión más reciente del modelo rápido de transferencia radiativa (RRTM, por sus siglas en inglés; ^{Iacono et al. 2000}) y con una mejor representación de la nubosidad. El RRTM también fue descrito por ^{Mlawer et al. (1997}) y utiliza el método k-correlacionado (^{Goody y Yung 1995}) para implementar un algoritmo que destaca frente a otros, por su rapidez.

Para la microfísica se escogió el esquema de ^{Lin et al. (1983}). Esta parametrización incluye el granizo blando, la sedimentación de hielo y el tiempo de disgregación del hielo.

El modelo de suelo utilizado es el conocido como Modelo de Superficie de Suelo Noah (Noah LSM, por sus siglas en inglés) y está basado en el Modelo de Superficie de Suelo de la Universidad del Estado de Oregón (OSU LSM, por sus siglas en inglés), desarrollado por ^{Pan y Mahrt (1987}) y empleado por ^{Chen et al. (1996}). Comparando con datos del primer experimento de campo realizado en la ejecución del Proyecto Internacional de Climatología de la Superficie Terrestre Mediante Satélites (conocido como experimento FIFE, por sus siglas en inglés), ^{Chen (1997)} encontró que el OSU LSM se comportaba a la altura de modelos más complejos. Su balance entre simplicidad y error lo llevaron a ser implementado en los modelos operacionales del NCEP y del NCAR. En ^{Chen y Dudhia (2001a} y ^2001b) puede encontrarse una descripción de la implementación del OSU LSM en el Modelo Meso Escalar de 5ta Generación (MM5, por sus siglas en inglés). El llamado Noah LSM es la evolución del OSU LSM implementada en WRF y en el Modelo Meso Escalar para Norteamérica (NAM, por sus siglas en inglés). Este modelo divide el suelo en 4 capas de grosores: 10, 30, 60 y 100 cm, hasta una profundidad total de dos metros.

El modelo químico utilizado está basado en los modelos Modelo Modal Dinámico de Aerosoles para Europa (^{Ackermann et al. 1998}) y el Modelo Secundario Orgánico de Aerosoles (^{Schell et al. 2001}) (MADE/SORGAM por sus siglas en inglés). MADE/SORGAM considera las interacciones de los aerosoles con la radiación solar y la microfísica de nubes. Se empleó el mecanismo fotoquímico RADM2 (^{NAPAP 1989}, ^{Stockwell et al. 1990}).

RESULTADOS Y DISCUSIÓN

En la figura 1 se muestran, con puntos, los focos de incendios en la región de estudio para los meses de junio, julio y agosto de 2015. Se aprecian puntos sobre Bolivia, Brasil y Perú en todos los meses, con un marcado incremento de un mes a otro. En el caso de Perú, la mayor parte de los focos se concentraron en la zona central del país, tanto zonal como meridionalmente. En el caso de Bolivia, los focos se concentran hacia la porción norte de la región representada en la figura y en el caso de Brasil hacia el suroeste del territorio enmarcado en el dominio de trabajo.

Fig. 1 Focos de incendios registrados en la región de estudio: (a) para junio, (b) para julio y (c) para agosto 

La figura 2 muestra la distribución de las concentraciones medias mensuales de PM₁₀, generadas a partir de las salidas del WRF-CHEM, para los tres meses analizados. Para el mes de junio se aprecian máximos sobre el suroeste de Brasil, norte de Bolivia y la zona central de Perú, en correspondencia con las regiones con mayor número de incendios mostrados en la figura 1. Para los meses de julio y agosto estos máximos se incrementan cuantitativamente, a la vez que sus zonas de impacto se amplían, lo cual se corresponde también con el incremento del número de incendios para estos meses.

Fig. 2 Distribución de las concentraciones medias mensuales de PM₁₀ (µg/m3), calculadas a partir de las salidas del modelo de pronóstico e investigaciones sobre el tiempo, acoplado con el modelo de química atmosférica (WRF-CHEM): (a) junio, (b) julio y (c) agosto 

La figura 3 muestra los campos de altura de geopotencial medios entre los 1000 y los 850 hPa (troposfera baja) para el período 1979-2013, correspondientes a los meses de junio, julio y agosto. Se puede apreciar un centro permanente de altas presiones en la porción sur del territorio de Bolivia ("A" en azul) y una zona de bajas presiones en el oeste de Brasil ("B" en rojo), que originan un corredor de aire orientado hacia las regiones central y norte de Perú. Lo anterior sugiere un transporte de aerosoles hacia Perú desde los puntos de incendios ubicados en Brasil y Bolivia.

Fig. 3 Campos de altura de geopotencial medios entre los 1000 y los 850 hPa. (troposfera baja) para el periodo 1979-2013: (a) corresponde a junio, (b) a julio y (c) a agosto. Las zonas de alta presión se indican con "A" en azul y las de baja presión con "B" en rojo 

Cuando se analiza el trimestre de forma integral (Fig. 4b) se confirma un flujo promedio desde el sureste hacia el noroeste, que sugiere un transporte promedio de partículas emitidas por la quema de biomasa en Bolivia y Brasil hacia territorio peruano. Aquí la letra "B" en rojo indica un centro de bajas presiones, y la "A" en azul, un centro de altas presiones.

Fig. 4 Análisis para el trimestre junio-agosto de 2015: (a) distribución de los focos de quema de biomasa en la región de estudio, (b) campo de altura de geopotencial medio entre los 1000 y los 850 hPa para el periodo 1979-2013 y (c) distribución de la concentración media trimestral (µg/m3). En (b) "A" en azul es una zona de alta presión y "B" en rojo, un centro de baja presión 

A modo de comprobación, se realizó una simulación para el 24 de agosto de 2015 a las 00 UTC (Plazo 54 h), incluyendo las emisiones producidas en Perú y fuera de este, considerando solamente las emisiones producidas en Brasil y Bolivia e incluyendo las emisiones producidas solamente en Brasil (Fig. 5). La figura 5b sugiere que cuando no se incluyen focos ubicados en Perú, se recibe la influencia de las emisiones de focos ubicados en Brasil y Bolivia, principalmente hacia el centro (y en menor medida hacia el norte) del territorio peruano. Esto se corresponde con la dirección del flujo predominante del Este-Sudeste mostrado en las figuras 3b y 4b. Sin embargo, en la figura 5c se aprecia que mayoritariamente la influencia se recibe de Brasil, mientras que desde Bolivia es poco significativa.

Fig. 5 Concentraciones de PM₁₀ (µg/m3) obtenidas con el modelo, inicializado el 24 de agosto de 2015 (plazo 54 h de pronóstico): (a) incluye los focos ubicados en Perú y fuera de este, (b) incluye solamente los focos ubicados en Brasil y Bolivia y (c) incluye solamente los focos ubicados en Brasil 

El cuadro III ilustra los porcentajes de aporte a las concentraciones de PM₁₀ en diferentes localidades, de los incendios producidos en Perú y de los que ocurrieron en Brasil. Como se puede apreciar, el mayor aporte lo realizaron los incendios ocurridos en Perú, aunque hay un determinado aporte de los que tienen lugar en Brasil, más significativo en localidades cercanas a las fronteras, como son Pucallpa y Puerto Maldonado. De este cuadro se deduce que el aporte de los incendios producidos en Bolivia es insignificante, con solo 0.5 % en Puerto Maldonado.

CUADRO III Aporte de las emisiones de quema de biomasa producidas en Perú y Brasil a las concentraciones de PM₁₀ para diferentes localidades según inicialización del 24 de agosto de 2015 a las 00 UTC* 

*Tiempo Universal Coordinado

En la figura 6 se representan las curvas normalizadas de los valores medios diarios de las concentraciones de PM₁₀ obtenidos con el modelo y del espesor óptico por aerosoles registrado para el período junio-agosto para la localidad de Huayao, donde se encuentra ubicado el observatorio "Huancayo". Se observa que de forma general existe una buena correspondencia entre ambas curvas para los tres meses, con un coeficiente de correlación del 72.4 % para el trimestre y de 76.8 % para el mes de agosto. Desde inicios del mes de junio hasta mediados de agosto se aprecia que tanto el EOA como las concentraciones tienen un comportamiento relativamente estable con un ligero crecimiento, sin embargo, a partir de ese momento comienza un incremento significativo de ambas variables.

Fig. 6 Curvas normalizadas del espesor óptico por aerosoles (EOA, azul) y de concentraciones de PM₁₀ (roja) obtenidas con el modelo para los meses de junio a agosto de 2015 

En la figura 7 se representan los porcentajes correspondientes a la clasificación de los diferentes tipos de aerosoles, para cada mes de los aquí analizados. En esta clasificación entran los aerosoles de tipo marítimo, continental, polvo, mezcla y contaminado. La mencionada clasificación se ha realizado de acuerdo con los criterios de ^{Hess et al. (1998}), Hoben et al. (2001) y ^{Dubovik et al. (2002}).

Fig. 7  Porcentaje de las clasificaciones de los diferentes tipos de aerosoles registrados por el fotómetro solar ubicado en el observatorio "Huancayo": (a) para junio, (b) para julio y (c) para agosto 

En correspondencia con los análisis anteriores se puede apreciar el incremento paulatino de los aerosoles de tipo contaminado a lo largo de los tres meses. Los aerosoles de tipo continental consideran la contribución de especies carbonáceas originadas fundamentalmente como consecuencia de contaminación de ciudades e industrias, básicamente, como consecuencia de la quema incompleta de combustibles fósiles. En el presente estudio los aerosoles de tipo continental constituyen básicamente el EOA de fondo en la zona de estudio, con un valor medio de 0.088 para el trimestre analizado. Sin embargo, el EOA del tipo contaminado alcanza un valor medio de 0.210 para el mismo periodo. Mientras que en los meses de junio y julio los mayores porcentajes corresponden a los aerosoles del tipo continental con valores superiores a 70 %, en el mes de agosto existe un predominio de los aerosoles del tipo contaminados, con un 48 %. Los aerosoles de tipo contaminado están compuestos básicamente por especies carbonáceas, las mismas que forman parte de los aerosoles generados como consecuencia de la quema de biomasa según ^{Dubovik et al. (2002}). Este análisis por tipos de aerosoles corrobora, por tanto, los resultados obtenidos a partir de las simulaciones realizadas con el WRF-CHEM, sobre todo si se tiene en cuenta que la información de entrada al modelo WRF-CHEM está constituida solamente por las emisiones de PM₁₀ en los puntos de incendio.

CONCLUSIONES

Los resultados del trabajo arrojaron que sobre toda la región de estudio tiene lugar un incremento progresivo, desde junio hasta agosto, del número de focos de incendios, lo cual dio lugar a un aumento, también progresivo, de las concentraciones de partículas PM₁₀ sobre todo el territorio y particularmente sobre Perú. Se pudo definir que el territorio peruano es afectado, tanto por partículas generadas por incendios dentro del país, como por partículas provenientes de incendios ocurridos en países vecinos, como Brasil y Bolivia, este último con una contribución muy baja en relación con Brasil y con el propio Perú.. El estudio de los patrones sinópticos medios de latitudes bajas para el período estudiado confirmó que en esta época del año predomina un flujo de aire orientado desde el sudeste hacia el noroeste, que garantiza el transporte de partículas desde los países antes mencionados hacia el Perú. Un estudio de caso simulado arrojó que Brasil realizó mayor aporte que Bolivia, mientras que las mayores concentraciones de PM₁₀ en Perú resultaron de incendios que tienen lugar en su el propio territorio.

El EOA registró un incremento progresivo a partir de mediados de julio, que fue más notable durante el mes de agosto, en correspondencia con el incremento de las concentraciones de PM₁₀. Los análisis realizados arrojaron un incremento de material carbonáceo, más notable en agosto, típico de las emisiones de partículas producidas por quema de biomasa. De todo lo anterior se infiere la necesidad de establecer una red de monitoreo de concentraciones de PM₁₀, fundamentalmente en las regiones central y norte de Perú, dada la importancia que su control significa para la salud de la población.