https://doi.org/10.29312/remexca.v15i6.3462

elocation-id: e3462

Almario, Mejia Salazar, and Meneses Buitrago: Production of forage leguminous species at different cutting times in Tolima, Colombia

Journal Metadata

Journal Identifier: remexca [journal-id-type=publisher-id]

Journal Title Group

Journal Title (Full): Revista mexicana de ciencias agrícolas

Abbreviated Journal Title: Rev. Mex. Cienc. Agríc [abbrev-type=publisher]

ISSN: 2007-0934 [pub-type=ppub]

Publisher

Publisher’s Name: Instituto Nacional de Investigaciones Forestales, Agrícolas y Pecuarias

Article Metadata

Article Identifier: 10.29312/remexca.v15i6.3462 [pub-id-type=doi]

Article Grouping Data

Subject Group [subj-group-type=heading]

Subject Grouping Name: Articles

Title Group

Article Title: Production of forage leguminous species at different cutting times in Tolima, Colombia

Contributor Group

Contributor [contrib-type=author]

Name of Person [name-style=western]

Surname: Almario

Given (First) Names: Nelson Perez

X (cross) Reference [ref-type=aff; rid=aff1]

Superscript: 1

Contributor [contrib-type=author]

Name of Person [name-style=western]

Surname: Mejia Salazar

Given (First) Names: Julian R.

X (cross) Reference [ref-type=aff; rid=aff1]

Superscript: 1

X (cross) Reference [ref-type=corresp; rid=c1]

Superscript: §

Contributor [contrib-type=author]

Name of Person [name-style=western]

Surname: Meneses Buitrago

Given (First) Names: Diego H.

X (cross) Reference [ref-type=aff; rid=aff2]

Superscript: 2

Affiliation [id=aff1]

Label (of an Equation, Figure, Reference, etc.): 1

Institution Name: in an Address: Centro de Investigación Nataima-Corporación Colombiana de Investigación Agropecuaria. Vía Espinal km 9, Chicoral, Tolima, Colombia. 730520-(1)4227300. (nperez@agrosavia.co). [content-type=original]

Institution Name: in an Address: Corporación Colombiana de Investigación Agropecuaria [content-type=normalized]

Institution Name: in an Address: Centro de Investigación Nataima- [content-type=orgdiv1]

Institution Name: in an Address: Corporación Colombiana de Investigación Agropecuaria [content-type=orgname]

Country: in an Address: Colombia [country=CO]

Email Address: nperez@agrosavia.co

Affiliation [id=aff2]

Label (of an Equation, Figure, Reference, etc.): 2

Institution Name: in an Address: Centro de Investigación Obonuco-Corporación Colombiana de Investigación Agropecuaria. Vía Pasto km 5. Obonuco, Colombia. 520001-(1)4227300. (dmeneses@agrosavia.co). [content-type=original]

Institution Name: in an Address: Corporación Colombiana de Investigación Agropecuaria [content-type=normalized]

Institution Name: in an Address: Centro de Investigación Obonuco [content-type=orgdiv1]

Institution Name: in an Address: Corporación Colombiana de Investigación Agropecuaria [content-type=orgname]

Country: in an Address: Colombia [country=CO]

Email Address: dmeneses@agrosavia.co

Author Note Group

Correspondence Information: [§] Autor para correspondencia: jmejia@agrosavia.co [id=c1]

Publication Date [date-type=pub; publication-format=electronic]

Day: 11

Month: 09

Year: 2024

Publication Date [date-type=collection; publication-format=electronic]

Season: Aug-Sep

Year: 2024

Volume Number: 15

Issue Number: 5

Electronic Location Identifier: e3462

History: Document History

Date [date-type=received]

Day: 01

Month: 01

Year: 2024

Date [date-type=accepted]

Day: 01

Month: 03

Year: 2024

Permissions

