Introducción
Los páramos son ecosistemas neotropicales entre el límite forestal superior y las nieves perpetuas; se encuentran distribuidos a lo largo de los Andes húmedos entre Perú, Ecuador, Colombia y Venezuela, con extensiones hasta Costa Rica y Panamá (García et al., 2019). Estos ecosistemas son fundamentales para la vida, dada su gran riqueza biológica y los numerosos servicios ecosistémicos que prestan a la sociedad (captación, almacenamiento y regulación hídrica, captura de carbono, regulación climática, formación y protección del suelo, entre otros) (Giné y Sánchez, 2016; Vargas, 2013). Sin embargo, hoy en día el cambio climático, el acelerado crecimiento poblacional y las diversas actividades económicas (agricultura, ganadería extensiva y minería) representan riesgos severos para la integridad de los páramos alrededor del mundo, ya que ejercen una enorme presión sobre los recursos naturales, lo que conduce a una pérdida de la calidad edáfica (Farfán et al., 2020; Fernández et al., 2019).
En Colombia, alrededor del 60 % de los páramos han sufrido procesos de transformación paisajística y ecosistémica como consecuencia de las actividades agropecuarias basadas en prácticas convencionales, entre ellas, los monocultivos de papa (Solanum tuberosum L.). El uso intensivo del suelo para la producción de papa (bajo el modelo de la revolución verde) ha ocasionado problemáticas ambientales como la erosión, la contaminación, la compactación, la pérdida de carbono orgánico y de biodiversidad, lo cual trae como consecuencia la degradación de los suelos (Avellaneda-Torres et al., 2018; Farfán et al., 2020; Fernández et al., 2019; Gutiérrez et al., 2020) y el deterioro de propiedades físicas, químicas y biológicas edáficas en agroecosistemas (Gutiérrez et al., 2020; Vallejo-Quintero et al., 2018).
La calidad del suelo (CS) en agroecosistemas ha sido evaluada y monitoreada a través de herramientas (cualitativas y cuantitativas) que incluyen, generalmente, la selección de un conjunto mínimo de indicadores (CMD) mediante modelos o procedimientos estadísticos multivariados o mediante la opinión de expertos, para finalmente realizar su respectiva puntuación e integración en índices de calidad (Maurya et al., 2020; Rezaee et al., 2020). Las evaluaciones de CS deben considerar indicadores físicos, químicos y biológicos, asociadas con importantes procesos del suelo, de tal modo que su interpretación se traduzca en una visión holística del estado general del suelo, permitiendo identificar fenómenos de degradación temprana, así como corroborar el mantenimiento o aumento de la CS en sistemas agrícolas sometidos a distintas prácticas de manejo (Afanador et al., 2020; Maurya et al., 2020; Yu et al., 2018).
Estudios llevados a cabo por Pinzón-Pinto et al. (2016) en el Páramo de Guerrero evidenciaron que la agricultura y la ganadería han contribuido a una pérdida y fragmentación de la vegetación natural y a la eliminación del horizonte genético O, lo que afecta la calidad física (retención de humedad, densidad aparente, temperatura del suelo y porosidad) y química de sus suelos (disminución del contenido de materia orgánica y nutrientes). Igualmente, estudios desarrollados por Álvarez-Yela et al. (2017) y Avellaneda-Torres et al. (2018) en páramos del Parque Nacional Natural de los Nevados (Colombia) demostraron cómo la agricultura intensiva afecta la dinámica y la calidad del suelo, generando un detrimento en la diversidad y función de grupos microbianos indispensables en la dinámica de los ciclos biogeoquímicos (C y N). Se evidenció que las actividades agropecuarias alteran la estructura de la comunidad microbiana y ocasionan una pérdida de la diversidad funcional metabólica, reduciendo su resistencia al estrés ambiental. Igualmente, se observó un aumento significativo de arqueas oxidadoras de amonio, bacterias copiotróficas y hongos ascomicetos, que afectan la respuesta fisiológica y funcional del ecosistema en las áreas intervenidas. Así mismo, se encontró impacto en indicadores químicos (carbono orgánico, Ca, K, capacidad de intercambio catiónico y amonio) y actividades enzimáticas (fosfatasas, fosfodiesterasas, β-glucosidasas y proteasas).
No obstante, a pesar de las importantes funciones ecológicas que ejercen este tipo de ecosistemas, continúan siendo limitados e insuficientes los estudios que evalúen el efecto de las actividades agropecuarias sobre la CS en zonas de páramo; este tipo de estudios resultan necesarios para plantear alternativas de mitigación, recuperación y conservación del recurso del suelo. Por otro lado, el Páramo de Guerrero es el segundo en extensión, en Colombia, después del de Sumapaz, es el más intervenido por actividades antrópicas (cultivos de papa, ganadería y explotación de carbón), y presenta distintos fenómenos de degradación que han puesto en peligro la prestación de los servicios ecosistémicos relacionados con la regulación del ciclo hidrológico y el ciclaje de nutrientes.
Por ello, el objetivo del presente estudio fue evaluar el efecto del monocultivo de papa sobre la CS en el Páramo de Guerrero, para lo cual se plantearon los siguientes objetivos específicos: (1) seleccionar un CMD y un ICS para evaluar el efecto generado por el monocultivo de papa sobre la CS y (2) evaluar la calidad física edáfica en agroecosistemas del Páramo de Guerrero a través de una guía diagnóstica de calidad en campo.
Materiales y métodos
Área de estudio
El Páramo de Guerrero se localiza en la cordillera Oriental de Colombia y ocupa un área de 15 000 ha, que comprenden una franja de montaña a una altitud de 3000 a 3900 m s.n.m. (Pinzón-Pinto, 2016). Esta investigación se llevó a cabo en el municipio de Zipaquirá (Cundinamarca) en la vereda Páramo de Guerrero oriental y occidental (05°07´34,3¨ N, 74°03´6,4¨ W).
La precipitación es de régimen bimodal, con dos periodos húmedos (abril, mayo y parte de junio) y dos secos, el primero en diciembre, enero, febrero y marzo, y el segundo en julio, agosto y parte de septiembre, con variaciones en los valores promedios multianual entre 604.9 y 1084.5 mm (Estación Carmen de Carupa) (Instituto de Investigación de Recursos Biológicos Alexander Humboldt, s. f.). El Páramo de Guerrero se encuentra entre un régimen semihúmedo y húmedo, y su temperatura oscila entre 7.2 y 11.1º C (Estación Guerrero). Se seleccionó una finca representativa (15 ha) constituida por monocultivos de papa, relictos de bosque con vegetación de páramo y zonas de descanso agrícola.
Diseño experimental y muestreo de suelos
Para el muestreo se seleccionaron tres agroecosistemas: (1) monocultivo convencional de papa (Solanum tuberosum L.) (SC) (05°07´34,3¨ N, 74°03´6,4¨ W) de 4 años, bajo labranza convencional y aplicación periódica de agroquímicos; (2) un suelo en descanso (18 meses) (SD) (05°07´33,6¨ N, 74°03´6,5¨ W), con adiciones de semilla de avena y; (3) un relicto de bosque con vegetación de páramo (BO) (05°06´32,7¨ N, 74°02´6,1¨ W). En cada agroecosistema se establecieron dos cuadrantes (3 x 3m), en los cuales se tomaron con barrenos metálicos dos muestras compuestas aleatorias (n=2) provenientes de 25 submuestras (profundidad de 0-15 cm). Cada muestra de suelo se almacenó en bolsas ziploc debidamente rotuladas, las cuales se dividieron en dos porciones: una parte fue tamizada (2 mm) y secada a temperatura ambiente, para análisis fisicoquímicos y la otra parte fue refrigerada (4° C) para análisis microbiológicos en los laboratorios de suelo de la Universidad Central.
Guía diagnostica in situ para la evaluación de la calidad física
Como parte de la evaluación de la calidad de suelos en el Páramo, se implementó la guía de calidad física del suelo desarrollada por Vallejo-Quintero et al. (2021). La guía considera indicadores en campo (velocidad de infiltración, resistencia a la penetración, estructura del suelo, consistencia del suelo, color, porosidad del suelo, resistencia al rompimiento), los cuales cuentan con una puntuación (menos deseable = 0, moderada = 1, preferido = 2) en función de su desempeño.
Métodos fisicoquímicos y microbiológicos
Se seleccionó un CTD de 14 propiedades físicas (textura, densidad real y aparente, porosidad, humedad y agregados estables en seco: índice de estabilidad de agregados, agregados estables y diámetro geométrico medio); químicas (pH, carbono orgánico y conductividad eléctrica); microbiológicas (densidad de hongos y heterótrofos totales); y bioquímicas (actividad de catalasa), cuyos métodos fueron descritos por Vallejo-Quintero et al. (2018) y por Afanador et al. (2020).
Análisis estadístico
Se realizó una prueba estadística de Shapiro-Wilk para determinar la distribución normal de los datos. Se empleó el análisis de varianza (ANOVA) (P≤0.05) para comparar las 14 variables fisicoquímicas y microbiológicas seleccionadas, así mismo, se realizó un análisis de comparaciones múltiples (Tukey) para identificar diferencias significativas (P≤ 0.05) entre los agroecosistemas mediante el software estadístico SPSS v.s 25 ®.
Selección del conjunto mínimo de datos (CMD)
Para la selección del CMD se realizó un análisis de componentes principales (ACP) con las 14 variables estudiadas, empleando el programa estadístico Paleontological Statistics (PAST), versión 2.16, con una matriz de datos estandarizada (matriz de correlación). Para la conformación del CMD, se seleccionaron los componentes principales (CP) que tuvieron valores propios (eigenvalues) ≥1, así como aquellos que explicaron al menos el 5 % de la varianza total del CTD. En cada CP se extrajeron las variables que presentaron un peso (valor de carga) ≥ 0.65. El CMD quedó conformado por aquellas variables que mejor explicaron la variabilidad total con respecto a todas las propiedades inicialmente seleccionadas (Afanador et al., 2020).
Índice de calidad de suelo (ICS)
Después de determinar las variables del CMD se calculó el ICS y cada variable fue transformada y normalizada en valores sin unidad entre 0 y 1, empleando funciones lineales. La técnica de indexación utilizada fue el IC aditivo ponderado propuesto por Karlen y Stott, (1994), la cual ha sido ampliamente utilizada por Chen et al. (2013); Nabiollahi et al. (2017); y Yu et al., (2018). El procedimiento detallado del cálculo del ICS es descrito por Afanador et al. (2020). Las calificaciones previamente establecidas para los indicadores se integraron en la siguiente ecuación (1):
Resultados y discusión
Efecto de las prácticas de manejo agrícolas en las propiedades de suelo
De acuerdo con los resultados obtenidos a partir de la guía diagnóstica de calidad física, se evidenció que el SC tuvo valores bajos con relación a las variables: resistencia al rompimiento, capas endurecidas y velocidad de infiltración (Tabla 1), y presentó, en términos generales, una menor calidad física del suelo (calificación de 22) en comparación con SD y BO (Tabla 1).
Los resultados evidenciados en el SC, con relación al comportamiento de dichas variables físicas en campo, se asocian a las prácticas de manejo implementadas por décadas en los cultivos de papa, principalmente, la labranza convencional y el limitado aporte de material orgánico al suelo, lo cual impacta la condición física de este. Lo anterior demuestra, al igual que investigaciones desarrolladas por Cherubin et al. (2019), la importancia de implementar e incluir este tipo de guías en campo, dada su fácil ejecución e interpretación y rápida determinación, pues permiten identificar señales tempranas de degradación edáfica y monitorear procesos de recuperación de suelos, cuyos efectos, en algunas ocasiones, no son rápidamente identificables por propiedades medibles en laboratorio.
Por otra parte, del CTD de 14 propiedades edáficas evaluadas, se evidenciaron diferencias significativas (ANOVA P≤0.05) con respecto a las variables: humedad (H), pH, carbono orgánico (CO), densidad de heterótrofos (HET) y actividad enzimática de catalasa (CAT) (Tabla 2).
Tabla 2 Propiedades fisicoquímicas y microbiológicas evaluadas en los agroecosistemas
![](/img/revistas/acag/v71n3//0120-2812-acag-71-03-280-gt2.png)
*IE= Índice de Estabilidad; EA= Estabilidad de Agregados; DGM= Diámetro Geométrico Medio; DPM= Diámetro Ponderado Medio; H= Humedad; DA= Densidad Aparente; PO= Porosidad; RP= Resistencia a la Penetración; CE= Conductividad Eléctrica; pH; CO= Carbono Orgánico; HET= Heterótrofos; HOT= Hongos Totales; CAT= Catalasa. Letras diferentes indican diferencias significativas entre los agroecosistemas (Tukey, P≤0.05).
En el presente estudio, la H fue significativamente mayor en el BO (32.68 %), en comparación con el SC y SD (Tabla 2) (Tukey P≤0.05). Los resultados concuerdan con Daza et al. (2014) y Camargo-García et al. (2012), quienes evidenciaron valores de H en un rango entre 26.8 % y 37.0 % en suelos del Páramo de Sumapaz con vegetación nativa y de 29.5 a 36.2 % en páramos del Parque Nacional Natural de Los Nevados. Los mayores porcentajes de H encontrados en el BO pueden asociarse al origen (glaciar y volcánico) de los suelos de páramos en Colombia, que están caracterizados por contener valores altos de materia orgánica (MO) que se mineraliza lentamente en bajas temperaturas. La combinación MO-cenizas volcánicas y el tipo de vegetación (musgos, líquenes y frailejones) que contribuye a la formación de una espesa y permanente cobertura orgánica del suelo, junto con las condiciones climáticas en estos ecosistemas, favorecen las condiciones de humedad del suelo y su capacidad de retención de agua (Vargas, 2013).
Por otra parte, los bajos valores de H del SC concuerdan con lo reportado por Farfán et al. (2020) en suelos bajo cultivo de papa en el páramo neotropical alto andino (Colombia); lo cual fue asociado al ingreso limitado de material orgánico al suelo en los monocultivos, de forma tal que el suelo permanece expuesto a la radiación solar, incrementando la evaporación en la superficie y limitando su capacidad para retener humedad.
Con relación a los valores de pH, estos fueron significativamente menores en el BO (3.66) (Tabla 2) (Tukey P≤0.05), y fueron similares a los reportados por Díaz-Granados et al. (2005) y Fernández et al. (2019). Los autores destacan que los suelos de páramos colombianos son generalmente de origen volcánico y se caracterizan por ser húmedos y ácidos, con valores de pH entre 3.9 y 5.4. El carácter ácido predominante en sus suelos favorece la formación de complejos de humus-Al; complejos que predominan en el horizonte A de la mayoría de los páramos colombianos de origen volcánico. Avellaneda-Torres et al. (2018) y Pinzón-Pinto (2016) reportaron valores de pH de 5.3 en suelos de ecosistemas neotropicales de páramos altoandinos en el Parque Natural Nacional de los Nevados y de 4.2 a 4.8 en horizontes genéticos O y A de suelos con vegetación natural en el Páramo de Guerrero (Colombia), respectivamente.
Con respecto al porcentaje de CO, este fue significativamente mayor en SD (9.32 %) (Tukey P≤0.05) en comparación con el SC y BO (Tabla 2). Los mayores valores de CO podrían asociarse con las actividades implementadas durante el periodo de descanso y recuperación del suelo, específicamente, la adición de la semilla de avena, lo cual contribuye a una mejoría en indicadores específicos del suelo. Los porcentajes de CO encontrados en el BO (6.61 %) son similares a los reportados por Fernández et al. (2019)) en suelos con vegetación nativa en el Páramo de Rabanal (Boyacá), con valores de CO entre 6.59 y 7.54 %.
Con relación a la densidad de HET, el SC presentó valores significativamente mayores (P≤0.05) en el SC (4.11 Log UFC/gps) (Tabla 2). Lo anterior podría asociarse con las prácticas de manejo implementadas en el monocultivo de papa, particularmente con las operaciones de labranza y los agroquímicos adicionados al suelo. Es bien conocido que la labranza a baja escala incrementa la aireación y acelera la descomposición de los residuos orgánicos, ocasionando un aumento en la disponibilidad de nutrientes que favorece la densidad de microorganismos edáficos. Avellaneda-Torres et al. (2020) evidenciaron cambios en la estructura y composición de las comunidades microbianas en suelos de páramos (Parque Nacional de los Nevados) cultivados con papa y bajo ganadería, sin embargo, los cambios dependieron de la elevación y las condiciones climáticas, siendo mayor la riqueza microbiana en los suelos cultivados.
Contrariamente, otros autores han reportado que la labranza convencional genera un detrimento en la abundancia de microorganismos del suelo, lo cual estaría asociado al disturbio físico del suelo, la compactación, la variación en las condiciones de humedad y la disponibilidad de nutrientes, así como al intercambio gaseoso; ello, conjuntamente, altera su densidad y actividad (Choudhary et al., 2018; Moreno et al., 2021). No obstante, en el presente estudio no hubo diferencias significativas en la mayoría de las variables físicas evaluadas en laboratorio (IE, EA, DGM, DPM, DA, RP y PO), las cuales se relacionan con una condición física y con el grado de compactación del suelo. Farfán et al. (2020) evidenciaron una disminución en la densidad de bacterias y hongos celulíticos en suelos de páramo neotropical alto andino (Colombia) disminución que está asociada con la aplicación de insumos químicos agrícolas.
Por su parte, la actividad enzimática de CAT presentó un comportamiento opuesto a la densidad de HET, siendo significativamente menor en el SC (0.16) (Tabla 2). Lo anterior podría relacionarse con la aplicación de pesticidas a largo plazo, lo que genera un impacto negativo en las actividades enzimáticas del suelo (Farfán et al., 2020). La actividad de CAT suele emplearse como un importante indicador de calidad edáfica debido a su alta sensibilidad a los cambios en el uso y manejo del suelo (Romeo et al., 2020; Vallejo-Quintero et al., 2018). En el presente estudio, los menores valores de H (20.81 %) y de CO (5.18 %) encontrados en el SC podrían haber afectado dicha actividad; este resultado es similar a lo reportado por Avellaneda-Torres et al. (2018).
Selección del CMD para la evaluación de la CS
La selección de indicadores para el CDM se realizó a partir de un CTD de 14 indicadores a través de PAST, de los cuales quedaron11 variables: DA, PO, H, CO, RP, IE, EA, DPM, DGM, HET y HOT (Tabla 3). El ACP indicó que el 72.381 % de las diferencias evidenciadas entre las propiedades estudiadas pueden ser explicadas en un 32.54 % por el componente uno, un 23.72 % por el componente dos y un 16.12 % por el componente tres. Los tres primeros CP se seleccionaron al presentar valores propios ≥1. De los tres CP seleccionados, solo los parámetros con mayores valores de carga (≥0.65) se reunieron para la indexación. Para el CP1 se seleccionaron: PO, H, EA, DPM. DGM y HOT; para el CP2: CO, RP y HET; y para el CP3: DA y EA (Tabla 3).
De los indicadores seleccionados para el CDM, específicamente, H, CO, RP, EA, DPM, DGM, IE, PO, HET y HOT concuerdan con los reportados por Afanador- Barajas et al. (2020) y Vallejo-Quintero et al. (2018) en agroecosistemas en el municipio de Cachipay (Cundinamarca), quienes seleccionaron 6 de los 11 indicadores trabajados en el presente estudio (EA, DPM, DGM, H y HET). El CMD seleccionado en este estudio puede ser utilizado para hacer un seguimiento periódico de los cambios en la CS en áreas protegidas, como páramos, que han sido afectadas por actividades antrópicas, en lugar de utilizar un CTD, ya que ahorra tiempo y gastos en programas de monitoreo.
Comparación de la calidad del suelo mediante el IC
Con relación al ICS, este fue mayor en el BO (1.11) y SD (1.10) y, por último, el SC (1.07), sin embargo, no hubo diferencias significativas entre ellos (Figura 1). Las variables que tuvieron mayor contribución al ICS fueron DPM (20.08 %), DGM (20.08 %), CO (18.31 %), EA (13.22 %), HOT (12.60 %), HET (9.23 %), PO (9.11 %), IE (7.45 %), RP (6.92 %), DA (5.61 %) y H (0.23 %) (Figura 1). La tendencia evidenciada en el ICS es consistente con Zou et al. (2021), quienes reportaron menores valores del ICS en monocultivos de caucho (Hevea brasiliensis) en comparación con sistemas agroforestales y bosques, lo cual está asociado a una significativa pérdida de nutrientes del suelo.
![](/img/revistas/acag/v71n3//0120-2812-acag-71-03-280-gf1.png)
Figura 1. Valores del índice de calidad del suelo en los agroecosistemas comparados. *SC= suelo cultivado, SD= Suelo en descanso, BO= Bosque
Conclusiones
El CMD seleccionado en el presente estudio permitió identificar 11 variables fisicoquímicas y microbiológicas (DA, PO, H, CO, RP, IE, EA, DPM, DGM, HET y HOT) altamente sensibles a las prácticas agrícolas del páramo. A su vez, el IC fue menor en el SC, en comparación con el SD y B, además, es relevante que estos dos últimos tuvieron índices muy similares, lo que refleja una progresiva recuperación del suelo previamente cultivado; posiblemente dicha recuperación se debe a la adición de materiales orgánicos al suelo, como la semilla de avena. Por su parte, la guía diagnóstica de calidad del suelo permitió identificar señales físicas de degradación del suelo, las cuales no fueron evidenciadas a través de los indicadores físicos ejecutados en laboratorio.