Research Article
Print
Research Article
Projected contraction of climatically suitable habitat for the forest-specialist anuran Quasipaa acanthophora under climate change in northern Vietnam
expand article infoLinh Thuy Thi Ha‡, Trang Hoang Nguyen‡, Thinh Van Nguyen§, Cuong The Pham|, Tien Quang Phan|, Truong Quang Nguyen|, Dung Van Tran#, Thomas Ziegler¤«, Dennis Rödder», Hai Ngoc Ngo˄˅
‡ Department of Environmental Science and Technology, School of Chemistry and Life Sciences, Hanoi University of Science and Technology, Hanoi, Vietnam
§ Central Institute for Natural Resources and Environmental Studies, Vietnam National University, Hanoi, Vietnam
| Institute of Biology, Vietnam Academy of Science and Technology, Hanoi, Vietnam
¶ Graduate University of Science and Technology, Vietnam Academy of Science and Technology, Hanoi, Vietnam
# Wildlife Department, Vietnam National University of Forestry, Hanoi, Vietnam
¤ Cologne Zoo, Cologne, Germany
« University of Cologne, Cologne, Germany
» LIB, Museum Koenig Bonn, Leibniz Institut zur Analyse des Biodiversitätswandels, Bonn, Germany
˄ Institute of Theoretical and Applied Research, Duy Tan University, Hanoi, Vietnam
˅ School of Medicine and Pharmacy, Duy Tan University, Da Nang, Vietnam

Corresponding author: Hai Ngoc Ngo ( ngongochai1@duytan.edu.vn )

Academic editor: Peter Mikulíček
© 2026 Linh Thuy Thi Ha, Trang Hoang Nguyen, Thinh Van Nguyen, Cuong The Pham, Tien Quang Phan, Truong Quang Nguyen, Dung Van Tran, Thomas Ziegler, Dennis Rödder, Hai Ngoc Ngo.
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: Ha LTT, Nguyen TH, Nguyen TV, Pham CT, Phan TQ, Nguyen TQ, Tran DV, Ziegler T, Rödder D, Ngo HN (2026) Projected contraction of climatically suitable habitat for the forest-specialist anuran Quasipaa acanthophora under climate change in northern Vietnam. Herpetozoa 39: 313-324. https://doi.org/10.3897/herpetozoa.39.e201361
ZooBank: urn:lsid:zoobank.org:pub:B8A8CF20-D798-4FC0-958F-A2AE7A8FD330
Open Access

Abstract

Anthropogenic impacts, including habitat loss and overexploitation, in combination with climate change, have resulted in rapid population declines and even local extirpations in many amphibians. As an endemic amphibian species in northern Vietnam, Quasipaa acanthophora is endangered by combined human impacts and potential effects of climate change, similar to other amphibian and reptile species occurring in the same habitats. In this study, we applied the Maximum Entropy algorithm in MaxEnt software to predict the current suitable habitats and their shifts in the future under scenarios of climate change. As the current prediction, the potential distribution of Q. acanthophora covers the known occurrence sites, including Lang Son, Bac Ninh, Quang Ninh provinces, in northeastern Vietnam, and extends to the southern border of Guangxi Province, China, with a total area of approximately 1,779 km2 within the study area, 1,567 km2 in the tropical evergreen rain forest, and only 231 km2 in the protected areas. Across all scenarios of climate change, future projections consistently indicate a substantial contraction in the potential distribution of Q. acanthophora. The study identifies priority areas, particularly forests in protected areas, for conservation measures targeting Q. acanthophora populations and their natural habitats.

Key Words

Mau Son Spiny Frog, MaxEnt, potential distribution, priority conservation areas, projection, tropical evergreen rain forest

Introduction

Under anthropogenic pressures, the global effects of climate change, including elevated CO2 levels, increasing temperatures, altered precipitation patterns, rising sea levels, and more frequent natural disasters (e.g., storms and floods), are exacerbating the biodiversity crisis (Bellard et al. 2012; IPBES 2019; IPCC 2021). These, in turn, threaten many species worldwide, driving their populations toward progressive decline and extinction risk (Schär et al. 2004; Diffenbaugh and Ashfaq 2010; Pimm et al. 2014; Ceballos et al. 2015). As ectotherms, amphibians are among the organism groups most strongly affected by climate change, due to the dependence of their physiology and ecology on environmental conditions (Carey and Alexander 2003), and these changes particularly affect species’ reproduction, metabolism and development (Deutsch et al. 2008). Luedtke et al. (2023) reported that 39% of amphibian species worldwide faced higher extinction risk due to climate change between 2004 and 2022.

The Mau Son Spiny Frog – Quasipaa acanthophora Dubois & Ohler, 2009, has been assessed as Vulnerable (VU) in the IUCN Red List (IUCN 2025) and listed as Endangered (EN) in the Vietnam Red Data Book (Nguyen et al. 2024), due to overexploitation for food consumption and habitat loss (Rowley et al. 2010; IUCN SSC Amphibian Specialist Group 2017; Nguyen et al. 2024; Pham et al. in press). The endemic species has only been recorded in the tropical evergreen rain forests, at the type locality of Mau Son Mountain, Lang Son Province, and two expanded provinces in northeastern Vietnam (i.e., Bac Ninh and Quang Ninh) (Dubois and Ohler 2009; Hecht et al. 2013; Nguyen et al. 2025). In these provinces, three sympatric tiger geckos of the genus Goniurosaurus (namely, G. huuliensis, G. lichtenfelderi and G. luii) were projected to be negatively affected by climate change under future scenarios (Ngo et al. 2021a, 2022, 2025). Pham et al. (in press) reported that forest degradation in coastal areas, including the natural habitats of Q. acanthophora in Quang Ninh Province, has been driven by extreme typhoons, which are becoming more intense and frequent as a consequence of climate change. The impacts of climate change on the target species, Q. acanthophora, may be further exacerbated under ongoing human impacts (Pham et al. in press).

