崔治家，雷豐順，盧有媛，晉玲，馬毅，王振恒

(甘肅中醫(yī)藥大學藥學院中(藏)藥資源研究所，甘肅 蘭州 730000)

·中藥農(nóng)業(yè)·

達烏里秦艽潛在適生分布區(qū)研究△

崔治家，雷豐順，盧有媛，晉玲*，馬毅，王振恒

(甘肅中醫(yī)藥大學藥學院中(藏)藥資源研究所，甘肅 蘭州 730000)

目的：研究影響達烏里秦艽適宜性分布的生態(tài)因子，預測達烏里秦艽在我國的潛在適生分布區(qū)。方法：利用達烏里秦艽在中國的154個地理分布記錄，結合28項環(huán)境因子，應用最大熵模型(MaxEnt)研究影響達烏里秦艽適宜性分布的生態(tài)因子，結合ArcGIS軟件預測達烏里秦艽全國及各省潛在適生分布區(qū)。結果：影響達烏里秦艽適宜性生長的主要環(huán)境因子為海拔、1月最高溫和1月降水量；達烏里秦艽在我國10個省的潛在適生性較高，適生區(qū)總面積達577 517.4 km2，占全國總面積的6.02%；甘肅省達烏里秦艽潛在適生區(qū)面積最大，青海次之。結論：研究結果對達烏里秦艽原產(chǎn)地資源保護和野生撫育具有指導作用，可以為達烏里秦艽栽培選地和生產(chǎn)區(qū)劃研究提供科學依據(jù)。

達烏里秦艽；潛在適生區(qū)；最大熵模型；ArcGIS

達烏里秦艽(GentianadahuricaFisch.)又稱小秦艽，其干燥根為中藥秦艽，具有祛風濕、清濕熱、止痹痛、退虛熱之功[1]；其花為藏藥解吉那保，具有清熱、消炎之功效[2]。達烏里秦艽喜涼爽、濕潤氣候，耐寒，主要分布于內(nèi)蒙古、陜西、甘肅、寧夏、山西、青海、新疆、四川、西藏及河北等地[3-4]。由于長期采挖，其野生資源日趨減少，不能滿足市場需求。達烏里秦艽地理分布及生境、形態(tài)特征、化學成分及藥理作用等方面已有相關研究報道[5-8]，但是長期依靠野生資源使得研究其生長的環(huán)境因子及適宜生長區(qū)工作開展較少。

隨著統(tǒng)計模型和GIS技術的興起，生態(tài)學原理、地理信息系統(tǒng)(Geographic Information System，GIS)結合環(huán)境變量進行物種潛在適生分布區(qū)預測得到了應用。目前，常用的生態(tài)位模型主要有6種，分別為GARP模型、ENFA模型、CLIMEX模型、MAXENT模型、BIOCLIM模型、DOMAIN模型[9-13]。其中，MaxEnt模型在農(nóng)作物適宜區(qū)預測及藥用植物潛在生境分布中得到了廣泛應用并取得了良好效果。5種生態(tài)位模型預測黃頂菊的研究證明MaxEnt模型的模擬精度最好[14]。本研究應用Maxent模型對達烏里秦艽的潛在適生分布區(qū)進行預測，研究影響其分布的主要環(huán)境因子，結合ArcGIS獲得達烏里秦艽的潛在適生區(qū)，為達烏里秦艽栽培生產(chǎn)選地和生產(chǎn)區(qū)劃研究提供科學依據(jù)。

1 材料與方法

1.1 達烏里秦艽分布數(shù)據(jù)的獲取與處理

