西藏高原冬蟲夏草資源適宜性區(qū)劃分析

周刊社,張建春,黃曉清,*,袁 雷,羅 珍,施 磊,魏克文

1 西藏自治區(qū)氣候中心,拉薩 850000 2 昌都市氣象局,昌都 854000 3 西藏自治區(qū)氣象信息與網(wǎng)絡中心,拉薩 850000 4 米林縣氣象局,林芝 860000

冬蟲夏草(Ophiocordycepssinensis(Berk.)G.H.Sung et al.)是蝙蝠蛾幼蟲被蟲草菌感染,死后尸體、組織與菌絲結成堅硬的假菌核[1]。從世界范圍看,冬蟲夏草主要分布在中國、尼泊爾、不丹和印度[2],中國主要分布在西藏的那曲地區(qū)、昌都市、林芝市,青海省的玉樹州、果洛州,四川省的阿壩州、甘孜州,甘肅省的甘南州,云南省的迪慶州等地[3- 4]。冬蟲夏草常見于海拔3500m以上的高寒山區(qū)和高山雪線附近,由于青藏高原生態(tài)地理環(huán)境的復雜性和多樣性,冬蟲夏草居群在區(qū)域尺度上呈不連續(xù)的斑塊狀分布[5]。不同地區(qū)由于土壤、植被、氣候狀況等的不同冬蟲夏草分布海拔高度有一定的差異[6- 11],總體上青藏高原主體分布海拔要高于周邊地區(qū)。袁峰模擬全球冬蟲夏草潛在適生分布區(qū)認為海拔2600—5300m是其主要潛在適生分布區(qū),海拔4200—4800m是集中分布區(qū)[12]。

生物一般均生活在較為適宜的生態(tài)地理環(huán)境中,為了某一生存目的,在可達的生境之間尋找一項最適宜生境,任何一種生物都是通過長期不斷的進化來適應棲息地的環(huán)境。在長期的進化過程中,冬蟲夏草寄主昆蟲與其生存的高寒草甸生態(tài)系統(tǒng)相互作用,長期適應,環(huán)境和氣候的變化會影響冬蟲夏草的產(chǎn)生和資源的可持續(xù)性。生態(tài)調(diào)查發(fā)現(xiàn)西藏高原冬蟲夏草一般生長在高山草甸、亞高山草甸和高山灌叢草甸之中,其分布與海拔、溫度、濕度、光照、土壤、植被等關系密切[13],降雨量和溫度是影響其產(chǎn)量的主要氣象因素[14]。冬蟲夏草生長地植被多為矮生墊狀灌叢,草本植物種類豐富,少有高大喬木,植被蓋度不高[15],朱印酒等[16]也對中國冬蟲夏草主要產(chǎn)區(qū)蝠蛾屬昆蟲生長發(fā)育的最適宜植被類型、土壤類型進行了調(diào)查。研究發(fā)現(xiàn)冬蟲夏草在生境平均氣溫達到2.6℃時開始生長,最適宜環(huán)境的氣溫為7.0—12.0℃[17],空氣相對濕度為80%—98%,土壤相對濕度為40%—60%[18]；云南西北部冬蟲夏草蝠蛾種類和種群主要分布在海拔3800m以上的高寒草甸區(qū)域,蝠蛾昆蟲種群的分布明顯受到海拔和草甸植物群體分布狀況的制約[19]；近幾十年青海玉樹、四川西部和云南西北部的冬蟲夏草分布海拔范圍均有所收窄[7,18,20- 21],據(jù)中國科學院的專家在冬蟲夏草研究基地的多年觀察發(fā)現(xiàn)部分破壞嚴重地區(qū)資源量不足30年前的1%—3%[22]。

因此,冬蟲夏草的產(chǎn)生與海拔高度、植被類型、土壤類型、氣候條件密切相關,其分布也與蝠蛾寄主昆蟲分布規(guī)律密切相關。許多學者雖然對不同地區(qū)冬蟲夏草生境做了調(diào)查,也對植被進行了分類,但大部分僅局限于單點觀測,缺乏全面性、系統(tǒng)性的研究。目前,冬蟲夏草資源適宜性區(qū)劃也只做了嘗試性研究[23- 24],沒有考慮當?shù)鼐唧w的氣候因素對冬蟲夏草產(chǎn)生的限制性影響。隨著氣候變暖,西藏高原冬蟲夏草資源數(shù)量及分布范圍也受到較大影響,如何利用現(xiàn)有資料做好冬蟲夏草資源適宜性區(qū)劃顯得尤為重要。本研究根據(jù)冬蟲夏草在西藏高原分布的地域性、海拔高度、土壤類型、植被類型及氣候因子等,結合冬蟲夏草資源調(diào)查基礎資料,利用地理信息系統(tǒng)(GIS)中的空間分析功能,研究西藏高原冬蟲夏草適宜性空間分布,以便為冬蟲夏草資源調(diào)查和保護提供方法,為當?shù)卣_展冬蟲夏草采挖等生態(tài)環(huán)境保護管理提供決策依據(jù)。

