Introducción

La gran expansión del ganado charolais se debe a su capacidad de adaptación, precocidad, rusticidad, robustez y propensión a ganar peso en todas las edades. Se encuentra en regiones de clima tropical, subtropical, templado y árido, lugares donde se logra su aclimatación (INRA, 2011). Uno de estos lugares es Cuba donde, en la actualidad, la hembra en desarrollo se incorpora a la reproducción en edades superiores a los 20 meses, ya que no alcanza ganancia de peso por encima de 400 g·animal·día^-1 (Vega et al., 2013), por lo que a esa edad, de acuerdo con los objetivos genéticos en Cuba, está por debajo de 300 kg de peso vivo necesarios para la incorporación a la reproducción debido, fundamentalmente, a los sistemas de crianza utilizados (Estadística Departamento de Genética Nacional, 2016).

El sistema de manejo y alimentación determina la velocidad de crecimiento (GMD) y la relación peso vivo por edad (Carriquiry et al., 2016), por lo que el estudio de estos indicadores debe garantizar la selección de las variantes tecnológicas efectivas de ceba que mejor satisfagan la demanda de hembras para la reproducción. En la empresa Genética y Cría “Manuel Fajardo”, en Granma, Cuba, la raza charolais no está exenta de lo planteado con anterioridad y esto determina que su incorporación a la reproducción no se realice antes de los 20 meses. Por lo que al analizar el sistema productivo se hizo ineludible la valoración del balance alimentario al incluir los árboles leguminosos, debido a sus ventajas de proporcionar alimentos nutritivos, con valores superiores de proteína, lo anterior favorecería las ganancias de peso y mejoraría la relación peso por edad a la incorporación de las hembras charolais en desarrollo (Vega et al., 2016). Además, proveen sombra, madera, conservan el suelo y el agua, ayudan a capturar carbono y nitrógeno del aire, entre otros factores (Milera, 2013).

Por ello, el objetivo de este trabajo fue valorar elementos nutricionales que definen la edad de incorporación a la reproducción de vaquillas charolais a sistemas de silvopastoreo y compararlos con el monocultivo de Cynodon nlemfuensis.

Materiales y métodos

En la finca San José del Retiro, Jiguaní, Granma, Cuba, que pertenece a la empresa Genética y Cría “Manuel Fajardo”, se realizó la investigación durante tres años, en un suelo pardo con carbonato (Hernández et al., 1999), ubicado a 115 m sobre el nivel del mar, con relieve ligeramente llano.

El comportamiento de las variables climáticas fue: temperatura media anual de 25.4 °C, con media de 23.9 °C y 26.8 °C en el período poco lluvioso y lluvioso, respectivamente. La precipitación media anual fue de 939 mm y la humedad relativa fue de 81.7%.

En un área total de 26.40 ha se establecieron dos tratamientos: silvopastoreo de Leucaena leucocephala-Cynodon nlemfuensis y monocultivo de Cynodon nlemfuensis. Se distribuyeron 12 ha para cada tratamiento en 24 cuartones de similar tamaño. Las restantes 2.4 ha se destinaron para las corraletas donde permanecían los animales durante la noche y para las mangas de acceso a los cuartones.

Para el establecimiento del silvopastoreo se realizó inicialmente la composición botánica del área donde las especies presentes fueron: 25% de estrella (C. nlemfuensis), 25% de romerillo (Bidens alba Lin.), 3% de Guinea (Megagthrysus maximus) y el resto era área despoblada. Después el suelo se preparó de forma convencional mediante arado, dos pases de grada y surcado con bueyes espaciando los surcos 5 m entre sí. La leucaena se sembró a una distancia de 3 m entre plantas con semilla de 85% de germinación y se logró una densidad promedio de 500 plantas ha^-1. Durante el establecimiento las malas hierbas se eliminaron manualmente y se consideró que la leucaena estaba establecida cuando alcanzó 2 m de altura a los 15 meses.

El monocultivo (C. nlemfuensis) se sembró a vuelta de arado y se consideró establecido a los 100 días de sembrado. Durante el período experimental se hizo una poda de la leucaena debido a la altura alcanzada (2.6 m) y coincidió con la entrada de los animales en el ciclo II. Además, se observó la presencia de Clitoria ternatea en 2%.