License Information [license-type=open-access; xlink:href=https://creativecommons.org/licenses/by-nc/4.0/; xml:lang=es]

Este es un artículo publicado en acceso abierto bajo una licencia Creative Commons

Abstract

Title: Abstract

Different studies have shown that leguminous forage species have a high impact and value in silvopastoral systems in terms of the contribution they represent for livestock feeding, based on the increase in both the potential of soils and the production of meat for livestock. The objective was to determine the productive potential in terms of dry matter in two contrasting seasons, dry and wet, for the leguminous species Albizia guachapele (Kunth) Dugand, Albizia niopoides (Benth.) Burkart, Albizia saman (Jacq.) Merr., and Leucaena leucocephala (Lam.) de Wit. This study was conducted at the Nataima Research Center of the Colombian Agricultural Research Corporation, located in the municipality of Espinal-Tolima, Colombia, between March and August 2018. A design of repeated measures over time and arrangement of randomized complete blocks, four treatments, and two evaluation periods were used, and eight interactions between the species and rainy season factors were worked on. The species A. guachapele presented the best productive potential compared to the others; there were statistical differences for cutting times in the species evaluated, and a high variability effect was observed due to the time factor. A. saman was the least productive species, while production stabilization was demonstrated for all those evaluated over time. There was no significant bidirectional interaction for the season factor, and the best performance corresponded to those carried out in the dry season. This allowed a greater expression of the productive potential as a food alternative for low rainfall areas.

Keyword Group [xml:lang=en]

Title: Keywords:

Keyword: adaptation

Keyword: interaction

Keyword: potential

Keyword: stabilization

Counts

Figure Count [count=4]

Table Count [count=1]

Equation Count [count=0]

Reference Count [count=24]

Page Count [count=0]

Abstract

Different studies have shown that leguminous forage species have a high impact and value in silvopastoral systems in terms of the contribution they represent for livestock feeding, based on the increase in both the potential of soils and the production of meat for livestock. The objective was to determine the productive potential in terms of dry matter in two contrasting seasons, dry and wet, for the leguminous species Albizia guachapele (Kunth) Dugand, Albizia niopoides (Benth.) Burkart, Albizia saman (Jacq.) Merr., and Leucaena leucocephala (Lam.) de Wit. This study was conducted at the Nataima Research Center of the Colombian Agricultural Research Corporation, located in the municipality of Espinal-Tolima, Colombia, between March and August 2018. A design of repeated measures over time and arrangement of randomized complete blocks, four treatments, and two evaluation periods were used, and eight interactions between the species and rainy season factors were worked on. The species A. guachapele presented the best productive potential compared to the others; there were statistical differences for cutting times in the species evaluated, and a high variability effect was observed due to the time factor. A. saman was the least productive species, while production stabilization was demonstrated for all those evaluated over time. There was no significant bidirectional interaction for the season factor, and the best performance corresponded to those carried out in the dry season. This allowed a greater expression of the productive potential as a food alternative for low rainfall areas.

Keywords:

adaptation, interaction, potential, stabilization.

Introduction

Implementing silvopastoral systems is very useful and favors the conservation of the environment; they involve tree and shrub species with high aptitude in biomass production suitable for bovine feeding and for other uses, such as shade (Uribe et al., 2011). One way to identify their capacity in the systems is to look at the productive capacity and the preferences of cattle for each of these in the grazing systems (Balehegn and Berhe, 2016).

There are species with good behavior and adaptation to different management techniques and adverse conditions present in the systems; this is reflected in the adequate productivity and supply of biomass per unit area, as is the case for Leucaena leucocephala (Francisco et al., 1998). Different studies attribute an important role in animal feeding to woody forage species (trees and shrubs).

The following species stand out: Albizia guachapele (Kunth) Dugand, Albizia niopoides (Benth.) Burkart, Albizia saman (Jacq.) Merr, and Leucaena leucocephala (Lam.) de Wit (Garcia et al., 2008; Uribe et al., 2011). Therefore, it is necessary to recognize their productive potential in terms of dry matter and incorporate their functional knowledge; in addition to providing biomass, they contribute to soil conservation, design efficiency, and maintenance or management of production systems that depend on the forage resource (Román-Miranda et al., 2017).

Additionally, it is of great importance to mention that many of these species with forage potential have excellent nutritional content, with high and low dry matter (DM) productions (Lombo, 2012). Some forage species of tree legumes have the ability to increase the level of nitrogen in the soil due to their ability to fix it from the atmosphere through the association that they have in the roots with nitrifying bacteria.

At the same time, they provide other ecosystem services due to nutrient cycling processes (Botero and Russo, 2002). This work aimed to determine the productive potential in terms of dry matter (DM) in two contrasting seasons for the leguminous species Iguá Albizia guachapele, Bayo Albizia niopoides, Saman Albizia saman, and Leucaena Leucaena leucocephala.

Materials and methods

Study site

Since 2013, production evaluations have been carried out for 22 tree and shrub species at the Nataima Research Center of the Colombian Agricultural Research Corporation (Agrosavia), located in the municipality of Espinal-Tolima, Colombia, Coordinates 4° 10’ 34” - 4° 11’ 26” north latitude and 74° 58’ 02” - 74° 50’ 44” west longitude at 391 masl, average temperature of 28 °C, average relative humidity of 65%, and annual rainfall of 1 400 m. This zone corresponds to Colombia’s tropical dry forest (TDF) (Serrano and Calderón, 2016).

Of the 22 species, the following were selected: Albizia guachapele, Albizia niopoides, Albizia saman, and Leucaena leucocephala, due to their excellent characteristics of adaptation to the area and high nutritional contribution. They were managed as shrubs of 1.8 to 2.0 m. They were used for the present research, carried out between March and August 2017.

A cutting or harvesting height of 80 cm was used, where a systematic random sampling (SRS) was carried out; it consisted of making daily homogenizations to groups of six plants (sampling unit), with systematic cuts and a daily cutting frequency until reaching forty days. All degraded leaves and branches less than 8 mm in diameter were considered. Each species had a population of 240 individuals of the same age planted in 1 x 1 m rows.

For the determination of dry matter (DM), 300 g of fresh fodder (FF) was taken from each species and packed in craft paper bags to facilitate the elimination of moisture; subsequently, by the traditional method, they were subjected to drying in an oven at temperatures between 60 and 65 °C for 24 to 72 h. Samples were weighed every 12 h until the time when the dry weight stabilized was determined. This process was carried out in the laboratories of Agrosavia, Espinal, Tolima.

The response variable was expressed in terms of dry matter (%) DM of forage, measured in grams (g), over 5.3 months in four production cycles of 40 days each (production times or replications).

Experimental design

The experiment used a design of repeated measures over time, with an arrangement of randomized complete blocks, four treatments (species) by two seasons, dry and wet, which resulted in eight interactions between the species and rainy season factors, 40 repetitions in 4 periods, giving 160 experimental units of six plants each.

During the species evaluation period from March to August for the year 2018, there were two well-defined and contrasting seasons in relation to the rainfall and temperatures recorded (Figure 1).

Figure 1

Figure 1. Accumulated precipitation and temperature from March to August 2018, Vintage station pro-2. Nataima Research Center, Colombian Agricultural Research Corporation (Agrosavia). Espinal-Tolima, Colombia.

2007-0934-remexca-15-05-e3462-gf5.jpg

Statistical analysis

The information obtained was analyzed through the R V.4.1.3 software (R Core Team, 2020) using the tidyverse (Wickham et al., 2019), ggpubr (Kassambara, 2020) and rstatix (Kassambara, 2021) packages. Assumptions were verified by identifying the presence of extreme outliers, normal distribution of the data with the Shapiro-Wilk statistic (p> 0.05), normal QQ plots for each time point, and the assumption of sphericity using the Mauchly statistic.

The effects of the treatments on the forage yield variable expressed in grams (g) of dry matter (DM) were analyzed through a unidirectional and bidirectional repeated measures analysis of variance and were considered significant at p< 0.05. In the presence of test significance, multiple paired t-tests were performed for the difference between treatments.

Results and discussion

For the species A. guachapele (Figure 2a), the data show that the dry matter production variable was statistically significant at the different times [F (1; 44; 56, 289) = 33.02; p= 0.0001, generalized eta squared= 0.46], which represents the high variability presented due to each of the cutting cycles over time. On the other hand, post hoc analyses with a Bonferroni adjustment revealed statistically significant differences for the DM variables between the pairwise comparisons for treatments (t), with p< 0.0001, between t1 and t2, t1 and t3, t1 and t4, t2 and t3, and t2 and t4, and highly significant between t3 and t4.

Figure 2

Figure 2. Effect of time on dry matter production (g) in A. guachapele (a); and A. niopoides (b), at the Nataima Research Center, Colombian Agricultural Research Corporation (Agrosavia). Espinal-Tolima, Colombia, 2018. The asterisks report the boundary of the significance value for which it is greater than the p-value. ****= probabilities less than 1%; ♦= it represents the arithmetic mean; ***= probability less than or equal to 1%; **= probability below 5%.

2007-0934-remexca-15-05-e3462-gf6.jpg

The fourth cut had the highest yields, averaging 205 g DM per plant. Treatments 1 and 2 presented yields of less than 200 g DM per plant, but with statistical differences between them due to the variability between cuts made for the four cycles.

For the species A. niopoides (Figure 2b), dry matter production was lower and exhibited statistically significant differences for the four times [F (1.64; 64.05) = 106.8, p= 0.0001 and a generalized eta squared= 0.73]. There was a variability effect due to cutting cycles. Time showed the highest value for eta in the four species evaluated, producing greater variability for species A. niopoides.

On the other hand, post hoc analyses with a Bonferroni adjustment revealed statistically significant differences in DM (g) between the comparisons, as those between t1 and t2, t1 and t3, t1 and t4, t2 and t3, t2 and t4, t2 and t4, and t3 and t4 had a p< 0.0001, where the fourth cut exhibited the highest yields and a mean of 160 g DM per plant, due to the variability that this species shows in relation to its growth.

For the species A. saman (Figure 3a), dry matter production was statistically significant at the different evaluation times [F (1.48; 57.63)= 42.34; p= 0.0001; generalized eta squared= 0.52]. There was high variability in the experiment due to the time factor (t) for the four moments evaluated.

Figure 3

Figure 3. Effect of time on dry matter production (g) in A. saman (a) and L. leucocephala (b) in the Nataima Research Center, Colombian Agricultural Research Corporation (Agrosavia). Espinal-Tolima, Colombia, 2018. ****= probabilities less than 1%; ♦= it represents the arithmetic mean; ***= probability less than or equal to 1%; **= probability below 5%.

2007-0934-remexca-15-05-e3462-gf7.jpg

Post hoc analyses with a Bonferroni adjustment revealed statistically significant differences in dry matter (g) between pairwise comparisons, p< 0.0001, between t1 and t2, t1 and t3, t1 and t4, t2 and t3, t2 and 4, and t3 and t4, where the first cut was below 150 g, and the fourth exhibited the highest yields with an average of 155 g DM per plant.

For the species L. leucocephala (Figure 3b), DM production was statistically significant at different evaluation times [F (1.66; 64.87)= 73.58; p= 0.0001 with generalized eta squared= 0.65], which indicated that there is a variability effect due to the t-factor. Post hoc analyses with a Bonferroni adjustment revealed statistically significant differences in DM (g) between pairwise comparisons, p< 0.0001, between t1 and t2, t1 and t3, t1 and t4, t2 and t3, t2 and t4, and t3 and t4, where the fourth cut showed the highest yields, with a mean of 158 g DM per plant.

All the species evaluated underwent a stabilization process over time. The availability of water resources and nutrients, among other production factors, was considered, which means that with longer recovery times, the species increases the production of DM (g). This is related to their physiology, morphology, and ability to adapt to the evaluated environmental conditions, represented by greater rooting and nutrient extraction (Cerón, 2013).

It was observed that there was no significant bidirectional interaction, which indicates that the impact that the species factor has on the dry matter variable does not depend on the levels for this (season) (Figure 4). Then, the main effects of each of the two factors (species and time) were interpreted.

Figure 4

Figure 4. Effect of the rainy season on dry matter (g) production in the legume species evaluated. Nataima Research Center, Colombian Agricultural Research Corporation (Agrosavia). Espinal-Tolima, Colombia, 2018.

2007-0934-remexca-15-05-e3462-gf8.jpg

The results of the analysis of variance were considered, and there were statistically significant main effects of the species factor [F (2.4; 189; 25)= 200.779, p= 0.0001], season factor [F (1.79)= 229.910, p< 0.0001] in the DM (g) variable. The cuts made for the dry season corresponded to those with the highest productivity recorded through the four evaluation moments.

Concerning the DM variable (g), the factors separately exhibited significant differences, where it was shown that for the species factor, A. guachapele had the highest yields in both the dry and wet seasons, and productions of 191.26 and 164.39 g, respectively, and statistically significant differences p< 0.0001 compared to the other treatments. These were not evident in the comparison made for the other treatments species A. niopoides, L. Leucocephala, and A. saman, and productions of 144.95, 142.67, and 140.39 g, respectively. The highest yields observed corresponded to the dry season, 164.39 g, and they were statistically different from the wet season, with 145.25 g harvested on average (Table 1).

Table 1

Table 1. Evaluation of dry matter yield in forage leguminous species for the independent factors, species and season.

Dry matter (g)
Species (legumes) A. guachapele 191.26a
A. niopoides 144.95 b
L. leucocephala 142.67 b
A. saman 140.39 b
Season (time) Dry 164.39 a
Wet 145.25 b

a , b, c= means with different letters in columns differ from each other according to Tukey’s test (p< 0.05)/ a, b, c. Nataima Research Center of the Colombian Agricultural Research Corporation (Agrosavia). Espinal-Tolima, Colombia. 2018.

Different trials have shown that soil and climatic conditions have an important effect on the survival and development of different plant species used in animal feeding (Lombo, 2012; Bueno et al., 2019). However, it is unclear how soil conditions, different rainfall, and the length of the dry period affect the individual behavior of each species (van Breugel et al., 2011).

In this sense, the species A. guachapele presented the highest DM yields for the four legume species evaluated; according to Lombo et al. (2012), this is due to the fact that this species, a pioneer in reforestation processes due to its rapid growth with great abundance in the secondary dry forest, showed a great capacity for adaptation and resistance to adverse climatic conditions, which makes it one of the main species to be used under conditions of water deficit and low fertility (Lombo, 2012; Guarachi, 2020).

It was shown that species development varies with edaphic conditions but that there are specific and differential responses to soil variables (Haggar et al., 1998). Therefore, species specialization to edaphic conditions not only shows adaptations to soil conditions but also the interaction of soil conditions with other covariates, such as pest pressure and phytosanitary disturbances (Fine et al., 2004).

There is evidence that sensitivity to drought varies greatly between species (Engelbrecht and Kursar, 2003; Poorter and Markesteijn, 2008); this was observed in the results obtained in this study developed in a dry environment, showing that the species presented similar survival but not the same production (DM) performance.

It is suggested that the history of the field regarding tillage practices affects the physicochemical properties of the soil by making these types of variables a more limiting factor than the variability in the distribution of precipitation itself (van Hughes et al., 2002; van-Breugel et al., 2011) All species showed a positive effect, with an increase in DM content through the two cutting moments, similar to what was found in other studies, where three cutting heights were evaluated for the species L. leucocephala, where there was an increase in its biomass as a function of the cutting height (Francisco et al., 1998).

The same study also showed that the higher the cutting height for the species, the higher its recovery rate; so, the cutting height of this evaluation, 40 cm, corresponded to half of the height to which the species were subjected in the present study (Francisco et al., 1998). This allowed these species to show greater and faster recovery and therefore, a progressive increase in DM content over time.

The biomass recorded for the wet season in this study was lower than that of the research carried out by (Francisco et al., 1998 ), where most of the yields for this trait occurred in the wet period; the above is due to the fact that the frequent and excessive rainfall of this period caused the cessation of metabolic and development processes, decreasing the production of biomass for these. Therefore, these species showed greater adaptation due to the fact that their origin presents extreme conditions for their adaptation (Pérez et al., 2015).

The lowest production corresponds to A. saman and L. leucocephala in the two periods evaluated, more for the wet season; it agrees with Olivarez et al. ( 2005) findings, where it is stated that the species L. leucocephala shows a reduced growth rate compared to other species, due to the fact that it has metabolisms that make it inefficient for the conversion of CO2 to carbohydrates, which are necessary for the formation of new tissues, limiting their development (Francisco, 1998).

At the same time, these species require well-drained soils without excess rainfall and good effective depth for their development, which leads to little development and potential for these species in the area (Botero and Russo, 2002). In contrast, lower yields have shown in the dry season (Latt et al., 2000). This is because there is less accumulation of carbohydrate reserves for these species for this dry season.

The species A. guachapele presented an excellent performance, especially for the dry season, with a DM weight of more than 200 g for the same period evaluated. This is due to its abundant growth, similar to Samán, and this agrees with what was reported by Pérez et al. (2015), where 19 species were evaluated, among which A. guachapele presented high productivity as well as survival with cutting periods of 30 days and 80 cm in height, it also showed the best adaptability behaviors to the conditions of Espinal-Tolima in Colombia.

Acknowledgments

The authors would like to thank the Colombian Agricultural Research Corporation (Agrosavia), where the Dynamic Agenda project, funded by the National Government, was carried out.

Bibliography

1 

Balehegn, M. and Berhe, K. 2016. Training reduced subjectivity of comparative yield method of estimation of grassland biomass. Grass and Forage Science. 71(3):482-489. 10.1111/gfs.12194.

2 

Botero, R. y Russo, R. 2002. Utilización de árboles y arbustos fijadores de nitrógeno en sistemas sostenibles de producción animal en suelos ácidos tropicales. Conferencia electrónica de la FAO sobre agroforestería para la producción animal en Latinoamérica. San José de Costa Rica. 121-143 pp. http://usi.earth.ac.cr/glas/sp/50000024.pdf.

3 

Bueno, A.; Pritsch, K. and Simon, J. 2019. Species-specific outcome in the competition for nitrogen between invasive and native tree seedlings. Frontiers in Plant Science. 10(337):1-17. 10.3389/fpls.2019.00337.

4 

Insuasty, S. E.; Apráez, G. y Gálvez, A. C. 2013. Caracterización botánica, nutricional y fenológica de especies arbóreas y arbustivas de bosque muy seco tropical. Ciencia Animal. (6):109-124.

5 

Engelbrecht, B. M. and Kursar, T. A. 2003. Comparative drought-resistance of seedlings of 28 species of co-occurring tropical woody plants. Oncología. 136(3):383-93. Doi: 10.1007/s00442-003-1290-8.

6 

Francisco, A. G. 1998. Efecto de tres alturas de corte en el rendimiento de biomasa de Leucaena leucocephala cv. CNIA-250. Pastos y Forrajes. 21(4):337-343.

7 

Fine, A. P.; Mesones, V. I. and Coley, P. D. 2004. Herbivores promote habitat specialization by trees in Amazonian forests. Science. 305(5684):663-665. https://doi.org/10.1126/science.1098982.

8 

Francisco, G.; Simon, L. y Soca. M. 1998. Efecto de tres alturas de corte en el rendimiento de biomasa de Leucaena leucocephala cv. CNIA-250. Pastos y Forrajes . 21(4):337-343.

9 

García, D. E.; Medina, M. G. L.; Cova, J.; Torres, A.; Soca, M.; Pizzani, P.; Baldizán, A. y Domínguez, C. E. 2008. Preferencia de vacunos por el follaje de doce especiescon potencial para sistemas agrosilvopastorilesen el Estado Trujillo, Venezuela. Pastos y Forrajes , 31(3):255-268 https://payfo.ihatuey.cu/index.php?journal=pasto&page=article&op=view&path%5B%5D=673&path%5B%5D=175.

10 

Guarachi, C. A. 2020. Efecto de la aplicación de dos tratamientos pre-germinativos y tres componentes de sustrato en la germinación de semillas de cedro amarillo (Albizia guachapele), en el municipio de el alto. Tesis profesional ingeniero agropecuario. Universidad Catolica Santiago de Guayaquil. Repositorio Rigital UCSG, Facultad de Educación Tecnica para el Desarrollo. https://repositorio.umsa.bo/bitstream/handle/123456789/25321/T-2790.pdf?sequence=1&isAllowed=y.

11 

Haggar, J. P.; Briscoe, C. B. and Butterfield, R. P. 1998. Native species: a resource for the diversification of forestry production in the lowland humid tropics. Forest ecology and management. 106(1-3):195-203. 10.1016/S0378-1127(97)00311-3.

12 

Hughes, R. F.; Kauffman, J. B. and Cummings, D. L. 2002. Dynamics of aboveground and soil carbon and nitrogen stocks and cycling of available nitrogen along a land-use gradient in Rondônia, Brazil. Ecosystems. 5(3):244-259. https://doi.org/10.1007/s10021-001-0069-1.

13 

Kassambara, A. 2020. ggpubr: ‘ggplot2’ Based Publication R Plots. R package version 0.4.0. https://CRAN.R-project.org/package=ggpubr.

14 

Kassambara, A. 2021. rstatix: pipe-friendly framework for basic statistical tests. R package version 0.7.0. https://CRAN.R-project.org/package=rstatix.

15 

Lombo-Ortiz, D. 2000. Evaluación de la disponibilidad de biomasa y capacidad de rebrote de leñosas forrajeras en potreros del trópico seco de Nicaragua. Tesis de maestría. Escuela de posgrados del CATIE. 1-25 pp. https://repositorio.catie.ac.cr/bitstream/handle/11554/8031/252.pdf?sequence=2&isAllowed=y

16 

Olivarez, P. J.; Jimenez, R. G.; Rojas, S. E. y Martinez, P. A. H. 2005. Uso de las leguminosas arbustivas en los sistemas de producción animal en el trópico. Redvet. Revista Electrónica de Veterinaria. 6(5):1-19. https://www.redalyc.org/articulo.oa?id=63617216006.

17 

Pérez, N.; Rojas. C.; Crillo, D.; Moreno, J. y Guzmán, L. 2015. Rasgos de especies forrajeras basados en diversidad funcional asociados a la producción de biomasa que contribuyen a la decisión de ganaderos en zonas secas. 317-322 pp. https://www.researchgate.net/publication/335472530-Rasgos-de-especies-forrajeras-basados-en-diversidad-funcional-asociados-a-la-produccion.

18 

Poorter. L., and MarkesteijnL. 2008. Seedling traits determine drought tolerance of tropical tree species. Biotropica. 40(3):321-331. https://doi.org/10.1111/j.1744-7429.2007.00380.x

19 

R Core Team. 2020. A language and environment for statistical computing. R Foundation for Statistical Computing.

20 

Roman-Miranda, M. L.; Palma-Garcia, J. M.; Zorrilla-Rios. J. M. y Mora-Santacruz, A. 2017. Producción de materia seca de Leucaena leucocephala y vegetación herbácea en un banco de proteína pastoreada por ovinos. Revista de Sistemas Esperimentales. 3(6):1-9.

21 

Serrano, P. R. y Calderón, F. 2016. Análisisn geoespacial del cambio de las zonas de vida de Holdridge en la provincia de Guayas. In: 14th LACCEI international multi conference for engineering, education, and technology: engineering innovations for global sustainability. 1-7 pp. 10.18687/LACCEI2016.1.1.185.

22 

Uribe, F. T. A. F.; Zuluaga, T. E.; Murgueitio, R. L. M.; Valencia, C. A.; Zapata, C. L. H.; Solarte, P. C. A. Cuartas Naranjo, C. F.; Galindo, W. F. S.; Gonzales, J. G S.; Sinisterra J. A. R.; Gomez, J. G. B.; Molina, C. E. D.; Molina, E. J. D.; Galindo, A. O.; Galindo, V. A. C. y Soto, B. R. 2011. Establecimiento y manejo de sistemas silvopastoriles. Proyecto Ganadería Sostenible Manual1. CIPAV. 1(1):1-78. http://ganaderiacolombianasostenible.co/web/wpcontent/uploads/2015/04/1.establecimiento-y-manejo-de-SSP.pdf.

23 

Van-Breugel, M.; Hall, J. S.; Craven, D.; Gregoire, T. G.; Park. A.; Dent, D.; Wishnie, M.; Mariscal, E.; Deago, J. Ibarra, D.; Cedeño, N. and Ashton, M. S. 2011. Early growth and survival of 49 tropical tree species across sites differing in soil fertility and rainfall in Panama. Forest ecology and management . 261(10):1580-1589. https://doi.org/10.1016/j.foreco.2010.08.019.

24 

Wickham, H. M.; Averick. J.; Bryan, W.; Chang, L.; D’Agostino, M. L.; François, R.; Grolemund, G.; Hayes, A.; Henry, L.; Hestre, J.; Kuhn, M.; Pedersen, T. L.; Miller, E.; Bache, S. M.; Muller, K.; Ooms, J.; Robinson, D.; Seide, D. P.; Spinu, V. and Yutani, H. 2019. Welcome to the tidyverse. Journal of Open-Source Software. 4(43):1-6 https://doi.org/10.21105/joss.01686.