1 材料與方法

1.1 研究區(qū)概況

西藏高原是青藏高原的主體,是一個集高原、山地、河谷、盆地、湖盆、冰川等多種地貌兼而有之的地區(qū),該區(qū)域大致可分為南部喜馬拉雅山區(qū),喜馬拉雅北麓湖盆區(qū),雅魯藏布江中游谷地、藏北高原湖盆區(qū)和藏東高山峽谷區(qū)[25]。植被水平地帶性從東南到西北,垂直地帶性從低海拔到高海拔依次分布半常綠雨林、常綠闊葉林、針闊混交林、亞高山暗針葉林、高山灌叢、高山草甸、高山草原、荒漠草原、草原化高山墊狀植被、高山亞冰雪稀疏植被等[26]。土壤從東南向西北水平分為亞熱帶山地森林土、高山氈土、高山草原土、高山漠土等[27]。冬蟲夏草廣泛分布于高原東部、中部地區(qū),地理位置為84°6′—99°6′E,26°40′—32°40′N之間(圖1)。這些縣主要為畜牧業(yè)縣,冬蟲夏草一般生長在這些縣的高山草甸或灌叢中。

圖1 西藏高原冬蟲夏草產(chǎn)區(qū)海拔分布Fig.1 Altitude distribution of Ophiocordyceps sinensis in Tibetan Plateau

1.2 數(shù)據(jù)來源

本研究的主要數(shù)據(jù)包括：高程數(shù)據(jù)空間分辨率為30m的DEM數(shù)據(jù),行政邊界數(shù)據(jù)采用國家基礎地理信息中心發(fā)布的1∶400萬西藏行政區(qū)劃數(shù)據(jù),植被數(shù)據(jù)為地球系統(tǒng)科學數(shù)據(jù)共享平臺發(fā)布的1∶100萬西藏植被圖,土壤數(shù)據(jù)為西藏1∶100萬數(shù)字化土壤圖；氣象數(shù)據(jù)為西藏38個氣象站1981—2010年年平均降水量和年平均氣溫,由西藏自治區(qū)氣象信息網(wǎng)絡中心提供；1951—2000年的年平均降水量、年平均氣溫柵格數(shù)據(jù)由月氣候柵格數(shù)據(jù)集計算獲得,由http://www.worldclim.org/version1提供；冬蟲夏草產(chǎn)地所在縣及鄉(xiāng)鎮(zhèn)資料由西藏自治區(qū)農(nóng)牧廳草原監(jiān)理站提供和實地調(diào)研獲得。利用ArcGIS9.3軟件實現(xiàn)數(shù)據(jù)投影轉(zhuǎn)換、圖層疊加運算、區(qū)劃等級劃分等空間分析功能。

1.3 研究方法

1.3.1 空間分析方法

通過調(diào)查發(fā)現(xiàn)西藏高原冬蟲夏草主要分布在39個縣區(qū),對西藏冬蟲夏草在主產(chǎn)區(qū)(那曲地區(qū)比如縣、索縣、嘉黎縣和昌都市的丁青縣、類烏齊縣、昌都縣)、次產(chǎn)區(qū)(林芝市高海拔縣區(qū))和非產(chǎn)區(qū)(那曲地區(qū)西部縣區(qū)和阿里地區(qū)各縣)的生境進行調(diào)查,利用手持GPS測量冬蟲夏草采挖地的經(jīng)緯度及海拔高度。在對西藏冬蟲夏草資源分布范圍調(diào)查的基礎上,結合38個氣象站30年氣象資料進行研究,初步確定了冬蟲夏草不同產(chǎn)區(qū)的年平均降水量和年平均氣溫,同時也對不同產(chǎn)區(qū)的海拔高度、植被類型和土壤類型進行了初步調(diào)查,確立了影響冬蟲夏草產(chǎn)生的幾個主要限制性因素。首先確立以海拔高度、植被、土壤、年平均降水量、年平均氣溫作為冬蟲夏草資源適宜性區(qū)劃主要依據(jù),并制定區(qū)劃指標(表1)；再利用ArcGIS 9.3軟件的空間分析模塊分析每個指標所確定的冬蟲夏草適宜性空間分布,通過重采樣技術將海拔高度、植被類型、土壤類型、年平均降水量和年平均氣溫柵格數(shù)據(jù)的分辨率均統(tǒng)一為90m,5個要素圖層疊加運算,根據(jù)疊加后值的大小重新分類,結果形成適宜區(qū)、次適宜區(qū)和不適宜區(qū)的冬蟲夏草綜合適宜性分布圖；最后將西藏冬蟲夏草所在區(qū)縣行政邊界疊加到適宜性分布圖上,形成西藏冬蟲夏草適宜性空間分布圖。