Se emplearon 30 novillas charolais de Cuba, 15 por tratamiento, repetido en dos ciclos (I y II) productivos para un total de 60 animales en un tiempo de tres años. Estos se agruparon en lotes homogéneos según el peso (168 kg en el ciclo I y 155 kg en el ciclo II) y edad de ocho meses, se desparasitaron antes del período experimental con 1 mL de Levamisol por cada 100 kg de peso vivo. Tuvieron libre acceso a la sal mineral compuesta por fosfato dicálcico 50.0%, 40.0% de sal común y 10.0% de una mezcla integrada por: maíz 19.0%, óxido de hierro 27.5%, sulfato de manganeso 23.0%, óxido de cobre 10.0%, óxido de zinc 20.0%, selenito de sodio 0.02% y sulfato de cobalto 0.1%. Los animales pastaron 14 horas diarias y se mantuvieron en las corraletas durante la noche con agua, pero sin alimentos. Se utilizó una carga de 1.1 animales·ha^-1. Los tiempos de reposo fueron 77 y 110 días por cuartón en el período lluvioso y poco lluvioso, respectivamente. Se cuantificaron los siguientes indicadores: peso vivo mediante el pesaje, en pesa de 500 kg, cada 28 días de todos los animales, en ayuno, para calcular la GMD entre pesadas y la ganancia acumulada del período mediante diferencia de peso inicial y final de cada ciclo según Preston y Willis (1970) y Menéndez (1984).

La disponibilidad de pasto se estimó según Martínez et al. (1990) mediante la toma de 80 muestras por réplica al inicio de cada rotación y paralelamente a los muestreos de disponibilidad, pero con una frecuencia bimestral, se tomaron muestras de pastos a 10 cm por encima del nivel del suelo (300 g) simulando con la mano la selección que hace el animal en pastoreo para determinar la PB (proteína bruta), FC, extracto etéreo (EE) y extracto libre de nitrógeno (ELN) de acuerdo con AOAC (1995).

La energía metabolizable (EM) se estimó mediante la ecuación de Ibrahim (1988) para calcular el total de nutrientes digestibles (TDN) según:

TDN = 37.937- 1.018FC – 4.886EE + 0.173ELN + 1.042PB + 0.015FB. – 0.058EE² + 0.008FB·ELN +0.038 EE·PB +0.003EE²·PB

Y

1kg TDN = 3.62 Mcal de EM.

Con posterioridad la Mcal (megacaloría) se convierten en MJ (megajoules).

La disponibilidad de L. leucocephala vc. Perú se cuantificó en 10 árboles al azar en cada réplica, simulando el ramoneo que realizan los animales. Se aplicó la técnica del ordeño de las partes más tiernas de las plantas, las hojas y los tallos tiernos hasta aproximadamente 2 mm de diámetro y se tomó una muestra cada dos meses para determinar la composición bromatológica (Lamela, 1998).

Para realizar el balance alimentario se estimó el consumo de MS de los animales por medio de las tablas presentadas por Martín (1981) y de los requerimientos de Martín y Palma (1999).

El análisis estadístico se realizó mediante la metodología de determinación del índice de impacto desarrollada por Torres et al. (2007). En la selección de los componentes principales se tuvo en cuenta que el valor propio fuera mayor que la unidad, se asumió que las variables tuvieran factor de preponderancia mayor a 0.60 y de acuerdo con las características de cada componente se les asignó una denominación (Torres et al., 2008). A partir del análisis de los gastos y el monto total de la inversión por tratamiento, se calcularon los principales indicadores económicos: costo·animal^-1, costo·kg^-1 de peso vivo, costo·peso cubano, costo·ha^-1 y la relación beneficio/costo. Se tomó como escala un rebaño de 120 animales, simulando un rebaño productivo, con carga de 1.10 animal·ha^-1 y los datos primarios puntualizados anteriormente en ambos ciclos.

Resultados