In this study, we employed the Maximum Entropy algorithm in MaxEnt software to predict the potential distribution of Q. acanthophora, enveloping suitable climates in the tropical evergreen rain forests, and to further assess its potential shifts under future scenarios of climate change. We hypothesize that suitable habitats of Q. acanthophora will contract across all scenarios, and the negative impacts could even be far more severe under human impacts. Based on the projections, we identify potential priority forests for conservation safeguarding wild populations of Q. acanthophora.

Methods

The study site

The study site for all climatic predictions was restricted to latitudes of 20°30’ – 22°30'N and longitudes of 106° – 108°30'E (Fig. 1), with a background area of 43,400 km2, covering all known distributions of Q. acanthophora in northeast Vietnam and its potential occurrence in southern China (Dubois and Ohler 2009; Hecht et al. 2013; Nguyen et al. 2025).

Figure 1. 

Map of study site and occurrences of Quasipaa acanthophora in Vietnam.

Data collection

Two coordinates of Q. acanthophora from the Global Biodiversity Information Facility (GBIF) were downloaded (GBIF.org 2026), using the “rgbif” package (Chamberlain et al. 2024) in R v4.6.0 (R Core Team 2026). Additionally, interviews with local people and forest rangers identified 32 streams likely occupied by Q. acanthophora, and 42 coordinates from 16 streams were subsequently verified by our direct observations during field surveys (Pham et al. in press). Only one occurrence locality of Q. acanthophora was randomly selected within each square kilometer, using the “dismo” (v1.3-8) package in R (Aiello-Lammens et al. 2015; Hijmans et al. 2022). Such spatial filtering can improve the quality of prediction models by decreasing geographical bias and autocorrelation effects (Radosavljevic and Anderson 2014; Kiedrzynski et al. 2017). Finally, 20 representative occurrence records of Q. acanthophora were employed for further analyses (Fig. 1).

To predict the current potential distribution of Q. acanthophora, we obtained 19 bioclimatic variables from the Worldclim 2.1 at 30 arc seconds resolution (Fick and Hijmans 2017). For future projections, we used three global circulation models (GCMs) from the sixth Coupled Model Intercomparison Project (CMIP6): Meteorological Research Institute Earth System Model (MRI-ESM2-0) developed by the Meteorological Research Institute in Japan (Eyring et al. 2016; Yukimoto et al. 2019), Interdisciplinary Research on Climate version 6 (MIROC6) developed by the Atmosphere and Ocean Research Institute, the National Institute for Environmental Studies, and the Japan Agency for Marine-Earth Science and Technology, Japan (Tatebe et al. 2019), and MR) developed by the Beijing Climate Center, China Meteorological Administration (Wu et al. 2019). These climate models are able to capture regional monsoon dynamics, to simulate ocean-atmosphere interactions, and to predict tropical cyclone behavior. Three climate change scenarios of the Shared Socioeconomic Pathways (SSPs) were downloaded to assess the different level of climate change impacts over 20-year periods (i.e. 2050s: 2041–2060; 2070s: 2061–2080; 2090s: 2081–2100): SSP2-4.5 (“middle of the road”: temperatures are projected to rise by 2.7 °C by the end of the century), SSP3-7.0 (“regional rivalry”: to rise by around 3.2 °C to 4.1 °C) and SSP5-8.5 (“fossil-fuelled development”: to rise by around 4.7 °C to 5.1 °C) (O’Neill et al. 2016; Riahi et al. 2017). In this study, we only presented one projected map for each future climate scenario, by averaging the raster cell values of suitability across the three global circulation models (GCMs).

To select climatic conditions for predictions, we first employed 19 variables to import into the MaxEnt software v. 3.4.4 (Phillips et al. 2006), and applied the default settings to identify the variables that contribute most to the model (Suppl. material 1: table S2). Next, a pairwise Pearson correlation analysis of selected variables was conducted to mitigate the influence of potential multicollinearity, using the “usdm” (v2.1-7) package in R to identify and eliminate highly correlated environmental variables with coefficient values larger than |0.75| (Dormann et al. 2013; Naimi 2015). Our final dataset comprised six variables, including Mean Diurnal Range (Bio2), Isothermality (Bio3), Temperature Seasonality (Bio4), Max Temperature of Warmest Month (Bio5), Mean Temperature of Coldest Quarter (Bio11), Precipitation of Driest Quarter (Bio17), which were finally selected for subsequent analyses (Suppl. material 1: table S3).

Species distribution modeling

Occurrence data of Q. acanthophora and six selected climate variables were imported into the MAXENT software to predict the potential distribution at present and to assess its shifts in all future scenarios.

For the species within a restricted distribution, defining the study background to include environmentally dissimilar areas beyond the accessible region of Q. acanthophora may artificially inflate the AUC (Area Under the Curve) and specificity values by increasing the contrast between presence and background locations. Thus, a buffer layer (18,492 km2) with a radius of 40 km around occurrence points with dissolved boundaries was created by using the R packages “dismo” and “raster” (v3.6-30). The buffer layer is expected to encompass an extent of occurrence (EOO) of Q. acanthophora within the study area (Fig. 2A) (IUCN SSC Amphibian Specialist Group 2017; Nguyen et al. 2024; Pham et al. in press). We randomly sampled 10,000 background points within the buffer area using the “randomPoints” function of the “dismo” package (Hijmans et al. 2022), and this approach can reduce the potential for artificial inflation of AUC values and provides a more realistic evaluation of model discrimination. Next, we applied the Jackknife method to generate a default model. We then used this default model to test model hyperparameters with six feature combinations (L, Q, H, LQ, QH and LQH; where L = linear, Q = quadratic, H = hinge) and regularization values ranging from 0.5 to 5 with a 0.5 increment, using the “gridSearch” function of the “SDMtune” package (v1.3.1) in R. A total of 60 candidate models were evaluated using the AIC test (Akaike Information Criterion), and the optimal model was selected with the lowest AICc value (i.e., delta. AICc = 0) (Warren and Seifert 2011). Based on these filtering criteria, the final SDM for Q. acanthophora was generated with the feature class of LQH and the regularization value of 1.5.

