基于MaxEnt和ArcGIS的定西市蒙古黃芪適宜性區(qū)劃研究

王莉 韋翡翡 成希 趙文龍 晉玲

中圖分類號 R-l

文獻標(biāo)志碼A

文章編號 1001-0408（2020）03-0321-04

DOI 10.6039/j.issn.1001-0408.2020.03.13

摘要 目的：對定西市蒙古黃芪進行生態(tài)適宜性區(qū)劃研究。方法：以定西市1001批蒙古黃芪為考察樣本（以每個有蒙古黃芪栽培的自然村為采集單位進行樣本采集），從定西市經(jīng)濟作物技術(shù)推廣站中獲取其經(jīng)、緯度信息，并從《中藥資源空間信息網(wǎng)格數(shù)據(jù)庫》中獲取其55項環(huán)境生態(tài)因子（包括氣候、地形、土壤等）信息。結(jié)合經(jīng)、緯度信息和環(huán)境生態(tài)因子信息，以75%樣本為訓(xùn)練集建立最大熵模型（MaxEnt模型），篩選出主要生態(tài)因子，并以另外25%樣本為驗證集進行模型驗證。隨后，應(yīng)用地理信息系統(tǒng)（Arc-GIS）對蒙古黃芪的適宜生長區(qū)進行劃分。結(jié)果：建立的MaxEnt模型預(yù)測效果良好（訓(xùn)練集和驗證集樣本受試者工作特征曲線下面積分別為0.970、0.968）。共篩選出海拔、降水量、溫度等8個主要生態(tài)因子（貢獻率總和為98.90%）;結(jié)合分析系統(tǒng)發(fā)現(xiàn)，海拔在1800～2650m、4月降水量在25～50mm、最冷月最低溫在-16～-8℃、最濕月降水量在95～110mm、溫度季節(jié)性變化標(biāo)準差在70～80、12月平均氣溫在-6～-3℃、10月降水量在30～50mm、12月降水量在0～10mm為定西市蒙古黃芪適宜生長的環(huán)境參數(shù)。結(jié)合ArcGIS分析結(jié)果，在定西市境內(nèi)蒙古黃芪普遍適宜生長;在臨洮縣西北部、安定區(qū)北部、通渭縣南部、隴西縣南部少量地區(qū)、岷縣南部少量地區(qū)以及渭源縣、漳縣、岷縣交界處不適宜生長。結(jié)論：本研究建立的模型分析結(jié)果與實際考察符合，可為定西市蒙古黃芪的種植區(qū)域規(guī)劃提供參考。

關(guān)鍵詞 蒙古黃芪;定西市;適宜性區(qū)劃;地理信息系統(tǒng)

黃芪為蒙古黃芪[Astragalus membranaceus （Fisch.）Bge.var.mongholicus（Bge.） Hsiao]或莢膜黃芪[Astraga-lus membranaceus （Fisch.）Bge.]的干燥根，具有補氣升陽、固表止汗、利水消腫、生津養(yǎng)血、行滯通痹、托毒排膿、斂瘡生肌的功效[1-2]。黃芪是中國傳統(tǒng)中藥，也是中醫(yī)最常用中藥之一，首載于《神農(nóng)本草經(jīng)》，并列為上品?，F(xiàn)代藥理研究表明，黃芪具有抗腫瘤、抗氧化、保肝等作用，對免疫系統(tǒng)、心血管系統(tǒng)、神經(jīng)系統(tǒng)等均有保護作用[3-4]，特別是在治療糖尿病和心血管等疾病上有顯著療效[5-7]。

