Next Article in Journal
The Impact of Environmental Enrichment on the Cortisol Level of Shelter Cats
Next Article in Special Issue
Analysis of Priority Conservation Areas Using Habitat Quality Models and MaxEnt Models
Previous Article in Journal
Physical and Chemical Characteristics of Droppings as Sensitive Markers of Chicken Health Status
Previous Article in Special Issue
Mitochondrial Variation of Bottlenose Dolphins (Tursiops truncatus) from the Canary Islands Suggests a Key Population for Conservation with High Connectivity within the North-East Atlantic Ocean
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Animals
Volume 14
Issue 9
10.3390/ani14091390
animals-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Prediction of Potential Suitable Distribution Areas for an Endangered Salamander in China

by
Jiacheng Tao
1,
Yifeng Hu
1,
Jianping Jiang
2,
Wanji Yang
3,
Tian Zhao
1,* and
Shengqi Su
1,*
1
College of Fisheries, Southwest University, Chongqing 400715, China
2
CAS Key Laboratory of Mountain Ecological Restoration and Bioresource Utilization & Ecological Restoration Biodiversity Conservation Key Laboratory of Sichuan Province, Chengdu Institute of Biology, Chinese Academy of Sciences, Chengdu 610041, China
3
Shengnongjia National Park Administration, Huibei Provincial Key Laboratory on Conservation Biology of the Shennongjia Golden Snub-Nosed Monkey, Shennongjia 442421, China
*
Authors to whom correspondence should be addressed.
Animals 2024, 14(9), 1390; https://doi.org/10.3390/ani14091390
Submission received: 23 March 2024 / Revised: 20 April 2024 / Accepted: 29 April 2024 / Published: 6 May 2024
(This article belongs to the Special Issue Protecting Endangered Species)
Browse Figures
Review Reports Versions Notes

Abstract

:

Simple Summary

Exploring species’ potential suitable habitats is crucial for endangered species conservation, in particular under future global climate change conditions. The Wushan salamander (Liua shihi) is an endangered salamander in China, which is a national protected species (level II). Based on the distribution records of L. shihi, the main objective of this study was to predict the distribution of suitable habitats under current and future climate conditions for L. shihi. Our results showed that precipitation, cloud density, vegetation type, and ultraviolet radiation were the main environmental factors affecting the distribution of suitable habitat for L. shihi. At present, the suitable habitats are mainly distributed in the Daba Mountain area. Under the future climate conditions, the area of suitable habitats increased, which mainly occurred in central Guizhou Province. These findings provided important information for the conservation of L. shihi.

Abstract

Climate change has been considered to pose critical threats for wildlife. During the past decade, species distribution models were widely used to assess the effects of climate change on the distribution of species’ suitable habitats. Among all the vertebrates, amphibians are most vulnerable to climate change. This is especially true for salamanders, which possess some specific traits such as cutaneous respiration and low vagility. The Wushan salamander (Liua shihi) is a threatened and protected salamander in China, with its wild population decreasing continuously. The main objective of this study was to predict the distribution of suitable habitat for L. shihi using the ENMeval parameter-optimized MaxEnt model under current and future climate conditions. Our results showed that precipitation, cloud density, vegetation type, and ultraviolet radiation were the main environmental factors affecting the distribution of L. shihi. Currently, the suitable habitats for L. shihi are mainly concentrated in the Daba Mountains, including northeastern Chongqing and western Hubei Provinces. Under the future climate conditions, the area of suitable habitats increased, which mainly occurred in central Guizhou Province. This study provided important information for the conservation of L. shihi. Future studies can incorporate more species distribution models to better understand the effects of climate change on the distribution of L. shihi.

1. Introduction

Understanding the effects of human-induced perturbations on biological diversity is one of the central concerns in modern ecology [1]. During the past hundreds of years, human activities have dramatically changed the environment on Earth, in particular the climate, which has strongly affected animals in recent decades [2]. Based on previous studies, over 27% of mammals, 13% of birds, 21% of reptiles, 41% of amphibians, and 37%of fishes are threatened due to climate change and other human disturbances [3]. For instance, although the Atlantic Forest network of protected areas (PAs) supports 18% of the amphibians in South America, the number of amphibian species in PAs is declining under changing climate conditions [4]. The reduction of precipitation will lead to a decrease in the reproductive success rate of shovel-nosed frog (Hemisus marmoratus) and natterjack toad (Bufo calamita), resulting in a decrease in their populations [5]. Climate change can also cause the alteration of amphibian phylogenetic composition and niche. For instance, most of the amphibians in PAs contracted their ranges and such responses are clade specific. Basal amphibian clades (e.g., Gymnophiona and Pipidae) were positively affected by climate change, whereas late-divergent clades (e.g., Cycloramphidae, Centrolenidae, Eleutherodactylidae, Microhylidae) were severely impacted [6]. From the wet season to the dry season, the vertical niche space of amphibians in Sierra Llorona has a clear downward trend in response to natural levels of climate variability [7]. In recent decades, increasing studies also indicated that climate change can lead to the shift of animals’ geographical distribution. For instance, Nottingham et al. showed that the suitable habitats of Del Norte salamander (Plethodon elongatus) and Siskiyou Mountains salamander (P. stormi) will shift to the coast and out of the valley with a move north into the mountains under future climate change conditions in the Pacific Northwest of the United States [8]. Duan et al. demonstrated that amphibians in China would lose 20% of their original distribution ranges on average, and over 90% of species’ suitable habitats will shift to the north when compared with their current distribution range. As a consequence, climate change can lead to significant changes in the spatial pattern of amphibian diversity in China [9].
Among all the vertebrates, amphibians are particularly sensitive to climate change as they cannot regulate their body temperature actively [10]. This is especially true for salamanders, which possess some specific traits such as cutaneous respiration and low vagility [11]. However, studies focused on the effects of climate change on salamanders are still limited (but see [12,13]). The Wushan salamander (Liua shihi) is a national protected (level II) salamander in China, which was classified as Near Threatened in the Red List of China. Although this species was listed as Least Concern (LC) in the International Union for Conservation of Nature (IUCN), the wild populations have declined continuously in recent years [14,15]. Based on the records, this species is widely distributed in montane streams of eastern Sichuan, Chongqing, western Hubei, and southern Shaanxi Provinces, with the elevation ranging from 900 to 2350 m [16], and it mainly feeds on aquatic insects and algae [17]. In recent years, the wild population of L. shihi has been decreasing continuously due to human-induced perturbations [18]. Therefore, it is urgent to understand the distribution of suitable habitats of this species, as well as how the suitable habitat will shift under future climate change conditions.
Species distribution models (SDMs) have been proved to be effective to predict the effects of climate change on species distribution patterns [19]. Based on species distribution points and environmental data, these models predict where species likely inhabit using approaches such as statistical and machine learning analyses [20]. Accordingly, SDMs are involved in several models such as Bioclim, random forest, maximum entropy, regression tree, and genetic algorithm [21,22]. Although none of the above models can be regarded as the best one, the maximum entropy model (MaxEnt) was considered to exhibit higher prediction accuracy, have a stronger ability to integrate multiple environmental variables, and provide more intuitive results [23,24]. Therefore, MaxEnt is increasingly used in ecological studies to investigate the responses of species distribution patterns to climate change. Using this model, Zhao et al. demonstrated that climate change can induce different effects on the evolutionarily significant units (ESUs) of Chinese giant salamander (Andrias davidianus) in China, with the northern ESU exhibiting more severe habitat loss [25]. Moreover, Zank et al. used MaxEnt to investigate the potential effects of climate change on 24 species of red-bellied toads (Melanophryniscus) in South America, and they found that 40% of the species may lose over 50% of their potential distribution area by 2080 [26]. However, most studies only used the default parameters provided by the MaxEnt model, despite the fact that MaxEnt is sensitive to sampling bias and prone to overfitting when using default parameters [23,27]. Therefore, it is essential to optimize the parameters of the MaxEnt before conducting the model analyses [27].
The main objective of the present study was to assess the effects of climate change on the distribution of suitable habitats for L. shihi. Specifically, we (1) investigated the distribution of suitable habitats for L. shihi under current climate conditions; (2) analyzed the key environmental factors affecting the distribution patterns of L. shihi; (3) revealed the shift of suitable habitats (i.e., the distribution patterns and the area) caused by climate change in the future. Based on previous studies (e.g., [28,29]), we predicted that the suitable habitats of L. shihi are mainly distributed in southwestern China at present. We also predicted that climate change will lead to the expansion of suitable habitats from the current distribution area to the southwest. In addition, the area of suitable habitats would decrease due to climate change.