Figure 2. 

A. The potential distribution at present; B. The prediction of multivariate environmental similarity surface (MESS) analysis.

Three types of replication procedures in MaxEnt with bootstrap (100 replicates), cross-validation and subsampling were tested (Phillips et al. 2006). The optimal procedure was selected based on its performance across multiple evaluation metrics of AUC, omission rate, sensitivity, and specificity. Consequently, the replication method of bootstrap was selected for all predictions (Suppl. material 1: table S1). Lastly, we selected the cloglog output format, in which values ranging from 0 to 1 represent cloglog-transformed suitability. The AUC score was used to evaluate the predictive performance, ranging from 0 to 1 (perfect performance). The contribution of each variable to the current climatic model was determined by measuring the permutation importance and contribution percentage (Phillips and Dudík 2008). The minimum training presence cloglog threshold and 10 percentile training presence cloglog threshold were both selected in Maxent to convert continuous predicted values to maps of suitable and highly suitable areas, respectively, which were separated to unsuitable areas (Phillips et al. 2006; Peterson 2011).

To assess the predictive capabilities, multivariate environmental similarity surface (MESS) analyses were further used to quantify areas characterized by potential extrapolation errors (Elith et al. 2010). MESS analysis compares the environmental similarity of any given grid cell projection area to the environmental background, which was used to train the models. We also presented one projected MESS map only for each future climate scenario, by averaging the raster cell values across the three global circulation models (GCMs). In particular, similarity/novelty was classified into three levels: “< 0” (Extrapolation), “=0” (Margin) and “>0” (Interpolation).

Priority conservation areas

The priority conservation areas for Q. acanthophora, which were identified within suitable habitats in terms of climates, were overlaid with the available layer of tropical evergreen rain forests within restricted latitudes of 21°–22°N and longitudes of 106°–108°E (JAXA EORC 2025). For assessing the conservation priorities, information on protected areas was included as shapefiles of nature reserves and national parks in Vietnam, which were further overlaid with highly suitable habitats of Q. acanthophora at present, using the higher threshold of 10% training presence cloglog. All maps were generated using Quantum GIS software (QGIS version 4.0.1), showing the study site, layers of suitable habitats, tropical evergreen rain forests and protected areas.

Results

The optimal model presented a very good prediction, with a mean AUC of 0.975 ± 0.011, a low omission rate (0.087), a high sensitivity (0.913) and specificity (0.936) (Suppl. material 1: table S1). The most important climate variable is the Max Temperature of Warmest Month (Bio5), followed by the Temperature Seasonality (Bio4), together accounting for 89.8% of model contribution and 75.3% of permutation importance (Suppl. material 1: table S3). According to the response curves of the climatic variables, the MaxEnt predicted the potential distribution of Q. acanthophora, in which the temperature of Bio5 ranges from 25.7 °C to 32.8 °C, the seasonal temperature of Bio4 (standard deviation x100) from 453.94 to 602.19 (corresponding to an actual standard deviation of 4.54 °C to 6.02 °C), and other climatic variables of the mean diurnal range (Bio2) from 4.5 °C to 8.23 °C, isothermality (Bio3, quantified as a percentage ratio of Bio 2 to Annual Temperature Range) from 26.42% to 37.86%, the mean temperature of the coldest quarter (Bio11) from 11.28 °C to 17.8 °C, and the precipitation of the driest quarter (Bio17) from 48 mm to 168 mm. These climatic conditions represent the suitable environmental ranges predicted for the species’ potential distribution (Suppl. material 2: fig. S1).

The potential distribution of Q. acanthophora was predicted covering the occurrence regions and surrounding areas in Lang Son, Bac Ninh and Quang Ninh provinces, northeast Vietnam (Fig. 2). According to the prediction, other areas may be prioritized for future field surveys to locate potentially unknown populations of Q. acanthophora: the north of Lang Son, Thai Nguyen and Cao Bang provinces in Vietnam, as well as southern Guangxi Province in China. It is noteworthy that we recorded islands in Hai Phong City and Quang Ninh Province, Vietnam, enveloping suitable conditions for the species (Fig. 2). Using the minimum threshold, the model predicted 5,038.3 km2 of climatically suitable habitat, which covers approximately 1,779.09 km2 (35%) under the 10% training presence threshold (Suppl. material 1: table S4). With the current MESS prediction, interpolated areas mostly envelop the buffer site (Fig. 2C).

Given future projections averaged across three GCMs among periods of the 2050s, 2070s and 2090s for each pathway (e.g., SSP2-4.5, SSP3-7.0 and SSP5-8.5), all climate change scenarios projected a significantly reduced extent of suitable habitats (more than 78%), compared to the current prediction (Figs 3, 4, Suppl. material 1: table S3, Suppl. material 2: fig. S2). Accordingly, the area of suitable habitat under future climate change will only cover a maximum of approximately 1.106.98 km2 (21.97%) in the SSP2-4.5 by the 2050s, and a minimum of 9.61 km2 (0.19%) in the SSP5-8.5 by the 2090s (Suppl. material 1: table S3, Suppl. material 2: fig. S2). MESS comparisons indicated that future suitable habitats retain only a very small overlapping area with climatic conditions analogous to the present (Suppl. material 2: figs S3, S4). In particular, interpolation was reduced to a maximum of 26.77% of the projected suitable area under SSP2-4.5 by the 2050s, whereas under SSP5-8.5 by the 2090s, as much as 96.64% of the analogous conditions within the buffer site were replaced by marginal or extrapolated conditions (Suppl. material 1: table S4). As a forest-specialist amphibian, the species’ suitable habitats will be even more restricted to a small area in all future projections (Fig. 5, Suppl. material 1: table S3).

