Biodiversity Science ›› 2018, Vol. 26 ›› Issue (9): 941-950.doi: 10.17520/biods.2018125

• Original Papers • Previous Article     Next Article

Amphibian species richness patterns in karst regions in Southwest China and its environmental associations

Bo Wang1, 4, Yong Huang2, Jiatang Li3, Qiang Dai3, Yuezhao Wang3, Daode Yang1, *()   

  1. 1 Institute of Wildlife Conservation, Central South University of Forestry and Technology, Changsha 410004
    2 Guangxi University of Chinese Medicine, Nanning 530200
    3 Chengdu Institute of Biology, Chinese Academy of Sciences, Chengdu 610041
    4 Guangxi Lujin Ecological Technology Company, Nanning 530028
  • Received:2018-04-22 Accepted:2018-07-05 Online:2019-01-05
  • Yang Daode
  • About author:

    # Co-first authors

Patterns in the distribution of species richness have always been a central theme in macroecology. The karst landforms in Southwest China (mainly Guangxi, Yunnan and Guizhou provinces) are among the largest of the global biodiversity hotspots. In this study, we sought to understand spatial patterns of amphibian species richness and its relationship with environmental factors. We compiled a large dataset of 18,246 records of point location data for 219 amphibian species occurring in China. We retrieved this data from published literature, Herpetology museums of Chengdu Institute of Biology and Kunming Institute of Zoology, Chinese Academy of Sciences, Guangxi Zhuang Autonomous Region Museum of Nature and the Central South University of Forestry and Technology, and published sources. We used this data to generate the potential distributions of each species using ecological niche modeling. We combined the potential distributions maps of all species into a composite map to describe species richness patterns on the grid cell of 10 km × 10 km, and then conducted multivariate regression and model selection. Our results showed that 12 species were distributed only in karst area, accounting for 5.48% of the total species pool, 104 species were found in non-karst area (47.49% of total species), and 103 species were found in both karst area and non-karst area (47.03% of total species). Based on the raw data of museum collections data and MaxEnt species distribution modeling, we found that amphibian species richness in the study area decreased at higher latitudes. Karst landforms and non-karst landforms differed in their distribution patterns of amphibian species richness (χ2 = 36.47, P < 0.0001), but the model was a poor fit to the data (McFadden’s Rho square = 0.0037). The most significant environmental predictors of species richness were mean annual rainfall (R2 = 0.232, P < 0.001) and precipitation of driest Month (R2 = 0.221, P < 0.001). The results based on model selection showed that underlying mechanisms related to landforms and different ecological hypotheses might simultaneously explain patterns of amphibian species richness in the study area. Future research should examine other biological factors such as interference, predation, and competition to understand the mechanisms controlling patterns of amphibian species richness.

Key words: biogeography, amphibia, species diversity, karst landforms, MaxEnt model

Table 1

Environmental variables represented by the four hypotheses"

