Research Article |
Corresponding author: Dung Van Tran ( trandungfuv@gmail.com ) Academic editor: Christoph Knogge
© 2023 Khoa Van Phung, Dung Van Tran, Hai Thanh Dong, Vinh Quang Luu, Van Bac Bui, Thinh Tien Vu.
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:
Van Phung K, Van Tran D, Thanh Dong H, Quang Luu V, Bac Bui V, Tien Vu T (2023) Changes in suitable habitat for the critically endangered Northern white-cheeked gibbon (Nomascus leucogenys) in the Western Nghe An Biosphere Reserve, Vietnam: Implication for conservation. Nature Conservation 51: 167-188. https://doi.org/10.3897/natureconservation.51.90373
|
Several recent studies have highlighted that change in land use and land cover (LULC) is the main threat causing the decline and extinction of certain species. Gibbons (Hylobatidae) could be excellent examples, on account of their potential for extinction in the near future under the effects of LULC changes due to their particular ecological traits. This study aims to model the current suitable habitat of the Northern white-cheeked gibbon (Nomascus leucogenys (Ogilby, 1840)) in the Western Nghe An Biosphere Reserve (BR), Vietnam and assess the changes in its suitable habitat following the changes in LULC from 1990 to 2020. The maximum entropy approach (MaxEnt) was used to predict the suitable habitat of the gibbon using its occurrence localities and environmental predictors. The model analysis showed that the “Distance to Agriculture” variable had the strongest impact on the gibbons’ suitable habitat. Our results predicted the present suitable habitat of the gibbon species at approximately 4,022.42 km2 (30.95% of the overall BR area) in three spatially separated areas inside the Western Nghe An BR. Furthermore, the suitable habitat areas of the gibbon in 1990, 2000, and 2010 were projected at roughly 4,347.68 km2, 4,324.97 km2, and 2,750.21 km2, respectively, following a decreasing trend from 1990 to 2010, but a gradual increase between 2010 and 2020. The suitable habitat of the gibbon inside three core protected areas (Pu Mat National Park, Pu Huong, and Pu Hoat Nature Reserves) showed a continually increasing trend from 1990 to 2020. Our results highlighted the influence of LULC changes and the role of the protected area network in gibbon conservation. The information from the study provides a quantitative baseline for the future conservation of the critically endangered gibbon in the Western Nghe An BR.
Gibbon, land use & land cover, MaxEnt, Nomascus, species distribution modeling, suitable distribution
Southeast Asia contains approximately 15% of the world’s tropical forests and is home to at least four globally important biodiversity hotspots (
The northern white-cheeked gibbon (NWCG) (Nomascus leucogenys) is classified as a Critically Endangered species on the International Union for Conservation of Nature (IUCN) Red List (
The largest population of NWCG was found in Pu Mat National Park (NP), with 22 confirmed groups; the site was also considered as the highest priority area for conserving the species (
In this study, we aimed to predict the suitable habitat of NWCG in the Western Nghe An BR by using species distribution modeling (MaxEnt) based on identified presence localities of the species and environmental predictors. We also attempted to project suitable habitat fluctuation in relation to the changes in LULC from 1990 to 2020. Based on the results, we discussed the impact of LULC changes and the role of the protected area network in gibbon conservation in Vietnam. Our results are a baseline for researchers, conservationists, and wildlife and habitat managers to aid decision-making and plan future conservation strategies for NWCG in the Western Nghe An BR, Vietnam.
The Western Nghe An Biosphere Reserve (BR) was recognized by UNESCO as the 6th Biosphere Reserve of Vietnam in 2007. The Western Nghe An BR (18°34'43"–19°59'44"N, 103°52'28"–105°30'07"E) is located in the western part of Nghe An province, central Vietnam, covering an area of 12,997.95 km2 with three functional zones: the core zone (1,683.01 km2), buffer zone (6,085.47 km2), and transition zone (5,229.47 km2). The Biosphere Reserve is expected to create a green corridor between the three protected areas: Pu Mat NP, Pu Huong NR, and Pu Hoat NR. The forest cover of the Western Nghe An BR is approximately 66.4% of the total area. The area is home to more than 3,000 vascular plants and more than 940 vertebrate animals, with several species listed as threatened, rare, and endemic. More than 900,000 people in five indigenous minority groups currently reside within Western Nghe An BR (
To predict the suitable habitat of NWCG in the Western Nghe An BR, we collected the occurrence of the endangered species through our field surveys in Pu Mat NP, Pu Huong NR, and Pu Hoat NR from March–May 2021 (
To predict the suitable habitat of NWCG in the Western Nghe An BR, we obtained variables on climate, topography, and LULC. The variables were selected based on our knowledge of the habitat of the gibbon and consulting the available research sources. For the climatic variables, we downloaded 19 predictors with the highest available resolution (30 arcseconds) from the World Clim database (http://www.worldclim.com/;
The LULC in the Western Nghe An BR in A 1990 B 2000 C 2010, and D 2020 extracted from
The predictor variables used for generating the habitat suitability of NWCG in the Western Nghe An BR.
No. | Name | Sources | Description |
---|---|---|---|
1 | Bio01 | WorldClim | Annual Mean Temperature |
2 | Bio03 | Isothermality (BIO2/BIO7) (×100) | |
3 | Bio07 | Temperature Annual Range (BIO5-BIO6) | |
4 | Bio12 | Annual Precipitation | |
5 | Bio13 | Precipitation of Wettest Month | |
6 | Bio15 | Precipitation Seasonality (Coefficient of Variation) | |
7 | Bio16 | Precipitation of Wettest Quarter | |
8 | Elevation | Earthexplorer.usgs.gov | Height above sea level |
9 | Slope | Degree of rise/run | |
10 | Aspect | Direction a slope face | |
11 | LULC | ALOS Research and Application Project | Land Use and Land Cover |
12 | Distance to settlement | Distance to residential land | |
13 | Distance to Agriculture | Distance to agricultural land |
To generate the suitable habitat model for NWCG in the Western Nghe An BR, we applied the Maximum Entropy approach (MaxEnt) version 3.4.4 available from https://biodiversityinformatics.amnh.org/open_source/maxent/. The MaxEnt model predicts the habitat suitability of species based on the probable occurrence of species in distinct localities finding the maximum entropy distribution of environmental predictors (
To assess the changes in the suitable habitat of NWCG from 1990 to 2020, we ran separately two models for 2010 and 2020 based on the presence localities and environmental data set of each year. For the model of 2020, we ran the model using 10-folds cross-validation to evaluate the model, while we applied the jackknife method due to the small sample size for the model 2010 (
The relationship between the predicted suitable habitat of species and environmental variables in the MaxEnt model was shown by response curves. We presented the response curves of the most contributing variables to identify the main impact of the variables on the change of suitable habitat of NWCG in the Western Nghe An BR.
The result of the MaxEnt model is presented in a logistic format ranging from 0 to 1, in which the higher values mean higher suitability. The threshold selection should be determined according to the objectives of the study (
The MaxEnt model predicted the habitat suitability of NWCG based on the available presence localities and predictor data set of LULC, climate, and topography with mean test AUC at 0.896 ± 0.043 and 0.936 ± 0.084 for the model of 2020 and 2010, respectively (Fig.
The results of AUC curves in developing the suitable habitat model of NWCG in A 2020, and B 2010. The red line showed the mean of AUC, and the blue area showed the standard deviation.
For the contribution of each environmental variable in the model of suitable habitat in 2020, the top three variables were “Distance to Agriculture” (39.2%), “Temperature Annual Range – Bio7” (25.3%), and “Precipitation of Wettest Quarter – Bio16” (23.3%). For the model of 2010, the highest contribution variables were “Distance to Agriculture” with 55.4%, followed by “Distance to Settlement” and “Elevation” with 29.3% and 6.6%, respectively (Table SI-II). The response curves in the MaxEnt model show the changes in habitat suitability in response to the changes in the predictors used in the model. The response curve of the highest contribution variables of both models (2020 and 2010) was mainly identical when the habitat suitability of NWCG increased significantly for the distance to agricultural land below around one km, and rose gradually for the distance above about one km (Fig.
The total suitable habitat of NWCG in Nghe An BR under the current land cover condition was estimated at around 4,022.42 km2 (30.95% of the overall BR area; Table
The suitable habitat of NWCG under present landcover condition (2020) predicted by the MaxEnt model.
Area of suitable habitat of NWCG projected by the MaxEnt model in 1990, 2000, 2010, 2020 (unit: km2).
Habitat suitability | Year | |||
---|---|---|---|---|
1990 | 2000 | 2010 | 2020 | |
Low | 3,091.51 | 3,221.81 | 1,719.21 | 2,488.18 |
Moderate | 871.28 | 793.30 | 708.17 | 1,173.40 |
High | 384.89 | 309.86 | 322.83 | 360.84 |
Total | 4,347.68 | 4,324.97 | 2,750.21 | 4,022.42 |
Using the MaxEnt model, we projected the suitable habitat of NWCG in 1990 and 2010 (Fig.
The suitable habitat of NWCG in the Western Nghe An BR in A 1990 B 2000 C 2010, and D 2020 (Pu Mat NR was established in 1995, and upgraded to be national park in 2001; Pu Huong, and Pu Hoat NRs were established 2003, and 2013, respectively).
The extracted LULC in the Western Nghe An BR from 1990 to 2020 constituted seven main types, including Settlement, Agriculture, Barren land, Broadleaf forest, Plantation, Water, and Bamboo forest (Table
Area and change of the LULC categories in the Western Nghe An BR from 1990 to 2020 extracted from
No | LULC categories | Year | Change | |||||
---|---|---|---|---|---|---|---|---|
1990 | 2000 | 2010 | 2020 | 1990–2000 | 2000–2010 | 2010–2020 | ||
1 | Settlement | 0.56 | 1.86 | 13.79 | 34.75 | 1.30 | 11.93 | 20.97 |
2 | Agriculture | 1,789.05 | 2,008.77 | 1,835.63 | 1,707.34 | 219.72 | -173.14 | -128.29 |
3 | Barren land | 791.70 | 869.97 | 899.11 | 549.95 | 78.27 | 29.13 | -349.16 |
4 | Broadleaf Forest | 6,366.36 | 7,145.05 | 6,902.46 | 7,695.72 | 778.69 | -242.59 | 793.26 |
5 | Plantation | 3,930.07 | 2,681.34 | 2,733.05 | 2,275.61 | -1,248.73 | 51.71 | -457.43 |
6 | Water | 68.27 | 83.86 | 69.91 | 125.89 | 15.59 | -13.94 | 55.96 |
7 | Bamboo forest | 889.56 | 1,044.72 | 1,381.62 | 1,446.31 | 155.16 | 336.90 | 64.68 |
The overlay between the existing protected areas within the Western Nghe An BR, including Pu Mat NP, Pu Huong, and Pu Hoat NRs, and the suitable habitat of NWCG predicted by the MaxEnt model was only approximately 868.05 km2 (25.11%) of suitable habitat lies within protected areas (PAs) in 2000. In 1990, three protected areas in the core zone of the BR were not yet established. In 2010, the areas of suitable habitat inside PAs increased to 1,110.13 km2. Then, the predicted suitable habitat rose to approximately 1,573.71 km2 in 2020 (Table
Area of suitable habitat inside and outside borders of protected areas (Pu Mat NP, Pu Huong NR, and Pu Hoat NR*) of NWCG projected by the MaxEnt model in 1990, 2000, 2010, 2020.
Suitability | 1990 | 2000 | 2010 | 2020 | ||||
---|---|---|---|---|---|---|---|---|
Inside | Outside | Inside | Outside | Inside | Outside | Inside | Outside | |
Low | 0 | 3,091.51 | 382.32 | 2,839.49 | 437.10 | 1,282.11 | 723.43 | 1,764.75 |
Moderate | 0 | 871.28 | 348.16 | 445.15 | 444.19 | 263.99 | 613.29 | 560.11 |
High | 0 | 384.89 | 137.57 | 172.28 | 228.84 | 93.99 | 236.99 | 123.85 |
Total | 0 | 4,347.68 | 868.05 | 3456.91 | 1,110.13 | 1,640.08 | 1,573.71 | 2,448.70 |
The AUC value of both models in 2010 and 2020 shows a high capacity to predict the suitable habitat of NWCG with a good performance. The MaxEnt model has been broadly employed to predict the potential distribution of primates (
In the present study, we tried to predict the suitable habitat changes of NWCG in the Western Nghe An BR. Modeling a partial distribution could not provide a complete species range because a part of the environmental conditions and presence records of species outside of the target area of the model are not included (
Our models indicated that the “Distance to Agriculture” was the most significant variable for both predicted models 2020, and 2010 of NWCG, with the highest contributions to the models. The habitat suitability was likely to increase when the distance to agricultural land value was high (Fig.
The suitable habitat generated by species distribution modeling could be very useful in suggesting areas to find new populations, especially for rare and threatened species (
Our current study indicates that LULC changes may significantly impact the distribution of gibbons in the Western Nghe An BR, central Vietnam, as projected from 1990 to 2020. Knowledge of habitat transformation plays a crucial role in making decisions regarding biodiversity conservation (
The Distance to Settlement variable was the second most important variable for the habitat suitability model of NWCG, although the Settlement was the smallest area among seven main LULC types in the Western Nghe An BR, probably indicating the high impact of human residence settlement areas to habitat suitability of NWCG. The settlement area in the transition zone of the Western Nghe An BR continually increased from 1990 to 2020, probably due to urbanization and population expansion. Between 1990 and 2000, the agricultural land area increased, leading to a slight decrease in the suitable habitat of NWCG. The broadleaf forest areas significantly decreased between 2000 and 2010, reflecting the ineffective activities for habitat conservation in the area. In contrast, there was an increase in the broadleaf forest from 2010 to 2020 that was probably facilitated by government policies, notably the closing of natural forest policy in the late 2010s, and the National Action Programme on REDD+ in the period 2011–2020 (
Although the Western Nghe An BR has been established since 2007 with only three core protected areas, including Pu Mat NP, Pu Huong, and Pu Hoat NR, most areas of the Biosphere Reserve have remained currently intact even with the increasing settlement, agricultural land, and the fluctuation of broadleaf forest areas. As a positive human policy, establishing protected areas for protecting all environmental components, including populations, habitats of wild species and natural ecosystems, is crucial because of its mitigations of negative impacts on biodiversity at a local scale (
In this study, land-use change has a strong relationship with the suitable habitat of NWCG, which resulted in the contraction of the suitable habitat of the endangered gibbon from 2000 to 2010. Interestingly, an expansion of suitable habitat of NWCG was revealed in the period 2010 – 2020, which provides a positive sign for habitat conservation in the study areas. Additionally, a limitation of the study is the lack of information on specific habitat parameters that might be critical for defining the selection habitat of gibbons, such as the height of tree canopy, food preference, and availability. We recommended that more intensive surveys should be focused on the suitable habitat areas suggested by our model to reveal the forest quality of these areas.
Vietnam is located in the Mekong region, which was assessed as one of the 11 areas experiencing the largest forest loss. The average rate of deforestation in the Mekong region was approximately 0.4% per year (
Finally, we highly recommend the implementation of the following actions for protecting existing populations of NWCG and its habitat in the Western Nghe An BR: 1) implement population surveys, monitoring and conservation awareness programmes not only in three core protected areas (Pu Mat NP, Pu Hoat, and Pu Huong NR) but also in other highly suitable habitats that are predicted by our model, especially for the northern Pu Hoat NR, northern and southern Pu Mat NP; 2) upgrade the Pu Hoat watershed protection forest, the highly suitable habitat in the north of Pu Hoat NR, to be a Nature Reserve (Fig.
There is currently a lack of research on the changes in suitable habitats for gibbons, leading to limited efforts to protect and conserve their habitats. Here, we applied the species distribution modeling to predict the suitable habitat of NWCG in the Western Nghe An BR, and assessed the impact of changes in land use and land cover from 1990 to 2020 on the suitable habitat of this species. The current suitable habitat for the gibbon in the Western Nghe An BR was estimated at approximately 4,022.42 km2, mainly concentrated in three distinctive areas (Pu Mat NP, Pu Huong, and Pu Hoat NR), while the predicted suitable habitat area lies outside the existing protected areas at around 2,448.70 km2. Our result also indicated that due to the changes in LULC, the predicted suitable habitat of the gibbon decreased from 1990 to 2010 but gradually increased from 2010 to 2020. Based on our findings, we suggested that more survey efforts should be focused on areas of predicted suitable habitat regions outside existing protected areas, in addition to continuing efforts for suitable habitat areas within protected areas. We also highly recommend enriching young forests and protecting the preferred habitat for the gibbon, specifically broadleaf forests. Finally, we proposed several actions to safeguard NWCG and its habitat in the Western Nghe An BR.
We would like to thank the “Mainstreaming Natural resource management and Biodiversity conservation objectives into socio-economic development planning and management of biosphere reserve in Viet Nam” project (BR project) for their supports. We are grateful to Dr. Greg Nagle for proofreading this manuscript. We wish to thank Management Broad of Western Nghe An BR, Pu Mat NP, Pu Hoat, and Pu Huong NR, and our field assistants for their support during our field surveys. Many thanks to anonymous reviewers for improving a previous version of the manuscript.
Supporting information
Data type: figures and tables
Explanation note: SI–I. Map of each variable that was used for the final model. SI–II. The percent contribution of environmental variables. SI–III. The suitable habitat of N. leucogenys were predicted from occurrence localities in 2020. SI–IV. Area of suitable habitat of N. leucogenys projected by MaxEnt model in 1990, 2000, 2010, 2020 from the model of 2020 (unit: km2).