Área de estudio

La zona de estudio se ubicó en los municipios de Irapuato, Silao, Abasolo, Cuerámaro y Salamanca, que corresponden a la parte suroeste del estado de Guanajuato, México (Figura 1). El estado de Guanajuato cuenta con una extensión territorial de 30 608 km² que representa 1.6% del territorio nacional mexicano. De acuerdo a la información del Consejo de Recursos Minerales (1992), el área bajo estudio pertenece a la provincia fisiográfica conocida como Eje Neovolcánico, situada en la parte sur del estado de Guanajuato.

Fuente: Elaboración propia.

Figura 1 Mapa del área de estudio y localización de los sitios de muestreo. Muestra agrícola Salamanca espárrago (MASE); Muestra agrícola Salamanca fresa (MASFE); Muestra agrícola lo de Juárez (MAJ); Muestra agrícola Copal (MAC); Muestra agrícola Abasolo (MAA); Muestra agrícola San Andrés de la Baraña (MASA); Muestra agrícola Silao (MASI); Muestra agrícola Cuerámaro Jitomate (MACJI); Muestra agrícola Cuerámaro con sodio (Na) (MACNA); Muestra agrícola Silao (MASIM); Muestra Irapuato (IRA); Muestra Irapuato 1 (IRP1); Muestra Irapuato 2 (IRP2); Muestra Irapuato 3 (IRP3); Muestra Irapuato 4 (IRP4); Muestra Lo de Juárez (JUA2); Muestra Cuerámaro 3 (CUE3); Muestra Abasolo 5 (ABA5). 

El Eje Neovolcánico se caracteriza por un típico paisaje volcánico conformado por sierras que presentan alturas medias de 1800 m.s.n.m. y extensos valles intermontanos rellenos de sedimentos continentales modernos. Con base en la carta estatal de suelos elaborada por la Secretaría de Programación y Presupuesto (1980), la zona de estudio ubicada al suroeste del estado de Guanajuato presenta predominantemente suelos tipo Vertisoles (Vertisol pélico [Vp]). Los Vertisoles son suelos muy arcillosos (> 30%), con alta proporción de arcillas expansibles (Grupo de las Esmectitas). Según Köppen, modificado por E. ^{García (1973)}, en el área donde se localizan los suelos analizados (Figura 1), predomina el clima templado sub-húmedo con lluvias en verano (Cw).

Toma y preparación de las muestras

En la Figura 1 se ubican geográficamente los puntos de muestreo de los suelos. El mapa de los sitios de muestreo fue elaborado usando ArcMap^TM 10.2.2 de ArcGIS® software propiedad de Esri. El muestreo se realizó de acuerdo a la NOM-021-SEMARNAT-2000 (Diario Oficial de la Federación ^{[DOF], 2002}). La toma de muestras se llevó a cabo durante los meses de enero a mayo de 2014 y 2015, periodo que corresponde a las estaciones de invierno y primavera boreal. El método de muestreo empleado fue el de zig-zag a 20 cm de profundidad. El número de muestras individuales recolectadas para cada una de las muestras compuestas fue de 15 y se abarcó un área de muestreo de 1 ha. Las submuestras se mezclaron y homogeneizaron hasta obtener una muestra representativa de 2 kg. Posteriormente, se procedió al secado de las muestras por una semana a temperatura ambiente, luego a la molienda y, por último, al tamizado con malla N° 10.

Análisis físicos

La caracterización física y química de los suelos se llevó a cabo en el Laboratorio de Tecnologías para Sustentabilidad Ambiental y el Laboratorio de Biotecnología Alimentaria y Vegetal (División de Ciencias de la Vida, Universidad de Guanajuato). La metodología empleada fue la reportada por la norma oficial mexicana NOM-021-SEMARNAT-2000 (^{DOF, 2002}) y la precisión analítica fue avalada por tres repeticiones de cada muestra.

El contenido de humedad (H) de los suelos se determinó mediante el método gravimétrico, calculando la diferencia de peso entre la muestra de suelo húmeda y después de haberse secado en una estufa marca Binder a 105 ºC durante 24 h.