El análisis de componente principal evidenció la presencia de cuatro componentes con valores propios mayores que la unidad y de conjunto contribuyeron a explicar 82.22% de la variabilidad entre los sistemas (cuadro 1). Las variables con las cargas factoriales de importancia fueron: en la CP1 (nombrada animal) las ganancias por animal y el peso; la CP2 (llamada alimento animal) mostró a la disponibilidad y edad de incorporación a la reproducción de la hembra Charolais; la CP3 (designada clima) señaló a la temperatura y la precipitación, y la CP4 (denominada edad de inicio) mostró sólo la propia edad de inicio.

Al tener en cuenta los resultados anteriores se realizó un análisis de conglomerado jerárquicos (cuadro 2) que permitió conformar grupos similares de animales a través del método de vinculación inter-grupo y la medida de intervalo de distancia euclidiana al cuadrado. Se obtuvieron tres grupos conformados por: los animales en silvopastoreo (grupo 1), aquellos que perdieron peso (grupo 2) y los que estuvieron en el monocultivo (grupo 3).

Los animales se evaluaron durante tres años en dos ciclos productivos, siendo superior el tiempo para alcanzar los 300 kg de peso vivo para su incorporación a la reproducción en el sistema monocultivo, lo cual se puede apreciar en las curvas de figura 1; las primeras corresponden a la tecnología silvopastoril y seguidamente la del monocultivo.

Figura 1
Índice de impacto animal en cada ciclo de ambos tratamientos con el peso y la ganancia.

Los indicadores en el ciclo II siguieron el mismo patrón de respuesta que en el ciclo I donde sobresalió el silvopastoreo aunque se encontraron valores superiores de eficiencia en la PB y EM en ambos sistemas (cuadro 4).

El análisis económico de ambos sistemas aparece en el cuadro 5. Todos los indicadores mostraron las ventajas económicas del sistema silvopastoril sobre todo en el costo por kg de peso vivo y por ha. Además, estos indicadores fueron superiores en el ciclo II comparado con el ciclo I, en especial en el silvopastoreo.

En la CP1 (nombrada animal) las ganancias por animal y el peso son las variables con cargas factoriales de importancia, lo que indicó la magnitud que ellas ocupan en tecnologías relacionadas con la incorporación de las hembras de reemplazo charolaise y contribuyeron a explicar 33.83% de la variabilidad de los sistemas.

La CP2 (llamada alimento-edad), mostró la relación positiva de la disponibilidad e inversa de la edad de incorporación de la hembra charolais; esto explicó 18.72% de la variabilidad. Con estos resultados se demostró la estrecha relación entre la disponibilidad de alimentos, incluyendo la calidad, o sea los aportes para cubrir los requerimientos de las hembras, y su edad a la incorporación a la reproducción. Esto pudiera reflejar que a mayor disponibilidad disminuiría la edad de incorporación y viceversa.

En la CP3 (designada clima) se destacaron la precipitación y la temperatura como las de mayor preponderancia en una relación directa y que explicaron 15.81 % de la variabilidad. Lo anterior está relacionado con el efecto de dichas variables en la disponibilidad de alimentos y con el bienestar animal.

La CP4 (denominada edad de inicio), mostró sólo la propia edad de inicio que explicó 13.86% de la variabilidad vinculada a los procesos de crecimiento de los animales que, a su vez, están determinados por el comportamiento de las variables que integran las componentes antes descritas.

Al tener en cuenta las cargas factoriales altas en los ciclos, las variables relacionadas con la producción animal, el alimento, clima y la edad inicial, fueron elementos importantes en la medición del impacto, debido a que las cargas factoriales expresan las correlaciones entre cada variable y el componente principal (cuadro 1). Por lo tanto, indican el grado de correspondencia entre ellas, haciendo a una variable con mayor carga representativa del componente principal. Se evidenció que la mayoría de dichas variables mencionadas aportan a la explicación de la variabilidad de los ciclos de los animales estudiados y que la rotación varimax logró soluciones factoriales más simples y eliminó alguna de las ambigüedades que acompañan a las soluciones factoriales inicialmente rotadas, por lo que se determinaron cuatro componentes principales que explican 82.22% de la variabilidad.

