Forest Structure and Threatened Tree Species along Elevation and Disturbance Gradients in a Tropical River Watershed

Ruffy Rodrigo

Department of Forest Science, Biliran Province State University, Biliran, Philippines

Corresponding Author Email: rodrigoruffy@gmail.com

DOI : https://doi.org/10.51470/JOD.2025.4.2.83

Abstract

Understanding the relationship between forest structure, elevation, and species conservation status is vital for developing effective biodiversity and carbon management strategies. The current study examined 27 forest plots within the selected River Watershed in the Philippines to assess how disturbance and elevation influence tree diameter, height, basal area, aboveground biomass, and the occurrence of threatened species. Structural attributes varied widely across plots, with mean tree diameter ranging from 13.9 to 29.2 cm, and biomass spanning from approximately 3,100 kg to over 30,600 kg per plot. Undisturbed plots consistently exhibited higher basal area and biomass, reflecting greater structural integrity and ecological resilience. Conservation status data revealed that plots hosted up to six IUCN-listed threatened species and as many as 15 DAO 2017-11 classified species, with elevated and undisturbed plots supporting the greatest numbers. Regression models confirmed that elevation significantly influenced the richness of IUCN-listed species, while disturbance strongly predicted reductions in both biomass and basal area. Plots above 700 meters above sea level emerged as conservation hotspots, emphasizing the value of upland forests as refugia for threatened taxa.
Visual analyses reinforced these findings: undisturbed plots had double the biomass and nearly twice the basal area compared to disturbed counterparts. These results align with regional and global literature, highlighting the urgency of protecting intact forest landscapes. By linking forest structure with conservation indicators, the study provides empirical evidence to inform sustainable forest management and biodiversity conservation in tropical ecosystems under increasing human pressure.

Keywords

biomass, conservation, disturbance, forest structure, tropical river watershed

Download this article as:

  1. Introduction

Forest ecosystems in tropical nations, such as the Philippines, are pivotal for providing essential ecological services, including biodiversity conservation, climate regulation, and carbon sequestration. However, these invaluable ecosystems face significant threats from land conversion, logging, and various human disturbances that jeopardize their integrity [1, 2]. Understanding the interplay between forest structure and conservation status is vital for developing effective forest management and protection strategies. Forest structure, which is characterized by metrics such as tree diameter, height, and basal area, serves as a useful proxy for assessing ecological integrity and habitat quality [3].  Further, complex forest structures typically support greater biomass and provide critical niches for diverse species, including those threatened with extinction [4].  However, there remains a paucity of studies directly linking forest structure with the presence of threatened species across varying disturbance gradients, particularly within the Philippine context [5].

The Amambahag River Watershed, located in Biliran Island, exemplifies a crucial ecological corridor that sustains biodiversity and underpins local livelihoods [6]. Unfortunately, parts of this landscape are undergoing increasing disturbances, which threaten its conservation potential. This present study seeks to investigate how disturbance and elevation modulate forest structure, biomass, and the richness of conservation-listed species. By integrating field-based forest inventories with species-level conservation assessments, the research aimed to offer empirical evidence regarding the structural and biodiversity value of undisturbed forests [7]. Such research is crucial for assessing forest condition in relation to conservation importance at a local scale and will support informed decision-making focused on forest protection and sustainable management practices [8].  

Hence, this research will address the interactions between elevation, disturbance, and forest metrics, not only enhancing our understanding of tropical forest ecosystems but is essential for calibrating conservation efforts within ecologically sensitive areas like the Amambahag River Watershed. Through this investigation, the findings will contribute to the broader body of knowledge that informs forest conservation strategies, which are imperative for safeguarding the rich biodiversity and ecosystem services that these forests provide [9].

  1. Materials and Methods

3.1 Study Area

