中麻黃生態適宜性區劃研究

馬曉輝+盧有媛+黃得棟+朱田田+呂培霖+晉玲

[摘要] 利用最大熵Maxent模型并結合GIS技術對中麻黃在我國的生態適宜性進行預測，以尋找影響中麻黃適宜性分布的主要生態因子及適宜生長區域。對采集到38份中麻黃樣本和從CVH中獲得的116條中麻黃分布信息運用ArcGIS技術進行分析。利用MaxEnt模型對中麻黃在我國的生態適宜區進行預測。中麻黃MaxEnt模型ROC曲線訓練數據集和測試數據集的AUC依次分別是0.986，0.958，均大于0.9，趨近于1。中麻黃較高適宜度主要集中在甘肅西部、青海東部、新疆西北部、內蒙古中部和西藏東北部。模型運算得出的中麻黃生境適宜度具有很高的準確度和可信度，表明MaxEnt模型對中麻黃在中國潛在分布區的預測效果非常好。為野生資源的探測和栽培區域的選擇提供參考依據。

[關鍵詞] 中麻黃；MaxEnt模型；生態因子；生態適宜性

[Abstract] The study aims at predicting ecological suitability of Ephedra intermedia in China by using maximum entropy Maxent model combined with GIS，and finding the main ecological factors affecting the distribution of E. intermedia suitability in appropriate growth area. Thirty-eight collected samples of E. intermedia and E. intermedia and 116 distribution information from CVH information using ArcGIS technology were analyzed. MaxEnt model was applied to forecast the E. intermedia in our country′s ecology. E. intermedia MaxEnt ROC curve model training data and testing data sets the AUC value was 0.986 and 0.958，respectively，which were greater than 0.9，tending to be 1.The calculated E. intermedia habitat suitability by the model showed a high accuracy and credibility，which indicated that MaxEnt model could well predict the potential distribution area of E. intermedia in China.

[Key words] Ephedra intermedia；MaxEnt model；ecological factors；ecology suitability

中麻黃Ephedra intermedia Schrenk et C. A. Mey.為麻黃科Ephedraceae麻黃屬Ephedra多年生草本狀灌木，其干燥草質莖入藥，具有發汗散寒，宣肺平喘，利水消腫的功效，用于風寒感冒，胸悶喘咳，風水浮腫；其干燥根作麻黃根用，具有固表止汗的功效，用于自汗，盜汗[1]。均為常用中藥。由于長期處于野生狀態，經過無序濫采濫挖，中麻黃野生資源遭到極大破壞，被列入為第二批《國家重點保護野生藥材物種名錄》[2]。其在我國產于新疆、內蒙古、西藏、甘肅、青海等省區，以西北地區常見，生長在海拔數百米至2 000多米的干旱荒漠、沙灘地區及干旱的山坡或草地上[3]。目前，對于中麻黃區域性的分布概況及生態環境已有相關研究報道，但是其在全國的適宜生長概況未見相關研究[4-5]。若進行全國范圍的資源調查，耗時耗力，工作量難以估量。“3S”技術是一種現代化信息綜合處理技術，它能夠實時動態監測、評估產量、顯示地理分布，已被應用到物種分布和適宜性評價研究等方面[6]。其中，地理信息系統（GIS）與傳統區劃方法相結合在資源區劃、農業區劃、氣候區劃、災害區劃等研究中顯示出巨大優勢[7]。

中藥區劃是在中藥資源調查的基礎上正確評價影響中藥資源開發和中藥生產的自然條件及社會經濟條件，從而揭示中藥資源與中藥生產的地域分布規律。生態適宜性區劃作為中藥區劃的一部分，在評價物種的適宜分布及生長區時尤為重要[8]。而最大熵MaxEnt模型通過物種的分布數據和環境數據，找出物種分布的最大熵，從而對物種的分布進行預測[9]，在農作物適宜區預測、動物潛在生境評價、外來入侵物種風險評估和藥用植物潛在生境分布中得到廣泛應用并取得了良好效果[10-11]。