查閱全球生物多樣性信息平臺(GBIF，http：//www.gbif.org/)、中國數(shù)字植物標本館(CVH，http：//www.cvh.org.cn/cms/)、中國自然保護區(qū)資源平臺(http：www.papc.cn/)、教學標本標準化整理整合與資源共享平臺(http：//mnh.scu.edu.cn/)，并參考文獻資料，共獲得達烏里秦艽分布信息1114條，沒有經(jīng)緯度坐標的分布點以最小可查行政單位為準，剔除重復數(shù)據(jù)及無法查詢數(shù)據(jù)，最終獲得154個采樣點信息，見圖1。采樣點數(shù)據(jù)按種名、分布點經(jīng)度和緯度順序排列后以.csv格式保存。

圖1 達烏里秦艽的樣點分布圖

1.2 環(huán)境變量數(shù)據(jù)獲取

從世界氣候數(shù)據(jù)庫(http：//www.worldclim.org)中獲得的1950～2000年的28項環(huán)境因子，精度為2.5弧分×2.5弧分(≈25 km)，詳細信息見(表1)。環(huán)境因子數(shù)據(jù)用ArcMap中的格式轉(zhuǎn)換工具轉(zhuǎn)換為ESRI ASCII格式。

1.3 MaxEnt模型設定

參考川貝母潛在分布區(qū)預測方法[15]，測試數(shù)據(jù)比例為25%，最大迭代次數(shù)106，設置ROC曲線評價模型準確度，響應曲線和刀切法檢驗環(huán)境因子權重，其余為模型默認設置。

1.4 適宜區(qū)劃分

根據(jù)Maxent模型預測達烏里秦艽在中國的分布概況及資料數(shù)據(jù)，將達烏里秦艽適宜性分為四類，分別為：P<0.1不適生區(qū)，0.1≤P<0.4低適生區(qū)，0.4≤P<0.6中適生區(qū)，P≥0.6高適生區(qū)。將三級以上河流分布區(qū)域設為低適宜區(qū)，最終獲得達烏里秦艽在中國的潛在適生區(qū)分布圖；利用提取工具獲得其在各省的潛在適生區(qū)分布圖；利用重分類工具計算得到中適生區(qū)和高適生區(qū)面積。

表1 達烏里秦艽適生區(qū)研究所用環(huán)境變量

表1(續(xù))

2 結果與分析

2.1 MaxEnt模型預測精度

采用ROC曲線對Maxent模型預測結果進行精度檢驗，即用AUC值(ROC曲線下的面積)作為衡量指標[16]，AUC的取值范圍為0.5～1，其值越靠近1，說明預測結果越好，即MaxEnt模型預測結果越準確。AUC值的標準為：0.5～0.6(失敗)，0.6～0.7(較差)，07～0.8(一般)，0.8～0.9(好)，0.9～1(非常好)。本研究MaxEnt模型預測結果測試數(shù)據(jù)集和訓練數(shù)據(jù)集的AUC值分別為0.994和0.996(圖2)，說明Maxent模型預測達烏里秦艽潛在適生區(qū)分布結果具有較高的額準確度和可信度。

圖2 達烏里秦艽潛在分布預測的ROC曲線

2.2 主導環(huán)境因子的確定

根據(jù)刀切法的結果可知(圖3)，1月降雨量、海拔、最冷季度平均溫度、最冷季度降雨量、1月均溫、1月最高溫是影響其分布的主要環(huán)境因子，根據(jù)模型中各環(huán)境變量的貢獻率(表2)可知，海拔(38.7%)、1月降雨量(25.7%)、1月最高溫(13.9%)、最暖季度降水量(10.7%)是影響達烏里秦艽分布的主要環(huán)境因子，取兩者共有環(huán)境因子，即海拔、1月降雨量、1月最高溫是影響達烏里秦艽分布的主要環(huán)境因子。

2.3 達烏里秦艽適生區(qū)分布

由全國適宜分布圖(圖4)可知，達烏里秦艽的適宜分布區(qū)主要集中在甘肅中部及南部、青海及西藏東部、四川西北部及北部、寧夏南部、陜西北部、山西北部及中部地區(qū)、內(nèi)蒙古中南部及河北北部地區(qū)，新疆中西部等地也有低適宜區(qū)分布。各省適宜分布圖(圖5)顯示，甘肅省達烏里秦艽的適生區(qū)主要分布張掖、武威、天祝等地及中部和西南部的甘南藏族自治州、白銀、平?jīng)?、定西、蘭州、臨夏等地，東部慶陽環(huán)縣等地也有中適宜分布區(qū)；內(nèi)蒙古達烏里秦艽適生區(qū)主要集中在中南部烏蘭察布、呼和浩特等地；寧夏達烏里秦艽適生區(qū)主要集中在南部固原、海源等地和中部吳忠部分地區(qū)，沿賀蘭山也有少量適生區(qū)分布；青海省達烏里秦艽適生區(qū)分布在東部海北藏族自治州、海南藏族自治州、西寧、黃南藏族自治州、海東等地區(qū)和中部海西蒙古族藏族自治州；山西省除南部運城、晉城及臨汾部分地區(qū)外，其他地區(qū)均為達烏里秦艽適生區(qū)；陜西省達烏里秦艽適生區(qū)主要分布在北部榆林地區(qū)和延安地區(qū)，寶雞及西安部分地區(qū)有低適宜分布區(qū)；四川省達烏里秦艽高適生區(qū)主要分布在北部地區(qū)阿壩藏族羌族自治州；新疆中西部伊利哈薩克有達烏里秦艽中、低適宜分布區(qū)；河北省達烏里秦艽中、高適宜區(qū)主要分布在北部的張家口，承德北部的圍場等地為低適宜分布區(qū)；西藏達東部的昌都、妥壩、左貢等地區(qū)有烏里秦艽低適生區(qū)分布。

圖3 刀切法的環(huán)境變量重要性分析

環(huán)境因子貢獻率(%)環(huán)境因子貢獻率(%)海拔387Bio17最干季度降水量051月降雨量257Bio14最干月降水量031月最高溫139Bio1年均溫02Bio18最暖季度降水量107Bio19最冷季度降水量017月降雨量43Bio7年均溫變化范圍01Bio6最冷月最低溫25Bio3等溫性01bio4溫度季節(jié)性變化的標準差15Bio8最濕季度平均溫度01Bio11最冷季度平均溫度07Bio15降水量變異系數(shù)01Bio12年均降水量05

圖4 達烏里秦艽在中國適生區(qū)分布圖

注：A.甘省肅；B.內(nèi)蒙古自治區(qū)；C.寧夏回族自治區(qū)；D.青海省E.山西??；F.陜西??；G.四川??；H.新疆維吾爾自治區(qū)；I.河北??；J.西藏自治區(qū)。圖5 達烏里秦艽在各省區(qū)的潛在適生區(qū)分布圖

各省達烏里秦艽高適生區(qū)和中適生區(qū)面積及其占全國和全省面積比例統(tǒng)計見表3。由表3可知，達烏里秦艽在全國的適生區(qū)面積為577 517.4 km2，占國土總面積的6.02%；甘肅省達烏里秦艽適生面積最大為139 913.2 km2，占國土面積的1.46%，占甘肅省總面積的30.86%，其次為青海??；寧夏適生區(qū)面積為24 166.67 km2，占到全省總面積的36.4%，為達烏里秦艽適生區(qū)面積占有率最高的??；相比較而言，西藏和新疆達烏里秦艽適生區(qū)面積較少，在全國及全省土地適生區(qū)的占有率較低。

表3 各省達烏里秦艽主要適生分布區(qū)面積統(tǒng)計

3 討論

本研究利用Maxent模型和地理信息系統(tǒng)，并結合環(huán)境變量對達烏里秦艽在我國的適生分布區(qū)進行預測，獲得了達烏里秦艽在我國的適生分布區(qū)，ROC曲線顯示測試數(shù)據(jù)和訓練數(shù)據(jù)的AUC值均在0.9以上，說明模型預測結果的可信較高。預測結果表明，達烏里秦艽的高適生分布區(qū)主要分布在甘肅、青海、山西、內(nèi)蒙古、四川、陜西、寧夏、新疆、河北、西藏等地，這與文獻記載達烏里秦艽實際分布區(qū)相吻合。研究結果從理論上明確了影響達烏里秦艽分布的主要環(huán)境因子及在各省的適宜分布區(qū)，對于進一步開展達烏里秦艽原產(chǎn)地資源保護和野生撫育具有指導作用，為達烏里秦艽野生資源的收集及栽培選地和生產(chǎn)區(qū)劃研究提供了參考。同時，

研究中也存在不足，如僅考慮了氣候和地形因子的影響，沒有涉及土壤類型、植被類型及人類活動范圍以及社會的城鎮(zhèn)化發(fā)展對其生態(tài)適宜性分布的影響。

[1] 國家藥典委員會.中華人民共和國藥典：一部[S].北京：中國醫(yī)藥科技出版社，2015：270.

[2] 楊永昌.藏藥志[M].西寧：青海人民出版社，1991：11．

[3] 中國植物志編輯委員會.中國植物志[M].北京：科學出版社，1988：62-64.

[4] 吳興海，黃馨慧.中藥秦艽及其可持續(xù)利用途徑研究[J].資源保護和可持續(xù)發(fā)展,2002，4(6)：58-61.

[5] 司慶文，侯勤正，朱興福，等.達烏里秦艽的傳粉生物學研究[J].西北植物學報，2010，30(12)：2433-2436．

[6] 唐建寧，康建宏，許強，等.秦艽與小秦艽光合日變化的研究[J].西北植物學報，2006，26(4)：836-841.

[7] 陳千良，石張燕，張雅惠，等.小秦艽化學成分研究[J].中藥材，2011，34(8)：1214-1216.

[8] 鄔衛(wèi)東，楊來秀，王玉華，等.蒙藥小秦艽花總黃酮超聲波提取工藝及其抗氧化研究[J].內(nèi)蒙古師范大學學報，2016，45(4)：510-514.

[9] 張海龍.基于GARP生態(tài)位模型的珍稀植物桃兒七適生區(qū)與生境分析[J].江西農(nóng)業(yè)學報，2013，25(7)：112-115.

[10]楊瑞，張雅林，馮紀年.利用ENFA生態(tài)位模型分析玉帶鳳蝶和箭環(huán)蝶異地放飛的適生性[J].昆蟲學報，2008，51(3)：290-297.

[11]馬駿，胡學難，彭正強，等.基于CLIMEX模型的扶桑綿粉蚧在中國潛在地理分布預測[J].植物檢疫，2011，25(1)：5-8.

[12]何淑婷，白碧玉，但佳惠，等，基于MaxEnt的南丹參在中國的潛在分布區(qū)預測及適生性分析[J].安徽農(nóng)業(yè)科學，2014，42(8)：2311-2314.

[13]徐家文，史家浩，任強，等.基于BIOCLIM模型的扶桑綿粉蚧在中國的適生性分析[J].湖北農(nóng)業(yè)科學，2015，54(11)：2631-2633.

[14]曹向鋒，錢國良，胡白石，等.采用生態(tài)位模型預測黃頂菊在中國的潛在適生區(qū)[J].應用生態(tài)學報，2010，21(12)：3063-3069.

[15]王娟娟，曹博，白成科，等.基于Maxent和ArcGIS預測川貝母潛在分布及適宜性評價[J].植物研究，2014，34(5)：642-649.

[16]景鵬飛，武坤毅，龔曄，等.藥用植物細辛在中國的潛在適生區(qū)分布[J].植物分類與資源學報，2015，37(3)：349-356.

StudyofPotentialDistributionRegionsofGentianadahurica

CUIZhijia，LEIFengshun，LUYouyuan，JINLin*，MAYi，WANGZhenheng

(SchoolofPharmaceuticalSciences，InstituteofTraditionalChineseMedicineandTibetanMedicineResources，GansuUniversityofChineseMedicine，Lanzhou730000，China)

Objective：To study effect of ecological factors on distribution ofGentianadahuricaand predict the potential suitable distribution regions ofG.dahuricain China.Methods：By using 154 distribution recordings in China and 28 ecological factors and combined with ArcGIS，the effect of the ecological factors on the suitable distribution regions ofG.dahuricaand the potential suitable distribution regions in China and each province were studied by MaxEnt.Results：The main factors effect of the suitable growth ofG.dahuricawere altitude，the highest temperature and the precipitation in January.10 provinces in China existed higher suitable distribution areas forG.dahurica.The suitable distribution areas were 577 517.4 km2，accounting for 6.02% of the total area of the country.Gansu province had the largest areas for the suitable growth ofG.dahurica，and the second was Qinghai.Conclusion：The results could guide for resources native protection and wild nurture ofG.dahurica.Meanwhile，it could provide scientific basis for selecting cultivation fields and suitable production area ofG.dahurica.

GentianadahuricaFisch.；potential suitable distribution regions；MaxEnt；ArcGIS

甘肅中醫(yī)藥大學中青年基金(ZQN-2013-1)；中央本級重大增減支項目(20603020212，20603020110)；中醫(yī)藥公共衛(wèi)生專項 (財社[2011]76 號)；中醫(yī)藥行業(yè)科研專項(201207002)

] 晉玲，教授，研究方向：珍惜瀕危和大宗常用中藥資源可持續(xù)利用；E-mail：zyxyjl@163.com

10.13313/j.issn.1673-4890.2017.6.016

2016-11-10)

*[