表1 冬蟲夏草資源適宜性空間劃分依據(jù)

1.3.2 實地驗證方法

調(diào)查全區(qū)39個縣出產(chǎn)冬蟲夏草的鄉(xiāng)鎮(zhèn),獲得156個鄉(xiāng)鎮(zhèn)政府所在地的經(jīng)緯度,同時利用ArcGIS 9.3軟件在DEM數(shù)據(jù)中提取156個鄉(xiāng)鎮(zhèn)政府所在地的海拔高度。利用ArcGIS 9.3軟件將產(chǎn)區(qū)鄉(xiāng)鎮(zhèn)駐地點數(shù)據(jù)疊加在區(qū)劃圖上進行合理性判斷,再結合單位草原面積產(chǎn)量對區(qū)劃圖的代表性和科學性進行驗證。

2 結果與分析

2.1 冬蟲夏草適宜性分布單因素分析

2.1.1 冬蟲夏草海拔適宜性分析

西藏高原不同地方冬蟲夏草生長海拔與品質(zhì)差異較大,以西藏那曲地區(qū)東部和昌都市西部出產(chǎn)的冬蟲夏草為優(yōu),這些產(chǎn)地的海拔高度較其他產(chǎn)地要高。楊大榮等的研究表明,蟲草蝠蛾的分布具有垂直分布和區(qū)域水平分布的特點,蝠蛾屬每一個種群都有特定的地理位置和分布格局。雖然蝠蛾寄主昆蟲具有物種多樣性,但產(chǎn)區(qū)均具有氣溫低,晝夜溫差大,無絕對無霜期,凍土時間長,日照充足等氣候特征。蝠蛾寄主昆蟲在長期進化過程中形成了對低溫、大風、缺氧、氣壓低等高原氣候的長期適應性,而且主要生長在高山草甸和亞高山草甸中,生態(tài)位相對固定,所以根據(jù)海拔限制來研究冬蟲夏草的垂直分布是可行的。調(diào)查發(fā)現(xiàn)冬蟲夏草在那曲地區(qū)東部海拔4100—5000m的垂直地帶均有分布,海拔4300—4800m是冬蟲夏草分布最適區(qū)域,4900m以上冬蟲夏草數(shù)量較少,5000m以上、4100m以下未發(fā)現(xiàn)分布[14]。調(diào)查發(fā)現(xiàn)那曲地區(qū)冬蟲夏草產(chǎn)區(qū)海拔高度聶榮縣為4500—5100m、那曲縣為4500—4900m、嘉黎縣和索縣為4000—4800m,海拔具有北部縣高于南部縣的特征。金云翔等[23]在西藏高原那曲地區(qū)研究認為海拔4300—4700m為冬蟲夏草生長的適宜區(qū)；4701—5000m和3701—4300m為冬蟲夏草生長的次適宜區(qū)；5000m以上和3700m以下為冬蟲夏草生長的不適宜區(qū)。由于那曲地區(qū)東部為西藏高原冬蟲夏草的主產(chǎn)區(qū),該地冬蟲夏草產(chǎn)量高,品質(zhì)優(yōu),具有代表性。根據(jù)此標準,結合西藏高程數(shù)據(jù)將冬蟲夏草資源的海拔適宜性空間分布按海拔高度劃分為3個區(qū)域(圖2)。

圖2 西藏高原冬蟲夏草海拔適宜性分布Fig.2 Suitability distribution of Ophiocordyceps sinensis based on the altitude in Tibetan Plateau

2.1.2 冬蟲夏草植被類型適宜性分析

調(diào)查發(fā)現(xiàn)冬蟲夏草產(chǎn)地植被類型多為高山草甸、亞高山草甸和高山灌叢草甸,主要植物有高山嵩草(Kobresiapygmaea)、藏苔草(Carextristachya)、披堿草(Elymusnatans)、珠芽蓼(Polygonumviviparum)、圓穗蓼(Polygonummacrophyllum)、風毛菊(Saussureasalsa)、火絨草(Leontopodiumleontopodioides)、龍膽(GentianascabraBunge)、黃芪(AstragalusSP.)等,也發(fā)現(xiàn)伴生有金露梅(Potentillafruticosa)、土大黃(RumexmadaioMakino)、小大黃(RheumpumilumMaxim.)、輪葉黃精(Polygonatumverticillatum)和垂葉黃精(Polygonatumcurvistylum)等植物。調(diào)査研究認為冬蟲夏草適宜生長的植被類型是高山灌叢、高山草甸和亞高山草甸；冬蟲夏草次適宜生長的植被類型是高山草原；闊葉林、針葉林、荒漠、裸露巖石等無植被地段很難發(fā)現(xiàn)寄主昆蟲和冬蟲夏草分布,是不適宜分布的植被類型或地段。根據(jù)西藏植被類型圖及冬蟲夏草對植被的依賴性將冬蟲夏草資源的植被適宜性分布劃分為3個區(qū)域(圖3)。