The study was situated in the Amambahag River Watershed on Biliran Island, located in the Eastern Visayas region of the Philippines (Figure 1). This watershed encompasses a diverse range of montane and submontane forest ecosystems recognized for their rich biodiversity and ecological significance, particularly regarding flora and fauna typical of tropical mountainous regions [10]. Elevations within the watershed vary from approximately 500 to 900 meters above sea level (masl), resulting in a range of microclimates. The forest areas surveyed included both disturbed sites, characterized by logging or human access, and undisturbed patches, which provide critical habitats for various species. Such a mixed environment facilitates the evaluation of how anthropogenic pressures and elevation influence forest structure and conservation value, as human disturbances can significantly alter both biodiversity and ecosystem function [11].  

3.2 Data Collection

A 2-km transect served as the baseline, with 20 m x 20 m quadrats established every 250 meters along the transect. A total of three transects was laid out, with a total of nine plots per transect. In each plot, adult trees with a diameter equal to or greater than 10 centimeters had measurements of diameters at breast height (DBH), merchantable height (MH), and total height (TH) were recorded (Figure 6). The approach complemented the Biodiversity Monitoring System (BMS) transects currently used by the Biodiversity Management Bureau (BMB) for faunal assessments. Each plot was placed to ensure adequate representation of both disturbed and undisturbed forest types, allowing for a comprehensive understanding of forest dynamics.

Further, key data collected for each tree included the DBH, total height, and species identity. To inform conservation assessments, the conservation status of each tree species was documented by consulting the International Union for Conservation of Nature (IUCN) Red List and the Philippine Department of Agriculture’s DAO 2017-11 classification, which provides vital information regarding the risk of extinction and necessary conservation measures for the species [12].  

3.3 Data Analysis

The analysis of forest metrics involved calculating the basal area (BA) using the formula:  BA = pi (3.1416) x DBH/2)2, which provides a standard measure of tree stand density at the plot level (Estavillo et al., 2013).  The biomass of trees (measured in kilograms) was estimated using a simplified allometric equation that incorporates DBH, height, and wood density; this estimation facilitates comparability across plots and aids in understanding carbon storage potential [13]. Data from each plot were summarized by averaging DBH, mean tree height, total basal area, total biomass, and the number of species designated as threatened under both the IUCN and DAO classifications.

In addition, to quantitatively assess the influence of disturbance and elevation on biomass, basal area, and species richness, linear models were employed. Furthermore, pie charts were generated to visually represent the distributions of conservation statuses across the recorded species, enhancing the interpretability of the data. All statistical analyses and visualizations were conducted using R statistical software, which is well-suited for ecological data analysis and graphical representation, allowing for robust conclusions regarding the ecological significance of the studied forest areas [14].  

  • Results and Discussion

The forest structure and conservation value across the 27 analyzed plots in the Amambahag River Watershed exhibit a noteworthy variability that reflects both ecological and anthropogenic influences. The mean tree diameter at breast height (DBH) ranged from 13.9 to 29.2 cm with an overall average of 20.6 cm, while the mean tree height was approximately 9.15 meters, indicative of a relatively mature forest structure [15]. Basal area measurements spanned from 0.22 m² to 1.88 m² per plot (Table 1; Figure 2), supporting the notion of substantial ecological differences within the watershed [16]. Biomass estimates showed considerable variations, from about 3,125 kg to over 30,650 kg per plot. Importantly, undisturbed forest areas demonstrated higher basal area and biomass values due to reduced human impact, corroborating findings of studies that highlight intact forests as critical for biodiversity and carbon storage [17].

Moreover, conservation metrics reveal considerable variation in the number of IUCN-listed threatened species per plot, ranging from 0 to 6, with a mean of 2.63 [17]. The Philippine DAO 2017-11 classification, the plots supported up to 15 species classified as threatened or vulnerable, averaging 5.33 species per plot, aligning with patterns observed in other forest types where biodiversity declines with disturbance levels [18]. Higher elevation plots, peaking at 932 masl, reinforced the ecological importance of undisturbed mountainous forests, as they generally exhibited greater representation of conservation-relevant species. Such findings resonate with previous literature highlighting the correlation between elevation and biodiversity conservation potential [19, 20].  

