Introducción

El suelo es un recurso natural fundamental en los agroecosistemas, con múltiples funciones ecológicas dentro de las cuales destacan su rol conectivo con la atmosfera y su influencia sobre la biosfera; su capacidad de almacenamiento de energía y materia; en adición a su papel como un enorme reservorio de carbono (^{Hillel 1998}, ^{Lal 2012}, ^{Hartge et al. 2016}).

El suelo es recurso de suma importancia a nivel agroecológico y económico, el cual es altamente vulnerable a la degradación y la pérdida irreversible del mismo (^{Lal 2001}, ^{Bertsch 2006}), tanto por procesos naturales (como el potencial erosivo de la lluvia, la gravedad y el viento) (^{Vahrson 1990}, ^{Villatoro-Sánchez et al. 2015}) como por procesos asociados a factores antrópicos, entre ellos, la labranza mecanizada (Lal 2019). En ese sentido, la labranza mecanizada (producto de sistemas agrícolas de producción intensiva) es un factor de suma importancia debido a su influencia directa sobre la magnitud de la pérdida de suelo que se asocia a los otros factores naturales anteriormente mencionados (^{Mehuys et al. 2009}).

Dichos procesos erosivos ocurren comúnmente en suelos derivados de cenizas volcánicas (Andisoles), en especial en aquellos que se desarrollan en las posiciones de ladera y las cimas de los sitios elevados (^{Shoji et al. 1993}, ^{Alvarado et al. 2001}). Estos suelos se extienden a lo largo de los países de Centro América, principalmente los que se localizan en las cercanías de los numerosos volcanes de la región (Alvarado et al. 2001, Alvarado 2022).

De forma general, los Andisoles son descritos en la literatura como suelos resistentes a los procesos erosivos, debido a sus propiedades físicas como altos contenidos de materia orgánica y materiales nano cristalinos, altas porosidades y tasas de infiltración que evitan los flujos erosivos por precipitación pluvial (^{Shoji et al. 1993}, ^{Zehetner y Miller 2006}).

Sin embargo, bajo las condiciones climáticas tropicales imperantes en la región Centroamericana (^{Solano y Villalobos 2001}, ^{Alfaro 2002}) las tasas de erosión reportadas para los Andisoles (tanto cuantificadas de forma experimental como calculadas con la ecuación universal de pérdida de suelo) pueden llegar a valores de hasta 76 Ton ha^-1 año^-1 (^{Lianes et al. 2009}, ^{Cerretelli et al. 2023}); valores considerados como erosión severa según la clasificación propuesta por FAO (1979). Curiosamente, estas elevadas tasas de erosión incluyen principalmente el poder erosivo de la lluvia como agente causal, sin contemplar otros posibles mecanismos de erosión y remoción del suelo asociados a las actividades agrícolas, entre ellos, erosión por labranza mecanizada y la pérdida de suelo por cosecha del cultivo (Soil loss due to crop harvesting; SLCH).

La erosión por labranza mecanizada consiste en la masa de suelo que es removida y transportada por el implemento de labranza (apero) desde un sector del sembradío a otro durante las labores de mecanización y preparación del suelo (^{Govers et al. 1999}). Este tipo de erosión depende en gran medida del valor de la pendiente del terreno y la dirección de la mecanización (es decir si esta se hace de forma perpendicular o paralela o a la pendiente), donde valores elevados de pendiente (15-35%) y con mecanización realizada a favor de la pendiente pueden acarrear cantidades significativas de sedimento (^{Mehuys et al. 2009}). En Andisoles de la zona alta de Cartago, se han reportado tasas de erosión por labranza mecanizada de hasta 250 Ton ha^-1 ciclo^-1, cuando la mecanización se realiza a favor de pendiente con arado de disco (^{Tiessen et al. 2010}).

Por otro lado, la pérdida de suelo por cosecha del cultivo (SLCH) es probablemente el mecanismo con menor reconocimiento dentro de los estudios de erosión y pérdida de suelo (^{Kuhwald et al. 2022}). SLCH consiste en la pérdida de sedimentos, los cuales salen adheridos a la parte comercial del cultivo y salen de forma irreversible del sitio de siembra (^{Ruysschaert et al. 2004}, ^{Mwango et al. 2015}, ^{Oshunsanya et al. 2018}, Kuhwald et al. 2022). Este mecanismo puede generar tasas de erosión de hasta 22 Ton ha^-1 ciclo^-1, y depende en gran medida de la morfología del cultivo, la densidad de siembra y el tipo de suelo en el que se establece el cultivo (Ruysschaert et al. 2004, ^{Han-Qing et al. 2016}, Kuhwald et al. 2022).

