One Ecosystem :
Research Article
|
Corresponding author: Wilbert A. Aureo (wilbert.aureo@bisu.edu.ph)
Academic editor: Stoyan Nedkov
Received: 15 Jul 2020 | Accepted: 11 Dec 2020 | Published: 11 May 2021
© 2021 Wilbert Aureo, Tomas Reyes Jr., Francis Carlo Mutia, Danilo Tandang, Reizl Jose
This is an open access article distributed under the terms of the Creative Commons Attribution License (CC BY 4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Citation:
Aureo WA, Reyes Jr. TD, Mutia FCU, Tandang DN, Jose RP (2021) Floristic composition and community structure along the elevational gradient of Balinsasayao Twin Lakes Natural Park in Negros Oriental, Philippines. One Ecosystem 6: e56536. https://doi.org/10.3897/oneeco.5.e56536
|
|
Balinsasayao Twin Lakes Natural Park (BTLNP) is one of the protected areas on the Island of Negros Oriental which is enormously rich in biodiversity due to different Lowland types formed along its elevation gradient. This study was conducted to better understand the composition and diversity of plant species in the natural park to improve conservation and management efforts of these remaining forests which are currently under threat from eco-tourism and other anthropomorphic influences. Within the 18 randomly distributed nested plots, a total of 351 species of plants were recorded. Of these, 183 species were trees, 54 herbs, 51 shrubs, 41 pteridophytes and 22 vines. The result of hierarchical cluster analysis showed differences in plant composition along the elevation gradient. There were 30 (9%) threatened species (vulnerable, endangered and critically endangered) and most of which were found in the submontane Lowland. Furthermore, the species diversity increases from lowland to sub-montane and eventually decreases towards montane Lowland. These results not only indicate the importance of BTLNP, but also highlights the submontane as a special area of concern due to the higher concentration of threatened and endemic species.
plant conservation, species diversity, Central Visayas, community structure, threatened species
The Philippines is one of the megabiodiverse countries and is considered a priority conservation area (
Species diversity can be influenced by many types of environmental changes (
Study site
The Balinsasayao Twin Lakes Natural Park (BTLNP) is located in Negros Oriental and is surrounded by Lakes Balinsasayao and Danao (Fig.
The study was conducted along the elevational gradient of Balinsasayao Twin Lakes Natural Park from April to May 2019. Three study forests were pre-determined, based on the general description by
The sampling plots in Balinsasayao Twin Lakes Natural Park (BTLNP) with the plot codes (Plots) used in the cluster analysis. Geographic coordinates (latitude, longitude) in each plot are in decimal format and the elevation is in metres above sea level (m a.s.l.).
Plot | Forest type | Latitude | Longitude | Elevation (m a.s.l.) |
LLQ1 | Lowland |
|
|
701 |
LLQ2 | Lowland |
|
|
680 |
LLQ3 | Lowland |
|
|
656 |
SMQ1 | Sub-montane upper part of the lake |
|
|
993 |
SMQ2 | Sub-montane upper part of the lake |
|
|
1034 |
SMQ3 | Sub-montane upper part of the lake |
|
|
1070 |
SMQ4 | Sub-montane along the lake |
|
|
886 |
SMQ5 | Sub-montane along the lake |
|
|
906 |
SMQ6 | Sub-montane along the lake |
|
|
923 |
SMQ7 | Sub-montane along the lake |
|
|
895 |
SMQ8 | Sub-montane along the lake |
|
|
877 |
SMQ9 | Sub-montane along the lake |
|
|
873 |
SMQ10 | Sub-montane along the lake |
|
|
867 |
SMQ11 | Sub-montane along the lake |
|
|
905 |
SMQ12 | Sub-montane upper part of the lake |
|
|
935 |
MOQ1 | Montane |
|
|
1670 |
MOQ2 | Montane |
|
|
1699 |
MOQ3 | Montane |
|
|
1710 |
Data recorded in the field were:
Relative locations of trees were sketched for tracing and monitoring purposes. For small sized plants (understorey and ground vegetation), these data were obtained:
Identification and nomenclature were aided using the following strategies:
Herbarium specimens were not collected during the first year of the study because of the limited time and collection restrictions.
The status of plant diversity was assessed using parameters of species abundance, such as: frequency, density and dominance and summed as the importance value. The Importance Value (IV) is the summation of the relative frequency, relative density and relative dominance combined as one.
The formulae to derive the IV are based on
Density = total number of individuals of a species / area sampled
Relative Density = density of species / total densities of all species x 100
Dominance = basal area (DBH area) of species / total area sampled
Relative Dominance = dominance of species / total dominance of all species x 100
Occurrence = number of plots a species is observed / total number of plots established
Frequency = number of occurrences / total number of occurrences
Relative Frequency = frequency of species / total of frequencies x 100
Importance value (IV) = Relative Density + Relative Dominance + Relative Frequency
On the other hand, the index of diversity was determined using the equations of
Shannon Diversity Index (H') = -∑ pi(LNpi)
where pi is the proportion (n/N) of individuals of one particular species found (n), divided by the total number of individuals found (N)
Pielou's Evenness Index (E') = H’/ln(s)
where s = number of species
Simpson Diversity Index (D) = 1 - (∑n(n-1)/N(N-1)
where n = the total number of organisms of a particular species; N = the total number of organisms of all species.
Values were interpreted, using the descriptions (Table
Relative Value Rating |
Species Diversity (H’) |
Evenness (E’) |
Very High |
3.50 – above |
0.75 – 1.00 |
High |
3.00 – 3.49 |
0.50 – 0.74 |
Moderate |
2.50 – 2.99 |
0.25 – 0.49 |
Low |
2.00 – 2.49 |
0.15 – 0.24 |
Very Low |
0.00 – 1.99 |
0.05 – 0.14 |
Hierarchical cluster analysis (HCA) of plots was done using Jaccard's similarity index from Paleontological Statistics (PAST version 2.17c) (
The conservation status of species was based on local, ‘‘The National List of Threatened Philippine Plants and their Categories” (
Floristic composition and community structure
A total of 351 plant species which represents 230 genera belonging to 103 families were recorded. These species were further classified as 183 species of trees, 54 herbs, 51 shrubs, 41 ferns and 22 vines. Angiosperms dominated the area, comprising 308 different plant species, while only two species were recorded for gymnosperms (Agathis philippinensis and Dacrydium beccarii). At the family level, the top 10 most abundant taxa included Rubiaceae (20 species; 6%), Lauraceae (19 species; 5%), Moraceae (18 species, 5%), Fabaceae (14; 4%), Myrtaceae (14 species; 4%), Urticaeae (11 species; 3%), Araceae (10 species; 3%), Euphorbiaceae (10 species; 3%), Melastomataceae (10 species; 3%) and Arecaceae (8 species; 2%). Moreover, amongst the elevational gradients surveyed, sub-montane ecosystems had the highest species richness (n = 287), followed by lowland (n = 96) and, lastly, montane (n = 53).
To understand the similarity of species in response to Lowland types along the elevation gradient, a cluster analysis was performed (Fig.
Dendogram of 18 sampling plots generated through UPGMA using Jaccard Similarity Index. Bootstraping was done at n = 1000; cophenetic correlation is 0.94. SMQ1-SMQ3 and SMQ12 are plots located in submontane (the upper part of the lake); SMQ4-SMQ11 - plots located in submontane (along the lake); LLQ1-LLQ3 are plots located in the lowland; MOQ1-MOQ3 – plots in the montane Lowland.
Species diversity and importance value
Lowland forest
The diversity estimates along the elevation gradient tend to follow this trend, high > very high > low (Fig.
Sub-montane Lowland
Sub-montane Lowland has an average diversity of H'= 4.20 which was categorised as very high, based on
Montane Lowland
This Lowland is characterised by the presence of numerous species of mosses, lichens and epiphytes. The dominant tree species was Dacrydium beccarii, while common shrubs dominating were Tasmannia piperita, Elaeocarpus argenteus, Melastoma crinitum, Vaccinium microphyllum, Freycinetia lagenicarpa, Astronia cumingiana, Dimorphanthera apoana and Eurya buxifolia and species of tree ferns were Cyathea negrosiana and Cyathea contaminans (Suppl. material
Conservation status and endemicity
Out of 351 species recorded at the study site, 30 (9%) were classified as threatened. Of the 30 threatened species, six (6) were categorised as critically endangered, four (4) endangered and 20 vulnerable. Remarkably, twenty five of these threatened species were found thriving in the Lowland of the sub-montane (Table
The findings in this study showed that the composition of flora of "BTLNP" in Negros Island is structured by elevational gradient. The distribution and composition of species within the three major Lowland types in the study site are defined and restricted by the differences in elevation which have a direct influence on the species ability to survive, compete and reproduce. Several tree species from the family Lauraceae, such as Litsea, Cinnamomum, Neolitsea, Actinodaphne, Phoebe and Machilus, were found thriving in the sub-montane forests. In addition, a number of dipterocarp species belonging to genus Shorea were also observed in this Lowland type. These observations conformed with the findings of
A sharp decline of species richness and diversity as altitude increased beyond 1,450 m a.s.l. was observed. This decline is due to the extreme environmental conditions, such as very low temperature, high relative humidity and low decomposition rate in the montane Lowland. Low temperature affects the respiratory process of trees through high retention of water vapour (
Thirty threatened and 88 Philippine endemic species were observed in BTLNP. Most of the threatened and endemic species were found in the sub-montane Lowland along the lake and are less in the montane and lowland. Elevation is known to be the greatest influencing factor in changes on floristic composition and diversity in tropical mountains (
The present study highlights the floristic composition, community structure and diversity of BTLNP along its elevational gradient, from lowland, sub-montane and montane. Our findings revealed that lowland and sub-montane forests had higher number of species compared to the montane Lowland. Furthermore, we observed a change in plant species composition along elevation gradients with a trend of decreasing species diversity with increasing elevation. It is worth noting that most of the threatened species observed in this study are mostly found in the sub-montane Lowland and are more easily accessible to humans compared to the montane.
These results indicate that BTLNP should follow appropriate conservation and management strategies to protect all elevations, but have a special focus on the sub-montane level, as well as threatened species. Additional care should be taken in the development of ecotourism, as human distrubance has been shown to alter the floral composition. Lastly, this study not only provides a floral baseline for BTLNP, but can also be used as a reference for additional species composition and elevational research within the Central Visayas.
The authors wish to thank Philippine Council for Agriculture, Aquatic and Natural Resources Research and Development (PCAARRD), for the financial assistance of the project. The Department of Environment and Natural Resources (DENR) and the Protected Area Management Board (PAMB) also with the BTLFAI Peoples Organization (PO) members of Balinsasayao Twin Lakes Natural Park is also acknowledged for allowing and supporting the researchers in the conduct of their study. We would like to also convey deep gratitude and heartfelt appreciation to the Central Visayas biodiversity survey team, Rochelyn Parba, Jessie Josol, Jessica Josol, Arianne Pacarat and Oscar Ido.
Biodiversity Assessment for Sustainable Management in Key Biodiversity Areas of Central Visayas
Bohol Island State University
Specimen Collection Gratuitous Permit (GP) Number VII-2019-04 (DENR Region VII)
The authors declare that there is no conflict of interest.
RDen = relative density; RFreq = relative frequency; RDom = relative dominance; IV = importance value.