INTRODUCCIÓN
El efecto del clima en la producción animal ha sido estudiado desde hace mucho tiempo, principalmente para el ganado bovino lechero, lográndose avances en el conocimiento de aspectos fisiológicos y de comportamiento animal bajo condiciones termoneutrales y de estrés climático. El cambio climático es una realidad y las estrategias de desarrollo agropecuario deben estar respaldadas por el profundo conocimiento de las características agrometeorológicas de las áreas correspondientes, lo cual es un factor determinante debido a la vinculación que el propio proceso productivo tiene con las mismas (Domínguez-Hurtado, Moya-Álvarez y Estrada-Moreno (2010).
En la ganadería vacuna, el bienestar de los animales se ha convertido en un factor determinante y el índice de temperatura-humedad (ITH) es un indicador de dicho bienestar muy usado para evaluar el nivel de confort de los animales (Carabaño et al., 2016; Nguyen et al., 2016). Estos estudios son escasos para las razas carniceras.
El objetivo de este trabajo fue caracterizar empresas genéticas del país que tienen razas carniceras a través de 5 variables climáticas y su evolución en el tiempo.
MATERIALES Y MÉTODOS
Localización, duración y variables climáticas
Se dispuso de la información mensual referente a temperatura ambiental mínima (Tmin), media (Tmed) y máxima (Tmax) siempre en 0C, la humedad relativa (%), de empresas y granjas genéticas pecuarias situadas en el occidente, centro y la región oriental del país. Las empresas y la ubicación de la estación meteorológica más cercana se muestran en la Tabla 1.
Empresa | Estación | Latitud | Longitud | Altura (m) |
---|---|---|---|---|
Camilo Cienfuegos | Pinar del Río | 220 24’16” | 830 39’14” | 56,48 |
San Juan de los Ramos | Jovellanos | 220 47’45” | 810 11’07” | 25,33 |
Rodas | Aguada | 220 22’22” | 800 49’35” | 28,34 |
Abra | Cienfuegos | 220 11’25” | 800 26’39” | 42,40 |
Rescate de Sanguily | Camagüey | 210 25’21” | 770 50’59” | 118,00 |
Manuel Fajardo | Manzanillo | 200 10’45” | 770 09’55” | 20,17 |
Turiguanó | Cayo Coco | 22º 32’21” | 780 22’07” | 3,40 |
Vallina | Valle de Caujerí | 20008’51” | 740 50’29 | 160,00 |
Los años incluidos fueron desde 1980 a 2018 y los doce meses del año y en total se utilizaron 3 774 registros.
Se elaboró la índice temperatura humedad (ITH) según lo planteado por Ravagnolo, Mistzal y Hoogenboom (2000).
ITH = 0,81*Ta + (HR/100) * (Ta - 14,4) + 46,4, donde Ta es la temperatura ambiente media en oC y HR la humedad relativa media en %.
Análisis estadísticos
Los datos se tabularon y se analizaron obteniéndose los estadígrafos generales, se realizaron análisis de varianza para las variables climáticas incluyendo como fuentes de variación empresa, quinquenio y cuatrimestre de inicio de la prueba de comportamiento (PC), así como las interacciones dobles y triples entre ellos.
Para el análisis estadístico se agruparon los meses en cuatrimestres naturales (enero-abril, mayo-agosto y septiembre-diciembre) y los años se agruparon por quinquenios (1980-1984, 1985-1989, 1990-1994, 1995-1999, 2000-2004, 2005-2009, 2010-2014 y 2015-2018). De esta forma las combinaciones generadas en el análisis de las interacciones fueron menores y permitieron interpretar mejor los resultados.
De las variables escogidas las temperaturas máximas y medias no presentaron distribución normal, por lo que para los análisis los datos fueron transformados por LOG(Y). El resto presentaron distribución normal. Se utilizó un modelo lineal general a través del PROC GLM del SAS (2013) con el siguiente modelo matemático:
Yijkl = µ + Ei + Aj + Mk + (ExA)ij + (ExM)ik + (AxM)jk + (ExAxM)ijk + eijkl
Donde Yijkl es la variable climática correspondiente, µ es la media general común a todas las observaciones; Ei el efecto de empresa (i = 1, 2, ….8); Aj el efecto de quinquenio de inicio de la PC (j = 1, 2, ….8) ; Mk es el cuatrimestre (k = 1, 2, 3) ; (ExA)ij, (ExM)ik, (AxM)jk y (ExAxM)ijk son las interacciones entre los efectos principales y eijkl es el efecto residual aleatorio N (0, σe2). Posteriormente las interacciones que no fueron significativas fueron eliminadas de los modelos y se hicieron nuevas corridas. Se obtuvieron las medias mínimas cuadráticas y se aplicó la dócima de Tukey para la comparación múltiple de medias.
Mediante el procedimiento PROC CORR del SAS (2013) se estimaron las correlaciones lineales de Pearson y además se estimaron las ecuaciones de regresión bajo determinadas circunstancias.
RESULTADOS Y DISCUSIÓN
Caracterización general
En la Tabla 2 se muestran los estadígrafos generales para las principales variables estudiadas. Los resultados obtenidos plantean coeficientes de variación relativamente bajos, presentándose las mayores variaciones para la Tmin y la Tmed. Suárez Tronco et al. (2018) en un análisis con menos empresas y registros encontraron resultados similares pero la variabilidad de los datos ahora fue menor quizás debido al mayor volumen de información. La temperatura mínima fue la más variable en ambos resultados. Independiente de que se analizaron datos de diferentes regiones del país las diferencias no son tan notables. Las condiciones climáticas resultan ser un elemento importante del agroecosistema y debe ser tenido en cuenta a la hora de escoger los genotipos y/o razas a explotar, así como el sistema de producción. Muchas veces se hace lo contrario, se escogen los genotipos independientes de las condiciones climáticas lo que no es aconsejable.
Variables | Medias | DS | CV (%) |
Temperatura mínima (0C) | 21,91 | 1,46 | 6,66 |
Temperatura media (0C) | 27,29 | 1,52 | 5,57 |
Temperatura máxima (0C) | 30,57 | 1,68 | 5,49 |
Humedad relativa (%) | 76,19 | 4,18 | 5,49 |
ITH | 77,78 | 3,45 | 4,43 |
DS desviación estándar; CV coeficiente de variación.
Existen diferentes clasificaciones del ITH para valorar el nivel de stress de los animales (Mader, 2003; de Rensis, García-Ispierto y López-Gatius, 2015); Enríquez Regalado y Álvarez Adán, 2020), la mayoría de las cuales han sido utilizadas para bovinos lecheros.
El ITH de 77,78 cae en la clasificación de stress leve para Mader (2003), stress severo para de Rensis, García-Ispierto y López-Gatius (2015), y supera el umbral señalado por Collier et al. (2012) para razas lecheras. En la provincia de Mayabeque, Enríquez Regalado y Álvarez Adán (2020) reportaron un ITH de 78 en el período 2005-2016, que consideraron stress moderado y un contraste marcado en su distribución anual. Recalcamos que las diferentes clasificaciones están aplicadas principalmente para zonas templadas y en razas de leche. Este indicador ha sido utilizado en algunos países para alertar a los productores sobre las condiciones de tensión de calor que amenazan el bienestar animal y evaluar el estado de confort del ganado (Vega, Almanza y Abraham, 2014). Como se ha dicho la mayor parte de los resultados se encuentran relacionados con la producción de leche, tales como: Ruíz-García, Carcelén-Caceres y Sandoval-Monzón (2018) y Ruíz, Carcelén y Sandoval-Monzón (2019) con la producción de leche diaria en Perú; Ruíz-Jaramillo et al. (2019) en Costa Rica y Enríquez Regalado y Álvarez Adán (2020) en la provincia de Mayabeque.
La temperatura media del aire de 27,29 °C, está dentro del rango reportado y calculada a partir de datos trihorarios, en el municipio de Guáimaro, en la provincia de Camagüey, en el periodo de 2000 a 2005, en explotaciones de ganado lechero con valores que oscilaron entre 22,9 °C y 27,7 °C. (Vega, Almanza y Abraham, 2014) y superior a la reportada por Suárez Tronco et al. (2018).
La humedad relativa (76,19%) resultó inferior a la reportada por Vega, Almanza y Abraham (2014) que presentó un valor medio de 82% y osciló durante el año entre 77%-86%. De igual forma resultó más baja que la reportada por Suárez Tronco et al. (2018) y similar al 76% en el área de investigación de la provincia de Mayabeque (Enríquez Regalado y Álvarez Adán, 2020).
La temperatura y la humedad relativa también han sido utilizadas en estudios sobre su influencia en los parámetros de la curva de lactancia en ganado Brown Swiss en condiciones tropicales (Lucena, 2014), así como algunas respuestas fisiológicas al calor y la humedad sobre diferentes genotipos (Espinoza et al., 2011).
Resultados de los análisis de varianza
Las interacciones triples de empresa x cuatrimestre x quinquenio no resultaron significativas en ningún caso, la interacción doble de cuatrimestre x quinquenio no resultó significativa para Tmed, HR, ITH y Tmax; por lo que fueron eliminadas de los modelos definitivos. La interacción más importante de acuerdo a nuestro objetivo empresa x quinquenio no resultó significativa para la Tmax.
Los resultados de los análisis de varianza para las variables climáticas analizadas aparecen en la Tabla 3. Todas las fuentes de variación analizadas en los modelos definitivos fueron altamente significativas, excepto la interacción empresa x quinquenio para la Tmax. Las variables climáticas variaron como era de esperar, en función de la empresa, el quinquenio y el cuatrimestre, así como las combinaciones entre algunos de ellos. Los coeficientes de determinación fueron relativamente altos para el ITH (76,02%) mientras que la HR fue bajo 51,25%. Todos los coeficientes de determinación fueron inferiores a los reportados por Suárez Tronco et al. (2018), quienes consideraron menos registros y empresas, pero igualmente dispersos por todo el país.
Fuentes de variación | GL | Tmin | Tmed | Tmáx | HR | ITH |
Empresa | 7 | *** | *** | *** | *** | *** |
Quinquenio | 7 | *** | *** | *** | *** | *** |
Cuatrimestre | 2 | *** | *** | *** | *** | *** |
Empresa x Cuatrimestre | 14 | *** | *** | *** | *** | *** |
Empresa x Quinquenio | 49 | *** | *** | NS | *** | *** |
R2 (%) | 71,78 | 74,57 | 63,63 | 51,25 | 76,02 |
*** P (≤ 0.001) NS no significativo
En la Fig.1 se presentan las variaciones para el ITH, la temperatura mínima y media a través de los cuatrimestres, y se aprecia un gran paralelismo entre las variables, independiente del efecto de escala. Se tomaron estas variables al presentar los coeficientes de determinación más altos y estar relacionadas positivamente entre sí (Ver Tabla 4). Hubo diferencias significativas (P<0,05) en todos los casos entre cuatrimestres.
El ITH se incrementa desde enero-agosto y luego comienza a declinar, presentándose los mayores valores entre los meses de julio (80,42) y agosto (80,61) que se pueden considerar como de peligro para los animales. Los valores más altos en estos meses también los encontraron Vega, Almanza y Abraham (2014) en Guáimaro, aunque fueron más bajos (76,9). Enríquez Regalado y Álvarez Adán (2020) en la provincia de Mayabeque encontró los valores más altos entre los meses de junio y septiembre (83,0) superiores a los de este trabajo. Los resultados anteriores corroboran que a fines y principios de año se presentan las mejores condiciones climáticas.
Dado que las correlaciones entre la Tmin y Tmed fue de 0,53 es lógico el paralelismo entre ellas (Fig. 1). Las temperaturas se incrementaron desde enero hasta agosto y luego disminuyeron. Las tendencias para las 3 variables son similares a las reportadas por Suárez Tronco et al. (2018). En aquella ocasión se reportaron las medias aritméticas y ahora son las medias mínimas cuadráticas ajustadas a los efectos estudiados.
Tmin | Tmed | Tmax | HR | ITH | |
Tmin | - | 0,53*** | 0,45*** | 0,27*** | 0,38*** |
Tmed | - | 0,55*** | 0,002NS | 0,68*** | |
Tmax | - | -0,14*** | 0,31*** | ||
HR | - | 0,16*** |
*** P (≤ 0.001) NS no significativo
Más importante que las variaciones entre épocas o cuatrimestres son las variaciones ocurridas a través del tiempo como clara señal del cambio climático. En la Fig. 2 se muestran las variaciones por quinquenios y a través de un análisis de regresión lineal se evidenció que la Tmin se ha incrementado en 0,29oC y la Tmed en 0,10oC por quinquenio, encontrándose coeficientes de determinación (R2) de 89,09% y 78,47%, respectivamente. En todos los casos el último quinquenio (2015-18) presentó los peores indicadores.
Caracterización de las empresas
Los estadígrafos generales para las empresas analizadas se presentan en la Tabla 5. Se puede apreciar que las diferencias entre empresas fueron en muchos casos notables. El ITH como medida integradora de la temperatura y la humedad, es un indicador muy importante y desde este punto de vista Turiguanó presentó los valores más altos que difirieron significativamente de los restantes, presentando Manuel Fajardo las mejores condiciones.
El ITH osciló entre 71,06 para Manuel Fajardo a 84,62 en Turiguanó. Valores de ITH < 70 se pueden considerar normales, sin afectaciones para los animales, entre 74 y 76 se califica el régimen como de alerta, con abatimiento insignificante de los animales, donde se encontraron solamente 3 empresas (Manuel Fajardo, Abra y Camilo Cienfuegos) encontrándose el resto por encima de estos valores que se considera como régimen de peligro con abatimiento significativo en los animales, de acuerdo a las clasificaciones anteriores, pero no fueron comprobadas las posibles alteraciones en los animales. En la Fig. 3 aparece el comportamiento por empresa para el ITH donde se ven más claramente las diferencias.
Empresa | HR (%) | T. max. (oC) | T. med. (oC) | T. min. (oC) | ITH |
Camilo Cienfuegos | 76,95 b | 28,07 e | 25,73 e | 23,10 b | 75,67 e |
San Juan | 74,57 e | 30,56 c | 27,72 b | 21,04 d | 78,79 c |
Rodas | 76,04 c | 30,10 c | 29,50 b | 20,85 e | 81,70 b |
Abra | 72,68 f | 31,64 b | 25,64 c | 21,98 c | 75,34 e |
Rescate de Sanguily | 75,09 d | 32,10 a | 26,40 d | 22,65 b | 76,80 d |
Manuel Fajardo | 78,03 b | 31,58 b | 24,99 f | 19,65 f | 71,06 f |
Turiguanó | 75,09 d | 32,10 a | 31,43 a | 22,65 b | 84,62 a |
Vallina | 81,09 a | 28,43 d | 26,88 c | 23,34 a | 78,30 c |
Letras distintas en la misma columna significa diferencias significativas (P≤ 0,05) prueba de Tukey
La empresa “Manuel Fajardo”, contrario a lo que debía esperarse debido a que se encuentra en la región oriental, presentó los valores más bajos de ITH (71,06) que difirió significativamente (P< 0,05) del resto, ligeramente inferior al reportado por Suárez Tronco et al. (2020) para la propia empresa que tiene ganado Criollo y Charolais, mientras que la peor empresa fue Turiguanó con ganado Santa Gertrudis. Todas las empresas excepto la Manuel Fajardo tuvieron valores promedio superiores a 74 que se pueden considerar en condiciones de alerta y 4 estuvieron por encima de 77 que son condiciones de stress.
En la Fig. 4 aparece la representación gráfica de la interacción empresa x quinquenio para el ITH. Dentro de cada empresa aparecen los puntos que representan las medias mínimas cuadráticas del ITH por quinquenio. Independientemente de que hay pequeñas variaciones en el orden de mérito de las empresas por quinquenio hay una tendencia bastante clara de que el ITH se ha incrementado en el tiempo en todas las empresas y situación similar se presenta con las temperaturas.
La actividad agropecuaria impone la necesidad de manejar eficientemente los recursos naturales, entre ellos el clima, y establecer estrategias que permitan la adaptación ante aquellos que resulten negativos para la producción. El bienestar de los animales se ha convertido en un factor determinante para lograr una mejor expresión productiva.
Los ambientes tropicales ofrecen ventajas y desventajas para la producción animal; dentro de los aspectos positivos se pueden mencionar que la precipitación alta y la duración solar diaria casi constante, favorecen la producción de forraje durante gran parte del año. En cambio, las desventajas están asociadas a: la proliferación de parásitos internos y externos; las enfermedades; los contenidos altos de pared celular en los forrajes que favorecen la producción de calor interno; y las temperaturas y HR altas, que conducen a una situación de estrés calórico afectando la expresión del potencial genético al comprometerse su bienestar (Tapki y Şahin, 2006). Entre los eventos fisiológicos desencadenados por el estrés calórico se mencionan el aumento de la temperatura rectal, la frecuencia respiratoria y el jadeo, que se manifiestan para buscar mantener la temperatura corporal (Pragna et al., 2017). Estos cambios conllevan una alteración en el patrón de alimentación y función del rumen, con reducción de la ingesta de materia seca y, por consiguiente, de la productividad.
CONCLUSIONES
Estos resultados permiten caracterizar y/o clasificar agrometeorológicamente nuestras empresas desde el punto de vista climático. Esto facilita valorar y escoger las razas y/o genotipos adecuados a las condiciones ambientales imperantes a tono con la regionalización de la explotación según razas lecheras o de carne y sus cruzamientos. Igualmente se observa un empeoramiento de las variables climáticas en el tiempo afectando la explotación ganadera.