2. Materials and Methods

2.1. Study Area

L. shihi is an endemic amphibian species in China. Although its distribution records were concentrated in the Daba Mountains, its potential suitable habitats could be widely distributed in China. Therefore, and in order to better protect this endangered species, we considered the whole of China as the study area.

2.2. Species Occurrence Data

The occurrence data of L. shihi in this study were obtained from published literature (Table S1), the Global Biodiversity Information Facility website (http://www.gbif.org, accessed on 28 May 2023) (Table S2), and our original field survey (Table S3). In total, 89 occurrence records of L. shihi were collected. To avoid spatial autocorrelation, redundant records within 5 × 5 km grids were excluded using SDMToolbox (version 2.4; [30]). Finally, a total of 53 occurrence records were obtained for further analyses (Figure 1).

2.3. Environmental Variables

Environmental variables were selected based on previous studies demonstrating that they may potentially affect the distribution of amphibians (e.g., [13,31,32]). These variables can be divided into five categories, including bioclimate, meteorology, vegetation, human disturbance, and topography. In total, we obtained 31 environmental variable raster layers (Table 1). Specifically, bioclimatic data were composed of 19 climate factors at a resolution of 2.5 min, which were derived from the WorldClim climate database (http://www.worldclim.org/, accessed on 26 May 2023) [33]. Meteorological factors were composed of ultraviolet-B (UV-B) radiation and cloud cover, which were derived from EarthEnv (https://www.earthenv.org/cloud, accessed on 3 June 2023) and Helmholtz Centre for Environmental Research (https://www.ufz.de/gluv, accessed on 3 June 2023), respectively. Vegetation data contained the percentages of tree coverage and vegetation types, which were from Global Map Data Archives (https://globalmaps.github.io/ptc.html, accessed on 3 June 2023) and Resources and Environmental Science Data Center (https://www.resdc.cn/, accessed on 3 June 2023), respectively. Human disturbance data were represented by the population density, which were downloaded from the Socioeconomic Data and Applications Centre (https://sedac.ciesin.columbia.edu/, accessed on 5 June 2023). Finally, topographic data included elevation, slope, and aspect at a resolution of 90 m, which were obtained from the Geospatial Data Cloud (https://www.gscloud.cn/, accessed on 5 June 2023). We unified their coordinate system as GCS_WGS_1984 and resampled them to obtain a consistent spatial resolution.
In order to reduce the influence of spatial correlation, environmental variables with high correlation but low contribution rate were removed before the model analyses [34]. Correlation analysis was performed using SPSS26.0 software. A Shapiro test was conducted using R software version 4.3.2 (https://www.r-project.org/, accessed on 7 November 2023) before the correlation analysis to identify the distribution of each variable [35]. Variables with a normal distribution were tested by Pearson correlations, and others were tested using Spearman correlations [36]. For the contribution rate, we performed a pre-simulation test in MaxEnt 3.4.4 with the distribution data of L. shihi and the 31 environmental variables. The contribution rate of the variables was tested using the jackknife test [37]. After that, variables with too high correlations (|PCCs| ≥ 0.8; Figure 2) but a low contribution rate (<1%) were removed [38], and the rest of the variables were used for secondary simulation. Based on our results, 15 environmental variables were finally selected for constructing the final models, including five for bioclimate, five for meteorology, two for vegetation, one for human disturbance, and three for topography (Table 1).
The future climate data were obtained from the BCC-CSM2-MR climate system model [39]. These data contained two shared socioeconomic pathways (SSPs), SSP126 and SSP585, which are scenarios of global economic, demographic, and energy development in the future [40]. Specifically, SSP126 represents the combined effects of low vulnerability, mitigating stress, and radiative forcing. SSP585 represents the future socioeconomic path of high-emission, high-carbon (coal, oil, and natural gas) use [41]. In this study, two future climate scenarios (SSP126, SSP585) of three periods (2021–2040, 2041–2060, 2061–2080) were selected for projecting the future distribution area of L. shihi.

2.4. Parameter Optimization and Model Construction

There are five feature types in MaxEnt models, including linear (L), quadratic (Q), hinge (H), product (P), and threshold (T). For parameter adjustment, we computed the AICc values of the modeling parameters’ regularization multiplier (RM) and feature combination (FC; the combination of the above five feature types) in R software using the ENMeval package [42]. In this study, we considered the range of RMs from 0.5 to 4.0 and selected six FC types (i.e., L, LQ, H, LQH, LQHP, and LQHPT). Then, we used the parameters corresponding to the minimum information criterion AICc value to construct the species distribution models [27].
The distribution data, environmental variables, and the optimized model parameters were input into MaxEnt3.4.4 software (New York, NY, USA, https://biodiversityinformatics.amnh.org/open_source/maxent/, accessed on 17 May 2023). The importance of environmental variables to the distribution of L. shihi was evaluated according to the relative contributions of environmental variables and the results of the jackknife test [40]. Twenty-five percent of the distribution data were randomly selected as the test set, while the rest were considered as the training set. The maximum number of background points was 10,000. A total of 10 runs were set for model construction, and the replicated run type was cross-validation.
We used the receiver operating characteristic (ROC) curve and the area under the ROC curve (AUC) to evaluate the accuracy of the model. The range of the AUC values was 0–1. A larger value indicates higher model accuracy, as well as higher credibility of the model. Models can be considered as having high prediction accuracy when the AUC value is greater than 0.8, and then the prediction results can be adopted [43]. An AUC value greater than 0.9 indicates that the prediction accuracy of the model is extremely high [23].

2.5. Parameter Optimization and Model Construction

We imported the average value of MaxEnt output results into ArcGIS 10.8 software and used a conversion tool to convert layers from asc format to raster data. The habitat suitability degree was divided into four levels, including high suitability area, moderate suitability area, low suitability area, and unsuitable area by natural breaks (Jenks) [44]. Finally, we calculated the area and proportion of suitable areas for each level. Moreover, we analyzed the change trend from current to future scenarios.

3. Results

3.1. Model Optimization and Accuracy Evaluation

For the current distribution models, the ΔAICc exhibited the lowest value when feature combination (FC) = LQHP and regulation multiplier (RM) = 2.5, indicating that the model was optimal with these parameters (Figure 3). This best model showed that the AUC value of the working curve of the subjects was 0.992 ± 0.004 (mean ± standard deviation), indicating the extremely high accuracy of the model prediction, thus the overfitting phenomenon could be effectively avoided (Figure 4).
In terms of the future distribution models, the optimal parameters were FC = LQPH and RM = 2 for the SSP585 (2021–2040) scenario, while FC = LQ and RM = 0.5 for the rest of the scenarios. After applying the above parameters in MaxEnt to construct models, the results showed that the AUC values of the working curve of the subjects were all > 0.9.

3.2. The Importance of Environmental Variables

For the MaxEnt models constructed under the current climate scenario, the top five environmental variables accounted for 84.7% of the cumulative contribution, including precipitation of the driest month (Bio14, 29.7%), cloud cover seasonal concentration (Mseason, 28.6%), vegetation type (Veg, 15.8%), mean UV-B of the lowest month (UVB4, 5.6%), and slope (5%; Table 1). In terms of the permutation importance (the extent to which the model depends on the variable; [45]), the top five environmental variables were mean UV-B of the highest month (UVB3, 25.6%), precipitation seasonality (Bio15, 22%), mean temperature of the driest quarter (Bio9, 14%), mean UV-B of the lowest month (UVB4, 13.2%), and precipitation of the driest month (Bio14, 6.9%). For the jackknife test (Figure 5), the test gain value was 3.7 when considering all the environmental variables. When considering the variables individually, precipitation of the driest month (Bio14), vegetation type (Veg), precipitation seasonality (Bio15), cloud cover seasonal concentration (Mseason), and mean UV-B of the lowest month (UVB4) were the top five variables that exhibited the highest test gain values.

3.3. Current Potential Suitable Habitats for L. shihi

Based on the results of MaxEnt models (Figure 6, Table 2), the potential suitable habitat for L. shihi was widely distributed in southwestern China, including Chongqing, Hubei, Sichuan, Shaanxi, Hunan, and Guizhou Provinces. In total, the suitable distribution area under current climate conditions for L. shihi was 45.61 × 104 km2. Specifically, the high-suitability region was mainly concentrated in the Daba Mountains and Shennongjia National Park, which are located at the junction of Chongqing, Hubei, and Shaanxi Provinces. In addition, there were a small number of high-suitability regions scattered in central Sichuan Province. The size of the high-suitability area was 6.51 × 104 km2, accounting for 14.3% of the total suitable habitat. The moderate-suitability region included the eastern part of Sichuan, southern Shaanxi, western Hubei, and eastern Chongqing Provinces, showing a ring shape, and the area was 9.77 × 104 km2, accounting for 21.4% of the total suitable habitats. The low-suitability region was wrapped around the periphery of the moderate- and high-suitability areas, showing a strip shape. Moreover, the area was 29.31 × 104 km2, accounting for 64.3% of the total suitable habitats.

3.4. Future Distribution Patterns of the Suitable Habitats for L. shihi

In the 2021–2040 period under SSP126, the high-suitability area increased to 8.89 × 104 km2, which was mainly contributed by the expansion in central Hubei Province and the junction of Chongqing and Guizhou Provinces. However, the high-suitability habitats in central Sichuan Province disappeared. The area of low-suitability habitats also decreased (244 km2), associated with the loss in central Sichuan Province. In the 2041–2060 period under SSP126, the total area of suitable habitats increased to 62.31 × 104 km2. Specifically, more high-suitability habitats occurred in Guizhou Province, despite the concentrated area in the Daba Mountains decreasing. The moderate-suitability habitats in central Sichuan Province disappeared, while there was no obvious change for low-suitability regions. In the 2061–2080 period under SSP126, a continuous decrease in high-suitability region was observed in the Daba Mountains, with the area being about 8.14 × 104 km2. The moderate-suitability region in Guizhou Province was lost, while the low-suitability region can be only observed in central and south China (e.g., Henan, Hubei, and Guizhou Provinces; Figure 7, Table 2).
In the 2021–2040 period under SSP585, the high-suitability area increased to 8.16 × 104 km2, which was mainly contributed by the expansion in the middle and north of Hubei Province. The moderate-suitability area expanded to the south, mainly located in the east of Chongqing Province and the north of Guizhou Province. The area of low-suitability habitats decreased (2.5 × 104 km2), associated with the loss in the junction of Chongqing and Guizhou Provinces. In the 2041–2060 period under SSP585, the total area of suitable habitats increased to 53.50 × 104 km2. Specifically, the high-suitability area expanded to Guizhou Province. The moderate-suitability area in the south of Henan Province expanded to 16.26 × 104 km2, while the area of low-suitability habitats decreased (0.53 × 104 km2). In the 2061–2080 period under SSP585, the highly suitable areas in the Daba Mountains were more concentrated. There was no obvious shift in the distribution pattern of the low-suitability area (Figure 7, Table 2).

4. Discussion

In the present study, we used optimized MaxEnt models to predict the distribution patterns of suitable habitats for L. shihi in China under current and future climate conditions. Based on the high AUC values, our models can be considered to have high accuracy in prediction [46]. Many previous studies only used the default parameters when conducting MaxEnt models (e.g., [47,48]). However, the default parameters will lead to over-fitting and high omission rates of the model. The “ENMeval” package developed by Muscarella et al. based on the R language has been widely used for optimizing the regularization multiplier (RM) and feature combination (FC) in the MaxEnt model to balance the complexity and avoid those defects [42]. Recently, increasing numbers of researchers have argued that MaxEnt models should be optimized before conducting predictions, as the default parameters may cause some bias [27,49]. Our results supported this claim as we found that the types of FC and the values of RM could change in different models. However, more theoretical work and field work are still needed to verify the effectiveness of parameter optimization in MaxEnt models.
As ectothermic animals, amphibians’ growth and distribution are strongly affected by external environments, in particular the climate conditions [28,29]. This is especially true for salamanders, which are more sensitive to the change in climatic factors [13]. Among all the climatic variables, precipitation of the driest month was the most important one that determined the distribution of potential suitable habitats for L. shihi. Based on previous studies [17], the breeding period for this species is between March and April, associated with the dry season in the Daba Mountains. Therefore, sufficient precipitation can provide suitable spawning sites for L. shihi in montane streams, and permanent streams were critical for them to complete the life cycle [50]. High concentration of cloud cover, woody plant coverage, and low UV-B also were the main environmental variables that affect the distribution of suitable habitats for L. shihi. This is consistent with previous studies showing that ultraviolet light can cause oxidative stress, DNA damage, and egg death in salamanders [51,52]. Since an increase in ultraviolet rays may also lead to dramatic habitat reduction and connectivity fragmentation in other amphibian species that live in montane streams (e.g., spiny-bellied frog: Quasipaa boulengeri; [31]), low UV-B could be an important factor driving the survival and distribution of aquatic amphibians. In the present study, it was found that high cloud density and forest coverage can effectively reduce the damage of ultraviolet rays to L. shihi [53], supporting the survival and distribution of this species. In addition, the influence of slope cannot be ignored, which was associated with the water flow rate and sunshine angle of the habitat, and salamanders usually preferred to select places with low water flow rate and sufficient light to grow and reproduce [54,55].
Our results showed that the high-suitability habitat for L. shihi was concentrated in the junction of Chongqing, Shaanxi, and Hebei Provinces, suggesting that this species may prefer some specific ecological conditions in this area [56]. Therefore, this region should be paid more attention for the protection of this species. For protected animals, a concentrated distribution pattern means they may be more easily threatened by climate change, and regional natural disasters and disease transmission will put the entire population at risk of extinction [57]. Interestingly, there were no distribution points recorded in some high-suitability regions (e.g., central Sichuan and central Chongqing Provinces), indicating that further field investigations can be carried out in these areas. In addition, a small number of existing distribution points were located in low-suitability or even non-suitable areas, suggesting that these populations should be paid more attention.
In the future, the total area of the suitable habitats for L. shihi will increase, although the main spatial distribution patterns did not change dramatically. This may be due to the unique climatic conditions (cool and humid all year round) in the Daba Mountains and Shennongjia National Park, which are climate transition regions between subtropical and northern warm temperate zones [58,59]. In two periods (2021–2040 and 2061–2080), the area of the suitable habitats under SSP126 was smaller than that of SSP585. This shows that the high-emission and high-carbon use scenario (SSP585) may cause an increase in the area of suitable habitat for L. shihi, which is similar to the finding of Wider et al. showing that the suitable range of the blue-spotted salamander (Ambystoma laterale) and the red-backed salamander (P. cinereus) increases with the increase in greenhouse gas concentration [60]. From the time point of view, the suitable area of L. shihi in the future is larger than the current results, and the high-suitability area under the SSP585 scenario will gradually increase with time. This increase may be the cumulative effect of climate change. This is contrary to previous studies showing that the area of suitable habitats of some other salamanders (e.g., leprous false brook salamander: Pseudoeurycea leprosa, streamside salamander: A. barbouri, and Cheat Mountain salamander: P. nettingi) will significantly decrease in the future [45,61]. We speculated that under this scenario, climate change has just reached the suitable conditions for L. shihi in some areas. It is worth noting that in the next three periods, the distribution range will be more concentrated. It indicates that the concentrated areas may have more important protection significance, as this region should be the refuge for L. shihi under future climate change.

5. Conclusions

In conclusion, the present study predicted the potential suitable habitats for L. shihi using a MaxEnt model with optimized parameters under current and future climate change scenarios for three time periods (SSP126 and SSP585). Our results indicated that precipitation of the driest month (Bio14), cloud cover seasonal concentration (Mseason), vegetation type (Veg), mean UV-B of the lowest month (UVB4), and slope are important environmental variables that have a great impact on the habitat suitability. The suitable habitats under the current situation are mainly distributed at the junction of Chongqing, Shaanxi, and Hubei Provinces. Under future climatic conditions, the total suitable area increased. The new suitable habitats were concentrated in the central part of Guizhou and Hubei Provinces. However, suitable habitats located in the central part of Sichuan and Chongqing Provinces were lost. The results of this study can help us better understand the distribution of L. shihi and can provide important information for determining the suitable areas of this species in China. Future studies can incorporate more species distribution models to better understand the effects of climate change on the distribution of suitable habitats for L. shihi.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ani14091390/s1, Table S1: The distribution data of Liua shihi collected in published literatures; Table S2: The distribution data of Liua shihi collected on website; and Table S3: The distribution data of Liua shihi collected from field survey.

Author Contributions

Conceptualization, J.T. and T.Z.; data curation, J.T., Y.H., W.Y., and J.J.; formal analysis, J.T. and Y.H.; methodology, J.T. and Y.H.; resources, J.T. and Y.H.; software, J.T. and Y.H.; funding acquisition, J.J., T.Z. and S.S.; writing—original draft, J.T.; supervision, T.Z. and S.S.; validation, T.Z. and S.S.; visualization, W.Y.; writing—review and editing, J.T., T.Z. and S.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Project of Background Resources Survey in Shennongjia National Park (SNJNP2022007), the Open Project Fund of Hubei Provincial Key Laboratory for Conservation Biology of Shennongjia Snub-nosed Monkeys (SNJGKL2022007), and the Chongqing Aquatic Science and Technology Innovation Alliance Research Program (CQFTIU2022-18).

Institutional Review Board Statement

Animal experiments were performed following the Guide for Care and Use of Laboratory Animals approved by the Committee of Laboratory Animal Experimentation of Southwest University (SWU2011032201, 22 March 2011).

Informed Consent Statement

Not applicable.

Data Availability Statement

The datasets presented in this study are available in Table S1.

Acknowledgments

We thank Chunlin Zhao and Da Kang for their help in constructing the MaxEnt models. We also thank Da Kang, Tianpeng Huang, and Jun Jiang for their help in collecting the distribution data of L. shihi.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Marzin, A.; Archaimbault, V.; Belliard, J.; Chauvin, C.; Delmas, F.; Pont, D. Ecological assessment of running waters: Do macrophytes, macroinvertebrates, diatoms and fish show similar responses to human pressures? Ecol. Indic. 2012, 23, 56–65. [Google Scholar] [CrossRef]
  2. Pacifici, M.; Foden, W.B.; Visconti, P.; Watson, J.E.M.; Butchart, S.H.M.; Kovacs, K.M.; Scheffers, B.R.; Hole, D.G.; Martin, T.G.; Akcakaya, H.R.; et al. Assessing species vulnerability to climate change. Nat. Clim. Change 2015, 5, 215–225. [Google Scholar] [CrossRef]
  3. Luedtke, J.A.; Chanson, J.; Neam, K.; Hobin, L.; Maciel, A.O.; Catenazzi, A.; Borzee, A.; Hamidy, A.; Aowphol, A.; Jean, A.; et al. Ongoing declines for the world’s amphibians in the face of emerging threats. Nature 2023, 622, 308. [Google Scholar] [CrossRef] [PubMed]
  4. Lemes, P.; Melo, A.S.; Loyola, R.D. Climate change threatens protected areas of the Atlantic Forest. Biodivers. Conserv. 2014, 23, 357–368. [Google Scholar] [CrossRef]
  5. Banks, B.; Beebee, T.; Cooke, A.S. Conservation of the natterjack toad Bufo-calamita in britain over the period 1970–1990 in relation to site protection and other factors. Biol. Conserv. 1994, 67, 111–118. [Google Scholar] [CrossRef]
  6. Loyola, R.D.; Lemes, P.; Brum, F.T.; Provete, D.B.; Duarte, L.D.S. Clade-specific consequences of climate change to amphibians in Atlantic Forest protected areas. Ecography 2014, 37, 65–72. [Google Scholar] [CrossRef]
  7. Basham, E.W.; Scheffers, B.R. Vertical stratification collapses under seasonal shifts in climate. J. Biogeogr. 2020, 47, 1888–1898. [Google Scholar] [CrossRef]
  8. Nottingham, S.; Pelletier, T.A. The impact of climate change on western Plethodon salamanders’ distribution. Ecol. Evol. 2021, 11, 9370–9384. [Google Scholar] [CrossRef] [PubMed]
  9. Duan, R.; Kong, X.; Huang, M.; Varela, S.; Ji, X. The potential effects of climate change on amphibian distribution, range fragmentation and turnover in China. Peerj 2016, 4, e2185. [Google Scholar] [CrossRef]
  10. Baskiera, S.; Gvozdik, L. Thermal independence of energy management in a tailed amphibian. J. Vertebr. Biol. 2020, 69, 20057.1-10. [Google Scholar] [CrossRef]
  11. Barrett, K.; Guyer, C. Differential responses of amphibians and reptiles in riparian and stream habitats to land use disturbances in western Georgia, USA. Biol. Conserv. 2008, 141, 2290–2300. [Google Scholar] [CrossRef]
  12. Ashrafzadeh, M.R.; Naghipour, A.A.; Haidarian, M.; Kusza, S.; Pilliod, D.S. Effects of climate change on habitat and connectivity for populations of a vulnerable, endemic salamander in Iran. Glob. Ecol. Conserv. 2019, 19, e00637. [Google Scholar] [CrossRef]
  13. Zhang, P.; Dong, X.; Grenouillet, G.; Lek, S.; Zheng, Y.; Chang, J. Species range shifts in response to climate change and human pressure for the world’s largest amphibian. Sci. Total Environ. 2020, 735, 139543. [Google Scholar] [CrossRef]
  14. IUCN Red List of Threatened Species. Available online: https://www.iucnredlist.org/en (accessed on 21 June 2019).
  15. Redlist of China’s Biodiversity. Available online: https://www.mee.gov.cn/xxgk2018/xxgk/xxgk01/202305/W020230522536559098623.pdf (accessed on 16 July 2023).
  16. McDonough, K. Amphibian Species of the World: An Online Reference (Version 6). Ref. Rev. 2014, 28, 32. [Google Scholar] [CrossRef]
  17. Fei, L.; Ye, C.; Jiang, J. Colorful Illustrations of Amphibians in China, 1st ed.; Sichuan Publishing Group Sichuan Science and Technology Press: Chengdu, China, 2010; pp. 66–67. (In Chinese) [Google Scholar]
  18. Ma, Q.; Wan, L.; Shi, S.; Wang, Z. Impact of Climate Change on the Distribution of Three Rare Salamanders (Liua shihi, Pseudohynobius jinfo, and Tylototriton wenxianensis) in Chongqing, China, and Their Conservation Implications. Animals 2024, 14, 672. [Google Scholar] [CrossRef] [PubMed]
  19. Jimenez-Garcia, D.; Peterson, A.T. Climate change impact on endangered cloud forest tree species in Mexico. Rev. Mex. Biodivers. 2019, 90, e902781. [Google Scholar] [CrossRef]
  20. Li, G.; Liu, C.; Liu, Y.; Yang, J.; Zhang, X.; Guo, K. Progress in theoretical research of species distribution models. Acta Ecol. Sin. 2013, 33, 4827–4835, (In Chinese with English Abstract). [Google Scholar]
  21. Gama, M.; Crespo, D.; Dolbeth, M.; Anastacio, P. Predicting global habitat suitability for Corbicula fluminea using species distribution models: The importance of different environmental datasets. Ecol. Model. 2016, 319, 163–169. [Google Scholar] [CrossRef]
  22. Padalia, I.; Srivastava, V.; Kushwaha, S.P.S. Modeling potential invasion range of alien invasive species, Hyptis suaveolens (L.) Poit. in India: Comparison of MaxEnt and GARP. Ecol. Inform. 2014, 22, 36–43. [Google Scholar] [CrossRef]
  23. Phillips, S.J.; Anderson, R.P.; Schapire, R.E. Maximum entropy modeling of species geographic distributions. Ecol. Model. 2006, 190, 231–259. [Google Scholar] [CrossRef]
  24. Phillips, S.J.; Dudik, M. Modeling of species distributions with Maxent:: New extensions and a comprehensive evaluation. Ecography 2008, 31, 161–175. [Google Scholar] [CrossRef]
  25. Zhao, T.; Zhang, W.; Zhou, J.; Zhao, C.; Liu, X.; Liu, Z.; Shu, G.; Wang, S.; Li, C.; Xie, F.; et al. Niche divergence of evolutionarily significant units with implications for repopulation programs of the world’s largest amphibians. Sci. Total Environ. 2020, 738, 140269. [Google Scholar] [CrossRef] [PubMed]
  26. Zank, C.; Becker, F.G.; Abadie, M.; Baldo, D.; Maneyro, R.; Borges-Martins, M. Climate Change and the Distribution of Neotropical RedBellied Toads (Melanophryniscus, Anura, Amphibia): How to Prioritize Species and Populations? PLoS ONE 2014, 9, e94625. [Google Scholar] [CrossRef] [PubMed]
  27. Zhu, G.; Qiao, H. Effect of the Maxent model’s complexity on the prediction of species potential distributions. Biodivers. Sci. 2016, 24, 1189–1196, (In Chinese with English Abstract). [Google Scholar] [CrossRef]
  28. Wu, J. The hazard and unsureness of reducing habitat ranges in response to climate warming for 91 amphibian species in China. Acta Oecol 2020, 108, 103640. [Google Scholar] [CrossRef]
  29. Zhang, Z.; Mammola, S.; Liang, Z.; Capinha, C.; Wei, Q.; Wu, Y.; Zhou, J.; Wang, C. Future climate change will severely reduce habitat suitability of the Critically Endangered Chinese giant salamander. Freshwater Biol. 2020, 65, 971–980. [Google Scholar] [CrossRef]
  30. Brown, J.L.; Bennett, J.R.; French, C.M. SDMtoolbox 2.0: The next generation Python-based GIS toolkit for landscape genetic, biogeographic and species distribution model analyses. Peerj 2017, 5, e4095. [Google Scholar] [CrossRef]
  31. Yang, S.; Wang, X.; Hu, J. Mountain frog species losing out to climate change around the Sichuan Basin. Sci. Total Environ. 2022, 806, 150605. [Google Scholar] [CrossRef] [PubMed]
  32. Jiang, D.; Zhao, X.; Lopez-Pujol, J.; Wang, Z.; Qu, Y.; Zhang, Y.; Zhang, T.; Li, D.; Jiang, K.; Wang, B.; et al. Effects of climate change and anthropogenic activity on ranges of vertebrate species endemic to the Qinghai-Tibet Plateau over 40 years. Conserv. Biol. 2023, 37, e14069. [Google Scholar] [CrossRef]
  33. Hijmans, R.J.; Cameron, S.E.; Parra, J.L.; Jones, P.G.; Jarvis, A. Very high resolution interpolated climate surfaces for global land areas. Int. J. Climatol. 2005, 25, 1965–1978. [Google Scholar] [CrossRef]
  34. Harisena, N.V.; Groen, T.A.; Toxopeus, A.G.; Naimi, B. When is variable importance estimation in species distribution modelling affected by spatial correlation? Ecography 2021, 44, 778–788. [Google Scholar] [CrossRef]
  35. Aslam, M. Analysing Gray Cast Iron Data using a New Shapiro-Wilks test for Normality under Indeterminacy. Int. J. Cast. Metal. Res. 2021, 34, 1–5. [Google Scholar] [CrossRef]
  36. Rovetta, A. Raiders of the Lost Correlation: A Guide on Using Pearson and Spearman Coefficients to Detect Hidden Correlations in Medical Sciences. Cureus J. Med. Sci. 2020, 12, e11794. [Google Scholar] [CrossRef] [PubMed]
  37. Hosni, E.M.; Nasser, M.G.; Al-Ashaal, S.A.; Rady, M.H.; Kenawy, M.A. Modeling current and future global distribution of Chrysomya bezziana under changing climate. Sci. Rep. 2020, 10, 4947. [Google Scholar] [CrossRef] [PubMed]
  38. Yang, X.; Kushwaha, S.P.S.; Saran, S.; Xu, J.; Roy, P.S. Maxent modeling for predicting the potential distribution of medicinal plant, Justicia adhatoda L. in Lesser Himalayan foothills. Ecol. Eng. 2013, 51, 83–87. [Google Scholar] [CrossRef]
  39. Fick, S.E.; Hijmans, R.J. WorldClim 2: New 1-km spatial resolution climate surfaces for global land areas. Int. J. Climatol. 2017, 37, 4302–4315. [Google Scholar] [CrossRef]
  40. Jiang, C.; Zhang, X.; Xie, W.; Wang, R.; Feng, C.; Ma, L.; Li, Q.; Yang, Q.; Wang, H. Predicting the potential distribution of the fall armyworm Spodoptera frugiperda (JE Smith) under climate change in China. Glob. Ecol. Conserv. 2022, 33, e01994. [Google Scholar] [CrossRef]
  41. Riahi, K.; van Vuuren, D.P.; Kriegler, E.; Edmonds, J.; O’Neill, B.C.; Fujimori, S.; Bauer, N.; Calvin, K.; Dellink, R.; Fricko, O.; et al. The Shared Socioeconomic Pathways and their energy, land use, and greenhouse gas emissions implications: An overview. Global Environ. Chang. 2017, 42, 153–168. [Google Scholar] [CrossRef]
  42. Muscarella, R.; Galante, P.J.; Soley-Guardia, M.; Boria, R.A.; Kass, J.M.; Uriarte, M.; Anderson, R.P. ENMeval: An R package for conducting spatially independent evaluations and estimating optimal model complexity for MAXENT ecological niche models. Methods Ecol. Evol. 2014, 5, 1198–1205. [Google Scholar] [CrossRef]
  43. Fu, A.; Gao, E.; Tang, X.; Liu, Z.; Hu, F.; Zhan, Z.; Wang, J.; Luan, X. MaxEnt Modeling for Predicting the Potential Wintering Distribution of Eurasian Spoonbill (Platalea leucorodia leucorodia) under Climate Change in China. Animals-Basel 2023, 13, 856. [Google Scholar] [CrossRef]
  44. Zhang, Y.; Hughes, A.C.; Zhao, Z.; Li, Z.; Qin, Y. Including climate change to predict the global suitable area of an invasive pest: Bactrocera correcta (Diptera: Tephritidae). Glob. Ecol. Conserv. 2022, 34, e02021. [Google Scholar] [CrossRef]
  45. Sutton, W.B.; Barrett, K.; Moody, A.T.; Loftin, C.S.; DeMaynadier, P.G.; Nanjappa, P. Predicted Changes in Climatic Niche and Climate Refugia of Conservation Priority Salamander Species in the Northeastern United States. Forests 2015, 6, 1–26. [Google Scholar] [CrossRef]
  46. Bartlow, A.W.; Giermakowski, J.T.; Painter, C.W.; Neville, P.; Schultz-Fellenz, E.S.; Crawford, B.M.; Lavadie-Bulnes, A.F.; Thompson, B.E.; Hathcock, C.D. Modeling the distribution of the endangered Jemez Mountains salamander (Plethodon neomexicanus) in relation to geology, topography, and climate. Ecol. Evol. 2022, 12, e9161. [Google Scholar] [CrossRef] [PubMed]
  47. Li, R.; Xu, M.; Wong, M.H.G.; Qiu, S.; Li, X.; Ehrenfeld, D.; Li, D. Climate change threatens giant panda protection in the 21st century. Biol. Conserv. 2015, 182, 93–101. [Google Scholar] [CrossRef]
  48. Stickley, S.F.; Fraterrigo, J.M. Microclimate species distribution models estimate lower levels of climate-related habitat loss for salamanders. J. Nat. Conserv. 2023, 72, 126333. [Google Scholar] [CrossRef]
  49. Karuppaiah, V.; Maruthadurai, R.; Das, B.; Soumia, P.S.; Gadge, A.S.; Thangasamy, A.; Ramesh, S.V.; Shirsat, D.V.; Mahajan, V.; Krishna, H.; et al. Predicting the potential geographical distribution of onion thrips, Thrips tabaci in India based on climate change projections using MaxEnt. Sci. Rep. 2023, 13, 7934. [Google Scholar] [CrossRef] [PubMed]
  50. Connette, G.M.; Osbourn, M.S.; Peterman, W.E. The Distribution of a Stream-breeding Salamander, Desmognathus ocoee, in Terrestrial Habitat Suggests the Ecological Importance of Low-order Streams. Copeia 2016, 104, 149–156. [Google Scholar] [CrossRef]
  51. Lesser, M.P.; Turtle, S.L.; Farrell, J.H.; Walker, C.W. Exposure to ultraviolet radiation (290-400 nm) causes oxidative stress, DNA damage, and expression of p53/p73 in laboratory experiments on embryos of the spotted salamander, Ambystoma maculatum. Physiol. Biochem. Zool. 2001, 74, 733–741. [Google Scholar] [CrossRef] [PubMed]
  52. Blaustein, A.R.; Edmond, B.; Kiesecker, J.M.; Beatty, J.J.; Hokit, D.G. Ambient ultraviolet-radiation causes mortality in salamander eggs. Ecol. Appl. 1995, 5, 740–743. [Google Scholar] [CrossRef]
  53. Xu, R.; Li, Y.; Teuling, A.J.; Zhao, L.; Spracklen, D.V.; Garcia-Carreras, L.; Meier, R.; Chen, L.; Zheng, Y.; Lin, H.; et al. Contrasting impacts of forests on cloud cover based on satellite observations. Nat. Commun. 2022, 13, 670. [Google Scholar] [CrossRef]
  54. Huang, H.; Bu, R.; Xie, H.; Hou, S.; Wu, Z. Habitat selection by Hynobius maoershanensis during its breeding period. Acta Ecol. Sin. 2019, 39, 6443–6451, (In Chinese with English Abstract). [Google Scholar]
  55. Zhang, Q.; Fang, R.; Deng, C.; Zhao, H.; Shen, M.H.; Wang, Q. Slope aspect effects on plant community characteristics and soil properties of alpine meadows on Eastern Qinghai-Tibetan plateau. Ecol. Indic. 2022, 143, 109400. [Google Scholar] [CrossRef]
  56. Siemens, A.L.; Bogart, J.P.; Linton, J.E.; Norris, D.R. Predicting the occurrence of an endangered salamander in a highly urbanized landscape. Endanger. Species Res. 2023, 52, 81–95. [Google Scholar] [CrossRef]
  57. Bakare, A.G.; Kour, G.; Akter, M.; Iji, P.A. Impact of climate change on sustainable livestock production and existence of wildlife and marine species in the South Pacific island countries: A review. Int. J. Biometeorol. 2020, 64, 1409–1421. [Google Scholar] [CrossRef] [PubMed]
  58. Deng, Z.; Cao, J. Incorporating ecosystem services into functional zoning and adaptive management of natural protected areas as case study of the Shennongjia region in China. Sci. Rep. 2023, 13, 18870. [Google Scholar] [CrossRef] [PubMed]
  59. Yao, Y.; Hu, Y.; Kou, Z.; Zhang, B. Spatial patterns of Pinus tabulaeformis and Pinus massoniana forests in Qinling-Daba Mountains and the boundary of subtropical and warm temperate zones. J. Geogr. Sci. 2020, 30, 1523–1533. [Google Scholar] [CrossRef]
  60. Widmer, B.W.; Gehring, T.M.; Heumann, B.W.; Nicholson, K.E. Climate change and range restriction of common salamanders in eastern Canada and the United States. J. Wildlife Manage 2022, 86, e22235. [Google Scholar] [CrossRef]
  61. Parra-Olea, G.; Martínez-Meyer, E.; de León, G. Forecasting climate change effects on salamander distribution in the highlands of central Mexico. Biotropica 2005, 37, 202–208. [Google Scholar] [CrossRef]
Animals 14 01390 g001
Figure 1. Distribution points of L. shihi (red triangles; after excluding autocorrelation).
Figure 1. Distribution points of L. shihi (red triangles; after excluding autocorrelation).
Animals 14 01390 g001
Animals 14 01390 g002
Figure 2. Correlation matrix between environmental variables. Bio1–19 are bioclimatic variables obtained from WorldClim website. Asp: aspect; Ele: elevation; Slope: slope; Pdensity: density of population; Mmean: mean annual cloud cover; Mseason: cloud cover seasonal concentration; Mspatial: cloud cover spatial variability; UVB1: annual mean UV-B; UVB3: mean UV-B of highest month; UVB4: mean UV-B of lowest month; Plantcover: density of trees on the ground; Veg: vegetation type. Positive correlations are displayed in blue and negative correlations in a red color. The color intensity and the size of the circle are proportional to the correlation coefficients.
Figure 2. Correlation matrix between environmental variables. Bio1–19 are bioclimatic variables obtained from WorldClim website. Asp: aspect; Ele: elevation; Slope: slope; Pdensity: density of population; Mmean: mean annual cloud cover; Mseason: cloud cover seasonal concentration; Mspatial: cloud cover spatial variability; UVB1: annual mean UV-B; UVB3: mean UV-B of highest month; UVB4: mean UV-B of lowest month; Plantcover: density of trees on the ground; Veg: vegetation type. Positive correlations are displayed in blue and negative correlations in a red color. The color intensity and the size of the circle are proportional to the correlation coefficients.
Animals 14 01390 g002
Animals 14 01390 g003
Figure 3. ΔAICc of the MaxEnt models under different regularization multipliers (RMs) and feature combinations (FCs).
Figure 3. ΔAICc of the MaxEnt models under different regularization multipliers (RMs) and feature combinations (FCs).
Animals 14 01390 g003
Animals 14 01390 g004
Figure 4. Receiver operating characteristic (ROC) curve and AUC value.
Figure 4. Receiver operating characteristic (ROC) curve and AUC value.
Animals 14 01390 g004
Animals 14 01390 g005
Figure 5. Jackknife of test gain for environmental variables in L. shihi.
Figure 5. Jackknife of test gain for environmental variables in L. shihi.
Animals 14 01390 g005
Animals 14 01390 g006
Figure 6. Potential suitable habitat for L. shihi under current climatic conditions.
Figure 6. Potential suitable habitat for L. shihi under current climatic conditions.
Animals 14 01390 g006
Animals 14 01390 g007
Figure 7. Potentially suitable climatic distribution of L. shihi under different climate change scenarios: (a) SSP126 from 2021–2040; (b) SSP126 from 2041–2060; (c) SSP126 from 2061–2082; (d) SSP585 from 2021–2040; (e) SSP585 from 2041–2060; (f) SSP585 from 2061–2080.
Figure 7. Potentially suitable climatic distribution of L. shihi under different climate change scenarios: (a) SSP126 from 2021–2040; (b) SSP126 from 2041–2060; (c) SSP126 from 2061–2082; (d) SSP585 from 2021–2040; (e) SSP585 from 2041–2060; (f) SSP585 from 2061–2080.
Animals 14 01390 g007
Table 1. Contribution and permutation importance of environmental variables in MaxEnt models.
Table 1. Contribution and permutation importance of environmental variables in MaxEnt models.
CodeEnvirnonmental VariablePercentage Contribution (%)Permutation Importance (%)
Bio14Precipitation of driest month29.76.9
MseasonCloud cover seasonal concentration28.66.6
VegVegetation type15.80.5
UVB4Mean UV-B of lowest month5.613.2
SlopeSlope55
Bio2Mean diurnal range3.50.6
MspatialCloud cover spatial variability3.51.6
Bio9Mean temperature of driest quarter3.114
Bio15Precipitation seasonality 1.622
Bio4Temperature seasonality1.42.5
UVB3Mean UV-B of highest month1.225.6
PdensityDensity of population0.41.5
UVB1Annual mean UV-B0.40
AspAspect00.1
PlantcoverDensity of trees on the ground00
Table 2. The area of suitable habitats in different periods (values in the brackets indicate the variations when compared with the area of the current period) (×104 km2).
Table 2. The area of suitable habitats in different periods (values in the brackets indicate the variations when compared with the area of the current period) (×104 km2).
GradeCurrent2021–20402041–20602061–2080
SSP126SSP585SSP126SSP585SSP126SSP585
Low29.30629.281
(−0.025)		26.800
(−2.506)		36.037
(6.731)		28.774
(−0.532)		30.802
(1.496)		35.229
(5.923)
Moderate9.77416.646
(6.872)		20.535
(10.761)		17.336
(7.562)		16.262
(6.488)		9.889
(0.115)		19.076
(9.302)
High6.5268.887
(2.361)		8.162
(1.636)		8.931
(−2.405)		8.457
(1.931)		8.137
(1.611)		9.307
(2.781)
Total45.60954.814
(9.205)		55.497
(9.888)		62.306
(16.697)		53.493
(7.884)		48.828
(3.219)		63.613
(18.004)
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Tao, J.; Hu, Y.; Jiang, J.; Yang, W.; Zhao, T.; Su, S. Prediction of Potential Suitable Distribution Areas for an Endangered Salamander in China. Animals 2024, 14, 1390. https://doi.org/10.3390/ani14091390

AMA Style

Tao J, Hu Y, Jiang J, Yang W, Zhao T, Su S. Prediction of Potential Suitable Distribution Areas for an Endangered Salamander in China. Animals. 2024; 14(9):1390. https://doi.org/10.3390/ani14091390

Chicago/Turabian Style

Tao, Jiacheng, Yifeng Hu, Jianping Jiang, Wanji Yang, Tian Zhao, and Shengqi Su. 2024. "Prediction of Potential Suitable Distribution Areas for an Endangered Salamander in China" Animals 14, no. 9: 1390. https://doi.org/10.3390/ani14091390

APA Style

Tao, J., Hu, Y., Jiang, J., Yang, W., Zhao, T., & Su, S. (2024). Prediction of Potential Suitable Distribution Areas for an Endangered Salamander in China. Animals, 14(9), 1390. https://doi.org/10.3390/ani14091390

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop