Vol. 27 Núm. 1 (2024): Enero-junio


Artículos de investigación científica y tecnológica

Biomass Storage Potential and Improvement in Soil Properties under Different Bamboo Plantations in the Terai Region of Central Himalaya

Almacenamiento de biomasa y mejora del suelo en plantaciones de bambú en la región de Terai en el Himalaya central


Palabras clave:

Bamboo species, nutrient uptake, aboveground biomass, soil microbial population, soil properties, Terai region (en).

Palabras clave:

Especies de bambú, absorción de nutrientes, biomasa aérea, población microbiana del suelo, propiedades del suelo, región de Terai (es).


Akinlabi, E. T., Anane-Fenin, K., & Akwada, D. R. (2017). Bamboo: The multipurpose plant. Springer. DOI:

Alemayehu, A., Mulatu, Y., Eshete, N., & Terefe, M. (2015). Growth performance and biomass accumulation of four different introduced bamboo species in South-Western Ethiopia. Growth, 5(3), 5-10.

Arce, J. J. C. (2019). Forests, inclusive and sustainable economic growth and employment. United Nations.

Arunachalam, A., & Arunachalam, K. (2002). Evaluation of bamboos in eco-restoration of ‘jhum’fallows in Arunachal Pradesh: Ground vegetation, soil and microbial biomass. Forest Ecology and Management, 159(3), 231-239. DOI:

Bai, W., Kong, L., & Guo, A. (2013). Effects of physical properties on electrical conductivity of compacted lateritic soil. Journal of Rock Mechanics and Geotechnical Engineering, 5(5), 406-411. DOI:

Benbi, D. K., Brar, K., Toor, A. S., & Singh P. (2015). Total and labile pools of soil organic carbon in cultivated and undisturbed soils in northern India. Geoderma, 237, 149-158. DOI:

Blake, G. R., & Hartge, K. H. (1986). Methods of soil analysis. American Society of Agronomy.

Bower, C. A., & Wilcox, L. V. (1965). Methods of soil analysis. American Society of Agronomy.

Cao, S. X., Chen, L., Shankman, D., Wang, C. M., Wang, X. B., & Zhang, H. (2011) Excessive reliance on afforestationin China’s arid and semi-arid regions, Lessons in ecological restoration. Earth Science Review, 104, 240-245. DOI:

Chang, E. H., Tian, G., Shiau, Y. J., Chen, T. H., & Chiu, C. Y. (2019). Influence of thorny bamboo plantations on soil microbial biomass and community structure in subtropical badland soils. Forests, 10(10), 854. DOI:

Chen, X. (2021). Variations in patterns of internode and branch lengths for several bamboo species. Plant Biosystems, 155(6), 1088-1099. DOI:

Chen, M., Lin, G., Ramakrishnan, M., Fei, Z., Kunnummal, K. V., Yulong, D., & Chen, J. (2022). Rapid growth of Moso bamboo (Phyllostachys edulis): Cellular roadmaps, transcriptome dynamics, and environmental factors. The Plant Cell, 34(10), 3577-3610. 10.1093/plcell/koac193 DOI:

Chhonkar, P. K., Bhadraray, S., Patra, A. K., & Purakayastha T. J. (2007). Experiments in soil biology and biochemistry. Westville Publishing House.

Danielson, R. E., & Southerland, P. L. (1986). Porosity. American Society of Agronomy, Soil Science Society of America. DOI:

Dong, L., & Kou, M. (2022). Soil Aggregate stability and carbon density in three plantations in the Loess Plateau, China. Forests, 13(7), 1096. DOI:

El-Amin, E. A., Diab, I. E., & Ibrahim, S. (2001). Influence of Eucalyptus cover on some physical and chemical properties of a soil in Sudan. Communications in Soil Science and Plant Analysis, 32(13-14), 2267-2278. DOI:

Farooq, T. H., Shakoor, A., Wu, X., Li, Y., Rashid, M. H. U., Zhang, X., & Yan, W. (2021). Perspectives of plantation forests in the sustainable forest development of China. Forest-Biogeosciences and Forestry, 14(2), 166-174. DOI:

Franklin, D. C., Prior, L. D., Hogarth, N. J., & McMahon, C. R. (2010). Bamboo, fire and flood: Consequences of disturbance for the vegetative growth of a clumping, clonal plant. Plant Ecology, 208(2), 319-332. DOI:

Fu, Y., Horton, R., & Heitman J. (2021). Estimation of soil water retention curves from soil bulk electrical conductivity and water content measurements. Soil and Tillage Research, 209, 104948. DOI:

Gaikwad, A. S., Thite, M., & Rananavare, S. I. A. (2021). Effect of different bamboo species on growth attributing characters grown on Entisol of semi-arid climate. The Pharma Innovation Journal, 10(12), 1980-1985. DOI:

Goh, Y., Yap, S. P., & Tong, T. Y. (2020). Bamboo: The emerging renewable material for sustainable construction. Encyclopedia of Renewable and Sustainable Materials, 2, 365-376. https://doi:org/10.1016/B978-0-12-803581-8.10748-9 DOI:

Guo, Z., Zhuang, M., Yang, L., Li, Y., Wu, S., & Chen, S. (2021). Differentiated mineral nutrient management in two bamboo species under elevated CO2 environment. Journal of Environmental Management, 279, 111600. DOI:

Hao, M., Zhang, J., Meng, M., Chen, H. Y., Guo, X., Liu, S., & Ye, L. (2019). Impacts of changes in vegetation on saturated hydraulic conductivity of soil in subtropical forests. Scientific Reports, 9(1), 8372. DOI:

Hassan, N. H. M., Abdullah, N., Kelanda, D. N. A., & Perumal, M. (2022). Early field growth performance of ten selected bamboo taxa: The case study of Sabal bamboo pilot project in Sarawak, Malaysia. Biodiversitas, 23(6), 2882-2892. DOI:

Iyer, R. (2019). Sustainability-based evaluation of bamboo harvest in Pachgaon [Doctoral dissertation, Department of Biology, IISRE].

Jackson, M. L. (1958). Soil chemical analysis. Asia Publishing House.

Jackson, M. L. (1967). Soil chemical analysis. Prentice Hall.