定西市位于甘肅中部，地處黃土高原、青藏高原和西秦嶺交匯地帶，通稱“隴中”，是古絲綢之路重鎮(zhèn)和新歐亞大陸橋必經(jīng)之地，也是甘肅省內(nèi)蒙古黃芪主產(chǎn)地之一。其轄下有安定區(qū)及通渭、隴西、渭源、臨洮、漳縣、岷縣等1區(qū)6縣，并且其中的隴西縣素有“中國黃芪之鄉(xiāng)”的美稱。目前，蒙古黃芪從種子萌發(fā)到高產(chǎn)的栽培種植都有成熟的技術(shù)支持，但膜莢黃芪栽培質(zhì)量不好、產(chǎn)量不高[8-11]，故現(xiàn)代栽培主要以蒙古黃芪為主。楊長花等[12]對4個主要黃芪產(chǎn)地的黃芪質(zhì)量進行比較后發(fā)現(xiàn)，不同產(chǎn)地黃芪中有效成分含量有差異，其中以內(nèi)蒙古、甘肅等地產(chǎn)黃芪中有效成分含量較高。

最大熵模型（MaxEnt）可在最大熵原理的基礎(chǔ)上較為精確地推斷和預(yù)測不同物種的生態(tài)適宜區(qū)域[13]。地理信息系統(tǒng)（ArcGIS）利用的是生態(tài)學(xué)原理和數(shù)據(jù)分析，對物種的潛在分布區(qū)域進行預(yù)測并將成果展示出來，在中藥相關(guān)領(lǐng)域也得到了廣泛應(yīng)用[14-15]。故本研究利用MaxEnt和ArcGIS相結(jié)合的方式，對定西市蒙古黃芪的資源分布進行整合，并對其進行生態(tài)適宜區(qū)預(yù)測，為蒙古黃芪的合理種植提供理論支持。

l 儀器與信息來源

1.1 儀器

Rin0 530HCx型便攜式衛(wèi)星導(dǎo)航儀（北京佳明航電科技有限公司）。

1.2 樣本信息

蒙古黃芪樣本由定西市經(jīng)濟作物技術(shù)推廣站于2018年10月進行批量采集（原則是以每個有蒙古黃芪栽培的自然村為采集單位進行樣本采集），并利用便攜式衛(wèi)星導(dǎo)航儀記錄每個樣本的海拔和經(jīng)、緯度信息。共采集到1001個蒙古黃芪樣本，其中安定區(qū)34個、臨洮縣108個、隴西縣181個、岷縣346個、通渭縣35個、渭源縣219個、漳縣78個。所有樣本均由甘肅中醫(yī)藥大學(xué)藥學(xué)院晉玲教授鑒定為蒙古黃芪真品。

1.3 生態(tài)因子信息

從《中藥資源空間信息網(wǎng)格數(shù)據(jù)庫》中下載相關(guān)生態(tài)因子信息，共55項。其中氣候變量數(shù)據(jù)是根據(jù)我國氣象觀測站1950-2000年的氣候觀測數(shù)據(jù)插值成空間分辨率為1kmxlkm的格網(wǎng)數(shù)據(jù);土壤類型數(shù)據(jù)采用FA0-90土壤分類系統(tǒng)，是根據(jù)第二次全國土壤調(diào)查《1：100萬中華人民共和國土壤圖》制成;植被類型數(shù)據(jù)來源于中科院植物所《1：100萬中華人民共和國植被圖》中的植被亞類數(shù)據(jù)[13]。

2 方法

2.1 MaxEnt模型參數(shù)及生態(tài)因子選擇

將采集到的1001個蒙古黃芪樣本的經(jīng)、緯度信息和生態(tài)因子數(shù)據(jù)添加到MaxEnt對應(yīng)的對話框中，按照要求設(shè)定參數(shù)，隨機選取75%的蒙古黃芪樣本數(shù)據(jù)作為訓(xùn)練集，以另外25%的蒙古黃芪樣本數(shù)據(jù)作為測試集，設(shè)定最大迭代次數(shù)為106，其余采用模型的默認選項。運行程序，剔除結(jié)果中貢獻率為O的生態(tài)因子，并再次運行程序，最終篩選出對蒙古黃芪生長有影響的生態(tài)因子。利用刀切法（反映不同生態(tài)因子在影響蒙古黃芪的生長適宜度中所占的權(quán)重）檢驗權(quán)重，研究各生態(tài)因子對蒙古黃芪生長的貢獻率[16]。刀切圖結(jié)合MaxEnt預(yù)測結(jié)果綜合得出影響蒙古黃芪生長的生態(tài)因子。

