基于最大熵模型的藏藥翼首草生態(tài)適宜性研究

唐策+文檢+張?chǎng)?蘇錦松+謝彩香+張藝

[摘要] 研究翼首草的生態(tài)適宜性，為翼首草GAP種植選址和區(qū)域發(fā)展提供參考。通過(guò)中國(guó)數(shù)字標(biāo)本館、文獻(xiàn)和實(shí)地考察等收集67個(gè)翼首草分布點(diǎn)，利用Maxent和GIS技術(shù)對(duì)翼首草進(jìn)行適宜性分析。發(fā)現(xiàn)翼首草潛在分布區(qū)域主要分布在西藏東部、四川西部、青海南部、云南西北部、甘肅南部；對(duì)翼首草生長(zhǎng)和分布貢獻(xiàn)率大的主要生態(tài)因子為海拔（62%）、最暖季風(fēng)降水量（14.4%）、降水量變異系數(shù)（7.2%）、最干季度平均溫度（3.5%）、土壤導(dǎo)電率（3%）、陽(yáng)離子交換總量（2.4%）、溫度季節(jié)性變化標(biāo)準(zhǔn)差（2.2%）。應(yīng)用Maxent模型對(duì)翼首草進(jìn)行生態(tài)適宜性研究，其準(zhǔn)確度高，可為翼首草野生撫育及GAP種植選址提供科學(xué)依據(jù)和參考。

[關(guān)鍵詞] 藏藥； 翼首草； Maxent模型； 生態(tài)適宜性

[Abstract] In order to study the ecology suitability of Pterocephalus hookeri， and provide a reference for GAP planting location and regional development， the Maxent model and GIS technology were used to investigate ecology suitability regions for P. hookeri based on the distribution points collected from Chinese virtual herbarium， the references and field trips. The potential distribution areas mainly concentrated in the eastern Tibet， western Sichuan， southern Qinghai， northwest Yunnan， and southern Gansu. There were 7 major environmental factors to have obvious influence on ecology suitability distributions of P. hookeri， including altitude （contribution rate of 62%）， precipitation of warmest quarter （contribution rate of 14.4%）， coefficient of variation of precipitation seasonality （contribution rate of 7.2%）， mean temperature of driest quarter （contribution rate of 3.5%）， the electrical conductivity of top and sub-soil （contribution rate of 3%）， the total exchangeable bases in the top- and subsoil （contribution rate of 2.4%） and SD of temperature seasonality （contribution rate of 2.2%）. The study of the ecological suitability regionalization of P. hookeri based on Maxent model can provide scientific basis for the selection of artificial planting base and GAP planting location.

[Key words] traditional Tibetan herb； Pterocephalus hookeri； Maxent model； potential distribution areas

翼首草（藏語(yǔ)名：榜孜毒烏）為川續(xù)斷科匙葉翼首草Pterocephalus hookeri（C.B.Clarke）Heck 的干燥全草，是藏醫(yī)臨床常用藥材[1-2]。翼首草為多年生草本植物，主要分布于西藏東部、云南西北部、青海南部、甘肅南部及四川西部和北部等，生長(zhǎng)在海拔1 800～5 700 m的山坡、草地、草甸、林間等[3-5]。翼首草具有抗炎鎮(zhèn)痛[6]、抗類(lèi)風(fēng)濕性關(guān)節(jié)炎[7]、抗腫瘤[8]、抗菌[9]、保肝[10]、免疫調(diào)節(jié)[9-10]等多種藥理作用，是多種藏成藥的重要原料。

翼首草是2015年版《中國(guó)藥典》（一部）[1]及衛(wèi)生部《藏藥部頒標(biāo)準(zhǔn)》[11]中收載品種。近年來(lái)，由于醫(yī)藥市場(chǎng)藥材價(jià)格飛漲，對(duì)其采取無(wú)限制、掠奪式、毀滅性的濫采濫挖，生態(tài)環(huán)境遭到嚴(yán)重破壞，導(dǎo)致野生翼首草藥材出現(xiàn)瀕危[12-13]。課題組前期對(duì)青藏高原藏藥瀕危物種進(jìn)行調(diào)查，使用《中國(guó)珍稀瀕危植物》的等級(jí)劃分方法，建議將匙葉翼首花增加為青藏高原五省（區(qū)）藏區(qū)瀕危藏藥物種，建議瀕危等級(jí)列為二級(jí)（稀有）[14]。故需加強(qiáng)對(duì)翼首草生態(tài)保護(hù)與人工栽培相關(guān)研究，通過(guò)翼首草的人工引種栽培與野生撫育，以提高翼首草的產(chǎn)量，并緩解其資源匱乏現(xiàn)狀。