Figure 3. 

Projected potential distribution for Quasipaa acanthophora under different future climate scenarios.

Figure 4. 

Projected areas of suitable habitats for Quasipaa acanthophora under future climate change scenarios (Yellow: the minimum training presence cloglog threshold, and Green: 10 percentile training presence cloglog threshold).

Figure 5. 

Projections of suitable habitats in the tropical evergreen rain forest for Quasipaa acanthophora under different future climate scenarios (Yellow: Tropical evergreen rain forests with suitable climates; Green: Tropical evergreen rain forests with unsuitable climates; Red: Suitable climate regions outside tropical evergreen rain forests).

Discussion

Climatic conditions

Quasipaa acanthophora only inhabits evergreen forests of northern Vietnam within the tropical monsoon zone, characterized by a seasonal variation in climatic conditions (Dubois and Ohler 2009; Hecht et al. 2013; Nguyen et al. 2025). The MaxEnt model indicated that temperature-related climatic variables, particularly the Max Temperature of Warmest Month (Bio5) during the rainy season in northern Vietnam and Temperature Seasonality (Bio4), play an important role in predicting the potential distribution of Q. acanthophora, accounting for 90% of the contribution. Amphibians, including Q. acanthophora, in tropical forests are highly sensitive to temperature changes and seasonal variability that strongly influence their physiological performance and distribution patterns (Buckley and Jetz 2007; Sunday et al. 2012; Scheffers et al. 2016). Pham et al. (in press) demonstrated that species’ temperature is significantly positively correlated with environmental conditions, including ambient temperatures with a range from 18.6 °C – 28.0 °C across four surveyed months (i.e., March, May, June and November). Humidity along the streams where the species occurs was relatively high, ranging from 68% to 99% (mean ± SE: 86.0 ± 1.2%) (Pham et al. in press). These variables likely reflect the species’ dependence on mild climatic conditions associated with streams in the mountain forests.

Potential distribution

The climatic model predicted suitable habitats, covering the natural distribution of Q. acanthophora across the mountain forests of Bac Ninh, Quang Ninh and Lang Son provinces, northeast Vietnam. Based on the predicted distribution of optimal environmental conditions, forest regions in Cao Bang and Thai Nguyen provinces (Vietnam), as well as Guangxi Province (China) should be prioritized for future field surveys to identify new populations of Q. acanthophora. Based on previous predictions for the potential distribution of Goniurosaurus lichtenfelderi, a tiger gecko species occurring in stream-associated habitats in syntopy with Q. acanthophora, extending potentially to the northeastern Vietnam – southwestern China border region, Ngo et al. (2022) subsequently conducted field surveys and confirmed its presence in the predicted areas.

In addition, islands of Bai Tu Long National Park, Ha Long Bay (Quang Ninh Province), and Cat Ba National Park (Hai Phong City) are suggested as the potential distribution of Q. acanthophora. However, herpetofaunal surveys have not recorded the species in these insular habitats (Gawor et al. 2015; Ngo pers. obs.). Geographic isolation (i.e., mainland and islands) and particular restriction to stream microhabitats likely limit the species’ dispersal and colonization to these islands (Heaney 2007; Whittaker et al. 2008; Lomolino et al. 2010; Pham et al. in press).

Climate change

Future climate change scenarios project a significant contraction in suitable habitats of Q. acanthophora, restricted to a very small region along the China–Vietnam border. Under the most extreme scenario with an increase of the global temperature, approximately 4.7 °C –5.1 °C, during the period 2060–2100, climatically suitable habitats for the species are projected to shrink severely. The pattern of contraction under climate change was noted for all tiger geckos (Goniurosaurus spp.) in northeastern Vietnam as well (Le et al. 2017; Ngo et al. 2021a, 2022, 2025). These projections underscore the potentially negative impacts of climate change not only on Q. acanthophora and Goniurosaurus species, but also on other threatened herpetofauna (including Tylototriton spp. and Shinisaurus crocodilurus) occurring in the study site (Bernardes et al. 2013; Nguyen et al. 2022). A previous study on Asian warty newts (Paramesotriton spp.) in northern Vietnam has shown that climate change can produce contrasting responses among genetic lineages, with the spatial distribution of lineage-specific populations either contracting or expanding, and reduced gene flow increasing vulnerability in certain regions (Tran et al. 2025). Thus, integrating population genetic data with species distribution models should therefore be considered in future studies for the target species of Q. acanthophora.

Addressing model limitations

A total of 20 occurrence records were used to predict the potential distribution of Q. acanthophora and its projected changes under future climate scenarios. However, the number of occurrence records was very limited because the species inhabits a highly restricted distribution (Nguyen et al. 2024; Pham et al. in press), which may reduce the robustness of the model predictions. To improve the accuracy of models, we restricted the background points to a smaller buffer surrounding the occurrence records rather than using the entire study area. This approach reduces the artificial inflation of AUC and specificity values that can arise from the broad background extent and provides a more realistic assessment of model discrimination. We also averaged projections across three global circulation models (GCMs) to reduce uncertainty associated with a single climate model only. Additionally, three MaxEnt replication procedures—bootstrap (100 replicates), cross-validation, and subsampling were compared, and the best-performing procedure was selected for the final model. Finally, we selected the optimal feature class and regularization multiplier through model tuning rather than relying on the default MaxEnt settings.

Conservation