Los estudios sobre SLCH son bastante escasos (^{Kuhwald et al. 2022}), y la gran mayoría de estos trabajos se centran en la pérdida de sedimentos. Un componente de suma importancia de la remoción de sedimentos que ocurre por este mecanismo es la pérdida irreversible de materia orgánica y elementos esenciales para el crecimiento vegetal que salen del sitio de siembra cada vez que se cosecha el cultivo (^{Mwango et al. 2015}, ^{Han-Qing et al. 2016}).

A nivel Centroamericano, no se tienen reportes en la literatura de estudios sobre erosión de suelo que cuantifiquen las pérdidas asociadas a SLCH, y mucho menos las cantidades de nutrimentos y la materia orgánica que se pierden de forma irreversible del sitio de siembra. Para atender esta problemática se planteó la presente investigación, la cual tuvo como objetivo cuantificar la pérdida de suelo por cosecha del cultivo (SLCH) en dos ciclos de cultivo de cebolla (Allium cepa L.) bajo tres implementos de labranza mecanizada en Andisoles de las tierras altas de Cartago, Costa Rica.

Materiales y métodos

Sitio de estudio

El presente estudio se desarrolló en una finca dedicada a la producción de los cultivos de cebolla (Allium cepa L.) y papa (Solanum tuberosum L.) ubicada en el distrito de Tierra Blanca, cantón de Cartago (9°54’20” N; 83°52’24” O) a una elevación de 1973 msnm (Figura 1). Los suelos de la finca se clasificaron dentro del orden de los Andisoles, con un régimen de humedad ústico (Suborden Ustands), comúnmente clasificados dentro del Subgrupo taxonómico Typic Haplustands (^{Alvarado et al. 2001}, ^{Alvarado y Mata 2016}).

Figura 1 Vista panorámica de la finca dedicada al cultivo de cebolla (Allium cepa L.) donde se estableció el presente estudio. Flechas blancas señalan preparación de las camas mediante un arado de vertedera tirado por un caballo (A). Labores de mecanización con la rastra rotativa (RR), uno de los implementos evaluados (B). 

La pendiente promedio del terreno (medida con un clinómetro Suunto^® PM-5) fue de 10%, y el suelo en su horizonte superficial Ap presentó una textura franco-arcillosa, un contenido de materia orgánica promedio de 2,9%, un valor de densidad aparente de 1,13 Mg m^-3 y un 52% de espacio poroso total antes del inicio del experimento. Más detalles sobre el sitio experimental se pueden consultar en el trabajo de ^{Marín-Araya y Villatoro-Sánchez (2023}).

Mecanización, siembra y manejo del cultivo

El experimento inició en mayo del 2019, cuando se sembró el primer ciclo de cebolla variedad “Álvara” durante el primer semestre del año, mientras que el segundo ciclo se sembró en abril del 2020 manteniendo la misma variedad.

Las parcelas establecidas con el cultivo de la cebolla recibieron tres fertilizaciones con fórmula N - P₂O₅ - K₂O distribuidas según la etapa fenológica: 1) 12-24-12 30 días después del trasplante (ddt), 2) 12-11-18 a los 60 ddt, y 3) 19-4-19 a los 90 ddt, para un total de 173 kg de N, 59 kg P y 164 kg de K por hectárea por ciclo. También se realizaron prácticas de control de malezas mediante la aplicación de herbicidas premergentes: Pendimetalina (2 L ha^-1) aplicado 2 días después del trasplante; Linurón (1 kg ha^-1) y Oxiflouren (200 mL ha^-1) ambos aplicados a los 30 días después del trasplante.

La labranza convencional para el cultivo de la cebolla consistió en un subsolado a 50 cm de profundidad en dirección transversal a la pendiente. Dicho subsolador se compuso de 7 picos, con 320 mm de espaciado entre picos, 1920 mm de largo y un peso de 530 kg, el cual se pasó a una velocidad de avance de 6 km h^-1. Posteriormente se utilizó un arado de disco a favor de pendiente, con una velocidad de pase de 10 km h^-1, y seguidamente se utilizó el rotador, a favor de pendiente para afinar el suelo. Estas operaciones de labranza fueron realizadas con un tractor Lamborghini^® R3.105, con un peso de 4067 kg y una potencia máxima de 102 CV y 2300 rpm. El tractor no usaba ningún tipo de lastre en las llantas, taco de tipo R1 y presión de 40 psi. Finalmente se elaboraron los surcos por medio de un arado de vertedera de tipo universal tirado por un caballo (Figura 1A).

