INTRODUCCIÓN
Las partículas atmosféricas (PM, del inglés particulate matter) son una mezcla compleja de pequeñas partículas sólidas y líquidas que se encuentran suspendidas en el aire (Meszaros 1999). Las más preocupantes son las partículas menores de 2.5 µm (PM2.5) y menores de 10 µm (PM10) cuya concentración y composición química varían según la ubicación geográfica y condiciones meteorológicas (Pöschl 2005).
Los principales efectos de las PM están asociados al clima y la salud pública. En el clima, causan un efecto directo de enfriamiento del planeta debido a su capacidad para reflejar la radiación solar. Además tienen un efecto indirecto debido a que promueven la formación de los núcleos de condensación de las nubes, incrementando el reflejo de la radiación solar (Andreae y Crutzen 1997, Pöschl 2005, Querol 2008, Olivier y David 2013). En la salud pública, los efectos de la exposición a las PM están fuertemente asociados con el tamaño, la composición y la concentración de las partículas (Davidson et al. 2005), pudiendo provocar enfermedades respiratorias y cardiovasculares (Hamanaka y Mutlu 2018, Hime et al. 2018). La Organización Mundial de la Salud (OMS) estima que 9 de cada 10 personas respiran aire contaminado (OMS 2018) siendo las partículas atmosféricas las responsables del 7.6 % de las muertes mundiales (Cohen et al. 2017). En áreas urbanas, destacan las emisiones vehiculares (Pant y Harrison 2013) siendo responsables, por ejemplo, de 200 000 muertes prematuras en los Estados Unidos de América cada año (Caiazzo et al. 2013); asimismo, en Lima (Perú), se ha mostrado la relación de las emisiones del transporte con la frecuencia de asma en niños (Carbajal-Arroyo et al. 2007, Tapia et al. 2018).
Para el diseño del control y la mitigación de la exposición a PM, es importante la identificación y cuantificación del aporte de las fuentes de emisión; para ello, en la identificación se usan trazadores elementales (Chavent et al. 2009) como C/Fe/Ba/Zn/Cu, indicadores de fuentes vehiculares, Al/Si/Ca/Fe del suelo o la corteza terrestre, Na/Cl/Mg del océano y V/Ni/, de emisiones industriales o quemas de combustible, que también están asociados con aerosoles secundarios (Viana et al. 2008). Los elementos S/K/Cl/Ca/Br/Zn son considerados trazadores de quema de biomasa (Echalar et al. 1998, Gilardoni et al. 2011)cobrando importancia en la región tropical que posee una alta frecuencia de quemas e incendios forestales (Reddington et al. 2019).
Uno de los modelos receptores exploratorios recomendados para estudios de calidad del aire urbano (Belis et al. 2014) es el análisis de componentes principales absoluto (ACPA), el cual calcula las puntuaciones absolutas (aportes) de los componentes (fuentes) a través de la regresión de las masas diarias y las concentraciones químicas elementales (Thurston y Spengler 1985). En América del Sur, varios autores, como Artaxo et al. (1999), Echalar et al. (1998), Gidhagen et al. (2002) y Goix et al. (2011) han empleado este modelo con resultados favorables.
En la región tropical, sobre los Andes centrales del Perú, se ubica el Valle del Mantaro (VM), cuyo nombre proviene del principal río que recorre la llanura, dividiéndola en dos márgenes. En este valle se ubican ciudades de importancia económica para la Región Junín y la capital del Perú (Lima) debido a que posee una gran productividad agropecuaria y comercial. Por el incremento poblacional y la diversificación de actividades de origen antrópico, la calidad de aire se ha venido degradando (CONAM 2005, Suárez-Salas et al. 2017) a tal punto que el Ministerio de Ambiente de Perú (MINAM) tiene incluida a la capital comercial del valle (Huancayo) dentro de las 31 zonas con atención prioritaria para la implementación de planes de acción de mejora de la calidad de aire (MINAM 2010). Por ello, se han desarrollado una serie de estudios (Suárez-Salas et al. 2017, Suárez- Salas et al. 2017, Estevan et al. 2019 y Huamán De La Cruz et al. 2019) cuya información necesita ser complementada con la determinación y cuantificación de las fuentes de emisión de PM10 y PM2.5 en el VM. Este estudio procesa los datos obtenidos en los años 2007-2008 para brindar información técnica para el desarrollo e implementación de políticas, planes y programas en favor de la mejora de la calidad de aire y protección de la salud de la población.
MATERIALES Y MÉTODOS
Ubicación de las estaciones de monitoreo
La presente investigación se enfocó en las áreas urbanas del VM, que consiste en las ciudades de Jauja, Concepción y Huancayo, cuya población total es de aproximadamente 600 000 habitantes (INEI 2007). La figura 1 muestra las ubicaciones de las estaciones de monitoreo de PM10 y PM2.5, instaladas en Jauja (JA, 11.77 S 75.49 W, 3475 msnm), Concepción (CON, 11.9 S 75.31 W, 3300 msnm) y Huancayo (HYO, 12.07 S, 75.21 W, 3272 msnm), las cuales se encuentran ubicadas sobre la llanura del VM rodeada de montañas desde los 3200 hasta los 4100 msnm.
En el VM la actividad económica principal es la comercial seguida del transporte público. Rodeando a las urbes, existen grandes extensiones de terreno con uso agrícola y ganadero cuyas poblaciones se encuentran conectadas por la carretera central, principal vía de la zona.
Climatología
La figura 2 ilustra los climatogramas obtenidos de las estaciones de Jauja (11.74S 75.47 W) e Ingenio (11.88 S 75.26 W) del Servicio Nacional de Meteorología e Hidrología (SENAMHI) y la estación del Huayao (12.04 S, 75.32 W) del Instituto Geofísico del Perú (IGP) desde 1966 hasta 2002 (IGP 2005). Estas estaciones son las más cercanas a los sitios de muestreo.
En las tres estaciones la variación de la temperatura media describe dos períodos definidos durante el año, el primero de valores mínimos, centrados en junio y julio (invierno), y el segundo de valores máximos entre octubre y diciembre (primavera-verano). La temperatura media anual para Jauja es de 11.9 ± 0.9 ºC, para Ingenio de 11 ± 1 ºC y para Huayao de 11.9 ± 1.2 ºC.
Del mismo modo, la precipitación en las tres estaciones registra dos períodos bien definidos, el primero corresponde a la época seca (junio-septiembre), y el segundo a la época húmeda (octubre-mayo), siendo enero (para la estaciones JA y CON) y febrero (para HYO) los meses de máxima precipitación. El promedio anual de precipitación para JA fue de 704 ± 47.2 mm, para CON de 755 ± 50.2 mm y para HYO de 752 ± 44.3 mm.
Muestreo de partículas atmosféricas
La campaña de monitoreo de PM10 y PM2.5 se llevó a cabo con un muestreador de bajo volumen (Partisol FRM 2000 de Thermo Scientific, EUA) a 16.7 L/min por un tiempo de 24 horas. Se empleó un ciclón para la colecta de PM2.5 y un cabezal de tipo inlet para PM10 (USEPA). En cada estación, el equipo se instaló entre los 8 y 10 m del suelo, cerca de las intersecciones de calles de gran confluencia vehicular y comercial; complementariamente se instaló una estación meteorológica (modelo 232 de Weatherhawk, EUA).
La campaña de monitoreo se inició el 31 de julio de 2007 y terminó el 17 de octubre de 2008; durante este tiempo se realizaron colectas de PM10, con una frecuencia de monitoreo de cada tres días, y colectas de PM2.5, con una frecuencia de monitoreo de cada 15 días. Ambos tamaños de partículas fueron colectados en filtros de teflón con diámetro de 47 mm y tamaño de poro de 2.0 µm.
La colecta de PM10 y PM2.5 se desarrolló en forma secuencial, por ello, los equipos de monitoreo se trasladaron de una estación a otra, cada mes. El detalle de los meses de monitoreo y la cantidad de filtros colectados en cada estación del VM se describen en el cuadro I. Para el aseguramiento de la calidad del proceso de monitoreo se emplearon filtros testigo (10 % de toda la campaña), cuyos resultados en su mayoría fueron menores al límite de detección, a excepción de los elementos: Ca, Cr, Ge, Fe, V, Co, Cu y Ni, Ti y Zn para filtros testigo de PM10 y Co, Cu, Ge, Fe, Mn, Ni y Zn en filtros testigo de PM2.5. Los valores de los elementos detectados en los filtros en testigo fueron restados a los resultados de las muestras.
Estación | Descripción | Meses de monitoreo | Número muestras PM10 | Número muestras PM2.5 |
Jauja (JA) | Ubicada en la azotea de un hospedaje a 10 m de altura, a 5 cuadras de la plaza principal. | septiembre, noviembre (2007). febrero, junio (2008). | 41 | 13 |
Concepción (CON) | Ubicada en el auditorio municipal a 8 m de altura, al frente de la plaza principal. | octubre, diciembre (2007). marzo, julio, agosto, septiembre, octubre (2008). | 47 | 11 |
Huancayo (HYO) | Ubicada en el Colegio Rosario a 8 m de altura, a 2 cuadras de la plaza principal. | agosto (2007). enero, abril, mayo (2008) | 41 | 9 |
Análisis químico de las muestras
Los filtros obtenidos en toda la campaña fueron analizados por el Laboratorio Chester LabNet en EUA, certificado por el Programa Nacional de Acreditación de Laboratorios Ambientales (NELAC, por sus siglas en inglés). La masa se determinó por el método gravimétrico (EPA IO-3.1) mediante una microbalanza de ± 0.01 mg de precisión. Las concentraciones de los elementos químicos se determinaron mediante la técnica de fluorescencia de rayos X (XRF) establecida en el método analítico EPA IO-3.3. Este método consiste en bombardear con rayos X a los filtros con muestra, para que, con base en su excitación y liberación de energía, se determine su concentración según la huella digital de cada elemento. Los límites de detección de los 38 elementos analizados por XRF se muestran en el cuadro II.
Elemento | µg/m3 | Elemento | µg/m3 | Elemento | µg/m3 |
Al | 0.008 | Sr | 0.0022 | Ag | 0.0075 |
Sb | 0.011 | P | 0.0047 | Pb | 0.0039 |
As | 0.002 | Ga | 0.0021 | K | 0.0030 |
S | 0.004 | Ge | 0.0013 | Rb | 0.0014 |
Ba | 0.083 | Fe | 0.0008 | Se | 0.0014 |
Br | 0.001 | In | 0.0196 | Si | 0.0055 |
Cd | 0.008 | La | 0.0508 | Na | 0.1484 |
Ca | 0.002 | Mg | 0.0178 | Ti | 0.0014 |
Cl | 0.005 | Mn | 0.0015 | V | 0.0010 |
Co | 0.001 | Hg | 0.0035 | Y | 0.0021 |
Cu | 0.001 | Mo | 0.0035 | Zn | 0.0010 |
Cr | 0.001 | Ni | 0.0008 | Zr | 0.0025 |
Sn | 0.013 | Pd | 0.0071 |
Análisis estadístico
Los resultados del laboratorio (µg/filtro) fueron convertidos a µg/m3 por el método EPA IO-2.4, tomando en cuenta el flujo de aire que circuló durante las 24 h de muestreo, la temperatura y la presión atmosférica de cada estación de monitoreo.
Se calcularon los estadísticos descriptivos y pruebas de T, análisis de varianza y Kruskal-Wallis con nivel de significancia del 95 % (p < 0.05). Se compararon las concentraciones de PM10 y PM2.5 con los Estándares de Calidad Ambiental (ECA) establecidos por la normatividad peruana (MINAM 2017).
Factor de enriquecimiento
Para diferenciar a los elementos químicos procedentes de actividades de origen en fuentes naturales y antrópicas, se aplicó el factor de enriquecimiento (FE). Para ello se usó la concentración de referencia del aluminio (Al) de Wedepohl K H. (1995), bajo la siguiente ecuación (Rahn 1971):
Donde X es la concentración de cada elemento obtenida en cada muestra, Al es la concentración de referencia del Al. Cuando el FE de una muestra se aproxima a la unidad, la fuente dominante para ese elemento es la corteza terrestre; Si FE > 10, una fracción significativa del elemento es de una fuente de origen antrópico (Chester et al. 2000).
Análisis multivariado
Se aplicó el análisis de componentes principales (ACP) con método varimax a las concentraciones elementales estandarizadas a media cero y varianza unidad (Pison et al. 2003). Para el caso de las PM10 de cada estación, se consideró a los elementos químicos cuyos datos fueron mayores a 30 (Hopke 2010); para las PM2.5, dado el menor número de muestras, se consideró a los elementos químicos trazadores, cuidando que el número de datos fuera mayor al número de elementos (Husson et al. 2010). El número de componentes (factores) considerados para el análisis e interpretación, en cada estación, fueron los factores con autovalores mayores a 1 (Kaiser 1958). Los elementos químicos con cargas factoriales mayores a 0.5 fueron considerados influyentes para cada componente.
Para confirmar los factores seleccionados, se utilizó el análisis de conglomerado jerárquico (ACJ) con método Ward y distancia euclidiana al cuadrado para identificar similitudes y características de las asociaciones entre los elementos químicos. El número de grupos se realizó conforme al número de factores determinados.
Dado que el ACP proporciona un reparto de los elementos químicos por factores, cada factor fue considerado como una fuente de emisión de PM. Para estimar la contribución de cada fuente se realizó el análisis de componentes principales absoluto (ACPA) según la metodología de Thurston y Spengler (1985).
En todos los análisis realizados se empleó el programa libre R versión 3.5.1 (R Core Team 2015).
RESULTADOS Y DISCUSIÓN
Dirección y velocidad de los vientos
Como se observa en la figura 3, en los meses monitoreados de la estación de JA, los vientos predominantes provienen del nornoroeste y norte, con velocidades dominantes comprendidas entre 7 a 31.3 m/s; asimismo, en la estación de CON, los vientos predominantes provinieron del noreste con velocidades comprendidas entre 1 a 3 m/s; en el caso de la estación de HYO, los vientos predominantes provinieron del sursuroeste y sur con velocidades dominantes de 1 a 3 m/s. En las tres estaciones, se observan vientos provenientes de otras direcciones, lo cual indica la existencia de mezcla de masas de aire (Flores-Rojas et al. 2019).
Partículas atmosféricas
PM 10
En el cuadro III, se muestran las concentraciones promedio ± desviación estándar (prom ± ds), valor máximo (máx.) y cantidad de datos (N) de las PM10 y los elementos químicos de las estaciones JA, CON y HYO. La concentración de PM10 en el centro urbano de Huancayo (65 ±31 µg/m3) fue significativamente mayor (F = 44.25; p < 0.0001) con respecto a las demás estaciones. Por lo tanto, el área urbana de Huancayo posee la mayor problemática de partículas del VM. Huancayo es considerada la capital comercial y económica de la región Junín y alberga el 84 % de la población del VM (aproximadamente 500 000 habitantes). No se encontró diferencia significativa (T = 0.09; p = 0.93) entre las PM10 de JA (32 ± 12 µg/m3) y CON (28 ± 13 µg/m3), lo cual se debe a que son ciudades más pequeñas.
JA | CON | HYO | |||||||||
Prom ± de. | Máx | N | Prom ± de | Máx | N | Prom ± de | máx | N | |||
PM10 | 31.7 ± 11.6 | 52.9 | 41 | PM10 | 28.2 ± 13.1 | 71.3 | 47 | PM10 | 64.5 ± 30.9 | 143 | 41 |
Si | 1.8 ± 0.7 | 3.2 | 41 | Si | 1.8 ± 0.9 | 3.7 | 47 | Si | 4.8 ± 2.4 | 11 | 40 |
Ca | 1.3 ± 0.5 | 2.2 | 41 | Al | 0.9 ± 0.6 | 3.9 | 45 | Ca | 2.4 ± 1.2 | 5.8 | 41 |
S | 1 ± 0.8 | 3 | 41 | Ca | 0.7 ± 0.3 | 1.5 | 47 | Al | 2 ± 1.1 | 4.7 | 40 |
Al | 0.6 ± 0.3 | 1.2 | 41 | S | 0.7 ± 0.6 | 2.6 | 45 | Fe | 1.4 ± 0.7 | 3.2 | 40 |
Fe | 0.6 ± 0.2 | 1 | 41 | K | 0.6 ± 0.8 | 5.4 | 46 | K | 1.1 ± 0.6 | 3.1 | 40 |
K | 0.5 ± 0.2 | 0.8 | 41 | Fe | 0.6 ± 0.2 | 1.1 | 45 | S | 0.9 ± 0.6 | 2.2 | 40 |
Na | 0.2 ± 0.05 | 0.3 | 6 | Cl | 0.3 ± 1.4 | 9.5 | 45 | Cl | 0.3 ± 0.29 | 1.3 | 40 |
Pb | 0.2 ± 0.2 | 0.8 | 40 | Na | 0.3 ± 0.03 | 0.3 | 2 | Na | 0.2 ± 0.1 | 0.4 | 9 |
As | 0.1 ± 0.2 | 0.6 | 40 | Mg | 0.1 ± 0.05 | 0.2 | 21 | Mg | 0.2 ± 0.1 | 0.5 | 38 |
Cl | 0.1 ± 0.1 | 0.3 | 37 | Ti | 0.1 ± 0.02 | 0.1 | 46 | Ti | 0.1 ± 0.1 | 0.3 | 40 |
Sb | 0.1 ± 0.1 | 0.4 | 18 | Zn | 0.03 ± 0.01 | 0.1 | 47 | Zn | 0.1 ± 0.03 | 0.1 | 40 |
Zn | 0.1 ± 0.1 | 0.2 | 41 | Pb | 0.02 ± 0.02 | 0.09 | 43 | Mn | 0.04 ± 0.02 | 0.1 | 40 |
Ti | 0.05 ± 0.02 | 0.08 | 41 | Sb | 0.02 ± 0.004 | 0.02 | 5 | Pb | 0.03 ± 0.03 | 0.19 | 40 |
Mg | 0.04 ± 0.01 | 0.06 | 14 | Mn | 0.02 ± 0.01 | 0.03 | 45 | As | 0.02 ± 0.03 | 0.12 | 38 |
Cu | 0.02 ± 0.02 | 0.08 | 40 | As | 0.01 ± 0.01 | 0.08 | 42 | Sb | 0.02 ± 0.005 | 0.03 | 9 |
Mn | 0.01 ± 0.01 | 0.03 | 41 | P | 0.01 ± 0.004 | 0.02 | 17 | P | 0.02 ± 0.01 | 0.05 | 25 |
Cd | 0.01 ± 0.003 | 0.02 | 7 | Sr | 0.01 ± 0.02 | 0.1 | 31 | Cu | 0.01 ± 0.01 | 0.02 | 40 |
P | 0.01 ± 0.01 | 0.02 | 5 | Br | 0.004 ± 0.008 | 0.047 | 32 | Cd | 0.01 ± 0.0004 | 0.01 | 3 |
Ag | 0.01 ± 0.001 | 0.01 | 2 | Cu | 0.004 ± 0.002 | 0.011 | 45 | Sr | 0.01 ± 0.004 | 0.02 | 39 |
Se | 0.01 ± 0.002 | 0.01 | 12 | Zr | 0.003 ± 0.001 | 0.005 | 7 | Zr | 0.005 ± 0.002 | 0.01 | 26 |
Sr | 0.004 ± 0.001 | 0.007 | 26 | Rb | 0.002 ± 0.001 | 0.005 | 39 | V | 0.005 ± 0.003 | 0.012 | 16 |
Zr | 0.004 ± 0.001 | 0.005 | 6 | V | 0.002 ± 0.001 | 0.004 | 14 | Rb | 0.004 ± 0.002 | 0.009 | 39 |
Rb | 0.003 ± 0.001 | 0.003 | 19 | Cr | 0.002 ± 0.0004 | 0.002 | 9 | Br | 0.003 ± 0.002 | 0.01 | 35 |
Br | 0.002 ± 0.001 | 0.005 | 14 | Ge | 0.001 ± 0.0004 | 0.001 | 3 | Cr | 0.003 ± 0.002 | 0.008 | 32 |
V | 0.002 ± 0.001 | 0.004 | 10 | Ni | 0.001 ± 0.001 | 0.002 | 22 | Ni | 0.002 ± 0.002 | 0.008 | 15 |
Cr | 0.002 ± 0.001 | 0.003 | 11 | Co | 0.0003 ± 0.0004 | 0.0007 | 3 | Se | 0.002 ± 0.0001 | 0.002 | 2 |
Ni | 0.001 ± 0.001 | 0.004 | 18 |
Prom = promedio aritmético, de = desviación estándar, máx = valor máximo, N = número de datos, JA = Jauja, CON = Concepción, HYO = Huancayo.
Los tres elementos químicos con mayor concentración registrada en las PM10 en las tres estaciones de monitoreo fueron Si, Ca y Al, trazadores de fuentes geológicas o suelo. Este resultado es razonable por la mezcla de vientos que transporta materiales de áreas rurales y calles que no están pavimentados, por otro lado, en el tiempo de monitoreo se observaron diversas construcciones de edificios.
En la figura 4, se observa la variabilidad diaria de las PM10 en las estaciones de JA, CON y HYO, adicionalmente se muestra en línea roja horizontal el ECA peruano (100 µg/m3), el cual es el doble del valor recomendado por la OMS (50 µg/m3) (OMS 2006). Sin embargo, en tres días de monitoreo en el centro de HYO fue superado el ECA, correspondiendo estos valores a los días 09/08/2007, 24/08/2007 y 23/05/2008. Estos valores altos son un indicador inicial de que los valores de PM representarían un riesgo para la población.
PM 2.5
Al igual que en el caso de las PM10, la concentración de las partículas finas del centro urbano de HYO (34 ± 15 µg/m3) fue mayor a lo registrado en las demás estaciones (Cuadro IV); igualmente, el promedio de PM2.5 de JA (22 ± 14 µg/m3) fue superior al de CON (16 ± 8 µg/m3) existiendo una diferencia significativa entre las concentraciones registradas en las estaciones (Kruskal-Wallis = 8.6; p = 0.01). Los elementos químicos más abundantes en las PM2.5 de las tres estaciones fueron el S, Si y K. Asimismo, se observa mayor detección de elementos químicos en JA (21 elementos químicos) comparado con CON (18 elementos químicos) y HYO (16 elementos).
JA | CON | HYO | |||||||||
Prom ± de | Máx | N | Prom ± de | Máx | N | Prom ± de | Máx | N | |||
PM2.5 | 21.8 ± 13.8 | 51.6 | 13 | PM2.5 | 15.8 ± 7.8 | 26.6 | 11 | PM2.5 | 34.5 ± 14.7 | 67.8 | 9 |
S | 1.2 ± 0.8 | 2.3 | 13 | S | 0.9 ± 0.9 | 2.4 | 11 | S | 0.9 ± 0.7 | 1.9 | 9 |
Si | 0.4 ± 0.3 | 0.9 | 13 | Si | 0.3 ± 0.2 | 0.5 | 11 | Si | 0.7 ± 0.3 | 1.1 | 9 |
K | 0.3 ± 0.2 | 0.8 | 13 | K | 0.2 ± 0.1 | 0.3 | 11 | K | 0.4 ± 0.3 | 1 | 9 |
Ca | 0.2 ± 0.1 | 0.4 | 13 | Al | 0.2 ± 0.1 | 0.3 | 11 | Ca | 0.3 ± 0.1 | 0.5 | 9 |
Na | 0.2 ± 0.05 | 0.2 | 2 | Ca | 0.1 ± 0.1 | 0.3 | 11 | Cl | 0.3 ± 0.2 | 0.5 | 8 |
Al | 0.2 ± 0.1 | 0.4 | 13 | Fe | 0.1 ± 0.04 | 0.2 | 11 | Al | 0.3 ± 0.2 | 0.6 | 9 |
Cl | 0.1 ± 0.1 | 0.3 | 6 | Pb | 0.1 ± 0.05 | 0.1 | 6 | Fe | 0.2 ± 0.1 | 0.3 | 9 |
Pb | 0.1 ± 0.2 | 0.7 | 12 | Cl | 0.1 ± 0.1 | 0.2 | 7 | Mg | 0.04 ± 0.01 | 0.05 | 4 |
Fe | 0.1 ± 0.1 | 0.3 | 13 | Sb | 0.02 ± 0.01 | 0.03 | 2 | Zn | 0.02 ± 0.01 | 0.03 | 9 |
As | 0.1 ± 0.1 | 0.3 | 12 | As | 0.02 ± 0.02 | 0.07 | 11 | Pb | 0.02 ± 0.01 | 0.04 | 9 |
Sb | 0.04 ± 0.03 | 0.1 | 5 | Zn | 0.02 ± 0.01 | 0.04 | 11 | Ti | 0.02 ± 0.01 | 0.03 | 9 |
Mg | 0.04 ± 0.02 | 0.1 | 3 | Ti | 0.01 ± 0.01 | 0.02 | 11 | As | 0.01 ± 0.01 | 0.03 | 9 |
Zn | 0.02 ± 0.03 | 0.1 | 13 | Se | 0.004 ± 0.001 | 0.005 | 3 | Mn | 0.01 ± 0.002 | 0.01 | 9 |
Cd | 0.01 ± 0.005 | 0.02 | 2 | Mn | 0.004 ± 0.001 | 0.01 | 9 | Cu | 0.004 ± 0.005 | 0.01 | 8 |
Ti | 0.01 ± 0.01 | 0.02 | 13 | Cu | 0.003 ± 0.002 | 0.01 | 8 | Br | 0.003 ± 0.001 | 0.01 | 8 |
Cu | 0.01 ± 0.01 | 0.02 | 11 | Ni | 0.003 ± 0.001 | 0.004 | 4 | Ni | 0.003 ± 0.002 | 0.004 | 2 |
Se | 0.01 ± 0.003 | 0.01 | 4 | Rb | 0.002 ± 0.0003 | 0.002 | 2 | ||||
Mn | 0.003 ± 0.001 | 0.01 | 8 | Br | 0.002 ± 0.0004 | 0.002 | 7 | ||||
Br | 0.003 ± 0.002 | 0.01 | 6 | ||||||||
Rb | 0.002 ± 0.0003 | 0.002 | 2 | ||||||||
Ni | 0.001 ± 0.0003 | 0.002 | 4 |
Prom = promedio aritmético, de = desviación estándar, máx = valor máximo; N = Número de datos, JA = Jauja, CON = Concepción, HYO = Huancayo.
En la figura 5, se observa la variabilidad diaria de PM2.5 en las estaciones de JA, CON y HYO durante toda la campaña de monitoreo. Comparando los datos, con el ECA peruano de 24 horas (50 µg/m3), se comprueba que la norma se rebasó en dos ocasiones: una en HYO (67.8 µg/m3, 25/08/2007) y la otra en JA (51.7 µg/m3, 02/10/2007), las cuales han sido registradas en temporada seca, época donde aumenta la concentración de PM10 y PM2.5 (Suárez et al. 2017 y Huamán de la Cruz et al. 2019) y el reporte de quemas de vegetación que afectan al VM (Moya-Álvarez et al. 2017).
Factor de enriquecimiento (FE)
PM 10
En la figura 6 se observan los factores de enriquecimiento determinados para los elementos químicos de PM10 en las estaciones evaluadas. Para JA, el Si, Al, Ti y Fe mostraron valores de FE menores que 1, lo que sugiere estos elementos provienen de fuentes naturales; en cambio, los valores de FE del As, Pb, S, Zn, Cu y Cl fueron mayores que 10, indicándonos su origen antrópico. Para el caso de CON, el Si, Al, Fe y Ti estarían proviniendo de fuentes naturales y el As, Pb, S, Cl, Zn y Cu de fuentes antrópicas. Para la estación de HYO, sólo el Mg y el Si estarían proviniendo de fuentes naturales y el As, Br, Pb, S, Zn, Cl, y Cu de fuentes antrópicas. Según el FE se observa que mayor número de elementos químicos provienen de fuentes antrópicas, siendo el As, Pb, Br, S, Cl, Cu y Zn los elementos de mayor factor.
PM 2.5
En la figura 7 se observan los FE determinados para los elementos químicos registrados en las PM2.5. Para JA, el Si, Al, Ti, Fe, Ca y K registraron valores de FE < 1, lo que indica que provienen de fuentes naturales; en cambio, el As, Pb, S y Zn, mostraron FE > 10, siendo su procedencia antrópica. Para el caso de CON, el Si, Al, Fe, Ti, Ca y K tendrían procedencia natural y el As, Br, S y Zn de fuentes antrópicas. Para HYO, el Al y Fe estarían proviniendo de fuentes naturales y el As, Br, S y Pb de fuentes antrópicas.
Por lo tanto, al igual que en las PM10, el As, Pb, Br, Zn y S tendrían procedencia antrópica y el Si, Al, Ti y Fe procedencia natural.
Identificación cuantitativa de las fuentes de emisión
PM 10
En el cuadro V se presentan las cargas factoriales, autovalores y porcentaje de varianza del ACP para las PM10 de las estaciones de JA, CON y HYO, respectivamente.
JA | CON | HYO | ||||||||||||
Elem. | Polvo del suelo | Fundición | Combustible-Aceite | Quema de biomasa | Elem. | Polvo del suelo | Fundición | Quema de biomasa | Vehículos | Elem. | Polvo del suelo | Fundición | Vehículos | Quema de biomasa |
Al | 0.96 | 0.2 | -0.05 | 0.01 | Si | 0.98 | 0.09 | 0.12 | 0.02 | Al | 0.94 | 0.21 | 0.18 | 0.09 |
Ca | 0.95 | 0.21 | 0.07 | 0 | Fe | 0.97 | 0.09 | 0.17 | 0.02 | Ca | 0.94 | 0.15 | 0.12 | 0.17 |
Si | 0.94 | 0.31 | -0.05 | 0.01 | Ca | 0.96 | 0.09 | 0.16 | -0.08 | Fe | 0.94 | 0.23 | 0.2 | 0.11 |
Fe | 0.93 | 0.36 | 0.01 | 0 | Mn | 0.95 | 0.18 | 0.18 | 0.01 | Si | 0.94 | 0.22 | 0.19 | 0.09 |
Ti | 0.93 | 0.32 | -0.01 | -0.04 | Ti | 0.94 | 0.03 | 0.16 | 0.08 | Ti | 0.94 | 0.21 | 0.17 | 0.14 |
Mn | 0.89 | 0.37 | 0.08 | -0.02 | Zn | 0.7 | 0.47 | -0.07 | 0.39 | Mn | 0.93 | 0.27 | 0.14 | 0.15 |
K | 0.87 | 0.3 | 0.02 | 0.19 | Al | 0.57 | 0.15 | 0.79 | 0.07 | K | 0.93 | 0.19 | 0.12 | 0.2 |
Pb | 0.27 | 0.94 | -0.06 | -0.02 | As | 0.13 | 0.95 | -0.07 | -0.09 | Rb | 0.93 | 0 | 0.16 | 0.14 |
Zn | 0.37 | 0.91 | -0.09 | 0.02 | Pb | 0.1 | 0.93 | 0.16 | 0.18 | Sr | 0.9 | 0.13 | 0.11 | 0.17 |
As | 0.32 | 0.91 | 0.05 | -0.04 | S | 0.02 | 0.75 | 0.47 | 0.06 | Zn | 0.87 | 0.37 | 0.16 | 0.21 |
Cu | 0.34 | 0.88 | -0.2 | 0.03 | Cu | 0.35 | 0.55 | 0.26 | 0.4 | Mg | 0.83 | 0.15 | 0.06 | 0.23 |
S | 0.21 | 0.73 | 0.03 | 0.16 | K | 0.16 | 0.17 | 0.95 | 0.16 | Cl | 0.79 | 0.04 | 0.25 | 0.25 |
Ni | 0.05 | -0.11 | 0.98 | -0.11 | Cl | -0.04 | 0.08 | 0.14 | 0.95 | Cu | 0.73 | 0.51 | 0.23 | 0.08 |
Cl | 0.03 | 0.07 | -0.11 | 0.98 | As | 0.02 | 0.94 | 0.16 | 0 | |||||
Pb | 0.37 | 0.87 | 0.11 | 0.01 | ||||||||||
S | 0.16 | 0.85 | -0.12 | -0.01 | ||||||||||
Cr | 0.4 | 0.1 | 0.89 | 0.04 | ||||||||||
Br | 0.47 | -0.04 | 0.04 | 0.87 | ||||||||||
Cargas | 6.46 | 4.51 | 1.05 | 1.05 | 5.58 | 2.98 | 2 | 1.31 | 11.01 | 3.13 | 1.21 | 1.12 | ||
Var. propor. | 0.46 | 0.32 | 0.07 | 0.07 | 0.43 | 0.23 | 0.15 | 0.10 | 0.61 | 0.17 | 0.07 | 0.06 | ||
Var. acumul. | 0.46 | 0.78 | 0.86 | 0.93 | 0.43 | 0.66 | 0.81 | 0.91 | 0.61 | 0.79 | 0.85 | 0.91 |
JA = Jauja, CON = Concepción, HYO = Huancayo, Elem.= elemento.
Para JA, se identificaron cuatro componentes principales que explican el 93 % de la varianza, el primer factor (46 % de la variabilidad) agrupa a los elementos: Al, Ca, Si, Fe, Ti, Mn y K, trazadores de polvo del suelo (Echalar et al. 1998, Artaxo et al. 1999). El segundo factor (32 %) conformado por Pb, Zn, As y Cu, fue referido a fundición (Querol et al. 2001, Querol et al. 2002 y Goix et al. 2011), dado que en la provincia de Jauja no existen industrias ni actividades relacionadas a fundiciones, la fuente más probable se encontraría ubicada a 50 km al noroeste del Valle del Mantaro, en la ciudad de La Oroya donde se ubica la fundición del Complejo Metalúrgico de La Oroya (CMLO), como se muestra en la figura 8. Ahí se fundían y refinaban minerales polimetálicos provenientes de los centros mineros de la zona central del Perú hasta su paralización en el año 2009. Según un estudio de Klepel (2005) existen indicios de que los contaminantes vertidos a la atmosfera por la chimenea del CMLO, estarían llegando al VM. El tercer factor (7 %), conformado por el Ni, es atribuible al combustible-aceite (Viana et al. 2008) y el cuarto factor (7 %), conformado por el Cl, estaría proviniendo de la quema de biomasa (Santanna et al. 2016). Como se mencionó anteriormente, en la periferia de la ciudad de Jauja se desarrollan actividades agrícolas y pecuarias donde la quema de pastizales y rastrojos es una práctica común.
Para CON, se obtuvieron cuatro componentes principales que explican el 91 % de la variación, el primer factor (43 %) conformado por el Si, Fe, Ca, Al, Mn, Ti y Zn se atribuye al polvo del suelo. El segundo factor (23 %) conformado por As, Pb, S y Cu se atribuye a la fundición del CMLO, debido a que también en la provincia de Concepción no existen industrias ni actividades relacionadas a fundiciones. El tercer factor (15 %) conformado por el Al y el K se atribuye a la quema de biomasa (Echalar et al. 1998). En la ciudad de Concepción una práctica común de los pobladores es el uso de leña como combustible para cocinar sus alimentos, así como la quema de pastizales y rastrojos en los trabajos de campo. El cuarto factor (10 %) conformado por el Cl, también se agrupa con el Cu y el Zn (cargas factoriales ≈ 0.4) por la significativa correlación (p < 0.05), por lo que se atribuye a fuentes vehiculares (Stechmaun y Daunecker 1990, Viana et al. 2008), que en la ciudad están conformadas especialmente por motos y vehículos menores, que realizan servicio de transporte público en forma frecuente.
Para el caso de HYO, se obtuvieron cuatro factores principales que explican el 91% de la variabilidad. El primer factor (61 %) agrupa a los elementos Al, Ca, Fe, Si, Ti, Mn, K, Rb, Sr, Zn, Mg, Cl y Cu atribuibles al polvo del suelo, el cual estaría relacionado al parque automotor, por la presencia de Cl, Cu y Zn (trazadores antrópicos). El segundo factor (17 %) está conformado por el As, Pb y S atribuible a emisiones de fundición. En Huancayo existen pequeñas fundiciones de hierro-acero y joyerías. El tercer factor (7 %), se debe al Cr atribuible a vehículos (Srimuruganandam and Shiva-Nagendra 2012a, 2012b, Banerjee et al. 2015), el cual es esperado dado que el sitio de monitoreo se ubicó en una zona de intenso tráfico vehicular. El cuarto factor (6 %) conformado por el Br es atribuible a la quema de biomasa (Echalar et al. 1998, Santanna et al. 2016). Alrededor de la ciudad de Huancayo también se desarrollan actividades agrícolas donde se queman pastizales y rastrojos agrícolas.
En la figura 9 se muestran los porcentajes de emisión para cada fuente de PM10, determinada mediante la aplicación del ACPA. En todas las estaciones el polvo del suelo fue el principal emisor de PM10 con 66 %, 57 % y 48 % para JA, CON y HYO, respectivamente. Para JA (14 %) y HYO (23 %) la quema de biomasa fue el segundo en importancia de emisión. En Concepción (7 %) y Huancayo (10 %) se registraron aportes de fuentes vehiculares, ello debido a que poseen mayor tráfico. Se determinaron indicios de aportes del CMLO en las PM10 registradas en Jauja (4%) y Concepción (5%). En cambio, para Huancayo se registraron aportes de fundiciones locales en un 13 %. Llama la atención que en la ciudad de Concepción el porcentaje no identificado de fuentes (Otros) ocupe el 27 % de aporte, ello quizás se deba a que es una zona de convergencia de vientos para el valle del Mantaro (Flores et al. 2019).
Las fuentes determinadas fueron confirmadas por el agrupamiento de los elementos químicos en los dendrogramas que se muestran en la figura 10. Para el caso de JA se determinaron cuatro grupos: G1, con elementos trazadores de polvo del suelo (K, Al, Si, Fe, Ti, Ca y Mn), G2, con elementos trazadores de fundición (S, As, Cu, Pb y Zn), G3 con un elemento trazador de quema de biomasa (Cl) y G4, con un elemento trazador de combustible-aceite (Ni). En CON, se determinaron cinco grupos: G1, con elementos trazadores de fundición (S, As y Pb), G2, con un elemento trazador de quema de biomasa (Cl), G3, con elementos trazadores de vehículos (Cu y Zn), y G4 y GV con elementos trazadores de suelo (Al y K, y Ti, Ca, Mn, Fe y Si, respectivamente). En HYO se determinaron cuatro grupos: G1, con elementos trazadores de fundición (S, As y Pb), G2, con un elemento trazador de vehículos (Cr), G3, con un elemento trazador de quema de biomasa (Br), y el G4 con elementos trazadores de suelo y vehículos (Mg, Sr, K, Rb, Zn, Ca, Mn, Ti, Al, Fe, Si, Cl y Cu). Lo anterior confirma que existe un mayor número de fuentes antrópicas de PM10 en cada lugar estudiado.
PM 2.5
En el cuadro VI se presentan las cargas factoriales, autovalores y porcentaje de varianza del ACP para PM2.5 de las estaciones de JA, CON y HYO. Para Jauja, se identificaron dos componentes principales que explican el 90 % de la varianza, el primer factor (53 % de la variabilidad) agrupa a los elementos Al, Ca, Fe, Si, Ti y K, trazadores de polvo del suelo. El segundo factor (36 %) conformado por As, Zn, Pb y S, atribuibles a la fundición del CMLO. Para Concepción, también se obtuvieron dos componentes principales que explican el 95 % de la variación, el primer factor (59 % de la variabilidad) conformado por el Al, Si, Fe, Ca y K es atribuible al polvo del suelo; el segundo factor (36 %) conformado por S y As fue atribuido a la fundición del CMLO. Para el caso de Huancayo, se obtuvieron cuatro factores principales que explican el 99 % de la variabilidad. El primer factor (30 % de la variabilidad) agrupa a los elementos Fe, Al y K atribuibles al polvo del suelo. El segundo factor (29 %) está conformado por el As y Pb atribuibles a fundición. El tercer factor (21 %), donde tenemos al Br es atribuible a la quema de biomasa y el cuarto factor (18 %) conformado por el S atribuible a vehículos.
JA | CON | HYO | ||||||||
Elemento | Polvo del suelo | Fundición | Elemento | Polvo del suelo | Fundición | Elemento | Polvo del suelo | Fundición | Quema de biomasa | Vehículos |
Al | 0.99 | -0.14 | Al | 0.95 | 0.27 | Fe | 0.87 | 0.39 | 0.26 | -0.08 |
Ca | 0.99 | 0.03 | Si | 0.93 | 0.35 | Al | 0.81 | 0.36 | 0.45 | -0.08 |
Fe | 0.99 | 0.06 | Fe | 0.9 | 0.29 | K | 0.59 | 0.34 | 0.55 | 0.47 |
Si | 0.99 | -0.02 | Ca | 0.87 | 0.48 | Pb | 0.3 | 0.92 | 0.15 | 0.14 |
Ti | 0.99 | 0.07 | K | 0.73 | 0.57 | As | 0.36 | 0.85 | 0.23 | 0.25 |
K | 0.62 | 0.28 | S | 0.4 | 0.88 | Br | 0.35 | 0.19 | 0.91 | 0.08 |
As | -0.04 | 0.99 | As | 0.28 | 0.95 | S | -0.08 | 0.18 | 0.06 | 0.98 |
Zn | 0.01 | 0.98 | ||||||||
Pb | -0.03 | 0.98 | ||||||||
S | 0.19 | 0.79 | ||||||||
Cargas | 5.32 | 3.64 | 4.11 | 2.52 | 2.11 | 2.04 | 1.48 | 1.27 | ||
Var. proporción | 0.53 | 0.36 | 0.59 | 0.36 | 0.3 | 0.29 | 0.21 | 0.18 | ||
Var. acumulada | 0.53 | 0.90 | 0.59 | 0.95 | 0.3 | 0.59 | 0.80 | 0.99 |
JA = Jauja, CON = Concepción, HYO = Huancayo.
En la figura 11 se muestran los porcentajes de emisión para cada fuente de PM2.5, determinada mediante la aplicación del ACPA. En todas las estaciones el polvo del suelo fue el mayor emisor con 50 %, 73 % y 34 % para JA, CON y HYO, respectivamente. Para JA (33 %) y CON (23 %) la fundición fue el segundo en importancia de emisión, seguidos de fuentes indeterminadas. Para HYO la quema de biomasa fue la segunda en importancia de aporte (27 %) seguida de vehículos (26 %); el aporte de las fundiciones locales fue del 7 %.
Los dendrogramas obtenidos con los elementos químicos de las PM2.5 se muestran en la figura 12. Para el caso de JA se determinaron dos grupos: G1, con elementos trazadores de fundición (S, Pb, As, y Zn) y G2, con elementos trazadores de polvo del suelo (K, Al, Ca, Ti, Fe y Si). En CON, también se determinaron dos grupos: G1 con elementos trazadores de polvo del suelo (Br, K, Fe, Ca, Al y Si) y G2, con elementos trazadores de fundición (As, S, Zn y Ti). En HYO se determinaron cuatro grupos: G1 con un elemento trazador de vehículos (S), G2, con elementos trazadores de fundición (As y Pb), G3, con elementos trazadores de polvo del suelo (Al y Fe), y G4, con elementos trazadores de quema de biomasa (Br y K).
Estos resultados demuestran que a nivel de las PM10 y PM2.5 de los centros urbanos del VM se tienen aportes de fuentes antrópicas y naturales, siendo mayor el número de fuentes antrópicas, cuya contribución se puede controlar o disminuir a través de la implementación de planes y programas de gestión de calidad de aire.
CONCLUSIONES
Se determinó que las concentraciones PM10 y PM2.5 superan a las de la normatividad peruana vigente, lo que demuestra que las partículas atmosféricas representa un riesgo para los pobladores de los centros urbanos del VM. Siendo la ciudad de Huancayo la que registró mayor problemática por superar más veces los límites establecidos.
A través de la aplicación de FE se obtuvieron elementos trazadores de fuentes naturales (FE < 1): Si, Al, Ti y Fe y elementos de origen antrópico (FE > 10): As, Pb, Br, Zn y S, tanto en PM10 como en PM2.5.
En el caso de las fuentes de emisión de PM del VM en total se han registrado cinco fuentes: polvo del suelo (Al, Ca, Si, Fe, Ti, Mn y K), quema de biomasa (Cl, Br, K), vehículos (Cu, Zn, Cl, Cr), combustible-aceite (Ni) y fundición (Pb, Zn, As y Cu). Se logró determinar mayor número de fuentes de PM10 por el mayor número de datos registrados.
La principal fuente de emisión de PM10 y PM2.5 para todas las áreas urbanas evaluadas fue el polvo del suelo, seguida de las quemas de biomasa, lo que demuestra la influencia de áreas rurales circundantes y las quemas en zonas agrícolas y ganaderas. También las emisiones del parque automotor y la quema de combustible-aceite tienen gran importancia.
Otra fuente de importancia para las poblaciones del VM fue la fundición, a la que están asociados dos elementos tóxicos para la salud de la población como el Pb y el As.