Introducción
En el Soconusco, Chiapas, por muchos años la preparación del suelo para la siembra de soya ha sido intensa y en condiciones de suelo húmedo, dicha actividad ha comprendido entre uno y siete pasos de rastra antes de los años 90’s. Actualmente, estos suelos presentan niveles variados de compactación que, como lo reporta la literatura en este ámbito (Van Ouwerkerk y Soane, 1994; Cabria y Culot, 1999; Alakukku et al., 2003; Ahmad et al., 2007), el paso repetitivo de la maquinaria compacta el suelo, reduce la permeabilidad y el almacenamiento de la humedad el perfil radical del cultivo (Comia et al., 1994; Van Wie et al., 2013).
También, cuando se mulle el suelo con la rastra, continúa la siembra de la soya y sí inmediatamente llueve, se forma una ‘costra’ que causa taponamiento de la superficie e impide la germinación de la semilla. Cuando esto ocurre se realiza una segunda siembra y genera más gastos de producción. En el ámbito mundial, la compactación de los suelos se interpreta como una forma de su degradación física ocasionada por el uso intensivo de la maquinaria agrícola (Van Ouwerkerk y Soane, 1994; Kuiperd y van de Zande, 1994), este tráfico produce una capa dura (piso de arado) por debajo de la capa arable que reduce la infiltración, obstaculiza el movimiento interno del agua, inhibe el crecimiento radical y el rendimiento de los cultivos (Czarnes et al., 2000; Hagedorn y Bundt, 2002).
Una manera de incrementar la infiltración es a través del subsoleo, cuya ruptura del ‘piso de arado’, incrementa la infiltración, el movimiento interno del agua, su redistribución, su capacidad de almacenamiento y crecimiento radical del cultivo para explorar mayor volumen de la rizosfera (Lampurlanes et al., 2001; Sharma et al., 2004). También, la ruptura del ‘piso de arado’ conduce a la mejora estructural del suelo, ya que permite el intercambio libre de flujos hacia el subsuelo, incrementa el intercambio de gases y reduce los escurrimientos superficiales (Bertolino et al., 2010).
En general, se concibe que la sustentabilidad en la agricultura sólo es posible mediante la conservación del suelo y agua. En condiciones de temporal, el cuidado del agua redunda en la mejora sustentable de la producción siempre y cuando su manejo sea eficiente. De aquí la importancia de conocer el contenido y dinámica del agua en el suelo con fines de su almacenamiento y conservación para los cultivos, principalmente para la soya en condiciones de temporal. La dinámica del agua en el suelo es interferida por varios factores que limitan su aprovechamiento; entre estos, los sistemas de labranza en condiciones de suelo húmedo y por tiempos prolongados (Alakukku et al., 2003).
Durante la compactación, algunas características físicas del suelo cómo reducción de la porosidad, incrementa la densidad aparente y limita el flujo de agua durante su redistribución. Algunas investigaciones sobre las variables hidráulicas del suelo recomiendan el desarrollo de teorías, mediciones en campo y el uso de modelos fisicomatemáticos para comprender el fenómeno de la compresión del suelo sobre sus propiedades hidráulicas (Horton et al., 1994).
Es así qué, la ruptura de la compactación mejora las propiedades hidráulicas del suelo, y por esto son necesarios el uso de métodos fisicomatemáticos para la medición y ajuste de esta información que evidencie la mejora sustantiva de esta práctica. Toda vez qué, la infiltración es la única fuente de ingreso y recarga de agua en el suelo, su conocimiento conlleva a recomendar el subsoleo como medida de conservación del suelo y agua, en particular para el cultivo de soya en condiciones de temporal (Klute y Dirksen, 1986; Cabria y Culot, 2000).
El laboreo repetitivo para siembra de soya en el Soconusco, Chiapas; reveló que estos suelos mantienen de mediana a alta compactación observable mediante perfiles de suelo, lo cual originó esta investigación cuyo objetivo fue evaluar por tres años consecutivos la infiltración en tres sistemas de manejo del suelo y su conductividad hidráulica a saturación de un predio, en cuyo suelo año tras año se siembra la soya en las condiciones mencionadas, así como la mejora del rendimiento del cultivo.
Materiales y métodos
Localización del área de estudio
El experimento se ubicó en el municipio de Tapachula, Chiapas; entre las coordenadas 14° 45’ latitud norte y 92° 23’ longitud oeste a una altitud de 16 m. El clima es cálido subhúmedo, cuya temperatura promedio fue 28 ±1 °C y la precipitación acumulada promedio anual es de 1 110 mm. El suelo con base en su textura es un tipo franco, el pH fue de 6.6 ligeramente ácido y un contenido de materia orgánica de 2.5%. El periodo experimental comprendió tres años consecutivos y durante los ciclos de producción primavera-verano (PV) 2017, 2018 y 2019.
Tratamientos y diseño experimental
El experimento consistió en tres tratamientos de 1.5 ha (150 m de largo por 100 m de ancho) distribuidos como parcelas independientes. El subsoleo (SUB) más un paso de rastra y siembra mecanizada, barbecho (BAR) más un paso de rastra y siembra mecanizada, y rastreo (RAS), dos pasos de rastra y siembra mecanizada. Anualmente la preparación del suelo para la siembra se hizo en la segunda quincena de junio. Las fechas de siembra fueron variables, 10 de julio en 2017, 12 y 8 de julio en 2018 y 2019, respectivamente.
El material vegetal fue la variedad de soya Huasteca 100 generada por el Campo Experimental Sur de Tamaulipas dependiente del Instituto Nacional de Investigaciones Forestales, Agrícolas y Pecuarias (INIFAP). Esta variedad es un material de porte y ciclo intermedio. En relación con la infiltración del agua en el suelo, ésta se sustenta en toda una teoría descrita mediante modelos fisicomatemáticos que gobiernan su comportamiento.
Estudios preliminares
A través de perfiles de suelo abiertos en enero de 2017, se detectó el ‘piso de arado’ aproximadamente a ±35 cm de profundidad, lo que indicó una profundidad para la práctica del subsoleo a 70 cm durante el periodo seco (abril). Dichos perfiles fueron observados al interior de cada tratamiento (tres) para este propósito, de donde al inicio del ciclo experimental se extrajeron muestras de suelo a las profundidades de 0-30 cm y se determinaron sus características físicas.
A partir de éstas se estimaron los parámetros de humedad y la conductividad hidráulica a saturación inicial para su comparación después del periodo experimental (tres años). Dichas estimaciones se hicieron con base en su textura, mediante el software SPAW V. 6.02.75; cuyos parámetros de humedad fueron: contenido volumétrico de humedad a capacidad de campo (ƟCC), punto de marchitamiento permanente (ƟPMP) contenido de humedad a saturación (Ɵs), tasa de infiltración (q0) conductividad hidráulica a saturación (Ks) y su densidad aparente (Da).
La infiltración
Brevemente se describe la teoría de la infiltración del agua que, en este estudio fue observada en los diferentes sistemas de manejo del suelo donde se ha cultivado la soya por más de 40 años consecutivos en condiciones de temporal. En términos generales el fenómeno de la infiltración no es más que la entrada del agua al interior del suelo por su superficie, físicamente en condiciones de suelo seco y de esta manera, este fenómeno está gobernado por los potenciales o fuerzas mátrica y gravitacional. Sin embargo, a medida que el agua entra al suelo, este se humedece y el potencial mátrico desaparece debido a que los poros se llenan de agua (suelo saturado) y cesan las fuerzas capilares; entonces, cuando esto ocurre sólo participa el potencial gravitacional (Cabria y Culot, 1999).
La descripción de este fenómeno a través del tiempo ha sido analizada por algunos autores mediante expresiones matemáticas cuya aproximación cuasi analítica revelan la dinámica de la infiltración (Green y Ampt, 1911; Kostiakov, 1932; Swartzendruber, 1993; Horton, 1940; Holtan, 1961; Philip, 1957), pero actualmente se continúa en la interpretación fisicomatemática de estos principios para determinar otras funciones que mejor describan el fenómeno.
Entre los autores citados, Philip (1957) dedujo una expresión para la infiltración acumulada, I(t), que hasta ahora presenta buenos resultados, ya que a partir de esta ecuación truncada en el término cuadrático se deducen los coeficientes mediante los cuales se estima la conductividad hidráulica a saturación (Ks) por Kutilek y Krejca (1987). En relación con los antecedentes señalados, es evidente que la fuente de ingreso del agua en el suelo es la infiltración, ya que mediante ésta se determina la cantidad de agua que ingresa al suelo por su superficie en función del tiempo.
En este trabajo de investigación se midieron la tasa de infiltración, q0(t) y la lámina de infiltración acumulada, I(t). En relación con el flujo de agua (q0), el subíndice 0, se refiere a la entrada de agua por la superficie del suelo, Z= 0. Existen métodos para medir la infiltración como el de doble cilindro (cilindros concéntricos), cuya información se ajusta por modelos fisicomatemáticos que permiten analizar su comportamiento.
De este modo, la tasa de infiltración, q0= dI/dt, (ecuación 1), es una magnitud (cm min-1), cuyo comportamiento es una función monótona potencial que disminuye rápidamente a través del tiempo (q0) y se aproxima a un valor constante; es decir, t= 0; cuando q0 → ∞ (ecuación 2); y cuando t → ∞, q0= constante; es decir, en teoría; q0= Ks (conductividad hidráulica a saturación; cm min-1); cuando t → ∞ (ecuación 3
Sucesivamente la infiltración del agua en el suelo es una variable del volumen de agua que ingresa en el suelo por unidad de área a través de un tiempo determinado, y se denota cómo infiltración acumulada, I(t), misma que se expresa en unidades lineales como lámina acumulada en el suelo en función del tiempo. Ésta es una curva continua ascendente expresada por la ecuación 4.
1).
2),
3)
4)
La tasa de infiltración, q0(t) = dI/dt, como se señaló, disminuye exponencialmente a través del tiempo y se aproxima a un valor constante después de un tiempo largo de infiltración. Entonces, si t → ∞, q0 → Ks, de tal manera que la tasa de infiltración alcanza un valor constante y se aproxima a la conductividad hidráulica a saturación, Ks (Chow et al., 1988). En general, las mediciones de la infiltración son ampliamente utilizadas para estimar la conductividad hidráulica a saturación (Ks), ya que esta propiedad física del suelo es clave para describir sus propiedades hidráulicas y el balance hidrológico sobre la superficie terrestre (Campbell, 1985; Van Looy et al., 2002; Hillel, 2003; Lal y Shukla, 2004; Morbidelli et al., 2011).
Además, su valor es un parámetro fundamental para el diseño correcto de sistemas de riego. Algunas propiedades del suelo, tales como su textura, densidad aparente, contenido de humedad inicial, capas o estratos en el perfil del suelo, pendiente del terreno, cobertura del suelo (residuos de cosecha y pedregosidad), patrón de la lluvia y el tiempo influyen en la variabilidad de la infiltración (Smith et al., 2002). En la tercera etapa de esta investigación se realizaron dos mediciones de la infiltración al interior de cada tratamiento.
La tasa de infiltración, q0(t) = dI/dt, cuyo comportamiento es similar a la ecuación de Kostiakov (1932), muestra un gráfico para un tiempo largo de la infiltración y experimenta una magnitud que decrece exponencialmente, y cuando este valor se aproxima a la asíntota con la horizontal, dicho valor es la infiltración básica (IB) del suelo particular. Asimismo, la infiltración acumulada, I(t) es la suma de la infiltración en el tiempo y representada por la ecuación de Philip (1957), reducida a los primeros tres términos (ecuación 5), misma información que fue sometida a su ajuste mediante la ecuación 6 (Kutílek y Krejča, 1987). I(t)= C1t1/2 + C2t + C3t3/2 + C4t2 + Cmtm/2 5). Donde: C1; C2, C3, Cm; son parámetros de la ecuación; y t es el tiempo.
Kutilek y Krejca (1987) sugirió los primeros tres términos de la ecuación de Philip (1957), de cuyos coeficientes de ajuste se estima la conductividad hidráulica a saturación (Ks) mediante la ecuación 6. Ks = (3C1*C3)1/2 + C2 6). Donde: C1 estima la sortividad, C3, y C2 son parámetros de la ecuación. A partir de la información de campo, la infiltración y otras características hidráulicas del suelo objeto de esta investigación se realizaron las estimaciones y ajustes correspondientes por los modelos señalados y procesados mediante el Software Curve Expert V. 2.6.
Resultados y discusión
En el Cuadro 1 se presentan la textura del suelo y los parámetros de humedad en relación con los diferentes tratamientos de manejo del suelo. Se observó que para el mismo suelo con base en su textura; después del subsoleo (abril de 2017), estos parámetros y su densidad aparente sólo mostraron ligera variación, denotados como valores iniciales. Sobre estos resultados, otras investigaciones reportaron similar variabilidad durante el primer año de observación y señalan que esto se debe a la misma heterogeneidad del suelo (Kutilek y Nielsen, 1994).
Textura del suelo | Variables hidráulicas | Subsoleo | Barbecho | Rastreo | |||||
---|---|---|---|---|---|---|---|---|---|
2017 | 2019 | 2017 | 2019 | 2017 | 2019 | ||||
Franco | ƟPMP (%) | 15.4 | 15.4 | 15.1 | 15.1 | 15.8 | 15.8 | ||
ƟCC (%) | 28.2 | 28.2 | 27.7 | 27.7 | 28.9 | 28.9 | |||
Ɵs (%) | 40.8 | 40.8 | 39.5 | 39.5 | 39.4 | 39.4 | |||
Da (g cm-3) | 1.04 | 1.04 | 0.97 | 0.97 | 0.94 | 0.94 | |||
q0 (cm min-1) | 0.027 | 0.06 | 0.011 | 0.044 | 0.019 | 0.033 | |||
Ks (cm h-1) | 1.27 | 2.34 | 1.32 | 2.1 | 1.57 | 2.22 | |||
Δ de Ks (%) | 184 | 159 | 141 |
ƟPMP= contenido volumétrico de humedad a punto de marchitamiento permanente en porciento; Ɵcc= contenido volumétrico de humedad a capacidad de campo en porciento; ƟS= contenido volumétrico de humedad a saturación en porciento; Da= densidad aparente del suelo; q0= tasa de infiltración; Δ= incremento y Ks= conductividad hidráulica a saturación.
En el Cuadro 1, también se muestra que el subsoleo en relación con los demás sistemas de manejo del suelo (barbecho y pasos de rastra) después de tres años, éste reportó un cambio significativo en la tasa de infiltración, cuyo valor se incrementó 184% con respecto a su valor inicial. Asimismo, el tratamiento de barbecho, su tasa de infiltración experimento un incremento con respecto a su valor inicial de 159% y el tratamiento de pasos de rastra se incrementó 141% durante este periodo experimental.
Aun cuando los resultados corresponden a un periodo de mediano plazo, se infiere qué entre los sistemas de manejo del suelo evaluados, el subsoleo evidencia que su práctica promueve el incremento de la infiltración y otras propiedades físicas e hidráulicas del suelo. Esta afirmación puede observase en la Figura 1, donde se muestra la dinámica de la tasa de infiltración en el tiempo. En tanto que en la Figura 2 se muestra la infiltración acumulada en el tiempo observado, donde el valor superior de la lámina infiltrada corresponde al subsoleo.
Con respecto a la lámina infiltrada acumulada, I(t), en la Figura 2 se muestra evidencia del impacto del sistema de subsoleo en los suelos productores de soya en condiciones de temporal, toda vez que en el corto plazo bajo el sistema de manejo del suelo evaluados, la lámina infiltrada acumulada, I(t), en el perfil radical del cultivo fue de 38 cm, seguida por los sistemas de manejo del suelo de sólo pasos de rastra y barbecho del suelo, cuyas láminas infiltradas acumuladas fueron de 28 cm y 21 cm, respectivamente, y para tiempos equitativos de evaluación de ocho y media horas en todos los sistemas de manejo del suelo evaluados.
Con base en estos resultados, las láminas de infiltración acumuladas revelaron ventaja potencial del sistema de subsoleo para el manejo del suelo en la producción sustentable de la soya en condiciones de temporal en Chiapas. Asimismo, otras ventajas además de las propiedades hidráulicas señaladas contribuyen a mejorar las propiedades físicas en suelos compactados. Algunas investigaciones en este ámbito señalan que el subsoleo con cincel, como fue el caso en esta investigación, produjo incremento en las tasas de infiltración (Álvarez et al., 2006; Soracco, 2009).
En el marco de investigaciones similares sobre el subsoleo, los análisis estadísticos no reportaron diferencia significativa sobre las propiedades físicas del suelo, en particular sobre la densidad aparente; mientras que, la tasa de infiltración fue altamente significativa sobre los tratamientos sin subsoleo (Solhjou y Niazi Ardekani, 2001).
En el mismo sentido, Desale et al. (2012) reportó altas tasas de infiltración, de infiltración acumulada y mayor retención de humedad mediante el subsoleo comparado con la labranza tradicional. También, se remarcan efectos atribuibles al subsoleo del cual se reportó incremento del contenido de humedad en la zona radical de los cultivos por la mejora asociada a la macroporosidad e inducción de mayores tasas de infiltración e incremento de la redistribución radical en el perfil del suelo (Mohanty et al., 2007; Zibilske and Bradfort, 2007).
Otras investigaciones remarcan ventajas sustantivas del subsoleo para mejora de la infiltración y redistribución del agua al interior del perfil del suelo (Sushil et al., 2018), derivado de la infiltración, su conductividad hidráulica a saturación ésta contribuye a la recarga de los acuíferos. En general el efecto del subsoleo sobre los suelos afectados por la labranza repetitiva y por periodos prolongados de tiempo como ocurre en el Soconusco, Chiapas, la presente investigación reveló resultados preliminares aceptables para la mejora de la infiltración y la conductividad hidráulica a saturación.
En relación con la conductividad hidráulica a saturación, en el Cuadro 1, hilera 8 se observó que el porcentaje superior de esta variable en el mediano plazo correspondió al subsoleo sobre el resto de los sistemas de manejo del suelo cultivados consecutivamente con soya en condiciones de temporal. Asimismo, otras investigaciones señalan que el movimiento interno del agua en el suelo no puede detectarse de un ciclo a otro y dicho fenómeno solo es observable a través de un perfil de suelo después de un periodo largo de tiempo (Klute y Dirksen, 1986; Cabria y Culot, 1999).
Sin embargo, esta investigación permitió observar que, en el corto plazo (tres años) el subsoleo incrementó la conductividad hidráulica a saturación. En lo sucesivo, investigaciones similares indican que el subsoleo en el cultivo de soya experimentó ventaja superior sobre la labranza convencional, toda vez que mediante esta práctica se incrementa la tasa de infiltración y el almacenamiento de agua en el perfil radical del cultivo, ya que dicha condición superior de humedad, minimiza el riesgo de déficit hídrico y favorece la productividad del sistema de producción soya en condiciones de temporal (Heatherly and Spurlock, 2001; Wesley et al., 2001).
Rendimiento y componentes de rendimiento en función del subsoleo
El análisis de resultados en relación con los sistemas de manejo del suelo en esta investigación, principalmente el subsoleo, éste en primer lugar, mostró clara evidencia de su efecto, ya que en este sistema se observó mayor capacidad de almacenamiento de agua en la zona radical del cultivo de soya en condiciones de temporal, en función de una lámina acumulada superior de 38 cm, comparada con los sistemas de barbecho y rastreo del suelo, respectivamente.
Por otra parte, se analizó el impacto de los sistemas de manejo del suelo sobre el rendimiento de soya y de sus variables agronómicas de interés. Estos resultados se presentan en el Cuadro 2, cuyos análisis de varianza de los años analizados (2017, 2018 y 2019) y su posterior separación de medias (p≤ 0.05), revelaron que el rendimiento, número de vainas, ramas y entrenudos; el subsoleo marcó una tendencia superior sobre los demás tratamientos. Sin embargo, en relación con las variables altura de planta y peso de 100 granos, el subsoleo mantuvo un comportamiento similar a los demás tratamientos (barbecho y rastreo).
Tratamientos | Rendimiento (kg ha-1) | Altura-planta (cm) | Núm. vainas (unidades) | Peso de 100 granos (g) | Núm. ramas (unidades) | Entrenudos (unidades) |
---|---|---|---|---|---|---|
Subsoleo | 2 634 a | 76.6 a | 85.3 a | 17.5 a | 5.8 a | 13.7 a |
Barbecho | 1 873 b | 71.1 a | 50.8 b | 18.9 a | 3.4 b | 12.6 a b |
Rastreo | 1 772 b | 68.7 a | 35.9 b | 18.7 a | 5.2 a | 10.1 b |
Medias de variables de respuesta. Tukey (p≤ 0.05)
En resumen, el subsoleo en esta investigación confirmó una mejora sustantiva para el manejo del suelo en el cultivo de soya en condiciones de temporal, ya que es un sistema de manejo del suelo efectivo para romper la capa de arado, reduce la densidad aparente y promueve la expansión de raíces en el subsuelo. Éste también ha demostrado mejorar las operaciones de labranza para la mejora de la rizófora, la diversidad microbiana y la capacidad de almacenamiento de agua en el suelo (Desale et al., 2012). Asimismo, el subsoleo establece condiciones ideales para el desarrollo radical y retarda la senescencia de la planta; aumenta el rendimiento de grano y hace posible el incremento de la densidad de población en este cultivo (Mohanty et al., 2007; Zibilske y Bradfort, 2007).
Conclusiones
Se determinó que la ruptura del piso de arado a través del subsoleo en suelos productores de soya de El Soconusco, Chiapas, en el corto plazo, por una parte, incrementó la tasa de infiltración y su conductividad hidráulica a saturación y por otra parte, incrementó el rendimiento y sus componentes de rendimiento.