Like many ectothermic species with a restricted distribution, Q. acanthophora is predicted to be particularly vulnerable to climate change (Luedtke et al. 2023; Hayden and Blom 2024; Li et al. 2024; Ngo et al. 2025). Extinction risk to the species may take place, even more rapidly under direct human impacts, such as overharvesting for food consumption and forest degradation (e.g., conversion to agriculture, mining, tourism development) (Pham et al. in press), and such impacts have been recorded in other endemic species in northeast Vietnam as well (e.g., Tylototriton vietnamensis, Goniurosaurus spp., Shinisaurus crocodilurus) (Bernardes et al. 2013; van Schingen et al. 2015, 2016; Ngo et al. 2021b, 2022, 2025). In fact, the current climate change under these impacts has increased the frequency, severity, and damage of natural disasters (e.g., storms, floods) on these species in northeast Vietnam (OCHA 2024; Kieu and Tran 2025; World Bank 2025; Pham et al. in press).

Most amphibian species, particularly narrowly distributed endemics (including Q. acanthophora), are unlikely to cope with climate change by shifting their distribution to newly suitable stream habitats within evergreen forests under future climate scenarios. This limitation is primarily due to narrow distribution ranges, specialized habitats and limited dispersal ability (Menin et al. 2011; Sunday et al. 2011). Furthermore, future models have not incorporated land-use change, yet continued human activities are likely to cause further declines in the natural forest in Vietnam. As a result, even areas projected to remain climatically suitable may no longer provide optimal forest conditions for the target species in the future. The preservation of forest habitats is therefore critical, as they maintain stable microclimatic conditions and buffer Q. acanthophora against broader macroclimatic changes projected under future climate change scenarios. As noted by Ngo et al. (2025), forest protection is also suggested as a conservation measure for the persistence of threatened herpetofauna in northeastern Vietnam, particularly under the combined pressures of future climate change and ongoing anthropogenic impacts.

Identifying core habitats is a crucial step to enhance the effectiveness of conservation measures to protect wild populations and natural habitats of the target species – Q. acanthophora. Based on the current prediction of potential distribution combined with layers of tropical evergreen rain forests and protected areas, we propose two areas, including Tay Yen Tu Nature Reserve (Bac Ninh) and Dong Son–Ky Thuong Nature Reserve (Quang Ninh), serving as the priority conservation areas of Q. acanthophora (Fig. 6).

Figure 6. 

Priority conservation areas for Quasipaa acanthophora.

Conclusion

In this study, species distribution modeling revealed that Q. acanthophora is strongly associated with stable temperature conditions and humid evergreen forest streams in northeastern Vietnam. The model identified several potentially suitable habitats beyond the currently known distribution, including areas along the Vietnam–China border, highlighting priorities for future field surveys. However, future climate scenarios projected a severe contraction of suitable habitats, with the potential complete loss of climatically suitable areas under extreme warming conditions. Combined with ongoing anthropogenic pressures, such as habitat degradation and overexploitation, these findings indicate that Q. acanthophora may face a high risk of population decline in the future.

Therefore, the protection of evergreen forest ecosystems and stream habitats, particularly within key conservation areas such as Tay Yen Tu and Dong Son–Ky Thuong NRs, is critical for the long-term persistence of the species. Population monitoring, habitat protection, assessments of harvesting pressure, forest quality, stream habitat conditions and habitat connectivity among occupied and potentially suitable sites are needed to inform conservation priorities and develop effective management strategies for the target species and other threatened herpetofauna in northeastern Vietnam.

Acknowledgments

We are grateful to the directorates of the Tay Yen Tu and Dong Son-Ky Thuong nature reserves, Forest Protection Departments of Bac Ninh, Lang Son, and Quang Ninh provinces for valuable support for our field work. We thank Do, Q.H., Luong, A.M. and Nguyen, N.T. (IB, Ha Noi) for their assistance in the field. We cordially thank Phi, T.Q. and Le, A.H. (IB, Hanoi), as well as T. Pagel and C. Landsberg (Cologne Zoo) for the support for our research work. The field equipment was provided by the Idea Wild (project code HA_TVIET0126-00).