2.1.3 冬蟲夏草土壤類型適宜性分析

土壤是蝠蛾昆蟲一生中最主要的活動場所,土壤質(zhì)地、土壤溫濕度是否適合幼蟲生活十分關鍵。同時,土壤類型的分布無形中也影響著蝠蛾幼蟲取食植物的分布,這些植物的嫩根或塊莖是蝠蛾幼蟲生存的食物基礎,土壤的適宜性影響寄主幼蟲的生長發(fā)育、種群數(shù)量、食物來源,甚至影響冬蟲夏草的產(chǎn)生。在西藏那曲地區(qū)調(diào)查發(fā)現(xiàn)冬蟲夏草多分布于比如縣、索縣、聶榮縣、嘉黎縣大部、那曲縣大部、巴青縣大部的亞高山草甸土和高山草甸土中。有研究認為最適宜蝠蛾幼蟲生長的土壤為高山草甸土和高山灌叢草甸土；其次為草甸混合土和流石灘,高山棕色和暗棕色林土也發(fā)現(xiàn)有少量寄主分布；其他土壤不適宜蝠蛾幼蟲生長。高山草甸土在西藏高原分布廣泛,主要發(fā)育在森林郁閉線以上的高山帶或無林山原高寒草甸植被下。這些土壤類型主要包括黑氈土、草氈土,生長在黑氈土和草氈土上的植被主要是高山嵩草、矮生嵩草等多種嵩草和薹草組成的高寒草甸,群落結構簡單,草層低矮,劃為適宜土壤類型。薄草氈土、黑褐土、赤紅壤、黃壤、褐土、灰褐土、赤紅壤等土壤類型中有少量寄主昆蟲幼蟲和冬蟲夏草資源分布,劃為次適宜土壤類型；而寒鈣土、暗寒鈣土、沼澤土、寒凍土、冷鈣土、實質(zhì)土等不適宜寄主昆蟲幼蟲生長和冬蟲夏草產(chǎn)生,因此劃為不適宜土壤類型。通過對冬蟲夏草分布區(qū)土壤類型及冬蟲夏草對土壤的適應性綜合分析,將冬蟲夏草資源的土壤類型適宜性空間分布劃分為3 個區(qū)域(圖4)。

圖4 西藏高原冬蟲夏草土壤適宜性分布Fig.4 Suitability distribution of Ophiocordyceps sinensis based on the soil types in the Tibetan Plateau

2.1.4 冬蟲夏草年平均降水量適宜性分析

在自然狀況下,濕度和降水是影響冬蟲夏草的主要氣候因子,特別是雪水對冬蟲夏草生長發(fā)育影響很大,干旱更不利于冬蟲夏草的生長發(fā)育。空氣濕度影響蝠蛾昆蟲的發(fā)育期長短,過低的空氣濕度會抑制蝠蛾昆蟲幼蟲發(fā)育,較高的空氣濕度促進生長發(fā)育,因此較多的降水和保持較高的空氣濕度有利于冬蟲夏草的形成。據(jù)調(diào)查,在3—4月降雪量較大的年份,冬蟲夏草產(chǎn)量相對較高；降雪少的年份,當年的冬蟲夏草產(chǎn)量相對較低。在降雪少的年份,如果空氣濕度低于50%時,會導致冬蟲夏草生長發(fā)育不良,子實體瘦小干癟。因此,降水量是影響冬蟲夏草產(chǎn)量最重要的限制性氣候因子之一,甚至降水的形態(tài)也決定著冬蟲夏草的生長發(fā)育及產(chǎn)量。根據(jù)對冬蟲夏草主產(chǎn)區(qū)、次產(chǎn)區(qū)和非產(chǎn)區(qū)年平均降水量、冬蟲夏草單位草原面積產(chǎn)量及各縣分布狀況的分析,結合38個氣象站1981—2010年年平均降水量數(shù)據(jù),初步將年平均降水量為451—600mm的區(qū)域劃為冬蟲夏草生長的適宜區(qū)；年平均降水量為301—450mm和601—750mm的區(qū)域劃為冬蟲夏草生長的次適宜區(qū)；年平均降水量≤300mm和≥750mm的區(qū)域劃為冬蟲夏草生長的不適宜區(qū)。依據(jù)以上標準,將冬蟲夏草資源降水量適宜性空間分布劃分為3 個區(qū)域(圖5)。

圖5 西藏高原冬蟲夏草年平均降水量適宜性分布Fig.5 Suitability distribution of Ophiocordyceps sinensis based on annual average precipitation in the Tibetan Plateau

2.1.5 冬蟲夏草年平均氣溫適宜性分析

