¹Profesor Auxiliar Universidad del Magdalena. Colombia. jose.quevas@yahoo.es ²Profesores Asociados. Universidad Nacional de Colombia sede Palmira, A.A. 237. Palmira, Valle del Cauca, Colombia. Autor para correspondencia: dbaenag@palmira.unal.edu.co, jcmenjivarf@palmira.unal.edu.co

En suelos destinados a la actividad agrícola de la granja experimental de la Universidad del Magdalena, Santa Marta - Colombia (74° 07' Y 74° 12' longitud oeste y 11° 11' y 11° 15' latitud), se estudió la variabilidad espacial de las propiedades físicas: humedad gravimétrica, textura, densidad aparente (Da), resistencia a la penetración (RP) y conductividad hidráulica (KH) y químicas: Materia orgánica (MO), pH, CIC, N, P, K, Ca⁺², Mg⁺², Na⁺, Cu, Zn, Fe, Mn, B), utilizando una grid o malla de 29m*29m con 184 puntos de muestreo. El análisis de los datos se realizó empleando enfoques univariado, geoestadístico y multivariado. El análisis estadístico de los resultados mostró gran variabilidad de las propiedades; la densidad aparente y el pH presentaron baja variación; los componentes texturales, la RP, Ca⁺², Mg⁺², Cu y B mostraron variación media, mientras que la humedad gravimétrica, KH, MO, K⁺, Na⁺, CIC, Zn, P, Mn y Fe, expresaron alta variación. Con el análisis Geoestadístico se visualizó y graficó la variabilidad de cada propiedad evaluada y se evidenció la existencia de moderada dependencia espacial en las propiedades excepto en los contenidos de arena, Ca⁺², CIC, P y Cu que revelan fuerte dependencia espacial. El índice de homogeneidad multivariado para los suelos en estudio fue bajo, siendo textura, Na⁺, B, Fe, K, pH, Ca y MO las propiedades con mayor contribución a la heterogeneidad. El análisis de clasificación permitió agrupar cuatro zonas con condiciones edáficas similares, las que se podrán utilizar para determinar el sistema de muestreo, número de muestras, tamaño y distribución de unidades experimentales y con ello facilitar la implementación de prácticas de manejo agronómico acordes con las características de cada zona.

Palabra clave: Variabilidad espacial, propiedades físicas y químicas, geoestadística, análisis multivariado, zonas homogéneas.

We have studied the spatial variability of agricultural soils in the experimental farm at Magdalena University (Santa Marta, Colombia; 74° 07' – 74° 12' W, 11° 11' – 11° 15' N). Analyzed soil properties included: gravimetric soil moisture content, grain size analysis, bulk density, resistance to root penetration, hydraulic conductivity, organic matter content, pH, cation exchange capacity, total N and P, exchangeable cations (Ca⁺², Mg⁺², K⁺, Na⁺), metals (Cu, Zn, Fe , Mn) and B content. Sampling was based on a 29 x 29 m grid, resulting in 184 sampling points. Data analysis was done using univariate, multivariate and geostatistical approaches. The analyzed soil properties showed different degrees of variability: bulk density and pH low variability; grain size, resistance to root penetration, divalent exchangeable cations (Ca⁺² and Mg⁺²), Cu and B intermediate variability; while moisture content, hydraulic conductivity, organic matter content, monovalent exchangeable cations (K⁺, Na⁺), cation exchange capacity, Zn, P Mn and Fe a large variability. The geostatistical analysis revealed the presence of moderate spatial dependency, except for sand content, exchangeable Ca⁺², cation exchange capacity, P and Cu, which showed strong spatial dependency. The multivariate homogeneity index was low, being exchangeable Na⁺, B, Fe, K, pH and organic matter content the soil properties with lager contributions to heterogeneity. Classification analysis enabled to indentify four areas with similar edaphic conditions, which will allow to stablish the sampling strategy, number of samples, size and distribution of experimental units and thus to improve the implementation of agronomic management practices in agreement with the characteristics of the soils.

Key words: Spatial variability, physical and chemical properties, geostatistics, multivariate analysis and homogeneous zones.

El objetivo del presente trabajo fue caracterizar la variabilidad espacial de las propiedades físicas y químicas de los suelos de la Granja Experimental de la Universidad del Magdalena destinados a la docencia, la investigación y la extensión.