La clase textural de los suelos se calculó por el método de ^{Bouyoucos (1962}), removiendo previamente la materia orgánica de los suelos con peróxido de hidrógeno al 30%. Luego, tres agentes dispersantes se agregaron, oxalato de sodio, metasilicato de sodio y hexametafosfato de sodio, para individualizar las partículas de arcilla y posteriormente se realizaron lecturas de la temperatura y del hidrómetro a 40 s y 2 h. La determinación de la clase textural se realizó con base en el triángulo de texturas del United States Department of Agriculture (USDA).

La densidad aparente (DA) se calculó por el método de la probeta a partir del conocimiento de la masa del suelo y el volumen total medido en probetas de 50 ml y 100 ml. La densidad real (DR) se determinó mediante el método del picnómetro, con la masa y el volumen total de suelo. La masa se determinó pesando directamente el suelo y el volumen de manera indirecta por el cálculo de la masa y la densidad del agua desplazada por la muestra de suelo. La porosidad total (Pt) de las muestras de suelo se calculó con base en la siguiente ecuación:

Análisis químicos

El pH fue determinado por el método potenciométrico en una suspensión suelo-agua en una relación 1:2 usando un potenciómetro Fisher Scientific AB15, previamente calibrado con soluciones reguladoras a pH 4, 7 y 10.

Para la cuantificación del contenido de materia orgánica (MO) se utilizó la técnica de Pérdida por Ignición (PPI₅₅₀; ^{Heiri, Lotter & Lemcke, 2001}), donde las muestras de suelo fueron calcinadas a una temperatura de 550 °C en una mufla marca Felisa durante tres horas y, por combustión húmeda, con la técnica de ^{Walkley & Black (1934)}. Esta última se basa en la oxidación del carbono orgánico (CO) del suelo por medio de una disolución de dicromato de potasio y el calor de reacción que se genera al mezclarla con ácido sulfúrico concentrado. El dicromato de potasio residual se valoró con sulfato ferroso 6-hidratado. Para calcular el porcentaje de MO se utilizó el factor de corrección de 1.724 (Factor de Van Benmelen).

La determinación de nitrógeno total (NT) se realizó siguiendo el método Micro-Kjeldahl. En una primera fase la muestra se digirió por calentamiento con ácido sulfúrico concentrado y una mezcla de sulfato de sodio y sulfato cúprico. En una segunda fase la muestra se expuso a hidróxido de sodio y agua destilada en un destilador marca Büchi modelo K-350. En la última fase se realizó la valoración del amoniaco recuperado en un volumen conocido de solución valorada de ácido bórico con ácido sulfúrico al 0.25 N (^{Ma & Zuazaga, 1942}).

Análisis de elementos traza

El análisis de elementos traza (As, Cd, Cr, Cu, Pb y Zn) fue realizado en el Laboratorio de Geoquímica Ambiental del Centro de Geociencias-Universidad Nacional Autónoma de México (UNAM) (Juriquilla, Querétaro), de acuerdo al método Environmental Protection Agency^{(EPA) (2001)}: “Trace elements in water, solids, and biosolids by inductively coupled plasma atomic emission spectrometry”. La medición se realizó mediante Espectrometría de Emisión Óptica con Plasma Acoplado Inductivamente Thermo iCAP 6500 Dúo, utilizando un patrón interno de itrio de 1 ppm para evaluar que la señal del instrumento se encontrara dentro de un margen de ± 10% de error. Los estándares de calibración fueron preparados a partir de soluciones de referencias certificadas de 100 ppm en siete diferentes concentraciones. Los coeficientes de correlación para las curvas de calibración fueron mayores a 0.9990, los porcentajes de recuperación fueron mayores al 90% y evaluados mediante un Material de Referencia Estándar el cual corresponde al National Institute of Standards & Technology (NIST) de Estados Unidos con la clave 2710a.

Cálculo del I_geo

Se calculó el I_geo para cada muestra de suelo con base en el indicador propuesto por ^{Müller (1969}) mediante la siguiente ecuación:

donde C _n es la concentración medida del elemento en el sedimento o suelo, B _n es el valor del elemento en la corteza continental. Los valores resultantes fueron contrastados con las clases correspondientes (Tabla 1).

Análisis estadísticos