氣溫直接影響冬蟲夏草的生長和發(fā)育,寒冷而濕潤,晝夜溫差大有利于冬蟲夏草繁衍。產(chǎn)區(qū)冬蟲夏草適應的氣溫有細微的差別,對西藏那曲蟲草進行生態(tài)調(diào)查發(fā)現(xiàn)產(chǎn)區(qū)平均氣溫為-1.9—4.0℃[6]。根據(jù)對冬蟲夏草主產(chǎn)區(qū)、次產(chǎn)區(qū)和非產(chǎn)區(qū)冬蟲夏草單位草原面積產(chǎn)量及各縣分布狀況的分析,結合38個氣象站1981—2010年年平均氣溫數(shù)據(jù),將年平均氣溫為-2.9—5.0℃的區(qū)域劃為冬蟲夏草生長的適宜區(qū)；-5.9—3.0℃和5.1—8.0℃的區(qū)域劃為冬蟲夏草生長的次適宜區(qū)；≤-6.0℃和≥8℃的區(qū)域劃為冬蟲夏草生長的不適宜區(qū),因此,將冬蟲夏草資源年平均氣溫適宜性空間分布按以上條件劃分為3 個區(qū)域(圖6)。

圖6 西藏高原冬蟲夏草年平均氣溫適宜性分布Fig.6 Suitability distribution of Ophiocordyceps sinensis based on annual average temperature in the Tibetan Plateau

2.2 冬蟲夏草綜合適宜性區(qū)劃

宏觀上分析得出決定冬蟲夏草分布的客觀條件主要為海拔高度、植被類型、土壤類型、年平均降水量及年平均氣溫,經(jīng)過對以上柵格數(shù)據(jù)重新分類將不適宜區(qū)賦值為0,將次適宜區(qū)賦值為1,將適宜區(qū)賦值為2,利用ArcGIS 軟件在同一投影下將所有柵格數(shù)據(jù)像元設為一致大小重采樣,用柵格計算工具將以上5個刪格數(shù)據(jù)圖層疊加。本文依據(jù)海拔高度確定冬蟲夏草海拔適宜性分布區(qū)域,作為第一圖層；在此圖層的基礎上,疊加植被類型確定的冬蟲夏草植被類型適宜性分布圖,兩圖層疊加時,取兩個圖層同一位置像元中的最小值作為新疊加后圖層的像元值,以此類推決定新圖層的原則進行,最終得出冬蟲夏草的綜合適宜性空間分布圖(圖7),該區(qū)劃結果與海拔高度、植被類型、土壤類型3要素的區(qū)劃結果相比,有了很大改進,更能客觀的反映冬蟲夏草的實際分布狀況,說明引進氣候因素對冬蟲夏草區(qū)劃是十分必要的。研究認為那曲地區(qū)的比如縣、索縣,昌都市的邊壩縣、丁青縣、類烏齊縣、江達縣、洛隆縣、察雅縣、貢覺縣等區(qū)域多處在適宜區(qū),該區(qū)主要處在那曲地區(qū)東部和昌都市大部；那曲地區(qū)的那曲縣、嘉黎縣、巴青縣、聶榮縣,拉薩市的當雄縣及其他產(chǎn)冬蟲夏草縣的部分區(qū)域為較適宜區(qū)；日喀則市的大部分地區(qū)、林芝市和山南市的低海拔地區(qū)為不適宜區(qū)。將西藏地區(qū)(市)行政邊界疊加到冬蟲夏草適宜性區(qū)劃圖上,分別提取6個地市冬蟲夏草不同適宜性等級分布面積(表2),發(fā)現(xiàn)在出產(chǎn)冬蟲夏草的地區(qū)(市)中,那曲地區(qū)、昌都市分布面積占當?shù)禺a(chǎn)冬蟲夏草縣域國土面積的52.1%—59.3%,山南市、日喀則市、拉薩市和林芝市分布面積占當?shù)禺a(chǎn)冬蟲夏草縣域國土面積的13.7%—26.1%,整個西藏高原的冬蟲夏草分布面積平均占當?shù)禺a(chǎn)冬蟲夏草縣域國土面積的34.5%。適宜區(qū)面積比例在各個地區(qū)差異也比較大，其中昌都市比例最大，其次是那曲地區(qū)和林芝市，山南市、日喀則市和拉薩市這一比例極低；山南市、日喀則市、拉薩市和林芝市不適宜區(qū)面積比例都特別高,這與該地區(qū)海拔低適合森林等茂密植被生長或主要為農(nóng)區(qū)有關,多因為海拔、降水量或植被所限制,很難產(chǎn)生冬蟲夏草。

圖7 西藏高原冬蟲夏草綜合適宜性空間分布Fig.7 Comprehensive suitability distribution of Ophiocordyceps sinensis in the Tibetan Plateau