Maxent模型可以根據(jù)物種當(dāng)前存在的不完整數(shù)據(jù)信息和環(huán)境變量間緊密的關(guān)系，找出物種分布規(guī)律的最大熵，從而對(duì)物種的潛在分布進(jìn)行預(yù)測(cè)，是目前預(yù)測(cè)物種潛在分布比較好的模型[15]。因此，本文以藏藥翼首草為研究對(duì)象，通過(guò)收集翼首草分布點(diǎn)的經(jīng)緯度信息，綜合氣候、地形、土壤等相關(guān)生態(tài)因子，運(yùn)用Maxent模型研究其在我國(guó)的適宜性區(qū)域，分析影響翼首草生長(zhǎng)的主要生態(tài)因子，確立翼首草適宜分布區(qū)域和生態(tài)條件，為藏藥翼首草的引種栽培以及產(chǎn)業(yè)發(fā)展提供參考。

1 數(shù)據(jù)來(lái)源

1.1 藏藥翼首草分布點(diǎn)經(jīng)緯度數(shù)據(jù)的獲取 本文收集翼首草在我國(guó)的分布數(shù)據(jù)主要來(lái)源于文獻(xiàn)檢索及中國(guó)數(shù)字植物標(biāo)本館（http：//www.cvh.org.cn/），以及對(duì)翼首草的實(shí)地采樣調(diào)查。標(biāo)本信息包括物種名、采樣點(diǎn)經(jīng)緯度，野外調(diào)查信息記錄包括經(jīng)緯度、海拔等。

1.2 生態(tài)因子數(shù)據(jù)的獲取 本文使用35個(gè)生態(tài)因子，各因子具體信息見(jiàn)表1。其中包括19個(gè)生物氣候變量（Worldclim； http：//www.worldclim.org/），1950—2000年監(jiān)測(cè)數(shù)據(jù)的平均值，精度2.5弧分；15個(gè)土壤數(shù)據(jù)來(lái)自世界土壤數(shù)據(jù)庫(kù)HWSD（http：//webarchive.iiasa.ac.at/Research/LUC/External-World-soil-database/HTML/）；1個(gè)地形數(shù)據(jù)和地理基礎(chǔ)數(shù)據(jù)來(lái)源于地理空間數(shù)據(jù)云（http：//www.gscloud.cn/sources/？cdataid=302& pdataid=10）。上述數(shù)據(jù)坐標(biāo)系為WGS84，圖層?xùn)鸥翊笮〖s為1 km2，溫度數(shù)值（℃）為實(shí)際數(shù)值的10倍。

2 方法

2.1 數(shù)據(jù)預(yù)處理 采用ENVI 4.8和Arcmap 10.0對(duì)土壤數(shù)據(jù)的柵格大小，地理坐標(biāo)進(jìn)行處理，得到與氣候數(shù)據(jù)相同的空間參數(shù)，同時(shí)對(duì)土壤數(shù)據(jù)進(jìn)行格式轉(zhuǎn)換，得到與氣候數(shù)據(jù)和海拔數(shù)據(jù)相同的數(shù)據(jù)格式（.asc），并直接加載于MaxEnt軟件（Version 3.3.3 k，http：//www.cs.princeton.edu/～schapire/ maxent/）。

2.2 Maxent模型參數(shù)設(shè)置 將生態(tài)環(huán)境圖層和翼首草分布點(diǎn)經(jīng)緯度數(shù)據(jù)（.csv格式）加載到Maximum entropy modeling of species geographic distributions，Version 3.3.3k版軟件中，設(shè)置參數(shù)運(yùn)行建模，設(shè)分布數(shù)據(jù)的25%被隨機(jī)抽取作為測(cè)試集（testdata），其余作為訓(xùn)練集，最大迭代次數(shù)是105。設(shè)置刀切法檢驗(yàn)權(quán)重，Maxent自定義設(shè)置ROC（receive operating characteristic）評(píng)價(jià)曲線(xiàn)，輸出文件類(lèi)型選為.asc，其他參數(shù)為軟件默認(rèn)值。

