Agustín Navarro Bravo^1§, Benjamín Figueroa Sandoval², Dora Ma. Sangerman-Jarquín³ y Esteban S. Osuna Ceja³

¹ Centro Nacional de Estandarización de Maquinaria Agrícola (CENEMA). INIFAP. Carretera Los Reyes- Lechería, km 18.5. Texcoco, Estado de México. C. P. 56230. ^§ Autor para correspondencia: navarro468@yahoo.com.mx.

² Campus San Luis Potosí. Colegio de Postgraduados en Ciencias Agrícolas. Iturbide Núm. 73, Salinas de Hidalgo, San Luis Potosí. C. P. 78600. México. (benjamin@colpos.mx).

³ Campo Experimental Valle de México. INIFAP. Carretera Los Reyes-Texcoco, km 13.5. Coatlinchán, Texcoco, Estado de México, C. P. 56250. (sangerman.dora@inifap.gob.mx).

³ Campo Experimental Pabellón. INIFAP. km 32.5. Carretera Aguascalientes-Zacatecas, Mpio. Pabellón de Arteaga, Aguascalientes. C. P. 20660. Tel: 01 465 958 01 86. (osuna.salvador@inifap.gob.mx).

El manejo físico del suelo implica optimizar, en la zona de las raíces de las plantas, las propiedades físicas, tales como: estructura, porosidad, infiltración, retención de agua, entre otras; que son importantes para el manejo sustentable del recurso e identificación de indicadores físicos del suelo. Por ello, se evaluó un agroecosistema de haba-maíz-frijol en rotación con labranza de conservación, para conocer los indicadores físicos a través de las propiedades físicas y químicas del suelo que se relacionen con el rendimiento de los cultivos, en el Colegio de Postgraduados en Ciencias Agrícolas, Texcoco, Estado de México. Las propiedades de suelo medidas fueron resistencia mecánica, densidad aparente, materia orgánica, conductividad eléctrica, porosidad, distribución del tamaño de partículas, humedad aprovechable, saturada y residual, conductividad hidráulica saturada, infiltración acumulada, sorbilidad; asimismo, se evaluó el rendimiento de haba, maíz y frijol. Los datos se analizaron mediante análisis de regresión simple y múltiple. Los resultados mostraron que la materia orgánica y la conductividad hidráulica saturada, son los indicadores físicos del suelo que definen el rendimiento de los tres cultivos bajo labranza de conservación.

Palabras clave: estructura del suelo, rendimiento de grano, sostenibilidad.

Introducción

El suelo es un recurso natural y un sistema muy complejo; que permite el sostenimiento de las actividades productivas, pero por su amplia variabilidad, es difícil establecer una sola medida física, química o biológica que refleje su calidad (Bandick y Dick, 1999).

Los métodos agrícolas de producción de cultivos se han concentrado en la labranza intensiva, suministro de insumos externos como estrategia para incrementar la fertilidad del suelo y el rendimiento de la cosecha. Sin embargo, las consecuencias de la labranza intensiva son la compactación del suelo, la pérdida de la biodiversidad, la contaminación del suelo y el agua por el uso excesivo de agroquímicos, así como el incremento de la erosión. Por el contrario, la labranza de conservación es un esquema de producción de cultivo, diferente a los métodos convencionales y es una alternativa de manejo productivo del suelo. Sus principales ventajas son: mejorar el contenido de materia orgánica en la superficie del suelo, mejorar la estabilidad estructural, incrementar la retención de agua, y reducir el escurrimiento y producción de sedimentos. Lo anterior promueve aumento en la calidad física de los suelos agrícolas y del ambiente en general (Gregorich y Carter, 1997).

Para cuantificar la calidad física del suelo, los indicadores necesitan ser medidos espacialmente a través de la evaluación de las propiedades del suelo (físicas, químicas y biológicas), las cuales deben ser fáciles de medir, aun los más sensibles a cambios generados por las prácticas de manejo, cuyos valores indicarían la calidad (Larson y Pierce, 1991; Dexter, 2004).

La calidad física del suelo se asocia con el uso eficiente del agua, nutrientes, y pesticidas; lo cual reduce las emisiones de gases que generan el efecto invernadero. Como consecuencia, se presenta un incremento en la producción agrícola (Lal et al, 1998). La calidad física del suelo no se puede ser medir directamente, pero se infiere a través de indicadores de calidad del suelo estáticos o dinámicos, y medición de atributos del mismo, que se influencian por el uso y prácticas de manejo (Sanchez-Maranon et al., 2002; Carter, 2002; Dexter, 2004).