地區(qū)Region產(chǎn)蟲草縣的國土面積Growth of Ophiocordyceps sinensis County area of land分布面積Areas/km2 /占該蟲草區(qū)面積比Percent/%適宜區(qū)Suitable areas次適宜區(qū)Sub-suitable areas不適宜區(qū)Unsuitable areas那曲地區(qū)6584034985.313076746.733157547.96昌都市1077152153519.994234439.314383640.70林芝市11321214171.252432321.488747277.26山南市59038340.06820213.895080286.05日喀則市53181400.08725213.644588986.29拉薩市1965800.00513326.111452573.89分布區(qū)域Distribution area418644265246.3411802128.1927409965.47

2.3 對比驗證分析

根據(jù)農(nóng)牧廳草原監(jiān)理站近11年(2005—2015年)調(diào)查資料分析,西藏的冬蟲夏草產(chǎn)量中那曲地區(qū)約占45.6%,昌都市約占33.2%,其他4個地市合計約占21.2%。利用西藏自治區(qū)草原監(jiān)理站2005—2015年各縣冬蟲夏草產(chǎn)量和草原面積數(shù)據(jù),計算獲得各縣單位草原面積產(chǎn)量(圖8),發(fā)現(xiàn)那曲地區(qū)東部的比如縣、索縣、巴青縣,昌都市的丁青縣、類烏齊縣,山南市的加查縣和林芝市的朗縣單位草原面積產(chǎn)量高；邊壩縣、嘉黎縣、洛隆縣、昌都縣等單位草原面積產(chǎn)量較高；其他縣域單位草原面積產(chǎn)量相對較低。這一分布事實與冬蟲夏草綜合適宜性空間分布較為一致,證明了區(qū)劃具有較強的科學性和客觀性,適宜性空間分布未能反映出加查、朗縣這一適宜分布小區(qū)域,可能是由于這兩縣地理位置偏南,海拔高度平均為3200m左右,海拔限制性區(qū)劃使得該區(qū)域部分被劃分為不適宜區(qū),實際上西藏高原南北適宜性海拔高度是有差別的,如果在藏東南地區(qū)和藏北高原采取不同的海拔區(qū)劃標準就可能解決這個問題。利用經(jīng)緯度提取156個冬蟲夏草采挖鄉(xiāng)鎮(zhèn)政府駐地高程,研究發(fā)現(xiàn)最低鄉(xiāng)鎮(zhèn)海拔為2304m,最高鄉(xiāng)鎮(zhèn)海拔為5074m,其中海拔為2304—3000m的鄉(xiāng)鎮(zhèn)5個,占3.2%；3001—3500m的鄉(xiāng)鎮(zhèn)18個,占11.5%；3501—4000m的鄉(xiāng)鎮(zhèn)45個,占28.8%；4001—4500m的鄉(xiāng)鎮(zhèn)58個,占37.2%；4501—5000m的鄉(xiāng)鎮(zhèn)29個,占18.6%；5000m以上的鄉(xiāng)鎮(zhèn)只有1個,占0.7%,海拔在3000—5000m之間冬蟲夏草產(chǎn)區(qū)鄉(xiāng)鎮(zhèn)占96.1%,可以初步判斷5000m和3000m是西藏高原產(chǎn)生冬蟲夏草的海拔上下限,但西藏地域遼闊,小氣候環(huán)境復雜而多樣,會出現(xiàn)各地海拔上下限有所差異,需要和實際調(diào)查結合。研究認為那曲地區(qū)東部產(chǎn)區(qū)大部分處在海拔4000—5000m之間,與陳仕江等[14]調(diào)查結果一致。整個西藏高原來看,冬蟲夏草產(chǎn)區(qū)海拔呈向東、向南降低趨勢,那曲地區(qū)產(chǎn)區(qū)海拔最高,其次為昌都市、拉薩市、山南市和日喀則市,而林芝市海拔最低(圖9)。

圖8 各縣單位面積草原冬蟲夏草產(chǎn)量Fig.8 The yield of Ophiocordyceps sinensis in the area per unitgrassland of each county

圖9 冬蟲夏草采挖點鄉(xiāng)鎮(zhèn)政府所在地海拔高度Fig.9 Altitude of Town government seat of Ophiocordyceps sinensisin the Tibetan Plateau

3 結論與討論

研究認為西藏高原那曲地區(qū)的比如縣、索縣,昌都市的邊壩縣、丁青縣、類烏齊縣、江達縣、洛隆縣、察雅縣、貢覺縣部分產(chǎn)區(qū)處在適宜區(qū)；那曲地區(qū)的那曲縣、嘉黎縣、巴青縣、聶榮縣,拉薩市的當雄縣及其他產(chǎn)區(qū)處在次適宜區(qū)；日喀則市的大部分地區(qū)、林芝市和山南市的低海拔地區(qū)為不適宜區(qū)?？傮w上西藏高原冬蟲夏草資源主要分布在那曲地區(qū)東部和昌都市西部的高海拔地區(qū),在林芝市、山南市、拉薩市、日喀則市的部分高海拔地區(qū)也有分布,各生長區(qū)域適宜區(qū)、次適宜區(qū)和不適宜區(qū)均呈斑塊狀鑲嵌分布。該適宜性區(qū)劃能比較真實的反應西藏冬蟲夏草分布格局,通過冬蟲夏草采挖縣產(chǎn)量和適宜度的驗證很好地證明了該方法的可行性；通過在那曲縣、比如縣、索縣冬蟲夏草產(chǎn)地實地調(diào)研,也印證了這一分布事實,認為這一區(qū)劃方法比較科學。