La comparación entre sitios de muestreo para las variables dependientes (H, clase textural, DA, DR, Pt, pH, CO, MO, NT, C/N y PPI₅₅₀) fue realizada mediante un diseño de una vía (prueba de F) o la prueba H de Kruskal-Wallis, dependiendo del cumplimiento de los supuestos para el análisis paramétrico en cada caso. Adicionalmente se obtuvieron los coeficientes de correlación (Pearson o Spearman) para determinar el grado de asociación entre los pares de variables (H, clase textural, DA, DR, Pt, pH, CO, MO, NT, relación C/N, PPI₅₅₀, As, Cd, Cr, Cu, Pb y Zn). Los análisis se ejecutaron con el paquete estadístico Statistical Analysis System (^{SAS, 2001}).

Análisis físicos

Las medianas de las propiedades físicas de las muestras de suelos se presentan en la Tabla 2. A excepción de la DR (prueba F), las otras características (arena, limo, arcilla, H, DA y P_t) fueron analizadas con la prueba H de Kruskal-Wallis, mostrando en todos los casos diferencias altamente significativas; lo cual permitió concluir en cada variable dependiente que la mediana de un sitio de muestreo es diferente a las demás. Enseguida se aplicó una prueba de F a la variable generada con los rangos asignados por esta prueba no paramétrica, confirmando las diferencias significativas ahora con los promedios de los rangos asignados, a los cuales además se les realizó la prueba de Tukey (p = 0.05), observándose que en cada una de estas variables se presentaron alrededor de 13 grupos estadísticos (Tabla 2)

El suelo con mayor porcentaje de H resultó ser la muestra IRP2, seguido de IRP1 y CUE3 (Tabla 2 y Figura 1). Estos 3 suelos tienen la particularidad de presentar los contenidos de arcilla más altos (máximo 83.80%; Tabla 2), lo que les confiere una retención de humedad mayor (^{Jury, Gardner & Gardner, 1991}). Por otra parte, los suelos JUA2, MAJ y MASIM (Figura 1) presentaron contenidos de arena comparativamente más altos (> 45%) en relación con los contenidos de H más bajos (H máxima: 2.22%; Tabla 2). De la misma manera, la matriz de correlación de Spearman (r _s ; Tabla 3) evidenció una correlación fuertemente positiva (r _s = 0.93753, valor-p = < 0.0001) entre el % arcilla y la H, y contrariamente, una correlación negativa entre el % arena y la H (r _s = -0.69659, valor-p = 0.0013; Tabla 3).

La mayoría de los suelos mostraron valores de DA comprendidos entre 1.00 g cm^-3 y 1.19 g cm^-3 (Tabla 2) y según la NOM-021-SEMARNAT-2000 (^{DOF, 2002}), estas DA son comunes en suelos minerales con textura arcilla. Los valores de DA de 1.00 g cm^-3 a 1.19 g cm^-3 indicarían que la mayor parte de los suelos agrícolas analizados no presentan problemas de compactación lo que conlleva a una buena capacidad de infiltración y aireación (^{Sumner, 1999}). La muestra con la DA más alta fue JUA2 (1.32 g cm^-3; Figura 1) acorde con un suelo de clase textural francosa. De la misma manera, la DA mostró una correlación positiva con el % arena (r _s = 0.61348, valor-p = 0.0068) y negativa con el % arcilla (r _s = -0.84967, valor-p = < 0.0001).

Según la clasificación de DR, reportada por ^{Porta, López-Acevedo & Roquero (2003)}, las muestras de suelo indicaron un predominio de minerales de arcilla (DR: 2.00 g cm^-3 - 2.65 g cm^-3) dentro de las partículas inorgánicas del suelo, densidades comunes en suelos con contenidos altos de filosilicatos y tectosilicatos. Los suelos estudiados en este trabajo derivarían de la meteorización de rocas predominantemente basálticas, los cuales se caracterizan por el predominio de estos minerales silicatados ricos en Fe y Mg como componentes sólidos del suelo (^{Tarbuck & Lutgens, 2005}).

Según ^{Oades (1984)}, la estabilidad de los poros asociada a los macro y micro agregados es esencial para el crecimiento óptimo de las plantas donde su estructura se ve favorecida por cationes multivalentes que actúan como puentes entre los coloides orgánicos y arcillas. La dispersión de las partículas de arcilla de los microagregados se promueve por la adsorción de los ácidos orgánicos producidos por plantas, bacterias y hongos que aumentan su carga negativa. La capacidad de dispersión de la arcilla se compensa por la acción unificadora de polisacáridos (principalmente mucílagos producidos por las bacterias), las raíces de las plantas e hifas de los hongos. Los suelos agrícolas con una adecuada actividad biológica generalmente contienen el 60% de espacio poroso y es indicativo del grado de aireación del suelo, retención de humedad y enraizamiento. En este trabajo el valor más bajo encontrado de Pt fue de 44.05% (JUA2: clase franco arenosa) y el más alto de 58.09% (IRP3: clase arcilla; Tabla 2 y Figura 1).

