El análisis multivariado de componentes principales y de agrupamiento jerárquico se hizo con el fin de identificar relaciones entre las variables estudiadas en los sistemas evaluados. En el análisis multivariado, se realizó análisis de factores por componentes principales (CP) (Manly, 1997). Para el análisis de agrupamiento se utilizaron las distancias de Ward, lo cual permitió identificar los grupos con base en el comportamiento de las variables evaluadas.

La tabla 3 presenta una descripción estadística de las variables químicas seleccionadas en el análisis multivariado relacionadas con los sistemas productivos de la Altillanura y Piedemonte.

Tabla 3


Fuente: Esta investigación.

Granada: S2= Monocultivo de soya en rotación con maíz; S3= Cultivo de maíz en rotación con yuca; S4= Monocultivo de plátano en rotación con yuca; S5= Monocultivo de arroz en rotación con maíz y soya; Villavicencio: S7= Sistema agroforestal de cacao y Acacia mangium; S9= Sistema agroforestal de café con plátano y árboles de leguminosas; S10= Monocultivo de arroz; Puerto López: S11= Bosque secundario; S12= Sistema silvopastoril de Acacia mangium y pastura mejorada; S13= Sistema agroforestal de caucho y leguminosas de cobertura; S14= Monocultivo de piña; S15= Monocultivo de cacao.

El pH del suelo varió entre 3.90 y 5.10, valores que correspondieron a los SAFs de cacao y Acacia mangium de Villavicencio (S7) y al SAF de café asociado de Villavicencio (S9), respectivamente. Según Rivera y Amézquita (2013), el pH presentó valores entre 4.0 y 4.4, siendo más alto en sabana nativa hasta 20 cm de profundidad.

Se comprobó que los diferentes manejos y usos del suelo presentaron influencia en la MOS, variando desde 1.60 a 4.20% en el sistema de monocultivo de piña de Puerto López (S14) y en el silvopastoril de Acacia mangium y pastura mejorada de Puerto López (S12).

El promedio de P disponible en el suelo varió entre 1.70 y 67.50 mg/kg de suelo; el menor valor correspondió al sistema agroforestal de caucho y leguminosas de cobertura (S13) y monocultivo de soya en rotación con maíz (S2), siendo más altos que los reportados por Rao et al. (1999).

Al respecto de los valores mínimos de P encontrados, "los suelos tropicales altamente meteorizados (Oxisoles, Ultisoles) se caracterizan por niveles bajos de P total y disponible y por la alta capacidad de fijación de P” (Friesen et al., 2013, p. 148).

En comparación con los otros nutrientes principales, el P es el menos móvil y disponible para las plantas en las condiciones de la mayoría de los suelos, particularmente en Oxisoles y, por lo tanto, es probable que se vea muy afectado por la labranza (Basamba et al., 2006; 2007). “Los valores de P reflejan el efecto de la intervención agrícola, con contenidos mayores a 5 mg/kg” (Camacho, 2010, p. 279). Sin embargo, estos mismos autores reportaron que Jaimes et al. (2003) encontraron contenidos inferiores a 1 mg/kg en sabana nativa de los Llanos Orientales.

Las leguminosas favorecen la disponibilidad del P en los suelos (Jaramillo, 2002), el crecimiento sustancialmente mayor de las raíces de las pasturas mejoradas (Brachiaria humidicola asociada con leguminosas) en las capas inferiores de suelo en condiciones ácidas con bajos niveles de P, refleja la mayor capacidad de los forrajes para movilizar y tomar P del suelo (Rao et al., 1997).

Según Friesen et al. (2013), durante los últimos 20 años se han identificado numerosos germoplasmas de pastos y forrajes, adaptados a bajos niveles de P en Oxisoles de los trópicos, pero las características de las plantas y los mecanismos que contribuyen a una mayor adquisición de P y/o utilización en estos genotipos son poco conocidos.

El Al intercambiable varió de 0.05 a 5.70 cmol (+)/kg de suelo; el menor valor correspondió al cultivo de soya en rotación con maíz de Granada (S2) y al SAF de cacao y Acacia mangium de Villavicencio (S7), siendo el máximo más alto que el reportado por Rivera y Amézquita (2013), donde el Al intercambiable varió entre 1.7 y 2.9 cmol (+)/kg de suelo en sistemas de maíz en monocultivo y en rotación con soya entre 0 y 10 cm de profundidad y en sabana nativa entre 0 y 20 cm de profundidad, respectivamente. El valor más alto encontrado estuvo posiblemente asociado a un pH más acido en los sistemas con mayores aportes de MOS. El Al intercambiable en suelos de carga variable, dependiente del pH como en Oxisoles, se debe neutralizar con altas aplicaciones de cal (Flores et al, 2008).

