nueva página del texto (beta)
Español (pdf)
Artículo en XML
Referencias del artículo
Traducción automática
Enviar artículo por email
Citado por SciELO 
Similares en
SciELO 
Permalink
Introducción
La ganadería es actualmente el tipo de uso de suelo de mayor extensión en el mundo, ocupando aproximadamente un 26 % de la superficie terrestre libre de hielo (Steinfeld et al. 2006). En las regiones áridas, debido a su escasa precipitación y a la incertidumbre espacio-temporal de la misma, la ganadería ha sido históricamente la principal opción de aprovechamiento antropico de los recursos naturales (Le Houérou 2006).
Alrededor del 40 % del territorio mexicano son regiones áridas o semiáridas, y más del 60 % del mismo está dedicado a la ganadería (SEMARNAT 2013). En los 1970s, el gobierno federal promovió un ambicioso plan de desmonte y establecimiento de praderas, debido a la escasez de forraje natural como la principal limitante de la producción ganadera (Brenner 2011). La estrategia se basó en la introducción de especies de pastos de mayor productividad, como el zacate buffel (Cenchrus ciliaris L.) (Moreno-Unda 2011).
La baja precipitación, es el principal factor limitante de la producción primaria de las regiones áridas, la cual tiene relación con la producción (Golodets et al. 2013). Como consecuencia de ello, la variabilidad interanual en la precipitación, está correlacionada con variabilidad interanual de la producción primaria (Hsu et al. 2012). Por tanto, la limitación productiva de estas regiones no sólo es consecuencia de la escasa precipitación, sino también del efecto de la variación interanual de la misma (Le Houérou et al. 1988, Swemmer et al. 2007).
Se ha observado que la productividad de las praderas establecidas en desmontes con el paso del tiempo se reduce, particularmente en regiones áridas (Radford et al. 2007). Esta disminución temporal determina la sostenibilidad de la producción, las causas de la disminución productiva son complejas, pero se han relacionado con el declive productivo intrínseco de la reducción progresiva de los nutrientes disponibles en el suelo, especialmente el nitrógeno (Dalai et al. 2013, Morales-Romero et al. 2015) y el fósforo (Morales-Romero et al. 2015), y por la degradación productiva debido al manejo, especialmente por sobre-pastoreo (Búrquez y Martínez-Yrízar 2006). La disminución de cobertura vegetal reduce la materia orgánica del suelo, lo que a su vez limita la retención y disponibilidad de los nutrientes y aumenta el riesgo de erosión, especialmente en suelos pobres en materia orgánica como los de zonas áridas (Loveland y Webb 2003).
Las estimaciones de la productividad de praderas de zacate buffel se basa con frecuencia en la relación de la precipitación - producción en condiciones experimentales durante breves periodos de tiempo, normalmente no más de dos años (Martin-R. et al. 1995, Velásquez-Valle et al. 2014). Estos experimentos generan una relación producción -precipitación que no incluye la variación de la relación causada por la variabilidad interanual de la precipitación o por el potencial decrecimiento productivo, debido al paso del tiempo. En ambientes áridos, que poseen elevada variabilidad interanual de precipitación, el uso de esta relación promedio como referencia para inferir la producción potencial del establecimiento de una pradera de zacate buffel podría resultar en una sobre-estimación. Por lo anterior el objetivo del presente trabajo fue analizar el efecto que tiene la inclusión de la variación interanual de la precipitación y el declive productivo con la edad de establecimiento sobre la estimación de la productividad potencial de las praderas de zacate buffel en ambientes áridos.
Materiales y métodos
Zona de estudio
El área de estudio se localiza en el extremo meridional del estado de Baja California Sur, México (Figura 1), en la transición entre la vegetación desértica y la árida subtropical (González-Abraham et al. 2010). La precipitación promedio anual varía entre 100 y 450 mm, con variación interanual en ciclos de 5.9 a 7.4 años (Díaz et al. 2001). A pesar de que Baja California Sur no se encontraba entre los 13 estados elegidos para implementar el PRONADE (Moreno-Unda 2011), se desmontaron y establecieron alrededor de 17 000 ha de praderas de zacate buffel para uso ganadero extensivo (Guillén-Trujillo et al. 2007, SEMARNAT 2011). Para luego establecer alrededor de 5 000 ha de praderas (Fundación Produce Baja California Sur 2011).
Figura 1 Localización geográfica de las 12 praderas de zacate buffel estudiados (círculos negros) y de las estaciones meteorológicas (círculos grises). La curva de nivel corresponde a los 600 m de altitud.
Con base en cartas de uso de suelo y vegetación del INEGI (series I, II y III) y por confirmación mediante visitas en campo en el 2012, se seleccionaron 12 praderas de entre 20 y 35 años de edad, con extensión superior a 100 ha distribuidas en un gradiente de 100 a 450 mm de precipitación (Tabla 1). Para cada una de las praderas se eligió la estación meteorológica de la Comisión Nacional del Agua (CONAGUA) más cercana que contase con al menos con 20 años de registros.
Modelo de estimación productiva
Se analizó el efecto de incorporar la variación interanual de la precipitación y el tiempo desde el establecimiento de la pradera sobre la estimación de la productividad, por medio de tres modelos en los que se incorporaron de forma escalonada los factores de: (1) precipitación anual constante, (2) variabilidad interanual de la precipitación, y (3) efecto del tiempo desde el establecimiento.
Modelo 1: No variación, precipitación anual constante. En este modelo la biomasa producida es función de la precipitación promedio del sitio. La única variación productiva es entre sitios, y esta causada por la heterogeneidad geográfica de la precipitación. En ambientes áridos, donde el agua es el principal factor limitante de la producción, se acepta la relación lineal entre la biomasa anual producida (B1) y la precipitación anual (P) (Hsu et al. 2012), donde a es la intersección de la línea con el eje de ordenadas y b la pendiente de la relación entre ambas variables, con el modelo: B1 = a + b · P.
En el presente estudio se utilizó la relación entre producción anual de zacate buffel y precipitación obtenida por Martin-R et al. (1995), con valores de B1 = -1 527.9 + 19.8 · P v , donde B1 indica la biomasa anual producida en kg ha-1 y P v la precipitación promedio de verano en mm. Considerando como verano aquellos meses en los que el zacate buffel puede tener crecimiento activo, con la temperatura mínima mensual promedio superior a 15 °C (Cox et al. 1988). Para luego calcular la precipitación promedio de verano (P v ) y la producción promedio (B1). Cuando la producción promedio estimada de una pradera fue inferior a cero, la pradera se consideró inviable y se le adjudicó una producción de cero.
Modelo 2: Variación interanual de la precipitación. La biomasa producida está relacionada no sólo con el promedio de la precipitación, sino también con la frecuencia de la misma (Le Houérou et al. 1988, Guevara et al. 1996). La posibilidad de incorporar la variabilidad interanual de la precipitación se realiza con la estimación de la productividad anual promedio mediante un umbral de precipitación confiable (Pc), con modelo: B2 = a + b · Pc.
Para la cual se utiliza una frecuencia umbral de 80 % (Le Houérou et al. 1988), es decir, que en cuatro de cinco años la precipitación será igual o superior a un determinado valor, que corresponde al percentil 20, denominado precipitación confiable. El índice de confiabilidad es el cociente entre la precipitación confiable y la precipitación promedio (Guevara et al. 1996), por lo que varia con el gradiente de precipitación. Se utilizaron registros históricos de las estaciones correspondientes para calcular la precipitación de verano confiable al 80 % (Pcv_0.8) de cada pradera y estimar su producción anual promedio confiable (B2).
Modelo 3: Variación interanual de la precipitación y efecto de la edad de las praderas. Se ha observado de forma recurrente la disminución productiva de praderas de zacate buffel al incrementarse el tiempo de establecimiento de la misma (Radford et al. 2007). Por lo que la producción esperada en el año t desde el establecimiento de la pradera es: B3(t) = B2 - D(t) donde B2 corresponde a la biomasa promedio estimada según el modelo dos, y D(t) representa la función que expresa la disminución productiva en función de la edad de la pradera (t). La pradera deja de ser viable productivamente cuando la reducción productiva por el paso del tiempo (D(t)) sea igual a la capacidad productiva potencial (B2). Por lo que se espera que la función D(t) varíe entre sitios de acuerdo a las condiciones ambientales y el manejo productivo de cada pradera. Sin embargo, dado el carácter exploratorio del presente estudio se consideró de manera conservadora que D(t) tuviera el mismo valor para todas las praderas y fuera igual a ln(t) (Radford et al. 2007). Por lo que se estimó la producción promedio de cada pradera con base en la ecuación del modelo para dos edades de establecimiento: (a) 10 años, y (b) la edad mínima estimada de cada pradera.
Para caracterizar la condición productiva actual de las praderas, se estimó su densidad mediante el establecimiento en cada una de las 12 praderas seleccionados dos parcelas de 200 χ 500 m, en las que se seleccionaron cinco puntos al azar. En cada punto se dispuso un transecto de 50 x 2 m con orientación aleatoria, y se establecieron 6 cuadros de 2 x 1 m a intervalos de 10 m en el interior del transecto. En total se mu est rea ron 12 praderas, por medio de 120 transectos y 720 cuadros. En cada cuadro se registró el número de individuos (macollos) de zacate y se estimó la densidad promedio de cada pradera.
Resultados
Efecto de la variación interanual de la precipitación en la estimación productiva
La extensión de las 12 praderas estudiadas varió entre 103 ha y 1 648 ha, con altitud entre 55 m y 470 m (Tabla 1). Las estaciones meteorológicas seleccionadas se encuentran a menos de 10 km de las praderas, excepto en el sitio 12 que está a 24 km (Figura 1). La precipitación anual promedio varió entre 107 y 465 mm, con cinco estaciones por debajo de 250 mm anuales. El número de meses de crecimiento activo potencial entre las distintas praderas fue de 3 a 8 meses, y la precipitación promedio de verano observada en los 12 sitios varió entre 65 y 413 mm (Tabla 1). La variación interanual de la precipitación de verano a nivel de sitio fue elevada con CV entre 58 y 241 %. La precipitación de verano confiable (P cv0.8 ) entre los sitios varió entre 8 y 276 mm y el índice de confiabilidad se redujo al disminuir la precipitación promedio de verano (Figura 2).
Figura 2 Relación entre la precipitación promedio de verano de las 12 praderas de zacate buffel estudiadas y su variabilidad interanual expresada en el índice de confiabilidad (cociente entre la precipitación de verano confiable al 80 % y la precipitación de verano promedio).
El escenario de productividad estimada cambió al incorporar la variabilidad interanual de la precipitación en el modelo (Figura 3). En el modelo uno la precipitación de verano que marca el umbral de viabilidad productiva es de 77 mm, tuvo una pradera no viable y producción promedio de las once praderas restantes de 790 a 6 653 kg ha-1. La inclusión de la variabilidad interanual de la precipitación de verano en el modelo 2 incrementó el umbral de viabilidad productiva a 171 mm, ya que la regresión lineal de la producción respecto a la precipitación promedio de verano fue B2 = -2663 + 15.57 · P v . En este modelo las praderas inviables fueron cinco y la producción de las siete praderas restantes fue de 712 a 3 946 kg ha-1. La producción estimada a nivel de pradera en el modelo 2 fue en promedio un 65 % menor de la obtenida con el modelo 1.
Figura 3 Producción de las praderas de zacate buffel estimada por los tres modelos en relación a la precipitación promedio de verano. Los modelos se diferencian en el grado de variabilidad temporal que incorpora en sus variables. Círculos blancos: modelo 1, no considera variación temporal, productividad es función de la precipitación promedio de verano; círculos grises: modelo 2, considera variación interanual de la precipitación (precipitación de verano confiable al 80 %); y círculos negros: modelo 3, igual al modelo 2 más la variación productiva por la edad de establecimiento.
Efecto de la edad de establecimiento en la estimación productiva
La condición productiva observada fue muy precaria. La densidad observada en 11 de las 12 praderas varió entre 0.17 y 2.7 individuos m-2, presentando una productividad de 5.38 individuos m-2 la pradera restante. El modelo 3 tiene el mismo número de praderas productivas que el modelo 2, y la viabilidad temporal de estas praderas varió entre 2 y 40 años. Con base en las fechas estimadas de establecimiento solo dos praderas siguen siendo viables. La producción promedio a 10 años de su establecimiento varió entre 73 y 2 435 kg ha-1, y se observó una reducción promedio a nivel de pradera respecto al modelo 2 del 63 %. Si se considera el tiempo desde su establecimiento hasta la actualidad, la producción promedio anual entre las praderas varía entre 19 y 1 236 kg ha-1, y la reducción promedio por pradera alcanza el 82 %.
Discusión
La región de estudio presenta los atributos característicos de las regiones áridas con bajo promedio anual de precipitación y elevada variación espacial y temporal. Si consideramos que la precipitación de verano es la funcionalmente importante para el zacate buffel, la variación espacial y temporal de la misma aumenta más. La producción anual promedio observada como función de la precipitación de verano promedio (modelo 1) es moderada y presenta una elevada heterogeneidad entre las praderas viables, de 790 a 6 653 kg ha-1; con CV1 de 69 %. Los resultados muestran que la inclusión de la variabilidad interanual de la precipitación en el modelo 2 tiene tres consecuencias generales sobre la productividad estimada. En primer lugar, aumenta el valor umbral de la precipitación de verano necesaria para obtener productividad viable, que pasa de 77 a 171 mm en la región de estudio. Esto implicaría una reducción a nivel regional del área con suficiente precipitación para garantizar la viabilidad productiva. En segundo lugar, se reduce el valor estimado de la producción anual promedio de cada pradera y, finalmente, aumenta la variabilidad geográfica de la producción estimada. El incremento de la heterogeneidad geográfica de la precipitación al incorporar la variabilidad temporal de la misma en el modelo 2 incrementa del 51 al 68 %, mientras que en la heterogeneidad geográfica de la producción pasa del 69 al 108 %. Este traslado amplificado de la variabilidad de la precipitación en ambientes áridos hacia la variabilidad espacial y/o temporal de la producción de forraje coincide con lo reportado por Lauenroth y Sala (1992) quienes encontraron que el coeficiente de variación temporal de la producción de forraje en una pradera natural del norte de Colorado fue mayor que la variación de la precipitación asociada. Le Houérou et al. (1988) reportan variabilidad de la producción entre 20 y 50 % superior a la variabilidad en la precipitación en las áreas de pastoreo del Mediterráneo árido.
La estimación productiva del modelo 2 incluye el efecto de la variabilidad productiva debido a la variación interanual de la precipitación. Sin embargo, hay que considerar que la estimación de este modelo se refiere a la producción promedio para un periodo determinado, en el que hay años con valores de producción por arriba y por debajo del promedio. Ante este reto de variabilidad productiva de las regiones áridas se ha sugerido como respuesta incorporar a la estrategia de manejo del ganado la movilidad, la carga de pastoreo variable y el manejo adaptativo (Vetter 2005, Pinedo et al. 2013). Entre las respuestas que se han venido estableciendo ante los periodos de sequía, habría que revisar el efecto a largo plazo del suplemento alimenticio de animales en pastoreo libre, ya que implica mantener la presión sobre la pradera deteriorada por la sequía (Von Wehrden et al. 2012). Cuanto más se deteriore un lugar, más le costará recobrar su producción cuando vuelva a existir precipitación suficiente, pues el deterioro incrementa los riesgos de erosión (Búrquez y Martínez-Yrízar 2006) y reduce la eficiencia en el uso de la precipitación por parte del pasto (Blanco et al. 2005). Por ejemplo, Pinedo et al. (2013) estimaron que después de un año de intensa sequía, como el 2011 en Chihuahua, el sobre-pastoreo acentuó el efecto de la misma y los pastizales naturales requerirán de 4 a 6 años para recuperar su productividad normal. Por otro lado, los bajos valores de densidad observados sugieren que el segundo factor de variación temporal productiva considerado en el análisis, el tiempo de establecimiento de las praderas, ha tenido un efecto significativo en la capacidad productiva de las mismas. Mientras que Guillén-Trujillo et al. (2007) encontraron baja densidad en cinco praderas que en el presente estudio mostraron valores incluso inferiores de densidad. Sin embargo, las diferencias se pueden deber al efecto de los diferentes muéstreos en ambos estudios.
Este declive productivo se ha observado frecuentemente en praderas establecidas después del desmonte y se ha explicado en parte como resultado de la disminución progresiva del Ν disponible en el suelo después del desmonte (Dalai et al. 2013, Morales-Romero et al. 2015, Celaya-Michel et al. 2015). Algunos autores sugieren que el declive productivo es la tendencia esperada hacia un nuevo estado de equilibrio en la dinámica de nutrientes de la pradera tras el desmonte de vegetación leñosa, y los altos valores productivos en los primeros años estarían relacionados con el incremento temporal en la disponibilidad de nitrógeno (Dubeux et al. 2007, Radford et al. 2007). También se ha señalado que una alta intensidad de pastoreo combinado con ausencia de fertilización externa puede producir una extracción de nutrientes del sistema, especialmente de nitrógeno, a través de la biomasa animal (Haynes y Williams 1993). El sobre-pastoreo puede producir también reducción de la cobertura vegetal que al disminuir la tasa de infiltración provoca un aumento del riesgo de pérdida de suelo, y por tanto de nutrientes, por erosión (Loveland y Webb 2003). La siembra de leguminosas que fijen N2 se ha propuesto corno una opción para paliar este declive productivo asociado a la disminución de Ν disponible (Dalai et al. 2013), idealmente debe combinarse el manejo de la densidad de animales en el tiempo con el objeto de regular la extracción de nitrógeno del sistema (Silcock y Johnston 1993).
Conclusiones
Los resultados sugieren reducción en la estimación de la producción promedio proyectada en praderas de zacate buffel, si se incluyen los factores de variabilidad temporal respecto al valor obtenido solo como función del promedio de la precipitación. La producción a futuro de las praderas de zacate buffel, actuales o por establecerse en las regiones áridas, debería considerarse con base no sólo en el valor promedio de la precipitación, sino también en la variabilidad interanual de la precipitación y el efecto de la edad de establecimiento.