2.3 主要生態(tài)因子的選擇 35個(gè)生態(tài)因子和翼首草67個(gè)分布點(diǎn)的經(jīng)緯度數(shù)據(jù)導(dǎo)入Maxent模型軟件，按照2.2項(xiàng)設(shè)置進(jìn)行數(shù)據(jù)處理，見(jiàn)表1。根據(jù)Maxent模型多次迭代計(jì)算結(jié)果，選取總貢獻(xiàn)率≥90%的生態(tài)因子作為主要生態(tài)因子。根據(jù)各主要生態(tài)因子的響應(yīng)曲線(xiàn)，獲得各生態(tài)因子的適宜值范圍。

2.4 適宜性區(qū)域的劃分 應(yīng)用ArcGIS 10.0軟件，對(duì)Maxent模型運(yùn)行結(jié)果圖進(jìn)行處理。采用人工分級(jí)的方法，劃分翼首草的適宜性分布區(qū)域。

3 結(jié)果與分析

3.1 翼首草分布點(diǎn)數(shù)據(jù) 對(duì)中國(guó)數(shù)字植物標(biāo)本館查詢(xún)，和對(duì)高原地區(qū)翼首草實(shí)地調(diào)查，收集到翼首草分布點(diǎn)104個(gè)，主要分布于四川、西藏、青海、云南等省區(qū)，見(jiàn)圖1。其中野外采集的標(biāo)本經(jīng)成都中醫(yī)藥大學(xué)民族醫(yī)藥學(xué)院張藝研究員鑒定為川續(xù)斷科匙葉翼首草P. hookeri植物。

3.2 模型的準(zhǔn)確度 Maxent模型一般根據(jù)ROC曲線(xiàn)下圍成的AUC值對(duì)預(yù)測(cè)結(jié)果進(jìn)行精度評(píng)測(cè)，AUC的取值范圍為0.5～1，值越大表示模型的準(zhǔn)確性越高。一般來(lái)說(shuō)AUC大于0.85的模擬驗(yàn)證結(jié)果就可以采納。當(dāng)AUC大于0.90時(shí)，模型模擬價(jià)值較高，可以較準(zhǔn)確的反應(yīng)物種的潛在分布[16]。根據(jù)Maxent模型預(yù)測(cè)結(jié)果的ROC曲線(xiàn)，見(jiàn)圖2，訓(xùn)練集和測(cè)試集的AUC均為0.995，表明Maxent模型預(yù)測(cè)結(jié)果準(zhǔn)確度高，可用于進(jìn)一步分析翼首草的分布區(qū)域。

3.3 主要生態(tài)因子 根據(jù)翼首草在我國(guó)的分布數(shù)據(jù)和環(huán)境數(shù)據(jù)，利用MaxEnt 軟件分析，獲得各生態(tài)因子對(duì)翼首草生長(zhǎng)適宜度的貢獻(xiàn)率，其中貢獻(xiàn)率大于0的有17個(gè)，其他生態(tài)因子貢獻(xiàn)率為0，見(jiàn)表1。對(duì)翼首草生長(zhǎng)貢獻(xiàn)率較大的因子為海拔（62%）、最暖季節(jié)降水量（14.4%）、降水量變異系數(shù)（7.2%）、最干季度平均溫度（3.5%）、土壤導(dǎo)電率（3%），總的貢獻(xiàn)率為90.1%。這5個(gè)主要生態(tài)因子中影響最大的為海拔，因此，海拔可能是影響翼首草分布的限制性因子；該結(jié)果也證實(shí)了藏藥翼首草為一種生長(zhǎng)在缺氧、高海拔高原環(huán)境的一種草本植物。

3.4 主要生態(tài)因子適宜值范圍 Maxent模型軟件運(yùn)行結(jié)果中各生態(tài)因子響應(yīng)曲線(xiàn)，即存在概率與生態(tài)

因子閾值范圍的關(guān)系，見(jiàn)圖3。根據(jù)翼首草的響應(yīng)曲線(xiàn)，劃分主要生態(tài)因子的適宜值范圍（存在概率大于0.1的范圍）：海拔高度范圍2 300～5 600 m，最適宜值在4 200 m左右；最暖季度降水量范圍：200～1 400 mm，最適宜值在300 mm左右；降水量變異系數(shù)范圍75～160，最適宜值在95左右；最干季度平均溫度范圍-9～10.5 ℃，最適宜值在-4 ℃左右；土壤導(dǎo)電率大于1%均可，最適宜值在2%左右；陽(yáng)離子交換總量范圍50～460 cmol·kg^-1，最適宜值在85 cmol·kg^-1左右。