Environmental variable
年均温度 Annual mean temperature (℃) BIO1
平均日温差 Mean diurnal range (℃) BIO2
等温性 Isothermality (℃) BIO3
季节性温度变化 Temperature seasonality BIO4
Max temperature of warmest month (℃)
Min temperature of coldest month (℃)
年均温差 Temperature annual range (℃) BIO7
Mean temperature of wettest quarter (℃)
Mean temperature of driest quarter (℃)
Mean temperature of warmest quarter (℃)
Mean temperature of coldest quarter (℃)
Mean annual sunshine (percent of daylength)
Mean annual frost-day frequency (days)
Mean annual potential evapotranspiration (mm/yr)
年均风速 Mean annual wind speed (m/s) WIND
Annual precipitation (mm/yr)
Precipitation of wettest month (mm/yr)
Precipitation of driest month (mm/yr)
季节性降雨量 Precipitation seasonality BIO15
Precipitation of wettest quarter (mm/yr)
Precipitation of driest quarter (mm/yr)
Precipitation of warmest quarter (mm/yr)
Precipitation of coldest quarter (mm/yr)
Mean annual relative humidity (%)
Productive energy
Normalized difference vegetation index
Mean annual actual evapotranspiration (mm/yr)
Mean annual solar radiation (W/m2)
海拔 Elevation ELE
Vegetation (number of vegetation
地貌类型 Landform LANDF

Table 2

Standardized regression coefficients of linear regression, t statistics, R2 and associated P-values for amphibian species richness against environmental predictors"

环境变量 Enviromental predictor AICc ΔAICc Wi t R2 P
年均降雨量 Annual precipitation (mm/yr) 41,710 0 1 48.465 0.232 < 0.001
最干月平均降雨量 Precipitation of driest month (mm/yr) 41,820 110 0 46.963 0.221 < 0.001
日照时数(白昼长的百分比) Mean annual sunshine (percent of daylength) 43,021 1,311 0 -27.903 0.091 0
年均实际蒸散量 Mean annual actual evapotranspiration (mm/yr) 43,071 1,361 0 26.895 0.085 < 0.001
年均相对湿度 Mean annual relative humidity (%) 43,479 1,769 0 17.021 0.036 0
年均风速 Mean annual wind speed 43,666 1,956 0 -9.861 0.012 < 0.001
年均潜在蒸散量 Mean annual potential evapotranspiration (mm/yr) 43,753 2,043 0 3.129 0.001 0.002
年均温度 Annual mean temperature (℃) 43,568 1,858 0 14.043 0.025 < 0.001
植被类型数 Vegetation (number of vegetation classes/quadrat) 43,758 2,048 0 -2.251 < 0.001 0.024
归一化植被指数 Normalized difference vegetation index 43,761 2,051 0 1.38 < 0.001 0.168
年均太阳辐射量 Mean annual solar radiation (W/m2) 43,761 2,051 0 -1.453 < 0.001 0.146

Fig. 1

Species distribution modeling of amphibians in Southwest China, based on (A) the raw data of museum collections data and (B) MaxEnt species distribution modeling. The area of each cell is 10 km × 10 km. Blank cells have no specimen based on the major collections."

Fig. 2

The relationships between amphibian richness and longitude (a) and latitude (b) based on the raw data of museum collections"

Fig. 3

The relationships between amphibian richness and longitude (a) and latitude (b) based on MaxEnt species distribution modeling"

Fig. 4

Moran’s I index correlogram for amphibian species richness and the residuals of multiple regression with environmental predictors"

Table 3

AICc value and adjusted R2 value for each hypothesis based on model selection approach"

假说 Hypothesis R2 AICc ΔAICc K Wi
Productivity energy
0.004 6,908 4,268 4 0
Habitat heterogeneity
0.005 6,898 4,258 3 0
能量假说 Ambient energy 0.173 5,462 2,822 5 0
Water-energy balance hypothesis
0.34 3,706 1,066 4 0
混合模型 Mixed model* 0.425 2,640 0 12 1
[1] Ahmed SE, McInerny G, O’Hara K, Harper R, Salido L, Emmott S, Joppa LN (2015) Scientists and software—surveying the species distribution modelling community. Diversity and Distributions, 21, 258-267.
[2] Ahn CH, Tateishi R (1994) Development of a global 30-minute grid potential evapotranspiration data set. Journal of the Japan Society of Photogrammetry, 33, 12-21.
[3] Anderson DR, Burnham KP, Thompson WL (2000) Null hypothesis testing: Problems, prevalence, and an alternative. Journal of Wildlife Management, 64, 912-923.
[4] Andrews P, O’Brien EM (2000) Climate, vegetation, and predictable gradients in mammal species richness in southern Africa. Journal of Zoology, 251, 205-231.
[5] Badgley C, Fox DL (2000) Ecological biogeography of North American mammals: Species density and ecological structure in relation to environmental gradients. Journal of Biogeography, 27, 1437-1467.
[6] Bohning-Gaese K (1997) Determinants of avian species richness at different spatial scales. Journal of Biogeography, 24, 49-60.
[7] Buckley LB, Jetz W (2007) Environmental and historical constraints on global patterns of amphibian richness. Proceedings of the Royal Society B: Biological Sciences, 274, 1167-1173.
[8] Burnham KP, Anderson DR (1998) Model Selection and Inference: A Practical in Formation Theoretic Approach. Springer-Verlag, New York.
[9] Costa GC, Nogueira C, Machado RB, Colli GR (2007) Squamate richness in the Brazilian Cerrado and its environmental-climatic associations. Diversity and Distributions, 6, 714-724.
[10] Currie DJ (1991) Energy and large-scale patterns of animal- and plant-species richness. The American Naturalist, 137, 27-49.
[11] Diniz-Filho JAF, Bini LM, Hawkins BA (2003) Spatial autocorrelation and red herrings in geographical ecology. Global Ecology & Biogeography, 12, 53-64.
[12] Fei L, Ye CY, Jiang JP (2012) Colored Atlas of Chinese Amphibians and Their Distributions. Sichuan Science and Technology Press, Chengdu. (in Chinese)
[费梁, 叶昌媛, 江建平 (2012) 中国两栖动物及其分布彩色图鉴. 四川科学技术出版社, 成都.]
[13] Ford DC, Williams PW (1989) Karst Geomorphology and Hydrology. Unwin Hyman, London.
[14] Gaston KJ (2000) Global patterns in biodiversity. Nature, 405, 220-227.
[15] Guo Q, Kelt DA, Sun Z, Liu H, Hu L, Ren H, Wen J (2013) Global variation in elevational diversity patterns. Scientific Reports, 3, 3007.
[16] Hanley JA, Mcneil BJ (1982) The meaning and use of the area under a Receiver Operating Characteristic (ROC) curve. Radiology, 143, 29-36.
[17] Hawkins BA, Diniz-Filho JA, Jaramillo CA, Soeller SA (2007) Climate, niche conservatism, and the global bird diversity gradient. The American Naturalist, 170, S16-S27.
[18] Hawkins BA, Field R, Cornell HV, Currie DJ, Guegan JF, Kaufman DM, Kerr JT, Mittelbach GG, Oberdorff T, O’Brien EM, Porter EE, Turner JRG (2003b) Energy, water, and broad-scale geographic patterns of species richness. Ecology, 84, 3105-3117.
[19] Hawkins BA, Porter EE, Diniz JAF (2003a) Productivity and history as predictors of the latitudinal diversity gradient of terrestrial birds. Ecology, 84, 1608-1623.
[20] Hou MF, López-pujol J, Qin HN, Wang LS, Liu Y (2010) Distribution pattern and conservation priorities for vascular plants in Southern China: Guangxi Province as a case study. Botanical Studies, 51, 377-386.
[21] Huang Y, Dai Q, Chen Y, Wan HF, Li JT, Wang YZ (2011) Lizard species richness patterns in China and its environmental associations. Biodiversity and Conservation, 7, 1399-1414.
[22] Jiménez-Valverde A, Peterson AT, Soberón J, Overton JM, Arago?n P, Lobo JM (2011) Use of niche models in invasive species risk assessments. Biological Invasions, 13, 2785-2797.
[23] Kerr JT, Packer L (1997) Habitat heterogeneity as a determinant of mammal species richness in high-energy regions. Nature, 385, 252-254.
[24] Legendre P (1993) Spatial autocorrelation: Trouble or new paradigm. Ecology, 74, 1659-1673.
[25] Li JJ, Zhou CH, Cheng WM (2009) 1 : 1000000 Geomorphological Atlas of the People’s Republic of China. Science Press, Beijing. (in Chinese)
[李吉均, 周成虎, 程维明 (2009) 中华人民共和国1 : 100万地貌图. 科学出版社, 北京.]
[26] Luo ZH, Tang SH, Jiang ZG, Chen J, Fang HX, Li CW (2016) Conservation of terrestrial vertebrates in a global hotspot of karst area in Southwestern China. Scientific Reports, 6, 25717.
[27] Mo YM, Wei ZY, Chen WC (2014) Color Atlas of Amphibians in Guangxi. Guangxi Science and Technology Press, Nanning. (in Chinese)
[莫运明, 韦振逸, 陈伟才 (2014) 广西两栖动物彩色图鉴. 广西科学技术出版社, 南宁.]
[28] New M, Lister D, Hulme M, Makin I (2002) A high-resolution data set of surface climate over global land areas. Climate Research, 21, 1-25.
[29] Phillips SJ, Anderson RP, Schapire RE (2006) Maximum entropy modeling of species geographic distributions. Ecological Modelling, 190, 231-259.
[30] Phillips SJ, Dudík MM (2008) Modeling of species distributions with Maxent: New extensions and a comprehensive evaluation. Ecography, 31, 161-175.
[31] Powney GD, Grenyer R, Orme CDL, Owens IPF, Meiri S (2010) Hot, dry and different Australian lizard richness is unlike that of mammals, amphibians, and birds. Global Ecology & Biogeography, 19, 386-396.
[32] Rahbek C, Graves GR (2001) Multiscale assessment of patterns of avian species richness. Proceedings of the National Academy of Sciences, USA, 98, 4534-4539.
[33] Rangel TF, Diniz-Filho JAF, Bini LM (2010) SAM: A comprehensive application for spatial analysis in macroecology. Ecography, 33, 46-50.
[34] Ricklefs RE (1987) Community diversity: Relative roles of local and regional processes. Science, 235, 167-171.
[35] Ruggiero A, Kitzberger T (2004) Environmental correlates of mammal species richness in South America: Effects of spatial structure, taxonomy and geographic range. Ecography, 27, 401-417.
[36] Schall JJ, Pianka ER (1978) Geographical trends in numbers of species. Science, 201, 679-686.
[37] Schoener TW (1976) The species-area relationship within archipelagoes: Models and evidence from island birds. Proceedings of XVI International Ornithological Congress, Australian Academy of Science, Canberra, Australia.
[38] Shen YH, Yang DD, Mo XY (2014) Fauna of Hunan: Amphibian. Hunan Science and Technology Press, Changsha. (in Chinese)
[沈猷慧, 杨道德, 莫小阳 (2014) 湖南动物志: 两栖纲. 湖南科学技术出版社, 长沙.]
[39] Shmida A, Wilson MV (1985) Biological determinants of species diversity. Journal of Biogeography, 12, 1-20.
[40] Stuart SN, Chanson JS, Cox NA, Young BE, Rodrigues AS, Fischman DL, Waller RW (2004) Status and trends of amphibian declines and extinctions worldwide. Science, 306, 1783-1786.
[41] Su ZM, Li XK (2003) The types of natural vegetation in karst region of Guangxi and its classified system. Guihaia, 23, 289-293. (in Chinese with English abstract)
[苏宗明, 李先琨 (2003) 广西岩溶植被类型及其分类系统. 广西植物, 23, 289-293.]
[42] Terribile LC, Olalla-Tárraga MA, Morales-Castilla I, Rueda M, Vidanes RM, Rodríguez MA, Diniz-Filho JA (2009) Global richness patterns of venomous snakes reveal contrasting influences of ecology and history in two different clades. Oecologia, 159, 617-626.
[43] Tognelli MF, Kelt DA (2004) Analysis of determinants of mammalian species richness in South America using spatial autoregressive models. Ecography, 27, 427-436
[44] Wang SJ, Lu HM, Zhou YC, Xie LP, Xiao DA (2007) Spatial variability of soil organic carbon and representative soil sampling method in Maolan karst virgin forest. Acta Pedologica Sinica, 44, 475-483. (in Chinese with English abstract)
[王世杰, 卢红梅, 周运超, 谢丽萍, 肖德安 (2007) 茂兰喀斯特原始森林土壤有机碳的空间变异性与代表性土样采集方法. 土壤学报, 44, 475-483.]
[45] Wright DH (1983) Species-energy theory: An extension of species-area theory. Oikos, 41, 496-506.
[46] Wu L, Dong Q, Xu RH (1987) Amphibians in Guizhou. Guizhou People’s Press, Guiyang. (in Chinese)
[伍律, 董谦, 须润华 (1987) 贵州两栖类志. 贵州人民出版社, 贵阳.]
[47] Xiong KN, Xiao SZ, Liu ZQ, Chen PD (2008) Comparative analysis on world natural heritage value of South China Karst. Engineering Science, 10, 17-28. (in Chinese with English abstract)
[熊康宁, 肖时珍, 刘子琦, 陈品冬 (2008) “中国南方喀斯特”的世界自然遗产价值对比分析. 中国工程科学, 10, 17-28.]
[48] Yang DT, Rao DQ (2008) Amphibians and Reptiles in Yunnan. Yunnan Science and Technology Press, Kunming. (in Chinese)
[杨大同, 饶定齐 (2008) 云南两栖爬行动物. 云南科技出版社, 昆明.]
[49] Yuan TX, Zhang HP, Ou ZY, Tan YB (2014) Effects of topography on the diversity and distribution pattern of ground plants in karst montane forests in Southwest Guangxi, China. Chinese Journal of Applied Ecology, 25, 2803-2810. (in Chinese with English abstract)
[袁铁象, 张合平, 欧芷阳, 谭一波 (2014) 地形对桂西南喀斯特山地森林地表植物多样性及分布格局的影响. 应用生态学报, 25, 2803-2810.]
[50] Zheng Z, Gong DJ, Sun CX, Li XJ, Li WJ (2014) Elevational pattern of amphibian and reptile diversity in Qinling Range and explanation. Biodiversity Science, 22, 596-607. (in Chinese with English abstract)
[郑智, 龚大洁, 孙呈祥, 李晓军, 李万江 (2014) 秦岭两栖、爬行动物物种多样性海拔分布格局及其解释. 生物多样性, 22, 596-607.]
[1] Hu Yifeng, Yu Wenhua, Yue Yang, Huang Zhenglanyi, Li Yuchun, Wu Yi. Species diversity and potential distribution of Chiroptera on Hainan Island, China [J]. Biodiv Sci, 2019, 27(4): 400-408.
[2] Yan Wenbo,Ji Shengnan,Shuai Lingying,Zhao Leigang,Zhu Dapeng,Zeng Zhigao. Spatial distribution patterns of mammal diversity in Yangxian County of Shaanxi Province on the southern slope of the Qinling Mountains [J]. Biodiv Sci, 2019, 27(2): 177-185.
[3] Chen Zuoyi, Xu Xiaojing, Zhu Suying, Zhai Mengyi, Li Yang. Species diversity and geographical distribution of the Chaetoceros lorenzianus complex along the coast of China [J]. Biodiv Sci, 2019, 27(2): 149-158.
[4] Hongfei Zhuang,Yinbo Zhang,Wei Wang,Yueheng Ren,Fangzheng Liu,Jinhong Du,Yue Zhou. Optimized hot spot analysis for probability of species distribution under different spatial scales based on MaxEnt model: Manglietia insignis case [J]. Biodiv Sci, 2018, 26(9): 931-940.
[5] Dexin Sun, Xiang Liu, Shurong Zhou. Dynamical changes of diversity and community assembly during recovery from a plant functional group removal experiment in the alpine meadow [J]. Biodiv Sci, 2018, 26(7): 655-666.
[6] Yu Zhang, Gang Feng. Distribution pattern and mechanism of insect species diversity in Inner Mongolia [J]. Biodiv Sci, 2018, 26(7): 701-706.
[7] Xiaqiu Tao,Shaopeng Cui,Zhigang Jiang,Hongjun Chu,Na Li,Daode Yang,Chunwang Li. Reptilian fauna and elevational patterns of the reptile species diversity in Altay Prefecture in Xinjiang, China [J]. Biodiv Sci, 2018, 26(6): 578-589.
[8] Cheng Tian,Junqing Li,Xuyu Yang,Lin Yu,Dan Yuan,Yunxi Li. Preliminary surveys of wild animals using infrared camera in Wanglang National Nature Reserve, Sichuan Province [J]. Biodiv Sci, 2018, 26(6): 620-626.
[9] Qian YANG, Wei WANG, Hui ZENG. Effects of nitrogen addition on the plant diversity and biomass of degraded grasslands of Nei Mongol, China [J]. Chin J Plan Ecolo, 2018, 42(4): 430-441.
[10] Han-Dong WEN, Lu-Xiang LIN, Jie YANG, Yue-Hua HU, Min CAO, Yu-Hong LIU, Zhi-Yun LU, You-Neng XIE. Species composition and community structure of a 20 hm2 plot of mid-mountain moist evergreen broad-leaved forest on the Mts. Ailaoshan, Yunnan Province, China [J]. Chin J Plan Ecolo, 2018, 42(4): 419-429.
[11] Xiaoyu Wu,Shikui Dong,Shiliang Liu,Quanru Liu,Yuhui Han,Xiaolei Zhang,Xukun Su,Haidi Zhao,Jing Feng. Identifying priority areas for grassland endangered plant species in the Sanjiangyuan Nature Reserve based on the MaxEnt model [J]. Biodiv Sci, 2018, 26(2): 138-148.
[12] Yiming Hu,Jianchao Liang,Kun Jin,Zhifeng Ding,Zhixin Zhou,Huijian Hu,Zhigang Jiang. The elevational patterns of mammalian richness in the Himalayas [J]. Biodiv Sci, 2018, 26(2): 191-201.
[13] LIU Hai-Yue, LI Xin-Mei, ZHANG Lin-Lin, WANG Jiao-Jiao, HE Xue-Li. Eco-geographical distribution of arbuscular mycorrhizal fungi associated with Hedysarum scoparium in the desert zone of northwestern China [J]. Chin J Plan Ecolo, 2018, 42(2): 252-260.
[14] Tianqi Shang, Nuonan Ye, Haiqing Gao, Hongdi Gao, Lita Yi. Community Structure Analysis of a Public Welfare Forest Based on Multivariate Regression Trees [J]. Chin Bull Bot, 2018, 53(2): 238-249.
[15] Zihong Chen, Xiaona Yang, Ningjing Sun, Ling Xu, Yuan Zheng, Yuming Yang. Species diversity and vertical distribution characteristics of Metarhizium in Gaoligong Mountains, southwestern China [J]. Biodiv Sci, 2018, 26(12): 1308-1317.
Full text



[1] Yili Guo,Dongxing Li,Bin Wang,Yunlin He,Wusheng Xiang,Yuliang Jiang,Xiankun Li. Composition and spatio-temporal dynamics of litter fall in a northern tropical karst seasonal rainforest in Nonggang, Guangxi, southern China[J]. Biodiv Sci, 2017, 25(3): 265 -274 .
[2] Li Jian-dong, Zheng Hui-ying. A Study of the Characteristics and the Eco-geographical Law of Stipa Steppe on the Song-Nen Plain[J]. Chin J Plan Ecolo, 1987, 11(3): 183 -192 .
[3] Ke Cao,Mide Rao,Jianping Yu,Xiaojuan Liu,Xiangcheng Mi,Jianhua Chen. The phylogenetic signal of functional traits and their effects on community structure in an evergreen broad-leaved forest[J]. Biodiv Sci, 2013, 21(5): 564 -571 .
[4] Zheng Yuanrun. The Applicability of Various Methods in Analysis of Picea mongolica Population Spatial Distribution Pattern[J]. Chin J Plan Ecolo, 1997, 21(5): 480 -484 .
[5] Zhu Hui-Zhong. Materials for the Oedogoniaceae of China (I)[J]. J Syst Evol, 1964, 9(3): 203 -212 .
[6] Jie-Lin GE, Gao-Ming XIONG, Jia-Xiang LI, Wen-Ting XU, Chang-Ming ZHAO, Zhi-Jun LU, Yue-Lin LI, Zong-Qiang XIE. Litter standing crop of shrubland ecosystems in southern China[J]. Chin J Plan Ecolo, 2017, 41(1): 5 -13 .
[7] Shi Zhi-xin. New Taxa of Chrysophyta from China[J]. J Syst Evol, 1997, 35(3): 269 -272 .
[8] Shi Chen, Zhonglai Luo, and Dianxiang Zhang. Pre- and post-zygotic reproductive isolation between co-occurring Mussaenda pubescens var. alba and M. shikokiana (Rubiaceae)[J]. J Integr Plant Biol, 2014, 56(4): 411 -419 .
[9] Li Jiandong. The 4th Congress of Northeast Grassland Was Held in Changchun City[J]. Chin J Plan Ecolo, 1983, 7(1): 88 .
[10] HUANG Zhong-Liang, KONG Guo-Hui, WEI Ping. Plant species diversity dynamics in Dinghu Mountain forests[J]. Biodiv Sci, 1998, 06(2): 116 -121 .