Se establecieron tres tratamientos de labranza mecanizada, donde se mantuvo el subsolado y arado inicial, pero reemplazando el rotador (RD) por una rastra rotativa (RR) o un palín mecánico (PM). Se mantuvo un tratamiento con rotador (el cual sirvió como control) para comparar con la labranza común que desarrolla el productor. Una descripción más detallada sobre los implementos de labranza empleados se puede encontrar en el trabajo realizado por ^{Marín-Araya y Villatoro-Sánchez (2023}).

El área de cada repetición (una parcela de 240 m²) fue cubierta por 2 a 3 pases adyacentes del implemento respectivo, la cantidad dependió del ancho de pase de cada uno de los implementos. En este caso, PM solamente necesitó 2 pases, mientras que RR y RD requirió de 3 pases.

Diseño experimental y muestreo

El diseño experimental utilizado en el presente ensayo fue de bloques completos al azar, con tres repeticiones para cada tratamiento de mecanización, con una parcela de 240 m² como una unidad experimental. Las muestras para cuantificar la pérdida de sedimento como consecuencia de la cosecha del cultivo fueron colectadas al final de cada ciclo del cultivo de la cebolla (en el momento de la cosecha), y consistió en tomar el sedimento pegado alrededor de los bulbos de las cebollas cosechadas (Figura 2). Para esto se cosecharon 12 cebollas a lo largo de 1,25 m lineales, se sacudieron y el sedimento remanente fue colectado en una bolsa con la ayuda de una brocha (Figura 2). En este caso, se tomaron tres puntos al azar por parcela (para un total de 36 cebollas evaluadas) y se confeccionó una muestra compuesta por cada repetición (para un total de 9 muestras compuestas por ciclo). Para cada muestra recolectada, se determinó su peso en campo (g) mediante una balanza granataria portátil, y posteriormente se determinó su porcentaje de humedad gravimétrica (w) siguiendo los lineamientos de ^{Topp y Ferré (2002}) en el laboratorio de Recursos Naturales del Centro de Investigaciones Agronómicas de la Universidad de Costa Rica (CIA-UCR).

Figura 2 Proceso de cuantificación de los sedimentos que se pierden por labores de cosecha en el cultivo de cebolla (Allium cepa L.) obtenidos en un Andisol de Tierra Blanca, Cartago. 

Las muestras compuestas fueron procesadas en el laboratorio de suelos y foliares del CIA-UCR. Se realizó un análisis químico completo (pH, acidez, Ca, Mg, K, P, Zn, Cu, Fe, Mn) donde se utilizó la solución extractora KCl - Olsen modificado (^{Díaz-Romeu y Hunter 1978}), y se determinó el contenido de nitrógeno y carbono orgánico (CO) mediante un auto analizador de C por combustión seca (^{Horneck y Miller 1998}).

El dato de materia orgánica (MO) se calculó a partir de los valores de C, multiplicando dichos valores por una constante de 1,43 para obtener el resultado de MO.

Adicionalmente, mediante muestras no alteradas (colectadas en cilindros) se determinó la densidad aparente del suelo en cada una de las unidades experimentales siguiendo los lineamientos de ^{Grossman y Reinsch (2002}).

Con los datos de peso de muestra en campo (M _w ) y su correspondiente porcentaje de humedad gravimétrica (w), se determinó la masa de sólidos (M _s ) mediante la siguiente ecuación:

Dichos valores de masa de sólidos permitieron calcular la masa de suelo que se pierde por la cosecha de la cebolla, y tomando en cuenta la densidad de plantas y el área de las parcelas, se pudo calcular la tasa de pérdida de sedimentos por ciclo (Ton ha^-1 ciclo^-1).

También, se calculó la cantidad de elementos (N, P, K, Ca, Mg, Fe, Cu, Zn, Mn) y MO que salen del sitio de siembra como sedimentos adheridos al bulbo, a partir de los valores de obtenidos de análisis químicos y la cantidad de sedimento mediante los siguientes procedimientos:

1) Se transformó la masa de sedimento (kg ha^-1) en volumen (L ha^-1) con la correspondiente densidad aparente del suelo (Mg m^-3 o 1000 kg m^-3) mediante la siguiente igualdad:

Esto con el objetivo de facilitar los cálculos, ya que los análisis químicos de suelo que se hacen de rutina en Costa Rica (^{Corrales et al. 2005}) expresan los contenidos de elementos por unidad de volumen (cmol₍₊₎ L^-1, mg L^-1) o bien como un porcentaje total.

