1. Introducción
En los últimos años, se ha observado un incremento en el número de trabajos orientados a identificar y detectar sitios con problemas de contaminación por concentración de metales pesados mediante el uso de parámetros magnéticos, esto se debe mayormente a que las técnicas magnéticas suelen ser de bajo costo, rápidas y no generan residuos tóxicos (Beckwith et al., 1986; Chaparro et al., 2004; Aguilar et al., 2011; Aguilar et al., 2013; Bautista et al., 2014; Limo et al., 2018; Declerq et al., 2019).
Los estudios de magnetismo ambiental emplean suelo superficial, material partículado (polvos urbanos o sedimentos) o material botánico (hojas, líquenes, etc.) para determinar el valor de su susceptibilidad específica (Xlf), magnetización remanente isotermal de saturación (MRIS) (Evans y Heller, 2003) y concentración de metales pesado (Cr, Cu, Ni, Pb, V, Zn y otros), con el objetivo de establecer relaciones entre estas variables y su uso como un método de monitoreo ambiental proxy (Maher et al., 2008; Qiao et al., 2013; Szczepaniak et al., 2016; Rachwal et al., 2017; Limo et al., 2018; Declercq et al., 2019).
Las relaciones que se han encontrado entre los parámetros magnéticos y la concentración de metales pesados son de tipo lineal, con coeficientes de correlación (r) de Pearson entre: 0.2 a 0.90 (Zhang et al, 2012; Zhu et al., 2013; Dytlow et al., 2019; Yang et al., 2019; Wang et al., 2020). Las relaciones entre los parámetros magnéticos y la concentración de metales pesados son distintas en cada área de estudio, por lo que no es posible establecer grados de comparación en lugares diferentes (Kim et al., 2009; Morton et al., 2009; Wei y Yang, 2010; Qiao et al., 2013; Limo et al., 2018; Declercq et al., 2019).
Los modelos magnéticos que predicen la concentración de metales pesados pueden llegar a presentar variaciones menores al 30% (entre el valor medido y el valor estimado de la concentración), por lo que pueden ser métodos proxy confiables (Kim et al., 2009; Qian et al., 2011; Jordanova et al., 2013; Jones et al., 2015).
El material partículado urbano o polvo urbano contiene diferentes materiales (magnéticos y no magnéticos) que provienen de fuentes antrópicas (emisiones vehiculares, emisiones de fabricas, entre otros) y naturales (suelo, caída de ceniza volcánica, entre otros), el material magnético de origen antrópico incrementa la señal magnética de cualquier muestra ambiental debido a la adición de material ferrimagnético de baja coercitividad, por lo que es posible identificar los sitios con mayor o menor actividad antrópica a partir de la acumulación de material magnético o del incremento en la señal magnética (Thompson y Oldfield, 1986; Evans y Heller, 2003).
La hipótesis principal de los estudios de magnetismo ambiental en entornos urbanos propone que, los minerales ferrimagnéticos son portadores de iones de metales pesados, por lo que, si el polvo urbano tiene minerales ferrimagnéticos, es probable que también contenga metales pesados (Jordanova et al., 2003; Gautam et al., 2004; Jelenska et al., 2008; Kim et al., 2009; Qian et al., 2011; Rai, 2013).
Los modelos que estiman las concentraciones de elementos mediante los parámetros magnéticos suelen omitir su grado de confiabilidad, además de que, no existe un valor de umbral magnético que permita examinar rápidamente si un área contienen concentración alta de metales pesados. Por otra parte, la distribución espacial de los parámetros magnéticos en las ciudades permite identificar los sitios con mayor concentración de materiales magnéticos y en consecuencia de metales pesados (Davila et al., 2006; Maher et al., 2008; Halsall et al., 2008; Mitchell et al., 2009; Morton et al., 2009;Jordanova et al., 2012; Qiao et al., 2013; Cejudo et al., 2015).
Este estudio tuvo como objetivo identificar los parámetros magnéticos que manifiestan una relación con las concentraciones de metales pesados para polvos urbanos, con la finalidad de establecer modelos matemáticos predictivos, así como proponer valores de umbrales magnéticos que permitan identificar áreas con una concentración baja, media y alta de metales pesados mediante los parámetros magnéticos.
2. Área de estudio
Las muestras de polvo urbano fueron recolectadas en la Ciudad de México (CDMX), la cual tiene cerca de 9 millones de habitantes, y una superficie de 1,485 km2 (INEGI, 2020). Existen estudios previos sobre propiedades magnéticas, metales pesados y su distribución espacial que sugieren la existencia de relaciones estadísticamente significativas entre los parámetros magnéticos, metales pesado y su distribución geográfica (Morton et al., 2009; Rodríguez et al., 2011; Cejudo et al., 2015).
La CDMX en el sistema de coordenadas Universal Transversal de Mercator (UTM WGS84) Zona 14 N se encuentra entre 473,351 a 501,401 E y 2,128,050 a 2,164,290 N, en la región del Valle de México y rodeada por la cordillera de Guadalupe, a una altitud de 2,800 metros sobre el nivel del mar. En la parte SE y a 89.9 km de la CDMX se encuentra el volcán Popocatépetl, el cual se encuentra activo y emite ceniza volcánica que puede llegar a la ciudad. El clima es templado con inviernos secos (clasificación del clima Köppen-Geiger: Cwb), con temperatura media de 17.3 °C, precipitación acumulada mensual de 10 mm (Conagua 2017), la dirección del viento de N-S y NW-SE durante Abril y Mayo (SEDEMA, 2018), tiempo durante el cual se llevó a cabo la campaña de muestreo.
La muestra de polvo urbano (PU) se obtuvo a partir del material recolectado dentro de un área de un metro cuadrado sobre el pavimento, el cual fue levantado por medio de un cepillo y recogedor plástico, los fragmentos de fracción gruesa (> 4 mm) fueron retirados en el sitio. Posteriormente, la muestra se colocó dentro de una bolsa de polipropileno para su almacenamiento y traslado al laboratorio, las coordenadas del sitio fueron tomadas en el sistema UTM WGS1984. Se recolectaron 374 muestras de polvo urbano distribuidas sobre la superficie de la Ciudad de México, la distancia entre cada sitio fue de 500 m aproximadamente (Figura 1).
En el laboratorio, las muestras se tamizaron en malla del número N°35 (0.5 mm), el material obtenido se empleó para los análisis magnéticos y la determinación de la concentración de Al, Cr, Cu, Fe, Ni, Pb, V y Zn.
2.1. MEDICIONES MAGNÉTICAS
Las mediciones de susceptibilidad magnética y magnetización remanente se hicieron con material tamizado que se colocó dentro de una capsula cúbica de acrílico de 8 cm3.
2.2. SUSCEPTIBILIDAD MAGNÉTICA
La susceptibilidad magnética (K) fue medida a 0.46 kHz (Klf) y 4.6 kHz (Khf) en un medidor de susceptibilidad Bartington MS3 con sensor dual MS2B a temperatura ambiente. Se determinó el parámetro de susceptibilidad específica (Xlf=Klf/ρ, ρ es la densidad del material en kg m-3); el porcentaje de la susceptibilidad magnética dependiente de la frecuencia (Xlf% = [(Klf-Khf)/ Klf]*100), el cual permite identificar muestras con minerales ferrimagnéticos superparamagnéticos (SP) ultra finos (< 0.3 μm). La concentración de partículas SP se define como: Baja si, Xlf% < 2.0; mixta (hay mezcla de partículas SP con partículas más gruesas entre 0.2 a 110 µm), si Xlf% ∈ [2.0,10.0]; alta (>75% de partículas SP), si Xlf% ∈ [10.0,14.0] y para Xlf% > 14 se les considera valores extraños pocas veces reportados (Dearing, 1999).
2.3. MAGNETIZACIÓN REMANENTE ISOTERMAL
Se hizo la medición de la magnetización remanente adquirida a 1 T, la cual fue considerada como el valor de la magnetización remanente isotermal de saturación (MRIS). La aplicación del campo magnético se hizo con un magnetizador de pulso IM-10 ASC Scientific y la magnetización remanente se midio con un magnetómetro de giro JR6A (AGICO). Para determinar el valor de cociente S-₃₀₀, a la muestra se le aplicó un campo magnético en dirección opuesta de la MRIS de 300 mT. Posteriormente se midió su magnetización remanente (IRM-₃₀₀), el valor de S-₃₀₀ fue calculado como: S-₃₀₀=MRI-₃₀₀/MRIS. Si la muestra contienen una mayor concentración de minerales de baja coercitividad o ferrimagnéticos, entonces S-₃₀₀ exhibirá valores entre 0.7 a 1.0; si la muestra contiene una mayor concentración de minerales de alta coercitividad magnética o antiferromagnéticos, entonces el valor de S-₃₀₀ exhibirá valores entre 0.0 a 0.7 (Evans y Heller, 2003; Frank y Nowaczyk, 2008).
2.4. CURVAS TERMOMAGNÉTICAS (kT)
La identificación de los portadores magnéticos en el PU se determinaron por medio del cambio de la susceptibilidad magnética en función de la temperatura. Se utilizó un medidor de susceptibilidad magnética AGICO Kappabridge modelo KLY5. Las muestras fueron calentadas de 35 °C a 600 °C, con un ritmo de calentamiento de 20° min-1 y enfriadas con la misma razón de cambio. Las fases magnéticas se determinaron por el método diferencial (d2M/dT2) descrito por Tauxe et al., (2018).
2.5. COMBINACIÓN DE PARÁMETROS MAGNÉTICOS
La combinación de los parámetros magnéticos proporciona información adicional de las características de los minerales magnéticos presentes en la muestra de PU. La combinación de Xlf vs MRIS permite estimar la concentración del material ferrimagnético (magnetita) y el tamaño de grano magnético, donde, las muestras con mayor concentración de materiales magnéticos de baja coercitividad se encuentran en la parte superior derecha del gráfico Xlf vs MRIS, mientras que, las muestras con menor concentración de ferrimagnéticos tienden a ubicarse en la región inferior izquierda, el tamaño de la partícula magnética crece de derecha a izquierda (Thompson y Oldfield, 1986).
2.6. ANÁLISIS GEOQUÍMICO
La determinación de la concentración de Al, Cr, Cu, Fe, Ni, Pb V y Zn en las 374 muestras de PU, se realizó con un espectrómetro de fluorescencia de rayos X de energía dispersa (FRX-ED) portátil: Genius 9000 Skyray (Skyray, 2019), el cual tiene un detector de SDD y resolución mínima del detector < 139 eV, límite de detección 10 mg kg-1. La curva de calibración fue hecha con los estándar de serie IGL: GLS-1, IGLa-1, IGLbb-3, obteniendo un coeficiente de correlación (r) de 0.98 (Lozano y Bernal, 2005).
Las mediciones se hicieron en material tamizado con malla N° 230, en seco y colocado en una copa plástica. El espectrómetro se orientó perpendicularmente a la muestra a una distancia de 2 mm, cada medición se hizo por triplicado a temperatura ambiente y con una exposición de 30 segundos. El protocolo de medición se basa en la recomendación sugerida en el método 6200 propuesto por EPA (US-EPA, 2007) y el manual de operación Skyray (Skyray, 2019).
2.7. ANÁLISIS ESTADÍSTICO
La estadística descriptiva, el análisis exploratorio de datos (AED), los coeficientes de correlación y los modelos de estimación se realizaron con el software Statgraphics© V18.
Con la finalidad de identificar las concentraciones de los elementos (Al, Cr, Cu, Fe, Ni, Pb, V y Zn) que muestran relación con los parámetros magnéticos (Xlf, y MRIS), se llevó a cabo un análisis de componente principales (ACP). Debido a que el ACP puede verse afectado por los valores atípicos, se hizo un AED y se hicieron prueba de normalidad. El ACP fue hecho con el software R.
El polvo urbano es una mezcla de suelo y material antrópico, y actualmente, no hay valores de fondo que determinen su contaminación. En este trabajo, se exploró el uso de los valores de fondo que hay para suelo de CDMX propuestos por Morton et al., (2009b), y el valor de la mediana obtenido de 25 muestras de usos de suelo de baja actividad antrópica y la metodología propuesta por Rueda et al., (2011). Los valores de fondo son propios para cada área de estudio, debido a la diversidad geológica y ambiental que existe. Por lo tanto, resulta difícil de aplicar o establecer un criterio general para decretarlos.
El factor de concentración (CF, siglas en inglés) se determinó por medio de los valores de fondo obtenidos de los usos de suelo con baja actividad antrópica y se calculó con (Ecuación 1):
Donde Cm es la concentración medida del m elemento de la muestra y Cm fondo es el valor de fondo para el elemento m.
El índice de carga contaminante (PLI , siglas en inglés) se determinó mediante (Ecuación 2):
Donde, CF es el factor de contaminación, n representa el número de metales analizados. Las categorías de PLI fueron: PLI < 1, entonces la concentración de metales pesados se encuentra en estado natural; PLI ∈ [1,2], entonces la concentración de metales es moderada y si PLI > 2, entonces la concentración de metales pesados es significativa (Tomlinson et al., 1980; Angulo, 1996).
2.8. ANÁLISIS ESPACIAL
El análisis estructuran fue hecho con el software GS+, el cual permitió obtener el semivariograma teórico para poder aplicar un modelo de interpolación espacial para los parámetros magnético, concentración de elementos y PLI de tipo Kriging ordinario (KO) (Webster y Oliver, 2007; Robertson, 2015).
La representación geográfica se hizo con el software Arcgis© V.9.
3. Resultados y discusión
3.1. RESULTADOS MAGNÉTICOS
La estadística descriptiva de los parámetros mag néticos se presenta en la Tabla 1.
Xlf | Xdf % | MRIS | S-300 | |
Máximo | 6.52 | 7.09 | 81.43 | 1.0 |
Mínimo | 1.51 | 0.39 | 16.40 | 0.7 |
Promedio | 3.31 | 2.31 | 44.28 | 0.9 |
Mediana | 3.18 | 2.14 | 43.08 | 0.9 |
Desviación | 0.95 | 0.90 | 13.52 | 0.1 |
Coeficiente | 0.29 | 0.39 | 0.31 | 0.1 |
Sesgo | 5.52 | 11.77 | 4.51 | -10.5 |
Curtosis | 1.16 | 15.57 | 0.70 | 3.5 |
Fondo | 1.50* | |||
Limite | 4.00* | 56.00** |
Al | Cr | Cu | Fe (mg kg-1) |
Ni | Pb | V | Zn | |
Máximo | 27336 | 117 | 405 | 39217 | 66 | 366 | 52 | 537 |
Mínimo | 3325 | 17 | 15 | 653 | 21 | 27 | 11 | 102 |
Promedio | 10577.4 | 47.4 | 94.7 | 4020.9 | 37.3 | 115.2 | 26.5 | 255.0 |
Mediana | 9194.3 | 43.7 | 81.8 | 3256.4 | 36.3 | 102.5 | 26.2 | 238.6 |
Desviación | 5522.0 | 16.9 | 57.2 | 3197.3 | 9.5 | 63.7 | 6.3 | 95.4 |
Coeficiente | 0.5 | 0.4 | 0.6 | 0.8 | 0.3 | 0.6 | 0.2 | 0.4 |
Sesgo | 8.3 | 9.9 | 13.2 | 34.9 | 2.7 | 10.2 | 4.1 | 5.2 |
Curtosis | 2.0 | 8.4 | 17.5 | 156.4 | -1.7 | 7.6 | 3.2 | -0.2 |
Fondo | 105 | 32 | 56 | 19 | 87 | 76 | ||
Límite | 280i | 3100ii | 1600i | 400i | 78ii | 23000ii |
i) Valor recomendado por el gobierno de México (SEMARNAT, 2007). ii) Valor recomendado por la Agencia de Protección Ambiental de Estados Unidos (US-EPA Region 9, 2012). Valores de fondo: reportados por Morton et al., (2009). Xlf: (µm3 kg-1); Xdf% en %; MRIS: (mAm2 kg-1) y S-300 adimensional.
Los valores obtenidos de Xlf y MRIS de las 374 muestras del PU de CDMX fueron: la Xlf entre 1.51 a 6.52 µm3 kg-1, con valor de mediana de 3.18 µm3 kg-1; la MRIS con valores entre 16.40 a 81.43 mAm2 kg-1, con un valor de mediana de 43.08 mAm2 kg-1 que indicó una variación en el contenidode material magnético en el PU de CDMX (Tabla 1).
Los resultados de Xlf y MRIS mostraron un comportamiento lineal, en donde el aumento de los valores de Xlf y MRIS se debe al incremento del material ferrimagnético o magnetita en el PU. Se observó una dispersión de los valores de los parámetros magnéticos con respecto a la tendencia lineal debido a la presencia de partículas ferrimagnéticas de diversos tamaños (1 µm a 256 mm) (Thompson y Oldfield 1986; Evans y Heller 2003, Figura 2).
El Xlf% mostró valores entre 0.39 a 7.09 % y valor de mediana de 2.14 % , indicando concentración baja y media de partículas SP con mezclas de partículas de fracciones más gruesas (Dearing, 1999). El cociente S-300 mostró valores entre 0.7 y 1.0, indicando que el contenido de material ferrimagnético fue predominante sobre el material antiferromagnético en el PU. Los datos formaron dos conglomerados que sugiere la presencia de sitios con material ferrimagnético antrópico mezclado con materiales antiferromagnéticos. (Evans y Heller, 2003)(Figura 2 y Tabla 1).
La comparación de los valores de Xlf y MRIS de este estudio con los valores umbrales propuestos por Morton et al., (2009) para Xlf (4.0 µm3 kg-1) y por Cejudo et al., (2015) para MRIS (56 mAm2 kg-1) indican que hay sitios con alta concentración de elementos en CDMX (Figura 2).
3.2. CURVAS DE SUSCEPTIBILIDAD MAGNÉTICA
El análisis de las curvas kT del PU mostró que el principal portador magnético es magnetita (material ferrimagnético), con una temperatura de Curie cercana a 580 °C (Tauxe et al., 2018). Algunas curvas presentaron varias fases magnéticas durante el calentamiento (120 y 350 °C), que indican la presencia de minerales magnéticos, tales como: goethita o greigita, y cuyo origen puede ser de fuentes naturales (Evans y Heller, 2003; Tauxe et al., 2018) (Figura 3). En varias muestras, se observó un incremento de susceptibilidad en forma de cresta alrededor de 580 °C, comportamiento asociado a la neoformación de magnetita, cuyo origen puede ser de emisiones vehiculares (Figura 3).
3.3. RESULTADOS GEOQUÍMICOS
La estadística descriptiva de las concentraciones de Al, Cr, Cu, Fe, Ni, Pb, V y Zn se muestran en la Tabla 1.
La concentración de los elementos en las 374 muestras de PU fueron diferentes. La concentración de cada elemento fue: Al entre 3,325 a 27,336 mg kg-1; Cr entre 17 a 117 mg kg-1; Cu entre 15 a 405 mg kg-1; Fe entre 653 a 39,217 mg kg-1; Ni entre 21 a 66 mg kg-1; Pb entre 27 a 366 mg kg-1; V entre 11 a 52 mg kg-1 y Zn entre 102 a 537 mg kg-1.
Las concentraciones de los elementos se encontraron por debajo del valor límite permisible recomendado por las entidades gubernamentales y ambientales en CDMX (SEMARNAT 2007; USA-EPA 2012). Sin embargo, algunos sitios mostraron una concentración mayor al valor de fondo (VF) propuesta para suelo, como en: Cr (VF: 105 mg kg-1), Cu (VF: 32 mg kg-1), Ni (VF: 56 mg kg-1), Pb (VR: 19 mg kg-1), V (VF: 87 mg kg-1) y Zn (VF: 76 mg kg-1) (Morton et al., 2009; Rodríguez et al., 2011; Ihl et al., 2015). Actualmente, no existen valores de fondo y de referencia para el PU, por lo que se usa el valor de fondo del suelo para su evaluación, aún cuando el PU tiene una matriz diferente (Morton et al., 2009; Rodríguez et al., 2011; Ihl et al., 2015; Cejudo et al., 2015).
La comparación de valores de referencias de los elementos geoquímicos con otros lugares alrededor del mundo puede provocar errores de interpretación, de tal forma que solo son útiles como un marco global y no para indicar si hay contaminación en una área específica, (Tabla 2).
Ciudad | País | Cr | Cu (mg kg-1) |
Ni | Pb | Zn | Referencia |
Holanda | 100 | 36 | 35 | 85 | 140 | Brus et al., 2009 | |
Málaga | España | 132 | 65 | 58 | 69 | 132 | Castillo et al., 2002 |
Granada | España | 66 | 26 | 20 | 36 | 76 | Díez et al., 2009 |
Almería | España | 25 | 394 | Gil et al., 2002 | |||
Tirana | Albania | 114 | 36 | 42 | 86 | 151 | Gjoka et al., 2010 |
China | 59 | 32 | 28 | 38 | 118 | Wei y Yang, 2010 | |
Bejing | China | 30 | 19 | 27 | 25 | 58 | Chen et al., 2005 |
Brasil | 40 | 35 | 13 | 17 | 60 | Fadigas et al., 2006 | |
México | 105 | 32 | 56 | 19 | 76 | Morton et al., 2009 | |
México | 39 | 64 | 34 | 100 | 25 | Este estudio | |
Suelos del mundo | Promedio | 7-221 | 6-80 | 4-55 | 10-84 | 17-125 | McBride, 1994 |
En este estudio, se usaron como valores de fondo las concentraciones de suelo reportados por Morton et al., (2009) para CDMX y los valores de concentración de PU de las zonas de más baja actividad antrópica (Tabla 2).
3.4. ANÁLISIS DE COMPONENTES PRINCIPALES
Los resultados del ACP se presentan en la Tabla 3, los datos transformados por logaritmo natural presentaron una distribución normal. El ACP fue usado para expresar las concentraciones de elementos y parámetros magnéticos en un conjunto de componentes lineales (Figura 4). El ACP mostró que, la componente 1 y 2 representan el 60% de la varianza total, por lo que la información de las concentraciones y parámetros magnéticos queda representada con dos componentes. En el plano principal y el círculo de correlación (Figura 4) se observaron cuatro agrupaciones de datos del PU de CDMX. El Grupo I mostró una relación fuerte entre Al, Fe y V, y una relación inversamente proporcional con los parámetros magnéticos y Ni; el Grupo II presentó una correlación inversamente proporcional con el contenido de Cr, Cu, Pb y Zn; el Grupo III exhibió una relación entre Xlf, MRIS y Ni, pero inversamente proporcional con Al, Fe y V; el grupo IV presentó una fuerte relación entre Cr, Cu, Pb y Zn, elementos que se asocian con emisiones de tipo antrópico. Dado que, las concentraciones de Al, Fe y V mostraron una relación inversamente proporcional con los parámetros magnéticos, estos no se consideraron para la evaluación de PLI, debido a que se buscó estimar el índice de carga contaminante con el crecimiento de los valores de los parámetros magnéticos.
Componente | Varianza | % Var | % Var Acumulada |
Dim.1 | 3.26 | 32.61 | 33 |
Dim.2 | 2.70 | 27.01 | 60 |
Dim.3 | 1.19 | 11.92 | 72 |
Dim.4 | 0.77 | 7.74 | 79 |
Dim.5 | 0.64 | 6.35 | 86 |
Dim.6 | 0.46 | 4.56 | 90 |
Dim.7 | 0.35 | 3.45 | 94 |
Dim.8 | 0.30 | 2.97 | 97 |
Dim.9 | 0.25 | 2.52 | 99 |
Dim.10 | 0.09 | 0.87 | 100 |
Componente | Cos2 | ||||
Parámetro | Dim.1 | Dim.2 | Dim.1 | Dim.2 | Suma |
ln MRIS | 0.745 | -0.157 | 0.55 | 0.02 | 0.76 |
ln Cr | 0.583 | 0.564 | 0.34 | 0.32 | 0.81 |
ln Cu | 0.638 | 0.420 | 0.41 | 0.18 | 0.76 |
ln Ni | 0.549 | 0.081 | 0.30 | 0.01 | 0.55 |
ln Pb | 0.403 | 0.566 | 0.16 | 0.32 | 0.69 |
ln V | -0.426 | 0.690 | 0.18 | 0.48 | 0.81 |
ln Zn | 0.467 | 0.562 | 0.22 | 0.32 | 0.73 |
ln Al | -0.615 | 0.679 | 0.38 | 0.46 | 0.92 |
ln Fe | -0.449 | 0.773 | 0.20 | 0.60 | 0.89 |
3.5. CORRELACIÓN ENTRE PARÁMETROS MAGNÉTICOS Y METALES PESADOS
Los trabajos previos de magnetismo ambiental en suelo urbano de la Ciudad de México, mostraron una relación entre parámetros magnéticos y metales pesados, los cuales exhibieron coeficiente de correlación lineal de Pearson (r) entre 0.3 a 0.5 (Morton et al., 2009; Cejudo et al., 2015). En este trabajo, las correlaciones entre los parámetros magnéticos y elementos en PU se muestra en la Tabla 4, donde, se observó r positivos (entre 0.25 a 52) para Xlf y MRIS y las concentraciones de Cr, Cu, Ni, Zn y PLI, resultados similares a los reportados en los trabajos de Morton et al. (2009) y Cejudo et al., (2015).
ln Xlf | ln MRIS | ln Cr | ln Cu | ln Ni | ln Pb | ln V | ln Zn | ln Al | ln Fe | PLI | |
ln Xlf | 1.00 | ||||||||||
ln MRIS | 0.69 | 1.00 | |||||||||
ln Cr | 0.34 | 0.26 | 1.00 | ||||||||
0.00 | 0.00 | ||||||||||
ln Cu | 0.43 | 0.27 | 0.49 | 1.00 | |||||||
0.00 | 0.00 | 0.00 | |||||||||
ln Ni | 0.36 | 0.52 | 0.32 | 0.29 | 1.00 | ||||||
0.00 | 0.00 | 0.00 | 0.00 | ||||||||
ln Pb | 0.16 | 0.03 | 0.69 | 0.42 | 0.15 | 1.00 | |||||
0.00 | 0.55 | 0.00 | 0.00 | 0.01 | |||||||
ln V | -0.18 | -0.26 | 0.06 | -0.04 | 0.00 | 0.06 | 1.00 | ||||
0.00 | 0.00 | 0.28 | 0.50 | 0.97 | 0.28 | ||||||
ln Zn | 0.26 | 0.19 | 0.53 | 0.65 | 0.14 | 0.49 | 0.14 | 1.00 | |||
0.00 | 0.00 | 0.00 | 0.00 | 0.01 | 0.00 | 0.01 | |||||
ln Al | -0.31 | -0.47 | -0.04 | -0.11 | -0.27 | 0.09 | 0.68 | 0.07 | 1.00 | ||
0.00 | 0.00 | 0.42 | 0.04 | 0.00 | 0.09 | 0.00 | 0.23 | ||||
ln Fe | -0.09 | -0.30 | 0.14 | 0.00 | -0.09 | 0.10 | 0.75 | 0.10 | 0.92 | 1.00 | |
0.07 | 0.00 | 0.01 | 0.94 | 0.09 | 0.05 | 0.00 | 0.07 | 0.00 | |||
PLI | 0.39 | 0.25 | 0.81 | 0.81 | 0.39 | 0.78 | 0.05 | 0.78 | -0.05 | 0.10 | 1.00 |
0.00* | 0.00 | 0.00 | 0.00 | 0.00 | 0.00 | 0.30 | 0.00 | 0.31 | 0.06 |
* valor de P.
La concentración de Al, Fe y V mostraron una relación inversamente proporcional con Xlf y MRIS (Tabla 4), este resultado concuerda con los datos obtenidos del ACP (Figura 4). La correlación inversa del Al y V se debe probablemente a las características paramagnéticas de estos elementos. Por otra parte, los minerales antiferromagnéticos con contenido Fe pueden provocar una baja señal magnética (Evans y Heller, 2003; Rengifo et al., 2013, Tauxe et al., 2018).
3.6. MODELO MATEMÁTICOS Y UMBRALES MAGNÉTICOS
Los modelos matemáticos para estimar la concentración de elementos (Cr, Cu, Ni, Pb y Zn) y el valor del PLI en el PU, por medio de los parámetros magnéticos, se muestran en la Tabla 5.
Parámetros | R | R2 | Modelo matemático | n obs |
n pred |
Precisión |
ln Xlf vs ln MRIS | 0.7 | 47.1 | ln MRIS = 2.89426 + 0.737451*ln Xlf | 370 | 289 | 78% |
ln Xlf vs ln Ni | 0.4 | 12.7 | ln Ni = 3.1819 + 0.344504*ln Xlf | 349 | 259 | 74% |
ln Xlf vs ln Zn | 0.3 | 6.6 | ln Zn = 5.07524 + 0.339069*ln Xlf | 343 | 211 | 62% |
ln Xlf vs ln Cr | 0.3 | 11.2 | ln Cr = 3.34864 + 0.391169*ln Xlf | 372 | 228 | 61% |
ln Xlf vs PLI | 0.4 | 14.9 | PLI = 0.802279 + 0.747927*ln Xlf | 374 | 225 | 60% |
ln Xlf vs ln Cu | 0.4 | 18.4 | ln Cu = 3.35683 + 0.886519*ln Xlf | 372 | 151 | 41% |
ln Xlf vs ln Al | -0.3 | 9.6 | ln Al = 9.77 - 0.542919*ln Xlf | 364 | 147 | 40% |
ln Xlf vs ln Pb | 0.2 | 2.5 | ln Pb = 4.256 + 0.298531*ln Xlf | - | - | - |
ln Xlf vs ln V | -0.2 | 3.3 | ln V = 3.42789 - 0.15255*ln Xlf | - | - | - |
ln MRIS vs ln Ni | 0.5 | 27.2 | ln Ni = 1.79952 + 0.474246*ln MRIS | 347 | 265 | 76% |
ln MRIS vs ln V | -0.3 | 6.7 | ln V = 4.02728 - 0.206784*ln MRIS | 370 | 281 | 76% |
ln MRIS vs PLI | 0.3 | 7.4 | PLI = -0.337863 + 0.536596*ln MRIS | 370 | 228 | 62% |
ln MRIS vs ln Cr | 0.3 | 6.8 | ln Cr = 2.7151 + 0.290241*ln MRIS | 368 | 209 | 57% |
ln MRIS vs ln Zn | 0.2 | -3.5 | ln Zn = 4.58362 + 0.236994*ln MRIS | 339 | 170 | 50% |
ln MRIS vs ln Cu | 0.3 | 7.5 | ln Cu = 2.36553 + 0.538765*ln MRIS | 368 | 126 | 34% |
ln MRIS vs ln Al | -0.5 | 22.4 | ln Al = 12.1375 - 0.798593*ln MRIS | 361 | 170 | 47% |
ln MRIS vs ln Fe | -0.3 | 9.2 | ln Fe = 10.5159 - 0.649695*ln MRIS | 370 | 124 | 34% |
Xlf: (µm3 kg-1); MRIS: (mAm2 kg-1).
Los modelos propuestos para estimar la concentración de Cr, Ni, Zn y PLI, por medio de Xlf y MRIS, mostraron una precisión de 57 a 74%. El PLI en PU puede ser estimado con Xlf o MRIS con una precisión del 60% (Figura 5).
En 2009, Morton et al., (2009) propusieron un valor umbral magnético de Xlf: 4.0 µm3 kg-1 y Cejudo et al., (2015) un valor de MRIS de 46 mAm2 kg-1 para identificar sitios con altas concentraciones de metales pesados. Para proponer, los nuevos valores umbrales, se definió los límites con el valor de PLI; si PLI ≤1, la concentración de elementos es baja; si PLI muestra valores entre 1 y 2, entonces la concentración de elementos es moderada y si PLI >2, entonces la concentración de elementos es alta. En cada modelo se determinó el valor medio del intervalo de PLI para determinar un valor umbral en termino del parámetro magnético. Los resultados se muestran en la Tabla 6.
Parámetro | PLI <1 |
PLI [1-2] |
PLI >3 |
Modelo | Parámetro magnético | ||
ln Cr | 3.44 | 3.69 | 4.15 | Xlf = exp (0.0629155 + 0.287546*ln Cr) | 2.86 | 3.08 | 3.51 |
ln Cu | 3.44 | 4.28 | 4.93 | Xlf = exp (0.247879 + 0.20769*ln Cu) | 2.62 | 3.12 | 3.56 |
ln Ni | 3.40 | 3.56 | 3.69 | Xlf = exp ( -0.143548 + 0.36757*ln Ni) | 3.02 | 3.21 | 3.36 |
ln Pb | 3.82 | 4.46 | 5.10 | Xlf = exp ( 0.789394 + 0.0820755*ln Pb) | 3.01 | 3.18 | 3.35 |
ln Zn | 4.79 | 5.36 | 5.85 | Xlf = exp ( 0.110147 + 0.193529*ln Zn) | 2.82 | 3.15 | 3.46 |
Mediana | 2.86 | 3.15 | 3.46 | ||||
ln Cr | 3.44 | 3.69 | 4.15 | MRIS = exp( 2.85859 + 0.234366*ln Cr) | 39.1 | 41.40 | 46.1 |
ln Cu | 3.44 | 4.28 | 4.93 | MRIS = exp(3.13983+0.139792*ln Cu) | 37.4 | 101.6 | 276.1 |
ln Ni | 3.40 | 3.56 | 3.69 | MRIS = exp( 1.71692 + 0.572697*ln Ni) | 39.0 | 42.8 | 46.1 |
ln Zn | 4.79 | 5.36 | 5.85 | MRIS = exp (2.84432 + 0.162912*ln Zn) | 37.5 | 41.2 | 44.6 |
Mediana | 38 | 42 | 46 |
Xlf: (µm3 kg-1); MRIS: (mAm2 kg-1).
Los valores umbrales que se proponen para estiman la concentración de Cr, Cu, Pb y Zn son: para Xlf: Baja <2.86 µm3 kg-1; Moderada: 2.86 a 3.46 µm3 kg-1 y Alta >3.46 µm3 kg-1; para MRIS: Baja <38 mAm2 kg-1; Moderada: entre 38 a 46 mAm2 kg-1 y Alta >46 mAm2 kg-1. Se encontró que entre el 36 a 40% de las muestras de PU de CDMX presentaron valores mayores al valor umbral de Xlf y MRIS. Por lo tanto, hay un número considerables de sitios que tienen una concentración alta de elementos (Tabla 6). Por otro lado, los cálculos con los valores umbrales magnéticos de Morton et al., (2009) y Cejudo et al., (2015), mostraron que entre el 20% y 17% de los sitios tienen probablemente un problema de concentración metales pesados.
3.7. DISTRIBUCIÓN ESPACIAL
Con la finalidad de establecer la relación espacial que hay entre la acumulación de material magnético y la concentración de elementos, se hizo un análisis geoestadístico para cada parámetro magnético, concentración de elemento y valor de PLI. Los resultados del análisis estructural indicaron : semivariogramas teóricos de tipo esférico y exponencial, con valores de nugget entre 0.003 y 0.013 y sill entre 0.8 a 0.14, con una varianza estructural > 91% y un rango entre 1.5 a 1.6 km, que indica una correlación espacial alta (Tabla 7, Figura 6). La representación espacial de los parámetros magnéticos muestra que, el mineral magnético se distribuye preferentemente en la dirección NO-SE, con una mayor acumulación en la parte SE de la ciudad (180 km2 aprox.), la cual puede estar influenciada por la caída de ceniza volcánica, ya que en esa dirección se encuentra un volcán activo (Popocatépetl). La parte central de la ciudad mostró pequeñas áreas (12 km2), en donde se acumula material magnético (Figura 6). La acumulación del material magnético de este estudio fue similar a la distribución espacia reportado por Cejudo et al., (2015).
Parámetro | Modelo | Nugget | Sill | C | varianza estructural % |
R2 | Rango (km) |
Xlf | Exponencial | 0.004 | 0.08 | 0.08 | 0.95 | 0.86 | 1.49 |
MRIS | Esférico | 0.003 | 0.07 | 0.07 | 0.96 | 0.96 | 1.56 |
Cr | Esférico | 0.003 | 0.10 | 0.10 | 0.97 | 0.91 | 1.42 |
Zn | Esférico | 0.003 | 0.14 | 0.13 | 0.98 | 0.97 | 0.99 |
Ni | Esférico | 0.000 | 0.06 | 0.06 | 1.00 | 0.97 | 1.40 |
PLI | Exponencial | 0.013 | 0.14 | 0.13 | 0.91 | 0.86 | 1.68 |
Los modelos de distribución espacial del Cr y Zn indicaron distribuciones similares en la ciudad, y algunas áreas de mayor acumulación en la parte norte de la ciudad (Figura 7). Mientras que, el Ni se concentró en la región centro y sur de la ciudad.
La distribución espacial del PLI mostró que, la mayor parte de la ciudad presenta valores de PLI entre 1 y 2, y una zona pequeña de 5 km2 al SE registro valores de PLI < 1. Las zonas con valores de PLI > 2, se ubican a lo largo de la ciudad, en dirección NO-SE y algunas zonas en dirección SO-NE (Figura 7).
4. Conclusiones
El análisis de propiedades magnéticas hecho en polvo urbano de la Ciudad de México mostró la presencia de material ferrimagnético mezclado con material antiferromagnético, siendo el principal portador magnético la magnetita, se detectó la presencia de partículas superparamagnéticas mezcladas con partículas de tamaño más grueso. El estudio magnético del polvo urbano de CDMX contiene material magnético de fuentes antrópicas y de fuentes naturales, la presencia de partículas ferrimagnéticas metaestables indica la adición de material magnético de fuentes vehiculares. Por otra parte, la detección de material antiferromagnético muestra la adición probable de polvo proveniente del suelo de la Ciudad de México.
El polvo urbano de CDMX tiene concentraciones variables de elementos mayores y metales pesados. El contenido de los elementos del polvo urbano de mayor a menor, fue: Al > Fe > Zn > Pb > Cu > Cr > Ni > V.
El análisis multivariante mostró que, hay una relación entre parámetros magnéticos y la concentración de elementos. Sin embargo, no todas las relaciones fueron del tipo directamente proporcional y positiva: La concentración de Al y V mostraron un comportamiento inverso con los parámetros magnéticos debido a su comportamiento paramagnéticos, mientras que, la concentración de Cr, Cu, Ni y Zn se logró estimar con los parámetros magnéticos con una precisión entre 60 a 70%. Se puede estimar el índice de carga contaminante con parámetros magnéticos. El modelo de regresión lineal es capaz de estimar el valor del índice de carga contaminante con una probabilidad del 60%.
Se proponen nuevos valores umbrales magnéticos para identificar sitios con problemas de acumulación de elementos como Cr, Cu, Ni, Pb y Zn en polvo urbano para la Ciudad de México, los cuales proporcionan un grado de eficacia del 60%, como: Xlf: Baja < 2.86 µm3 kg-1; Moderada: 2.86 a 3.46 µm3 kg-1 y Alta > 3.46 µm3 kg-1; para MRIS: Baja < 38 mAm2 kg-1; Moderada: entre 38 a 46 mAm2 kg-1 y Alta > 46 mAm2 kg-1.
El modelo de distribución espacial de susceptibilidad magnética, magnetización remanente isotermal y el índice de carga contaminante mostraron que la dirección noreste-sureste de la Ciudad de México es más propensa a exhibir concentraciones altas de metales pesados. La influencia de la ceniza del volcán Popocatépetl se ve reflejada en la acumulación de material magnético en la zona sureste y fue detectada con los parámetros magnéticos.
Este estudio, proponen nuevos modelos matemáticos que estiman la concentración de elemento e índice de carga contaminante a través de los parámetros magnéticos con un grado de precisión del 60%. Los umbrales magnéticos que se proponen para polvo urbano podrían usarse como proxy ambientales.
Contribuciones de los autores
Conceptualización: R. Cejudo (RC), F. Bautista (FB) y A. Goguitchaichvili (AG); Análisis o adquisición de datos: RC, FB, AG y M.Cervantes (MC); Desarrollo metodológico/técnico: RC, FB, AG y MC; Redacción del manuscrito original: RC, AG, FB; Redacción del manuscrito corregido y editado: RC, FB y AG Diseño gráfico: RC, FB y AG; Trabajo de campo: RC; Interpretación: RC, FB, AG y MC.
Conflictos de interés
No existen conflictos de interés con otros autores, instituciones u otros terceros sobre el contenido del presente artículo.