References

  • Aiello-Lammens ME, Boria RA, Radosavljevic A, Vilela B, Anderson RP (2015) spThin: An R package for spatial thinning of species occurrence records for use in ecological niche models. Ecography 38: 541–545. https://doi.org/10.1111/ecog.01132
  • Bernardes M, Rödder D, Nguyen TT, Pham TC, Nguyen QT, Ziegler T (2013) Habitat characterization and potential distribution of Tylototriton vietnamensis in northern Vietnam. Journal of Natural History 47(17–18): 1161–1175. https://doi.org/10.1080/00222933.2012.743611
  • Buckley LB, Jetz W (2007) Environmental and historical constraints on global patterns of amphibian richness. Proceedings of the Royal Society B: Biological Sciences 274(1614): 1167–1173. https://doi.org/10.1098/rspb.2006.0436
  • Ceballos G, Ehrlich PR, Barnosky AD, García A, Pringle RM, Palmer TM (2015) Accelerated modern human-induced species losses: Entering the sixth mass extinction. Science Advances 1(5): e1400253. https://doi.org/10.1126/sciadv.1400253
  • Chamberlain S, Barve V, McGlinn D, Oldoni D, Desmet P, Geffert L, Ram K (2024) rgbif: Interface to the Global Biodiversity Information Facility API. Version 3.8.4.3. https://CRAN.R-project.org/package=rgbif
  • Deutsch CA, Tewksbury JJ, Huey RB, Sheldon KS, Ghalambor CK, Haak DC, Martin PR (2008) Impacts of climate warming on terrestrial ectotherms across latitude. Proceedings of the National Academy of Sciences of the United States of America 105(18): 6668–6672. https://doi.org/10.1073/pnas.0709472105
  • Dormann CF, Elith J, Bacher S, Buchmann C, Carl G, Carré G, Marquéz JRG, Gruber B, Lafourcade B, Leitão PJ, Münkemüller T, McClean C, Osborne PE, Reineking B, Schröder B, Skidmore AK, Zurell D, Lautenbach S (2013) Collinearity: A review of methods to deal with it and a simulation study evaluating their performance. Ecography 36: 27–46. https://doi.org/10.1111/j.1600-0587.2012.07348.x
  • Eyring V, Bony S, Meehl GA, Senior CA, Stevens B, Stouffer RJ, Taylor KE (2016) Overview of the Coupled Model Intercomparison Project Phase 6 (CMIP6) experimental design and organization. Geoscientific Model Development 9: 1937–1958. https://doi.org/10.5194/gmd-9-1937-2016
  • Fick SE, Hijmans RJ (2017) WorldClim 2: New 1-km spatial resolution climate surfaces for global land areas. International Journal of Climatology 37: 4302–4315. https://doi.org/10.1002/joc.5086
  • Gawor A, Schmitz A, Ziegler T (2015) The herpetofauna of Bai Tu Long National Park, northeastern Vietnam. Salamandra 51(1): 23–44.
  • GBIF.org (2026) GBIF occurrence download. Global Biodiversity Information Facility. https://www.gbif.org
  • Hayden BSI, Blom MPK (2024) Climate change from an ectotherm perspective: Evolutionary consequences and demographic change in amphibian and reptilian populations. Biodiversity and Conservation 33: 905–927. https://doi.org/10.1007/s10531-024-02785-3
  • Hecht VL, Pham TC, Nguyen TT, Nguyen QT, Bonkowski M,, Ziegler T (2013) First report on the herpetofauna of Tay Yen Tu Nature Reserve, northeastern Vietnam. Biodiversity Journal 4(4): 507–522.
  • IPBES (2019) Global assessment report on biodiversity and ecosystem services of the Intergovernmental Science-Policy Platform on Biodiversity and Ecosystem Services. Eds: Brondízio ES, Settele J, Díaz S, Ngo HT. IPBES Secretariat, Bonn, Germany. https://doi.org/10.5281/zenodo.3831673
  • IPCC (2021) Climate Change 2021: The Physical Science Basis. Contribution of Working Group I to the Sixth Assessment Report of the Intergovernmental Panel on Climate Change. Eds: Masson-Delmotte V, Zhai P, Pirani A, Connors SL, Péan C, Berger S, Caud N, Chen Y, Goldfarb L, Gomis MI, Huang M, Leitzell K, Lonnoy E, Matthews JBR, Maycock TK, Waterfield T, Yelekçi O, Yu R, Zhou B. Cambridge University Press, Cambridge. https://doi.org/10.1017/9781009157896
  • Kiedrzynski M, Zielińska KM, Rewicz A, Kiedrzyńska E (2017) Habitat and spatial thinning improve the Maxent models performed within complete data. Journal of Geophysical Research: Biogeosciences 122(6): 1359–1370. https://doi.org/10.1002/2016JG003629
  • Le QT, Le MH, Le XT, Tran AT, Chu TH, Nguyen QT, Ngo NH (2017) Existing status and impact of climate change to the distribution of Goniurosaurus catbaensis. In: Proceedings of the 7th National Scientific Conference on Ecology and Biological Resources. 1034–1040.
  • Lomolino MV, Riddle BR, Whittaker RJ, Brown JH (2010) Biogeography. 4th edn. Sinauer Associates, Sunderland.
  • Luedtke JA, Chanson J, Neam K, Hobin L, Maciel AO, Catenazzi A, Borzée A, Hamidy A, Aowphol A, Jean A, Sosa-Bartuano Á, Fong GA, De Silva A, Fouquet A, Angulo A, Kidov AA, Muñoz Saravia A, Diesmos AC, Tominaga A, Wangyal JT, Rowley JJL, Measey J, Vasudevan K, Chan KO, Gururaja KV, Ovaska K, Warr LC, Canseco-Márquez L, Toledo LF, Díaz LM, Khan MMH, Meegaskumbura M, Acevedo ME, Napoli MF, Ponce MA, Vaira M, Lampo M, Yánez-Muñoz MH, Scherz MD, Rödel MO, Matsui M, Fildor M, Kusrini MD, Ahmed MF, Rais M, Kouamé NG, García N, Gonwouo NL, Burrowes PA, Imbun PY, Wagner P, Kok PJR, Joglar RL, Auguste RJ, Brandão RA, Ibáñez R, Von May R, Hedges SB, Biju SD, Ganesh SR, Wren S, Das S, Flechas SV, Ashpole SL, Robleto-Hernández SJ, Loader SP, Incháustegui SJ, Garg S, Phimmachak S, Richards SJ, Slimani T, Osborne-Naikatini T, Abreu-Jardim TPF, Condez TH, De Carvalho TR, Cutajar TP, Pierson TW, Nguyen TQ (2023) Ongoing declines for the world’s amphibians in the face of emerging threats. Nature 622: 308–314. https://doi.org/10.1038/s41586-023-06578-4
  • Menin M, Waldez F, Lima AP (2011) Effects of environmental and spatial factors on the distribution of anuran species with aquatic reproduction in central Amazonia. Austral Ecology 36: 940–948. https://doi.org/10.1111/j.1442-9993.2010.02226.x
  • Ngo NH, Nguyen QH, Phan QT, Tran MH, Nguyen QT, Ziegler T, Rödder D (2021a) Vulnerability of an endemic tiger gecko (Goniurosaurus huuliensis) to climate change: Modeling environmental refugia and implications for in-situ conservation. Salamandra 57(4): 464–474.
  • Ngo NH, Nguyen QH, Tran MH, Ngo TH, Le DN, Gewiss LR, van Schingen-Khan M, Nguyen QT, Ziegler T (2021b) A morphological and molecular review of the genus Goniurosaurus, including an identification key. European Journal of Taxonomy 751: 38–67. https://doi.org/10.5852/ejt.2021.751.1379
  • Ngo NH, Nguyen QH, Phan QT, Nguyen QT, Gewiss LR, Rödder D, Ziegler T (2022) Modeling the environmental refugia of the endangered Lichtenfelder’s tiger gecko (Goniurosaurus lichtenfelderi) towards implementation of transboundary conservation. Frontiers of Biogeography 14: e51167. https://doi.org/10.21425/F5FBG51167
  • Ngo NH, Nguyen TT, Nguyen HH, Nguyen QT, Nguyen QH, Ngo TTH, Nguyen TTH, Ha TTL, Grismer LL, Rödder D, Ziegler T (2025) Slipped through the narrow door: Effects of climate change on the distribution and conservation of Goniurosaurus luii. Nature Conservation 59: 21–41. https://doi.org/10.3897/natureconservation.59.151689
  • Nguyen THV, Pham VA, Nguyen QT, Tran TT, Ziegler T, Pham TC (2025) New record and diet of the endangered Mau Son spiny frog (Quasipaa acanthophora Dubois & Ohler, 2009) in Vietnam. Biodiversity Data Journal 13: e154204. https://doi.org/10.3897/BDJ.13.e154204
  • Nguyen TT, Ngo TH, Ha QQ, Nguyen QT, Le QT, Nguyen HS, Pham TC, Ziegler T, van Schingen-Khan M, Le DM (2022) Molecular phylogenetic analyses and ecological niche modeling provide new insights into threats to the endangered crocodile lizard (Shinisaurus crocodilurus). Frontiers of Biogeography 14: e54779. https://doi.org/10.21425/F5FBG54779
  • Nguyen QT et al. [Eds] (2024) Vietnam Red Data Book, Part 1. Animals. Publishing House for Science and Technology, Ha Noi, 1250 pp. [In Vietnamese]
  • O’Neill BC, Tebaldi C, van Vuuren DP, Eyring V, Friedlingstein P, Hurtt G, Knutti R, Kriegler E, Lamarque JF, Lowe J, Meehl GA, Moss R, Riahi K, Sanderson BM (2016) The ScenarioMIP project for CMIP6. Geoscientific Model Development 9(9): 3461–3482. https://doi.org/10.5194/gmd-9-3461-2016
  • Pham TC, Phan QT, Nguyen QT, Nguyen THV, Ziegler T, Nguyen VT, Ngo NH (in press) On the edge of survival: Population status, habitat preferences, and conservation challenges of an endangered amphibian – Quasipaa acanthophora in Vietnam. Salamandra.
  • Pimm SL, Jenkins CN, Abell R, Brooks TM, Gittleman JL, Joppa LN, Raven PH, Roberts CM, Sexton JO (2014) The biodiversity of species and their rates of extinction, distribution, and protection. Science 344(6187): 1246752. https://doi.org/10.1126/science.1246752
  • R Core Team (2026) R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna. https://www.r-project.org/
  • Radosavljevic A, Anderson RP (2014) Making better MaxEnt models of species distributions: Complexity, overfitting and evaluation. Journal of Biogeography 41: 629–643. https://doi.org/10.1111/jbi.12227
  • Riahi K, van Vuuren DP, Kriegler E, Edmonds J, O’Neill BC, Fujimori S, Bauer N, Calvin K, Dellink R, Fricko O, Lutz W, Popp A, Cuaresma JC, Leimbach M, Jiang L, Kram T, Rao S, Emmerling J, Ebi K, Hasegawa T, Havlík P, Humpenöder F, Da Silva LA, Smith S, Stehfest E, Bosetti V, Eom J, Gernaat D, Masui T, Rogelj J, Strefler J, Drouet L, Krey V, Luderer G, Harmsen M, Takahashi K, Baumstark L, Doelman J, Kainuma M, Klimont Z, Marangoni G, Lotze-Campen H, Obersteiner M, Tabeau A, Tavoni M (2017) The shared socioeconomic pathways and their energy, land use, and greenhouse gas emissions implications: An overview. Global Environmental Change 42: 153–168. https://doi.org/10.1016/j.gloenvcha.2016.05.009
  • Rowley J, Brown R, Bain R, Kusrini M, Inger R, Stuart B, Wogan G, Neang T, Tanya C, Trung TC, Diesmos A, Iskandar DT, Lau M, Ming LT, Makchai S, Nguyen QT, Phimmachak S (2010) Impending conservation crisis for Southeast Asian amphibians. Biology Letters 6(3): 336–338. https://doi.org/10.1098/rsbl.2009.0793
  • Schär C, Vidale PL, Lüthi D, Frei C, Häberli C, Liniger MA, Appenzeller C (2004) The role of increasing temperature variability in European summer heatwaves. Nature 427: 332–336. https://doi.org/10.1038/nature02300
  • Scheffers BR, De Meester L, Bridge TCL, Hoffmann AA, Pandolfi JM, Corlett RT, Butchart SHM, Pearce-Kelly P, Kovacs KM, Dudgeon D, Pacifici M, Rondinini C, Foden WB, Martin TG, Mora C, Bickford D, Watson JEM (2016) The broad footprint of climate change from genes to biomes to people. Science 354(6313): aaf7671. https://doi.org/10.1126/science.aaf7671
  • Sunday JM, Bates AE, Dulvy NK (2011) Global analysis of thermal tolerance and latitude in ectotherms. Proceedings of the Royal Society B: Biological Sciences 278: 1823–1830. https://doi.org/10.1098/rspb.2010.1295
  • Tatebe H, Ogura T, Nitta T, Komuro Y, Ogochi K, Takemura T, Sudo K, Sekiguchi M, Abe M, Saito F, Chikira M, Watanabe S, Mori M, Hirota N, Kawatani Y, Mochizuki T, Yoshimura K, Takata K, O’ishi R, Yamazaki D, Suzuki T, Kurogi M, Kataoka T, Watanabe M, Kimoto M (2019) Description and basic evaluation of simulated mean state, internal variability, and climate sensitivity in MIROC6. Geoscientific Model Development 12(7): 2727–2765. https://doi.org/10.5194/gmd-12-2727-2019
  • Tran VD, Suzuki T, Fukuyama I, Vera RJ, Nishikawa K (2025) Population genetics provides insights into the impact of future climate change on the genetic structure and distribution of Asian Warty Newts (Genus Paramesotriton). Ecology and Evolution 15(1): e71054. https://doi.org/10.1002/ece3.71054
  • van Schingen M, Pham TC, An TH, Nguyen QT, Bernades M, Bonkowski M, Ziegler T (2015) First ecological assessment of the endangered crocodile lizard, Shinisaurus crocodilurus Ahl, 1930 in Vietnam: Microhabitat characterization and habitat selection. Herpetological Conservation and Biology 10: 948–958.
  • van Schingen M, Ha QQ, Pham TC, Le QT, Nguyen QT, Bonkowski M, Ziegler T (2016) Discovery of a new crocodile lizard population in Vietnam: Population trends, future prognoses and identification of key habitats for conservation. Revue Suisse de Zoologie 123: 241–251. https://doi.org/10.5281/zenodo.155300
  • Warren DL, Seifert SN (2011) Ecological niche modeling in Maxent: The importance of model complexity and the performance of model selection criteria. Ecological Applications 21(2): 335–342. https://doi.org/10.1890/10-1171.1
  • Whittaker RJ, Triantis KA, Ladle RJ (2008) A general dynamic theory of oceanic island biogeography. In: Losos JB, Ricklefs RE (Eds) The Theory of Island Biogeography Revisited. Princeton University Press, Princeton, 88–115. https://doi.org/10.1515/9781400831920.88
  • World Bank (2025) Viet Nam 2045—Growing greener: Pathways to a resilient and sustainable future. World Bank, Washington, DC.
  • Wu T, Lu Y, Fang Y, Xin X, Li L, Li W, Jie W, Zhang F, Liu V, Zhang L, Zhang J, Zhang Y, Wu F, Rong X, Yan J, Wu X, Shi X, Chu M, Wei M, Wang Q, Che J, Fang C, Zhang L, Kang W, Bao J (2019) The Beijing Climate Center Climate System Model (BCC-CSM): the main progress from CMIP5 to CMIP6. Geoscientific Model Development 12(4): 1573–1600. https://doi.org/10.5194/gmd-12-1573-2019
  • Yukimoto S, Kawai H, Koshiro T, Oshima N, Yoshida K, Urakawa S, Hiroyuki T, Deushi M, Tanaka T, Hosaka M, Yabu S, Yoshimura H, Shindo E, Mizuta R, Obata A, Adachi Y, Ishii M (2019) The Meteorological Research Institute Earth System Model version 2.0, MRI-ESM2.0: description and basic evaluation of the physical component. Journal of the Meteorological Society of Japan 97(5): 931–965. https://doi.org/10.2151/jmsj.2019-051