Al tener en cuenta los resultados anteriores se realizó un análisis de conglomerado jerárquicos (cuadro 2) que permitió conformar grupos similares de animales a través del método de vinculación inter-grupo y la medida de intervalo distancia euclidiana al cuadrado. Se obtuvieron tres grupos conformados por: los animales en el silvopastoreo (grupo 1), aquellos que perdieron peso (grupo 2) y el tercer grupo corresponde al monocultivo (grupo 3).

El grupo 1 (sistema silvopastoril) contó con la mayor disponibilidad de biomasa por la inclusión de la arbórea y mejoró la calidad del alimento cosechado por los animales, lo que propició las mayores ganancias de peso (0.55 kg·animal·día^-1) que significó 16.72 kg·animal·mes^-1 y determinó la menor edad de comienzo a la vida reproductiva (19.93 meses).

El grupo 2 registró pérdidas de peso con el menor peso vivo de los tres grupos, y aun cuando la disponibilidad fue de 24.7 kg de MS, existió déficit en el suministro de agua y no permitió mejorar la edad de incorporación a la reproducción de las hembras de reemplazo.

El grupo 3 (monocultivo de pasto estrella) alcanzó el mayor peso con menor disponibilidad de pasto, pero la edad de incorporación fue de 25 meses con 9.22 kg·animal·mes^-1 de ganancia y 0.31 kg·animal·día^-1, lo que propició que la edad para comenzar la reproducción se retardara, debido a que los animales no pudieron cubrir sus requerimientos nutricionales.

La comparación de los valores obtenidos para cada variable identificó, claramente, la situación particular de cada grupo de animales, así como la repercusión de la introducción de la tecnología de los árboles leguminosos en las áreas de pastoreo. El grupo 1 alcanzó la mayor disponibilidad de alimento (29.8 kg MS·animal^-1) y calidad debido al efecto de asociar la leucaena al pasto estrella, lo que propició las mayores ganancias comparado con el monocultivo de C. nlemfuensis representado por el grupo 3 con una oferta de pasto de 23.0 kg, 9.22 kg·animal·mes^-1 y 0.31 kg·animal·día^-1, respectivamente. El grupo 2 presentó los peores valores lo cual estuvo motivado, principalmente, por las dificultades presentadas con el abasto de agua y es conocido el importante papel que desempeña el adecuado balance hídrico del animal en la eficiente utilización de los nutrientes disponibles en la dieta consumida.

Las evaluaciones mostraron que con la carga utilizada (1.1 animales·ha^-1) y sal mineral a voluntad sólo se alcanzó por encima 500 g·animal^-1 de GMD de peso vivo (grupo 1) en la tecnología silvopastoril (figura 1), aunque al inicio del ciclo I por estrés hídrico transitorio sólo se llegó a 69 g·animal^-1 de GMD , pero al ser eliminado el referido estrés se logró alcanzar la edad y peso a la incorporación para la reproducción, con hembras charolais de reemplazo y esto correspondió con la mayor disponibilidad. Sin embargo, el monocultivo de C. nlemfuensis (grupo 3) no superó los 310 g·animal^-1 de GMD en ambos ciclos, lo que evidenció la superioridad de la introducción de los árboles leguminosos en los sistemas de producción de hembras de reemplazo.

Se logró demostrar que es posible incorporar los animales a la reproducción con más de 300 kg de PV entre 21 y 17 meses de edad, en 273 y 366 días de permanencia en el sistema silvopastoril y GMD total de 416 y 567 g·animal^-1 en cada ciclo correspondiente. Este resultado es relevante ya que se obtuvo sin gasto de suplemento·animal^-1, en comparación con los que pastorearon en el C. nlemfuensis donde la edad de incorporación estuvo entre 25 y 28 meses, en 578 y 548 días de permanencia y una GMD total inferior de 241 y 278 g·animal^-1 en cada ciclo.