2.2 模型精確性驗證

受試者工作特征曲線（ROC曲線）分析法在物種潛在分布預(yù)測模型評價中得到了廣泛應(yīng)用，用ROC曲線下面積（AUC，其大小可作為預(yù)測模型準確度的衡量指標(biāo)）進行模型精確度預(yù)測[17]。AUC取值范圍為0.5～1，數(shù)值越接近1說明該預(yù)測的結(jié)果越好，其模型預(yù)測的結(jié)果就越準確。當(dāng)AUC在0.6以下時，說明該模型預(yù)測失敗;當(dāng)AUC在0.6～0.9之間時，說明該模型預(yù)測結(jié)果一般;當(dāng)AUC在0.9以上時，說明該模型預(yù)測結(jié)果非常好。

2.3 定西市蒙古黃芪適宜區(qū)域劃分標(biāo)準

利用ArcGIS軟件Arctoolbox模塊中的ConversionTools將MaxEnt軟件迭代運行得到的最終結(jié)果轉(zhuǎn)換為柵格數(shù)據(jù)，用“按掩膜提取”工具將定西市的預(yù)測結(jié)果從全國預(yù)測結(jié)果中提取出來。在提取結(jié)果的屬性表中計算得到最大值（max）、最小值（min）、平均值（μ）、標(biāo)準差（σ）、平均值與標(biāo)準差相加（μ+σ）、平均值與標(biāo)準差相減（μ-σ）的值。以[0，μ-σ]為不適宜區(qū)，[μ-σ，μ]為次適宜區(qū)，[μ-σ，μ+σ]為適宜區(qū)，可得到其區(qū)劃等級分布。

3 結(jié)果

3.1 MaxEnt模型的預(yù)測結(jié)果分析

按上述方法建立MaxEnt模型重復(fù)運行5次后，得到的訓(xùn)練集樣本的AUC為0.970，測試集樣本的AUC為0.968，表明MaxEnt模型對定西市蒙古黃芪生態(tài)適宜區(qū)的預(yù)測效果良好。刀切圖中顯示，有16個生態(tài)因子貢獻率不為O，對蒙古黃芪生長有影響見圖1（只變量：生態(tài)因子;不變量：生長適宜度;全部變量：有貢獻率的生長因子）。定西市蒙古黃芪ROC曲線見圖2。

3.2 主要生態(tài)因子的確定

生態(tài)因子貢獻率大小會影響物種分布，本次參與模型建立的生態(tài)因子有55個，MaxEnt運行結(jié)束后，有16個生態(tài)因子貢獻率不為O。其中8個生態(tài)因子的貢獻率總和占比達到98.90%，包括海拔（38.07%）、4月降水量（32.58%）、最冷月最低溫（6.03%）、最濕月降水量（5.79%）、溫度季節(jié)性變化標(biāo)準差（5.71%）、12月平均氣溫（4.79%）、10月降水量（3.50%）、12月降水量（2.43%），為蒙古黃芪適宜區(qū)的主要影響生態(tài)因子。上述生態(tài)因子對蒙古黃芪生長影響的適宜范圍分別為海拔1800～2650m、4月降水量25～50mm、最冷月最低溫-16～-8℃、最濕月降水量95～110mm、溫度季節(jié)性變化標(biāo)準差70～80、12月平均氣溫-6～-3℃、10月降水量30～50mm、12月降水量0～10mm。影響定西市蒙古黃芪的8種主要生態(tài)因子響應(yīng)圖見圖3。

3.3 定西市蒙古黃芪適宜區(qū)域劃分

對定西市蒙古黃芪生態(tài)適宜區(qū)進行綜合評價后將區(qū)劃等級劃分為3類：以黃芪采樣點生境適宜度值O～0.42為不適宜區(qū)，0.42～0.50為次適宜區(qū)，0.50--0.58為適宜區(qū)。結(jié)果顯示，蒙古黃芪在定西市普遍適宜生長，在臨洮縣西北部、安定區(qū)北部、通渭縣南部、隴西縣南部少量地區(qū)、岷縣南部少量地區(qū)以及渭源縣、漳縣、岷縣交界處為蒙古黃芪不適宜區(qū)，具體分布見表1。