Additional information

Conflict of interest

The authors have declared that no competing interests exist.

Ethical statement

No ethical statement was reported.

Artificial Intelligence (AI) use

The authors accept full responsibility for the content of the manuscript, including the disclosure of any use of AI.

No AI tools were used in the preparation of this manuscript.

Funding

This research was supported by the Vietnam Academy of Science and Technology (Project code THTETN.03/25-26), Ho Chi Minh Communist Youth Union Central Committee (Ben Dam Me – Award) for Ngo, N.H., and Stiftung Artenschutz (Species Conservation Foundation) (Project code 2023-71A).

Author contributions

All authors have contributed equally.

Author ORCIDs

L.T.T. Ha https://orcid.org/0009-0002-5344-6092

C.T. Pham https://orcid.org/0000-0001-5158-4526

T.Q. Phan https://orcid.org/0000-0002-2738-5364

T.Q. Nguyen https://orcid.org/0000-0002-6601-0880

D.V. Tran https://orcid.org/0000-0002-6211-9360

T. Ziegler https://orcid.org/0000-0002-4797-609X

D. Rödder https://orcid.org/0000-0002-6108-1639

H.N. Ngo https://orcid.org/0000-0001-5471-965X

Data availability

All of the data that support the findings of this study are available in the main text or Supplementary Information.