Análisis químicos

En la Tabla 4 se muestran las medias o en su caso las medianas de las variables químicas en los sitios evaluados. A las variables: pH, CO, MO y NT se les aplicó la prueba H de Kruskal-Wallis, y a la relación C/N y PPI₅₅₀ la prueba de F.

De acuerdo a la NOM-021-SEMARNAT-2000 (^{DOF, 2002}), el pH en los diferentes suelos mostró una gran variabilidad de resultados, desde moderadamente ácido hasta fuertemente alcalino (6.37 - 8.69; Tabla 4). Según ^{Porta et al. (2010)}, este parámetro es un buen indicador de la presencia o ausencia de determinados constituyentes en el suelo, de los cationes que habrá en el complejo de cambio, de la disponibilidad de nutrientes para los cultivos y de la actividad biológica en el suelo. En general, se considera que el rango de pH obtenido para los suelos agrícolas resultaría adecuado para lograr una disponibilidad óptima de macro y micro elementos, sin sufrir deficiencias de nutrimentos esenciales para los cultivos (^{Núñez, 1985}).

Para el caso particular de las muestras de suelo localizadas en los municipios de Cuerámaro y Abasolo (i.e., MACJI, MACNA, ABA5), las cuales presentaron pH relativamente más alcalinos (entre 8.14 y 8.69; Tabla 4), los resultados probablemente se expliquen por la geología y topografía de esta zona. Los municipios de Abasolo y Pénjamo actúan como cubetas que reciben las escorrentías y el lavado del material rocoso y edafológico de las zonas circundante más elevadas (Servicio Geológico Mexicano, 1999), confiriéndole a los suelos de ciertas partes de estos municipios una basicidad mayor.

En este estudio las muestras MASE, MASFE y MAC (Figura 1) presentaron las concentraciones más altas de MO (6.68% - 14.39%; tabla 4), lo que indicaría suelos con una gran cantidad de sustancias húmicas y una alta reactividad debido a sus superficies específicas elevadas (^{Porta et al., 2003}; ^{Sumner, 1999}). Estos suelos particularmente fueron manejados mediante labranza de conservación con cubiertas de paja, restos de rastrojos y regados con lixiviados de lombricomposta. Asimismo, con base en lo propuesto por ^{Oades (1984}), los componentes orgánicos de las muestras analizadas les otorgarían estabilidad a los agregados del suelo, mejorando la estructura de los mismos y contribuyendo a la disminución del riesgo de sellado y encostramiento superficial. Otra ventaja de las partículas orgánicas del suelo, es que constituyen una fuente de alimentos y energía metabólica para los microorganismos, la meso y macro fauna del suelo, liberando nutrientes al mineralizarse y neoformando sustancias húmicas al biodegradarse. Por otra parte, la muestra que presentó una categoría baja en MO fue la muestra IRP1 (0.88%; Tabla 4 y Figura 1), un suelo con manejo convencional y de temporal.

Los suelos extremadamente ricos en NT fueron MAJ, MAC, MASFE, MASE, MASIM, MACJI e IRP3 (NT > 0.221%; Tabla 4 y Figura 1). Estos suelos se caracterizan por contener restos de paja en descomposición de sus cultivos anteriores incorporados al suelo. La matriz de correlación mostró que el CO tiene correlación positiva con el NT (r _s = 0.67980, valor-p = 0.0019). La paja al contener una alta proporción de C/N favorece el crecimiento y multiplicación de organismos en el suelo. Cuando la paja se descompone y la mayor parte de los materiales ricos en energía han sido utilizados, la actividad microbiana disminuye y el nitrógeno nuevamente se acumula en el suelo (^{Miller, Graves, Williams & Madson, 1989}). Estos suelos se comportarían como reservas importantes de nitrógeno y como fuentes potenciales que favorecen la transformación de compuestos solubles o gaseosos de distintas formas nitrogenadas. La muestra MAA, (Figura 1) con una concentración de NT de 0.07%, correspondió únicamente a la categoría medianamente pobre. Este es un suelo agrícola con 10 años de producción con dos temporadas de siembra por año con uso de fertilizantes químicos.