本文基于中麻黃的實際分布信息，結合氣候因子數據、土壤數據和植被類型，利用最大熵Maxent模型并結合GIS技術對中麻黃在我國的生態適宜性進行預測，以尋找影響中麻黃適宜性分布的主要生態因子及適宜生長區域，為野生資源的探測和栽培區域的選擇提供參考依據。

1 材料與方法

1.1 儀器

GPS（Rino 530HCx 便攜式衛星導航儀，北京佳明航電科技有限公司）。

1.2 樣本收集

查閱相關文獻明確中麻黃分布范圍，依據分布范圍，于2012—2016年采集中麻黃樣本，記錄經緯度、海拔和生境；查閱中國數字標本館（CVH） （http：//www.cvh.org.cn./cms），收集中麻黃的分布信息，沒有記錄經緯度坐標的分布點以最小可查行政單位為準。根據MaxEnt軟件要求，將采樣點按物種名、分布點經度和緯度順序儲存成后綴名為.csv文件。

1.3 中麻黃生長分布預測

1.3.1 生態因子數據的獲取 本研究所使用生態因子數據庫來源于“中藥資源空間信息網格數據庫”，包含氣候、土壤、地形和植被類型，涵蓋了影響中藥資源分布的所有生態因子數據。

1.3.2 MaxEnt模型參數設置 將中麻黃的采樣信息和環境數據按要求格式轉換后分別導入MaxEnt軟件，設置參數運行建模，具體參數如下：隨機測試比例為15%，最大迭代次數為106，設置響應曲線，MaxEnt自定義設置ROC評價曲線和刀切法（jackknife）檢驗權重，其他參數為軟件默認值。

1.3.3 模型預測評價采用 ROC工作曲線和曲線下面積（AUC值）進行模型精度評測。

1.3.4 主生態因子的選擇 第一步，根據中麻黃分布的GPS數據和15個生態因子分別在MaxEnt中的迭代運算結果，選取總貢獻率≥95%的生態因子作為供篩選的生態因子；第二步，分別提取中麻黃分布點的供篩選的生態因子數值，進行Spearman相關性分析，相關系數大于0.8的生態因子棄去貢獻率較小的一個；第三步，將第二步剩余的生態因子利用MaxEnt模型分別計算分析，得到生態因子對中麻黃4個種生態適宜性分布的貢獻率和刀切法獲得的主要生態因子權重分析，選取貢獻率大權重高的生態因子作為主生態因子，根據響應曲線獲取適宜性。

1.4 中麻黃在我國的生態適宜性區劃分析

利用GIS軟件MaxEnt模型計算結果進行疊加和地圖制作，生態適宜性區劃根據采樣點處生境適宜度值，按正態分布參數μ，δ進行劃分，[0，μ-δ]為不適宜區，[μ-δ，μ]為次適宜區，[μ，1]為適宜區，繪制出中麻黃的生態適宜性區劃分布圖。

2 結果與分析

2.1 樣本收集

本次研究共采集到中麻黃樣本38份；從CVH中獲得的中麻黃分布信息116條。

2.2 模型預測評價

ROC曲線下面積（即AUC）的大小作為模型預測準確度的衡量指標，其取值范圍為[0，1]，值越大表示模型判斷力越強[11]。AUC在0.5～0.6為失敗，0.6～0.7為較差，0.7～0.8為一般，0.8～0.9為好，0.9～1.0為非常好[10-11]。中麻黃MaxEnt模型ROC曲線訓練數據集和測試數據集的AUC值依次分別是0.986，0.958，均大于0.9，趨近于1。表明MaxEnt模型對中麻黃在中國潛在分布區的預測效果非常好，由模型運算得出的中麻黃生境適宜度具有很高的準確度和可信度。

2.3 主要生態因子的選擇

15個生態因子中6個對中麻黃生態適宜性有貢獻率，按1.3.4進行主生態因子的選擇結果見表1。

2.4 中麻黃在我國的生態適宜性區劃