本文首次將氣候因子引入到冬蟲夏草適宜性區(qū)劃,是對西藏高原地形復雜,地貌組合區(qū)域差異大,氣候差異大的情況下進行冬蟲夏草資源適宜性區(qū)劃的有效嘗試；初步明確了西藏高原冬蟲夏草分布的年平均降水量和年平均氣溫適宜性的閾值范圍,如果經(jīng)過多點驗證,可更加精確冬蟲夏草分布的氣候因子閾值。由于本研究使用的環(huán)境因子圖層數(shù)據(jù)(氣候因子、植被、土壤等)的空間分辨率較粗,疊加后產(chǎn)生一定的誤差,對區(qū)劃圖的準確性有一定影響。應該縮小研究范圍,應用分辨率更高的數(shù)據(jù),采取更加科學的數(shù)據(jù)處理方法減小誤差,提高制圖的準確度。另外,由于生態(tài)因子的重要性有區(qū)別,對不同生態(tài)因子分級處理,賦予不同的權重,加權運算獲得的區(qū)劃結果可能會更加科學。同時考慮寄主昆蟲的生長發(fā)育、取食植物的豐富度以及冬蟲夏草產(chǎn)生易受到坡度、坡向等地貌因子,日照、空氣濕度、土壤溫濕度等氣候因子的影響,將這些因素考慮進區(qū)劃會更加科學；同樣冬蟲夏草分布也受到河流、湖泊、凍土、冰蓋等分布的影響,加強這方面的研究也十分必要。

本區(qū)劃標準與李仕江等[6]等在西藏高原那曲地區(qū)的生態(tài)調(diào)查較為一致,在那曲地區(qū)海拔下限要比大部分產(chǎn)地高100—600m,海拔上限大多數(shù)產(chǎn)區(qū)均為5000m。隨著經(jīng)濟的發(fā)展,道路修建、牧區(qū)居民點建設、過度放牧、人為破壞草皮、濫砍濫伐森林灌木、亂扔生活垃圾、不合理采挖等人為干擾均會影響冬蟲夏草資源的分布；氣溫升高,雪線抬升,冬蟲夏草核心分布區(qū)海拔明顯上升,冬蟲夏草產(chǎn)區(qū)面積明顯變小或發(fā)生位置遷移。近年來對冬蟲夏草寄主蝠蛾種類的多樣性進行了系統(tǒng)研究[28- 31],由于冬蟲夏草蝠蛾昆蟲的多樣性,不同寄主的生態(tài)特性不同,會存在不同的適應性和分布。對于某一種蝠蛾昆蟲分布,必須從該蝠蛾昆蟲的生物學和生態(tài)學特性,結合外界生態(tài)環(huán)境開展研究。今后研究冬蟲夏草生長適宜性應從冬蟲夏草蝠蛾昆蟲和蟲草菌生長發(fā)育的各個環(huán)節(jié)本身來研究,加強各種環(huán)境因子對冬蟲夏草產(chǎn)生和不同發(fā)育階段影響的研究,加強氣候變化對冬蟲夏草產(chǎn)生的影響研究；同時加強利用遙感技術對冬蟲夏草資源分布區(qū)的植被等生態(tài)環(huán)境監(jiān)測,有效保護冬蟲夏草資源。

參考文獻(References):

[1] 梁宗琦. 中國真菌志第三十二卷蟲草屬. 北京: 科學出版社, 2007: 1- 1.

[2] Winkler D. Caterpillar fungus (Ophiocordycepssinensis) production and sustainability on the Tibetan Plateau and in the Himalayas. Asian Medicine, 2009, 5(2): 291- 316.

[3] Li Y, Wang X L, Jiao L,Jiang Y, Li H, Jiang S P, LhosumtseiringN, Fu S Z, Dong C H, Zhan Y, Yao Y J. A survey of the geographic distribution ofOphiocordycepssinensis. The Journal of Microbiology, 2011, 49(6): 913- 919.

[4] 中國科學院西北高原生物研究所, 青海省藥品檢驗所. 中國蟲草:歷史·資源·科研[M]. 陜西科學技術出版社, 2008: 54- 58.

[5] 董應才, 羅都強. 我國的冬蟲夏草研究. 國土與自然資源研究, 1996, (3): 47- 51.