Esto significó que los animales entran a los sistemas (silvopastoril y monocultivo) con 155 a 170 kg y siete a nueve meses de edad, se incorporan a la reproducción con 17 meses de edad y más de 300 kg de PV solo en el silvopastoreo (figura 1), aunque en el primer ciclo se obtuvieron valores por debajo de cero debido al estrés hídrico mencionado anteriormente; los valores positivos demostraron un comportamiento que no se había reportado para la hembra de reemplazo en pastos tropicales en áreas de Latinoamérica al considerar el pasto (sin riego ni fertilizantes) como único alimento. Lo anterior permite potenciar el incremento de la GMD, desde el nacimiento hasta el destete, en las fincas de ganado de cría con ganado charolais de Cuba, en pastoreo asociado de gramíneas y leguminosas.

El uso de Cynodon en monocultivo no produjo alta GMD, tercera y cuarta curvas (figura 1), pero su evaluación fue necesaria, porque representa una estrategia, para comparar la oferta de MS, dado que, en el período poco lluvioso, en esta zona climática, los animales pierden peso a causa de la baja disponibilidad de alimento voluminoso y falta de suplementos. Además, este trabajo constituyó un aporte en la evaluación del silvopastoreo, con el mejor ganado de carne de Cuba, con pasto y sal mineral (Díaz, 2008). Iglesias (2003) reportó GMD acumuladas similares a las del presente trabajo, de 524 y 441 g·animal^-1, pero con F1 (1/2 Holstein x 1/2 Cebú) y 5/8 H x 3/8 C, donde el sistema de alimentación estaba basado en pastoreo restringido de mezclas de pastos naturales con estrella, concentrado comercial y forraje de caña con urea a voluntad, aunque con hembras charolais en desarrollo reportaron ganancias similares (Vega et al., 2016).

La GMD, el crecimiento de los animales, la conversión de la EM y la PB se pueden relacionar con la suficiente disponibilidad de pasto en el silvopastoreo que se expresó en función de la presión de pastoreo (26.40 y 18.40 kg de MS·100 kg PV^-1), la composición proteico-energética del pasto y la leucaena permitió mejor ajuste de la relación N total: EM y la conversión de la PB en GMD (cuadros 1 y 2). Esto se debe, según Preston y Willis (1970), al menor gasto energético de mantenimiento de este genotipo, porque es de maduración tardía. Similar resultado obtuvo Fernández (2017) al estudiar los sistemas silvopastoriles.

Las superiores GMD de los animales del silvopastoreo, con respecto a los del pastoreo Cynodon en monocultivo se puede explicar por la mayor población de bacterias celulolíticas y viables totales que aparecen cuando los animales consumen leucaena según lo reportado para esta misma raza en similares condiciones de alimentación (Díaz, 2008). Galindo et al. (2014) determinaron que el aumento de estas poblaciones microbianas permite mejor degradación de la MS del pasto y la asimilación de nutrientes.

La nutrición desempeña un papel muy importante en el comportamiento de la reproducción del vacuno, lo que es, tal vez, uno de los factores que más contribuyen a la expresión del potencial genético de los animales e incluso de su progreso (Vega, 2012).

En el silvopastoreo en ambos ciclos el consumo de materia seca del pasto, (cuadro 3 y 4), mantuvo una presión de pastoreo (kg MS·100 kg^-1 de peso vivo) superior a 5 kg, este suministro garantizó la selección del material a consumir y en consecuencia, mejor comportamiento animal como reflejo de la presencia de la leguminosa y de la contribución de ésta al sistema.

La disponibilidad de materia seca y los nutrientes de los sistemas de pastoreo no fue igual si se considera por épocas dentro de los años. El período poco lluvioso mostró ser el más crítico, con un comportamiento estacional marcado en el caso del Cynodon en monocultivo, aunque logró cubrir las recomendaciones, según Martín (1981), para la oferta de materia seca y en el caso del ciclo II (cuadro 4) hubo un déficit de -15 g de proteína en junio. Sin embargo, fue el déficit de agua la principal limitante en dos meses del ciclo I para ambos sistemas (figura 1) alcanzando valores menores de cero.