2) Si la concentración del elemento fue expresada en cmol₍₊₎ L^-1 (Ca, Mg, K) se calculó la cantidad del elemento mediante su correspondiente masa equivalente (ME _(g) ), que resultó de dividir su masa molar por la carga correspondiente. El cálculo final se realizó de esta forma:

3) Si la concentración (P, Fe, Cu, Zn, Mn) se expresó en unidades de ppm (mg L^-1) se calculó la cantidad de elemento mediante la siguiente igualdad:

4) Si la concentración (N, MO) se expresó en unidades de porcentaje (%) se calculó la cantidad de elemento mediante la siguiente igualdad:

Análisis estadístico

Los resultados obtenidos fueron analizados mediante un modelo lineal generalizado (GLM), con una estructura factorial que incluyó los factores: Ciclo de cultivo (C), Implemento de labranza (I) y una interacción entre factores (C x I), con el bloque como factor aleatorio. Para la separación de medias se utilizó el método HSD-Tukey con un p-value<0,05. Se verificaron los supuestos de normalidad y homocedasticidad de varianzas a través de la inspección visual de la distribución de los residuos obtenidos de los modelos ajustados (normal QQ Wally plots y Wally plots) para las variables de respuesta evaluadas siguiendo los lineamientos descritos por ^{Kozak y Piepho (2018}).

Todos estos procedimientos anteriores se realizaron mediante el programa estadístico R Studio (R v. 4.0.4, 2021-02-15) “Lost Library Book” (^{R Studio Team 2015}).

Resultados y discusión

Pérdida de suelo por cosecha del cultivo (SLCH) de cebolla en Andisoles

Un componente que se ha estudiado muy poco dentro de los procesos de erosión y pérdida de suelo en es la masa de sedimentos que se sale del suelo por la cosecha de cultivos hortícolas. Los resultados obtenidos para la pérdida de suelo por cosecha del cultivo (SLHC por sus siglas en inglés) en función de los implementos de labranza (Figura 3), muestran una tendencia a la reducción en la cantidad de sedimento con la rastra rotativa (0,60 ±0,47 Ton ha^-1), seguida por el palín mecánico (0,87 ± 0,55 Ton ha^-1), y el rotador (1,05 ± 0,46 Ton ha^-1) siendo el tratamiento que registró mayor cantidad de sedimento adherido al bulbo, sin diferencias significativas entre los tres tratamientos evaluados (Tabla 1; Figura 3).

Figura 3 Pérdida de sedimentos por la cosecha del cultivo de cebolla, como respuesta al efecto de tres implementos de labranza en ciclos de producción no continuos en un Andisol de Tierra Blanca, Cartago. Barras de error como error estándar. * Sin diferencias significativas (p-value>0,05) entre las medias por los factores evaluados en el modelo (Tabla 2). 

Una posible causa de esta tendencia fue un problema biótico que se presentó en un sector de la plantación. Debido al diseño del ensayo, una repetición con la rastra rotativa y otra repetición con el palín mecánico coincidieron con un sector donde las plantas presentaron clorosis, un crecimiento reducido y mortalidad (Figura 4). Esta condición afectó el tamaño de los bulbos y de forma indirecta, la cantidad de sedimento colectado en dichas unidades experimentales, en concordancia con lo reportado por ^{Isabirye et al. (2007}) donde menores rendimientos del camote (Ipomoea batatas L.) representaron menores cantidades de sedimentos perdidos como SLCH.

Figura 4 Daño abiótico en el cultivo de cebolla (Allium cepa L.), que generó reducción de crecimiento, mortalidad de las plantas y espacios vacíos en las parcelas evaluadas en el presente estudio. Tierra Blanca, Cartago. 

En el segundo ciclo de cultivo, la menor cantidad de sedimento colectado se obtuvo con el palín mecánico (0,94 ± 0,23 Ton ha^-1), seguido por la rastra rotativa (0,99 ± 0,26 Ton ha^-1) y el rotador (1,08 ± 0,08 Ton ha^-1). Tampoco se presentaron diferencias significativas entre los implementos de labranza (p > 0,05) para este ciclo de cultivo (Tabla 1; Figura 3).