4 討論

通過MaxEnt運算，此次模型的AUC大于0.9，可信度高。在55個生態(tài)因子中運算得到16個貢獻率不為O的生態(tài)因子，并篩選對模擬結(jié)果影響較大的8個生態(tài)因子，分別是海拔（38.07%）>4月降水量（32.58%）>最冷月最低溫（6.03%）>最濕月降水量（5.79%）>溫度季節(jié)性變化標(biāo)準差（5.71%）>12月平均氣溫（4.79%）>10月降水量（3.50%）>12月降水量（2.43%），累積貢獻率達98.90%。結(jié)合已知蒙古黃芪生長所需條件發(fā)現(xiàn)，海拔對蒙古黃芪生長影響最大，需要保障在1800～2650m;在4月蒙古黃芪長苗期、最濕月（即花果期）以及10、12月地上部分枯萎時期需要較多水分，而在12月、最冷月需要提供適宜溫度以使其正常發(fā)育。這與實際調(diào)查的蒙古黃芪生長所需條件相吻合。

定西市臨洮縣西北部（紅旗鄉(xiāng)大部、中鋪鎮(zhèn)西南部、太石鎮(zhèn)西部、辛店鎮(zhèn)西部及東北部、新添鎮(zhèn)西部、上營鄉(xiāng)北部），安定區(qū)北部（門碌鄉(xiāng)、石峽溝鄉(xiāng)、魯家溝鎮(zhèn)交界處大部），通渭縣南部（常家河鎮(zhèn)中部至南部、李家店中部至南部、襄南鄉(xiāng)南部、碧玉鄉(xiāng)中部至南部少量地區(qū)、雞川鎮(zhèn)南部、隴川鄉(xiāng)北部），隴西縣南部少量地區(qū)（文峰鎮(zhèn)中南部至鞏昌鎮(zhèn)交界處），渭源縣南部少量地區(qū)（峽城鄉(xiāng)南部、會川鎮(zhèn)南部與五竹鎮(zhèn)交界處），漳縣西北部和東部少量地區(qū)（金鐘鎮(zhèn)西部、南部和殪虎橋鄉(xiāng)、大草灘鄉(xiāng)交界處、武陽鎮(zhèn)南部、武當(dāng)鄉(xiāng)南部、新寺鎮(zhèn)北部），岷縣南部少量地區(qū)（秦許鄉(xiāng)西部至中部、寺溝鄉(xiāng)南部至中部、麻于川鄉(xiāng)南部至中部、茶埠鎮(zhèn)西南部、寺溝鄉(xiāng)、禾馱鄉(xiāng)交界處、禾馱鄉(xiāng)南部和申都鄉(xiāng)交界處、申都鄉(xiāng)南部、閭井鎮(zhèn)南部至鎖龍鄉(xiāng)交界處、鎖龍鄉(xiāng)南部與馬塢鄉(xiāng)交界處）為蒙古黃芪的不適宜區(qū)，主要是由于該地區(qū)內(nèi)有高山，海拔過高不適宜蒙古黃芪生長。臨洮縣西北部（紅旗鄉(xiāng)東部、中鋪鎮(zhèn)大部、太石鎮(zhèn)中部、辛店鎮(zhèn)中部、新添鎮(zhèn)大部、八里鋪鎮(zhèn)與洮陽鎮(zhèn)交界處），安定區(qū)大部（除內(nèi)官鎮(zhèn)、香泉鎮(zhèn)、高峰鄉(xiāng)），隴西縣中部少量地區(qū)（鞏昌鎮(zhèn)、和平縣、云田鄉(xiāng)、渭陽鄉(xiāng)、首陽鎮(zhèn)東部、永吉鄉(xiāng)北部），通渭縣大部（除華家?guī)X鄉(xiāng)、北城鋪鄉(xiāng)大部、隴陽東西部、隴山東西部、馬營鎮(zhèn)南北部、什川鄉(xiāng)南北部、第三鋪鄉(xiāng)北部、榜羅鎮(zhèn)中西部），岷縣東部少量地區(qū)（鎖龍鄉(xiāng)北面、馬塢鄉(xiāng)和閭井鎮(zhèn)交界處）為蒙古黃芪次適宜區(qū)，其他地區(qū)都適宜蒙古黃芪生長。隴西縣大部分地區(qū)都為蒙古黃芪適宜生長區(qū)，這也與其實際地理考證相吻合。