De hecho, el peso es más importante que la edad en la manifestación y posterior mantenimiento de celos. Lo anterior resalta, en consecuencia, la importancia de la alimentación para mantener ritmos de crecimiento adecuado para lograr la manifestación del ciclo estral a partir de los dos años en todo tipo de ganado.

Los resultados obtenidos en las hembras charolais con solo pasto estrella comparadas con el sistema silvopastoril, permitió demostrar en este último la capacidad de hacer más eficiente la conversión alimenticia en períodos de restricción y de abundancia continuos, el efecto en la ganancia de peso y el inicio temprano de la pubertad, el cual se reflejó en mejores indicadores productivos y reproductivos del rebaño en general.

Se deben aprovechar las ventajas de la charolais, que mantenidas en pasturas de regular calidad pueden igualar o superar la ganancia si tienen acceso a pastos de mejor calidad (Unduganado, 2015), ya que con su alto poder de conversión alimenticia (1 kg de peso vivo x 7.26 kg de alimento), le permite alta velocidad de crecimiento (hasta 2.5 kg·día^-1) por su buena capacidad de adaptación a diferentes condiciones de cría y elevada ganancia de peso con forrajes bastos como los resultados del grupo 1 (silvopastoreo), aunque el grupo 2 debido a la falta de agua presentó pérdidas de peso de 0.71 kg·animal·dia^-1 (cuadro 2).