3.5 翼首草適宜性區(qū)域區(qū)劃 根據(jù)翼首草適宜性區(qū)域區(qū)劃結(jié)果，見(jiàn)圖4。采用人工分級(jí)的方法，將翼首草分布區(qū)域劃分為5個(gè)區(qū)間，即相似度<10%、相似度10%～30%、相似度30%～50%、相似度50%～70%、相似度>70%。翼首草潛在分布區(qū)域主要在青藏高原及其周邊區(qū)域，以橫斷山脈附近最佳。翼首草潛在分布區(qū)域以西藏、四川、青海、云南、甘肅5個(gè)省區(qū)分布面積最大，陜西、山西、寧夏、貴州4省區(qū)也有少量分布。

相似度10%～30%的區(qū)域，面積為70.9萬(wàn)km2，分布在四川的德陽(yáng)市、廣元市、雅安市、攀枝花市，西藏的昂仁縣、措勒縣、班戈縣、申扎縣、謝通門(mén)縣、安多縣、聶榮縣、薩嘎縣、嘉黎縣，青海的雜多縣、稱(chēng)多縣、甘德縣、瑪沁縣、同德縣、門(mén)源縣、樂(lè)都縣、貴德縣，云南的楚雄縣、大理、昆明、曲靖，甘肅的天水、成縣、定西、臨夏、平?jīng)隹h、蘭州市、武威，貴州畢節(jié)市、六盤(pán)水，寧夏固原，山西太原、晉中等地區(qū)。

相似度30%～50%的區(qū)域，面積為17.9萬(wàn)km2，主要在甘肅的夏河縣、卓尼縣較為集中，其他省份分布較為分散。

相似度50%～70%的區(qū)域，面積為47.8萬(wàn)km2，分布在四川的阿壩藏族羌族自治州的大部分區(qū)域；西藏林芝地區(qū)的察隅縣、米林縣、墨脫縣，日喀則地區(qū)的崗巴縣、亞?wèn)|縣、薩迦縣、江孜縣，昌都地區(qū)丁青縣、八宿縣；青海的斑馬縣、久治縣、玉樹(shù)縣，甘肅的碌曲縣、瑪曲縣、舟曲縣；云南的德欽縣、麗江市。

相似度>70%的區(qū)域是翼首草最適宜生長(zhǎng)區(qū)域，面積為20.2萬(wàn)km2主要分布在四川、西藏、云南、青海4個(gè)省（自治區(qū)）。主要縣城有：四川的甘孜藏族自治州的大部分區(qū)域，以及涼山彝族自治州的木里縣和九龍縣；西藏昌都地區(qū)的貢覺(jué)縣、芒康縣、江達(dá)縣、察雅縣、妥壩縣、類(lèi)烏齊縣，林芝地區(qū)波密縣，日喀則地區(qū)定結(jié)縣、定日縣、拉孜縣、日喀則市；青海的囊謙縣；云南的中甸縣。

4 討論

本文基于氣候、土壤和地形相關(guān)的35個(gè)生態(tài)因子，運(yùn)用Maxent模型對(duì)翼首草進(jìn)行了潛在分布預(yù)測(cè)。從預(yù)測(cè)結(jié)果顯示，翼首草的潛在適宜區(qū)主要分布在我國(guó)的青藏高原及其周邊區(qū)域，以橫斷山脈附近最佳。翼首草潛在分布區(qū)域主要分布在西藏東部、四川西部、青海南部、云南西北部、甘肅南部，與目前翼首草實(shí)際分布區(qū)域一致。同時(shí)各適宜指數(shù)區(qū)間劃分也基本符合翼首草在我國(guó)的分布現(xiàn)狀，因此，本次研究結(jié)果準(zhǔn)確度較高。

藏藥材是藏醫(yī)藥發(fā)展的物質(zhì)基礎(chǔ)，近年來(lái)隨著藏醫(yī)藥市場(chǎng)的開(kāi)拓和生產(chǎn)的工業(yè)化，藏醫(yī)藥產(chǎn)業(yè)得到了迅猛發(fā)展，藏藥材的需求量越來(lái)越大，但目前90%以上的藏藥材供給僅靠野生資源[17]。隨著采集量的增大，不僅破壞了藏藥材種群的自然更新，也破壞著野生藏藥材的生長(zhǎng)環(huán)境，由此引發(fā)的“掠奪式”采挖，致使一些常用藏藥材原料陷入瀕危或滅

