nueva página del texto (beta)
Español (pdf)
Artículo en XML
Referencias del artículo
Traducción automática
Enviar artículo por email
Citado por SciELO 
Similares en
SciELO 
Permalink
Introducción
La contaminación por metales pesados es un problema que se presenta en prácticamente todas las ciudades del mundo, debido al desarrollo económico e industrial. En los últimos años, se ha observado un aumento en las emisiones de partículas suspendidas que contienen metales pesados (MP) por parte de fuentes antrópicas (emisiones vehiculares e industriales) y, en menor medida, por fuentes naturales (actividad volcánica, tormentas de arena, incendios forestales, entre otros). Las partículas en suspensión transportan MP y otras sustancias tóxicas que, al precipitarse sobre la superficie, generan un incremento en las concentraciones naturales del suelo, sedimentos y cuerpos de agua (superficiales y subterráneos), que pueden alcanzar niveles de toxicidad perjudiciales para el ecosistema y, en el caso de suelos urbanos, pueden representar una amenaza para la población (Yap et al. 2019, Wang et al. 2021, Salazar-Rojas et al. 2023).
Los elementos como arsénico (As), cadmio (Cd), cobre (Cu), cromo (Cr), hierro (Fe), manganeso (Mn), níquel (Ni), plomo (Pb), vanadio (V) y zinc (Zn), conocidos como metales pesados y metaloides o elementos potencialmente tóxicos, se encuentran de forma natural en la corteza terrestre (suelo y rocas) en concentraciones que varían de un lugar a otro (Mirsal 2004, Kabata-Pendias y Mukherjee 2007, US-EPA 2024). En suelos urbanos, las concentraciones naturales de elementos se alteran (aumentan) a causa de la precipitación de MP derivada de actividades antrópicas (Vareda et al. 2019, Ahmad et al. 2021, Aguilera et al. 2022).
Recientemente, se ha empleado el término “elementos potencialmente tóxicos” (EPT) en la literatura para referirse a un grupo de metales pesados, considerando los elementos esenciales y no esenciales que requieren los organismos para el desarrollo de sus funciones vitales, tales como Cr, Cu, Fe, Mn y otros. Sin embargo, concentraciones elevadas de estos elementos en el organismo resultan dañinas y el Pb, un metal pesado que no se incluye en estos EPT, es altamente tóxico a cualquier concentración (WHO 2023). Para la cuantificación de concentraciones de elementos como Cr, Cu, Ni, Pb, V y Zn en el suelo es preferible el empleo del término metales pesados (Vareda et al. 2019, Zhang et al. 2022, WHO 2023).
El crecimiento urbano acelerado característico de muchas ciudades del mundo limita la capacidad de identificar con precisión la concentración de fondo total (CFT) para evaluar el enriquecimiento de MP en suelos. Varios estudios evalúan el enriquecimiento de elementos en suelo urbano mediante valores de referencia globales (VRG), valores de referencia de unidades geológicas o edafológicas similares o con valores de concentración especifica total (CET), lo cual puede causar ambigüedad en la evaluación del enriquecimiento de elementos del suelo (Morton-Bermea et al. 2009, Aguilera et al. 2022, 2024).
La información disponible sobre el valor de concentración de fondo total (CFT) de MP en suelos urbanos es muy escasa debido a la diversidad edafológica y geológica o la ausencia de ésta, dificultando la evaluación precisa de la amenaza o el riesgo a la salud y al ambiente. La propuesta para cuantificar el valor de CFT de elementos en suelo consiste en obtener muestras de suelo dispersas por la ciudad que provengan de usos de suelo con baja actividad antrópica (áreas verdes o residenciales), con el fin de medir la concentración de elementos mediante un modelo estadístico. Se toma como CFT la media de la población cuando los valores obtenidos presenten una distribución normal; en caso de presentarse una distribución diferente se considera el valor de la mediana (SEMARNAT 2007).
En México hay muy pocos datos sobre la concentración de elementos en suelos, lo que dificulta evaluar la contaminación en diferentes partes del país. La presente investigación se orientó a realizar un análisis de la concentración de elementos y de magnetismo ambiental en suelos superficiales de la Zona Metropolitana de Guadalajara, con el fin de identificar y cuantificar el estado de enriquecimiento de elementos potencialmente tóxicos en el suelo por medio de la concentración de fondo local. Se busca también evaluar la relación entre parámetros magnéticos y concentración de elementos en el suelo, para establecer un método de monitoreo mediante proxys que permita identificar suelos contaminados mediante el uso de mediciones magnéticas.
MATERIALES Y MÉTODOS
Área de estudio
El área de estudio comprende la superficie urbana de la Zona Metropolitana de Guadalajara (ZMG), la cual ocupa seis municipios: San Pedro Tlaquepaque, Tonalá, Zapopan, Tlajomulco de Zúñiga, El Salto, Acatlán de Juárez y Zapotlanejo, ubicados en la región central del estado de Jalisco dentro del Valle de Atemajac, entre las coordenadas 20º45’-20º28’ N y 103º31’-103º13’ W, con una altitud promedio de 1566 msnm, un área de 1564 km2 y población de 5 266 642 habitantes (INEGI 2022; Fig. 1). La superficie urbana está edificada sobre rocas (brecha y basalto) de los periodos Neógeno y Cuaternario, y sepulta suelos del grupo Regosol, Feozem, Fluvisol y Litosol. En la ZMG se presenta un clima templado del grupo semicálido subhúmedo (ACw en la clasificación de Koppen), una precipitación < 35 mm y temperatura promedio de 22 ºC en el mes de marzo (INEGI 2022a, b, GDL 2024, SGM 2024a, b).
Fig. 1 Sitios de muestreo de suelo urbano en la Zona Metropolitana de Guadalajara y población en las áreas geográficas básica (AGEB) en 2015.
El levantamiento de suelo urbano (SU) de la ZMG se realizó el 23 de marzo de 2015, empleando un muestreo sistemático por rejilla rectangular, con separación de 600 m entre sitios (Fig. 1). Se recolectaron 140 muestras de SU superficial, que fueron obtenidas con un nucleador cilíndrico de plástico (10 cm de largo y 5 cm de diámetro), el cual se introdujo a 5 cm del suelo mediante percusión, obteniéndose dos muestras por sitio (~ 200 cm3). In situ, se retiró la materia vegetal y elementos no pertenecientes al suelo; cada muestra se colocó en una bolsa de polipropileno con cierre para su traslado. Adicionalmente, se recolectaron datos de uso de suelo y ubicación geográfica mediante un instrumento geoposicionador (GPS) en el sistema de coordenadas UTM WGS 84.
La clasificación del uso de suelo se elaboró de acuerdo con la norma urbanística del plan de ordenación municipal (IMEPLAN 2016) en cinco categorías: áreas verdes (áreas con extensión > 200 m2 con cubierta vegetal nativa), comercio y servicios (áreas que tienen mercados, plazas comerciales o venta de consumo), equipamiento (escuelas, centros de recreación, áreas deportivas, iglesias y hospitales), industrial (áreas con fábricas) y habitacional (viviendas).
En el laboratorio, las muestras fueron secadas a temperatura ambiente (20-22 ºC), humedad relativa de (30 a 50 %) y en ausencia de luz por dos semanas en charolas plásticas. Una vez secada la muestra, se destruyen los agregados con mazo de madera y se retiran los fragmentos de piedra y plantas. Seguidamente, la muestra fue tamizada con redes plásticas: ASTM núm. 10 (ISO: 2 mm), ASTM núm 40 (ISO: 0.425 mm) y ASTM núm 200 (0.075 mm), con el fin de homogenizar el tamaño de partícula a 0.075 mm, procedimiento basado en la Norma Mexicana NMX-AA-132-SCFI-2006 (SE 2017).
Análisis elemental
La determinación del contenido de elementos se fecetuó con espectroscopia de fluorescencia de rayos X (FRX), usando un espectrómetro portátil FRX Genius 9000, marca Skyray. Se colocaron 3 g de material tamizado dentro de un portamuestras cilíndrico hueco con ventana de mylar (portamuestras espacial para el análisis). Las mediciones se realizaron por triplicado y se reportó el valor promedio como el dato de concentración. El espectrómetro se configuró para determinar la concentración de Cr, Cu, Pb, Ni, V y Zn en suelo, con 30 s de exposición, voltaje de 45 keV, corriente de 40 μA y temperatura de 32.7 ºC. Las curvas de calibración fueron construidas con cuatro estándares (IGLs-1, IGLa-1, IGLsy-2 y IGLgb-3) del Instituto de Geología de la Universidad Nacional Autónoma de México (Lozano y Bernal 2005).
Análisis de propiedades magnéticas
Cada muestra tamizada se colocó dentro de un contenedor cúbico de acrílico de 8 cm3, el cual no interfirió con las mediciones magnéticas. Todas las muestras fueron pesadas en una balanza con susceptibilidad de 0.01 g. Seguidamente, se midió la susceptibilidad magnética a baja (κlf: 0.7 kHz) y alta frecuencia (κhf: 4.7 kHz) con un equipo Bartington MS3 con sensor dual MS2B a temperatura ambiente (22 ºC). Estos valores se normalizaron con la densidad (ρ en kg/m3) para determinar los parámetros magnéticos de susceptibilidad específica (χlf = κlf/ρ) y el porcentaje de susceptibilidad dependiente de la frecuencia (χdf % = ([κlf - κhf]/κlf) × 100), parámetros que proporcionan información sobre el contenido de material magnético y la presencia de partículas finas superparamagnéticas (< 30 nm), respectivamente (Dearing 1999, Evans y Heller 2013). Cada muestra se midió por triplicado, con un error < 1 %.
La magnetización remanente isotérmica de saturación (MRIS) se obtuvo mediante la aplicación de un pulso de campo magnético de corriente continua (CC) de 1000 mT con un magnetizador de pulso IM-10-30, ASC Scientific. La magnetización remanente se midió con un magnetómetro de giro JR6 Agico. Posteriormente, se aplicó un campo magnético de CC invertido de 300 mT y se midió la magnetización de remanencia (IRM-300). Con esta medición, se determinó el parámetro S-300 (S-300 = IRM300mT /MRIS), el cual permite identificar la presencia de materiales de baja coercitividad (S-300 ≥ 0.7) o ferrimagnéticos (MFM), y de alta coercitividad (S-300 < 0.7) o antiferromagnéticos (MFAM; Thompson y Oldfield 1986). Todas las mediciones se realizaron en el Laboratorio Universitario de Geofísica Ambiental de la Universidad Nacional Autónoma de México.
Análisis estadístico
El análisis estadístico se hizo con Statgraphics Centurion XVI y el programa R. El valor de CFT se obtuvo de los datos de concentración clasificados por uso de suelo de tipo residencial (la distribución de datos de suelo de área verde está sesgada al suroeste de la ZMG) y usando el criterio de la Norma Oficial Mexicana NOM-147-SEMARNAT/SSA1-2004, que considera concentración de fondo la media de la población cuando se presenta distribución normal, o la mediana para una distribución diferente (SEMARNAT 2007). A cada grupo de uso de suelo se le aplicó una prueba de normalidad Kolmogorov-Smirnov (KS) con corrección de Lilliefors para comprobación de supuestos estadísticos (linealidad, normalidad, homocedasticidad e independencia). Para encontrar diferencias entre grupos, se usaron una prueba no paramétrica (prueba de Kruskal-Wallis) y una prueba ad-hoc (prueba de Dunn).
Evaluación del enriquecimiento de elementos
El factor de enriquecimiento (FE) se usó para evaluar el incremento de concentración de cada elemento en el suelo urbano, y se calculó de acuerdo con McKenzie et al. (2008):
donde Cm i es el valor de la concentración del elemento en la muestra i y C FT representa la concentración de fondo total del elemento. Se empleó la escala propuesta por Sutherland (2000) para la evaluación: FE < 2: enriquecimiento bajo; FE ϵ [2, 5): enriquecimiento moderado; FE ϵ [5, 20): enriquecimiento alto y FE > 20: enriquecimiento muy alto.
El índice de carga de contaminación (PLI) de Tomlinson et al. (1980), que permite evaluar la contaminación de un sitio como contaminado o no contaminado, se calcula en función de los FE como:
donde PLI es igual a la raíz n del producto de los FE analizados. La interpretación del índice según Tomlinson et al. (1980) indica que si PLI ≤ 1 no hay contaminación, pero si PLI > 1 hay una degradación del sitio o está contaminado. En este estudio se amplía la clasificación y se consideran tres niveles para valores de PLI > 1, con el fin de evaluar mejor el incremento de elementos en el suelo: si PLl ϵ (1, 2]: moderadamente contaminado; si PLI ϵ (2, 3]: altamente contaminado y Si PLI > 3: muy altamente contaminado.
Análisis espacial
Para determinar un modelo de distribución espacial estacionario con predictor de estimación lineal insesgado se utilizó un conjunto de datos normalizados. El análisis geoestadístico consta de tres etapas: análisis exploratorio de datos (AEO), análisis estructural (AE) y un modelo de estimación o interpolación (en este estudio se usó un modelo de interpolación de tipo Kriging ordinario [KO]). En el AEO se verifican los parámetros de tendencia (media, mediana y coeficiente de asimetría), se identifican los valores atípicos y se evalúa la normalidad de la variable. Cuando no se cumplen con los supuestos de normalidad se procede a un cambio de variable de tipo log, inverso u otro, y se repite el AEO. Las variables que no cumplen con los criterios de normalidad son rechazadas para el análisis geoestadístico.
En el AE, se determina la correlación espacial que hay entre los datos de una variable para el área específica y el grado de influencia que existe entre los datos cercanos. Se determina el semivariograma, que es una media aritmética de la diferencia cuadrática de la variable Z mediada en un punto y a una distancia determinada (h), evaluada como:
donde γ(h) representa el valor de la función de semivarianza experimental, N(h) los pares de observación entre ubicaciones a una distancia h, Z(X i ) el valor medido de la variable Z en X i y Z(X i + h) el valor de la variable X i a una distancia h. Los valores de X son las posiciones georreferenciadas de los sitios de muestreo (Webster y Oliver 1990, Hernández-Stefanoni y Ponce-Hernández 2006)
Los parámetros estructurales se determinan a partir del modelo de semivariograma teórico que mejor se ajusta al variograma experimental. El valor de Nugget o varianza pepita (Co) es intercepto con el eje de la semivariograma; representa la variabilidad de las variables a distancias que no presentan relación con la distancia y designa la varianza no explicada. El sill o meseta (Co + C) es el límite del semivariograma cuando la varianza estructural alcanza su valor máximo; se puede presentar como meseta finita o no finita. El rango (Ao) corresponde a la distancia a partir de la cual dos observaciones son independientes, donde las variables pierden su relación espacial. La varianza estructural relativa (VE = C/[Co + C]) se define como la proporción de la varianza no explicada (nugget) entre la varianza estructural (sill). Cuando el modelo se ajusta a los datos de semivariograma teórico, la suma de los residuales cuadráticos (RSS) se reduce y los coeficientes de determinación (R2) se incrementan (Isaaks y Srivastava 1989, Webster y Oliver 1990, Delgado et al. 2010). En la tercera etapa, estiman los datos mediante interpolación Kriging ordinaria, una técnica que proporciona el mejor estimador lineal insesgado.
La comparación entre las capas o modelos de interpolación se efectuó mediante un análisis de geoprocesamiento (álgebra de mapas) para identificar las zonas que presentan una vinculación espacial entre parámetros magnéticos y concentración de elementos (Tomlin 1994). Los modelos y geoprocesamientos se obtuvieron con el software GS+ 9.0 (GD 2015) y la representación espacial fue editada en arcMAP de ArcGIS (ESRI 2004).
RESULTADOS Y DISCUSIÓN
Parámetros magnéticos
Los resultados de los parámetros magnéticos, clasificados por uso de suelo, se presentan en el cuadro I. Los parámetros de dispersión (rango, desviación estándar) indican fluctuaciones entre los datos de los parámetros magnéticos de χlf, χdf% y MRIS, que revelan que la concentración de material magnético no es homogénea en el suelo de la ZMG (Cuadro I).
CUADRO I RESUMEN ESTADÍSTICO DE PARÁMETROS MAGNÉTICOS DE 140 MUESTRAS DE SUELO DE LA ZONA METROPOLITANA DE GUADALAJARA.|
| Parámetro | Uso de suelo | ||||||
| Industrial | Comercial | Equipamiento | Residencial | Área verde | Total | ||
| (20) | (15) | (16) | (78) | (11) | (140) | ||
| χlf | Rango | 0.19-4.22 | 0.6-3.34 | 0.54-3.29 | 0.18-4.8 | 0.52-2.37 | 0.18-4.8 |
| (µm3/kg) | Prom ± SD | 1.71 ± 0.98 | 1.37 ± 0.68 | 1.11 ± 0.72 | 1.09 ± 0.71 | 1.01 ± 0.64 | 1.21 ± 0.77 |
| Mediana | 1.39 | 1.32 | 0.83 | 0.90 | 0.84 | 0.95 | |
| χdf% | Rango | 0.3-4.4 | 0.5-5.2 | 1.2-4.4 | 0.6 - 6.0 | 1.3-4.3 | 0.3-6.0 |
| (%) | Prom ± SD | 2 ± 1.1 | 2.7 ± 1.3 | 2.6 ± 0.9 | 2.5 ± 1.0 | 3.1 ± 0.9 | 2.5 ± 1.1 |
| Mediana | 1.9 | 2.7 | 2.4 | 2.5 | 3.4 | 2.5 | |
| MRIS | Rango | 4.4-78.01 | 5.38-32.38 | 7.93-48.95 | 4.32-55.34 | 5.97-40.01 | 4.32-78.01 |
| mA m2/kg | Prom ± SD | 21.72 ± 16.22 | 16.92 ± 8.86 | 17.77 ± 12.87 | 15.81 ± 10.8 | 15.17 ± 11.25 | 16.95 ± 11.82 |
| Mediana | 17.81 | 15.65 | 12.49 | 12.78 | 9.44 | 13.14 | |
| S-300 | Rango | 0.28-1.00 | 0.87-0.99 | 0.85-1.00 | 0.42-1.00 | 0.89-0.95 | 0.28-1.00 |
| AD | Prom ± SD | 0.89 ± 0.15 | 0.93 ± 0.03 | 0.92 ± 0.04 | 0.93 ± 0.13 | 0.92 ± 0.02 | 0.92 ± 0.12 |
| Mediana | 0.92 | 0.93 | 0.91 | 0.92 | 0.91 | 0.92 | |
χlf: susceptibilidad específica; χdf%: porcentaje de susceptibilidad dependiente de la frecuencia; MRIS: magnetización remanente isotermal de saturación: Prom: promedio; SD: desviación estándar; AD: adimensional.
Los valores de χlf en varios suelos de tipo industrial y comercial son altos (mayores que 1.39 y 1.32 µm3/kg, respectivamente), en comparación con los suelos de tipo equipamiento, residencial y de área verde (< 0.9 µm3/kg), indicando un aumento importante de material magnético para este tipo de suelos, diferencia que se comprobó posteriormente mediante la prueba de Kruskal-Wallis (KW) (Cuadro I). Se observaron valores atípicos en suelos de áreas verdes, comerciales, de equipamiento y residenciales, que sugieren un enriquecimiento de material magnético por una o varias fuentes antrópicas (Cuadro I, Fig. 2).
Fig. 2 Gráfico box-plot de valores de los parámetros magnéticos: susceptibilidad específica (χlf), porcentaje de susceptibilidad magnética dependiente de la frecuencia (χdf%), magnetización remanente isotermal de saturación (MRIS) y cociente S-300 de muestras de suelo superficial de la Zona Metropolitana de Guadalajara, agrupados por uso de suelo. Las barras indican la variación de los datos fuera de los cuartiles superior e inferior, en tanto que los valores atípicos se muestran en círculos negros fuera de las barras. La línea roja punteada en χdf% indica el límite para las muestras que contienen concentraciones bajas de partículas superparamagnéticas (< 2 %; Dearing 1999). La línea azul punteada en S-300, indica el límite que clasificar muestras con mayor contenido de minerales ferrimagnéticos (S-300 > 0.7) o antiferromagnéticos (S-300 < 0.7; Aguilar et al. 2013).
Para identificar las diferencias estadísticas entre los grupos de suelo y, en vista de que no se presenta normalidad entre los datos, se procedió con un análisis no paramétrico (prueba KW). Para los valores de χlf, el análisis KW mostró un valor de p < 0.05 (p = 0.006), indicando una diferencia estadísticamente significativa entre los grupos, con un nivel del 95 % de confianza. La prueba de Dunn (análisis post-hoc) reveló que las diferencias se presentan en los grupos de suelo industrial-área verde (valor p = 0.03) e industrial-residencial (p = 0.01; Fig. 2).
El registro de χdf% mostró valores de 2 a 3 %, indicando una baja concentración de minerales ferrimagnéticos ultrafinos (< 50 nm) o superparamagnéticos (SP) en el suelo de la ZMG, con presencia de granos gruesos magnéticos, de acuerdo con Dearing (1999). Los suelos del área verde presentan mayor concentración de partículas SP, valor esperado por los procesos de intemperismo o bioquímicos que se producen en el suelo (presencia de bacterias magnetotácticas), aunque la adición de partículas finas suspendidas que provienen de fuentes antrópicas (partículas ultra finas < 1 µm) puede mostrar valores altos de χdf% (Dearing 1999, Mistry et al. 2024). Se observaron valores atípicos de χdf% en suelos de tipo residencial (Cuadro I). El análisis KW para χdf% mostró un valor de p > 0.05 (p = 0.07), por lo que no existe una diferencia estadísticamente significativamente entre los grupos de suelo para este parámetro (Fig. 2).
Los resultados de la MRIS (concentración de material magnético de baja coercitividad o ferrimagnéticos [MFM]) indicaron varios suelos de tipo comercial e industrial con alto contenido de MFM (> 17.81 mA m2/kg), en comparación con los suelos del área verde y residencial (< 2.5 mA m2/kg). El aumento de la MRIS en el suelo urbano está relacionado con la adición de partículas de material procedente de fuentes antrópicas, debido a que se observa una mayor concentración de MFM en su composición que en las fuentes naturales (Evans y Heller 2013). Se observó la presencia de valores atípicos en los suelos de área verde, residencial, de equipamiento e industrial (MRIS > 32 mA m2/kg), los cuales exhibieron un enriquecimiento de material magnético de baja coercitividad importante en algunos sitios de la ZMG (Cuadro I). El análisis KW para MRIS mostró un valor de p > 0.05 (p = 0.14), por lo que no existe una diferencia estadísticamente significativamente entre los grupos de suelo (Fig. 2).
Los suelos de la ZMG mostraron valores de S-300 entre 0.70 y 1.00 que revelaron la presencia de minerales de baja coercitividad (ferrimagnéticos), responsables de la señal magnética, cuyo origen es probablemente antrópico. Se reportan dos muestras de suelo (residencial e industrial) con valores de S-300 < 7.00, valor indicativo de suelos que presentan un alto grado de intemperismo y baja perturbación antrópica debido a una elevada concentración de minerales antiferromagneticos (MAFM) o de alta coercitividad (Evans y Heller 2003; Cuadro I, Fig. 2). El análisis KW para S-300, mostró un valor de p > 0.05 (p = 0.82); por lo tanto, no existe una diferencia estadísticamente significativamente entre los grupos de suelo (Fig. 2).
Comparando el valor de la mediana de χlf y MRIS de este estudio con los resultados reportados de suelos y polvo urbanos (PU) de ciudades mexicanas, se encontró que la χlf de los suelos urbanos de la ZMG (0.95 µm3/kg) fue similar a la de suelos de Mérida (Aguilar et al. 2011) y el PU de Mexicali (Sánchez-Duque et al. 2015), y menor para suelos de Morelia (Aguilar et al. 2013) y Ciudad de México (Cejudo et al. 2015a), y PU de Ciudad de México (Cejudo et al. 2015a, 2022), Villahermosa (Cejudo et al. 2017) y Ensenada (Cortés et al. 2017). Esta comparación sugiere que el contenido de material magnético varía entre las ciudades de México, en que la combinación de factores geológicos y antrópicos dan características particulares de acumulación (Fig. 3 y Cuadro II).
Fig. 3 Valores promedios de susceptibilidad específica (χlf) y magnetización remanente isotermal de saturación (MRIS) en suelos y polvo, reportados en diferentes ciudades de México: Guadalajara (GUAD; este estudio), Mérida (MERI), Ciudad de México (CDMX), Morelia (MOR), Mexicali (MEXI), Ensenada (ENSE) y Villahermosa (VH). En la figura se representa el valor medio y su desviación estándar para cada parámetro magnético del suelo de la Zona Metropolitana de Guadalajara con líneas de colores azules y verdes, respectivamente. (CDMXa y CDMXb: resultados reportados en 2015; CDMX: resultados reportados en 2022.)
CUADRO II RESUMEN ESTADÍSTICO DE PARÁMETROS MAGNÉTICOS REPORTADO EN VARIAS CIUDADES DE MÉXICO.
| Parámetro | Ciudad | Tipo | n | Prom | Min | Max | Referencia |
| χlf µm3/kg | ZMG | Suelo | 140 | 1.2 | 0.2 | 4.8 | Este estudio |
| Mérida | Suelo | 90 | 1.0 | 0.0 | 2.0 | Aguilar et al. 2011 | |
| Morelia | Suelo | 98 | 5.2 | 0.3 | 10.0 | Aguilar et al. 2013 | |
| CDMXa | Suelo | 88 | 2.7 | 2.0 | 4.6 | Cejudo et al. 2015a | |
| CDMXb | Polvo | 89 | 4.9 | 3.0 | 10.2 | Cejudo et al. 2015b | |
| CDMX | Polvo | 374 | 3.3 | 1.5 | 6.5 | Cejudo et al. 2022 | |
| Mexicali | Polvo | 108 | 1.4 | 0.2 | 6.9 | Sánchez-Duque et al. 2015 | |
| Ensenada | Polvo | 86 | 4.3 | 1.1 | 10.1 | Cortés et al. 2017 | |
| Villahermosa | Polvo | 82 | 4.8 | 3.3 | 6.2 | Cejudo et al. 2017 | |
| CDMX* | Suelo | 135 | 3.38 | 0.49 | 1.21 | Morton-Bermea et al. 1999* | |
| San Luis Potosí** | Polvo | 100 | 0.0003 | 0.000001 | 0.0034 | Aguilar et al. 2020** | |
| MRIS mA m2/kg | ZMG | suelo | 140 | 16.95 | 4.32 | 78.1 | Este estudio |
| Mérida | Suelo | 90 | 4.00 | 0.00 | 8.00 | Aguilar et al. 2011 | |
| Morelia | Suelo | 98 | 45.00 | 0.00 | 90.00 | Aguilar et al. 2013 | |
| CDMX a | Suelo | 88 | 41.10 | 19.90 | 81.10 | Cejudo et al. 2015a | |
| CDMX b | Polvo | 89 | 62.90 | 33.20 | 144.90 | Cejudo et al. 2015b | |
| CDMX | Polvo | 374 | 44.28 | 16.40 | 81.43 | Cejudo et al. 2022 | |
| Mexicali | Polvo | 108 | 12.51 | 1.61 | 87.27 | Sánchez-Duque et al. 2015 | |
| Villahermosa | Polvo | 82 | 28.30 | 14.30 | 42.30 | Cejudo et al. 2017 | |
| Ensenada | Polvo | 86 | 31.50 | 7.70 | 95.00 | Cortés et al. 2017 | |
| San Luis Potosí** | Polvo | 100 | 0.002 | 0.0006 | 0.04 | Aguilar et al. 2020** |
n: número de muestras; Prom: promedio; Min: mínimo; Max: máximo; χlf: susceptibilidad específica; MRIS: magnetización remanente isotermal de saturación; ZMG: Zona Metropolitana de Guadalajara; CDMX: Ciudad de México.
*Susceptibilidad reportada en valores × 10-8 de acuerdo con el Sistema Internacional de Unidades; **valores no comparables; CDMXa y CDMXb, valores reportados en 2015.
Los valores de MRIS en los suelos de la ZMG son mayores que en suelos de Mérida (Aguilar et al 2011), los cuales exhibieron mayor contenido de material ferromagnético. Sin embargo, el factor geológico puede ser un indicador de esta diferencia. En Mérida, la mayor parte de las rocas son calizas (bajo contenido de MFM), a diferencia de la ZMG, que está asentada sobre roca volcánica (rica en MFM). Por otra parte, la comparación de valores de MRIS en la ZMG con suelos y polvo de Morelia (Aguilar et al. 2013), Ciudad de México (Cejudo et al. 2015a, b, 2022) y Ensenada (Cortés et al. 2017), indica que éstos son menores, aun cuando las condiciones geológicas son similares (Cuadro II). Por lo tanto, el aporte de minerales ferrimagnéticos en las ciudades se debe a otro tipo de fuentes (antrópicas o ambientales) que provocan variaciones en el contenido de los minerales magnéticos en el suelo.
La comparación de valores medios de MRIS en suelos de la ZMG, y PU en Mexicali (Sánchez-Duque et al. 2015) y Villahermosa (Cejudo et al. 2017), indica que no hay diferencias significativas. Sin embargo, esta igualdad puede ser sólo coincidente, ya que las matrices de las muestras provienen de condiciones geológicas y ambientales diferentes, por lo que es necesario hacer más estudios para verificar este resultado.
Los parámetros magnéticos reportados en polvo y suelo en diferentes ciudades en México presentan una gran variabilidad, lo que dificulta establecer un valor estándar para determinar la concentración de material magnético como referencia o contaminante (Cuadro II y Fig. 3). Sin embargo, existen propuestas en biomonitores para generar modelos basados en aprendizaje automatizado, que podrían utilizarse tanto en suelo como en polvo para estandarizar el uso de parámetros magnéticos como proxys de contaminación (Salazar-Rojas et al. 2023, Chaparro et al. 2024). Hay otros trabajos que reportan datos magnéticos para la Ciudad de México (Morton-Bermea et al. 2009) y San Luis Potosí (Aguilera et al. 2020). No obstante, una falta de estandarización en el uso de unidades dificulta su comparación con la literatura existente (Cuadro II).
Concentración de elementos
Los resultados de la concentración de elementos (Cr, Cu, Pb, V y Zn) y el factor de enriquecimiento de los usos de suelo en la ZMG se presentan en el cuadro III. Los valores de concentración de fondo total (CFT) se obtuvieron del suelo de tipo residencial (58 muestras), debido a su amplia distribución en la ZMG y con base en el criterio de la Norma Mexicana NMX-AA-132-SCFI-2006 (SE 2017). Los valores de Cr, Cu, Ni, Pb, V y Zn, fueron 13, 18, 10, 26, 36 y 269 mg/kg, respectivamente (Cuadro III).
CUADRO III RESUMEN ESTADÍSTICO DE CONCENTRACIÓN DE ELEMENTOS (mg/kg) Y PLI EN MUESTRAS DE SUELO DE LA ZONA METROPOLITANA DE GUADALAJARA.
| Elemento | n | Uso de suelo | Total | FE | |||||
| Industrial | Comercial | Equipamiento | Residencial | Área verde | |||||
| Cr | Rango | 10-41 | 7-32 | 10-55 | 7-63 | 8-22 | 7-63 | 0.5-4.8 | |
| 106 | Prom ± SD | 19 ± 9 | 17 ± 8 | 19 ± 13 | 16 ± 9 | 15 ± 6 | 17 ± 9 | 1.2 ± 0.4 | |
| Mediana | 16 | 18 | 14 | 13 (58) | 14 | 15 | 1.1 | ||
| Cu | Rango | 11-682 | 7-80 | 7-166 | 7-159 | 7-28 | 7-682 | 0.4-38 | |
| 110 | Prom ± SD | 126 ± 185 | 36 ± 23 | 32 ± 44 | 30 ± 31 | 18 ± 7 | 48 ± 90 | 2.6 ± 5.0 | |
| Mediana | 44 | 34 | 18 | 18 (58) | 19 | 22 | 1.2 | ||
| Ni | Rango | 7-20 | 8-15 | 7-46 | 7-22 | 7-10 | 7-46 | 0.7-4.6 | |
| 103 | Prom ± SD | 12 ± 4 | 11 ± 3 | 14 ± 10 | 11 ± 3 | 9 ± 2 | 11 ± 5 | 1.1 ± 0.5 | |
| Mediana | 12 | 10 | 10 | 10 (58) | 9 | 10 | 1.0 | ||
| Pb | Rango | 15-531 | 14-292 | 14-330 | 8-305 | 11-59 | 8-531 | 0.3-20.5 | |
| 138 | Prom ± SD | 94 ± 124 | 60 ± 72 | 68 ± 89 | 46 ± 56 | 31 ± 14 | 56 ± 75 | 2.2 ± 2.9 | |
| Mediana | 46 | 37 | 28 | 26 (76) | 31 | 28 | 1.1 | ||
| V | Rango | 10-86 | 9-51 | 15-68 | 8-123 | 11-11 | 8-123 | 0.2-3.4 | |
| 42 | Prom ± SD | 47 ± 27 | 36 ± 19 | 29 ± 21 | 43 ± 28 | 11 ± 0 | 41 ± 26 | 0.3 ± 0.7 | |
| Mediana | 48 | 43 | 19 | 36 (19) | 11 | 35 | 1.0 | ||
| Zn | Rango | 178-1768 | 140-716 | 123-1348 | 79-1215 | 109-316 | 79-1768 | 0.3-6.6 | |
| 140 | Prom ± SD | 641 ± 364 | 363 ± 161 | 337 ± 301 | 325 ± 213 | 220 ± 75 | 367 ± 264 | 1.4 ± 1.0 | |
| Mediana | 556 | 378 | 244 | 269 (78) | 203 | 283 | 1.1 | ||
| PLI | Rango | 0.8-5.2 | 0.8-2.0 | 0.4-5.3 | 0.5-3.4 | 0.6-1.3 | 0.4-5.3 | ||
| 140 | Prom ± SD | 2.0 ± 1.1 | 1.3 ± 0.3 | 1.3 ± 1.1 | 1.2 ± 0.6 | 0.9 ± 0.2 | 1.3 ± 0.8 | ||
| Mediana | 1.7 | 1.3 | 1.0 | 1.0 | 0.9 | 1.0 | |||
n: número de muestras; Prom: promedio; SD: desviación estándar; FE: factor de enriquecimiento; PLI: índice de carga contaminante. Las cifras en negritas representan los valores de concentración de fondo total propuestos para la Zona Metropolitana de Guadalajara.
Las pruebas de normalidad (Kolmogorov-Smirnov con corrección de Lilliefors) no se satisficieron para los resultados de la concentración de elementos (Cr, Cu, Ni, Pb, V y Zn) clasificados por uso de suelo. La comparación entre los valores de las medianas indicó que el contenido de elementos en suelos de la ZMG disminuye en el siguiente orden: Zn > V > Pb > Cu > Cr > Ni (Cuadro III).
Contrastando los valores de las medianas para Ni, Pb y V (14, 10, 27 mg/kg, respectivamente) con la concentración de referencia total (CRT) de suelos contaminados para uso agrícola/residencial/comercial de la Norma Oficial Mexicana NOM-147-SEMARNAT (Ni, Pb y V: 1600, 400 y 78 mg/kg, respectivamente), se observó que, la mayor parte de suelos de la ZMG se clasifican como no contaminados (SEMARNAT 2007). Sin embargo, se observó que algunos suelos de tipo industrial y residencial presentan valores superiores a los valores de CRT para Pb y V que requieren de atención (Cuadro III).
Los FE obtenidos con las CFT de este trabajo indicaron un incremento en la concentración de varios elementos, principalmente Cu, Pb y Zn, cuyos valores de FE estuvieron entre 5 y 20 y, fueron categorizados como enriquecimiento medio, alto y muy alto. Se observaron muestras de suelos con FE > 5 para Pb, en suelo de tipo residencial, comercial y equipamiento, lo que representa un peligro para las personas que habitan cerca de estos lugares (Fig. 4 y Cuadro III).
Fig. 4 Diagramas de caja para factores de enriquecimiento (FE) de Cr, Cu, Ni, Pb, V y Zn (eje x) para los distintos usos de suelo (eje y) en la Zona Metropolitana de Guadalajara. Las líneas punteadas representan los valores de la escala propuesta por Sutherland (2000) para evaluar el factor de enriquecimiento (FE < 2: sin enriquecimiento; FE ϵ [2,5): enriquecimiento moderado; FE ϵ [5,20): alto enriquecimiento; FE > 20: muy alto enriquecimiento). Las barras indican la variación de los datos fuera de los cuartiles superior e inferior, en tanto que los valores atípicos se muestran en círculos negros fuera de las barras.
Las pruebas de KW y Dunn mostraron que los FE de Cu, Pb y Zn presentan diferencias entre grupos en el siguiente orden: industrial > comercial-equipamiento > residencial-área verde (Fig. 4); por lo tanto, suelos de uso industrial-comercial-equipamiento requieren mayor vigilancia sobre el monitoreo de elementos.
La estimación de valores de PLI en los suelos de la ZMG mostró valores de 0.4 a 5.3; por lo tanto, hay suelos considerados no contaminados, con contaminación moderada a muy altamente contaminados (residencial, equipamiento e industrial; Cuadro III).
Correlaciones de los FE de los elementos, índice de carga contaminante y propiedades magnéticas
Se hizo un análisis exploratorio de datos y aplicación de transformación logarítmica a los datos de χlf, MRIS y FE de Cr, Cu, Pb, Zn y PLI, los cuales satisfacen la prueba de normalidad de Kolmogorov-Smirnov con corrección de Lilliefors, con valores de p ≥ 0.05 para obtención de correlaciones significativas (Cuadro IV). Los valores de V y Ni no cumplen con las pruebas de normalidad y se excluyeron de los análisis.
CUADRO IV CORRELACIONES ENTRE FACTORES DE ENRIQUECIMIENTO DE LA CONCENTRACIÓN DE ELEMENTOS, ÍNDICE DE CARGA CONTAMINANTE Y PARÁMETROS MAGNÉTICOS DE SUELOS URBANOS DE LA ZONA METROPOLITANA DE GUADALAJARA.
| Log(FE Cr) | Log(FE Cu) | Log(FE Pb) | Log(FE Zn) | Log(PLI) | Log(χlf) | Log(MRIS) | |
| Log(FE Cr) | |||||||
| Log(FE Cu | 0.25 | ||||||
| Log(FE Pb) | 0.39 | 0.48 | |||||
| Log(FE Zn) | 0.26 | 0.68 | 0.62 | ||||
| Log(PLI) | 0.56 | 0.82 | 0.80 | 0.80 | |||
| Log(χlf) | 0.37 | 0.56 | 0.53 | 0.48 | 0.70 | ||
| Log(MRIS) | 0.25 | 0.50 | 0.43 | 0.36 | 0.54 | 0.80 |
Log: logaritmo; PLI: índice de carga contaminante; χlf: susceptibilidad específica; MRIS: magnetización remanente isotermal de saturación.
Las cifras en negritas indican correlación valida con valor p < 0.05.
Las variables transformadas de los parámetros magnéticos (χlf, MRIS) con los factores de enriquecimiento de Cr, Cu, Pb, Zn y PLI mostraron coeficiente de correlación bajo (R2 ϵ (0.3, 0.5]) y medio (R2 ϵ (0.5,0.7]; Cuadro IV). La transformación de χlf mostró una baja correlación con los factores de enriquecimiento de Cr y Zn, y media con Cu, Pb y PLI, en tanto que la trasformación de MRIS exhibió una baja correlación con Cu, Pb y Zn, y media con PLI (Cuadro IV).
Dada la correlación observada entre los parámetros magnéticos con el PLI, se proponen dos modelos, para χlf: (Log(PLI) = 0.94 + 0.89 × χlf, R2 = 0.57); y MRIS: Log(PLI) = -1.0 + 0.42 × Log(MRIS, R2 = 0.35) como proxy para la evaluación o monitoreo de la contaminación en el suelo de la ZMG. Estos modelos satisfacen los supuestos de linealidad (media de los residuos: -8.8 × 10-8 y 3.1 × 10-18), normalidad (p = 0.09 y 0.59), homocedasticidad (p = 0.35, 0.14) e independencia (estadístico Duvin-Watson: 1.75, 1.21).
Los valores equivalentes de χlf y PLI para categorizar el incremento de elementos fueron los siguientes: sin contaminación, < 1.20 µm3/kg; contaminación moderada, de 1.20 a 1.69 µm3/kg; contaminación alta, de 1.70 a 2.30 µm3/kg; y contaminación muy alta, > 2.30 µm3/kg. En el caso del MRIS: sitios sin contaminación, < 13.0 mA m2/kg; con contaminación moderada, de 13.0 a 20.0 mA m2/kg; con contaminación alta, de 20.1 a 25.0 mA m2/kg; y con contaminación muy alta > 25.0 mA m2/kg. Los valores de clasificación de sitios contaminados para los parámetros magnéticos en la ZMG son diferentes a los propuestos para la Ciudad de México por Morton-Bermea et al. (2009) y por Cejudo et al. (2022) (χlf: 4.0 µm3/kg; MRIS: 38 mA m2/kg), debido a las diferencias geológicas (tipos de suelos, fisiología y clima) y sociales (variación de la densidad de población y fuentes antrópicas), que aportan cantidades variables de materiales magnéticos al suelo y polvo urbano.
Análisis espacial
Modelos de distribución espacial
Se construyó una base de datos georreferenciados para los parámetros magnéticos (χlf y MRIS), factores de enriquecimiento (FE) de elementos (Cr, Cu, Pb y Zn) e índice de carga contaminante (PLI). Para el análisis geoestadístico (AG) se utilizó el software GS+ (GD 2015). El AG se realizó en tres etapas. La primera consistió en un análisis exploratorio de datos (AEO), cuyo objetivo fue verificar la normalidad de los datos y determinar parámetros de tendencia central, coeficientes de asimetría y valores atípicos. Para este estudio se retiraron los valores atípicos y se usaron los datos trasformados por logaritmo de χlf, MRIS, FE de Cu, Cr, Pb, Zn y PLI. Se excluyeron los datos que no presentaron normalidad (V y Ni). En la segunda etapa, se analizaron la correlación espacial y la influencia de los datos cercanos, se usó una función experimental de semivariograma, y se ajustó un modelo teórico. Los resultados de los semivariogramas para los factores de concentración de elementos, parámetros magnéticos y PLI se presentan en el cuadro V.
CUADRO V SEMIVARIOGRAMAS EXPERIMENTALES PARA LA CONCENTRACIÓN DE ELEMENTOS DEL SUELO DE LA ZONA METROPOLITANA DE GUADALAJARA E ÍNDICE DE CARGA CONTAMINANTE.
| Variable | Modelo | Varianza pepita | Meseta | Rango | Varianza estructural relativa | Modelo |
| Co | Co + C | (m) | C/(C + Co) (%) | R2 | ||
| FE (Cr) | Exponencial | 0.0064 | 0.0675 | 1770 | 91 | 0.93 |
| FE (Cu) | Exponencial | 0.0840 | 0.5050 | 9990 | 83 | 0.90 |
| FE (Ni)* | - | - | - | - | - | - |
| FE (Pb) | Esférico | 0.0010 | 0.5950 | 2530 | 99 | 0.95 |
| FE (V)* | - | - | - | - | - | - |
| FE (Zn) | Exponencial | 0.0157 | 0.1174 | 5310 | 86 | 0.82 |
| χlf | Exponencial | 0.0063 | 0.0461 | 5040 | 86 | 0.85 |
| MRIS | Exponencial | 0.0340 | 0.3450 | 3840 | 90 | 0.84 |
| PLI | Exponencial | 0.0191 | 0.0798 | 6840 | 76 | 0.85 |
χlf: susceptibilidad específica; MRIS: magnetización remanente isotermal de saturación; PLI: índice de carga contaminante.
*Variables que no cumplen con los criterios de normalidad para interpolación espacial.
La evaluación de los modelos de semivariogramas teóricos se hizo con los parámetros del nugget o varianza pepita (Co), sill o meseta (Co+C), rango, varianza estructural y coeficiente de determinación (R2).
Los factores de enriquecimiento de Cr, Cu, Pb, Zn, PLI y los parámetros magnéticos son representados por semivariogramas de tipo esférico y exponencial (Cuadro V), con un valor de R2 de 0.82 a 0.95 y una varianza estructural relativa entre 76 y 99 %. Esto indicó que la varianza no explicada o varianza pepita es pequeña comparada con el valor de la meseta o sill. El rango indica la distancia en que las variables presentan una correlación espacial; en este estudio, este parámetro varió de 1770 a 9990 m (Cuadro V). La distancia de muestreo de este estudio fue de 600 m y el rango obtenido de los semivariogramas tuvo un valor mínimo de correlación de 1770 m, lo cual sugiere un muestreo representativo para las variables medidas en suelos de la ZMG, indicando la viabilidad de usar la interpolación KO para Cr, Cu, Pb, Zn, χlf, MRIS y PLI (Webster y Oliver 1990, Delgado et al. 2010).
Los modelos de distribución espacial permitieron identificar las zonas de acumulación para los elementos de Cr, Cu, Pb y Zn, la concentración de material magnético (χlf y MRIS) y las zonas que requieren mayor vigilancia o estudios complementarios para determinar el riesgo potencial debido a la acumulación conjunta de elementos (PLI).
Los mapas de distribución del factor de enriquecimiento de Cr y Zn mostraron que la mayor parte de los suelos urbanos de la ZMG mantiene una baja acumulación de Cr y Zn. En la parte central de la ZMG hay una zona con factor de enriquecimiento medio, muy próxima a una zona industrial, donde el tránsito vehicular es intenso. Estas actividades pueden ser las responsables del enriquecimiento del suelo urbano con Zn, ya que este elemento está asociado con emisiones vehiculares (Huang et al. 1994, Wang et al. 2021; Fig. 5)
Fig. 5 Distribución de los factores de enriquecimiento de Cr, Cu, Pb y Zn en suelos de la Zona Metropolitana de Guadalajara. (GDL: Gobierno de Guadalajara.)
Los mapas de distribución del FE de Cu y Pb mostraron tres zonas: baja, media y alta. En ambos mapas se observó un gradiente de enriquecimiento alto en el centro de la ciudad que disminuye radialmente hacia la zona periurbana. Las zonas que presentan enriquecimiento medio a alto se expanden sobre polígonos que tienen población alta (> 10000 habitantes, Fig. 1) en la parte central de la ZMG, por lo que ésta se identificó como área de riesgo (Figs. 1 y 5). Las zonas con factor de enriquecimiento alto se encuentran próximas a una zona industrial y un patio de maniobras ferroviario, condición que puede incrementar el contenido de Cu y Pb (Fig. 5).
El Pb es un elemento altamente tóxico para el ser humano. Este estudio evidencia varias zonas de enriquecimiento alto en suelos de la ZMG, por lo que es necesaria una evaluación mayor para identificar y delimitar el área de potencial riesgo para la salud humana (Fig. 5).
Los modelos de distribución espacial de los parámetros magnéticos obtenidos en suelo urbano de la ZMG mostraron un gradiente de concentración radial del contenido de material magnéticos (Fig. 5).
Los mapas de χlf y MRIS mostraron que la zona central (5 % del área de la ZMG) contiene una alta concentración de material magnético, en tanto que el 80 % del área de la ZMG exhibió una concentración moderada de material magnético (Fig. 5). Se observó que la ubicación de industrias, áreas de maniobras del ferrocarril y aeropuertos muestran una relación con la adición de material magnético en el suelo, ya que en la proximidad de estos sitios (~ 5 km) la concentración de material magnético fue moderada a alta (Fig. 6).
Fig. 6 Distribución espacial de parámetros magnéticos e índice de carga contaminante (PLI, por su sigla en inglés) en suelos de la Zona Metropolitana de Guadalajara. (GDL: Gobierno de Guadalajara.)
Las áreas de la ZMG con mayor concentración de elementos se determinaron con el modelo de distribución del PLI (Fig. 6). Los suelos de la zona centro fueron los más contaminados (PLI > 3). Las zonas de contaminación moderada (PLI ϵ (2,3]) se extienden desde la parte central en dirección noreste-sureste (~ 80 % del área de la ZMG). Se observa que los parques industriales, así como las instalaciones aeroportuarias y del ferrocarril, son un factor influyente (Fig. 6). Las zonas sin contaminación (PLI ≤ 1) se encuentraron en la periferia de la ZMG (Fig. 6).
De acuerdo con censo de población de 2020, la zona con alta contaminación registra una población de casi 300 000 habitantes, que equivale al 6.3 % de la población total de la ZMG (~ 5 millones de habitantes), por lo que hay un riesgo potencial para la población que habita en esta área (INEGI 2022a; Fig. 6).
La comparación de superficies de zonas de clasificación de PLI con los parámetros magnéticos en los mapas de distribución mostraron que hay una correspondencia espacial, que para la χlf fue del 79 % del área. Para la MRIS la correspondencia fue del 71 % del área. Por otra parte, la comparación espacial de las zonas de distribución de Pb con los parámetros magnéticos fue de 75 % para la χlf y de 72 % para la MRIS; para el Cu fue < 60 % para ambos parámetros (Figs. 5 y 6). Por lo tanto, los parámetros magnéticos de χlf y MRIS pueden usarse como proxys de monitoreo para estimar valores de PLI y contenido de Pb en suelos urbanos de la ZMG.
CONCLUSIÓN
Se logró determinar el contenido de elementos en suelos urbanos de la Zona Mteropolitana de Guadalajara (ZMG) y establecer valores de concentración de referencia local para Cr, Cu, Ni, Pb, V y Zn.
El contenido de elementos de los suelos urbanos de Guadalajara disminuye de la siguiente forma: Zn > V > Pb > Cu > Cr > Ni. Se observaron varios sitios en las zonas centro y noroeste de la ciudad con áreas densamente pobladas que exhiben concentraciones altas de Cu, Pb y Zn, por lo cual se requieren mayores estudios y vigilancia en estas áreas, ya que existe un riesgo potencial para la salud de la población.
Hay indicios de que las altas concentraciones de Cu, Pb y Zn están relacionadas con los usos de suelo en el siguiente orden: industrial > comercial-equipamiento > residencial-área verde en la ZMG, siendo necesarios más estudios para identificar las fuentes.
Los suelos de la ZMG presentan concentraciones variables de material magnético, siendo los suelos industrial y comercial los que contienen mayor cantidad de minerales ferrimagnéticos de baja coercitividad.
La susceptibilidad magnética y magnetización remanente isotermal de saturación mostraron una relación con la concentración de elementos de forma espacial y estadística, por lo que es viable su uso como indicadores proxy de contaminación, principalmente para acumulación de Pb e identificación de sitios contaminados por acumulación de elementos.
Los resultados de este estudio ayudarán a mejorar las evaluaciones ambientales para las concentraciones de Cr, Cu, Ni, Pb, V y Zn en suelos de la ZMG y apuntan al uso de los parámetros magnéticos como alternativa de monitoreo en suelos, lo cual permitiría disminuir costos o priorizar áreas de estudio para el beneficio de la población.