El trabajo se desarrolló en la Granja Experimental de la Universidad del Magdalena, Santa Marta D.T.C.H, Colombia (74° 07' y 74° 12' de longitud oeste y 11° 11' y 11° 15' de latitud norte), a 15 m.s.n.m, clima estepario tropical cálido y formación vegetal de bosque tropical seco, vegetación xerofítica, precipitación media anual de 674mm, de carácter unimodal distribuida entre abril y noviembre, temperatura promedio de 28.5°C, humedad relativa 75% y vientos alisios del hemisferio norte (IDEAM, 2001). Los suelos de la granja son de origen aluvial, poco evolucionados con predominio de texturas gruesas, clasificados en el orden Entisoles y suborden Psamments (Lobato, 2003). Se realizó un muestreo preliminar de densidad aparente, pH y materia orgánica (propiedad de mayor variaciación) para determinar el tamaño de muestra, estableciendo dos transectos de 40m de longitud en direcciones ortogonales e intervalos de 5m. Con los valores obtenidos, se estimó el tamaño de muestra de 184 sitios, los cuales fueron distribuidos en el área de forma sistemática, utilizando un intervalo de muestreo de 29m*29m, sobre el mapa base del Campo Experimental, realizando observaciones en las intersecciones de la cuadrícula y georreferenciando cada punto con GPS 12 XL GARMIN (Figura 1) En los puntos de muestreo se midieron las propiedades físicas (humedad, textura, densidad aparente, conductividad hidráulica y resistencia a la penetración) y químicas (MO, pH, CIC, N, P, ,K⁺, Ca,⁺², Mg⁺², Na⁺, Cu, Zn, Fe, Mn, B) a 25 cm. de profundidad. Los análisis se realizaron según las metodologías descritas en el Manual de laboratorio de suelos y tejidos vegetales del Instituto Geográfico Agustín Codazzi (IGAC) 2006.

El análisis de la información se realizó con enfoque univariado, geoestadístico y multivariado. Para el análisis geoestadístico, se corrieron los datos en el programa GS+ (Geostatistics for the Enviromental Sciences Versión 5.1 de Gamma Design Software), se obtuvo el semivariograma experimental para cada propiedad y se estableció el modelo teórico de mejor ajuste, teniendo en cuenta como parámetros de decisión el coeficiente de determinación (R²) y la suma de cuadrados de los residuales (RSS). Posteriormente se realizó el análisis de validación cruzada, mediante el método de interpolación Kriging puntual, atendiendo la metodología sugerida por Moral 2004, que fue base para la construcción de mapas temáticos que permitió establecer la variabilidad espacial de las propiedades evaluadas. Para el análisis multivariado, se realizó análisis factorial de componentes principales a partir de los cuales se seleccionaron aquellos con valores propios mayores o iguales a uno, explicando así el mayor porcentaje de la variación total de los datos y análisis de clasificación o Clúster Análisis (método jerárquico) con utilización del algoritmo de Ward (mínima inercia), a fin de conformar grupos de suelos con características semejantes y delimitar zonas homogéneas.

Enfoque univariado. De acuerdo con el grado de variabilidad de cada propiedad, definido en términos del coeficiente de variación (CV%), se establecieron cuatro categorías de variables: relativamente homogéneas con CV%<20% (Da y pH), moderadamente heterogéneas con 20<CV%<40% (textura y Ca⁺²), normalmente heterogéneas 40<CV%<60% (RP, KH, MO, Mg⁺², Cu y B) y extremadamente heterogéneas CV%>60% (humedad, K⁺, Na⁺, CIC, P, Zn, Mn e Fe). De lo anterior se deduce que con excepción del Ca⁺² y el pH, las propiedades químicas de los suelos de la granja presentan niveles de variación altos, asociados con una evidente heterogeneidad en toda el área de estudio. Los valores de densidad aparente determinados por el método del terrón parafinado (IGAC, 2006) en promedio pueden considerarse de normales a ligeramente altos para suelos de textura media; sin embargo, la propiedad no tiene variación significativa entre sitios de muestreo, posiblemente debido al efecto del uso y manejo y a las condiciones ambientales (baja precipitación y altas temperaturas) de la época seca en que se realizó el muestreo. El 69.6% de los sitios presentan predominio de textura franco arcillo arenosa, evidenciando la intergradación de texturas gruesas a medias, propiedad de los suelos de origen aluvial (Lobato, 2003). El promedio de humedad (3%), evidencia las condiciones climáticas imperantes en la época de muestreo. La resistencia a la penetración con promedio de 2.1 MPa se presenta en niveles relativamente limitantes para los cultivos (Salamanca, 2000), mientras que la conductividad hidráulica con promedio de 1.58cm/h se califica como lenta a moderada. De acuerdo con las consideraciones de Ortega y Flores, 1999 y del Manual de Laboratorio del ICA 1992, el pH osciló entre ligeramente alcalino y alcalino, con máximo de 9.7 muy alcalino, situación extremadamente limitante para el desarrollo adecuado de los cultivos y mínimo de 6.8. De acuerdo con los valores promedio, las propiedades químicas se pueden estratificar de la siguiente manera: niveles de bajo a medio: MO, Fe, medio: Cu, Zn, CIC y alto: bases, P, Mn y B (Cuadro 1).