Los contenidos de Ca variaron de 0.05-5.70 cmol (+)/kg suelo considerados bajos y medios; siendo el más bajo en SAF de cacao y Acacia mangium de Puerto López (S12) y el valor máximo correspondió a SAFs de café asociado de Villavicencio (S9); estos valores estuvieron asociados con el pH ácido del suelo en estos dos sistemas; el Mg varió de 0.01 a 0.80 cmol (+)/kg, considerados bajos en los sistemas monocultivo de arroz (S10) de Villavicencio y monocultivo de soya en rotación (S2) de Granada; K vario entre 0.07 a 0.38 cmol (+)/kg de suelo, contenidos que son considerados bajos y medios en el monocultivo de arroz en rotación con soya y maíz de Granada (S5) y en monocultivo de plátano en rotación con yuca de Granada (S4), respectivamente.

Al respecto, Kirkby (1979) menciona que los Oxisoles y Ultisoles predominantes en los Llanos Orientales de Colombia se caracterizan por un alto contenido de arcillas caolinitas tipo 1:1, que requieren menos nivel de saturación de bases para una adecuada disponibilidad de Ca y Mg para las plantas.

Kamprath (1971) menciona que la disponibilidad de Ca y Mg en Oxisoles y Ultisoles es afectada por prácticas de manejo como labranza, rotación de cultivos y aplicación de enmiendas o fertilizantes, al igual que por los niveles de la forma intercambiable, la capacidad de intercambio catiónico efectiva, el Al intercambiable, la textura y la mineralogía de las arcillas.

La aplicación de Ca como enmienda, ayuda a elevar el nivel de este y del Mg en el suelo (Fassbender, 1972), favoreciendo el movimiento de ambos hacia otras partes del suelo y aumenta el pH (Morelli et al., 1971), sin afectar mucho al K.

El Na no es una limitante en estos suelos (Jaramillo, 2002); en cuanto a los elementos menores, Fe y Mn fueron elevados en cultivo de soya en rotación con maíz (S2); y en cultivo de maíz en rotación con yuca de Granada (S3); B fue de bajo a alto en sistema agroforestal de caucho y leguminosas de cobertura de Puerto López (S13) y en monocultivo de soya en rotación de Granada (S2). La fijación biológica del N ayuda a aumentar los contenidos de B en los suelos (Moreira y Siqueira, 2006).

El análisis de correlaciones entre las variables asociadas a la fertilidad química de los suelos de los sistemas productivos de la Altillanura y de Piedemonte se observa en la tabla 4.

Tabla 4


Fuente: Esta investigación.

Los resultados obtenidos confirman que la materia orgánica del suelo (MOS), correlacionó con todas las características químicas del suelo, menos con B (Tabla 4), como reportado por West et al (2010). El P estuvo altamente correlacionado con la MOS (Tabla 4).

La correlación negativa y positiva encontrada entre el pH con el Al intercambiable y con Ca, Mg, K y Na (Tabla 4), estaría indicando que a un menor pH es mayor el Al y menor el contenido de bases (Fassbender, 1972). El Ca correlacionó con el K, el Mg con el Zn y el K con el Mn (Tabla 4), existiendo sinergia entre nutrientes y su influencia en el mejoramiento de la fertilidad de los suelos (Moreira y Siqueira, 2006). El Na no es una limitante en estos suelos, como si lo es el Al (Jaramillo, 2002).

Stemberg (1998)) al evaluar aspectos físicos, químicos y biológicos de los suelos destacó la importancia de las técnicas de análisis multivariado (componentes principales), las cuales facilitaron la clasificación de la calidad de los suelos en tres categorías (buena, normal y baja). En el análisis de componentes principales se seleccionaron los primeros tres componentes, que explicaron el 79.63% de la variabilidad de los datos.

Tabla 5


Fuente: Esta investigación.

La conformación del CP1 estuvo relacionado principalmente con la MOS, pH, Al, todas las bases y micronutrientes del suelo, menos el Fe; incluido el Boro. El P no guardo relación. El primer componente explicó el 50.54% de la variabilidad total (Tabla 5).

El CP2 explicó el 17.25% de la variabilidad total y se relacionó principalmente con la acidez del suelo, pH, Al, Na, Fe y B en los suelos (Tabla 5). Las variables de mayor peso en la conformación del CP3 fueron Cu, Zn y P en los sistemas evaluados, explicando el 11.84% de la variabilidad total de las propiedades químicas (Tabla 5).

El análisis de Cluster estuvo conformado de la siguiente manera (Figura 1).

Figura 1
Análisis de clúster asociados a la fertilidad química de los suelos en sistemas productivos de la Altillanura y de Piedemonte, Meta, Colombia, 2015.
Fuente: Esta investigación.

Cluster 1: Sistemas productivos con suelos medianamente fértiles. En este cluster se agruparon los sistemas de monocultivo de plátano en rotación con yuca de Granada (S4), Monocultivo de arroz de Villavicencio (S10), sistema agroforestal de caucho y leguminosas de cobertura (S13), monocultivo de piña (S14) y monocultivo de cacao de Puerto López (S15); parece ser que la alta heterogeneidad de los suelos y la alta variabilidad de las propiedades químicas no determinó un patrón específico en cuanto al tipo de uso del suelo; hubo predominio de suelos con bajo pH y altos contenidos de Al, altos contenidos de Fe y bajos contenidos de B en los suelos (Figura 1).