En un estudio realizado en las tierras altas de Tanzania (1400-1600 msnm), ^{Mwango et al. (2015}) reportaron pérdidas de sedimentos por la cosecha de cebolla de hasta 3,0 Ton ha^-1 ciclo^-1 cultivadas en Entisoles y Ultisoles, mientras que ^{Antoneli et al. (2021}) obtuvieron valores de 3,35 Ton ha^-1 ciclo^-1 en Ultisoles de las montañas de Papuã en Paraná, Brasil. Todos estos valores mencionados fueron superiores a los máximos obtenidos en el presente estudio (1,05-1,08 Ton ha^-1 ciclo^-1) en los Andisoles de Tierra Blanca de Cartago. Estas diferencias se deben probablemente a un menor contenido de arcilla en texturas francas o franco arenosas características de los Andisoles (^{Alvarado et al. 2001}), las cuales exhiben una menor expresión de las fuerzas de adhesión y cohesión de las partículas del suelo (^{Nanzyo et al. 1993}), y por ende una menor cantidad de sedimento adherido al bulbo de la cebolla en comparación con suelos de texturas arcillosas (Antoneli et al. 2021).

Aunque no se encontró diferencias significativas entre los implementos de labranza evaluados, los valores obtenidos de SLCH fueron similares a los resultados de otros trabajos en donde la labranza y el cultivo presentan un efecto importante sobre la cantidad final de sedimento perdido por labores de cosecha, incluyendo el sedimento adherido al producto (^{Ruysschaert et al. 2004}, ²⁰⁰⁶, y 2007, ^{Thomaz y Bereze 2021}, ^{Kuhwald et al. 2022}).

Una mayor alteración del suelo por labranza mecanizada podría tener un efecto en la cantidad de sedimento adherido al bulbo de la cebolla, el cual se podría atribuir a dos posibles causas: 1) un incremento en la densidad de los pelos radicales, y 2) un aumento en el rendimiento del cultivo (^{Oshunsanya et al. 2018}). En ese sentido, la labranza modifica la estructura del suelo e incrementa la porosidad al separar los agregados, creando espacios que facilitan el crecimiento y desarrollo radicular (^{Büchi et al. 2017}). Dicho incremento en la cantidad de raíces permite una mayor cantidad de partículas de suelo adheridas al cultivo (Oshunsanya et al. 2018); lo cual no se vio reflejado en los resultados obtenidos en presente estudio.

Pérdida de nutrimentos por la cosecha del cultivo de cebolla

Uno de los mayores impactos de la erosión del suelo es la pérdida irreversible de elementos esenciales para el crecimiento vegetal (^{Bertsch 1995}, ^{Alewell et al. 2020}). Los sedimentos que salen del sembradío llevan consigo una cantidad importante de nutrimentos, los cuales pueden encontrarse en distintas formas: 1) nutrimentos intercambiables (adsorbidos a las partículas coloides) y 2) nutrimentos en cantidades totales (en este caso, las intercambiables y no intercambiables).

Los valores obtenidos para N, P, y K se detallan en la Figura 5. No se encontró diferencias significativas asociadas a los factores Ciclo (C), Implemento (I), o su interacción (C x I) para los contenidos de N, mientras que los elementos P y K presentaron diferencias asociadas solamente al factor C (Tabla 2; p-value <0,05). Estas diferencias se pudieron asociar con la variabilidad en la temperatura, la evapotranspiración y la precipitación pluvial en ambos ciclos de cultivo, donde observó un incremento en la magnitud de las lluvias durante los meses de mayo, junio y septiembre del segundo ciclo de cultivo (Figura 6), lo cual pudo haber tenido un efecto en el crecimiento de la cebolla y la dinámica de los elementos en el suelo, lo que de forma indirecta podría afectar la cantidad de sedimento adherido al bulbo.

En un ciclo de cultivo de cebolla, se pueden perder hasta 3,04 ± 0,27 kg de N, 0,15 ± 0,02 kg de P, y 1,16 ± 0,12 kg de K por hectárea (Figura 5). Un aspecto importante para considerar es la naturaleza de dichas cantidades, donde el N corresponde a una cantidad total, mientras que el P y K corresponden solamente a la fracción intercambiable (extraída con KCl-Olsen; ^{Bertsch et al. 2005}), por lo que las pérdidas de P y K como SLCH podrían ser significativamente mayores.

Figura 5 Pérdida de macronutrientes (N, P, K) por la cosecha del cultivo de cebolla como respuesta al efecto de tres implementos de labranza en ciclos de producción no continuos, en un Andisol de Tierra Blanca, Cartago. Barras de error como error estándar. * Sin diferencias significativas (p-value >0,05) entre las medias por los factores evaluados (Tabla 2). Letras distintas denotan diferencias significativas (p-value <0,05) entre las medias para el factor Ciclo (C) evaluado en el modelo (Tabla 2). Los valores de P y K corresponden a cantidades intercambiables, mientras que el N corresponde a N total. P analizado como PO₄ ^-. 