Análisis de la variabilidad enfoque geoestadístico. Teniendo en cuenta lo enunciado por Cambardella et al., 1994, citados por Oliveira et al., 1999 y Mejía (2007), el contenido de arena presentó una fuerte dependencia espacial lo cual es válido si se tiene en cuenta que los cambios texturales entre muestras vecinas son mínimos. Con excepción de las propiedades Ca, CIC, P y Cu con débil dependencia espacial, las propiedades restantes presentaron porcentajes de nugget relativo entre 25%-75% que corresponde a una moderada dependencia espacial. El alcance o rango efectivo, dentro del cual se mantiene la dependencia espacial varió de una propiedad a otra. La densidad aparente, la resistencia a la penetración, limo, conductividad hidráulica, pH, MO, CIC, Zn, Mn y el Fe presentan un alcance mayor a la máxima longitud del lote y por tanto, la dependencia espacial para estas propiedades se manifiesta a lo largo y ancho del mismo. Los semivariogramas para las propiedades físicas y químicas presentan estructuras diferentes de dependencia espacial, puesto que en ningún caso la semivarianza es constante en función de la distancia. Se observan además diferencias en los parámetros geoestadísticos y comportamiento característico de poblaciones entremezcladas (Reichart, 1985), debido a la presencia de condiciones de suelo diferentes a nivel de detalle en el área de estudio. Los modelos teóricos de mejor ajuste en la mayoría de los casos fueron el esférico y el exponencial (Cuadro 2), lo cual permite afirmar que las variables estudiadas cumplen con la hipótesis de fuerte estacionalidad, dado que sus modelos son acotados (Samper y Carrera, 1990).

Los niveles de pH predominantes en las cuatro condiciones de suelo limitan la absorción de nutrientes por los cultivos, en este sentido los valores de los macro y micro nutrientes (considerados altos en la descripción de cada clúster), se deben analizar con detenimiento y mucho cuidado al realizar cualquier recomendación ya que solo representan los contenidos de la fase intercambiable del suelo y podrían enmascarar posibles deficiencias o toxicidades que se manifiesten en los cultivos que se proyecten, por esta razón, para evaluar la fertilidad física y química de las cuatro condiciones de suelo agrupadas es necesario complementar la información determinando los cationes y aniones en solución, carbonatos, bicarbonatos, porosidad, infiltración, contenidos de humedad a capacidad de campo y punto de marchitez permanente.

La variabilidad espacial de los atributos del suelo en el lote experimental es alta, siendo mayor en las propiedades químicas que en las físicas; con excepción del Ca⁺² y el pH, las propiedades químicas de los suelos de la granja presentan niveles de variación altos, asociados con una evidente heterogeneidad en toda el área de estudio. El separado arena presentó una fuerte dependencia espacial lo cual es válido si se tiene en cuenta que los cambios texturales entre muestras vecinas son mínimos. Con excepción de las propiedades Ca⁺², CIC, P y Cu con débil dependencia espacial, las propiedades restantes mostraron una moderada dependencia espacial.

La densidad aparente, la resistencia a la penetración, limo, conductividad hidráulica, pH, MO, CIC, Zn, Mn y el Fe presentan un alcance mayor a la máxima longitud del lote (entre 594 m y 1010 m) y por tanto la dependencia espacial para estas propiedades se manifiesta a lo largo y ancho del lote.

Con el análisis de componentes principales de las variables físicas y químicas evaluadas en el presente estudio, se conformaron tres y cuatro componentes principales respectivamente que explicaron el 72% y 66% de la variación, siendo las variables textura, Na⁺ y pH las que más influyen en la diferenciación de las condiciones de suelos.

Los enfoques univariado, multivariado y geoestadístico para el análisis de información, son herramientas que utilizadas individualmente o en conjunto ayudan en el análisis de la variabilidad espacial de los atributos del suelo y así, sus potencialidades o limitantes. Además, el mapa de zonas homogéneas debe constituirse entonces en instrumento de ayuda que guíe la toma de decisiones a fin de implementar prácticas acordes con las características físicas y químicas específicas (Figura 2).

Burrough, P. A. 1986. Soil variability: a late 20th century view. Soils and Fertilizers. 56 (5): 529-562.        [ Links ]