The multiple linear regression analyses conducted in this study elucidate significant relationships between forest disturbance, elevation, and key ecological metrics (Table 3). Notably, elevation positively influenced the number of IUCN-listed threatened species (Estimate = 0.00, 95% CI: 0.00–0.01, p = 0.022), signifying that higher elevations correlate with greater conservation value [21].  While the effect size was minimal due to the scale of elevation, the statistical significance highlights its ecological implications. Disturbance conditions emerged as strong predictors of total biomass and basal area. Specifically, undisturbed plots outperformed disturbed plots by an average of 7,522 kg in biomass (95% CI: 2,921–12,123 kg, p = 0.002) and 0.49 m² in basal area (95% CI: 0.22–0.77 m², p = 0.001), echoing findings from regional studies that connect forest disturbance with decreased structural integrity [22].

Illustrative evidence from boxplots indicates a marked biomass disparity between disturbed and undisturbed plots (Figure 3). The undisturbed forests frequently surpassed 30,000 kg in biomass, while disturbed areas typically remained below 15,000 kg [23]. This pattern emphasizes the significance of preserved forests for biomass storage, crucial for carbon sequestration and habitat maintenance, while also highlighting the detrimental impact of disturbance on forest structure. Basal area metrics further reinforced this argument, with undisturbed plots boasting a mean of approximately 1.1 m² compared to 0.6 m² in disturbed plots, highlighting degradation in structural capacity due to human interventions [24, 25].

Lastly, the positive correlation between elevation and the number of IUCN-listed threatened species per plot (Figure 4), observed through scatter plot analyses, affirms the role of montane forests as critical refugia for species of conservation concern within the Amambahag River Watershed. Higher elevation zones, particularly over 700 masl, showed increased counts of threatened species in less disturbed contexts, which aligns with ecological patterns documented in various tropical and subtropical forest studies [26, 27].

Multiple linear regression models revealed important relationships between forest disturbance, elevation, and key ecological metrics. Elevation was significantly associated with the number of IUCN-listed threatened species (Estimate = 0.00, 95% CI: 0.00–0.01, p = 0.022), suggesting that higher elevation plots support greater conservation value. While the effect size is small due to the unit scale of elevation (masl), the trend remains statistically meaningful.

Disturbance condition was a significant predictor of both total biomass and basal area. Undisturbed plots had, on average, 7,522 kg more biomass (95% CI: 2,921–12,123 kg, p = 0.002) and 0.49 m² more basal area (95% CI: 0.22–0.77 m², p = 0.001) compared to disturbed plots. These findings highlight the structural degradation associated with disturbance. While disturbance had a positive effect on IUCN species count (Estimate = 1.03), it was not statistically significant (p = 0.135). The models explained 28–36% of the variance in the response variables, with the basal area model showing the highest adjusted R² (0.333).

  • Conclusion

The study highlights the critical role of undisturbed and high-elevation forest areas in maintaining structural integrity and conservation value within the Amambahag River Watershed. Plots with minimal human disturbance exhibited significantly higher biomass and basal area, reinforcing the ecological function of intact forests as key carbon sinks and biodiversity reservoirs. Elevation also emerged as a strong ecological determinant, with higher-altitude plots supporting more IUCN-listed threatened species, particularly in less disturbed settings. These findings highlight the compounded effect of disturbance and elevation on forest health and conservation outcomes. Maintaining and protecting undisturbed upland forests is essential not only for sustaining carbon storage but also for safeguarding species of conservation concern. The integration of forest structure and species status data provides valuable insights for formulating targeted conservation strategies in tropical landscapes under increasing anthropogenic pressure.

  1. Acknowledgment

The author gratefully acknowledges the Provincial Environment and Natural Resources Office (PENRO) in Naval, Biliran, for providing access to essential field data. Special thanks are extended to Forester Jufamar Fernandez, MSc. for the support and logistical assistance. Appreciation is also given to Biliran Province State University for the institutional support and for providing internal funding assistance for the publication of this study. Their contributions were instrumental in the successful completion and dissemination of this research.

  • References