Las muestras MASIM, IRA, MASE y MAC presentaron la relación C/N más aproximada a la recomendable para los suelos de uso agrícola (C/N = 10 - 12; ^{Kass, 1996}; Tabla 4 y Figura 1), logrando un equilibrio óptimo entre los procesos de descomposición, mineralización y humificación. La relación C/N del suelo MAA (C/N = 29.06; Tabla 4) resultó más alta que el valor límite óptimo recomendable por la norma mexicana de suelos. Este valor de C/N indica que el aporte de los nutrimentos al suelo como nitrógeno, fósforo y azufre, por la descomposición de la materia orgánica, es más lento. Además, la muestra agrícola MAA presentó un contenido bajo en NT total (0.07%), lo que indicaría una posible competencia del nitrógeno orgánico entre los microorganismos y los cultivos del suelo.

Análisis de elementos traza e I_geo

^{Rudnick & Gao (2003}) reportaron las concentraciones de ET de la corteza continental (B_n; Tabla 5). En este estudio se encontró que la concentración de As en los suelos analizados es superior al valor reportado para la corteza (Figura 2), principalmente en la muestra MASE (14.97 mg kg^-1), aunque no rebasa el valor de la norma mexicana de metales en suelos (Tabla 5). En cuanto al I_geo, este suelo perteneció a la clase 2 (moderadamente contaminado; Figura 3). Muchos factores ambientales, como el pH, potencial redox (Eh), la presencia de óxidos de Fe y Al, la textura, la capacidad de intercambio catiónico y la materia orgánica, controlan la concentración del As en los suelos (^{Chen, Ma & Harris, 1999}; ²⁰⁰²). Según ^{Varsányi & Kovács (2006)} y ^{García, Scracek, Fernández & Del Valle (2007)}, la movilidad del As se encuentra ligada fundamentalmente a los procesos de adsorción-desorción con los minerales de Fe, Mn y Al, y con la materia orgánica, ante condiciones de pH y Eh específicas.

Fuente: Elaboración propia.

Figura 2 Índice de geoacumulación (I_geo) de los suelos analizados. 

Fuente: Elaboración propia

Figura 3 Índice de geoacumulación (I_geo) de los suelos analizados. 

En el caso especial del suelo MAC, ubicado en la comunidad de “El Copal” (municipio de Irapuato; Figura 1), si bien ha sido regado en los últimos años con agua procedente de pozos profundos (más de 100 m) con concentraciones elevadas de As (0.33 mg L^-1; ^{Rodríguez, Morales-Arredondo & Rodríguez, 2016}) que superan la norma mexicana de aguas para consumo humano (0.025 mg L^-1; NOM-127-SSA1-1994 (^{DOF, 2000})), este metaloide no fue detectado en la muestra de suelo (Tabla 5). Es importante destacar que en el área de estudio (principalmente el noreste del municipio de Irapuato) existen rocas basálticas constituidas por silicatos de Al y Fe y óxidos de Fe, por lo que se esperaría la presencia de cantidades importantes de estos componentes inorgánicos en el suelo. Es probable que las fases portadoras fundamentales del metaloide en la muestra MAC lo constituyan los oxihidróxidos de Fe, que bajo condiciones oxidantes y de pH elevadas puedan ocasionar la desorción del elemento al medio. Esto sería avalado por la escasa o nula correlación encontrada entre el As y la materia orgánica como así también, la correlación negativa entre el As y el pH (Tabla 3).

Con base en la norma mexicana para suelos de uso agrícola, el elemento Cd se encontró en valores muy por debajo de la concentración de referencia total (CR_T) (Tabla 5). Sin embargo, comparando los valores reportados en la corteza continental (Figura 2) y las concentraciones cuantificadas de Cd para IRA e IRP4 (Figuras 1 y 3), estos resultaron con clase de I_geo 1 (no contaminado a moderadamente contaminado). A los suelos agrícolas MASE, IRP1, IRP2, IRP3, JUA2, CUE3 y ABA5 les correspondió una clase de I_geo 2 mientras que MAJ y MAC resultaron con clase de I_geo 3 (de moderado a fuertemente contaminado; tabla 5; figuras 1 y 3). Esto indicaría que su acumulación procede de fuentes antropogénicas. Como fue reportado por ^{Alloway (2012)}, la aplicación de biósolidos, estiércol y las partículas sólidas presentes en el humo de los incendios y las emisiones de diversas industrias son una fuente de aporte de este contaminante en los suelos.