Cambardella, C. A.; Moorman T. B.; Konopka, A. E. 1994. Field scale variability in central Iowa soils. Soil Sci. 58:1501-1511.        [ Links ]

Giraldo, S.; Mendoza, A. 2003 Caracterización de la variabilidad espacial de las propiedades físicas de los lotes experimentales destinados a los cultivos de hortalizas en el CEUNP. 57p.        [ Links ]

Instituto colombiano agropecuario. 1992. Manual del Laboratorio de Suelos y Aguas. 350p.        [ Links ]

Instituto colombiano de hidrología y meteorología, IDEAM. 2001. Sección de meteorología sinóptica aeronáutica. Datos diarios. Estación Simón Bolívar Santa Marta, Magdalena: IDEAM (mimeografiado).        [ Links ]

Instituto geografico agustin codazzi. 2006. Métodos analíticos del Laboratorio de Suelos. 400p.        [ Links ]

Jaramillo, D. 2003. Introducción a la Ciencia del Suelo. Universidad Nacional de Colombia. Sede Medellín. p 553-572        [ Links ]

Lobato, R. 2003. Mapa taxonómico de los suelos de la Universidad del Magdalena. 36p.        [ Links ]

Mejia, A. 2007. Variabilidad espacial de las Propiedades Físicas y Químicas en un área productiva de palma africana (Elaeis guineensis L). Tesis de grado para optar al título de Magíster en Ciencias Agrarias. Universidad Nacional de Colombia. Pasto. 171p        [ Links ]

Moral, F. J. 2004. Aplicación de la geoestadística en las ciencias ambientales. Ecosistemas. Vol 13, No. 1, pp. 78-86. Disponible en Internet: http://www.aeet.org/ecosistemas/041/documentos/Revision3.doc [Consultada: 10 Enero. 2006.]        [ Links ]

Nielsen, D. R.; Tillotson, P. M.; Vieira, S. R. 1983. Analysing field measured soil water properties. Agr.Water.Manag. Amsterdam. v 6, p.93-109.        [ Links ]

Obando Moncayo, A. Et al. 2006. Variabilidad Espacial de Propiedades Químicas Y Físicas En Un Typic Udivitrands, Arenoso de la Región Andina Central Colombiana. En revista Facultad de Agronomía. Universidad Nacional sede Medellín. Vol 59, No. 1 p. 3217-3235.        [ Links ]

Oliveira, J. J. Et al. 1999. Variabilidad espacial de propiedades químicas en un solo salino-sódico. Revista Brasileira de Ciencia do Solo. 23: 783-789.        [ Links ]

Ortega, R. Flores, L. 1999. Introducción al manejo sitio - específico. Ministerio de Agricultura. Instituto de Investigaciones Agropecuarias. (ed.) En Agricultura de precisión. CRI Quilamapu-Chile. p13-46.        [ Links ]

Ovalles, F. 1992. Metodología para determinar la superficie representada por muestras tomadas con fines de fertilidad. FONAIAP-CENIAP-IIAG. Maracay. 44 p.        [ Links ]

Ponce, D., et al. 1999. Variabilidad espacial del pH, fósforo y potasio en muestreos de suelos con fines de fertilidad en plantaciones de caña de azúcar en Cuba. En Bioagro. Vol 11, No.1 p. 3-11.        [ Links ]

Reichart, C. 1985. Variabilidad espacial (e Temporal) de solos. Procesos de transferencia no Sistema Solo-Planta-Atmosfera. Cuarta edición. Fundacao Cargill. p. 391-417.        [ Links ]

Salamanca, A. 2000. Influencia de las características físicas del suelo y su interacción con las condiciones climáticas en el comportamiento y calidad del Desmodium heteracarpon subsp. Ovalifolium (OHASHI) en tres regiones de Colombia. Tesis de grado para optar el título de ingeniero agrónomo. Universidad Nacional de Colombia. Palmira. 170p        [ Links ]

Samper, J., Carrera, L. 1990. Geoestadistica. CIMNE. Barcelona. España. 484p.        [ Links ] ]]> 1986 56 5 5 529-562 1994 58 1501-1511 2003 57p. Instituto colombiano agropecuario 1992 350p. Instituto colombiano de hidrología y meteorología 2001 Instituto geografico agustin codazzi 2006 400p. 2003 p 553-572 2003 36p. 2007 171p 2004 Vol 13 No. 1 No. 1 pp. 78-86 1983 v 6 p.93-109 2006 Vol 59 No. 1 No. 1 p. 3217-3235 1999 23 783-789 1999 13-46 1992 44 p. 1999 Vol 11 No.1 No.1 p. 3-11 1985 Cuarta edición p. 391-417. 2000 170p 1990 484p.