中國(guó)蛇類(lèi)物種豐富度地理格局及其與生態(tài)因子的關(guān)系

蔡 波, 黃 勇, 陳躍英, 胡軍華, 郭憲光, 王躍招,*

(1. 中國(guó)科學(xué)院成都生物研究所, 四川 成都 610041; 2. 中國(guó)科學(xué)院研究生院, 北京 100049; 3. 廣西藥用植物園, 廣西 南寧 530023)

中國(guó)蛇類(lèi)物種豐富度地理格局及其與生態(tài)因子的關(guān)系

蔡 波1,2, 黃 勇3, 陳躍英1, 胡軍華1, 郭憲光1, 王躍招1,*

(1. 中國(guó)科學(xué)院成都生物研究所, 四川 成都 610041; 2. 中國(guó)科學(xué)院研究生院, 北京 100049; 3. 廣西藥用植物園, 廣西 南寧 530023)

物種豐富度地理格局成因是生態(tài)學(xué)和生物地理學(xué)研究重要目標(biāo)之一。生態(tài)假說(shuō)在解釋物種豐富度地理格局的成因上受到廣泛關(guān)注。該文基于100 km×100 km空間分辨率研究中國(guó)蛇類(lèi)物種豐富度的地理分布格局, 并結(jié)合生態(tài)假說(shuō)探討影響蛇類(lèi)分布格局的因素。該研究采用主軸鄰距法獲得基于特征值的空域數(shù)據(jù), 并同生態(tài)因子進(jìn)行多元回歸分析, 結(jié)果表明: (1) 中國(guó)蛇類(lèi)物種豐富度在經(jīng)、緯度上呈現(xiàn)多峰分布格局, 物種豐富度最高的地區(qū)位于東洋界亞熱帶、熱帶, 豐富度較低的地區(qū)位于青藏高原、北方草原荒漠、黃淮平原、兩湖平原及鄱陽(yáng)湖平原等; (2) 多元回歸分析能解釋56.5%的蛇類(lèi)物種豐富度變化，分析得出蛇類(lèi)物種豐富度格局的主要影響因子是歸一化植被指數(shù)、最冷季降水量和年溫差。(3) 模型選擇結(jié)果顯示, 在多元回歸分析中，P<0.05的變量 (歸一化植被指數(shù)、最冷季降水量和年溫差) 組成的模型是解釋蛇類(lèi)物種豐富度格局的最優(yōu)模型。這說(shuō)明蛇類(lèi)物種豐富度格局是由不同生態(tài)因子共同作用的結(jié)果?；谥袊?guó)蛇類(lèi)物種豐富度地理格局成因研究的復(fù)雜性, 該文提出在進(jìn)一步研究中，需重視各假說(shuō)中影響因子的選擇和人類(lèi)活動(dòng)的影響, 并在不同空間尺度上對(duì)蛇類(lèi)物種豐富度地理格局進(jìn)行綜合分析。

宏生態(tài)學(xué); 生態(tài)假說(shuō); 生物地理學(xué); 空間分布格局; 能量假說(shuō); 生境異質(zhì)性假說(shuō)

