Biodiversity Data Journal :
Research Article
|
Corresponding author: Wilbert A. Aureo (wilbert.aureo@bisu.edu.ph)
Academic editor: Anatoliy Khapugin
Received: 24 Jun 2020 | Accepted: 25 Sep 2020 | Published: 29 Dec 2020
© 2020 Wilbert Aureo, Tomas Reyes, Francis Carlo Mutia, Reizl Jose, Mary Beth Sarnowski
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 TD, Mutia FCU, Jose RP, Sarnowski MB (2020) Diversity and composition of plant species in the forest over limestone of Rajah Sikatuna Protected Landscape, Bohol, Philippines. Biodiversity Data Journal 8: e55790. https://doi.org/10.3897/BDJ.8.e55790
|
Rajah Sikatuna Protected Landscape (RSPL), considered the last frontier within the Central Visayas region, is an ideal location for flora and fauna research due to its rich biodiversity. This recent study was conducted to determine the plant species composition and diversity and to select priority areas for conservation to update management strategy. A field survey was carried out in fifteen (15) 20 m x 100 m nested plots established randomly in the forest over limestone of RSPL from July to October 2019. Three hundred and sixty eight (368) species of plants were identified up to species level. This represented 327 angiosperms, one gymnosperm and 40 pteridophytes. Common plant families with more than 10 representative species were Moraceae, Meliaceae, Lauraceae, Dipterocarpaceae, Rubiaceae, Myrtaceae, Phyllanthaceae, Annonaceae, Araceae and Lauraceae. There were 93 (28%) endemic and 46 (14%) threatened species (vulnerable to critically endangered) observed. The cluster analysis and species accumulation curve suggests that plant species are not homogeneously distributed which implies that different management and conservation strategies should be implemented across RSPL. These results not only indicate the importance of RSPL, but also highlights areas with higher diversity and concentration of threatened and endemic species as a special area of concern. Furthermore, areas with high biodiversity value were recommended for immediate protection, while areas with low biodiversity value were recommended for reforestation programmes using species with high importance value.
Central Visayas, endemic, plant habit, species diversity, threatened species
The Philippines is one of the mega diverse (
The RSPL is one of the last parts in Bohol with relatively-undisturbed forest with a high diversity of plant species (
RSPL was declared a national park in 1987 and a protected landscape in 2000 by the Department of Environment and Natural Resources (DENR) (
Meanwhile, there is growing emphasis on the importance of conservation planning which identifies priority areas and directs limited conservation resources in a strategic manner to help address challenges to biological diversity. Approaches in conservation planning have evolved over the past few decades from focusing mainly on species to encompassing broader aspects of biodiversity.
Rajah Sikatuna Protected Landscape is located in the six municipalities of Batuan, Bilar, Carmen, Dimiao, Garcia Hernandez, Sierra Bullones and Valencia (Fig.
The general and specific locality of sampling sites in Rajah Sikatuna Protected Landscape (RSPL) with the plot codes (Plot) used in the hierarchical cluster analysis. Geographic coordinates (Lat: latitude, Long: longitude) in each plot are in decimal format.
Plot | Municipality | Barangay | Total land area (ha) | Forest cover (%) | Topography | Lat / Long |
1 | Bilar | Zamora | 459.87 | 80 | Steep to rolling |
|
2 | Valencia | Omjon | 714.63 | 60 | Sloping |
|
3a | Carmen | Montehermoso | 350.69 | 50 | Sloping |
|
3b | Carmen | Montehermoso | Rolling |
|
||
3c | Carmen | Montehermoso | Sloping |
|
||
4a | Sierra Bullones | Nan-od | 812.91 | 85 | Relatively flat to rolling |
|
4b | Sierra Bullones | Nan-od | Relatively flat |
|
||
4c | Sierra Bullones | Nan-od | Relatively flat |
|
||
4d | Sierra Bullones | Nan-od | Flat to rolling |
|
||
5a | Garcia Hernandez | Cambuyo | 403.12 | 65 | Rolling |
|
5b | Garcia Hernandez | Cambuyo | Rolling |
|
||
6a | Garcia Hernandez | Datag | 252.55 | 75 | Flat to rolling |
|
6b | Garcia Hernandez | Datag | Flat to rolling |
|
||
7 | Sierra Bullones | Bugsoc | 190.56 | 50 | Steep to rolling |
|
8 | Batuan | Cabacnitan | 306.46 | 50 | Relatively flat |
|
This study was conducted from 22 July- 7 October 2019. We followed the random nested sampling plot design by
Data recorded in the field were: (i) plant names from family down to species level; (ii) diameter at breast height (cm) and total height (m) of species in the canopy layer; (iii) plant groups of observed plants; and (iv) GPS coordinates of all corners of each segment and nested plots. Relative locations of trees were sketched for tracing and monitoring purposes. For small-sized plants (understorey and ground vegetation), these data were obtained; (i) number of individuals and (ii) estimated percent crown cover.
Identification and nomenclature were aided using the following strategies: (i) expert determination; (ii) use of flora databases (
Importance value
The relative density, relative frequency and relative dominance for each tree species in all plots were determined to obtain their importance value (IV), a standard measure in ecology that determines the rank relationships of species. A high importance value of species indicates a composite score for high relative species dominance, density and frequency and provides a basis on what species can be used for restoration.
To compute for the relative density, relative dominance and relative frequency, the following formulae were used (
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 times a species is encountered / 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
Diversity index
Shannon-Weiner (H') and Simpson (D) indices and eveness were 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)
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
The description of the diversity index value is presented in Table
Classification of diversity indices interpreted using the descriptions proposed by
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 |
Cluster analaysis and species accumulation curve
The species accumulation curve (SAC) was done using specaccum (vegan) method = random and fitted values, residuals, non-linear model coefficients using the fitspecaccum and hclust() functions in R package version 3.6.3 (
The conservation status of plant species was based on local, ‘‘The National List of Threatened Philippine Plants and their Categories” (
We identified six criteria which relate to biodiversity value (e.g. computed diversity, richness, abundance, importance value, basal area and number of threatened and endemic species) to select priority areas for conservation. Some of these criteria (e.g. endemicity, conservation status, richness) were used by
Floristic composition
Rajah Sikatuna Protected Landscape was recorded to have 368 plant species belonging to 223 genera from 101 families. Angiosperms were the most diverse represented by 327 species, while there were 40 species of pteridophytes and a single gymnosperm species (Gnetom gnemon) (Table
Plant groups of the species in each family recorded in RSPL. The entire list of plant species is available in Suppl. material
Family | Angiosperm | Gymnosperm | Pteridophyte | Total species |
Acanthaceae | 2 | 2 | ||
Aceraceae | 1 | 1 | ||
Actinidiaceae | 1 | 1 | ||
Anacardiaceae | 6 | 6 | ||
Annonaceae | 12 | 12 | ||
Apocynaceae | 8 | 8 | ||
Araceae | 11 | 11 | ||
Araliaceae | 5 | 5 | ||
Arecaceae | 7 | 7 | ||
Asparagaceae | 1 | 1 | ||
Aspleniaceae | 7 | 7 | ||
Asteraceae | 1 | 1 | ||
Athyriaceae | 1 | 1 | ||
Begoniaceae | 1 | 1 | ||
Bignonaceae | 1 | 1 | ||
Blechnaceae | 1 | 1 | ||
Boraginaceae | 2 | 2 | ||
Brownlowiaceae | 1 | 1 | ||
Burseraceae | 7 | 7 | ||
Calophyllaceae | 2 | 2 | ||
Cannabaceae | 1 | 1 | ||
Casuarinaceae | 1 | 1 | ||
Celastraceae | 4 | 4 | ||
Clusiaceae | 2 | 2 | ||
Combretaceae | 2 | 2 | ||
Connaraceae | 1 | 1 | ||
Cornaceae | 3 | 3 | ||
Cunoniaceae | 1 | 1 | ||
Davalliaceae | 1 | 1 | ||
Dilleniaceae | 1 | 1 | ||
Dipterocarpaceae | 17 | 17 | ||
Dryopteridaceae | 2 | 2 | ||
Ebenaceae | 3 | 3 | ||
Elaeocarpaceae | 2 | 2 | ||
Euphorbiaceae | 10 | 10 | ||
Fabaceae | 15 | 15 | ||
Fagaceae | 1 | 1 | ||
Flagellariaceae | 1 | 1 | ||
Gesneriaceae | 2 | 2 | ||
Gnetaceae | 1 | 1 | ||
Hymenophyllaceae | 3 | 3 | ||
Hypericaceae | 2 | 2 | ||
Hypoxidaceae | 1 | 1 | ||
Juglandaceae | 1 | 1 | ||
Lamiaceae | 2 | 2 | ||
Lauraceae | 10 | 10 | ||
Lecythidaceae | 1 | 1 | ||
Leeaceae | 1 | 1 | ||
Lindsaeaceae | 1 | 1 | ||
Lomariopsidaceae | 2 | 2 | ||
Lygodiaceae | 1 | 1 | ||
Magnoliaceae | 3 | 3 | ||
Malvaceae | 1 | 1 | ||
Maranthaceae | 1 | 1 | ||
Marattiaceae | 1 | 1 | ||
Meliaceae | 17 | 17 | ||
Menispermaceae | 5 | 5 | ||
Monimiaceae | 1 | 1 | ||
Moraceae | 21 | 21 | ||
Myristicaceae | 3 | 3 | ||
Myrsinaceae | 1 | 1 | ||
Myrtaceae | 13 | 13 | ||
Nephrolepidaceae | 3 | 3 | ||
Nyctaginaceae | 1 | 1 | ||
Oleaceae | 1 | 1 | ||
Orchidaceae | 3 | 3 | ||
Pandanaceae | 2 | 2 | ||
Pentapetaceae | 2 | 2 | ||
Pentaphylaceae | 1 | 1 | ||
Phyllanthaceae | 15 | 15 | ||
Piperaceae | 2 | 2 | ||
Pittosporaceae | 1 | 1 | ||
Poaceae | 2 | 2 | ||
Polygalaceae | 1 | 1 | ||
Polypodiaceae | 3 | 3 | ||
Primulaceae | 2 | 2 | ||
Proteaceae | 1 | 1 | ||
Pteridaceae | 2 | 2 | ||
Rhamnaceae | 1 | 1 | ||
Rosaceae | 3 | 3 | ||
Rubiaceae | 18 | 18 | ||
Rutaceae | 7 | 7 | ||
Salicaceae | 2 | 2 | ||
Sapindaceae | 6 | 6 | ||
Sapotaceae | 8 | 8 | ||
Schizaeaceae | 1 | 1 | ||
Selaginellaceae | 4 | 4 | ||
Simaroubaceae | 1 | 1 | ||
Smilacaceae | 1 | 1 | ||
Sparmanniaceae | 1 | 1 | ||
Stemonuraceae | 1 | 1 | ||
Sterculiaceae | 5 | 5 | ||
Strombosiaceae | 1 | 1 | ||
Symplocaceae | 2 | 2 | ||
Tectariaceae | 4 | 4 | ||
Thelypteridaceae | 3 | 3 | ||
Thymelaeaceae | 4 | 4 | ||
Urticaceae | 7 | 7 | ||
Vitaceae | 3 | 3 | ||
Zingiberaceae | 2 | 2 | ||
Total | 327 | 1 | 40 | 368 |
In the analysis of similarity of vegetation structure and composition of each plot, three main floristic groups were identified (Fig.
Plant diversity
High diversity was observed in all established plots in RSPL (Fig.
Species importance value
Based on the computed IV (Table
Plot | Latin name | Family name | Common name | IV |
1 | Trigonostemon philippinensis Stapf in Elmer | Euphorbiaceae | Croton | 44.01 |
Mallotus cumingii Müll.Arg. | Euphorbiaceae | Apanang | 16.02 | |
Shorea contorta S.Vidal, | Dipterocarpaceae | White lauan | 15.84 | |
2 | Symplocos odoratissima (Blume) | Symplocaceae | Bayokbok | 33.02 |
Hopea acuminata Merr. | Dipterocarpaceae | Manggachapui | 31.41 | |
Tarennoidea wallichii (Hook.f.) | Rubiaceae | Pototan gubat | 21.83 | |
3a | Symplocos odoratissima (Blume) | Symplocaceae | Bayokbok | 93.16 |
Diplodiscus paniculatus Turcs. | Brownlowiaceae | Balobo | 49.36 | |
Aglaia exstipulata (Griff.) W.Theob. | Meliaceae | Aglia gagmay | 15.12 | |
3b | Chisocheton cumingianus (C.DC.) | Meliaceae | Balukanag | 38.70 |
Pterocymbium tinctorium (Blanco) Merr | Sterculiaceae | Taluto | 37.93 | |
Trigonostemon philippinensis Stapf in Elmer | Euphorbiaceae | Croton | 37.48 | |
3c | Ficus ampelas Burm.f. | Moraceae | Upling gubat | 31.37 |
Trigonostemon philippinensis Stapf in Elmer | Euphorbiaceae | Croton | 28.49 | |
Chisocheton cumingianus (C.DC.) | Meliaceae | Balukanag | 17.04 | |
4a | Shorea squamata (Turcz.) Benth. & Hook.f. ex DC | Dipterocarpaceae | Mayapis | 34.31 |
Symplocos odoratissima (Blume) | Symplocaceae | Bayokbok | 22.66 | |
Diplodiscus paniculatus Turcs. | Brownlowiaceae | Balobo | 20.18 | |
4b | Alangium meyeri Merr. | Cornaceae | Putian | 56.36 |
Pterospermum celebicum Miq. | Pentapetaceae | Bayokbayokan | 11.14 | |
Litsea cordata (Jack) Hook.f. | Lauraceae | Marang timber | 10.20 | |
4c | Shorea squamata (Turcz.) Benth. & Hook.f. ex DC | Dipterocarpaceae | Mayapis | 34.13 |
Pterocymbium tinctorium (Blanco) Merr | Sterculiaceae | Taluto | 27.27 | |
Chisocheton cumingianus (C.DC.) | Meliaceae | Balukanag | 21.25 | |
4d | Lithocarpus coopertus (Blanco) | Fagaceae | Ulaian | 26.56 |
Shorea squamata (Turcz.) Benth. & Hook.f. ex DC | Dipterocarpaceae | Mayapis | 24.37 | |
Hopea acuminata Merr. | Dipterocarpaceae | Manggachapui | 19.46 | |
5a | Trigonostemon philippinensis Stapf in Elmer | Euphorbiaceae | Croton | 32.70 |
Aglaia edulis (Roxb.) Wall. | Meliaceae | Malasaging | 29.19 | |
Chisocheton cumingianus (C.DC.) | Meliaceae | Balukanag | 22.88 | |
5b | Gomphandra mappioides Valeton, Crit. Overz | Stemonuraceae | Taguibokbok | 32.84 |
Chisocheton cumingianus (C.DC.) | Meliaceae | Balukanag | 24.57 | |
Hopea acuminata Merr. | Dipterocarpaceae | Manggachapui | 17.71 | |
6a | Pterocymbium tinctorium (Blanco) Merr | Sterculiaceae | Taluto | 43.46 |
Symplocos odoratissima (Blume) | Symplocaceae | Bayokbok | 23.02 | |
Lithocarpus coopertus (Blanco) | Fagaceae | Ulaian | 20.07 | |
6b | Shorea contorta S.Vidal, | Dipterocarpaceae | White lauan | 34.20 |
Trigonostemon philippinensis Stapf in Elmer | Euphorbiaceae | Croton | 18.09 | |
Baccaurea lanceolata (Miq.) | Phyllanthaceae | Mala Ulalian | 14.62 | |
7 | Symplocos odoratissima (Blume) | Symplocaceae | Bayokbok | 38.35 |
Lepiniopsis ternatensis Valeton | Apocynaceae | Paginga | 19.80 | |
8 | Lithocarpus coopertus (Blanco) | Fagaceae | Ulaian | 15.89 |
Trigonostemon philippinensis Stapf in Elmer | Euphorbiaceae | Croton | 15.59 | |
Shorea guiso (Blanco) Blume | Dipterocarpaceae | Guijo | 15.49 |
Tree stucture and density
Tree species with diameter class (DBH) of 10-20 cm had the highest proportion of individuals across all size classes (57%), followed by 21-30 cm (24%), 31-40 cm (11%) and 40+ cm (7%) (Table
Diameter class (cm) | Plot | ||||||||||||||||
1 | 2 | 3a | 3b | 3c | 4a | 4b | 4c | 4d | 5a | 5b | 6a | 6b | 7 | 8 | Total | % | |
10 to 20 | 69 | 35 | 23 | 31 | 76 | 102 | 83 | 64 | 61 | 66 | 60 | 86 | 93 | 95 | 69 | 1013 | 57 |
21 to 30 | 20 | 18 | 20 | 33 | 24 | 30 | 42 | 38 | 33 | 31 | 26 | 22 | 29 | 32 | 31 | 429 | 24 |
31 to 40 | 5 | 16 | 13 | 13 | 12 | 15 | 11 | 21 | 20 | 11 | 10 | 9 | 19 | 14 | 11 | 200 | 11 |
41 above | 11 | 13 | 8 | 14 | 10 | 7 | 6 | 5 | 11 | 5 | 2 | 8 | 4 | 10 | 12 | 126 | 7 |
Threatened and endemic species
Forty-six (46) species were recorded as threatened (critically endangered, endangered and vulnerable) under either the Philippine Red List (
Conservation status and endemicity of plant species recorded in each plot.
Conservation status | Plot | ||||||||||||||
1 | 2 | 3a | 3b | 3c | 4a | 4b | 4c | 4d | 5a | 5b | 6a | 6b | 7 | 8 | |
Threatened | 25 | 36 | 21 | 19 | 20 | 26 | 23 | 33 | 37 | 20 | 27 | 23 | 36 | 29 | 27 |
Endemic | 22 | 27 | 28 | 24 | 29 | 24 | 34 | 31 | 37 | 30 | 27 | 29 | 29 | 25 | 31 |
Some threatened tree species of Rajah Sikatuna Protected Landscape. Dipterocarpus grandiflorus (A), Hopea philippinensis (B), Shorea guiso (C), Hopea acuminata (D), Shorea squamata (E), Vatica mangachapoi (F), Pterocarpus indicus (G), Anisoptera thurifera (H), Hopea quisumbingiana (I), Shorea polysperma (J), Shorea contorta (K) and Shorea astylosa (L).
The biodiversity value of an area is always measured in terms of species richness and the number of endemic and threatened species present (
Plot | Endemic | Richness | Abundance | Threatened | Basal area | Diversity | Rank |
6b | 29 | 134 | 331 | 36 | 145 | 3.68 | 1 |
4d | 37 | 123 | 325 | 37 | 125 | 3.77 | 2 |
4b | 34 | 133 | 332 | 23 | 142 | 3.56 | 3 |
4c | 31 | 109 | 307 | 33 | 128 | 3.56 | 4 |
6a | 29 | 124 | 337 | 23 | 125 | 3.54 | 5 |
7 | 25 | 101 | 325 | 29 | 151 | 3.64 | 6 |
4a | 24 | 95 | 272 | 26 | 154 | 3.64 | 7 |
2 | 27 | 105 | 326 | 36 | 82 | 3.45 | 8 |
8 | 31 | 101 | 248 | 27 | 123 | 3.65 | 9 |
5b | 27 | 102 | 268 | 27 | 98 | 3.67 | 10 |
1 | 22 | 107 | 317 | 25 | 105 | 3.43 | 11 |
3a | 28 | 114 | 307 | 21 | 64 | 3.42 | 12 |
3c | 29 | 100 | 285 | 20 | 122 | 3.26 | 13 |
5a | 30 | 100 | 302 | 20 | 113 | 3.13 | 14 |
3b | 24 | 93 | 255 | 19 | 91 | 3.33 | 15 |
The study provides baseline data on plant composition in the forest over limestone of the Rajah Sikatuna Protected Landscape. A total of 368 plant species, belonging to 223 genera from 101 families, were recorded. This result is higher compared to similar studies conducted in forest over limestone of Dinagat Island which accounted for 144 plant species, belonging to 50 families and 88 genera (
Species were also observed to form their own phytosociological group, distinct from each other. Symplocos odoratissima, in particular, was found common on top of the mountain hills though was also present in relatively flat terrain. Trigonostemon philippinensis and Pterocymbium tinctorium, on the other hand, formed thickets along valleys and flat rocky substrates in between mountain hills which may occasionally be submerged in flowing water during rainy season, while Shorea contorta and Shorea squamata extended their dominance from hilly areas to relatively rolling undisturbed closed canopy forests. These observations conformed to the findings of
All plots showed high diversity in both Shannon (H’) and Simpson’s (D) Diversity Indices, with the plot in Brgy. Nan-od exhibited the highest value of 3.77, followed by Plot 6b in Datag, while the lowest diversity was recorded in Plot 5a in Cambuyo. It is important to note that, in terms of the Shannon Diversity Index, the ordering of the plots was likely affected by the topography, dense forest cover and maturity of the forest (
Of the 368 plant species encountered, 93 species were found endemic to the Philippines and 46 were listed as threatened, based on
Areas with the highest totals, Datag, Nan-od, Bugsoc and Omjon, summed from the identified criteria, need the most immediate conservation actions. The species in these areas are highly threatened due to their economic value, being locally accessible and often used for construction and as firewood. Barangays Bugsoc, Datag and Nan-od exhibited very high species diversity, especially in the canopy layer (tree species). These Barangays are similarly the areas displaying high basal area of tree species mostly belonging to the family Dipterocarpaceae - a very good source of timber - highlighting the risk of future exploitation. In contrast, Barangays Montehermoso, Bilar and Cambuyo exhibit relatively-lower species diversity. These areas are located near communities where we can expect high human encroachment in the forest. In addition, it is worth mentioning that despite low diversity in Barangay Montehermoso, Shorea negrosensis, a critically-endangered species was only observed in this area, thereby presenting a unique conservation priority. While this present study does not account for all species within RSPL, these results provide a data-driven basis for prioritisation of areas for conservation.
RSPL has a very high diversity value and is a habitat to at least 368 flora species. A total of 93 plants are exclusively found in the Philippines and a significant number of forty six (46) threatened species were observed in the area. Many of these threatened species like White lauan (Shorea contorta), Mayapis (Shorea squamata) and Mangachapui (Hopea acuminata), although listed as critically endangered, were commonly found within its natural forest. Therefore, it is recommended that immediate conservation and management activities should be conducted to parts of RSPL with high diversity and number of threatened and endemic species like Barangays Bugsoc, Datag, Nan-od and Omjon to save these plants from extirpation. However, for areas with low diversity and number of threatened and endemic plants, we suggest that reforestation initiatives, like assisted natural regeneration (ANR) and rain-forestation, be implemented. Additionally, species with high importance value in each area are also recommended to be used in the said reforestation programme. Lastly, future studies should also consider plant associations and environment interactions.
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) members of Rajah Sikaituna Protected Landscape are also acknowledged for allowing and supporting the researchers in the conduct of their study. We would like to also convey deep gratefulness and heartfelt appreciation to the Central Visayas biodiversity survey team, Rochelyn Parba, Jessie Josol, Jessica Josol, Arianne Pacarat and Oscar Ido and the National Museum of the Philippines - Herbarium and Botany Section consultants: Danilo Tandang, John Rey Callado, and Edwin Tadiosa and Fulgent P. Coritico of CMU-CEBREM for their help identifying the seed plants and fern species. Special thanks to BS Forestry students who joined the survey and utilised study data for their coursework.
Biodiversity Assessment for Sustainable Management in Selected Key Biodiversity Areas of Central Visayas
Bohol Island State University (BISU)
Specimen Collection Gratuitous Permit (GP) Number VII-2019-04 (DENR Region VII) and Prior Informed Consent (PIC) were secured before conducting the research. The research proposal was also endorsed by the Central Visayas Agriculture, Aquatic and Natural Resources Research and Development Consortium (CVAARRDEC) to PCAARRD for funding after passing ethics and security standards of the university.
The authors declare that there is no conflict of interest.
This contains all the raw data and analyses performed in this article