Kanwal, M. S. (2011). Influence of different bamboo species on growth, carbon sequestration potential and soil properties in tarai region of northern india [Master’s thesis, G. B. Pant University of Agriculture and Technology, Department of Agronomy.

Kaushal, R., Singh, I., Thapliyal, S. D., Gupta, A. K., Mandal, D., Tomar, J. M. S., & Durai, J. (2020 a). Rooting behaviour and soil properties in different bamboo species of Western Himalayan Foothills, India. Scientific Reports, 10(1), 4966. DOI:

Kaushal, R., Tewari, S., Banik, R. L., Thapliyal, S. D., Singh, I., Reza, S., & Durai, J. (2020 b). Root distribution and soil properties under 12-year old sympodial bamboo plantation in Central Himalayan Tarai Region, India. Agroforestry Systems, 94, 917-932. DOI:

Kaushal, R., Islam, S., Tewari, S., Tomar, J. M. S., Thapliyal, S., Madhu, M., & Durai, J. (2022). An allometric model-based approach for estimating biomass in seven Indian bamboo species in western Himalayan foothills, India. Scientific Reports, 12(1), 7527. DOI:

Kim, C., Baek, G., Yoo, B. O., Jung, S. Y., & Lee, K. S. (2018). Regular fertilization effects on the nutrient distribution of bamboo components in a moso bamboo (Phyllostachys pubescens (Mazel) Ohwi) stand in South Korea. Forests, 9(11), 671. DOI:

Kittur, B. H., Sudhakara, K., Mohan Kumar, B., Kunhamu, T. K., & Sureshkumar, P. (2016). Bamboo based agroforestry systems in Kerala, India: performance of turmeric (Curcuma longa L.) in the subcanopy of differentially spaced seven year-old bamboo stand. Agroforestry Systems, 90, 237-250. 10.1007/s10457-015-9849-z DOI:

Kodešová, R., Kodeš, V., Žigová, A., & Šimůnek, J. (2006). Impact of plant roots and soil organisms on soil micromorphology and hydraulic properties. Biologica, 61(19), 339-343. DOI:

Kumar, B. M., Rajesh, G., & Sudheesh, K. G. (2006). Aboveground biomass production and nutrient uptake of thorny bamboo [Bambusa bambos (L.) Voss] in the homegardens of Thrissur, Kerala. Journal of Tropical Agriculture, 43, 51-56.

Kumar, R., Thangaraju, M. M., Kumar, M., Thul, S. T., Pandey, V. C., Yadav, S., & Kumar, S. (2021). Ecological restoration of coal fly ash–dumped area through bamboo plantation. Environmental Science and Pollution Research, 28, 33416-33432. 10.1007/s11356-021-12995-7 DOI:

Kumar, P. S., Shukla, G., Nath, A. J., & Chakravarty, S. (2022). Soil properties, litter dynamics and biomass carbon storage in three-bamboo species of Sub-Himalayan region of Eastern India. Water, Air, & Soil Pollution, 233(1), 12. DOI:

Kumari, Y., & Bhardwaj, D. R. (2017). Distribution of nutrition elements in different bamboos in bamboo based agroforestry plantation: A comparative study. Journal of Pharmacognosy and Phytochemistry, 6(3), 556-561.

Kumari, Y., & Bhardwaj, D. R. (2019). Bioproduction and nutrient cycling in seven bamboo species in subtropical Indian Himalayas. Journal of Tree Sciences, 38(2), 1-12. DOI:

Lei, Z., Yu, D., Zhou, F., Zhang, Y., Yu, D., Zhou, Y., & Han, Y. (2019). Changes in soil organic carbon and its influencing factors in the growth of Pinus sylvestris var. mongolica plantation in Horqin Sandy Land, Northeast China. Scientific Reports, 9(1), 16453. DOI:

Li, C., Cai, Y., Xiao, L., Gao, X., Shi, Y., Zhou, Y., & Zhou, G. (2021). Rhizome extension characteristics, structure and carbon storage relationships with culms in a 10-year moso bamboo reforestation period. Forest Ecology and Management, 498, 119556. DOI:

Liu, R., Yang, X., Gao, R., Hou, X., Huo, L., Huang, Z., & Cornelissen, J. H. (2021). Allometry rather than abiotic drivers explains biomass allocation among leaves, stems and roots of Artemisia across a large environmental gradient in China. Journal of Ecology, 109(2), 1026-1040. DOI:

Lobovikov, M., Paudel, S., Ball, L., Piazza, M., Guardia, M., Wu, J., & Ren, H. (2007). World bamboo resources: A thematic study prepared in the framework of the global forest resources assessment 2005 (No. 18). Food & Agriculture Organization.

Luo, X., Keenan, T. F., Chen, J. M., Croft, H., Colin Prentice, I., Smith, N. G., & Zhang, Y. (2021). Global variation in the fraction of leaf nitrogen allocated to photosynthesis. Nature Communications, 12(1), 4866. DOI:

Mensah, S., Kakaï, R. G., & Seifert, T. (2016). Patterns of biomass allocation between foliage and woody structure: the effects of tree size and specific functional traits. Annals of Forest Research, 59(1), 49-60. 10.15287/afr.2016.458 DOI:

Nawaz, M. F., Bourrie, G., & Trolard, F. (2013). Soil compaction impact and modelling. A review. Agronomy for Sustainable Development, 33, 291-309. DOI:

Noguchi, M., & Yoshida, T. (2005). Factors influencing the distribution of two co-occurring dwarf bamboo species (Sasa kurilensis and S. senanensis) in a conifer-broadleaved mixed stand in northern Hokkaido. Ecological Research, 20, 25-30. DOI:

Oli, B. N., & Kandel, C. M. (2005). Biomass estimation of Bambusa nutans subspecies cupulata grown at Eastern Terai, Nepal. Banko Janakari, 15(2), 34-37. DOI:

Ovington, J. D. (1965). Organic production, turnover and mineral cycling in woodlands. Biological Reviews, 40(3), 295-336. DOI:

Pandey, C. S. (2019). Influence of tree cropping systems on physical properties of salt affected soil. International Journal of Current Microbiology and Applied Science, 8(4), 739-747. DOI:

Pattanaik, S., & Hall, J. B. (2014). Patterns of morphometric variability in Dendrocalamus hamiltonii Munro populations across East Khasi Hills, Northeast India. Indian Forester, 140(9), 868-874.

Patra, S., Kaushal, R., Singh, D., Kumar, R., Gadedjisso-Tossou, A., & Durai, J. (2022). Surface soil hydraulic conductivity and macro-pore characteristics as affected by four bamboo species in North-Western Himalaya, India. Ecohydrology & Hydrobiology, 22(1), 188-196. DOI:

Puangchit, L., Hnin, S. M., & Sungkaew, S. (2019). Allometric equations for estimating the aboveground biomass of a 14-Year-old bamboo plantation at moeswe research station, Myanmar. Journal of Tropical Forest Research, 3(1), 1-19.

Qian, Z., Sun, X., Gao, J., & Zhuang, S. (2021). Effects of bamboo (Phyllostachys praecox) cultivation on soil nitrogen fractions and mineralization. Forests, 12(8), 1109. DOI:

Qiao, Y., Miao, S., Yue, S., Wu, H., & Han, X. (2016). How 23-year continuous soybean cultivation led to more SOC and thermal energy stored in Mollisol micro-aggregates. Polish Journal of Environmental Studies, 25(3), 1215-1221. DOI:

Robinson, D. A., Thomas, A., Reinsch, S., Lebron, I., Feeney, C. J., Maskell, L. C., & Cosby, B. J. (2022). Analytical modelling of soil porosity and bulk density across the soil organic matter and land-use continuum. Scientific Reports, 12(1), 7085. DOI:

Seobi, T., Anderson, S. H., Udawatta, R. P., & Gantzer, C. J. (2005). Influence of grass and agroforestry buffer strips on soil hydraulic properties for an Albaqualf. Soil Science Society of America Journal, 69(3), 893-901. DOI:

Selassie, Y. G., & Ayanna, G. (2013). Effects of different land use systems on selected physico-chemical properties of soils in Northwestern Ethiopia. Journal of Agricultural Science, 5(4), 112. DOI:

Shanmughavel, P., & Francis, K. (2003). Biomass accumulation and nutrient distribution in Dendrocalamus hamiltonii [Conference presentation]. XII World Forestry Congress, Quebec, Canada.

Sharma, M. L., & Nirmala, C. (2015). Bamboo diversity of India: An update [Conference presentation]. 10th World Bamboo Congress, Seoul, South Korea.

Shiau, Y. J., Wang, H. C., Chen, T. H., Jien, S. H., Tian, G., & Chiu, C. Y. (2017). Improvement in the biochemical and chemical properties of badland soils by thorny bamboo. Scientific Reports, 7(1), 40561. DOI:

Silva, L. N., Freer-Smith, P., & Madsen, P. (2019). Production, restoration, mitigation: a new generation of plantations. New Forests, 50(2), 153-168. DOI:

Singh, J., Sharma, R., Dhakad, A. K., & Chauhan, S. K. (2018). Defining growth, quality and biomass production of different bamboo species in central plains of Punjab. Journal of Pharmacognosy and Phytochemistry, 7(5), 1328-1332.

Singh, K. A., & Kochhar, S. K. (2005). Effect of clump density/spacing on the productivity and nutrient uptake in Bambusa pallida and the changes in soil properties. Journal of Bamboo Rattan, 4(4), 323-334. DOI:

Singh, A. N., & Singh, J. S. (2006). Experiments on ecological restoration of coal mine spoil using native trees in a dry tropical environment, India: A synthesis. New Forests, 31, 25-39. DOI:

Singh, K. A., & Rai, A. K. (2012). Nutrient accumulation, distribution and use efficiency in different bamboo plant species in north-eastern hills region of India. Indian Journal of Agricultural Sciences, 82(2), 134. DOI:

Subbaiah, B. V., & Asija, G. L. (1956). A rapid method for estimation of available N in soil. Current Science, 25, 259-260.

Tian, G., Justicia, R., Coleman, D. C., & Carroll, C. R. (2007). Assessment of soil and plant carbon levels in two ecosystems (woody bamboo and pasture) in montane Ecuador. Soil Science, 172(6), 459-468. DOI:

Trivedi, P., Leach, J. E., Tringe, S. G., Sa, T., & Singh, B. K. (2020). Plant-microbiome interactions: From community assembly to plant health. Nature Reviews Microbiology, 18(11), 607-621. DOI:

Tu, Z., Chen, L., Yu, X., & Zheng, Y. (2013). Effect of bamboo plantation on rhizosphere soil enzyme and microbial activities in coastal ecosystem. Journal of Food, Agriculture and Environment, 2(3), 2333-2338.

Udawatta, R. P., Kremer, R. J., Garrett, H. E., & Anderson, S. H. (2009). Soil enzyme activities and physical properties in a watershed managed under agroforestry and row-crop systems. Agriculture, Ecosystems & Environment, 131(1-2), 98-104. DOI:

Walkley, A., & Black, I. A. (1934). An examination of the Degtjareff method for determining soil organic matter, and a proposed modification of the chromic acid titration method. Soil Science, 37(1), 29-38. DOI:

Wang, Y., Qiao, G., Xu, J., Jin, K., Fan, M., Ding, Y., & Zhuo, R. (2022). Anatomical characteristics and variation mechanisms on the thick-walled and dwarfed culm of Shidu bamboo (Phyllostachys nidularia f. farcta). Frontiers in Plant Science, 13, 1-11. DOI:

Watanabe, F. S., & Olsen, S. R. (1965). Test of an ascorbic acid method for determining phosphorus in water and NaHCO3 extracts from soil. Soil Science Society of America Journal, 29(6), 677-678. DOI:

Wu, J., Xu, Q., Jiang, P., & Cao, Z. (2009). Dynamics and distribution of nutrition elements in bamboos. Journal of Plant Nutrition, 32(3), 489-501. DOI:

Xie, J., Huang, X., Qi, J., Hse, C. Y., & Shupe, T. F. (2014). Effect of anatomical characteristics and chemical components on microwave-assisted liquefaction of bamboo wastes. BioResources, 9(1), 1-9. DOI:

Xue, L., Ren, H., Li, S., Leng, X., & Yao, X. (2017). Soil bacterial community structure and co-occurrence pattern during vegetation restoration in karst rocky desertification area. Frontiers in Microbiology, 8, 2377. DOI:

Zhan, H., Zhang, L. Y., Niu, Z. H., Wang, C. M., & Wang, S. G. (2016). Chemical properties and fiber morphology of Dendrocalamus hamiltonii as potential pulp material. European Journal of Wood and Wood Products, 74, 273-276. DOI:

Zhang, Y. M., Yu, Y. L., & Ji, Y. W. (2013). Effect of thermal treatment on the physical and mechanical properties of Phyllostachys pubescen bamboo. European Journal of Wood and Wood Products, 71(1), 61-67. 10.1007/s00107-012-0643-6. DOI:

Zhang, M. M., Fan, S. H., Guan, F. Y., Yan, X. R., & Yin, Z. X. (2020). Soil bacterial community structure of mixed bamboo and broad-leaved forest based on tree crown width ratio. Scientific Reports, 10(1), 6522. DOI:

Zhao-Hua, L. I., Denich, M., & Borsch, T. (2005). Growth behavior of Phyllostachys nigra var. henonis (Bambusoideae) in Central China. Journal of Forestry Research, 16, 163-168. DOI:

Zheng, Y., & Lin, X. (2020). Niche specialization and functional overlap of bamboo leaf and root microbiota. Frontiers in Microbiology, 11, 571159. DOI:

Zhou, Z., & Shangguan, Z. (2007). Vertical distribution of fine roots in relation to soil factors in Pinus tabulaeformis Carr. forest of the Loess Plateau of China. Plant and Soil, 291, 119-129. DOI:

Cómo citar


Sharma, N., Kumar, D., N.R., S., Yeasin, M., y ., B. (2024). Biomass Storage Potential and Improvement in Soil Properties under Different Bamboo Plantations in the Terai Region of Central Himalaya. Colombia forestal, 27(1), e20898.


Sharma, N. et al. 2024. Biomass Storage Potential and Improvement in Soil Properties under Different Bamboo Plantations in the Terai Region of Central Himalaya. Colombia forestal. 27, 1 (ene. 2024), e20898. DOI:


Sharma, N.; Kumar, D.; N.R., S.; Yeasin, M.; ., B. Biomass Storage Potential and Improvement in Soil Properties under Different Bamboo Plantations in the Terai Region of Central Himalaya. Colomb. for. 2024, 27, e20898.


SHARMA, Nnidhi; KUMAR, Dhirender; N.R., Sudhakara; YEASIN, Mohammad; ., Bharti. Biomass Storage Potential and Improvement in Soil Properties under Different Bamboo Plantations in the Terai Region of Central Himalaya. Colombia forestal, [S. l.], v. 27, n. 1, p. e20898, 2024. DOI: 10.14483/2256201X.20898. Disponível em: Acesso em: 16 abr. 2024.


Sharma, Nnidhi, Dhirender Kumar, Sudhakara N.R., Mohammad Yeasin, y Bharti . 2024. «Biomass Storage Potential and Improvement in Soil Properties under Different Bamboo Plantations in the Terai Region of Central Himalaya». Colombia forestal 27 (1):e20898.


Sharma, N. (2024) «Biomass Storage Potential and Improvement in Soil Properties under Different Bamboo Plantations in the Terai Region of Central Himalaya», Colombia forestal, 27(1), p. e20898. doi: 10.14483/2256201X.20898.


N. Sharma, D. Kumar, S. N.R., M. Yeasin, y B. ., «Biomass Storage Potential and Improvement in Soil Properties under Different Bamboo Plantations in the Terai Region of Central Himalaya», Colomb. for., vol. 27, n.º 1, p. e20898, ene. 2024.


Sharma, Nnidhi, et al. «Biomass Storage Potential and Improvement in Soil Properties under Different Bamboo Plantations in the Terai Region of Central Himalaya». Colombia forestal, vol. 27, n.º 1, enero de 2024, p. e20898, doi:10.14483/2256201X.20898.


Sharma, Nnidhi, Dhirender Kumar, Sudhakara N.R., Mohammad Yeasin, y Bharti . «Biomass Storage Potential and Improvement in Soil Properties under Different Bamboo Plantations in the Terai Region of Central Himalaya». Colombia forestal 27, no. 1 (enero 1, 2024): e20898. Accedido abril 16, 2024.


Sharma N, Kumar D, N.R. S, Yeasin M, . B. Biomass Storage Potential and Improvement in Soil Properties under Different Bamboo Plantations in the Terai Region of Central Himalaya. Colomb. for. [Internet]. 1 de enero de 2024 [citado 16 de abril de 2024];27(1):e20898. Disponible en:

Descargar cita






Los datos de descargas todavía no están disponibles.

Recibido: 6 de marzo de 2023; Aceptado: 8 de agosto de 2023


El concepto de productividad forestal es biológicamente complejo. En masas forestales, se estima indirectamente mediante modelos de crecimiento en altura dominante e índice de sitio. El objetivo de este estudio fue desarrollar un modelo de altura dominante para Gmelina arborea con datos provenientes de 129 parcelas de monitoreo. Se generaron 9 ecuaciones dinámicas de altura dominante a partir de los modelos teóricos de Hossfeld y Bertalanffy-Richards, ajustadas individualmente y sujetas a validación cruzada. De los modelos evaluados, sobresalió el modelo polimórfico de Bertalanffy-Richards (parámetro c). Este modelo presentó una tendencia biológica realista, con curvas polimórficas y una única asíntota. El modelo seleccionado es una herramienta importante para el manejo de plantaciones forestales de G. arborea en zonas que presenten características de sitio semejantes a las de este estudio.

Palabras clave:

curvas de índice de sitio, ecuaciones dinámicas, polimorfismo, productividad forestal, validación cruzada.


The concept of forest productivity is biologically complex. In forest masses, it is indirectly estimated via dominant height growth models and site indices. The aim of this study was to develop a dominant height model for Gmelina arborea using data from 129 monitoring plots. 9 dynamic equations were generated from the Hossfeld and Bertalanffy-Richards theoretical models, which were individually fitted and cross-validated. Among the evaluated models, the Bertalanffy-Richards model (parameter c) stood out. This model showed a logical biological trend, with polymorphic curves and a unique asymptote. The selected model is an important tool for the management of G. arborea forest plantations in areas with site characteristics similar to those of this study.


site index curves, dynamic equations, polymorphism, forest productivity, cross-validation.


La estimación del crecimiento futuro es una cuestión fundamental en la planificación de las inversiones y la ordenación forestal. Estas proyecciones permiten proyectar la respuesta de un sitio frente a una especie en particular, disminuyendo el sesgo en los supuestos, lo que las hace una herramienta indispensable para la planificación forestal (Aguirre & Aguirre, 2014). Los proyectos forestales, al ser afectados de forma directa por la variación del clima (Sabatia & Burkhart, 2014), necesitan de herramientas que permitan proyectar tasas diferenciales de productividad para realizar prescripciones silvícolas más precisas.

El potencial productivo de las masas forestales habitualmente se determina evaluando la calidad del sitio (Burkhart & Tomé, 2012), que a su vez se evalúa mediante un indicador conocido como el índice de sitio (IS), el cual se ve poco afectado por la densidad de plantación. Este índice consiste en unas curvas de crecimiento donde se relaciona una altura dominante proyectada hasta una edad base definida (Tamarit-Urias et al., 2022). Estas curvas muestran la respuesta de la especie al sitio, relacionando directamente los niveles de productividad de una especie con una ecorregión (Akbas & Senyurt, 2018).

Los métodos para la construcción de las curvas de IS se clasifican en cuatro grupos: el método de la curva guía, el enfoque por diferencias algebraica (EDA) desarrollado por Bailey & Clutter (1974), el método de predicción de parámetros (Clutter et al, 1983) y el enfoque por diferencias algebraicas generalizadas (EDAG) desarrollado por Cieszwski & Bailey (2000). Para el modelado dinámico, es necesario disponer de datos medidos en árboles o rodales al menos en dos ocasiones utilizando parcelas de intervalo, parcelas permanentes o análisis troncal (Hernández-cuevas et al., 2018). El método EDA involucra esencialmente la sustitución de un parámetro del modelo base para expresarlo como una función del sitio. Cieszwski & Bailey (2000) introdujeron una generalización a la metodología EDA (i.e., enfoque EDAG), considerando que la ecuación base se expande para permitir que más de un parámetro dependa de la calidad de sitio y que las familias de curvas obtenidas incrementen su flexibilidad. Una de las ventajas de la metodología EDAG es su capacidad para ampliar la base de ecuaciones de acuerdo con varias teorías sobre las características del crecimiento (asíntota, tasa de crecimiento), incluyendo la posibilidad de simular polimorfismos concurrentes y múltiples asíntotas, que son vistos como propiedades importantes de las ecuaciones de sitio (Guerra-Hernández et al., 2021).

En Colombia, las plantaciones forestales cubren 538 000 ha, localizadas en la Región Andina, la Costa Atlántica y la Orinoquía (MADR, 2022), estimando tan solo un 2.2 % de la superficie potencial para el establecimiento de plantaciones (Gómez et al., 2012). Algunas especies forestales tardan menos tiempo en crecer y ser productivas en comparación con otros países. Para el caso de la Gmelina arborea, se reporta un rendimiento de 20-25 m3.ha.año-1, mientras que, en promedio, países como Bolivia, Brasil, Costa Rica, Filipinas, Indonesia, Malawi y Malasia reportan entre 12 y 19 m3.ha.año-1. G. arborea se encuentra dentro de las diez especies con alto potencial de incrementar sus rendimientos por aumento en la producción, lo que a su vez se asocia con la disminución en los tiempos de rotación, que varían entre los 10 y 14 años (Escobar Quemba, 2013). Esta especie es plantada especialmente en ecosistemas de bosque seco tropical, alcanzando un área estimada de 21 000 ha para la industria de la madera en Colombia (MADR, 2022), lo que la convierte en una de las especies más promisorias para el desarrollo forestal del país (Melo Cruz, 2015). Estudios previos han modelado el crecimiento de Gmelina arborea, proyectando índices de sitio a diferentes intervalos. Esto permite conocer la respuesta de la especie a las condiciones de sitio. Este tipo de investigaciones ha sido reportado para países como Filipinas, México y Nigeria (Vanclay et al., 2008; Martínez-Zurimendi et al., 2015; Nurudeen et al., 2018). En Colombia se desarrolló un modelo empírico para predecir la altura dominante (H) y el IS de plantaciones de G. arborea, relacionando variables biofísicas y de rodal (Barrios et al., 2021). López et al. (2011) modelaron la altura dominante en la costa Atlántica y la región Andina de Colombia bajo el EDA, en el cual los modelos de crecimiento de Chapman-Richards y Korf presentaron el mejor ajuste para la especie. Así, este modelo se convirtió en un predictor caracterizado por la generación de curvas polimórficas con única asíntota.

El objetivo principal del presente estudio es generar un modelo predictor de productividad de sitio para Gmelina arborea en Colombia a partir de datos recolectados en los departamentos de Magdalena y Tolima en múltiples temporalidades. Para ello, se evalúan la calidad de ajuste y la capacidad predictiva en modelos dinámicos derivados de los enfoques EDA y EDAG, seleccionando la ecuación que mejor represente el crecimiento de la especie. Con esta ecuación, se clasifican los niveles de productividad en bajo, medio y alto.


Área de estudio y base de datos

Se empleó una red de 129 parcelas permanentes de monitoreo instaladas durante el periodo 2008-2013. Las parcelas se establecieron en 5 municipios del departamento del Magdalena, ubicado en la región Caribe de Colombia, y en 4 municipios del departamento del Tolima, el cual se encuentra en la región Andina del centro-oeste del país (Figura 1). Se registraron mediciones dasométricas desde el primer año de establecimiento hasta los 16 años, con un total de 2 a 5 remediciones en cada parcela. En cada parcela y medición, se registró el diámetro a la altura del pecho (DAP), tomado a 1.30 m del suelo (DAP) con cinta diamétrica, y la altura total (HT), medida con un hipsómetro laser.

Distribución espacial de las unidades muestrales en los departamentos de Magdalena y Tolima

Figura 1: Distribución espacial de las unidades muestrales en los departamentos de Magdalena y Tolima

En el departamento del Magdalena, la altura dominante media fue de 23.28 m, con un rango de 6.9-30.28 m. En el departamento del Tolima, la altura dominante varió entre 5.72 y 36.0 m, con una media de 18.6 m (Tabla 1). La base de datos presenta mediciones a edades tempranas (mayores a un año del establecimiento de la plantación), intermedias (antes de la edad de corta) y mayores a la edad de aprovechamiento, lo que permite afirmar la idoneidad de los datos para el ajuste del modelo predictor.

Tabla 1: Resumen de las mediciones realizadas en parcelas permanentes de G. arborea en cada departamento y municipio

Altura dominante

Existen diferentes criterios para determinar la altura dominante de un rodal. Para este estudio, la altura dominante se considera como la altura media aritmética de los 100 árboles.ha-1 con los diámetros de mayor tamaño, tal como lo explica Johansson (2013). Una correcta predicción de la altura dominante basada en modelos de IS no sesgados y precisos es esencial para modelar el crecimiento y la producción (Vargas et al., 2013).

La Tabla 1 resume las características temporales y los rangos de altura dominante de la base de datos a nivel de departamento y municipio.

Modelos de altura dominante evaluados

Se usó un total de 9 funciones derivadas de 2 modelos de crecimiento (Hossfeld y Bertalanffy-Richards), los cuales han sido reportados por Burkhart & Tomé (2012), Panik (2014) y Tamarit-Urias et al. (2022). Se generaron 4 ecuaciones dinámicas del modelo Hossfeld y 5 ecuaciones del modelo de Bertalanffy-Richards, usando el EDA y el EDAG. En la Tabla 2 se presenta el modelo en su forma base, los parámetros seleccionados para el despeje de cada ecuación y la respectiva forma del modelo dinámico que se utiliza para modelar el IS.

Tabla 2: Modelos de crecimiento ajustados para la predicción del desarrollo en altura dominante

El EDA solo admite una hipótesis de crecimiento respecto a los parámetros, es decir, se tiene una sola asíntota común (modelos polimórficos) o se tienen asíntotas variables (modelos anamórficos), pero no es posible incluir ambas hipótesis en un modelo (Pacheco-Aguilar et al., 2016). Esta metodología, formalizada por Bailey & Clutter (1974), involucra esencialmente la sustitución de un parámetro del modelo base para expresarlo como una función del IS o de una condición inicial de la variable de interés. Estas ecuaciones poseen una asíntota horizontal y generan curvas anamórficas o polimórficas siguiendo un patrón de crecimiento sigmoide y un punto de inflexión correspondiente a un comportamiento esperado para cada especie (Talavera García, 2013).

El enfoque por diferencias algebraicas generalizadas (EDAG), desarrollado por Cieszewski & Bailey (2000), considera dos parámetros como dependientes de la calidad de sitio, lo que permite la generación de curvas polimórficas con múltiples asíntotas. Este procedimiento se considera como un método que produce estimaciones imparciales para ecuaciones invariantes a la edad base (Nunes et al., 2011). Esta es una propiedad deseable al momento de describir el crecimiento en altura dominante y clasificar la productividad de sitio (Pacheco-Aguilar et al., 2016).

Las ecuaciones dinámicas generadas son de la forma Y = f(t, Y 0, t 0, β i ), donde Y y Y 0 indican la altura dominante futura y la altura dominante de referencia respectivamente, t y t 0 son las edades respectivas para cada altura dominante y β i representa los parámetros a estimar.

Aquí, a, b y c son los parámetros propios en cada ecuación; b1, b2 y b3 son los parámetros de la ecuación dinámica; X 0 y L 0 son los parámetros dependientes de la calidad de sitio; ln es el logaritmo natural; y e es la función exponencial.

Ajuste y validación

Los parámetros de estos modelos se estimaron mediante el procedimiento Model en la interfaz del Statistical Analysis System (SAS) 9.4 (Statistical Analysis System, 2021), con el algoritmo de minimización de Marquardt y utilizando el método de máxima verosimilitud con información completa (FIML). Se empleó regresión no lineal ponderada asumiendo que la varianza de los residuales puede ser modelada como una función potencial de la altura dominante predicha w i = ŷ i m (Tamarit-Urias et al., 2022). La potencia (m) se determinó empleando la metodología descrita por Barrios et al. (2014). Adicionalmente, para evitar problemas de autocorrelación, los residuales fueron modelados a través de una estructura autorregresiva en tiempo continuo AR(1), empleando el macro %(ar) del procedimiento Model.

Se usó la técnica de validación cruzada para valorar la precisión de los modelos, la cual permite evaluar la generalización de los resultados empleando un conjunto independiente de datos (Refaeilzadeh et al., 2008). Para ello, los datos se dividieron aleatoriamente en 10 grupos de igual tamaño que contenían 13 parcelas (10 % de los datos totales). El procedimiento de validación cruzada (k-fold cross-validation) consistió en una división aleatoria en k número de grupos. Para este estudio, el total de subconjuntos fue 10; cada grupo contenía 13 parcelas (aproximadamente 10 % de los datos totales). Posteriormente, se realizaron k iteraciones (proceso de validación), de tal manera que, en cada iteración, se mantuvo un conjunto diferente de los datos, mientras que el resto (k-1) se utilizó para ajuste (Krkač et al., 2020). La bondad de ajuste se evaluó considerando el coeficiente de determinación ajustado (Burkhart & Tomé, 2012):

Como medidas de bondad de predicción, se utilizaron el sesgo promedio (Sesgo), el error absoluto medio (EAM) y la raíz del error medio cuadrático (REMC) (Barrios-Trilleras et al., 2021):

donde y i es la altura dominante observada, ŷ i es la altura dominante predicha, ȳ es la altura dominante media, n es el número de observaciones y p es el número de parámetros del modelo.

La evaluación de la autocorrelación se realizó con el análisis estadístico de Durbin-Watson (DW) (Panik, 2014), y los residuos de cada modelo se graficaron frente a los valores predichos para evaluar su distribución y heterocedasticidad (Barrios-Trilleras et al., 2021).

A partir de los análisis estadísticos de bondad de ajuste y la validación, se construyó un sistema de clasificación donde se establece un orden jerárquico de calidad y desempeño estadístico, en el cual se selecciona el modelo que obtenga la calificación más baja (Tewari & Singh, 2018), que finalmente se utiliza para predecir el IS de cada parcela, considerando una edad base igual a 10 años, usualmente empleada para determinar el IS de G. arborea (López et al., 2011; Martínez-Zurimendi et al., 2015).


Modelos ajustados de crecimiento de altura dominante

La Tabla 3 presenta el valor de los parámetros estimados, la bondad de ajuste y la bondad de predicción para los modelos evaluados, así como la calificación total de cada uno. En términos generales, los modelos restantes evaluados presentaron buenos ajustes. Todos los modelos presentaron R2 aj ≥ 0.914, a excepción de los modelos 3 y 4, los cuales presentaron R2 aj = 0.88.

Tabla 3: Parámetros estimados y bondad de ajuste y predicción para los modelos de crecimiento de altura dominante evaluados

En paréntesis: el intervalo de confianza para cada parámetro estimado.

Los datos estadísticos de bondad de predicción se consideran aceptables. Los valores de sesgo indican que los modelos tendieron a subestimar los valores de altura dominante con sesgos entre 0.009 y 0.651 m. Los modelos presentaron EAM entre 1.343 y 1.582 m, y el REMC presentó valores entre 1.746 y 2.099 m.

De acuerdo con los resultados del sistema de clasificación implementado, el modelo 7 (Bertalanffy-Richards), el modelo 2 y el modelo 4 sobresalieron con los mejores datos estadísticos. En relación con el análisis estadístico DW, los valores oscilaron entre 1.996 y 2.223, en general fueron próximos a 2 y, de acuerdo con Trim et al. (2020), se evidencia la inexistencia de autocorrelación.

La Figura 2 muestra el comportamiento de los residuales del modelo 7. Como se aprecia, no hay tendencias anómalas, indicando la ausencia de heterocedasticidad. Al graficar los valores observados frente a los predichos, se observa que el modelo 7 no presenta tendencias significativas de sobrestimación o subestimación, lo que lo vuelve una ecuación idónea para la proyección del IS, facilitando la construcción de planes precisos para la planificación de la especie.

Dispersión de residuales y valores observados frente a valores predichos para el modelo 7 (modelo de Richards, parámetro c).

Figura 2: Dispersión de residuales y valores observados frente a valores predichos para el modelo 7 (modelo de Richards, parámetro c).

Finalmente, se analiza la familia de curvas generadas por la ecuación de Richards. Se presentan curvas ajustadas a IS con edad de referencia de 10 años superpuestas a las trayectorias de las parcelas observadas a través de la edad. De esta manera, se buscó evaluar el comportamiento biológico de los modelos, representado en una familia de tres curvas generadas para clasificar la productividad de los sitios como alto, medio y bajo, que corresponden al percentil 95, la mediana y el percentil 5 de los datos observados respectivamente.

Las curvas de IS generadas por los modelos presentan diferentes tasas de crecimiento a edades tempranas (de 0 a 5 años). En edades medias o mayores de 5 años, se ve un crecimiento continuo hacia la formación de una asíntota horizontal, lo cual indica un punto de inflexión en los niveles de desarrollo de la especie. Esta información generalmente se utiliza para la planificación de prácticas silvícolas y la proyección de la productividad.

Para el modelo de Richards, parámetro c (modelo 7), se observa en la Figura 3 un ajuste claro para las tres curvas evaluadas, donde, a edades tempranas, se presenta un crecimiento acelerado, marcado por una deflexión donde la tasa de crecimiento se reduce en el tiempo. Esto es bastante acertado para un comportamiento biológico promedio. Este modelo mostró un rango de IS entre los 15.5 y los 26.4 m a una edad base de 10 años. Para clasificar la productividad de plantaciones de Gmelina arborea, se generaron tres rangos basados en las curvas de la Figura 3: sitios con bajas productividades, o aquellos con IS ≤ 18.6 m; sitios con productividades medias entre 18.6 m < IS ≤ 24.1 m; y sitios con productividad alta, con IS > 24.1 m.

Familia de curvas de índice de sitio para G. arborea obtenidas a partir del modelo de Richards, parámetro c (modelo 7)

Figura 3: Familia de curvas de índice de sitio para G. arborea obtenidas a partir del modelo de Richards, parámetro c (modelo 7)

Finalmente, se contrastan los valores de IS para el modelo 7 para cada uno de los municipios muestreados en ambos departamentos (Tabla 4).

Tabla 4: Índice de sitio promedio para plantaciones de G. arborea por municipio y departamento

Como se evidencia en la Tabla 4, los departamentos presentan el mismo nivel de productividad. Todos los municipios fueron categorizados como de productividad media, excepto el municipio de Coyaima en el departamento del Tolima, que presenta un valor más bajo de productividad. La Ecuación (5) de IS para G. arborea según el modelo 7 se muestra a continuación:

donde t0 hace referencia a la edad base o edad clave (10 años) y H y t son la altura dominante y edad del rodal.


Este estudio presenta una evaluación de modelos de crecimiento en altura dominante para Gmelina arborea, los cuales son de utilidad para clasificar la productividad de plantaciones de la especie. A partir de dos funciones de crecimiento base, se generaron 9 modelos dinámicos de altura dominante empleando las metodologías EDA y EDAG. Estudios previos han demostrado que el uso de estos enfoques ofrece estimaciones de altura dominante sustancialmente más precisas en comparación con el enfoque de curva guía debido a su naturaleza polimórfica (Scolforo et al., 2020).

Los más destacados fueron los modelos 7, 8 y 2, los dos primeros originados a partir del modelo de Bertalanffy-Richards y el tercero a partir del modelo de Hossfeld. Estudios previos en Colombia (López et al., 2011; Barrios et al., 2021) habían coincidido en que el modelo de crecimiento base de Bertalanffy-Richards representa mejor el crecimiento en altura dominante de Gmelina arborea. Sin embargo, en un estudio realizado para plantaciones forestales de Gmelina arborea en Tabasco, México, se seleccionó el modelo base de Hossfeld para describir el crecimiento de altura dominante (Martínez-Zurimendi et al., 2015). Estos resultados reflejan la complejidad del modelamiento de la calidad de sitio, debido a que el modelo refleja el desarrollo local o regional en altura dominante de los árboles.

En este estudio se selecciona el modelo dinámico 7, el cual se origina a partir del modelo base de Bertalanffy- Richards, dejando el parámetro c dependiente del sitio. El modelo de crecimiento seleccionado para esta especie muestra las características típicas de un modelo de crecimiento biológico, que, según Scolforo et al. (2013) inicia en el punto cero, presenta crecimiento acelerado a edades juveniles, tiene un punto de inflexión donde el crecimiento se desacelera e incluye una asíntota que refleja el valor máximo que puede alcanzar el organismo. Los resultados obtenidos concuerdan con los de López et al. (2011) en cuanto a que el desarrollo de altura dominante de Gmelina arborea se caracteriza por tener diferentes tasas de crecimiento, lo que da lugar a polimorfismos. No obstante, los rodales finalmente alcanzan una altura dominante máxima (asíntota).

La productividad de G. arborea fue categorizada en tres rangos: sitios con bajas productividades con IS ≤ a 18.6 m, sitios con productividades medias con IS entre 18.6 y 24.1 m y sitios con productividad alta con IS > 24.1 m. Estos intervalos se encuentran en los rangos propuestos por diversos autores que evaluaron la productividad de Gmelina arborea (Hughell, 1991; Vásquez & Ugalde, 1995; Cifuentes, 2010; Martínez- Zurimendi et al., 2015). En todos los casos se empleó una edad de referencia de 10 años. En términos generales, la productividad de las plantaciones forestales de G. arborea evaluadas puede ser considerada como media. Sin embargo, es factible que esta incremente a causa del desarrollo genético y la silvicultura. En todo caso, el modelo desarrollado tiene la capacidad de adaptarse a estos cambios y mantener su capacidad predictiva.


Este estudio evaluó una amplia gama de modelos de proyección de altura dominante, lo que permitió disminuir el error de selección, reduciendo considerablemente las posibilidades de sobreestimar o subestimar la calidad de sitio, pues esto afectaría directamente los supuestos de planificación y liquidez durante el desarrollo de la plantación.

El modelo que mejor se ajustó a los datos y al comportamiento de la especie fue el número 7 (Bertalanffy- Richards, parámetro c), presentando una tendencia biológica acorde con tasas de crecimiento acelerado en edades juveniles y un punto de inflexión donde el crecimiento se desacelera. También e incluye una asíntota que refleja la altura dominante máxima que puede alcanzar el rodal. Por lo tanto, se considera Gmelina arborea como una especie que se ajusta mejor a modelos que generen familias de curvas polimórficas y que tiendan a alcanzar una única asíntota en edades avanzadas. El modelo seleccionado es un modelo integral y funcional que puede ser utilizado para evaluar la calidad de sitio en los departamentos de Magdalena y Tolima.



Los autores agradecen a la Oficina de Investigación y Desarrollo Científico de la Universidad del Tolima por el apoyo financiero otorgado para desarrollar la investigación, a través del proyecto no. 340130516.


Aguirre-Mendoza Z, Aguirre-Mendoza C. Especies leñosas y cultivos objetivos para sistemas agroforestales en zonas semiáridas del sur del Ecuador. Bosques. Latitud Cero, 4, 21-30. [Link]

Akbas U, Senyurt M. Site quality estimations based on the generalized algebraic difference approach a case study in Çankiri forests. Revista Árvore, 42(3), e420311. [Link]

Bailey R, Clutter JL. Base-age invariant polymorphic site curves. Forest Science, 20(2), 155-159.

Barrios-Trilleras A, López-Aguirre AM, Báez-Aparicio CA. Modelamiento de la productividad de Gmelina arborea Roxb con base en variables biofisicas y de rodal. Colombia Forestal, 24(1), 71-87. [Link]

Barrios A, López A.M, Nieto V. Predicción de volúmenes comerciales de Eucalyptus grandis a través de modelos de volumen total y de razón. Colombia Forestal, 17(2), 137-149. [Link]

Burkhart HE, Tomé M. Modeling forest trees and stands. Springer Science & Business Media; 2012. [Link]

Cieszewski CJ, Bailey RL. Generalized algebraic difference approach: Theory-based derivation of dynamic site equations with polymorphism and variable asymptotes. Forest Science, 2000;461116-126.

Cieszewski, Harrison, Martin. Practical methods for estimating non-biased parameters in self-referencing growth and yield models.[Link]

Cifuentes JG. Evaluación y monitoreo de plantaciones forestales en Guatemala. Instituto Nacional de Bosques. 2010.

Clutter JL, Forston JC, Pienaar LV, Brister GH, Bailey RL. Timber Management: A quantitative approach. John Wiley & Sons. 1983.

Escobar Quemba L. Relacion de parámetros de fertilidad del suelo con el Índice de sitio determinado para plantaciones forestales de melina (Gmelina arbórea) y ceiba (Pachira quinata) en Zambrano-Bolivar (Colombia). Tesis de maestría, Universidad Nacional de Colombia.[Link]

Espitia M, Murillo O, Castillo C. Ganancia genética esperada en melina (Gmelina arborea Roxb.) en Córdoba (Colombia). Revista Árvore, 2016;40171-80.[Link]

Gómez EA, Rios LA, Peña JD. Madera, un potencial material lignocelulósico para la producción de biocombustibles en Colombia. Informacion Tecnológica, 2012;23673-86.[Link]

Guerra-Hernández J, Arellano-Pérez S, Gonzáles-Ferreiro E, Pascual A, Altelarrea VS, Ruiz-González AD, Álvarez-Gonzáles JG. Developing a site index model for Pinus pinaster stands in NW Spain by combining bi-temporal ALS data and environmental data. Forest Ecology and Management, 481, 118690. 2021.[Link]

Hernández-Cuevas M, Santiago-Garcia W, De los Santos-Posadas HM, Martinez-Antúnez P, Ruiz-Aquino F. Modelos de crecimiento en altura dominante e índices de sitio para Pinus ayacahuite Ehren. Agrociencia, 2018;523437-453.[Link]

Hughell D. Modelo preliminar para la predicción del rendimiento de Gmelina arborea Roxb. en América Central. CATIE. 1991.

Johansson T. A site dependent top height growth model for hybrid aspen. Journal of Forestry Research, 24, 691-698. 2013.[Link]

Krkač M, Bernat-Gazibara S, Arbanas Ž, Sečanj M, Mihalić-Arbanas S. A comparative study of random forests and multiple linear regression in the prediction of landslide velocity. Landslides, 17, 2515-5131. 2020.[Link]

López AM, Barrios A, Nieto V. Monitoreo y modelamiento del crecimiento para el manejo de plantaciones forestales comerciales. Intersolugráficas. 2011.

Martinez-Zurimendi P, Dominguez-Dominguez M, Juárez-Garcia A, López-López LM, De la Cruz V, Álvarez J. Índice de sitio y producción maderable en plantaciones forestales de Gmelina arborea en Tabasco. Revista Fitotecnia Mexicana, 38, 415-425. 2015.[Link]

Melo Cruz OA. Modelación del crecimiento, acumulación de biomasa y captura de carbono en árboles de Gmelina arborea Roxb., asociados a sistemas agroforestales y plantaciones homogéneas en Colombia. Tesis de doctorado, Universidad Nacional de Colombia.[Link]

Ministerio de Agricultura y Desarrollo Rural (MADR). Boletín estadístico forestal. 2022.

Nunes L, Patricio M, Tomé J, Tomé M. Modeling dominant height growth of maritime pine in Portugal using GADA methodology with parameters depending on soil and climate variables. Annals of Forest Science, 68, 311-323. 2011.[Link]

Nurudeen TA, Odeleye S, Akinyemi OD, Salami KD, Lawal AA. Index age for estimating site index of Gmelina arborea (Roxb.) stands in ioluwa forest reserve, South Western Nigeria. Fudma Journal of Sciences, 2, 154-158. 2018.[Link]

Pacheco-Aguilar G, Santiago-Juárez W, Martínez-Santiago D, Ortiz-Barrios R. Analisis del crecimiento e incremento y estimacion de índice de sitio para Pinus montezumae Lamb. en Santiago Textitlán, Sola de Vega, Oaxaca. Foresta Veracruzana, 2, 21-28. 2016.[Link]

Panik MJ. Growth curve modeling: Theory and applications. John Wiley & Sons. 2014.[Link]

Quiñonez Barraza G, De los Santos Posadas HM, Cruz Cobos F, Velázquez Martínez A, Ángeles Pérez G, Ramírez Valverde G. Índice de sitio con polimorfismos complejos para masas forestales en Durango, México. Agrociencia, 49, 439-454. 2015.

Refaeilzadeh P, Tang L, Lui H. Cross Validation. Liu L, Özsu MT. Encyclopedia of Database Systems, 677-684. Springer. 2008.[Link]

Richards FJ. A flexible growth function for empirical use. Journal of Experimental Botany, 10, 290-301. 1959.[Link]

Sabatia CO, Burkhart HE. Predicting site index of plantation lobolly pine from biophysical variables. Forest Ecology and Management, 326, 142-156. 2014.[Link]

Scolforo HF, McTague J, Burkhart H, Roise J, Alvares CA, Stape JL. Site index estimation for clonal eucalypt plantations in Brazil: A modeling approach refined by environmental variables. Forest Ecology and Management, 466, 118079. 2020.[Link]

Scolforo J, Maestri R, Ferraz A, Mello J, Oliveira A, Assis A. Dominant height model for site classification of Eucalyptus grandis incorporating climatic variables. International Journal of Forestry Research, 2013, 139236. 2013.[Link]

Statistical Analysis System. Statistical Analysis System-SAS. 2021.[Link]

Talavera García J. Modelización de crecimiento de Pinus pinaster. Tesis de doctorado. Universidad Autónoma de Nuevo León.[Link]

Tamarit-Urias JC, Rodríguez-Acosta M, Quiñonez-Barraza G, Santiago-Garcia W. New dynamic equation of dominant height and site index for Pinus chiapensis (Martínez) Andresen in Puebla and Veracruz, Mexico. Revista Chapingo Serie Ciencias Forestales y del Ambiente, III(28), 331-347. 2022.[Link]

Tewari VP, Singh B. Total wood volume equation for Tectona grandis linn F. stands in Gujarat, India. Journal of Forest and Environmental Science, 34(4), 313-320. 2018.[Link]

Trim KR, Coble DW, Weng Y, Stovall JP, Hung IK. A new site index model for intensively managed loblolly pine (Pinus taeda) plantations in the west Gulf Coastal Plain. Forest Science, 66(1), 2-13. 2020.[Link]

Vanclay JK, Baynes J, Cedamon E. Site index equation for smallholder plantations of Gmelina arborea in Leyte Province. Small-Scale Forestry, 7, 87-93. 2008.[Link]

Vargas B, Aguirre OA, Corral JJ, Crecente F, Diéguez U. Modelo de crecimiento en altura dominante e índice de sitio para Pinus pseudostrobus Lindl. en el noreste de México. Agrociencia, 47(1), 91-106. 2013.[Link]

Vásquez W, Ugalde L. Rendimiento y calidad de sitio para Gmelina arborea, Tectona grandis, Bombacopsis quinatum y Pinus caribaea en Guanacaste, Costa Rica. CATIE. 1995

Collazos-Avendaño, J. C., Barrios-Trilleras, A., & Díaz-Sierra, D. R. (2024). Modelos de altura dominante para Gmelina arborea Roxb. en Colombia. Colombia Forestal, 27(1), e20569