En el caso estudiado de la charolais, el costo·peso^-1 se comportó inferior en el silvopastoreo y no superó los 0.42 centavos por cada peso y en monocultivo los 0.52. El indicador beneficio/costo en el sistema silvopastoril supera en 26% al monocultivo (23 y 28 para los ciclos I y II respectivamente (cuadro 5). Este análisis permite recomendar la tecnología silvopastoril garantizando que, si se cumple con la disciplina tecnológica, el sistema es sostenible desde el punto de vista económico.

Conclusiones

Las valoraciones de elementos nutricionales como la PB y EM durante los ciclos productivos de investigación argumentaron la superioridad del sistema silvopastoril (Leucaena-pasto estrella) con respecto al monocultivo de pasto estrella, expresado por la menor edad de incorporación, mayor peso en vaquillas charolais, no empleo de riego, fertilización, ni suplemento y mejor relación costo/beneficio.

Recomendaciones

Extender este sistema a otras unidades de producción comercial y desarrollar similar investigación en otras condiciones ambientales y razas vacunas.

Agradecimientos

A la empresa de Genética y Cría “Manuel Fajardo” en el municipio Jiguaní, Granma, Cuba y a sus trabajadores, por las facilidades brindadas para la realización de esta investigación.

Literatura citada

AOAC (1995). Official methods of analysis. Ass. Off. Agric. Chem. 16^th ed. Washington, DC.

Carriquiry, E.J.; Fernández, J.A. y Yáñez, E.A. (2016). Influencia del peso vivo y la alimentación sobre el desempeño de terneros en su primer invierno. Rev. Vet. 27(1): 14-20.

Díaz, A. (2008). Producción de carne bovina en pastoreo con gramíneas y leguminosas. Tesis presentada en opción al grado científico de Doctor en Ciencias Veterinarias. Instituto de Ciencia Animal. MES. La Habana, Cuba.

Estadística Departamento de Genética Nacional (2016). Ministerio de la Agricultura. La Habana. Cuba. 100 p.

Fernández, A. (2017). Producción de carne y leche bovina en sistemas silvopastoriles. Ediciones INTA. Argentina. 195 p.

Galindo, J.; González, N.; Iraola, J.; Muñoz, E.; Aldana. A.I.; Moreira, O.; Hernández, J.L. y Moreira, E. (2014). Efecto de un sistema silvopastoril con Leucaena leucocephala vc Perú y la época de año en la población microbiana ruminal de toros mestizos de Cebú en crecimiento-ceba. III Convención Internacional Agrodesarrollo 2014. EEPF Indio Hatuey. ISBN 978-959-7138-18-1. Matanzas, Cuba. Pp. 1053-1057.

Hernández, A.; Pérez, J.M. y Boch, O. (1999). Nueva versión de la clasificación genética de los suelos de Cuba. Instituto de suelos. AGROINFOR-MINAG. Cuba. 64 p.

Ibrahim, M.M.M. (1988). Feeding tables for ruminants in Sri Lanka. Animal Feed Act No. 15 1st edition. Costa Rica. 156 p.

Iglesias, J. (2003). Los Sistemas Silvopastoriles, una alternativa para la crianza de bovinos jóvenes en condiciones de bajos insumos. Tesis presentada en opción al grado científico de Doctor en Ciencias Veterinarias. Instituto de Ciencia Animal - Universidad Agraria de La Habana. La Habana, Cuba.

INRA. (2011). Origen de la raza Charolaise (Multipropósito). Disponible en: http://hdavidgarciam.blogspot.com/2011/04/origen-de-la-raza-charolaise.html. (Consultado en 25 enero de 2013).

Lamela, L. (1998). Técnica de muestreo. Conferencia del curso de posgrado “Manejo de los pastos y forrajes para la producción animal”. Maestría en Pastos y Forrajes. Estación Experimental de Pastos y Forrajes “Indio Hatuey”. Matanzas, Cuba. 12 p.

Martín, P.C. (1981). Metodología de balance alimentario y formulario de raciones para el ganado bovino en Nicaragua. Ed. Ministerio de Desarrollo Agropecuario y Reforma Agraria. 177 p.

Martín, P.C. y Palma, J.M. (1999). Manual Para Fincas y Ranchos Ganaderos. Colima, México. Ed. AgroSystems.120 p.

Martínez, J.; Milera, M.; Remy, V.; Yepes, I. y Hernández, J. (1990). Un método ágil para estimar la disponibilidad de pasto en una vaquería comercial. Pastos y Forrajes 14: 101-109.

Menéndez, A. (1984). Un método simple para evaluar el crecimiento y desarrollo de nuestras hembras lecheras. ACPA 3: 13-16.

Milera, M. (2013). Contribución de los sistemas silvopastoriles en la producción y el medio ambiente. Avances en Investigación Agropecuaria 17(3):7-24.

Unduganado, A. (2015). Historia y origen de la raza Charolaise. Publicada en Unduganado com. El portal de Honduras para ganaderos. Disponible en: http://www.honduganado.com/historia-y-origen-de-la-raza-charolais/ (Consultado el 14 de noviembre 2017).

Preston, T. y Willis, M. (1970). Producción intensiva de carne. B. SC. (Dunelm). PH. (EDIN). División de Ciencia Animal. ICA. Habana. Cuba. Pp. 176-507.

Torres, V.; Benítez, D.; Lázaro, D. y Álvarez, A. (2007). Modelo estadístico para la medición del impacto de la innovación o transferencia tecnológica en la rama agropecuaria. XI Conferencia Española y Primer Encuentro Iberoamericano de Biomatemática. Universidad de Salamanca. España. Pp. 25-32.

Torres, V.; Ramos, N.; Lizazo, D.; Monteagudo, F. y Noda, A. (2008). Modelo estadístico para la medición del impacto de la innovación o transferencia tecnológica en la rama agropecuaria. Revista Cubana de Ciencia Agrícola 42(2): 133-139.

Vega, A. (2012). Crianza de hembras de reemplazo Charolaise de Cuba, en silvopastoreo con Leucaena leucocephala y Cynodon nlemfuensis. Tesis presentada en opción al grado de Doctor en Ciencias Veterinarias. Instituto de Ciencia Animal. La Habana, Cuba.

Vega, A.; Herrera, R.S.; Lamela, L.; Santana, A. y Rodríguez, G.A. (2013). Evaluación en silvopastoreo y un monocultivo de hembras en crecimiento Charolaise en el oriente cubano. REDVET. Revista Electrónica de Veterinaria 14(2).

Vega, A.; Herrera, R.S.; Torres, V.; Lamela, L.; Montejo, I.; Santana, A. y Cino, D. (2016). Comportamiento de las hembras de reemplazo Charolais de Cuba en un sistema silvopastoril comparada con un monocultivo. Revista Cubana de Ciencia Agrícola 50 (1): 51-59.