綜上，蒙古黃芪在定西大部分地區(qū)適宜種植，適宜區(qū)研究可為提升當(dāng)?shù)刂兴幉漠a(chǎn)業(yè)結(jié)構(gòu)、人工種植栽培以及可持續(xù)利用研究提供依據(jù)。本研究可為后期定西市各縣（區(qū)）蒙古黃芪質(zhì)量評價及其品質(zhì)區(qū)劃研究提供基礎(chǔ)。

參考文獻

[1]中國科學(xué)院中國植物志編輯委員會，中國植物志[M]北京：科學(xué)出版社，1993：133.

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

[3]林紅強，楊娜，王涵，等.黃芪的化學(xué)成分、藥理活性及臨床應(yīng)用研究進展[J]世界最新醫(yī)學(xué)信息文摘，2018，18（38）：45-47、49

[4]聶娟，謝麗華，馬港圓，等，中藥黃芪的化學(xué)成分及藥理作用研究進展[J]湖南中醫(yī)雜志，2018，34（7）：228-231.

[5]尤良震，林逸軒，方朝暉，等，黃芪甲苷治療糖尿病及其并發(fā)癥藥理作用研究進展[J]中國中藥雜志，2017，42（24）：4700-4706

[6]張瑜，中藥黃芪的心血管藥理作用及臨床應(yīng)用研究[J].世界最新醫(yī)學(xué)信息文摘，2017，17（74）：47-48

[7]田楓，李佶操，唐紹微，等，基線等比增減設(shè)計法優(yōu)選黃芪注射液聯(lián)合燈盞細辛注射液抗大鼠腦缺血再灌注損傷的最佳配比[J]中國藥房，2019，30（14）：1885-1889

[8]程瑤蒙古黃芪硬實種皮解剖特征及休眠破除方法研究[D]蘭州：甘肅農(nóng)業(yè)大學(xué)，2017.

[9]張訓(xùn)中草藥黃芪栽培技術(shù)[J]南方農(nóng)業(yè)，2017，11（2）：9-10

[10]黃文全黃芪高產(chǎn)栽培技術(shù)[J]農(nóng)村科技，2017（9）：65-66.

[11]譚根堂黃芪栽培技術(shù)[J].農(nóng)村百事通，2018（9）：30-32

[12]楊長花，程芬，李慧，等4個主要黃芪產(chǎn)地黃芪質(zhì)量比較[J]西部中醫(yī)藥，2019，32（1）：27-31.

[13]邵揚，葉丹，歐陽臻，等薄荷的生境適宜性區(qū)劃及品質(zhì)區(qū)劃研究[J]中國中藥雜志，2016，41（17）：3169-3175

[14]孫成忠，郝振同，朱壽東，等，地理信息技術(shù)在中藥資源生產(chǎn)區(qū)劃研究中的應(yīng)用[J]測繪通報，2016（12）：96-99、123

[15]李足意，謝宗明，孑L祥柏，等，瀏陽河中上游地區(qū)中藥資源區(qū)劃研究[J]中國衛(wèi)生產(chǎn)業(yè)，2019，16（19）：126-127

[16]嚴輝，張小波，朱壽東，等當(dāng)歸藥材生產(chǎn)區(qū)劃研究[J]中國中藥雜志，2016，41（17）：3139-3147.

[17]王茹琳，高曉清，王閆利，等.基于MaxEnt的非洲橘硬薊馬在全球及中國的潛在分布區(qū)預(yù)測[J].中國農(nóng)學(xué)通報，2014，30（28）：315-320.