中麻黃生態適宜性區劃見圖1～3，由圖可知，中麻黃較高適宜度主要集中在甘肅西部、青海東部、新疆西北部、內蒙古中部和西藏東北部；與實際資源調查來源的中麻黃實際分布位置相一致。

3 結論

中麻黃是瀕危物種，其野生資源破壞嚴重，雖然在甘肅、青海等省區也有栽培[12-13]，但其資源量依然不容樂觀，其野生資源保護和資源開發勢在必行。本研究通過中麻黃專項調查及網絡共享平臺數據收集了中麻黃154份分布信息，利用ArcGIS技術和MaxEnt模型對中麻黃在我國的生態適宜區進行預測，預測結果經ROC曲線評價AUC均在0.9以上，模型預測結果的準確度和可信度非常高。

試驗利用MaxEnt模型和ArcGIS技術相結合，較全面的考慮了各環境因子的影響，最終得出對中麻黃生長影響較大的環境因子及其適宜值范圍。綜合分析中麻黃的主環境因子，海拔、溫度和降水量對中麻黃生態適宜性分布均產生主要影響，從主環境因子生態適宜值分析，中麻黃適生海拔較高，對降水量的需求較少。

研究得出中麻黃在我國的生態適宜性區劃分布圖，其中中麻黃生長適宜區主要集中在甘肅西部、新疆西北部、西藏東南部、寧夏全境和內蒙古中部。試驗所得分布區與《中國植物志》中的中麻黃的分布區域基本吻合。

綜合整個研究，樣點信息基本覆蓋了中麻黃分布區，考慮了大多數中麻黃分布實際生態位，預測得出了影響生態適宜區的主要生態因子，并從理論上預測了中麻黃的生態適宜區和次適宜區，可以為中麻黃野生資源的收集、保護和栽培區域的選擇提供參考。

[參考文獻]

[1] 中國藥典.一部[S].2015：320.

[2] 洪浩，陳虎彪，徐風，等. 麻黃藥材原植物資源和市場品種調查[J]. 中國中藥雜志，2011，36（9）：1129.

[3] 中國科學院中國植物志編輯委員會.中國植物志.第7卷[M].北京：科學出版社，1997：474.

[4] 劉運東，齊妍婷，邱遠金，等. 麻黃屬的地理分布與起源演化[J]. 干旱區資源與環境，2009，23（6）：120.

[5] 滿多清，寥空太.中麻黃生境及栽培因子研究[J].甘肅農業大學學報，2003，38（1）：84.

[6] 張海龍.基于生態位模型的傳統中藥秦艽潛在地理分布研究[D]. 西安：陜西師范大學，2014.

[7] 郭蘭萍，黃璐琦，蔣有緒，等.基于3S的蒼術道地藥材研究（二）——蒼術揮發油組分的氣候主導因子篩選及氣候適宜性區劃[C].北京：全國中藥資源生態學學術研討會，2006.

[8] 張小波，郭蘭萍，黃璐琦.中藥區劃研究進展[J].中國農業與區劃，2010，31（3）：64.

[9] 張海娟，陳勇，黃烈健，等.基于生態位模型的薇甘菊在中國適生區的預測[J].農業工程學報，2011，27（1）：413.

[10] 何淑婷，白碧玉，但佳惠，等.基于MaxEnt的南丹參在中國的潛在分布區預測及適宜性分析[J].安徽農業科學，2014，42（8）：2311.

[11] 王茹琳，高曉清，王閆利，等.基于MaxEnt的非洲橘硬薊馬在全球及中國的潛在分布區預測[J].中國農學通報，2014，30（28）：315.

[12] 王運生，謝丙炎，萬方浩，等.ROC曲線分析在評價入侵物種分布模型中的應用[J].生物多樣性，2007，15（4）：365.

[13] 楊自輝，趙翠蓮.中麻黃的適應性及其抗逆性栽培[J].干旱區資源與環境，2003，17（1）：119.

[責任編輯 呂冬梅]