Supplementary materials

Supplementary material 1 

Additional tables

Linh Thuy Thi Ha, Trang Hoang Nguyen, Thinh Van Nguyen, Cuong The Pham, Tien Quang Phan, Truong Quang Nguyen, Dung Van Tran, Thomas Ziegler, Dennis Rödder, Hai Ngoc Ngo

Data type: docx

This dataset is made available under the Open Database License (http://opendatacommons.org/licenses/odbl/1.0/). The Open Database License (ODbL) is a license agreement intended to allow users to freely share, modify, and use this Dataset while maintaining this same freedom for others, provided that the original source and author(s) are credited.
Download file (21.47 kb)
Supplementary material 2 

Additional figures

Linh Thuy Thi Ha, Trang Hoang Nguyen, Thinh Van Nguyen, Cuong The Pham, Tien Quang Phan, Truong Quang Nguyen, Dung Van Tran, Thomas Ziegler, Dennis Rödder, Hai Ngoc Ngo

Data type: docx

This dataset is made available under the Open Database License (http://opendatacommons.org/licenses/odbl/1.0/). The Open Database License (ODbL) is a license agreement intended to allow users to freely share, modify, and use this Dataset while maintaining this same freedom for others, provided that the original source and author(s) are credited.
Download file (2.72 MB)
login to comment