Figura 6 Dinámica de las variables climáticas determinadas para el sitio de estudio en Tierra Blanca de Cartago, durante los dos ciclos del cultivo de la cebolla (2019-2020).  

Al respecto, ^{Fallas-Corrales (2013}) encontró contenidos totales de P en cinco Andisoles de Costa Rica entre 760 - 2009 mg kg^-1 P, lo que podría representar una pérdida total entre 1,48 - 3,93 kg ha^-1 ciclo^-1 de P (10 - 26 veces más que los obtenidos en el presente estudio) si se considera en el cálculo los valores máximos de sedimentos obtenidos en las parcelas de cebolla evaluadas en ambos ciclos. De forma similar al P, las cantidades de K que se podrían perder como sedimentos durante la cosecha de cebolla podrían ser hasta 3 - 8 veces mayores (3,07 - 9,58 kg ha^-1 ciclo^-1 de K, si consideramos los contenidos totales de K en los Andisoles, como los reportados por ^{Molina et al. (1986}) para seis Andepts de Costa Rica (actualmente Andisoles) los cuales variaron entre 1575 - 4250 mg k kg^-1.

Los resultados obtenidos para el Ca y el Mg se detallan en la Figura 7. No se encontraron diferencias significativas asociadas al implemento de labranza, el ciclo de cultivo o su interacción para el Ca, mientras que en el caso del Mg se observó diferencias significativas asociadas al factor ciclo de cultivo (Tabla 2; p-value <0,05). Se encontró valores máximos de pérdida de nutrimentos por cosecha para estos dos elementos de hasta 2,28 ± 0,07 kg de Ca y 0,33 ± 0,02 kg de Mg por hectárea por ciclo.

Figura 7 Pérdida de cationes (Ca y Mg) por la cosecha del cultivo de cebolla como respuesta al efecto de tres implementos de labranza en ciclos de producción no continuos, en un Andisol de Tierra Blanca, Cartago. Barras de error como error estándar. * Sin diferencias significativas (p-value>0,05) entre las medias por los factores evaluados (Tabla 2). Letras distintas denotan diferencias significativas (p-value<0,05) entre las medias para el factor Ciclo (C) evaluado en el modelo (Tabla 2). Los valores de Ca y Mg corresponden a cantidades intercambiables. 

Finalmente, las cantidades de Fe, Cu, Zn y Mn que se pierden por labores de cosecha se resumen en la Figura 8. No se encontraron diferencias significativas asociadas al implemento de labranza, el ciclo de cultivo o su interacción para los cuatro microelementos (Tabla 3; p-value >0,05). Se encontró valores máximos de pérdida de hasta 0,23 ± 0,10 kg de Fe, 0,03 ± 0,01 kg de Cu, 0,02 ± 0,01 kg de Zn y 0,07 ± 0,04 kg de Mn por hectárea por ciclo de cultivo de la cebolla. Dichos valores resultaron significativamente más bajos comparación con los obtenidos para el N, K, y Ca (Figuras 5 y 7). Sin embargo, estas cantidades representan la fracción intercambiable de los elementos en el suelo (^{Bertsch et al. 2005}), por lo que es esperar que la pérdida total sea aún mayor.

Figura 8 Pérdida de micronutrientes (Fe, Cu, Zn, y Mn) por la cosecha del cultivo de cebolla como respuesta al efecto de tres implementos de labranza en ciclos de producción no continuos, en un Andisol de Tierra Blanca, Cartago. Barras de error como error estándar. *Sin diferencias significativas (p-value>0,05) entre las medias por los factores evaluados en el modelo (Tabla 2). 

Pérdida de materia orgánica por labores de cosecha del cultivo de cebolla

Comúnmente, la pérdida de suelo por procesos erosivos (naturales y antrópicos) se asocia con la pérdida de la capa más fértil que es capaz de sustentar la producción agrícola. Sin embargo, un componente normalmente olvidado en los estudios de erosión y pérdida de suelo es el contenido de materia orgánica (MO) que se pierde por estos fenómenos, promovido en gran medida por la conversión de ecosistemas a agroecosistemas (^{Lal 2019}).

En el presente estudio, se cuantificó las pérdidas totales de la materia orgánica que salen del sembradío por la cosecha. Las cantidades de materia orgánica obtenidos se detallan en la Figura 9. No se encontró diferencias significativas asociadas a los dos factores evaluados y su interacción (Tabla 4; p-value > 0,05).

Figura 9 Pérdida de materia orgánica del suelo por la cosecha del cultivo de cebolla en respuesta al efecto de tres implementos de labranza durante ciclos de producción no continuos, en un Andisol de Tierra Blanca, Cartago. Barras de error como error estándar. *Sin diferencias significativas (p-value >0,05) entre las medias por los factores evaluados en el modelo (Tabla 2). 

En un ciclo de cebolla se puede perder valores de MO entre 20,43 ± 16,61 y 39,06 ± 4,21 kg de materia orgánica por hectárea (Figura 9), donde dichos valores representan respectivamente 0,02 y 0,04% de la MO total cuantificada en los primeros 20 cm de profundidad de este suelo. Pese a que estas cantidades se pudiesen considerar como muy bajas, dicha pérdida repercute no solamente en la pérdida de materia total del suelo, sino también en el componente fertilidad de suelo (^{Bertsch 1995}), así como en propiedades físicas como el desarrollo y estabilidad de los agregados estructurales del suelo, su porosidad y subsecuentes fenómenos de retención y conducción de agua, así como la aireación del suelo (^{Hartge et al. 2016}).

Muchos de los estudios sobre pérdida de suelo por cosecha del cultivo (SLCH) se enfocan en la pérdida de masa como sedimentos que salen adheridos al cultivo (^{Ruysschaert et al. 2007}, ^{Mwango et al. 2015}, ^{Oshunsanya et al. 2018}, ^{Antoneli et al. 2021}, ^{Thomaz y Bereze 2021}, ^{Kuhwald et al. 2022}), sin contemplar la posible salida no solamente de nutrimentos sino de materia orgánica de los sitios de siembra y sus implicaciones en la calidad y salud del suelo, por lo cual no se tienen valores de comparación con los obtenidos en el presente estudio.

Sin embargo, las pérdidas de carbono orgánico del suelo (SOC) por procesos erosivos son considerables, alcanzando valores globales de hasta 5,7 Pg año^-1 (5,7 x 10¹⁵ g año^-1) los cuales contribuyen de forma significativa tanto a la degradación del suelo como a las emisiones de CO₂ a la atmósfera (^{Lal 2019}).

Un aspecto importante por considerar es el papel de la materia orgánica en la génesis de los Andisoles, donde la interacción entre la materia orgánica y minerales nano cristalinos como la alofana, imogolita y ferrihidrita produce complejos órgano-minerales resistentes a la degradación y descomposición, lo que resulta en la estabilización de dicha materia orgánica (^{Shoji et al. 1993}, ^{Buol et al. 2011}). Por tanto, la pérdida de estos materiales orgánicos (sujeta a cambios en otros factores formadores como el clima) podría repercutir en los procesos pedogenéticos y la evolución del suelo hacia taxones más alterados como Inceptisoles, Alfisoles, Ultisoles y Oxisoles (Shoji et al. 1993), reduciendo gradualmente la capacidad potencial de producción hortícola en las tierras altas de Cartago.

Implicaciones del SLCH en el manejo y conservación del recurso suelo en Costa Rica

Los Andisoles son el tercer grupo taxonómico más abundante a nivel nacional, con una cobertura del 15,9% de la superficie del país (^{Mata et al. 2020}). Debido a sus propiedades física particulares (texturas francas, baja densidad aparente) son suelos muy fáciles de mecanizar, lo cual incrementa la susceptibilidad a la pérdida por erosión hídrica y eólica (^{Buol et al. 2011}). Aparte del poder erosivo de la lluvia (^{Vahrson 1990}), la erosión por labranza es un fenómeno de suma importancia en suelos derivados de cenizas volcánicas, con un potencial de remoción del suelo considerable (hasta 192 kg m^-1) en sitios cuya pendiente se encuentre entre 5 - 35% (^{Mehuys et al. 2009}).

En el presente trabajo se determinaron cantidades significativas de suelo removido del sembradío por cosecha de cebolla (entre 0,59 ± 0,47 y 1,08 ± 0,078 Ton ha^-1 ciclo^-1), las cuales podrían representar un aporte importante a las actuales tasas de erosión reportadas para Andisoles de la zona alta de Cartago, las cuales varían entre 0,7 y 14 Ton ha^-1 ciclo^-1 (^{Castillo 2012}) considerándose como erosión ligera a moderada según la clasificación propuesta por ^{FAO (1979)}. No obstante, en esta región se han obtenido tasas elevadas de erosión pluvial, superiores a 50 Ton ha^-1 año^-1 (^{Lianes et al. 2009}), y valoraciones de erosión moderadas a severas (^{Ramírez et al. 2008}). Sin embargo, dichas tasas no incluyen ni el aporte de la erosión por labranza mecanizada (^{Mehuys et al. 2009}, ^{Tiessen et al. 2010}) o el sedimento exportado por la cosecha del cultivo de la cebolla (^{Antoneli et al. 2021}).

Los resultados de SLCH obtenidos en el presente estudio son pioneros no solamente en Costa Rica, sino también a un nivel regional más amplio (América Central y el Caribe), en los cuales se cuantifica no solamente la pérdida de sedimento adherido al bulbo de la cebolla, sino también la pérdida de materia orgánica y elementos esenciales que se pierden de forma irreversible del sembradío; una información novedosa y útil en comparación con los estudios tradicionales de SLCH (^{Oshunsanya et al. 2018}, ^{Antoneli et al. 2021}, ^{Thomaz y Bereze 2021}, ^{Kuhwald et al. 2022}).

Estos valores tienen serias implicaciones en el potencial agrícola de los Andisoles, ya que la pérdida de sedimentos por procesos erosivos no solamente genera una pérdida de la masa de la capa arable, sino también la salida irreversible de enormes cantidades de nutrimentos tanto adicionados como fertilizante como los contenidos naturales en el suelo, visibilizando la problemática de los procesos de erosión en la degradación del suelo y sus implicaciones en la agricultura y el medio ambiente (^{Lal 2001}). Esto anterior resulta de suma importancia para el manejo de los Andisoles, ya que estos suelos que cubre una gran cantidad del paisaje agrícola del territorio Centroamericano (^{Buol et al. 2011}, ^{Alvarado 2022}).

En adición a las prácticas recomendadas para el manejo y conservación del suelo en la zona alta de Cartago (^{Ramírez et al. 2008}, ^{Lianes et al. 2009}, ^{Mehuys et al. 2009}), se debe explorar soluciones para evitar la salida del sedimento de la finca por la cosecha del cultivo. Por ejemplo, evaluar la factibilidad del lavado del bulbo in situ previo al secado de la cebolla, impidiendo que los sedimentos salgan del sitio de siembra.

Esto anterior representa una gran ventada de oportunidad para el desarrollo de estudios sobre erosión por SLCH, en adición a estudios más complejos que involucren el aporte de la erosividad de la lluvia (^{Vahrson 1990}) o la pérdida de suelo por labranza mecanizada (^{Mehuys et al. 2009}, ^{Marín-Araya 2022}), los cuales no solamente visibilicen esta problemática sobre la degradación de los suelos en la zona alta de Cartago, sino también desarrollen soluciones prácticas para la mitigación de esta problemática, incluyendo a otros cultivos de importancia hortícola en la región como la papa (Solanum tuberosum L.) y la zanahoria (Daucus carota subsp. Sativus Hoffm.).

En adición, se debería considerar otras posibles estrategias de manejo y conservación de suelos, como la rotación de cultivos con especies de la familia Poaceae. Estos sistemas de rotación han demostrado tener efectos positivos en las propiedades físicas del suelo, mejorando la estructuración e incrementando el contenido de carbono y otros nutrimentos esenciales (^{Zani et al. 2020}), así como la mitigación de la erosión de suelos en sistemas agrícolas (^{Boincean y Dent 2019}).

Conclusiones

Bajo las condiciones agroclimáticas de las tierras altas de Cartago, se pueden perder cantidades de sedimentos (adheridos al bulbo) de entre 0,6 y 1,1 Ton ha^-1 en un ciclo de cultivo de cebolla, de forma indiferente del implemento de labranza que se utilice para preparar el terreno previo a la siembra del cultivo. Dichas cantidades podrían incrementar si se realiza más de dos ciclos de producción de cebolla por año.

En esos sedimentos que salen por cosecha del cultivo (SLCH), se pueden perder macronutrientes, microelementos, y materia orgánica (MO) los cuales salen del sitio de siembra de forma irreversible. En el caso del N y el K, pueden salir respectivamente cantidades de hasta 3 y 1,2 kg ha^-1 ciclo^-1, y la salida de MO puede llegar a valores de 39 kg ha^-1 ciclo^-1.

Un aspecto importante para la interpretación de estas cantidades perdidas es la naturaleza del análisis químico de laboratorio que se realizó al sedimento, ya que si se solamente se analiza cantidades intercambiables (con soluciones extractoras KCl-Olsen o Mehlich III), las pérdidas de nutrimentos por cosecha del cultivo pueden ser considerablemente más grandes (hasta 10 veces mayores).

Se recomienda continuar con este tipo de estudios sobre SLCH no solamente en el cultivo de cebolla, sino también en otros cultivos hortícolas como la papa y la zanahoria, explorando otras metodologías de análisis químico que permitan cuantificar las cantidades totales de elementos en los sedimentos colectados.