Eso muestra que los atributos mencionados poco contribuyeron para separar los sistemas productivos analizados dependientes de tipos de manejo y uso del suelo con características similares.

A pesar de que últimamente los sistemas de monocultivos en las tres zonas están incluyendo algún tipo de manejo sostenible de suelo como es la práctica de rotación de cultivos y la implementación de SAFs, la baja fertilidad química de estos suelos supera a las innovaciones tecnológicas de los sistemas. Según Vera (2004), White (1984) y Bernoux et al. (2006), la rotación de cultivos incrementa los contenidos de la MOS con resultados visibles a mediano y largo plazo, siendo que con los SAFs los resultados pueden darse en el corto plazo (Kirby y Potvin, 2007; Dossa et al., 2008; Somarriba et al., 2013; Jadan et al., 2015).

Cluster 2: Sistemas productivos con suelos de alta fertilidad. Agrupa los sistemas de pastura mejorada de Granada (S1); pastura mejorada de Villavicencio (S6); SAF de cacao y Acacia mangium de Villavicencio (S7); bosque secundario (S11) y sistema silvopastoril de Acacia mangium y pastura mejorada (S12) de Puerto López; la conformación de este clúster estuvo más asociado a la baja acidez de los suelos, menores contenidos de Al intercambiable, altos contenidos de bases y micronutrientes, pero principalmente a altos contenidos de MOS (Figura 1).

Doran y Parkin (1994) y Reeves (1997), han incluido a la materia orgánica del suelo (MOS) como la propiedad química más importante a seleccionar entre un grupo de datos mínimos y necesarios para definir la calidad del suelo en sistemas productivos.

En las pasturas mejoradas y los SAFs la materia orgánica del suelo (MOS) se acumula en los horizontes superiores (Bessam y Mrabet, 2003; Guimarães et al., 2004; Basamba et al, 2006), pudiendo explicar los resultados obtenidos, lo cual puede estar relacionado con los tipos de labranza reducida y/o cero y la alta entrada de residuos.

La recuperación de pasturas degradadas con la adopción de prácticas conservacionistas como las pasturas mejoradas y los SAFs, junto con el uso de correctivos y fertilizantes, evita la degradación del suelo y favorece incrementos en los contenidos de MOS y P principalmente (Nair, 1993; Carvalho et al., 2007; Giraldo et al., 2008; Latriglia y Vera, 2014).

Cluster 3: Sistemas productivos con suelos de baja fertilidad. Agrupa los sistemas de monocultivo de soya en rotación con maíz (S2), cultivo de maíz en rotación con yuca (S3) y monocultivo de arroz en rotación con maíz y soya (S5), todos éstos de Granada; monocultivo de piña (S8) y sistema agroforestal de café asociado de Villavicencio (S9); los suelos de estas dos zonas estuvieron asociados a alta acidez, que definió altos contenidos de Fe y Zn, y bajos contenidos de P-disponible (Figura 1).

Según Friesen et al. (2013), estudiando la adquisición y reciclaje de fósforo en sistemas de cultivos y pasturas en suelos tropicales de baja fertilidad,

Las diferencias en las tasas de liberación puede ser de menor importancia en las pasturas donde las raíces activas están siempre presentes para recuperar el P liberado de la hojarasca en descomposición. En los sistemas que involucran cultivos, sin embargo, el sistema radicular del cultivo no puede ser suficientemente desarrollado para recuperar el P liberado de los residuos del cultivo anterior. (p. 156)

De esta manera se pueden explicar las diferencias en el reciclaje de nutrientes de los sistemas productivos incluidos en los cluster 2 y 3.

Para el arroz, el patrón de liberación de P de los residuos aplicados a la superficie sugiere que <60% estará disponible para el siguiente cultivo de caupí (o componentes alternativos) en un sistema de rotación. Para el monocultivo de arroz, la liberación completa dentro de los 4 meses implica una interacción sin restricciones con el suelo de alrededor de 6 meses antes que el próximo cultivo presente un sistema radicular suficientemente desarrollado para adquirirlo. (Friesen et al., 2013, p. 156)

Por otra parte, en sistemas de monocultivo donde se realiza labranza intensiva es un factor que contribuye a la degradación de la fertilidad química de los suelos en el corto tiempo (Lal, 2008); debido principalmente a la liberación del carbono del suelo que sucede por una oxidación intensa de la MOS y las pérdidas de nutrientes por erosión (Chen et al., 2009; West et al., 2010); presentándose una disminución de los contenidos de P en los suelos de los sistemas productivos del clúster 3 y de micronutrientes en sistemas productivos del clúster 1.