El promedio del contenido de Cr en los suelos analizados fue de 21.50 mg kg^-1 (Tabla 5). En la corteza continental, el Cr alcanza los 92 mg kg^-1 por lo que ninguna concentración de Cr cuantificada en este estudio supera la de la corteza (Figura 2), lo cual sugiere que los suelos agrícolas no estarían contaminados por este metal.

Por su parte, el Cu se encuentra normalmente en el suelo en una concentración de 28 mg kg^-1 - 30 mg kg^-1, acumulado normalmente en los horizontes superficiales. El Cu puede llegar a la superficie del suelo por la aplicación de fertilizantes foliares para el crecimiento y salud de las plantas, como un complemento a las deficiencias de este micronutriente esencial en el suelo (^{Machender et al., 2011}; ^{Wuana & Okieimen, 2011}). El valor más elevado de Cu se encontró en el suelo IRA, localizado en el área NE del municipio de Irapuato (Figura 1) con un valor de 99.98 mg kg^-1 (Tabla 5). Este suelo es trabajado mediante agricultura convencional lo que implica una aplicación alta de agroquímicos que acumulan metales pesados en el suelo. La clase de I_geo para esta muestra de suelo resultó de 2 (moderadamente contaminado).

El contenido de Pb en el suelo se asocia principalmente con minerales de arcilla, óxidos de Mn, hidróxidos (Fe y Al) y con la materia orgánica. El contenido de Pb de los suelos superficiales a escala mundial se ha estimado en 25 mg kg^-1; valores por encima de este indican una influencia antropogénica derivada de la emisión aérea de la combustión de gasolina con tetraetilo de plomo (^{Wuana & Okieimen, 2011}). Los suelos agrícolas MASE, MAJ, IRA Y JUA2 mostraron una clase de I_geo 1 (no contaminado a moderadamente contaminado, Figuras 1 y 3), todos ubicados en las adyacencias a zonas urbanas y/o a las carreteras. El Pb mostró correlación negativa estadística significativa con el pH (r_s = -0.48244, valor-p = 0.0426). El valor del pH desempeña un importante papel en la disponibilidad del Pb, cuanto menor sea el pH mayor cantidad de Pb se encontrará disponible en la solución del suelo. A pH alto es un elemento que precipita como hidróxido, fosfato y carbonato, o puede ligarse también a la materia orgánica. El Pb también mostró una correlación positiva con el As (r_s= 0.51756, valor-p = 0.0278), con el Cd (r_s= 0.50828 valor-p = 0.0313) y con el Cu (r_s= 0.6099, valor-p = 0.0072; Tabla 3).

El Zn pertenece al grupo de elementos trazas esenciales para la vida y es muy peligroso en concentraciones altas, ya que influye negativamente en la actividad de los microorganismos y lombrices de tierra, retardando la descomposición de la materia orgánica (^{Wuana & Okieimen, 2011}). Las fuentes principales de contaminación de Zn son las industrias, los lixiviados de estiércol y compostas, fertilizantes y pesticidas (^{Gowd, Reddy & Govil, 2010}). De esta manera, el suelo MASE regado con lixiviados de lombricomposta es que el presentó el contenido más alto de Zn de todos los suelos analizados (109.44 mg kg^-1; Tabla 5). En las rocas de la corteza continental se reportó una concentración de 67 mg kg^-1 de Zn, resultando las muestras MASE e IRA en la categoría 1 de I_geo (no contaminado a moderadamente contaminado; Figuras 1 y 3), ambos suelos sometidos a fertilización por agroquímicos.

En general, las correlaciones positivas del Pb, Cd, Cu y Zn concuerdan con lo reportado en el trabajo de ^{Alloway (2012)}, donde se indica que estos elementos se encuentran asociados y su contaminación en los suelos proviene de los fertilizantes aplicados para dar crecimiento a los cultivos, así como de zonas urbanas, emisiones del tránsito vehicular, pinturas y muchas otras fuentes no específicas.