[6] 陳仕江, 尹定華, 丹增, 扎西, 扎瑪次仁. 中國西藏那曲冬蟲夏草的生態(tài)調(diào)查. 西南農(nóng)業(yè)大學學報, 2001, 23(4): 289- 292, 296- 296.

[7] 劉兆紅, 李玉玲. 玉樹州冬蟲夏草資源與分布. 草業(yè)與畜牧, 2006, (12): 34- 36.

[8] 黃天福, 陳仕江, 傅善全, 羅慶明. 康定冬蟲夏草菌優(yōu)勢寄主昆蟲的生態(tài)型的研究. 時珍國藥研究, 1996, 7(3): 178- 179.

[9] 肖生榮, 施至用, 陳慶濤. 蟲草的研究——Ⅰ. 生境考察和形態(tài)描述. 微生物學通報, 1983, (1): 5- 6.

[10] 郭相, 劉蓓, 馬紹賓, 桂明英. 云南冬蟲夏草生態(tài)環(huán)境調(diào)查及生物學特性分析. 中國食用菌, 2008, 27(6): 8- 11.

[11] 顧龍云, 安旺盛. 甘肅省甘南藏族自治州蟲草考察初報. 西北植物學報, 1987, 7(4): 261- 265.

[12] 袁峰. 冬蟲夏草居群譜系地理與適生區(qū)分布研究[D]. 昆明: 云南大學, 2015.

[13] 羅玉秀. 冬蟲夏草的研究現(xiàn)狀. 青海大學學報: 自然科學版, 2003, 21(2): 38- 40, 52- 52.

[14] 陳仕江, 尹定華, 李黎, 扎西, 孫健華, 扎瑪次仁. 西藏那曲地區(qū)冬蟲夏草資源及分布. 中藥材, 2000, 23(11): 673- 675.

[15] 高葵. 川西冬蟲夏草分布及生境中種子植物區(qū)系研究. 綿陽師范學院學報, 2008, 27(11): 81- 87.

[16] 朱印酒, 段雙全, 歐珠朗杰. 冬蟲夏草的研究進展. 中央民族大學學報: 自然科學版, 2009, 18(2): 27- 33.

[17] 雷豪清. 淺談玉樹州冬蟲夏草生長與水熱條件的關系. 青海草業(yè), 1995, (4): 19- 20, 36- 36.

[18] 刁治民. 青海冬蟲夏草資源及生物學特性的初步研究. 生物學雜志, 1996, 13(2): 20- 22.

[19] 楊大榮, 李朝達, 舒暢, 楊躍雄. 蟲草蝠蛾分布與高寒草甸植物的關系. 西南農(nóng)業(yè)學報, 1992, 5(2): 68- 73.

[20] 楊大榮, 彭艷瓊, 陳吉岳, 曹永強, 楊培. 中國冬蟲夏草分布格局與環(huán)境變化對其分布的影響. 中國草地學報, 2010, 32(S1): 22- 27.

[21] 沈發(fā)榮, 楊躍雄, 楊大榮, 春生, 魯自, 斯那都吉, 喻潤清. 云南的冬蟲夏草初步研究. 微生物學通報, 1988, 15(2): 49- 51.

[22] 周東平. 經(jīng)濟聚焦: 冬蟲夏草, 如何把“根”留住. 北京: 人民網(wǎng)-《人民日報》, [2010-07- 27].http://cpc.people.com.cn/GB/64093/64387/12256253.html.

[23] 金云翔, 徐斌, 楊秀春, 覃志豪, 高懋芳, 呂海燕, 朱立博. 青藏高原那曲地區(qū)冬蟲夏草資源分布空間分析方法. 生態(tài)學報, 2010, 30(6): 1532- 1538.

[24] 李芬, 吳志豐, 徐翠, 徐延達, 張林波. 三江源區(qū)冬蟲夏草資源適宜性空間分布. 生態(tài)學報, 2014, 34(5): 1318- 1325.

[25] 西藏自治區(qū)土地管理局. 西藏自治區(qū)土地利用. 北京: 科學出版社, 1992: 165- 198.

[26] 張新時. 西藏植被的高原地帶性. 植被學報, 1978, 20(2): 140- 149.

[27] 李明森. 西藏土壤資源的特點及其開發(fā)利用. 資源科學, 1984, (2): 1- 6.

[28] 王曉亮. 冬蟲夏草及其寄主昆蟲的多樣性研究[D]. 北京: 中國科學院研究生院, 2012.

[29] 張姝. 冬蟲夏草菌群體遺傳學及其與寄主昆蟲之間的進化關系研究[D]. 太原: 山西大學, 2013.

[30] 全慶梅. 冬蟲夏草菌與其寄主昆蟲的遺傳多樣性及協(xié)同進化研究[D]. 上海: 同濟大學, 2014.

[31] 吳長奎. 蘭坪蟲草及其寄主蝠蛾居群遺傳學和適生區(qū)分布研究[D]. 昆明: 云南大學, 2016.