Los cambios temporales en las propiedades del suelo con la labranza, son indicadores que se usan para determinar si mejora la calidad del suelo, si es estable o si disminuye con el tiempo (Lal, 1998). Lo anterior permite establecer niveles críticos de indicadores de calidad del suelo, asignar un factor de peso y relacionarlo con la productividad (Lal, 1994).

La hipótesis de este estudio fue que el rendimiento de los cultivos evaluados bajo el sistema de rotación que se implemento en labranza de conservación, se puede explicar mediante el comportamiento de las propiedades físicas y químicas del suelo. Por tanto, el objetivo del presente trabajo fue: encontrar los indicadores físicos del suelo a través de las propiedades físicas y químicas de éste, y que se relacionen con el rendimiento del haba maíz y frijol.

Materiales y métodos

El estudio se realizó en el Campus Montecillo del Colegio de Postgraduados, Texcoco, Estado de México, que se ubica en las coordenadas geográficas (19° 28' 4. 26" latitud norte y 98° 53' 42.18" longitud oeste con una altitud de 2 250 m). La precipitación media anual es de 550 mm, con una buena distribución (> 450 mm) durante la estación de crecimiento del haba, maíz y frijol, y temperatura media anual de 18 ^oC con un período libre de heladas de 150 días. El suelo tiene textura migajón arcillo limoso, de origen aluvial y lacustre, moderadamente drenado, clasificado como Fluvaquentic Endoaquolls (Gutiérrez-Castorena, 1999).

Durante el ciclo primavera-verano de 1995 se estableció un experimento de largo plazo con labranza de conservación, donde se realizan prácticas de rotación de cultivos de gramíneas y leguminosas para producción de grano. En esta área experimental (0.84 ha), se realizó un análisis geoestadístico de materia orgánica (MO) en la capa arable (5, 10, y 15 cm de profundidad) del suelo, porque en esta capa hay mayor acumulación de MO en un sistema de producción de cultivos con labranza de conservación, dando las 10 parcelas de escurrimiento (2 m * 7 m, 14 m²) con aplicación de lluvia simulada con diferentes niveles de MO, para medir propiedades físicas y químicas del suelo en tres puntos centrales equidistantes y lineales de las parcelas de escurrimiento; así mismo se evaluó el rendimiento de grano de los cultivos dentro y alrededor de dichas parcelas.

El contenido de humedad volumétrica (θ_v) en el perfil edáfico (15 cm), la resistencia mecánica (Rm) y la infiltración acumulada del agua (I_a), se determinaron con métodos no invasivos de exploración electromagnética: reflectometría del dominio del tiempo (TDR), penetrómetro de ultrasonido (PU) y un simulador de lluvia móvil tipo ARS/USDA respectivamente. Además se uso métodos tradicionales para medir densidad aparente (ρ_a), materia orgánica (MO), conductividad eléctrica (CE) y textura del suelo.

La determinación de ρ_a se realizó al igual que la θ_v y Rm en cada punto de muestreo a 15 cm de profundidad, se utilizó el método del cilindro de volumen conocido (Jury et al, 1991).

La porosidad total (ƒ) se determinó a partir de la densidad aparente (ρ_a) medida y se consideró una densidad real igual a 2.6 Mg m³ (Hillel, 1980).

La humedad aprovechable (θ_a) se estimó como la diferencia en humedad entre capacidad de campo (CC) y punto de marchitamiento permanente (PMP); la humedad residual (θ_r) fue predicha por el modelo propuesto por Rawls et al. (1991).

La infiltración del agua (I) se determino al usar un simulador de lluvia tipoARS/USDA (Norton et al., 1995), con el cual se aplicaron tres eventos de lluvia en cada área de muestreo respectivo con una intensidad de 20, 60 y 100 mm h^-1. La tasa de infiltración se determinó a partir de la diferencia entre la intensidad de la lluvia y el escurrimiento ocurrido (Q). La conductividad hidráulica saturada (K_s) se determino con los datos de la tasa de infiltración, los cuales se graficaron contra el tiempo para los eventos de 60 y 100 mm h^-1 de lluvia que se aplico en cada área de muestreo, y se calculó el promedio aritmético de los tres últimos valores de la tasa de infiltración. Estos valores corresponden a la porción de la curva de velocidad de infiltración donde ésta se hace constante. La infiltración acumulada (I_a) se obtuvo de la diferencia entre la lluvia total y el escurrimiento. La sortilidad (S_o) se calculó de acuerdo con Chong (1983).

La siembra de los cultivos de haba, maíz y frijol se realizó a equidistancias (entre hileras y entre plantas) de 30 cm en forma manual, para garantizar una densidad de 100 000 plantas ha^-1 en cada uno de los cultivos evaluados; como alternativa para obtener una alta población de plantas en un sistema de labranza de conservación que maneja residuos de cosecha sobre la superficie del suelo.

Se determinaron las estadísticas descriptivas básicas y posteriormente se hicieron análisis de correlación por pares de parámetros con los 14 atributos de suelo medidos en la capa de 0-15 cm de profundidad y el rendimiento de los tres cultivos, se uso el software SAS V8.0 (1999). Los atributos que se uso para el análisis de correlación fueron: Rm, ρ_a, MO, CE, ƒ, contenido de arena(Ar), limo (Li), arcilla (Ac), θ_a, θ_s, θ_r, K_s, I_a, S_o, rendimiento de haba (Ŷ_h), maíz (Ŷ_m) y frijol (Ŷ_f). También se realizó un análisis de regresión múltiple, con el rendimiento de grano de los cultivos (Ŷ_h, Ŷ_m, y Ŷ_f) como variables dependientes y los atributos del suelo como variables independientes.

Resultados y discusión

Las propiedades del suelo que presentaron coeficientes de variación (CV), mayores de 25% fueron la Rm, I_a, CE y K_s. Lo que indicó una mayor variación espacial de estos en el área donde se ubicaron las parcelas. Los coeficiente de variación de las variables físicas: θ_a, θ_s, θ_r, S_o, Ar, Lm y Ac, fueron menores de 20%; pero la ρ_a y ƒ fueron cercanos a cero o con una distribución espacial uniforme. Para el rendimiento de grano del haba (Ŷ_h), maíz (Ŷ_m,) y frijol (Ŷ_f), los coeficientes fueron de bajos a moderados, con un rango de 4.38 a 24.96%. Otras estadísticas tales como: mínimo, máximo, media, mediana y la desviación estándar se muestran en el (Cuadro 1).

El análisis de correlación por pares de parámetros de los 14 atributos de suelo y tres rendimientos de cultivos, mostró una correlación significativa con una probabilidad p≥ 0.05 en 21 de los 136 pares de atributos del suelo (Cuadro 2). La Rm se correlacionó negativamente con CE (R²= -0.85), pero con las demás variables disminuyó; la CE tuvo un R²= 0.75 con la Ar, y disminuyó estadísticamente con las otras variables. Por el contrario, la MO se correlacionó ficativamente con ρ_a, ƒ, K_s, S_o, Ŷ_h y Ŷ_m, (R²> 0.70), y correlación R² = 0.59 con el rendimiento del frijol. Una alta correlación positiva fue obtenida entre I_a con K_s y Ŷ_f· (R²> 0.77), pero negativamente con ρ_a (R²= -0.55). La K_s, tuvo una alta correlación positiva con Ŷ_h, Ŷ_m, Ŷ_f, I_a y MO (R²> 0.72) y una correlación negativa con Ac (R²= -0.7). La S_o se correlacionó positivamente con MO y ƒ (R²> 0.72), y negativamente con ρ_a (R²= -0.72).

La ƒ se correlaciona positivo con S_o y MO (R²= 0.72), pero con ρ_a tuvo una R²= -1, (dado que ésta se estimó a partir de la ρ_a). De la misma manera, existió una alta correlación positiva entre θ_a y θ_s (R²= 1), pero estas dos variables presentaron una correlación baj a (R²< 0.52) con las otras variables estudiadas. La θ_r se correlacionó positivamente con Lm (R²= 0.75) y negativamente con Ac (R²= -0.99). En este análisis de correlación las variables θ_a y θ_s fueron las únicas que no se correlacionaron significativamente con el total de variables estudiadas (Cuadro 2).

Al efectuar las regresiones lineales simples, la MO, K_s, I_a y S_o tuvieron una regresión lineal aceptable con el rendimiento de los tres cultivos. En la Figura 1 se ilustran que la MO y K_s se relacionaron con el rendimiento de haba, maíz y frijol (R²< 0.60; Figura 1a-f); S_o con el rendimiento de haba (R²= 0.42; Figura 1g); mientras que la I_a se relacionó con el rendimiento de maíz y frijol (R²< 0.66; Figura 1i-j).

El rendimiento del maíz se influenció por K_s, MO, I_a, S_o, ρ_a, ƒ, Ac, Rm, θ_r y CE; las cuales explicaron individualmente de 51.8 a 14.2% (Figura 2b). Además, se observó que la línea acumulativa inicia con K_s a S_o, con un ligero aumento en porcentaje de acumulación, pero de ƒ a Rm y de θ_r a CE, existe un mayor aumento que finaliza con un 100% de explicación. Por el contrario, en el tramo de S_o a ƒ y de Rm a θ_r la pendiente es cero, es decir, indica que en forma acumulativa estas variables no explican un aumento razonable como lo hacen individualmente.

El rendimiento del frijol se debió a la I_a, K_s, Lm, MO, Ac, θ_a, θ_s, θ_r, S_o y ρ_a, con un rango de 65.2 a 8.5% y las variables restantes explican individualmente menos del 8% (Figura 2c). En la misma figura se muestra la línea del porcentaje acumulativo de las variables del suelo que explicaron el 100% el rendimiento del cultivo del frijol, en este caso existen pendientes marcadas, por ejemplo, de I_a a Lm se encuentra el mayor porcentaje de explicación (85.5%); además, se observan tramos de MO a θ_a, y θ_r, aρ_a, con menor porcentaje de explicación, pero hay tramos de Lm a MO y θ_a a θ_r, donde no hay pendiente del rendimiento.

Con la inclusión de dos o más propiedades del suelo como variables independientes y el rendimiento de grano del haba, maíz y frijol como variable dependiente en las regresiones múltiples la R² mejoró de 0.64 a 0.97. La MO, K_s y Ac, explicaron 80% el rendimiento de haba, pero al incluir Rm y Ar se explicó 85% y al incluir la θ_r se alcanzó 96% (Cuadro 3).

El rendimiento de grano del cultivo de maíz se influenció por K_s, MO y Ar, las cuales explican 64% dicho rendimiento, pero al incluir la I_a, aumentó a 73% y posteriormente al adherir a S_o, se alcanzó 94% (Cuadro 3). De la misma manera, las propiedades del suelo que influenciaron el rendimiento de grano de frijol fueron I_a, K_s y MO, variables que explican 71% el rendimiento, pero al incluir θ_a a la ecuación, aumento a 88% y finalmente al adherir Ar se alcanzó 97% (Cuadro 3).

De nueve ecuaciones de regresión múltiple que se muestran en el Cuadro 3, en ocho de ellas el coeficiente de la MO fue positivo, lo cual indicó que el rendimiento de grano de los cultivos estudiados aumenta conforme el contenido de MO fue mayor en el suelo. Cada una de las variables del suelo tuvieron una influencia diferente para explicar el rendimiento de los cultivos, pero en general la K_s, MO y Ar fueron tres de 14 variables del suelo que estuvieron presentes en las regresiones múltiples con el rendimiento de grano de los tres cultivos (haba, maíz y frijol). Otras variables que influenciaron en menor importancia en las ecuaciones de regresión múltiple fueron: Ac, Rm, y θ_r para explicar el rendimiento del haba, así también I_a y S_o influenció en el rendimiento del maíz, del mismo modo I_a y θ_a en el rendimiento de fríjol (Cuadro 3).

Conclusiones

La materia orgánica, conductividad hidráulica saturada y contenido de arcilla fueron las variables que definieron el rendimiento del haba. El rendimiento de maíz influenciado por la conductividad hidráulica saturada, materia orgánica y contenido de arena. El rendimiento de frijol está asociado con la infiltración acumulada, conductividad hidráulica saturada y materia orgánica. Por lo tanto, la conductividad hidráulica saturada y la materia orgánica (agua y nutrimento) son los indicadores del suelo que definen el rendimiento de los tres cultivos manejados con labranza de conservación.

Bandick, A. K. and Dick, R. P. 1999. Field management effect on soil enzyme activities. Soil Biol. Biochem. 31(11):1471-1479.         [ Links ]

Brejada, J. J.; Moorman, T. B.; Karlen, D. L. and Dao, T. H. 2000. Identification of regional soil quality factors and indicators: I. Central and southern high plains. Soil Sci. Ame. J. 64:2115-2124.         [ Links ]

Carter, R. 2002. Soil quality for sustainable land management: Organic matter and aggregation interactions that maintain soil functions. Agron. J. 94:38-47.         [ Links ]

Chong, S. K. 1983. Calculation of sorptivity from constantrate rainfall is infiltration measurement. Soil Sci. Soc. Ame. J. 47:627-630.         [ Links ]

Dexter, A. R. 2004. Soil physical quality part I. Theory, effect of soil texture, density and organic matter, and effect on root growth. Article in press. Geoderma 120:201-214.         [ Links ]

Doran, J. W. and Parkin, T. B. 1994. Defining and assessing soil quality in defining soil quality for a sustainable environment. In: defining soil quality for a sustainable environment. Doran, J. W.; Coleman, D. C, Bezdicek, D. F. and Stewart, B. A. (Ed.). SSSA. Spec. Publ. Núm. 35. SSSA/ASA, Madison, Wis. USA. 3-21 pp.         [ Links ]

Gregorich, E. G. and Carter, M. R. 1997. Soil quality for crop production and ecosystem health. Develop. Soil Sci. 25:125-165.         [ Links ]

Gutiérrez-Castorena, Ma. C. y Ortiz-Solorio, C. A. 1999. Origen y evolución de los suelos del ex-lago de Texcoco, México. Agrociencia. 33:199-208.         [ Links ]

Hillel, D. 1980. Fundamentals of soil physics. Academic Press, Inc. New York, USA. 413 p.         [ Links ]

Jury, W. A.; Gardner, W. R. and Gardner, W. H. 1991. Soil physics. John, W. Sons. New York. 328 p.         [ Links ]

Lal, R. 1994. Methods and guidelines for assessing sustainable use of soil and water resources in the tropics. Soil management support services. USDA-NRCS, Washington. 57 p.         [ Links ]

Lal, R. 1998. Soil quality and agricultural sustainability. In: soil quality and agricultural sustainability. Lal, R. (ed.). Ann Arbor Press, Chelsea, MI. 3-12 pp.         [ Links ]

Lal, R.; Kimble, J. M.; Follett, R. F. and Cole, C. V. 1998. The potential of us cropland to sequester carbon and mitigate the greenhouse effect. Ann Arbor Press, Chelsea, MI. 457 p.         [ Links ]

Larson, W. E., F. J. Pierce. 1991. Conservation and enhancement of soil quality in evaluation for sustainable land management in the developing world. Vol. 2. IB SRAM Proc. 12(2). Int. Board for soil resource and management. Bangkok, Thailand. 37-51 pp.         [ Links ]

Li, Y.; and Lindstorm, M. J. 2001. Evaluating soil quality-soil redistribution relationship on terraces and steep hillslope. Soil Sci. Soc. Am. J. 65:1500-1508.         [ Links ]

Mendham, D. S.; Smethurst, P. J.; Holz, G. K.; Menary, R. C.; Grove, T. S.; Weston, C. and Baker, T. 2002. Soil analysis as indicators of phosphorus response in young eucalyptus plantations. Soil Sci. Soc. Am. J. 66:959-968.         [ Links ]

Norton, L. D.; Figueroa-Sandoval, B. and Oropeza-Mota, J. L. 1995. Comparison of runoffand erosion form natural and reclaimed tepetate from the semiarid area East of Mexico city. Water Research and Management in Semiarid Environments. Tucson, Arizona. 15 p.         [ Links ]

Rawls, W. J.; Gish, T. J. and Brakensiek, D. L. 1991. Estimating soil water retention from soil physical properties and characteristics. Adv. Soil Sci. 16:213-234.         [ Links ]

Salchow, E. R.; Lal, R.; Fausey, N. R. and Ward, R. 1996. Pedotransfer functions for variable alluvial soils in southern Ohio. Geoderma 73:165-181.         [ Links ]

Sanchez-Maranon, M.; Soriano, M.; Delgado, G. and Delgado, R. 2002. Soil quality in mediterranean mountain environments : Effects of land use change. Soil Sci. Soc. Am. J. 66:948-958.         [ Links ]

Statistical Analysis System (SAS V8.0). 1999. The SAS system for window. SAS Institute Inc. Cary, NC27513, USA. 100-120 pp.         [ Links ]