絕的困境。據(jù)了解，以翼首草為處方原料的藏藥材資源嚴(yán)重不足，翼首草的年消費(fèi)均在1 000 t左右。因此，本文基于Maxent模型和GIS技術(shù)對(duì)翼首草進(jìn)行適宜性區(qū)域劃分研究，分析影響其生長(zhǎng)的主要生態(tài)因子以及適宜范圍，對(duì)其進(jìn)行分布可能性進(jìn)行等級(jí)劃分，為藏藥翼首草野生撫育提供科學(xué)依據(jù)，也為GAP種植基地的選址提供參考。

[參考文獻(xiàn)]

[1] 中國(guó)藥典.一部[S].2015：381.

[2] 賈敏如，張藝.中國(guó)民族藥辭典[M].北京：中國(guó)醫(yī)學(xué)科技出版社，2016：678.

[3] 中國(guó)科學(xué)院西北高原生物研究所.藏藥志[M].西寧：青海人民出版社，1991：191.

[4] 青海省藥品檢驗(yàn)所，青海省藏醫(yī)藥研究所.中國(guó)藏藥.第1卷[M].上海：上海科技出版社，1996：273.

[5] 中國(guó)科學(xué)院中國(guó)植物志編輯委員會(huì).中國(guó)植物志.第73卷. 第1分冊(cè)[M].北京：科學(xué)出版社，1986：69.

[6] Zhang L， Hu J J， Lin J W， et al. Anti-inflammatory and analgesic effects of ethanol and aqueous extracts of Pterocephalus hookeri （C.B. Clarke） Hoeck[J]. J Ethnopharmacol，2009， 123（3）： 510.

[7] 王曉東，鄧志文，彭曉華，等.藏藥然降多吉膠囊對(duì)大鼠佐劑性關(guān)節(jié)炎細(xì)胞因子及炎癥介質(zhì)的影響[J].中國(guó)中藥雜志，2004，29（12）：1186.

[8] 郭晨旭，吳迎春，朱元章，等.藏族藥翼首草正丁醇部位提取物體外抑制人肝癌Hep3B細(xì)胞增殖和侵襲轉(zhuǎn)移[J].中國(guó)實(shí)驗(yàn)方劑學(xué)雜志，2015，21（3）：100.

[9] 甄梓娟，徐元江，廖志華，等.藏藥匙葉翼首草及其同屬植物的研究進(jìn)展[J].中藥材，2016，39（1）：223.

[10] 郭晨旭，朱國(guó)福.藏藥翼首草化學(xué)成分及藥理作用研究進(jìn)展[J].世界中醫(yī)藥，2015，10（9）：1440.

[11] 中華人民共和國(guó)衛(wèi)生部藥典委員會(huì).中華人民共和國(guó)衛(wèi)生部藥品標(biāo)準(zhǔn)·藏藥.第1冊(cè)[S].北京：化學(xué)工業(yè)出版社，1995：138.

[12] 李隆云，占堆，衛(wèi)瑩芳，等.瀕危藏藥資源的保護(hù)[J].中國(guó)中藥雜志，2002，27（8）： 561.

[13] 何淑玲，馬令法，常毓巍，等.甘南州道地藏藥植物資源區(qū)域分布與保護(hù)對(duì)策初探[J].廣東農(nóng)業(yè)科學(xué)，2012，（6）：45.

[14] 趙彩云，劉歡，蘇錦松，等.青藏高原藏藥瀕危現(xiàn)狀及其資源保護(hù)策略[J].中國(guó)中藥雜志，2016，41（23）：4451.

[15] 文檢，呂秀梅，洪道鑫，等.基于Maxent模型的青藏高原大花紅景天生態(tài)適宜性分析[J].中國(guó)中藥雜志，2016，41（21）：3931.

[16] 吳慶明，王磊，朱瑞萍，等.基于MAXENT模型的丹頂鶴營(yíng)巢生境適宜性分析——以扎龍保護(hù)區(qū)為例[J].生態(tài)學(xué)報(bào)，2016，36 （12）：1.

[17] 李軒豪，趙彩云，劉悅，等.青藏高原藏藥資源與品質(zhì)評(píng)價(jià)研究進(jìn)展[J].中國(guó)中藥雜志，2016，41（4）：562.

[責(zé)任編輯 呂冬梅]