Spanish (pdf)
Article in xml format
Article references
Automatic translation
Send this article by e-mail
Cited by SciELO 
Cited by Google 
Similars in
SciELO 
Similars in Google 
Permalink
Introducción
La ganadería de cría extensiva se distribuye a nivel mundial en zonas semiáridas y áridas, caracterizadas por temperaturas elevadas y de reducida calidad y disponibilidad de agua y forraje, donde no es posible desarrollar con éxito otras actividades agropecuarias, sin embargo, estos sistemas de producción están caracterizados por una baja carga animal y bajos índices de procreo (Viglizzo, 2014).
Algunas variables, como la carga ganadera, la ubicación de las aguadas, la sombra y las características del pastizal (terreno y vegetación), tienen un impacto significativo en el comportamiento del ganado en pastoreo. Cuando los pastizales son extensos y las fuentes de agua escasas, el ganado sobrepastorea las zonas más cercanas a las aguadas, lo que provoca un deterioro del estado del pastizal y un desaprovechamiento de las zonas alejadas de estas (Blanco et al., 2008; Rivero et al., 2021).
Aunque el ganado puede recorrer grandes distancias en busca de forraje, su eficiencia de cosecha disminuye notablemente más allá de los tres km desde la fuente de agua. Esto es especialmente notorio en razas británicas durante veranos muy calurosos (Herrera Conegliano et al., 2022; Hunt et al., 2014; Nyamuryekung’e et al., 2021).
El estrés térmico puede afectar severamente la tasa de fecundación y la supervivencia de embriones en razas no adaptadas a elevadas temperaturas. Esta condición compromete la esteroidogénesis, reduce la viabilidad de los ovocitos y favorece la mortalidad posnatal por hipertermia en genotipos no adaptados (Brown-Brandl, 2018; Kassab et al., 2020). Por su parte, el balance energético negativo agudo es una de las principales causas de mortalidad embrionaria en rodeos de cría.
El fuego es un importante factor ecológico en muchos ecosistemas de climas semiáridos-áridos, influyendo en la estructura y la composición de la vegetación; donde este suele ser utilizado para gestionar la vegetación de las zonas de pastoreo: controlar la vegetación leñosa y mantener un sotobosque productivo, eliminar el material herbáceo senescente para facilitar su rebrote con mayor calidad y reducir el riesgo de incendios forestales destructivos (Briske, 2017).
Se sabe que Argentina experimenta un proceso de intensificación del sistema productivo ganadero durante los últimos 30 años, por lo que los cambios ocurridos en su cría bovina son el resultado de la competencia por la tierra con otros productos agropecuarios, el deterioro de los recursos naturales, la lenta adopción de tecnología, la rentabilidad y la incapacidad de generar un conjunto de políticas útiles para sostener y mejorar la productividad (Arelovich et al., 2011). Desarrollada en tres cuartas partes del territorio nacional, bajo condiciones áridas o semiáridas, la ganadería de cría estancó su porcentaje de destete en un 60% a nivel nacional, valor que es superado por otros países en condiciones ambientales similares (Arceo, 2017; Viglizzo, 2014).
La provincia de la Pampa dispone de diez millones de hectáreas dedicadas principalmente a esta actividad, correspondientes a la región del Espinal (semiárido) y el Monte Occidental (árido). El porcentaje de destete provincial es similar a la media histórica nacional (60%), con valores promedios de 59% para el semiárido y 50% para el árido (Senasa, 2021), donde estas dos regiones se caracterizan por la utilización de grandes parcelas de pastoreo (> 500 ha), con escasos puntos de distribución de agua, pobres en cantidad y calidad (Umazano et al., 2005). Debido a ello, el pastoreo afecta la heterogeneidad espacial, generando extensas áreas subutilizadas y zonas cercanas a las aguadas muy degradadas por defoliación extrema. El proceso de arbustización, provocado en parte por el inadecuado manejo del pastoreo, afecta casi a cuatro millones de hectáreas (Vázquez et al., 2013) y esta situación es una severa limitación para el desarrollo de la ganadería de cría, con impacto directo sobre la capacidad de carga animal (Blanco et al., 2008).
La normativa provincial vigente sobre la conservación de bosques nativos, Ley 2624 de 2012, estableció estrategias de intervención controlada sobre todas las áreas definidas en la categoría II del ordenamiento territorial de los bosques nativos. La presencia de doce eventos climáticos adversos en La Pampa, entre 2008 y 2021 (Vázquez, 2021), está asociada a un proceso de cambio climático detectado en la región central de Argentina, caracterizado por un aumento significativo de la temperatura y precipitaciones extremas (Basile et al., 2021). Esto representa un importante reto a superar, ya que los modelos más conservadores pronostican, para el periodo 2015-2099, eventos de estrés calóricos en bovinos cada vez más frecuentes para la región (Armendano et al., 2020).
Por su parte, los acuíferos de la región son generalmente de poco espesor, con escaso rendimiento y de regular o mala calidad para la actividad pecuaria (Umazano et al., 2005), lo que limita severamente la posibilidad de aumentar la carga animal en un área de aproximadamente un millón de hectáreas (Adema, 2015).
Todos estos factores limitantes interactúan entre sí y condicionan el desarrollo del sistema criador, donde la productividad es controlada por la disponibilidad del recurso más escaso. Además, solo se lograrán progresos significativos cuando la tecnología aplicada considere la jerarquía de estas limitantes.
El objetivo de este trabajo fue analizar la asociación existente entre los factores: agua, forraje, carga animal, temperatura, arbustización y fuego, con la tasa de destete del sistema productivo a escala regional y cuantificar la magnitud de sus impactos en la región criadora de la provincia de La Pampa, Argentina. Resultado de la investigación fue notar que la producción de terneros es afectada por limitantes que actúan de manera jerárquica y estas deberán ser atendidas en ese orden para lograr impactos productivos positivos en la región.
Materiales y métodos
Selección de establecimientos
Se realizó una selección de 180 establecimientos al azar (15% del total) en la zona criadora de la provincia de la Pampa, Argentina (-37° latitud, -66° longitud), caracterizados por un ambiente semiárido-árido en el sentido este-oeste. Para representar correctamente los establecimientos con destetes promedios superiores al 65%, se estimó el tamaño muestral según el cálculo de razón de probabilidad con tamaño poblacional conocido (Lwanga & Lesmehow, 1991):
(1)
Donde:
n = tamaño de muestra mínimo (156)
N = tamaño población objetivo (5000)
Z = nivel de confianza (1,95)
P = probabilidad a favor (0,12) (estadísticas del Senasa 2004-2017)
1- P = probabilidad en contra (0,88)
e = error de estimación (0,05)
Una condición fue que todos los establecimientos estuvieran identificados con sus números de Registro Nacional Sanitario de Productores Agropecuarios (Renspa). Esta identificación facilitó la reconstrucción de la información climática y productiva seriada del periodo 2004-2017 (figura 1). Así, la caracterización climática regional se analizó a partir de series temporales completas de temperatura del aire del National Oceanic and Atmospheric Administration (NOAA, 2020) y pluviométricas del Global Precipitation Climatology Centre (gpcc) (Schneider et al., 2018; Ziese et al., 2014) en grillas de 0,5°, correspondientes al periodo 1980-2020.
Fuente: Elaboración propia
Figura 1. Distribución espacial de los establecimientos analizados en la zona criadora de la provincia de La Pampa
Las relaciones ternero/vaca registradas durante las campañas de vacunación antiaftósica (periodo 2004-2017) fueron corregidas a partir de las declaraciones de movimiento de hacienda (Senasa, 2021), por lo que se seleccionaron 10 variables relacionadas con la disponibilidad de agua, forraje, temperatura y arbustización.
Calidad de agua (Cal_agua)
Se generó un mapa de calidad de agua para bebida de ganado con base en la información físico- química de 1850 perforaciones realizadas en toda la región (Instituto de Hidrología de Llanuras & Secretaría de Recursos Hídricos, 2021) según normas de calidad establecidas por el Instituto Nacional de Tecnología Agropecuaria, INTA (Sager, 2001), donde:
Sales totales=residuo seco(gl−1) (2) y
Sales beneficiosas(gl−1)=Sales totales − Sulfatos(SO4−2) (3)
Así, se complementó la información con las concentraciones halladas de arsénico (As) y flúor (F) y el agua con valores superiores a 0,3 y 2,0 ppm, respectivamente, fueron consideradas no aptas.
Precipitación acumulada en otoño-invierno-primavera y verano-otoño (PPT OIP, PPT VO)
Se analizó la precipitación previa y durante la fecha de servicio provista por un modelo de precipitación global (Schneider et al., 2018) validado para la provincia (Vázquez et al., 2013).
Porcentaje del establecimiento incendiado (% de fuego)
Se cuantificó la superficie afectada por fuego durante la temporada de incendios (noviembre- marzo) en cada establecimiento, a partir de digitalización directa sobre imágenes “near real time” (NRT) Modis Terra y Aqua de 250 m de resolución espacial, bandas 7-2-1, suministradas por la NASA (https://lance.modaps.eosdis.nasa.gov/imagery/subsets/?subset=LaPampa).
Presencia de arbustización (arbustización)
Se utilizó un modelo desarrollado para la provincia de La Pampa que identificó procesos arbustivos en unidades de 260 hectáreas (Vázquez et al., 2016), basado en el algoritmo Seasonal Trend Decomposition by Loess (Lu, 2003), y se estimó el coeficiente de área arbustizada del establecimiento como:
ArbPred = arbSTL (ha)⁄establecimiento (ha) (4), donde:
ArbPred = coeficiente de arbustización del establecimiento.
arbSTL (ha) = hectáreas clasificadas con procesos arbustivos por algoritmo STL. Establecimiento (ha) = hectáreas del establecimiento.
Para ajustar estos resultados a la situación predial, se analizó mediante el procedimiento no paramétrico Mann-Kendall (Terrset, 18.3 ©Clark University, 2017) dentro de estas unidades y la tendencia de la cobertura de bosques a partir de la serie temporal (anual) del producto MOD44B del sensor Modis (Hansen et al., 2002), ajustada a los valores de los inventarios de bosque nativos realizados en 2018 (Roberto & Carreño, 2018). Se consideró como valor de corte Z ≥ |1,95| (α = 0,05). Finalmente, se definió arbustización como:
Arbustización=ArbPredZ⁄medio anual (5), donde:
Arbustización ≥ 1: hay proceso arbustivo (≈ ≥ 30% del establecimiento con arbustización). Arbustización < 1: no hay proceso arbustivo importante.
Carga térmica ambiental (estrés)
Se estimó el estrés en todos los establecimientos monitoreados durante los 14 años estudiados a partir de temperaturas de superficie (Ts), determinadas a partir de sensores remotos (Niclòs et al., 2010). La Ts posee la capacidad de reflejar la influencia de la humedad relativa, la radiación incidente y la velocidad del viento (Sánchez et al., 2007), parámetros involucrados en el cálculo de indicadores de estrés calórico, como el índice de Temperatura-Humedad (THI) y el índice de carga calórica (HLI). Ello se estimó como:
(6)
donde:
Tst (e-f): temperatura de superficie media del establecimiento, periodo enero-febrero, durante el año analizado.
Ts2002-2020 (e-f): temperatura de superficie media histórica del establecimiento, periodo enero- febrero (2002-2020).
Además, se utilizaron composiciones de ocho días de temperatura de superficie del sensor Modis Aqua (producto MYD11A2) con resolución espacial de 1 km², para el periodo 2002-2020.
Distancia a las aguadas (distancia)
Se censaron los pozos de agua (visualización de tanques australianos) y tajamares mediante imágenes de alta resolución espacial (0,5-2,0 m), disponibles entre los años 2010-2017 (©SAS.Planet, 2017) y se eliminaron del análisis los tanques secos. Una vez censadas las fuentes de agua (figura 5), se calculó la distancia media de cada establecimiento a las mismas mediante el procedimiento Varcost y Distance (Terrset, 18.3 ©Clark University, 2017) y se consideró que 2500 m es una distancia crítica a partir de la cual existe aprovechamiento mínimo del forraje para equilibrar posibles errores de omisión de aguadas (Bavera, 2011).
Balance hídrico (BalHid)
Se calculó la carga animal promedio anual, expresada en equivalente vaca (e.v.) (Cocimano et al., 1975) y se asumió un consumo de 70 l de agua diario por e.v., valor que refleja la mayor demanda diaria por animal adulto de razas británicas durante la lactancia o el periodo de calor intenso (Bavera, 2011). La demanda anual de agua del establecimiento se definió como:
Demanda de agua anual (l)=e.v.establecimiento∗70le.v⁄.∗365 días (7)
La oferta anual de agua se estimó a partir del número de pozos funcionales, acceso al agua de acueductos y tajamares censados en los establecimientos; mientras que la cantidad de agua aportada por día se estimó a partir de los caudales de extracción recomendados (Instituto de Hidrología de Llanuras & Secretaría de Recursos Hídricos, 2009). Por su parte, el aporte de los tajamares se calculó estimando el volumen de cada uno, la precipitación registrada, la eficiencia de captación y la evapotranspiración local estimada por sensores remotos (Rivas, 2004). En promedio, estos tienen una dimensión de 60 x 40 x 2 m. y posteriormente se estimó el volumen de agua aportado como:
Oferta agua anual (l)=Σ[pozos∗caudal(lh⁄)∗bombeo diario(h)]+aguatajamares(l) (8)
Finalmente se definió BalHid como:
BalHid=Oferta de agua anual−Demanda de agua anual (9)
Balance forrajero (BalFor)
Se consideraron los requerimientos energéticos de los bovinos presentes en cada establecimiento y la oferta forrajera anual estimada a partir de los aportes del campo natural, pasturas y verdeos, expresados en e.v. (Cocimano et al., 1975), con lo que se modeló la producción forrajera anual del campo natural (oferta forrajera 1 kg MS ha/año) del año correspondiente al momento de servicio (Vázquez et al., 2016) (figura 2).
Fuente: Elaboración propia
Figura 2. Esquema de cálculo de oferta forrajera del campo natural por establecimiento (CN) considerando el momento de servicio como periodo crítico para afectar el destete
El aporte forrajero de cada recurso implantado se estimó considerando la eficiencia del uso del agua de los principales recursos forrajeros de la zona (Moreyra et al., 2014) y la precipitación registrada durante el ciclo de crecimiento (Schneider et al., 2018). Sobre una máscara de actividad agrícola, se definieron parámetros biométricos en 14 series temporales anuales de 12 composiciones de NDVI (Vázquez et al., 2013). Con ellos se identificó la superficie ocupada por pasturas y verdeos de invierno y verano (oferta forrajera 2).
Oferta forrajera (2)=E.U.A.(kgMSha⁄∗mmagua−1)∗ppt(mm)∗área(ha) (10)
Donde:
E. U. A.= eficiencia del uso del agua de pastura y verdeos
ppt= precipitación
y
BalFor=[Oferta forrajera 1+2](e.v.)−Demanda forrajera (e.v.) (11)
Análisis estadístico
Se analizó el grado de asociación existente entre factores ambientales y de manejo con las tasas de destete de terneros obtenidas en la región, mediante el test no paramétrico de tablas de contingencia (Whitley & Ball, 2002). Las variables fueron categorizadas en 0 y 1, considerándose (1) a las relaciones ternero/vaca < 0,65 y a todas las condiciones de factores ambientales y de manejo que las favorecen (tabla 1). Su asociación fue probada mediante el test Chi-cuadrado de Pearson, (χ2, α = 0,10).
Tabla 1. Sistematización y codificación de variables
Nota aclaratoriacalidad de agua (Cal_agua); precipitación acumulada en otoño-invierno- primavera (PPT OIP); precipitación acumulada en verano-otoño (PPT VO); porcentaje del establecimiento incendiado (% de fuego); presencia de arbustización (arbustización); carga térmica ambiental (estrés); distancia a las aguadas (distancia); balance hídrico (BalHid) y balance forrajero (BalFor).
Fuente: Elaboración propia
Para estimar la probabilidad de que una condición esperada esté presente en un grupo de la población frente al riesgo que ocurra en otro, se calculó el odds ratio de cada combinación significativa obtenida en la tabla de doble entrada (Cerda et al., 2013) como:
Odds ratio = (a⁄c)⁄(b⁄d), (12) donde:
a = cantidad de resultados deseados obtenidos a partir de la aplicación de un tratamiento.
c = cantidad de resultados deseados obtenidos sin la aplicación de un tratamiento.
b = cantidad de resultados no deseados obtenidos a partir de la aplicación de un tratamiento.
d = cantidad de resultados no deseados obtenidos sin la aplicación de un tratamiento.
El intervalo de confianza se calculó como:
𝐼𝐶(95 %)=𝑒(𝑙𝑛(𝑂𝑑𝑑𝑠 𝑟𝑎𝑡𝑖𝑜)±1,96∗(1𝑎⁄+1𝑏⁄+1𝑐⁄+1𝑑⁄)0,5) (13)
La probabilidad se estimó como en Whitley & Ball (2002):
Probabilidad (%)=Oddsratio(Odds ratio+1)⁄ (14)
Resultados
Caracterización de variables explicativas
Los parámetros climáticos observados en la provincia de La Pampa presentaron variaciones durante los últimos 40 años. Las temperaturas del aire máximas medias mensuales se incrementaron en 2 °C, la precipitación media mensual disminuyó 20 mm, se duplicó la frecuencia de días con temperaturas superiores a 35 °C y los valores de estrés hídrico trimestrales pasaron de un valor medio de 0,5 (moderadamente húmedo) durante la primera mitad del ciclo evaluado a -0,5 (moderadamente seco) en la segunda mitad (figura 3).
Fuente: Elaboración propia
Figura 3. Valores medios mensuales observados y evolución de parámetros climáticos (valor entre paréntesis) en la zona de cría de la provincia de la Pampa, periodo 1980-2020.Notas aclaratorias: a) temperatura máxima del aire, b) precipitaciones, c) frecuencia de días con temperaturas > 35 °C y d) índice de sequía trimestral.
La carga térmica del periodo enero-febrero, analizado a partir de Ts, afectó a nueve de los diez millones de hectáreas de la región, con tendencia ascendente durante el periodo 2002-2020 (figura 4).
Factores limitantes de la relación ternero/vaca
La información obtenida a partir de los censos de fuentes de agua permitió reconstruir la calidad, la disponibilidad de agua y el balance hídrico a escala predial (figura 5).
La relación ternero/vaca < 0,65 estuvo asociada significativamente a cargas térmicas ambientales elevadas (estrés = 1), precipitaciones estivales por debajo de la media histórica (PPT VO = 1) y balances hídricos del rodeo negativos (BalHid = 1) (tabla 2).
Fuente: Elaboración propia
Figura 5. Reconstrucción de factores asociados a la disponibilidad de agua a escala predial en la provincia de La PampaNota aclaratoria: a) censo de tanques australianos y tajamares, b) caudales de pozos (l/h), c) calidad de agua para el ganado según la escala de Sager (2001) y d) distancia a las aguadas.
Tabla 2. Odds Ratio (OR), intervalo de confianza (IC) y probabilidad (p) de ocurrencia de una relación ternero/vaca (Te/vc) < 0,65 en la provincia de La Pampa, en función de las variables analizadas
Nota aclaratoriacalidad de agua (Cal_agua); precipitación acumulada en otoño-invierno- primavera (PPT OIP); precipitación acumulada en verano-otoño (PPT VO); porcentaje del establecimiento incendiado (% de fuego); presencia de arbustización (arbustización); carga térmica ambiental (estrés); distancia a las aguadas (distancia); balance hídrico (BalHid) y balance forrajero (BalFor).
*α=0,10
**α= 0,05
***α=0,001
nsno significativo
Fuente: Elaboración propia
La combinación estrés = 1, BalHid = 1 y BalFor =1 tuvo una probabilidad del 83 % de obtener relaciones ternero/vaca < 0,65 (tabla 3).
Tabla 3. Odds ratio (OR), intervalo de confianza (IC) y probabilidad (p) de ocurrencia de una relación ternero/vaca (Te/vc) < 0,65 en función de la carga térmica (estrés), balance hídrico (BalHid) y balance forrajero (BalFor) en la provincia de La Pampa
*α=0,10
**α= 0,05
***α=0,001
nsno significativo
Fuente: Elaboración propia
La variable distancia no estuvo asociada a relaciones ternero/vaca < 0,65, pero permitió explicar el proceso arbustivo presente en la zona de cría. Los establecimientos con distancias medias a las aguadas mayores a 2500 m tuvieron una probabilidad de 66 % de estar afectados por procesos arbustivos y de 80 % si hubo incendios (tabla 4).
Discusión
Las altas temperaturas estivales presentes en la región (T > 35 °C) afectaron significativamente el rendimiento productivo de los bovinos y la incorporación de genotipos más rústicos o la suplementación de animales con selenio y zinc son opciones que han permitido incrementar las tasas de destete en un 25 % en condiciones similares (Burns et al., 2010; Nyamuryekung’e et al., 2021; Kassab et al., 2020).
Por su parte, las precipitaciones estivo-otoñales afectan directamente la producción forrajera y la persistencia de los pastizales. La mayor frecuencia de periodos secos y días con temperaturas elevadas halladas en este trabajo son situaciones características del proceso de cambio climático instaladas en la región (Armendano et al, 2020; Basile et al., 2021), coincidentes con lo hallado en otras zonas de pastoreo a escala mundial (Briske, 2017).
El balance hídrico animal en la región de estudio está ajustado a una demanda promedio de 35 l de agua por animal por día, valor característico en razas británicas. Esto es insuficiente si se consideran los 65 l requeridos durante los periodos de máximo requerimiento, afectando negativamente la capacidad de servicio, fertilidad y lactancia (Wagner & Engle, 2021).
El análisis espacial indicó que siete de las diez millones de ha están a más de 2500 m de distancia de una fuente de agua. Esto significa que 70% de la superficie no es pastoreada eficientemente y el 30% restante es sometida a sobrepastoreo, favoreciendo los procesos arbustivos (Blanco et al., 2008; Vázquez et al., 2016). Esta característica es una condición frecuente en áreas de pastoreo de zonas áridas y semiáridas de América, África y Australia (Briske, 2017).
Las razas índicas y criollas pueden desplazarse a mayores distancias en busca de forraje durante los periodos de altas temperaturas, comparadas con las razas británicas. Su locomoción y su capacidad de desplazarse a mayores distancias es una propiedad heredable en los bovinos, característica que permitiría desarrollar mayores índices de aprovechamiento de las parcelas de pastoreo (Nyamuryekung’e et al., 2021; Rivero et al., 2021).
Es de considerar que las altas temperaturas estivales y la escases de agua imponen la necesidad de evaluar e incorporar tecnologías no consideradas aún en la región para incrementar la producción de carne. Estas limitantes presentan una jerarquía superior a otras más comúnmente intervenidas y su control permitirá evidenciar los beneficios de aplicar tecnologías enfocadas a la disponibilidad y la calidad forrajera, el manejo reproductivo o el control de enfermedades.
Conclusiones
Los resultados hallados indican que el sistema criador de la provincia de La Pampa debería ser replanteado si se desea incrementar sensiblemente su nivel productivo. Deberán considerarse la escasez de agua en la región y el aumento de las temperaturas estivales pronosticadas como limitantes a superar. La incorporación de rusticidad al ganado bovino ante las condiciones de cambio climático actual (y a futuro) debe ser considerado estratégicamente durante el proceso de reestructuración del sistema criador, similar a lo ocurrido en países con condiciones equivalentes.
Contribución de los autores
Pablo Mauricio Vázquez: definición de variables a analizar, construcción y análisis del sistema de información geográfico, análisis estadísticos, interpretación de datos y escritura. María del Carmen Rojas: análisis estadístico, construcción de las gráficas y figuras, escritura y revisión del artículo. Ezequiel Fernandez: adecuación y depuración de la base de datos histórica de existencias ganaderas provinciales y revisión del artículo, Ana Clara González: Adquisición de fotografías, definición de variables a analizar y revisión del artículo.
