Spanish (pdf)
Article in xml format
Article references
Automatic translation
Send this article by e-mail
Cited by SciELO 
Cited by Google 
Similars in
SciELO 
Similars in Google 
Permalink
1. INTRODUCCIÓN
Los aerosoles son partículas finas en estado sólido o líquido que se encuentran suspendidas en un gas; se originan de forma natural (biogénicos) o por actividades humanas (antropogénicos) y su rango de tamaño está en el orden de unos cuantos nanómetros hasta decenas de micras [1]. Un componente que puede dar origen a los aerosoles y que se emite a la atmósfera es el contaminante denominado Black Carbon (BC), que incluye componentes como Elemental Carbon (EC) y Organic Carbon (OC) [2], material carbonáceo producido por la combustión incompleta de biomasa, combustibles fósiles y biocombustibles [3],[4].
El BC se caracteriza por absorber radiación en todas las longitudes de onda y se ha establecido que puede causar calentamiento de la atmósfera cuando se encuentra suspendido en el aire, y derretimiento de la nieve o el hielo tras su deposición [5]. Varios autores han estudiado la contribución del BC al cambio climático, concluyendo que tiene un efecto neto de calentamiento, ocupando junto al CO2 y el metano, los primeros lugares entre los contaminantes que mayor aporte realizan al calentamiento global [5]-[6].
Los resultados acerca de la contribución del BC al calentamiento global se basan en su mayoría en cálculos de forzamiento radiactivo, es decir del balance entre la radiación solar entrante y la infrarroja saliente. El Panel Intergubernamental sobre el Cambio Climático (2007) (IPCC por sus siglas en inglés) reportó tres tipos de forzamientos que al combinarse dan como resultado un forzamiento neto, dichos tipos son: Forzamiento directo, que se refiere a la absorción directa de la radiación solar o terrestre; forzamiento albedo nieve/hielo, resultado del cambio en la reflectividad de la nieve o el hielo a causa de su oscurecimiento por el BC; y forzamiento indirecto, resultante del impacto sobre las características de las nubes, como tiempo de vida, composición, reflectividad, etc.
El IPCC [5], Bond [7], entre otros, reportan valores diferentes para el forzamiento radiactivo ejercido por el BC, evidenciando que existe aún mucha incertidumbre sobre la magnitud del mismo [2].
Además de su impacto sobre el clima y los fenómenos meteorológicos, se ha encontrado que el BC causa problemas en la salud de las personas, afectando principalmente el sistema respiratorio debido a que es capaz de alcanzar los alveolos pulmonares al ser inhalado [6],[8].
Existe también evidencia de que personas expuestas a emisiones provenientes de motores diésel, las cuales contienen BC, pueden desarrollar cáncer de pulmón [9].
De acuerdo con estudios realizados, se ha establecido que las emisiones de BC en países desarrollados han tenido una importante disminución en los últimos años, en parte debido a mejores sistemas de control y combustibles más limpios; mientras que en países en vía de desarrollo de América Latina, Asia y África ha ocurrido lo contrario. Se espera que estas tendencias continúen en las próximas décadas [10]. En el caso particular del Área Metropolitana del Valle de Aburrá, más de dos tercios del total de emisiones de PM2.5 (del cual forma parte el BC) provienen de fuentes móviles, siendo el empleo del combustible diésel el responsable del 80 % de las emisiones de este contaminante [11].
En vista de la importancia que reviste el BC tanto para la salud humana, como para el calentamiento global, y tomando en cuenta que se trata de un contaminante de corta vida, es decir, que sus efectos sobre el clima son casi inmediatos, así como los beneficios de su mitigación [7],[8], se hace necesario conocer sus principales fuentes, su comportamiento en la atmósfera y su tiempo de vida en la región del Área Metropolitana del Valle de Aburrá.
Con el fin de dar cumplimiento a estos propósitos, se decidió emplear el Modelo Integrado de Calidad del Aire con Extensiones (CAMx, por sus siglas en inglés). De acuerdo con [11], CAMx es un modelo fotoquímico de dispersión euleriano que posibilita la evaluación de la contaminación gaseosa y particulada en diferentes escalas. Es capaz de simular la emisión de contaminantes, su dispersión, reacciones químicas y remoción de la atmósfera en un sistema tridimensional de mallas anidadas, a través de la solución de la ecuación de continuidad para cada especie.
Para la simulación de resultados se aplicó el Sistema de Modelización Atmosférico Regional Brasilero (BRAMS) que simula la meteorología, y la actualización del inventario de emisiones para el Valle de Aburrá año 2015, elaborado por la Universidad Pontifica Bolivariana (U.P.B.) y el Área Metropolitana del Valle de Aburrá (autoridad ambiental de la región).
En el presente artículo se presenta la estimación de emisiones de BC en el Área Metropolitana del Valle de Aburrá a partir del inventario de emisiones de esta región (año 2015), la dispersión de este contaminante en la zona de estudio y la evaluación de los modelos (CAMx y BRAMS), comparando los resultados de simulación con los resultados medidos.
2. METODOLOGÍA
2.1 Configuración del modelo CAMx
2.1.1 Descripción del modelo
El campo de dispersión de las partículas de BC se simulan usando el Modelo Integrado de Calidad del Aire con Extensiones (CAMx) que dispone de un módulo de especiación de varias especies químicas, incluyendo las partículas finas, específicamente la denominada Primary Elemental Carbon (PEC).
El modelo está diseñado para reunir todas las características técnicas especificadas en el estado del arte de los modelos de calidad del aire, todo en un sistema simple open-source (libre) que es computacionalmente eficiente, flexible y que está disponible públicamente [12].
Adicionalmente, se empleó el mecanismo de reacción químico Carbon Bond 6 (CB06) el cual describe la química oxidante troposférica de una manera concisa y adecuada para el uso en modelos atmosféricos complejos de tres dimensiones y el esquema de aerosol CF (Coarse/Fine) para representar la distribución de tamaño de partícula.
El modelo CAMx fue ejecutado con un dominio de resolución horizontal de 60 x 60 km2 y 13 niveles verticales. Este dominio cubre la región del Valle de Aburrá y parte del oriente en el departamento de Antioquia, tal y como se muestra en la Figura 1.
2.1.2 Área de estudio
El área de estudio se encuentra ubicada en el departamento de Antioquia, al noroccidente de Colombia y abarca la región del Área Metropolitana del Valle de Aburrá con sus diez municipios (Caldas, La Estrella, Sabaneta, Envigado, Itagüí, Medellín, Bello, Copacabana, Girardota y Barbosa), así como algunos municipios de la región del Oriente Antioqueño (Figura 1).
Es un valle con una altura promedio sobre el nivel del mar de 1480 metros, rodeado por montañas que alcanzan los 2800 msnm y una topografía pendiente e irregular. Además, los municipios que conforman la región pertenecen a la cuenca hidrográfica del río Aburrá o Medellín, que atraviesa el valle de sur a norte.
2.1.3 Condiciones iniciales
Las condiciones iniciales para el modelo fueron las sugeridas por Christopher Emery de ENVIRON [12] y se presentan en la Tabla 1.
Los parámetros meteorológicos alimentados al modelo se obtienen del Sistema de Modelización Atmosférico Regional Brasilero (BRAMS) versión 5, simulando la circulación atmosférica sobre un área geográfica limitada [13].
2.1.4 Emisiones de fuentes contaminantes
El Área Metropolitana del Valle de Aburrá, viene desarrollando desde el año 2000 el inventario de emisiones en convenio con la Universidad Pontificia Bolivariana. En el presente estudio se utiliza la última versión reportada, la cual es con año base 2015 [11] y se usa como datos de entrada al modelo fotoquímico CAMx.
El inventario tiene en cuenta las emisiones de fuentes de área, fijas y móviles de la región, que son estimadas a partir de factores de emisión (FE) reportados en el AP-42 Compilation of Air Pollution Emission Factor de la Agencia de Protección Ambiental de Estados Unidos (EPA) [11] para las dos primeras fuentes de emisión, mientras que para las móviles, se utilizan los FE reportados en el modelo Internacional de Emisiones Vehiculares (IVE por sus siglas en inglés) [14].
2.2 Mecanismo de reacción
El mecanismo de reacciones químicas que usa el modelo es el Carbon Bond (CB) que describe la química oxidante troposférica de una manera concisa y adecuada para el uso en modelos atmosféricos complejos en tres dimensiones. Las versiones existentes (CB4 y CB05) han sido usadas en modelos fotoquímicos como CMAQ, CAMx y WRF-Chem. En este estudio, se empleó el mecanismo CB6, con el cual se mejora la formación de oxidantes con la inclusión de varios compuestos orgánicos que son de larga duración y relativamente abundantes como propano, acetona, benceno y etino (acetileno). Adicionalmente, se mejora el modelizado de Aerosoles Orgánicos Secundarios (SOA) con la inclusión de compuestos alfa-dicarbonilo (glioxal y análogos) y la formación de peróxido de hidrógeno (H2O2) y, por lo tanto, la formación de aerosoles de sulfato con la actualización en la química de los radicales peroxi [12].
Las 213 sustancias químicas con las que cuenta el inventario de emisiones del Valle de Aburrá son agrupadas por especies funcionales que producen 39 ecuaciones de especies reactantes, las cuales hacen finalmente el mecanismo de reacción del modelo CAMx.
Además de utilizar dicho mecanismo de reacción, se cuenta con la herramienta para la especiación de partículas finas obtenida de la documentación de desarrollo de bases de datos de especiación de la EPA - Speciate Version 4.2 (EPA), donde se presentan las distribuciones (porcentaje en peso) de especies químicas para material particulado.
2.3 Escenario modelizado
Se realizó una simulación empleando el modelo CAMx para conocer la dispersión de BC en el Valle de Aburrá empleando el inventario para el día 9 de noviembre de 2016, el cual tiene condiciones meteorológicas favorables (verano) para una adecuada dispersión de contaminantes con alta radiación y baja nubosidad a baja altura.
El día mencionado fue evaluado en el punto donde se encuentra la estación automática de monitoreo de la calidad del aire ubicada dentro de la Universidad Pontificia Bolivariana (75°35’27,08’’O, 6°14'36,85"N), debido a que es una zona poblacional expuesta en donde es posible de evaluar el desempeño del modelo pues se cuenta con una campaña de medición del contaminante estudiado.
2.4 Evaluación del escenario
El Grupo de Higiene y Gestión Ambiental (GHYGAM) perteneciente al Politécnico Colombiano Jaime Isaza Cadavid, en convenio con el Área Metropolitana del Valle de Aburrá, realizó una campaña de medición de BC en la zona Laureles, específicamente en la Universidad Pontificia Bolivariana, durante 11 días del mes de noviembre de 2016 (del 2 al 5, del 9 al 12 y del 16 al 18).
El monitoreo se llevó a cabo usando un medidor Automático de BC conocido en el mercado como Etalómetro Portable Magee Scientific Modelo AE42, el cual reporta en tiempo real a partir de la cuantificación de la tasa de cambio en la absorción de la luz transmitida, debido a la continua recolección del filtro de la deposición de aerosol. En el rango 880 nm es evaluada como medición directa de la concentración de BC pues representa la materia carbonácea tipo antracita que absorbe energía calórica [11]. Los datos recolectados en la campaña fueron usados para evaluar el desempeño del modelo CAMx.
En cuanto al análisis de desempeño del modelo meteorológico, se emplearon los datos de temperatura y humedad relativa proporcionados por la estación automática de monitoreo de la calidad del aire ubicada dentro de la Universidad Pontificia Bolivariana, zona Laureles, la cual, además, mide contaminantes como SOx, NOx, O3 y PM2.5.
3. RESULTADOS Y DISCUSIÓN
Una vez configurado el dominio computacional del modelo que contiene las variables topográficas y la definición del mecanismo de reacción química, se alimenta al modelo con las emisiones de BC con una desagregación espacial de 1 x 1 km2 y una resolución temporal de una hora.
3.1 Estimación de emisiones de BC
A partir de las emisiones de PM2.5 calculadas para fuentes fijas y móviles, provenientes del inventario de emisiones del Valle de Aburrá año 2015, se estimaron las respectivas emisiones de BC para cada una de las fuentes contaminantes, tal como se discute a continuación.
3.2 Fuentes móviles
Para las emisiones de BC provenientes de las fuentes móviles, se utilizó la documentación de desarrollo de bases de datos de especiación de la EPA - Speciate Version 4.2 (EPA), donde se presentan las distribuciones (porcentaje en peso) de especies químicas para material particulado. Es así como, a partir de esta distribución y de la emisión de PM2.5 calculado, es posible estimar las emisiones de BC.
En la Tabla 2 se presentan las emisiones anuales de BC discriminadas por categoría vehicular. Se observa que la categoría camiones seguido de las volquetas son los mayores emisores aportando el 43 % y 37 %, respectivamente, del total de las emisiones asociadas a fuentes móviles.
Las emisiones anuales por tipo de combustible se presentan en la Tabla 3, donde se observa que el diésel es el principal responsable por las emisiones de BC con un 95 % del total de las emisiones asociadas a fuentes móviles, el cual, a pesar de su versatilidad y su eficiencia, históricamente se ha evidenciado que producen altas emisiones de NOx y material particulado. Adicionalmente, se observa que la gran mayoría de las emisiones de BC proviene de las categorías Camiones y Volquetas, que principalmente utilizan el combustible diésel.
3.3 Fuentes fijas
Utilizando la clasificación establecida en la Tabla 4, se muestran los resultados en la Tabla 5 de las emisiones anuales de BC por actividad productiva asociada con las fuentes fijas propias de la región.
A partir de los resultados de emisiones obtenidos para las diferentes actividades industriales, se concluye que el sector Textil es el que más aporta en emisión de este contaminante, seguido de la industria maderera con un 59 % y 30 %, respectivamente, del total de las emisiones de BC asociadas a fuentes fijas. Los demás sectores productivos aportan muy poco a la emisión de este contaminante por lo que no se incluyeron en la Tabla 5.
En la Tabla 6 se muestran las emisiones anuales por municipio, en donde se observa que el municipio de Barbosa realiza los mayores aportes de BC seguido del municipio de Itagüí con un 32 % y 27 %, respectivamente, del total de las emisiones en el Valle de Aburrá. Los demás municipios del Valle de Aburrá aportan muy poco a la emisión de este contaminante.
La razón por la que el sector textil y maderero son los que más aportan a la emisión del BC, así como el municipio de Barbosa, se debe a la distribución de las fuentes por sector en cada municipio.
Debido a que en el municipio de Barbosa se encuentran las 3 empresas más grandes del sector maderero, éstas aportarían la mayor carga de BC en el municipio.
Mientras que, para el caso del sector textil, si bien emite el doble al maderero, esta industria se encuentra asentada en los diferentes municipios del Valle de Aburrá, principalmente en el municipio de Itagüí, siendo éste el segundo municipio emisor más importante.
3.4 Comparación inventario de BC
En total, en la región del Área Metropolitana del Valle de Aburrá para el año 2015, se emitieron 850 toneladas de BC, de las cuales el 97 % pertenecen a las fuentes móviles, especialmente los vehículos que usan diésel, como son los camiones y las volquetas.
En la Tabla 7 se muestran las estimaciones de emisión anuales de BC por diferentes inventarios globales y regionales para el año 2000, los cuales muestran que la variación de estas emisiones es considerablemente amplia, donde los rangos a nivel global varían de 4’500.000 a 12’600.000 Ton/año, mientras que los rangos a nivel regional varían de 260.000 a 1’800.000 Ton/año [15].
Tabla 7 Comparación de estimaciones de emisión anual de BC por diferentes inventarios
Fuente: elaboración propia.
En la Tabla 8 se muestra la comparación de emisiones de BC por cada 10000 habitantes para el Valle de Aburrá y los países de China e India, en donde se muestra que China tiene la mayor emisión de BC, la cual es aproximadamente cinco veces la emisión del Valle de Aburrá por cada 10000 habitantes.
Tabla 8 Comparación de estimaciones de emisión anual de BC por cada 10000 habitantes por diferentes inventarios
Fuente: elaboración propia.
Esto se debe al uso desmesurado del carbón como fuente de energía a causa de la alta industrialización, el transporte y el alto desarrollo económico de ese país.
3.5 Dispersión de BC en el Valle de Aburrá
En la Figura 2 se muestra el perfil horario de las concentraciones de BC para el día 9 de noviembre de 2016, en donde se puede observar un pico máximo de concentración de 8.47 µg/m3 a las 7 horas. En la Figura 3, se muestra la dispersión de BC resultado de la modelización en CAMx, para el día de estudio a las 7 horas, donde, de acuerdo con los resultados, se presentan las mayores concentraciones. En la dispersión de BC se puede evidenciar que el área de mayor impacto se encuentra localizada en el centro de Medellín y el norte de este municipio, la cual sigue la dirección de viento en esta zona de estudio.
3.6 Evaluación y verificación del modelo BRAMS
La evaluación del desempeño del modelo se realizó con la comparación de las variables temperatura y humedad relativa del aire como parámetros de comparación entre los datos simulados y los medidos para el día 9 de noviembre de 2016. Posteriormente, se realizó la verificación del modelo comparando dichos parámetros para el día 16 de noviembre de 2016.
Los dos días son evaluados en el punto donde se encuentra la estación automática de monitoreo de la calidad del aire ubicada dentro de la Universidad Pontificia Bolivariana. La evaluación y verificación del modelo meteorológico mostró que este representó la variabilidad diaria de la temperatura y humedad relativa en esta zona del Valle de Aburrá.
Como se observa en la Figura 4, el modelo BRAMS se ajusta a los datos reales, mostrando una tendencia similar para el día en el que se evaluó el modelo y para el día en que se realizó la verificación del mismo.
Fuente: elaboración propia.
Figura 4 Comparación del perfil horario de temperatura y humedad relativa, simulada y medida, para los días evaluados
El Error Porcentual Absoluto Medio (MAPE) de la temperatura y humedad relativa en los días evaluados, arrojó un valor de 0.8 % y 2.6 %, respectivamente, para el 9 de noviembre y 7.7 % y 3.4 %, respectivamente, para el 16 de noviembre de 2016, los cuales se consideran bajos (por debajo de 15 %) y, por lo tanto, de acuerdo con estos resultados, se corrobora que hay un buen ajuste del modelo BRAMS a las condiciones reales medidas en la estación automática de la Universidad Pontificia Bolivariana.
3.7 Evaluación del modelo CAMx
En primera instancia, en la Figura 5, se presenta un análisis cualitativo de la precisión del modelo CAMx al predecir las concentraciones horarias de BC para el 9 de noviembre de 2016. La Figura 5 indica que el modelo se encuentra sobre-estimando aproximadamente el 25 % de las concentraciones simuladas, mientras que el 67 % de éstas se encontraron en la zona 1:1, es decir, los valores se ajustaron muy bien a la realidad. El 8 % restante, muestra que el modelo subestimó las concentraciones, entregando valores menores a los medidos.
Fuente: elaboración propia.
Figura 5 Comparación de la concentración de BC, simulada y medida, para el 9 de noviembre de 2016
La inferencia difusa es una herramienta automática para la evaluación de los resultados simulados de calidad del aire [16], en la cual se obtiene un único indicador de confiabilidad o correlación (desempeño) para los datos, basándose en múltiples mediciones estadísticas indicadoras de confiabilidad/correlación.
El desempeño del modelo para predecir las concentraciones de BC en la estación automática U.P.B. se encuentra en el 66 % para el día de estudio.
A continuación, se presenta la comparación de los perfiles horarios, simulado y observado, de las concentraciones de BC para el día de estudio, el cual tiene condiciones meteorológicas favorables para la dispersión de contaminantes y se encuentra dentro de la campaña de medición realizada por el grupo GHYGAM.
Como se muestra en la Figura 6, las concentraciones de BC simuladas se encuentran muy cerca de las medidas observando una tendencia muy similar en el perfil horario, subestimando ligeramente los datos en la mayoría del día analizado.
Fuente: elaboración propia.
Figura 6 Comparación del perfil horario de la concentración de BC, simulada y medida, para el 9 de noviembre de 2016
Adicionalmente, se observan dos máximos notables, uno en la mañana y otro en la tarde, pero de menor intensidad entre las 5 y las 10 horas y entre las 16 y las 19 horas, respectivamente, obteniendo un pico máximo de concentración BC observado de 12.09 µg/m3 mientras que el pico máximo de concentración BC simulada fue de 8.47 µg/m3. El Error Porcentual Absoluto Medio (MAPE) del día evaluado, arroja un valor de 12.7 % por lo que se considera que hay un ajuste aceptable del modelo CAMx a las condiciones de BC medidas.
La eficiencia del modelo CAMx reportada en la literatura se ve influenciada de forma significativa por varios factores como el componente a modelizar, la meteorología y la locación, mostrando, en la mayoría de estudios, una subestimación con respecto a las concentraciones observadas [17].
En este estudio, las diferencias entre lo medido y lo observado en el modelo CAMx, se deben en gran parte a que en la malla de simulación no se encuentran aún representadas la totalidad de vías importantes en la zona de estudio (teniendo en cuenta que el 97 % de las emisiones de BC se deben a las fuentes móviles), lo cual es una oportunidad de mejora para este Sistema de Modelización Atmosférica. Además, los errores en las suposiciones del modelo, su parametrización, inexactitud en los datos de entrada e incertidumbres relacionadas con la naturaleza estocástica de los procesos atmosféricos, son aspectos que se deben tener en cuenta para minimizar dichas diferencias.
4. CONCLUSIONES
A partir de la estimación de emisiones de BC para el Valle de Aburrá, se concluye que la categoría camiones seguido de la categoría volquetas son los mayores emisores del contaminante BC aportando el 41 % y 38 %, respectivamente del total de las emisiones pertenecientes a fuentes móviles.
El combustible diésel es el principal responsable por las emisiones de BC seguido de la gasolina, con el 96 % y 4 %, respectivamente.
El sector textil es el que más aporta en emisión de este contaminante seguido de la industria maderera, con un 59 % y 30 %, respectivamente, del total de las emisiones de BC asociadas a fuentes fijas.
El municipio de Barbosa realiza los mayores aportes de BC seguido del municipio de Itagüí, con un 30 % y 27 %, respectivamente, del total de las emisiones en el Valle de Aburrá pertenecientes a fuentes fijas.
Para el día de estudio (9 de noviembre de 2016), se obtuvo un comportamiento del perfil horario de las concentraciones de BC con un pico máximo de concentración para el día de verano de 8,47 µg/m3 a las 7 horas.
La dispersión de BC para el día de estudio, muestra que el área de mayor impacto se encuentra localizada en el centro de Medellín y el norte de este municipio.
El modelo meteorológico BRAMS se ajusta a los datos reales de las variables meteorológicas estudiadas (temperatura y humedad relativa), mostrando una tendencia similar para el día en el que se evaluó el modelo y para el día en que se realizó la verificación del mismo, con un error MAPE de la variable temperatura de 0.8 % y 7.7 %, respectivamente, mientras que para la humedad relativa fue de 2.6 % y 3.4 %, respectivamente.
El desempeño del modelo CAMx (inferencia difusa) para predecir las concentraciones de BC, en la estación automática U.P.B. para el día de estudio, se encuentra en el 66 %. A su vez, el modelo se encuentra sobre-estimando aproximadamente el 25 % de las concentraciones simuladas, mientras que el 67 % de éstas se encontraron en la zona 1:1, es decir, los valores se ajustaron muy bien a la realidad. El 8 % restante, muestra que el modelo subestimó las concentraciones entregando valores menores a los medidos.
El día 9 de noviembre de 2016, las concentraciones de BC simuladas se encuentran cerca de las medidas observando una tendencia muy similar en el perfil horario, subestimando ligeramente los datos en la mayoría del día analizado, en el cual el error MAPE es del 12.7 %, el cual es considerado bajo (por debajo de 15 %), por lo que se considera que hay un ajuste aceptable del modelo CAMx a las condiciones de BC medidas.
En este estudio, las diferencias entre lo medido y lo observado en el modelo CAMx se deben en gran parte a que en la malla de simulación no se encuentran aún representadas la totalidad de vías importantes en la zona de estudio (teniendo en cuenta que el 97 % de las emisiones de BC se deben a las fuentes móviles), lo cual es una oportunidad de mejora para este Sistema de Modelización Atmosférica. Además, los errores en las suposiciones del modelo, su parametrización, inexactitud en los datos de entrada e incertidumbres relacionadas con la naturaleza estocástica de los procesos atmosféricos, son aspectos que se deben tener en cuenta para minimizar dichas diferencias adecuadas concordantes con los resultados del trabajo.