1 Pendrill, F., Persson, U., Godar, J., & Kästner, T. (2019). Deforestation displaced: Trade in forest-risk commodities and the prospects for a global forest transition. Environmental Research Letters, 14(5), 055003. https://doi.org/10.1088/1748-9326/ab0d41

2 Busch, J., & Engelmann, J. (2017). Cost-effectiveness of reducing emissions from tropical deforestation, 2016–2050. Environmental Research Letters, 13(1), 015001. https://doi.org/10.1088/1748-9326/aaa4c4

3 Igu, N. (2017). Species diversity and structure of an intact freshwater swamp forest in the Niger Delta. Open Journal of Forestry, 7(2), 242–254. https://doi.org/10.4236/ojf.2017.72014

4 Ferreira, J., Lennox, G., Gardner, T., Thomson, J., Berenguer, É., Lees, A., & Barlow, J. (2018). Carbon-focused conservation may fail to protect the most biodiverse tropical forests. Nature Climate Change, 8(8), 744–749. https://doi.org/10.1038/s41558-018-0225-7

6 Hansen, M., Wang, L., Song, X., Tyukavina, A., Turubanova, S., Potapov, P., & Stehman, S. (2020). The fate of tropical forest fragments. Science Advances, 6(11), eaax8574. https://doi.org/10.1126/sciadv.aax8574

7 Brodie, J., Williams, S., & Garner, B. (2021). The decline of mammal functional and evolutionary diversity worldwide. Proceedings of the National Academy of Sciences, 118(3), e1921849118. https://doi.org/10.1073/pnas.1921849118

8 Taylor, P., Asner, G., Dahlin, K., Anderson, C., Knapp, D., Martin, R., & Townsend, A. (2015). Landscape-scale controls on aboveground forest carbon stocks on the Osa Peninsula, Costa Rica. PLOS ONE, 10(6), e0126748. https://doi.org/10.1371/journal.pone.0126748

9 Liang, J., Crowther, T., Picard, N., Wiser, S., Zhou, M., Alberti, G., & Reich, P. (2016). Positive biodiversity–productivity relationship predominant in global forests. Science, 354(6309), aaf8957. https://doi.org/10.1126/science.aaf8957

10 Española, C., Collar, N., Mallari, N., & Marsden, S. (2016). Large avian frugivores in the Philippines show linear responses to improvements in forest quality. Forest Ecology and Management, 375, 127–133. https://doi.org/10.1016/j.foreco.2016.05.022

11 Edwards, D., Magrach, A., Woodcock, P., Ji, Y., Lim, N., Edwards, F., & Yu, D. (2014). Selective-logging and oil palm: Multitaxon impacts, biodiversity indicators, and trade-offs for conservation planning. Ecological Applications, 24(8), 2029–2049. https://doi.org/10.1890/14-0019.1

12 Kittelberger, K., Neate‐Clegg, M., Blount, J., Posa, M., McLaughlin, J., & Şekercioğlu, Ç. (2021). Biological correlates of extinction risk in resident Philippine avifauna. Frontiers in Ecology and Evolution, 9, 618441. https://doi.org/10.3389/fevo.2021.618441

13 Igu, N., Ezenwenyi, J., Ayogu, C., & Okolo, N. (2023). Carbon storage and environmental determinants in a tropical rainforest landscape. Open Journal of Ecology, 13(4), 229–240. https://doi.org/10.4236/oje.2023.134014

14 Mukul, S., Herbohn, J., & Firn, J. (2016). Co-benefits of biodiversity and carbon sequestration from regenerating secondary forests in the Philippine uplands: Implications for forest landscape restoration. Biotropica, 48(6), 882–889. https://doi.org/10.1111/btp.12391