物種豐富度大尺度地理格局及其成因是生態(tài)學(xué)和生物地理學(xué)研究的重要目標(biāo)之一(Gaston,2000)。目前已有多種假說(shuō)解釋物種豐富度地理格局的形成機(jī)制(Currie, 1991; Gaston, 2000; Hawkins et al, 2003a,b; Rahbek, 2005), 其中, 受到廣泛關(guān)注的生態(tài)假說(shuō)(ecological hypothesis)主要包括生產(chǎn)力假說(shuō)(productivity hypothesis)、環(huán)境能量假說(shuō)(ambient energy hypothesis)、寒冷忍耐假說(shuō)(freezing tolerance hypothesis)、水能動(dòng)態(tài)假說(shuō)(water-energy dynamic hypothesis)及生境異質(zhì)性假說(shuō)(habitat heterogeneity hypothesis)(Gaston, 2000; Terribile et al, 2009)。生產(chǎn)力假說(shuō)認(rèn)為能量的增加提高了一個(gè)地區(qū)的凈初級(jí)生產(chǎn)力, 導(dǎo)致流向食物網(wǎng)各營(yíng)養(yǎng)級(jí)的能量增加, 保證了更多個(gè)體的存活, 從而提高了物種多樣性(Currie et al, 2004; Hawkins et al, 2003a, b; Wang et al, 2009; Wright, 1983); 環(huán)境能量假說(shuō)認(rèn)為環(huán)境能量直接影響生物個(gè)體的生理調(diào)節(jié)機(jī)制, 從而影響種群增長(zhǎng)進(jìn)而影響物種豐富度( Currie, 1991; Huang et al, 2011; Turner, 2004; Turner et al, 1987; Wang et al,2009); 寒冷忍耐假說(shuō)認(rèn)為多數(shù)物種因受最冷季低溫忍受力影響而不能在更冷的地區(qū)生存(Currie et al,2004; Hawkins, 2001; Hawkins et al, 2003a; Wang et al, 2009); 水能動(dòng)態(tài)假說(shuō)基于水分對(duì)生物能量利用過(guò)程的影響, 認(rèn)為生物多樣性的大尺度分布格局由水能動(dòng)態(tài)決定( Hawkins et al, 2003a; O'Brien, 1993;Wang et al, 2009); 生境異質(zhì)性假說(shuō)認(rèn)為物理因子或生物因子的復(fù)雜性孕育較多生態(tài)位, 環(huán)境越復(fù)雜,空間異質(zhì)性程度越高, 小生境類(lèi)型越多樣, 所能支持的物種數(shù)也越多( Cox & Moore, 2010 ; Jiang et al,1997; Kerr & Packer, 1997; Pianka, 1966)。但越來(lái)越多的研究表明物種豐富度地理格局的形成機(jī)制很復(fù)雜, 不能僅由某一個(gè)假說(shuō)來(lái)解釋(Currie, 1991;Hawkins et al, 2003a; Rahbek, 2005)。

中國(guó)疆域廣闊, 跨越古北與東洋兩大動(dòng)物區(qū)系。第三紀(jì)晚期以來(lái)發(fā)生的地質(zhì)事件, 特別是印度板塊北移與青藏高原隆起, 使得區(qū)域內(nèi)地理環(huán)境復(fù)雜, 氣候多樣, 物種(包括蛇類(lèi))多樣性豐富(Zhang,1999)。截止2012年2月, 我國(guó)境內(nèi)記錄蛇類(lèi)～222種( Gong et al, 2010; Guo et al, 2008; Guo & Wang,2011; Guo et al, 2011; Huang et al, 2012; Jiang &Zhao, 2009; Luo et al, 2010; Malhotra et al, 2011; Rao& Zhao, 2005; Vogel & David, 2010; Vogel & Luo,2011; Yang et al, 2011a, b; Zhang et al, 2011a; Zhang et al, 2011b; Zhao, 2006; Zhao & Adler, 1993), 約占全球蛇類(lèi)物種數(shù)的 7% (http://www.reptile-database.org)。關(guān)于中國(guó)蛇類(lèi)物種多樣性, 早在三千年前的《山海經(jīng)》中就有記載, 其分布格局與生態(tài)因子關(guān)系則在20世紀(jì)80年代開(kāi)始受到關(guān)注(Shen,1985)。盡管在一些兩棲爬行動(dòng)物的分布格局研究中涉及蛇類(lèi), 但尚無(wú)針對(duì)整個(gè)中國(guó)區(qū)域蛇類(lèi)分布格局的相關(guān)研究, 這與該類(lèi)群所具有的獨(dú)特進(jìn)化地位很不相稱。

不同類(lèi)群物種具有不同的生理特性, 適應(yīng)不同的生態(tài)環(huán)境, 其分布格局及其成因也有所不同(Gaston, 2000; Lin et al, 2009), 特別在全球氣候變化背景下, 更需要了解某一類(lèi)群的物種多樣性分布格局與生態(tài)因子關(guān)系機(jī)制, 以利于生物多樣性保護(hù)。因此, 本文將探討中國(guó)蛇類(lèi)物種豐富度分布格局及其影響因素, 結(jié)合生態(tài)假說(shuō)分析其成因, 為蛇類(lèi)多樣性保護(hù)提供科學(xué)依據(jù)。

1 材料和方法

1.1 物種數(shù)據(jù)

以爬行動(dòng)物數(shù)據(jù)庫(kù)(http://www.reptiledatabase.org)和《中國(guó)蛇類(lèi)》(Zhao, 2006)的分類(lèi)系統(tǒng)為依據(jù), 本研究收集中國(guó)蛇類(lèi)物種信息, 包括采集地點(diǎn)經(jīng)緯度、生境與海拔等。數(shù)據(jù)主要來(lái)源于中國(guó)科學(xué)院兩棲爬行動(dòng)物標(biāo)本館[Museum of Herpetology,Chengdu Institute of Biology (CIB), Chinese Academy of Sciences], 以及《中國(guó)動(dòng)物志·有鱗目·蛇亞目》、《中國(guó)蛇類(lèi)》和在以下雜志公開(kāi)發(fā)表的文獻(xiàn), 如《Zootaxa》、《Bonn Zoological Bulletin》、《Asian Herpetological Research》、《兩棲爬行動(dòng)物學(xué)報(bào)》、《動(dòng)物學(xué)研究》、《四川動(dòng)物》、《動(dòng)物分類(lèi)學(xué)報(bào)》等。若資料中記錄有經(jīng)、緯度數(shù)據(jù), 則采用原資料的經(jīng)、緯度; 若無(wú), 則根據(jù)資料記載地點(diǎn)(至少到縣級(jí)行政區(qū)劃)使用Google Earth查找物種分布的經(jīng)、緯度。另外, 我們還邀請(qǐng)相關(guān)專(zhuān)家對(duì)原始數(shù)據(jù)審閱、更新和修正。數(shù)據(jù)整理截止2011年8月, 共獲得19 575條記錄(圖1), 最終構(gòu)建了中國(guó)219種蛇類(lèi)分布數(shù)據(jù)庫(kù)。

圖1 基于標(biāo)本館記錄和各種文獻(xiàn)收集的中國(guó)蛇類(lèi)分布點(diǎn)圖Fig. 1 Raw sample point of snakes based on museum collections and published literature

1.2 生態(tài)數(shù)據(jù)

本研究共采用31個(gè)生態(tài)因子作為自變量(表1),并基于生態(tài)假說(shuō)將其分為5組: 生產(chǎn)力假說(shuō)、環(huán)境能量假說(shuō)、寒冷忍耐假說(shuō)、水能動(dòng)態(tài)假說(shuō)和生境異質(zhì)性假說(shuō)(Currie, 1991; Hawkins et al, 2003a; Kerr &Packer, 1997; Terribile et al, 2009; Wang et al, 2009)。相關(guān)假說(shuō)和生態(tài)因子如下:

生產(chǎn)力假說(shuō)一般用實(shí)際蒸散量(mean annual actual evapotranspiration, AET) (Currie & Paquin,1987; Huang et al, 2011)和歸一化植被指數(shù)(normalized difference vegetation index, NDVI)(Evans et al, 2006; Huang et al, 2011)衡量。AET 的空間分辨率為30′, 數(shù)據(jù)來(lái)源于http://www.grid.unep.ch/data/grid/gnv183.html。NDVI數(shù)據(jù)來(lái)自 1998—2008年, 空間分辨率為 1 km2, 來(lái)源于 http://westdc.westgis.ac.cn。

環(huán)境能量假說(shuō)采用溫度及相關(guān)指標(biāo)衡量(Currie,1991; Huang et al, 2011; Turner et al, 1987; Wang et al, 2009)。本文采用年均溫(℃)(annual mean temperature,Tmean)、等溫性(℃) (isothermality, ISO)、溫度季節(jié)性變化(℃) (temperature seasonality, TS) 、風(fēng)速(m/s)(mean annual wind, WIND)衡量, 其中Tmean、ISO、TS數(shù)據(jù)來(lái)自1961—1990年, 空間分辨率為 30″(～1 km2), 來(lái)源于 http://www.worldclim.org。WIND來(lái)自1981—1996年, 空間分辨率為500 m2, 來(lái)源于 http://www.geodata.cn/Portal。

寒冷忍耐假說(shuō)常用最冷季均溫、最冷月均溫或年最低溫等極端低溫相關(guān)指標(biāo)來(lái)衡量(Hawkins,2001; Hawkins et al, 2003a; Wang et al, 2009), 又因物種對(duì)低溫的耐受能力受所處環(huán)境溫度及其變化影響, 因此, 同時(shí)也需要采用一般環(huán)境溫度、溫差指標(biāo)來(lái)衡量(Ricklefs & Schluter, 1993; Wiens &Donoghue, 2004)。本文采用年最高溫(℃)(mean annual highest temperature,Tmax)、年最低溫(℃)(mean annual lowest temperature,Tmin)、最冷月均溫(℃)(temperature of coldest month,TCMmean)、最熱月均溫(℃)(temperature of warmest month,TWMmean)、最潮季均溫(℃)(mean temperature of wettest quarter,TWetQmean)、最旱季均溫(℃)(mean temperature of driest quarter,TDQmean)、最冷季均溫(℃)(mean temperature of coldest quarter,TCQmean)、最熱季均溫(℃)(mean temperature of warmest quarter,TWarmQmean)、年溫差(℃)(temperature annual range,TAR)、日溫差(℃)[mean diurnal range-mean of monthly (max temp-min temp),TDR]衡量, 這些數(shù)據(jù)來(lái)自 1961—1990年, 空間分辨率為 30″(～1 km2),來(lái)源于http://www.worldclim.org。

表1 表征各個(gè)假說(shuō)的生態(tài)因子Tab. 1 Tested ecological factors associated with the five hypotheses

水能動(dòng)態(tài)假說(shuō)是基于水分對(duì)生物能量利用過(guò)程的影響, 認(rèn)為生物多樣性大尺度分布格局由水能動(dòng)態(tài)決定, 通常用蒸散量、溫度、降水量和濕度相關(guān)指標(biāo)來(lái)衡量(Hawkins et al, 2003a; Wang et al,2009)。本文采用年均溫(℃)(mean annual temperature,Tmean)、實(shí)際蒸散量(mm/y)(actual evapotranspiration,AET)、潛在蒸散量(mm/a)(potential evapotranspiration,PET)、年降水量(mm/a)(mean annual precipitation,AP)、最潮月降雨量(mm/a)(precipitation of wettest month,PWM)、最旱月降雨量(mm/a)(precipitation of driest month,PDM)、季節(jié)性降水量(mm/a)[precipitation seasonality (coefficient of variation),PS]、最潮季降水量(mm/a)(precipitation of wettest quarter,PWetQ)、最旱季降水量(mm/a) (precipitation of driest quarter,PDQ)、最冷季降水量(mm/a)(precipitation of coldest quarter,PCQ)、最熱季降水量(mm/a)(precipitation of warmest quarter,PWarmQ)、相對(duì)濕度(%)(relative humidity,REH)、潮濕天數(shù)(mm/a)(wet-days, Number days with>0.1 mm rain per month,WD)衡量。其中PET的空間分辨率 為 30′, 來(lái) 源 于 http://www.grid.unep.ch/data/grid/gnv183.html。AP、PWM、PDM、PS、PWetQ、PDQ、PCQ、PWarmQ數(shù)據(jù)來(lái)自1961—1990年, 空間分辨率為 30″(～1 km2), 來(lái)源于 http://www.worldclim.org。REH、WD數(shù)據(jù)來(lái)自1961—1990年,空間分辨率為 10′, 來(lái)源于 http://www.geodata.cn/Portal。

生境異質(zhì)性假說(shuō)常用生境類(lèi)型或者海拔指標(biāo)衡量(Currie, 1991; Huang et al, 2011; Kerr & Packer,1997)。本文采用每個(gè)柵格的土地覆蓋類(lèi)型數(shù) (land use and land cover change， LUCC)、平均海拔(mean elevation, ELEmean)、最高海拔(maximum elevation,ELEmax)、最低海拔(minimum elevation, ELEmin)。海拔數(shù)據(jù)來(lái)源http://eros.usgs.gov, 空間分辨率為30′。植被類(lèi)型數(shù)據(jù)來(lái)源國(guó)家自然科學(xué)基金委員會(huì)“中國(guó)西部環(huán)境與生態(tài)科學(xué)數(shù)據(jù)中心” http://westdc.westgis.ac.cn, 空間分辨率為1 km2。

1.3 統(tǒng)計(jì)分析

不同空間尺度或面積對(duì)物種分布格局的形成有重要影響(He & Legendre, 1996; Willis & Whittaker,2002)。為了消除空間尺度或面積對(duì)物種分布格局的影響, 本研究將中國(guó)國(guó)土劃分為100 km×100 km的柵格系統(tǒng)(Huang et al, 2011); 為了防止海岸線和邊界地區(qū)不完整柵格影響后續(xù)統(tǒng)計(jì)分析, 去除了不足完整柵格面積75%的柵格(Huang et al, 2011); 為了消除島嶼的影響(Connor & McCoy, 2001; Schall &Pianka, 1978), 去除所有島嶼柵格。本研究利用ArcGIS軟件將蛇類(lèi)物種分布數(shù)據(jù)生成每個(gè)物種的Arc/Info柵格圖, 有某一物種分布的柵格賦值為“1”,無(wú)某一物種分布的柵格賦值為“0”, 然后疊加所有物種的柵格圖以生成物種豐富度格局, 最后投影為Alberts Equal-Area Conic Projection (雙標(biāo)準(zhǔn)緯線:N25°和 N47°, 中央經(jīng)度 E105°)。同時(shí), 所有生態(tài)因子投影重采樣為100 km×100 km并提取數(shù)據(jù)。

由于生態(tài)因子間的強(qiáng)相關(guān)性會(huì)給后續(xù)多元回歸分析帶來(lái)偏差, 因此研究中盡可能選擇相關(guān)系數(shù)低的變量并且每個(gè)假說(shuō)至少有一個(gè)變量代表(Huang et al, 2011)。于是, 在SPSS軟件中對(duì)生態(tài)因子進(jìn)行兩兩之間相關(guān)分析后, 選擇相關(guān)系數(shù)<0.7的變量 (表2)。

表2 生態(tài)因子的相關(guān)矩陣 (r <0.7)Tab. 2 Correlation matrix of ecological factors (r <0.7)

統(tǒng)計(jì)分析在 SAM v4.0 (Spatial Analysis in Macroecology, http://www.ecoevol.ufg.br/sam)(Rangel et al, 2006, 2010)軟件中完成。為了消除回歸殘差空間自相關(guān), 本研究采用主軸鄰距法(PCNM,principal coordinates neighbor matrices)獲得基于特征值的空域數(shù)據(jù)(eigenvector-based filters)并與生態(tài)因子一起進(jìn)行多元回歸分析, 該方法能有效降低回歸殘差空間自相關(guān)(Diniz-Filho & Bini, 2005)。為降低空間自相關(guān)性, 本文選擇最小殘差的 Moran’sI指數(shù)的空域文件(filters)(Diniz-Filho & Bini, 2005),將其與生態(tài)因子組成多元回歸模型。Moran’sI指數(shù)取值范圍為?1～1, “?1”意味著強(qiáng)烈的負(fù)空間自相關(guān), “0”意味著無(wú)相關(guān)性, “1”意味存在強(qiáng)烈的正空間自相關(guān)。該多元回歸分析包括了空間因素,r2由三部分組成: 僅生態(tài)因子(predictors only)、僅空域(filters only)、兩者(生態(tài)因子和空域)(predictors+filters)。

為了探討5個(gè)假說(shuō)對(duì)中國(guó)蛇類(lèi)物種豐富度分布格局的解釋程度, 本研究對(duì)各假說(shuō)包含的所有變量采用主軸鄰距法分別進(jìn)行多元回歸分析, 同時(shí)另增加兩個(gè)混合模型: 大混合模型(所有變量), 小混合模型(多元回歸分析中P<0.05的變量)。采用樣本大小矯正過(guò)的赤池信息準(zhǔn)則(AICc)來(lái)衡量模型擬合度優(yōu)良性。最小AICc值的模型是最簡(jiǎn)約或優(yōu)先考慮

的模型, 是在模型構(gòu)建中偏差和精度之間最優(yōu)化的權(quán)衡(Burnham & Anderson, 1998)。當(dāng)所有模型的AICc值都很接近時(shí), 用候選模型和最優(yōu)模型之差(ΔAICc)衡量相關(guān)模型擬合度。ΔAICc越小, 候選模型越接近最優(yōu)模型。即, 當(dāng)ΔAICc≤2時(shí), 數(shù)據(jù)強(qiáng)力支持該模型; 當(dāng)4≤ΔAICc≤7時(shí), 數(shù)據(jù)一定程度支持該模型; 當(dāng)ΔAICc＞10時(shí), 數(shù)據(jù)不支持該模型。使用赤池權(quán)重(Akaike weights, wi)估計(jì)的模型平均數(shù)來(lái)解決模型選擇的不確定性(Burnham &Anderson, 1998)。

2 結(jié) 果

2.1 物種豐富度地理格局

由圖2可見(jiàn)中國(guó)蛇類(lèi)物種豐富度熱點(diǎn)地區(qū)在橫斷山脈南部、南嶺、武夷山脈、大婁山?巫山山脈、青藏高原東南緣、秦嶺?大巴山脈、藏南河谷(雅魯藏布江下游地區(qū)), 以及北方的長(zhǎng)白山脈、小興安嶺和西部的天山山脈西部等。中國(guó)蛇類(lèi)物種豐富度并不隨緯度升高而單調(diào)降低, 而是在 N22°～28°為最高峰, 在 N28°～37°隨緯度升高而顯著降低, 在～N37°有一個(gè)低谷(圖 3A); 在 E75°～90°變化不大,在 E90°～110°隨著經(jīng)度升高而上升, 在 E110°～128°隨著經(jīng)度升高而降低, 在 E128°之后隨著經(jīng)度升高而上升, 豐富度最高值位于～E110°處 (圖3B)。結(jié)合圖 2可知中國(guó)蛇類(lèi)物種豐富度最高的地區(qū)出現(xiàn)東洋界亞熱帶、熱帶地區(qū), 物種豐富度較低的地區(qū)在青藏高原、北方草原荒漠地區(qū)、黃淮平原、兩湖平原、鄱陽(yáng)湖平原等地。

圖2 中國(guó)蛇類(lèi)物種豐富度柵格圖 (100 km×100 km)Fig. 2 Grid map of snake species richness patterns in China at 100 km spatial resolution

圖3 中國(guó)蛇類(lèi)物種豐富度與緯度(A)、經(jīng)度(B)關(guān)系(100 km×100 km)Fig. 3 Relationship of snake richness (100 km×100 km) in China with (A) latitude and (B) longitude

2.2 物種豐富度地理格局與生態(tài)因子的關(guān)系

主軸鄰距法分析的截尾距離(truncation distance)為142.265 km。選擇最小殘差的Moran’sI指數(shù)的空域文件將其與生態(tài)因子組成多元回歸模型, 有效消除了空間自相關(guān)性(圖4)。

圖4 蛇類(lèi)物種豐富度和多元回歸殘差的Moran’s I 指數(shù)圖Fig. 4 Moran’s I correlogram for snake species richness and the residuals of multiple regression with ecological predictors

使用空間特征值模型的多元回歸分析能夠解釋蛇類(lèi)物種豐富度總體 56.5%的變化(r2=0.565;F=66.86;P=0), 僅生態(tài)因子解釋了 39.9%, 僅空域解釋了 8.6%, 兩者(生態(tài)因子和空域)共同解釋了8.1%。在所有生態(tài)因子中, 歸一化植被指數(shù)和最冷季降水量是影響中國(guó)蛇類(lèi)物種豐富度最大的因子,且它們均與中國(guó)蛇類(lèi)物種豐富度呈正相關(guān), 年溫差是影響中國(guó)蛇類(lèi)物種豐富度次大的因子, 與中國(guó)蛇類(lèi)物種豐富度呈負(fù)相關(guān)(表3)。

2.3 生態(tài)假說(shuō)對(duì)物種豐富度地理格局的解釋

基于赤池信息量準(zhǔn)則的模型選擇方法, 小混合模型(多元回歸分析中P<0.05的變量)的ΔAICc最低,赤池權(quán)重為1(表4), 是解釋中國(guó)蛇類(lèi)物種豐富度格局的最優(yōu)模型。其他模型(如: 大混合模型、生境異質(zhì)性假說(shuō)、生產(chǎn)力假說(shuō)、環(huán)境能量假說(shuō)、寒冷忍耐假說(shuō)、水能動(dòng)態(tài)假說(shuō)) ΔAICc均>10, 赤池權(quán)重也很低。

表 3 中國(guó)蛇類(lèi)物種豐富度與生態(tài)因子的多元回歸分析Tab. 3 Snake richness regressed against environmental variables in China

3 討 論

3.1 中國(guó)蛇類(lèi)物種豐富度地理格局

物種豐富度具有空間異質(zhì)性特征(Gaston,2000), 蛇類(lèi)也不例外。中國(guó)區(qū)域蛇類(lèi)物種豐富度最高的地區(qū)在東洋界亞熱帶、熱帶地區(qū), 如: 武夷山脈、橫斷山脈南部、南嶺、藏南河谷(雅魯藏布江下游地區(qū))、大婁山?巫山山脈、青藏高原東南緣和秦嶺?大巴山脈(圖 2)。這些地區(qū)均為南方山區(qū), 與中國(guó)最冷季降雨量、歸一化植被指數(shù)高的地區(qū)(圖 5)高度重疊。另外北方豐富度較高地區(qū)為長(zhǎng)白山脈、小興安嶺、天山山脈西部等地(圖2), 這些地區(qū)的共同特點(diǎn)也是山區(qū), 與北方地區(qū)歸一化植被指數(shù)、年均溫高的地區(qū)(圖5)重疊較多。青藏高原、北方草原荒漠地區(qū)、大興安嶺地區(qū)是蛇類(lèi)物種豐富度低的地區(qū), 極端嚴(yán)寒或干旱, 與最冷季降雨量、歸一化植被指數(shù)、年均溫低的地區(qū)(圖5)重疊程度高。最冷季降水量、歸一化植被指數(shù)與蛇類(lèi)物種豐富度呈正相關(guān)。最冷季降水量高能有效保證這些地區(qū)在蛇類(lèi)冬眠期間有足夠的水分支撐蛇類(lèi)捕食對(duì)象的生存, 并保證蛇類(lèi)冬眠所需環(huán)境濕度, 幫助幼蛇順利度過(guò)冬眠期; 歸一化植被指數(shù)高能促進(jìn)蛇類(lèi)食物鏈下端的動(dòng)物大量生長(zhǎng)繁殖, 增加蛇類(lèi)數(shù)量和種數(shù)。年溫差則與蛇類(lèi)物種豐富度呈負(fù)相關(guān)。雖然冬眠是蛇類(lèi)有效的體溫調(diào)節(jié)機(jī)制, 但作為冷血?jiǎng)游? 蛇類(lèi)的最適氣溫范圍是20～30 ℃(Zhao, 2006), 過(guò)大溫差會(huì)阻礙蛇類(lèi)生長(zhǎng)發(fā)育(Sun, 2001), 因此年溫差高的地區(qū)蛇類(lèi)物種豐富度低, 年溫差低的地區(qū)蛇類(lèi)物種豐富度高。

表4 模型選擇 (K表示模型變量數(shù), 不包括截距; wi表示權(quán)重)Tab. 4 Summary of the model selection approach based on AICc wi was the Akaike weight (K was the number of variables of the model, not including intercept)

圖5 中國(guó)部分生態(tài)因子分布圖Fig. 5 Maps of some ecological factors in China

采用主軸鄰距法的多元回歸分析結(jié)果與海拔、最冷季降水量分布圖、年溫差分布圖具有良好的對(duì)應(yīng)關(guān)系, 但局部區(qū)域除外, 如: 東洋界的黃淮平原、兩湖平原、鄱陽(yáng)湖平原等區(qū)域蛇類(lèi)物種豐富度很低。這些區(qū)域的蜥蜴類(lèi)、哺乳類(lèi)物種豐富度也低(Lin et al, 2009; Huang et al, 2011)。原因可能是這些區(qū)域歷史上受河患影響嚴(yán)重且農(nóng)業(yè)發(fā)達(dá)(Gong, 2002;Lin & Mo, 2011; Zhao et al, 2007)。中國(guó)爬行動(dòng)物瀕危的主要原因是食用和醫(yī)用造成過(guò)度利用，亂采濫伐導(dǎo)致生境破碎化，污水和農(nóng)藥污染(Li & Wllcove,2005)。因此這些地區(qū)蛇類(lèi)物種豐富度低，除了河患導(dǎo)致外, 還可能是由人類(lèi)影響所致。具體原因有待進(jìn)一步研究。

中國(guó)蛇類(lèi)物種豐富度分布格局與蜥蜴類(lèi)較為一致(Huang et al, 2011), 但也存在細(xì)微差異, 較明顯的地區(qū)在北方草原荒漠和天山山脈西部地區(qū)。中國(guó)蜥蜴類(lèi)物種豐富度格局與年均霜期呈負(fù)相關(guān), 與海拔、植被類(lèi)型數(shù)及潮濕天數(shù)頻率呈正相關(guān)(Huang et al, 2011), 而年均霜期、植被類(lèi)型數(shù)及潮濕天數(shù)頻率與本研究中所選的蛇類(lèi)各生態(tài)因子相關(guān)性很高,海拔相關(guān)性相對(duì)較低(表 2), 海拔變幅對(duì)蜥蜴影響顯著而對(duì)蛇類(lèi)不顯著, 與Shen (1985)的結(jié)論一致。這可能是蛇類(lèi)、蜥蜴類(lèi)在這些區(qū)域的物種豐富度存在微小差異的原因之一; 其次可能是由于兩類(lèi)群本底資料調(diào)查不一致所致, 因而有待進(jìn)一步研究。

3.2 中國(guó)蛇類(lèi)物種豐富度地理格局的成因及其復(fù)雜性

生態(tài)假說(shuō)中, 除去生境異質(zhì)性假說(shuō)之外的假說(shuō)均屬于能量假說(shuō)(energy hypothesis), 該假說(shuō)認(rèn)為氣候的能量變化在初期產(chǎn)生了全球豐富度梯度, 能量可利用性導(dǎo)致并維持物種豐富度梯度(Hawkins et al,2003a; von Humboldt, 1808)。Hawkins et al (2003a)研究得出能量組成的氣候模型能更好地解釋大尺度豐富度格局, Guisan & Hofer (2003)也認(rèn)為中等尺度能量組成的氣候模型比地形模型(地形為生境異質(zhì)性假說(shuō)的成分之一)更能解釋爬行動(dòng)物物種豐富度格局, 而生境異質(zhì)性假說(shuō)很難在大尺度上解釋爬行動(dòng)物空間分布格局(Kiester, 1971; Rodríguez et al,2005)。但也存在其他觀點(diǎn)。中國(guó)爬行動(dòng)物物種豐富度格局成因研究中, 有的既支持生境異質(zhì)性假說(shuō)又支持能量假說(shuō)(Huang et al, 2011; Qian et al 2007;Shen, 1985), 有的僅支持能量假說(shuō)(Fu et al, 2007;Qian et al, 2007; Wang, 2004)。

中國(guó)蛇類(lèi)物種豐富度格局的強(qiáng)影響因子為最冷季降水量和歸一化植被指數(shù)。最冷季降水量與蛇類(lèi)物種豐富度分布格局呈正相關(guān), 該研究結(jié)果支持水能動(dòng)態(tài)假說(shuō)。美洲的珊瑚蛇屬(Terribile &Diniz-Filho, 2009)和中國(guó)爬行動(dòng)物( Fu et al, 2007;Qian et al, 2007; Wang, 2004)相關(guān)研究也支持該假說(shuō)。該假說(shuō)主要基于水分對(duì)生物體能量利用過(guò)程的影響(Wang et al, 2009), 能有效解釋整個(gè)地球生物多樣性梯度(Hawkins et al, 2003a), 但對(duì)動(dòng)物的影響可能是通過(guò)對(duì)植物生產(chǎn)力和生物量的影響而間接實(shí)現(xiàn)(Hawkins et al, 2003a)。歸一化植被指數(shù)對(duì)蛇類(lèi)豐富度格局的影響則一定程度上證明了這種間接性。NDVI是生產(chǎn)力假說(shuō)的重要內(nèi)容, 本研究表明NDVI在一定程度支持生產(chǎn)力假說(shuō), 通過(guò)控制區(qū)域內(nèi)的蛇類(lèi)食物數(shù)量來(lái)控制蛇類(lèi)種群密度, 進(jìn)而影響其物種豐富度水平。而研究中代表生境異質(zhì)性假說(shuō)的土地覆蓋類(lèi)型數(shù)、平均海拔、最高海拔、最低海拔與蛇類(lèi)物種豐富度格局相關(guān)性很小, 因此本研究支持能量假說(shuō), 不支持生境異質(zhì)性假說(shuō)。

物種豐富度分布格局成因研究中雖然出現(xiàn)了多種假說(shuō)或者理論, 但仍然沒(méi)有單一假說(shuō)能充分解釋(Gaston, 2000; Cox & Moore, 2010)。本研究表明各種假說(shuō)并不相互排斥, 而代表各假說(shuō)的生態(tài)因子共同影響中國(guó)蛇類(lèi)物種豐富度格局。本研究使用空間特征值模型的多元回歸分析能夠解釋中國(guó)蛇類(lèi)物種豐富度總體 56.5%的變化, 尚有 43.5%不能解釋。其復(fù)雜性主要表現(xiàn)在以下三個(gè)方面。

首先, 各假說(shuō)中生態(tài)因子的復(fù)雜性。眾多模型中, 小混合模型最能解釋中國(guó)蛇類(lèi)物種豐富格局,因?yàn)樗嗽诙嘣貧w分析中所有P<0.05的變量。這說(shuō)明物種豐富度格局影響因子多樣, 而現(xiàn)有的假說(shuō)解釋度低。如果進(jìn)行赤池信息量準(zhǔn)則的模型選擇時(shí)去除小混合模型, 得出的解釋模型將依次是水能動(dòng)態(tài)假說(shuō)、寒冷忍耐假說(shuō)、環(huán)境能量假說(shuō)及生產(chǎn)力假說(shuō)。出現(xiàn)這一情況的原因就是各假說(shuō)中生態(tài)因子成分的影響。很多生態(tài)因子之間關(guān)系密切, 各假說(shuō)之間也緊密聯(lián)系(Allen et al, 2007; Wang et al,2009), 各假說(shuō)或模型受生態(tài)因子影響很大, 稍有不慎則結(jié)果迥異, 如本研究中的大、小混合模型。因此, 應(yīng)當(dāng)重視各假說(shuō)所包含影響因子的選擇。建議先參考前人的研究, 選擇相關(guān)假說(shuō)生態(tài)因子, 盡可能全面獲取生態(tài)因子數(shù)據(jù)及后各個(gè)生態(tài)因子的變量(如最低值、最高值、平均值、變異系數(shù)等), 在統(tǒng)計(jì)分析前根據(jù)需要去掉相關(guān)性高的因子(如本文去掉相關(guān)性>0.7的因子), 以免影響統(tǒng)計(jì)結(jié)果。

其次, 尺度復(fù)雜性。豐富度格局會(huì)隨尺度變化,局域尺度進(jìn)程也影響區(qū)域尺度豐富度的觀測(cè)(Gaston, 2000)。尺度具有變異性, 各生態(tài)因子在不同尺度表現(xiàn)出不同的特征(Fu et al, 2007; Gaston,2000; Lü & Fu, 2001), 即對(duì)物種多樣性分布格局的支持度, 也影響各假說(shuō)對(duì)分布格局的解釋(Hawkins et al, 2003a)。Hu et al (2005)分析了非環(huán)境類(lèi)因子對(duì)物種數(shù)目協(xié)同變化的影響發(fā)現(xiàn), 不同鳥(niǎo)獸間的協(xié)同變化不存在面積效應(yīng), 但其空間尺度間的相關(guān)性差異明顯。本研究與后續(xù)各尺度的研究(結(jié)果另文發(fā)表)表明在小尺度與較大尺度上物種豐富度格局影響因子不同但存在關(guān)聯(lián), 局部尺度影響因子也會(huì)對(duì)區(qū)域尺度產(chǎn)生影響(圖5)。因此將來(lái)需要研究不同尺度對(duì)物種豐富度格局的影響, 以及不同尺度上各假說(shuō)對(duì)物種豐富度格局成因的解釋。

最后, 人類(lèi)因素影響的復(fù)雜性。人類(lèi)因素影響屬于生物因子影響之一, 從生物干擾、食物競(jìng)爭(zhēng)、捕食等作用上直接或間接影響物種分布格局。Li &Wllcove (2005)研究得出資源過(guò)度利用、生境破碎化以及各種污染是我國(guó)兩棲爬行動(dòng)物瀕危的主要原因。另外, 人類(lèi)生產(chǎn)生活導(dǎo)致的氣候變化將影響動(dòng)物的行為和物候(Bull & Vogt, 1979; Gargaglioni &Milsom, 2007), 改變分布范圍, 繼而改變豐富度甚至導(dǎo)致物種滅絕(Pounds et al, 1999; Sun & Zhang,2000)。蛇類(lèi)具有重要的經(jīng)濟(jì)、科研價(jià)值, 盡管在2000年出版《國(guó)家保護(hù)的有益的或者有重要經(jīng)濟(jì)、科學(xué)研究?jī)r(jià)值的陸生野生動(dòng)物名錄》所列的中國(guó)395種爬行動(dòng)物中, 蛇類(lèi)就有207種, 是當(dāng)時(shí)我國(guó)有分布記錄的所有蛇類(lèi), 但我國(guó)蛇類(lèi)受威脅程度很高,幾乎所有蛇的整體和器官都能作為中藥、工藝、食品材料 (Zhao, 2006)。因此, 必須重視人類(lèi)因素對(duì)我國(guó)蛇類(lèi)豐富度格局的影響。

中國(guó)蛇類(lèi)物種豐富度地理格局成因研究有其復(fù)雜性, 因此, 在進(jìn)一步研究中需要重視各假說(shuō)中影響因子的選擇和人類(lèi)活動(dòng)的影響, 并在不同空間尺度上對(duì)蛇類(lèi)物種豐富度地理格局進(jìn)行綜合分析。

致謝：感謝中國(guó)科學(xué)院成都生物研究所趙爾宓院士提供中國(guó)早期蛇類(lèi)研究文獻(xiàn), 感謝宜賓學(xué)院郭鵬教授、中國(guó)科學(xué)院成都生物研究所丁利博士等對(duì)中國(guó)蛇類(lèi)數(shù)據(jù)的修正。

Allen AP, Gillooly JF, Brown JH. 2007. Recasting the Species-Energy Hypothesis: The Different Roles of Kinetic and Potential Energy in Regulating Biodiversity [M]// Storch D, Marquet PA, Brown JH.Scaling Biodiversity. Cambridge: Cambridge University Press,283-299.

Bull JJ, Vogt RC. 1979. Temperature-dependent sex determination in turtles[J].Science,206(4423): 1186-1188.

Burnham KP, Anderson DR. 1998. Model Selection and Inference: A Practical Information-Theoretic Approach[M]. New York:Springer-Verlag.

Connor EF, McCoy ED. 2001. Species-Area Relationships [M] // Levin SA.Encyclopedia of Biodiversity. New York: Academic Press, 397-411.

Cox CB, Moore PD. 2010. Biogeography: An Ecological and Evolutionary Approach. (8th ed)[M]. New Jersey: John Wiley & Sons.

Currie DJ. 1991. Energy and large-scale patterns of animal- and plant-species richness[J].Amer Nat,137(1): 27-49.

Currie DJ, Mittelbach GG, Cornell HV, Field R, Guégan JF, Hawkins BA,Kaufman DM, Kerr JT, Oberdorff T, O'Brien E, Turner JRG. 2004.Predictions and tests of climate-based hypotheses of broad-scale variation in taxonomic richness[J].Ecol Lett,7(12): 1121-1134.

Currie DJ, Paquin V. 1987. Large-scale biogeographical patterns of species richness of trees[J].Nature,329(6137): 326-327.

Diniz-Filho JAF, Bini LM. 2005. Modelling geographical patterns in species richness using eigenvector-based spatial filters[J].Glob Ecol Biogeogr,14(2): 177-185.

Evans KL, James NA, Gaston KJ. 2006. Abundance, species richness and energy availability in the North American avifauna[J].Glob Ecol Biogeogr,15(4): 372-385.

Fu CZ, Wang JX, Pu ZC, Zhang SL, Chen HL, Zhao B, Chen JK, Wu JH.2007. Elevational gradients of diversity for lizards and snakes in the Hengduan Mountains, China[J].Biod Conserv,16(3): 707-726.

Gargaglioni LH, Milsom WK. 2007. Control of breathing in anuran amphibians[J].Comp Biochem Physiol A,147(3): 665-684.

Gaston KJ. 2000. Global patterns in biodiversity[J].Nature,405(6783):220-227.

Gong SP, Wang HJ, Yang CT, Yu CS, Pan DS, Wang FM. 2010.Protobothrops cornutusdiscovered in Tianjingshan Forest, Guangdong Province, China[J].Chn J Zool,45(3): 170-173.[龔世平, 王海京, 楊昌騰, 余春勝, 潘崠生, 王付民. 2010. 廣東天井山林區(qū)發(fā)現(xiàn)角原矛頭蝮. 動(dòng)物學(xué)雜志, 45(3): 170-173.]

Gong SS. 2002. Historical variation and sustainable utilization of the Jianghan-Dongting Plain’s wetlands[J].Resour Environ Yangtze Basin,11(6): 569-574. [龔勝生. 2002. 江漢——洞庭湖平原濕地的歷史變遷與可持續(xù)利用. 長(zhǎng)江流域資源與環(huán)境, 11(6): 569-574.]

Guisan A, Hofer U. 2003. Predicting reptile distributions at the mesoscale:relation to climate and topography[J].J Bioger,30(8): 1233-1243.

Guo P, Liu SY, Feng JC, He M. 2008. The Descriptionof a new species ofThermophis(Serpentes: Colubridae)[J].Sichuan J Zool,27(3): 321. [郭鵬, 劉少英, 馮今朝, 何苗. 2008. 溫泉蛇屬Thermophis一新種——四川溫泉蛇. 四川動(dòng)物, 27(3): 321.]

Guo P, Wang YZ. 2011. A new genus and species of cryptic Asian green pitviper (Serpentes: Viperidae: Crotalinae) from southwest China[J].Zootaxa,2918: 1-14.

Guo YH, Wu YK, He SP, Shi HT, Zhao EM. 2011. Systematics and molecular phylogenetics of Asian snail-eating snakes[J].Zootaxa,3001:57-64.

Hawkins BA. 2001. Ecology's oldest pattern?[J].Endeavour,25(3):133-134.

Hawkins BA, Field R, Cornell HV, Currie DJ, Guégan JF, Kaufman DM,Kerr JT, Mittelbach GG, Oberdorff T, O'Brien EM, Porter EE, Turner JRG. 2003a. Energy, water, and broad-scale geographic patterns of species richness[J].Ecology,84(12): 3105-3117.

Hawkins BA, Porter EE, Diniz-Filho JAF. 2003b. Productivity and history as predictors of the latitudinal diversity gradient of terrestrial birds[J].Ecology,84(6): 1608-1623.

He FL, Legendre P. 1996. On species-area relations[J].Amer Nat,148(4):719-737.

Hu HJ, Jiang ZG, Wang ZW. 2005. Impacts of the non-environmental factors on covariance between avian and mammalian species richness[J].Zool Res,26(4): 358-364.[胡慧建, 蔣志剛, 王祖望. 2005. 非環(huán)境類(lèi)因子對(duì)物種數(shù)目協(xié)同變化的影響. 動(dòng)物學(xué)研究, 26(4): 358-364]

Huang S, Ding L, Burbrink FT, Yan J, Huang JT, Ling C, Chen X, Zhang YP.2012. A new species of the genusElaphe(Squamata: Colubridae) from Zoige County, Sichuan, China[J].Asia Herpetol Res,3(1): 38-45.

Huang Y, Dai D, Chen YY, Wan HF, Li JT, Wang YZ. 2011. Lizard species richness patterns in China and its environmental associations[J].Biod Conserv,20(7): 1399-1414.

Hutchinson GE. 1959. Homage to Santa Rosalia or why are there so many kinds of animals?[J].Amer Nat,93(870): 145-159.

Jiang F, Zhao EM. 2009.Gloydius lijianlii, a new species from the northern coastal islands along Shandong peninsula (Reptilia, Squamata,Viperidae)[J].Acta Zootax Sin,34(3): 642-646. [江帆, 趙爾宓. 2009.山東半島北部沿岸島嶼蝮屬一新種(爬行綱,有鱗目,蝰科). 動(dòng)物分類(lèi)學(xué)報(bào), 34(3): 642-646.]

Jiang ZG, Ma KP, Han XG. 1997. Conservation Biology[M]. Zhejiang:Zhejiang Science and Technology Press. [蔣志剛, 馬克平, 韓興國(guó).1997. 保護(hù)生物學(xué). 杭州: 浙江科學(xué)技術(shù)出版社.]

Kerr JT, Packer L. 1997. Habitat heterogeneity as a determinant of mammal species richness in high-energy regions[J].Nature,385(6613):252-254.

Kiester AR. 1971. Species density of North American amphibians and

reptiles[J].Syst Zool,20(2): 127-137.

Li YM, Wllcove DS. 2005. Threats to vertebrate species in China and the United States[J].BioScience,55(2): 147-153.

Lin X, Wang ZH, Tang ZY, Zhao SQ, Fang JY. 2009. Geographic patterns and environmental correlates of terrestrial mammal species richness in China[J].Biodiv Sci,17(6): 652-663. [林鑫, 王志恒, 唐志堯, 趙淑清,方精云. 2009. 中國(guó)陸棲哺乳動(dòng)物物種豐富度的地理格局及其與環(huán)境因子的關(guān)系. 生物多樣性, 17(6): 652-663.]

Lin ZH, Mo XG. 2011. Interactions between evolution of farming systems and environment for agricultural production in the Huang-Huai-Hai Plain during the historic period[J].Chn J Eco-Agricul,19(5):1072-1079. [林忠輝, 莫興國(guó). 2011. 歷史時(shí)期黃淮海平原農(nóng)作制度變遷與農(nóng)業(yè)生產(chǎn)環(huán)境演變. 中國(guó)生態(tài)農(nóng)業(yè)學(xué)報(bào), 19(5): 1072-1079.]

Luo J, Ryabov SA, Luo Y, Gao HY, Luo ZR, Hu XC. 2010. Classification and distribution of the GenusDinodon[J].Sichuan J Zool,29(4):579-582, 677. [羅鍵, Ryabov SA, 羅穎, 高紅英, 羅至仁, 胡先成.2010. 鏈蛇屬的分類(lèi)分布現(xiàn)狀. 四川動(dòng)物, 29(4): 579-582, 677.]

Lü YH, Fu BJ. 2001. Ecological scale and scaling[J].Acta Ecol Sin,21(12):2096-2105. [呂一河, 傅伯杰. 2001. 生態(tài)學(xué)中的尺度及尺度轉(zhuǎn)換方法. 生態(tài)學(xué)報(bào), 21(12): 2096-2105.]

Malhotra A, Dawson K, Guo P, Thorpe RS. 2011. Phylogenetic structure and species boundaries in the mountain pitviperOvophis monticola(Serpentes: Viperidae: Crotalinae) in Asia[J].Molec Phylo Evol,59(2):444-457.

O'Brien EM. 1993. Climatic gradients in woody plant species richness:towards an explanation based on an analysis of southern Africa's woody flora[J].J Bioger,20(2): 181-198.

Pianka ER. 1966. Latitudinal gradients in species diversity: a review of concepts[J].Amer Nat,100(910): 33-46.

Pounds JA, Fogden MPL, Campbell JH. 1999. Biological response to climate change on a tropical mountain[J].Nature,398(6728): 611-615.

Qian H, Wang XH, Wang SL, Li YL. 2007. Environmental determinants of amphibian and reptile species richness in China[J].Ecography,30(4):471-482.

Rahbek C. 2005. The role of spatial scale and the perception of large-scale species-richness patterns[J].Ecol Lett,8(2): 224-239.

Rangel TF, Diniz-Filho JAF, Bini LM. 2010. SAM: a comprehensive application for spatial analysis in macroecology[J].Ecography,33(1):46-50.

Rangel TFLVB, Diniz-Filho JAF, Bini LM. 2006. Towards an integrated computational tool for spatial analysis in macroecology and biogeography[J].Glo Ecol Biog,15(4): 321-327.

Rao DQ, Zhao EM. 2005. A new record from China-Protobothrops Kaulbacki(Reptilia, Serpentes, Viperidae)[J].Acta Zootax Sin,30(1):209-211. [饒定齊, 趙爾宓. 2005. 中國(guó)一新紀(jì)錄種——緬北原矛頭蝮 (爬行綱, 蛇亞目, 蝰科). 動(dòng)物分類(lèi)學(xué)報(bào), 30(1): 209-211.]

Ricklefs RE, Schluter D. 1993. Species Diversity in Ecological Communities: Historical and Geographical Perspectives[M]. Chicago:University of Chicago Press.

Rodríguez Má, Belmontes JA, Hawkins BA. 2005. Energy, water and large-scale patterns of reptile and amphibian species richness in Europe[J].Acta Oecol,28(1): 65-70.

Schall JJ, Pianka ER. 1978. Geographical trends in numbers of species[J].Science,201(4357): 679-686.

Shen GS. 1985. Fuzzy cluster analyses to the geographical distribution of Amphibians and Reptiles in China[J].Sichuan J Zool,4(2): 13-16. [沈國(guó)坤. 1985. 中國(guó)兩棲爬行動(dòng)物地理分布的模糊聚類(lèi)分析. 四川動(dòng)物, 4(2): 13-16.]

Sun QH, Zhang ZW. 2000. The impact of Climate warming on the distribution of Chinese Birds[J].Chn J Zool,35(6): 45-48. [孫全輝, 張正旺. 2000. 氣候變暖對(duì)我國(guó)鳥(niǎo)類(lèi)分布的影響. 動(dòng)物學(xué)雜志, 35(6):45-48.]

Sun RY. 2001. Principles of Animal Ecology. (4th ed)[M]. Beijing: Beijing Normal University Publishing Group. [孫儒泳. 2001. 動(dòng)物生態(tài)學(xué)原理. (第四版)[M]. 北京: 北京師范大學(xué)出版社.]

Terribile LC, Diniz-Filho JAF. 2009. Spatial patterns of species richness in New World coral snakes and the metabolic theory of ecology[J].Acta Oecol,35(2): 163-173.

Terribile LC, Olalla-Tárraga Má, Morales-Castilla I, Rueda M, Vidanes RM,Rodríguez Má, Diniz-Filho JAF. 2009. Global richness patterns of venomous snakes reveal contrasting influences of ecology and history in two different clades[J].Oecologia,159(3): 617-626.

Turner JRG. 2004. Explaining the global biodiversity gradient: energy, area,history and natural selection[J].Basic Appl Ecol,5(5): 435-448.

Turner JRG, Gatehouse CM, Corey CA. 1987. Does solar energy control organic diversity? Butterflies, moths and the British climate[J].Oikos,48(2): 195-205.

Vogel G, David P. 2010. A new species of the genusLycodon(Boie, 1826)from Yunnan Province, China (Serpentes: Colubridae)[J].Bonn Zool Bull,57(2): 289-296.

Vogel G, Luo J. 2011. A new species of the genusLycodon(Boie, 1826)from the southwestern mountains of China (Squamata: Colubridae)[J].Zootaxa,2807: 29-40.

von Humboldt A. 1808. Ansichten der Natur: Mit Wissenschaftlichen Erl?uterungen[M]. Tübingen: J. G. Cotta’schen Buchhandlung.

Wang CH. 2004. Reseach on the Biodiversity Distribution Patterns in Mainland of China[D]. Ph.D. dissertation, Shanxi University, Taiyuan.[王翠紅. 2004. 中國(guó)陸地生物多樣性分布格局的研究. 山西大學(xué)博十論文]

Wang ZH, Tang ZY, Fang JY. 2009. The species-energy hypothesis as a mechanism for species richness patterns[J].Biodiv Sci,17(6): 613-624.[王志恒, 唐志堯, 方精云. 2009. 物種多樣性地理格局的能量假說(shuō).生物多樣性, 17(6): 613-624.]

Wiens JJ, Donoghue MJ. 2004. Historical biogeography, ecology and species richness[J].Trends Ecol Evol,19(12): 639-644.

Willis KJ, Whittaker RJ. 2002. Species diversity-scale matters[J].Science,295(5558): 1245-1248.

Wright DH. 1983. Species-energy theory: an extension of species-area theory[J].Oikos,41(3): 496-506.

Yang JH, Orlov N, Wang YY. 2011a. A new species of pitviper of the genusProtobothropsfrom China (Squamata: Viperidae)[J].Zootaxa,2936:59-68.

Yang JH, Wang YY, Zhang B, Lau MW, Chou WH. 2011b. Revision of the diagnostic characters ofOpisthotropis maculosaStuart and Chuaynkern, 2007 with notes on its distribution and variation, and a key to the genusOpisthotropis(Squamata: Natricidae)[J].Zootaxa,2785: 61-68.

Zhang J, Jiang K, Vogel G, Rao DQ. 2011a. A new species of the genusLycodon(Squamata, Colubridae) from Sichuan Province, China[J].Zootaxa,2982: 59-68.

Zhang L, Jiang K, Hu P, Yu BC, Peng YB, Tang XP, Hu HJ. 2011b.Lycodon futsingensis: a new snake record in Guangdong priovince, China[J].Chn J Zool,46(1): 128-130. [張亮, 蔣珂, 胡平, 遇寶成, 彭永波, 唐小平, 胡慧建. 2011. 廣東省蛇類(lèi)新紀(jì)錄——福清白環(huán)蛇. 動(dòng)物學(xué)雜志, 46(1): 128-130.]

Zhang RZ. 1999. Zoogeography of China[M]. Beijing: Science Press. [張榮祖. 1999. 中國(guó)動(dòng)物地理. 北京: 科學(xué)出版社.]

Zhao EM. 2006. Snake of China[M]. Hefei: Anhui Science and Technology Publishing House. [趙爾宓. 2006. 中國(guó)蛇類(lèi). 合肥: 安徽科學(xué)技術(shù)出版社.]

Zhao EM, Adler K. 1993. Herpetology of China[M]. Oxford, Ohio: Society for the Study of Amphibians and Reptiles in cooperation with Chinese Society for the Study of Amphibians and Reptiles.

Zhao QG, Huang GQ, Qian HY. 2007. Ecological environment and Sustainable development of poyang lake[J].Acta Pedol Sin,44(2):318-326.[趙其國(guó), 黃國(guó)勤, 錢(qián)海燕. 2007. 鄱陽(yáng)湖生態(tài)環(huán)境與可持續(xù)發(fā)展. 土壤學(xué)報(bào), 44(2): 318-326.]

Geographic patterns and ecological factors correlates of snake species richness in China

CAI Bo1,2, HUANG Yong3, CHEN Yue-Ying1, HU Jun-Hua1, GUO Xian-Guang1,WANG Yue-Zhao1,*

(1.Chengdu Institute of Biology,the Chinese Academy of Science,Chengdu610041,China; 2.Graduate School of Chinese Academy of Science,Beijing100049,China; 3.Guangxi Botanical Garden of Medicinal Plants,Nanning530023,China)

Understanding large-scale geographic patterns of species richness as well its underlying mechanisms are among the most significant objectives of macroecology and biogeography. The ecological hypothesis is one of the most accepted explanations of this mechanism. Here, we studied the geographic patterns of snakes and investigated the relationships between species richness and ecological factors in China at a spatial resolution of 100 km×100 km. We obtained the eigenvector-based spatial filters by Principal Coordinates Neighbor Matrices, and then analyzed ecological factors by multiple regression analysis. The results indicated several things: (1) species richness of snakes showed multi-peak patterns along both the latitudinal and longitudinal gradient. The areas of highest richness of snake are tropics and subtropical areas of Oriental realm in China while the areas of lowest richness are Qinghai-Tibet Plateau, the grasslands and deserts in northern China, Yangtze-Huai Plain, Two-lake Plain, and the Poyang-lake Plain; (2) results of multiple regression analysis explained a total of 56.5% variance in snake richness. Among ecological factors used to explore the species richness patterns, we found the best factors were the normalized difference vegetation index,precipitation in the coldest quarter and temperature annual range ; (3) our results indicated that the model based on the significant variables that (P<0.05) uses a combination of normalized difference vegetation index, precipitation of coldest quarterand temperature annual range is the most parsimonious model for explaining the mechanism of snake richness in China. This finding demonstrates that different ecological factors work together to affect the geographic distribution ofsnakes in China. Studying the mechanisms that underlie these geographic patterns are complex, so we must carefully consider the choice of impact-factors and the influence of human activities.

Macroecology; Ecological hypothesis; Biogeography; Spatial distribution patterns; Energy hypothesis;Habitat heterogeneity hypothesis

Q14

A

0254-5853-(2012)04-0343-11

10.3724/SP.J.1141.2012.04343

2012-03-20；接受日期：2012-05-05

國(guó)家自然科學(xué)基金面上項(xiàng)目(NSFC31071892); 中國(guó)科學(xué)院知識(shí)創(chuàng)新工程方向項(xiàng)目(KSCX2-EW-J-22); 科技部基礎(chǔ)專(zhuān)項(xiàng)(2006FY110500-3)

?通信作者(Corresponding author), E-mail: arcib@cib.ac.cn

蔡波(1986-), 男, 碩士研究生, 研究方向?yàn)樯鷳B(tài)生物地理學(xué)、動(dòng)物生態(tài)學(xué)、保護(hù)生物學(xué)。E-mail: caibo.cd@gmail.com