15 Osuri, A., Machado, S., Ratnam, J., Sankaran, M., Ayyappan, N., Muthuramkumar, S., & Naeem, S. (2020). Tree diversity and carbon storage cobenefits in tropical human-dominated landscapes. Conservation Letters, 13(2), e12704. https://doi.org/10.1111/conl.12704

16 Holdaway, R., Easdale, T., Carswell, B., Richardson, S., Peltzer, D., Mason, N., & Coomes, D. (2016). Nationally representative plot network reveals contrasting drivers of net biomass change in secondary and old-growth forests. Ecosystems, 20(5), 944–959. https://doi.org/10.1007/s10021-016-0072-0

17 Magnago, L., Magrach, A., Laurance, W., Martins, S., Meira-Neto, J., Simonelli, M., & Edwards, D. (2015). Would protecting tropical forest fragments provide carbon and biodiversity cobenefits under REDD+? Global Change Biology, 21(9), 3455–3468. https://doi.org/10.1111/gcb.12937

18 Baraloto, C., Rabaud, S., Molto, Q., Blanc, L., Fortunel, C., Hérault, B., & Fine, P. (2011). Disentangling stand and environmental correlates of aboveground biomass in Amazonian forests. Global Change Biology, 17(8), 2677–2688. https://doi.org/10.1111/j.1365-2486.2011.02432.x

19 Qi, G., Peng, S., Chen, C., Zhao, G., Liang, Y., Li, P., & Wen, Z. (2013). Biomass carbon storage of three forests at different altitudes in Baotianman Mountain, Central China. Advanced Materials Research, 864–867, 2463–2468. https://doi.org/10.4028/www.scientific.net/AMR.864-867.2463

20 Zhou, L., Dai, L., Wang, S., Huang, X., Wang, X., Qi, L., & Shao, G. (2011). Changes in carbon density for three old-growth forests on Changbai Mountain, Northeast China: 1981–2010. Annals of Forest Science, 68(5), 953–958. https://doi.org/10.1007/s13595-011-0097-1

21 Saatchi, S., Harris, N., Brown, S., Lefsky, M., Mitchard, E., Salas, W., & Morel, A. (2011). Benchmark map of forest carbon stocks in tropical regions across three continents. Proceedings of the National Academy of Sciences, 108(24), 9899–9904. https://doi.org/10.1073/pnas.1019576108

22 Hanberry, B., & He, H. (2015). Effects of historical and current disturbance on forest biomass in Minnesota. Landscape Ecology, 30(8), 1473–1482. https://doi.org/10.1007/s10980-015-0190-y

23 Munteanu, C., Senf, C., Niță, M., Sabatini, F., Oeser, J., Seidl, R., & Kuemmerle, T. (2021). Using historical spy satellite photographs and recent remote sensing data to identify high-conservation-value forests. Conservation Biology, 36(2), e13724. https://doi.org/10.1111/cobi.13724

24 Holdaway, R., Easdale, T., Carswell, B., Richardson, S., Peltzer, D., Mason, N., & Coomes, D. (2016). Nationally representative plot network reveals contrasting drivers of net biomass change in secondary and old-growth forests. Ecosystems, 20(5), 944–959. https://doi.org/10.1007/s10021-016-0072-0

25 Estavillo, C., Pardini, R., & Rocha, P. (2013). Forest loss and the biodiversity threshold: An evaluation considering species habitat requirements and the use of matrix habitats. PLOS ONE, 8(12), e82369. https://doi.org/10.1371/journal.pone.0082369

26 Padgurschi, M., Vieira, S., Stefani, E., Nardoto, G., & Joly, C. (2018). Bamboos and a new perspective on nitrogen input in tropical forests. PeerJ Preprints. https://doi.org/10.7287/peerj.preprints.26458

27 Sutton, L., Ibañez, J., Salvador, D., Taraya, R., Opiso, G., Senarillos, T., & McClure, C. (2023). Priority conservation areas and a global population estimate for the critically endangered Philippine eagle. Animal Conservation, 26(5), 684–700. https://doi.org/10.1111/acv.12900