基于MaxEnt模型附子的全國藥材生態(tài)適應(yīng)性區(qū)劃研究

崔 顥，景志賢，周 華,3

(1. 澳門科技大學(xué)中醫(yī)藥學(xué)院，澳門特別行政區(qū)氹仔 999078；2. 中國中醫(yī)科學(xué)院中藥資源中心道地藥材國家重點(diǎn)實(shí)驗(yàn)室，北京 100700；3.中藥質(zhì)量研究國家重點(diǎn)實(shí)驗(yàn)室(澳門科技大學(xué))，澳門特別行政區(qū)氹仔 999078)

【研究意義】附子(Aconiti Lateralis Radix Praeparata,Fuzi)為毛茛科植物烏頭(AconitumcarmichaeliiDebeaux)子根的加工品，是我國傳統(tǒng)名貴中藥材，具有回陽救逆、補(bǔ)火助陽、逐寒除濕等功效[1]，其主要成分為烏頭堿、次烏頭堿和中烏頭堿等生物堿類物質(zhì)。四川江油是公認(rèn)的附子道地產(chǎn)區(qū)，有1300余年栽培歷史[2]；陜西漢中也有300余年栽培歷史，是附子的主產(chǎn)區(qū)之一[3]；云南目前也發(fā)展成為附子的主產(chǎn)區(qū)[4]，2020年栽培面積超過1667 hm2(數(shù)據(jù)來源云南省農(nóng)業(yè)廳)。但由于附子生態(tài)環(huán)境適宜性相關(guān)基礎(chǔ)性研究工作不足[5]，導(dǎo)致各產(chǎn)區(qū)附子藥材質(zhì)量參差不齊，給我國附子的規(guī)范化生產(chǎn)和質(zhì)量控制造成了較大的不利影響。因此，開展附子藥材的生態(tài)適宜性區(qū)劃研究，對提高附子規(guī)范化種植水平，保障附子藥材質(zhì)量具有重要意義。【前人研究進(jìn)展】附子產(chǎn)量及品質(zhì)不僅與種質(zhì)資源有關(guān)[6]，還受生態(tài)環(huán)境制約[7]。從采收期看，附子采收期也隨著氣候的變遷發(fā)生變化[8]。從產(chǎn)地來看，一般認(rèn)為四川江油與陜西漢中栽培的附子品質(zhì)最佳[9]。研究表明，生態(tài)因子對附子的生物產(chǎn)量及有效成分次烏頭堿含量均有影響，其中日照時(shí)數(shù)、降雨量及土壤pH等因素顯著影響附子的道地性[10]。 MaxEnt模型 (最大熵模型) 可以根據(jù)已知的物種分布點(diǎn)和氣候變量，對物種潛在的分布區(qū)進(jìn)行預(yù)測。該模型具有運(yùn)行時(shí)間短，操作簡單，精確度高等優(yōu)點(diǎn)[11]，因此得到廣泛應(yīng)用。MaxEnt模型被用于農(nóng)作物、經(jīng)濟(jì)作物和林業(yè)物種的潛在分布區(qū)預(yù)測，如未來氣候變化情景下中國區(qū)域水稻潛在適生區(qū)變化的預(yù)測[12]、經(jīng)濟(jì)作物獼猴桃潛在種植區(qū)的預(yù)測[13]、薄殼山核桃當(dāng)前和未來中國適宜種植區(qū)的預(yù)測[14]、小黃花茶在中國的適生區(qū)域研究[15]等。目前，利用中藥資源空間信息網(wǎng)絡(luò)數(shù)據(jù)庫，基于MaxEnt模型和ArcGIS軟件分析進(jìn)行中藥材生態(tài)適應(yīng)性區(qū)劃已經(jīng)成為研究熱點(diǎn)。張琴等[16]采用最大熵模型(MaxEnt)預(yù)測黃花蒿在中國的生態(tài)適宜區(qū)，應(yīng)用R語言構(gòu)建青蒿素含量與環(huán)境因子的關(guān)系模型，用ArcGIS軟件預(yù)測青蒿素空間分布以進(jìn)行黃花蒿生產(chǎn)區(qū)劃；嚴(yán)輝等[17]利用GIS技術(shù)和最大信息熵模型分析了當(dāng)歸適宜性分布；趙文龍等[18]采用同樣的方法結(jié)合全國范圍內(nèi)55個(gè)環(huán)境氣候因子的空間數(shù)據(jù)和中藥大黃3種基原植物的實(shí)際地理分布信息對中藥大黃3種不同基原植物的生境適宜性進(jìn)行評價(jià),預(yù)測了中藥大黃全國生境適宜性空間分布并進(jìn)行等級劃分；王浩等[19]采用該方法分析了野生防風(fēng)分布信息及多個(gè)環(huán)境因子，并進(jìn)行了防風(fēng)生態(tài)適宜性區(qū)劃，為栽培防風(fēng)的擴(kuò)大生產(chǎn)提供基礎(chǔ)研究數(shù)據(jù)。在附子適宜性研究方面，趙潤懷等[20]采用TCMGIS分析了我國附子道地藥材產(chǎn)地適應(yīng)性，發(fā)現(xiàn)土壤可能是影響附子分布的主導(dǎo)因素；魏玲等[15]對陜西漢中市附子種植氣候區(qū)劃進(jìn)行了研究，認(rèn)為氣候條件為漢中是否適合大規(guī)模發(fā)展商品附子的首要因素；肖小河等[21]首次運(yùn)用模糊數(shù)學(xué)法將四川附子的生態(tài)適宜區(qū)劃分為最適宜區(qū)、適宜區(qū)、較不適宜區(qū)、很不適宜區(qū)4類，并提出了烏頭及附子合理生產(chǎn)布局的初步設(shè)想。但未見MaxEnt 模型的附子生態(tài)適應(yīng)性分析的文獻(xiàn)報(bào)道?！颈狙芯壳腥朦c(diǎn)】本論文基于MaxEnt模型，對我國附子生態(tài)適應(yīng)性進(jìn)行研究?！緮M解決的關(guān)鍵問題】明確影響附子適宜性分布的主要生態(tài)因子，為我國附子道地產(chǎn)區(qū)的合理規(guī)劃及優(yōu)質(zhì)栽培提供科學(xué)依據(jù)。

1 材料與方法

1.1 樣品采集

2018—2019年于采收期7—9月在四川、陜西、云南等附子產(chǎn)區(qū)(圖1)采集附子樣品，記錄海拔、經(jīng)緯度、樣品數(shù)量(表1)。

表1 附子樣品采集信息

續(xù)表1 Continued table 1

1.2 分布位點(diǎn)數(shù)據(jù)收集

查詢中國數(shù)字植物標(biāo)本館(http://www.cvh.ac.cn)，記錄附子基原植物烏頭AconitumcarmichaeliiDebeaux項(xiàng)下所有標(biāo)本的采集地信息1294條，結(jié)合全國第四次中藥資源普查結(jié)果、文獻(xiàn)查閱、經(jīng)緯度定位、行政區(qū)劃變遷信息等，整理獲得以我國現(xiàn)行行政區(qū)劃下，以縣為單位的烏頭分布點(diǎn)287個(gè)(圖1)，作為生態(tài)區(qū)劃的基點(diǎn)。使用中國地圖對來自自然資源部門戶網(wǎng)站“標(biāo)準(zhǔn)地圖服務(wù)”(http://www.mnr.gov.cn/sj/sjfw/ch/)制圖。

1.3 環(huán)境及氣候數(shù)據(jù)

本研究采用全國各地1990—2019年常規(guī)氣象資料(包括平均氣溫、最高氣溫、最低氣溫、≥10 ℃積溫、相對濕度、風(fēng)速、日照時(shí)數(shù)等對附子分布影響較大的主要?dú)庀笠蜃覽21-23])，氣象數(shù)據(jù)來源于中國氣象數(shù)據(jù)網(wǎng)共享數(shù)據(jù)庫(http://data.cma.cn/)提供的中國地面氣候資料日值數(shù)據(jù)集(江油、漢中、布施、麗江4個(gè)地區(qū)氣象站)。其它生態(tài)因子數(shù)據(jù)來源于中藥資源空間信息網(wǎng)絡(luò)數(shù)據(jù)庫。

1.4 生態(tài)適應(yīng)性分析

1.4.1 模型設(shè)置 參考文檢等的分析方法[24]，將環(huán)境氣候數(shù)據(jù)和樣點(diǎn)坐標(biāo)數(shù)據(jù)導(dǎo)入MaxEnt模型，以85%的樣點(diǎn)數(shù)據(jù)作為訓(xùn)練樣本集，設(shè)置1×104的迭代次數(shù)用于估算各特征在訓(xùn)練數(shù)據(jù)中的分布，直到訓(xùn)練集樣本的特征分布和模型特征分布相同時(shí)給出模型最優(yōu)參數(shù)；剩余15%的樣點(diǎn)數(shù)據(jù)作為測試數(shù)據(jù)集，對預(yù)測結(jié)果進(jìn)行隨機(jī)測試，同時(shí)評價(jià)各環(huán)境氣候因子的貢獻(xiàn)率。MaxEnt模型預(yù)測精度以ROC曲線下的面積(AUC)為衡量指標(biāo)，其評價(jià)標(biāo)準(zhǔn)為：0.5～0.6為模型預(yù)測失敗，0.6～0.為7較差，0.7～0.8為一般，0.8～0.9為好，0.9～1.0為非常好[23]。

1.4.2 附子生態(tài)適應(yīng)性區(qū)劃 采用MaxEnt 模型預(yù)測不同環(huán)境氣候因子對附子生境適應(yīng)性，使用正態(tài)分布參數(shù)劃分法，通過ArcGIS 10.0軟件空間分析工具對MaxEnt預(yù)測圖進(jìn)行區(qū)劃，采用人工分級法將附子生態(tài)適宜性產(chǎn)區(qū)的適宜指數(shù)劃分4個(gè)區(qū)間。最適宜區(qū)：適宜指數(shù)[0.50，1.00]；適宜區(qū)：適宜指數(shù)[0.45，0.65)；次適宜區(qū)：適宜指數(shù)[0.25，0.45)；不適宜區(qū)：適宜指數(shù)[0，0.25)。適宜指數(shù)是根據(jù)最大熵原理由軟件直接計(jì)算獲得，其范圍為0～1.00。

2 結(jié)果與分析

2.1 MaxEnt模型預(yù)測精度及主導(dǎo)環(huán)境因子分析

對74個(gè)生態(tài)因子進(jìn)行4次運(yùn)算，除去貢獻(xiàn)率為0的生態(tài)因子，最終選擇26個(gè)生態(tài)因子進(jìn)行分析(表2)。附子樣品信息及生態(tài)因子數(shù)據(jù)通過MaxEnt模型分析后，訓(xùn)練數(shù)據(jù)AUC為0.995，測試數(shù)據(jù)AUC為0.985，說明MaxEnt模型對附子潛在適宜分布預(yù)測的準(zhǔn)確度和可信度非常高(圖2)。結(jié)合MaxEnt模型計(jì)算和Jackknife檢驗(yàn)進(jìn)行主導(dǎo)環(huán)境因子分析。其中，MaxEnt 模型計(jì)算得到附子各主要潛在環(huán)境因子貢獻(xiàn)率(表2)大于1的從大到小依次是：10月月降水量(35.6%)>海拔(24.9%)>濕潤指數(shù)(10.0%)>晝夜溫差與年溫差比值(5.6%)>土壤有效含水量等級(3.6%)>植被類型(2.2%)=11月月降水量(2.2%)>最干月降雨量(1.3%)>5月月降水量(1.1%)；Jackknife 檢驗(yàn)結(jié)果(圖3)顯示，各環(huán)境因子對附子的分布影響排序?yàn)椋?月月降水量>海拔>12月月降水量>4月月降水量=植被類型>7月月降水量>2月月降水量>5月、9月、10月月降水量>5月月均溫>3月月平均氣溫>11月月平均氣溫>11月月降水量=6月月平均氣溫。根據(jù)模型和刀切法檢驗(yàn)結(jié)果，選擇二者共有的環(huán)境因子，得出10月月降水量、海拔、植被類型、11月月降水量、5月月降水量，該5個(gè)環(huán)境因子是影響附子分布的主要環(huán)境因子。

表2 MaxEnt模型運(yùn)算的附子生態(tài)因子貢獻(xiàn)率和重要性

2.2 附子的生態(tài)適應(yīng)性分析及種植區(qū)劃

2.2.1 附子的生態(tài)適應(yīng)性分析 利用ArcGIS軟件對MaxEnt模型導(dǎo)出的生態(tài)適應(yīng)性進(jìn)行重分類，得到附子的生態(tài)適應(yīng)性區(qū)劃圖(圖4)。附子在我國的最大生態(tài)相似區(qū)域主要分布在云南、四川、西藏三省區(qū)結(jié)合部，也是目前附子的主要產(chǎn)區(qū)；適宜分布區(qū)包括貴州、重慶、陜西的秦嶺一帶，湖南、河南與陜西交界部分區(qū)域以及遼寧、吉林與黑龍江東部的少部地區(qū)；次適宜區(qū)包括西藏東南部，云南、貴州大部，四川西部及北部、重慶東部；其它地區(qū)為非宜分布區(qū)。分析結(jié)果與文獻(xiàn)記載的我國西南橫斷山區(qū)南部(四川西部、云南西北部和西藏東部)是國產(chǎn)烏頭屬植物的重要分布區(qū)高度吻合[16]，也與現(xiàn)階段附子的主要生產(chǎn)區(qū)分布一致。目前，我國附子種植面積約2667 hm2，其中，云南產(chǎn)區(qū)1333 hm2、四川產(chǎn)區(qū)667 hm2、陜西產(chǎn)區(qū)333 hm2，均在附子的最適宜和適宜分布區(qū)域。

2.2.2 附子產(chǎn)區(qū)種植區(qū)劃 我國附子種植區(qū)域的形成和發(fā)展與自然資源特點(diǎn)、生態(tài)環(huán)境因素、社會(huì)經(jīng)濟(jì)發(fā)展、科學(xué)技術(shù)進(jìn)步密切相關(guān)。本文在劃分我國附子產(chǎn)區(qū)種植區(qū)劃時(shí)，主要考慮以下幾個(gè)原則：一是附子藥材的道地性；二是附子的生態(tài)適應(yīng)性，產(chǎn)區(qū)的種植區(qū)應(yīng)集中在最適宜區(qū)和適宜區(qū)域；三是附子在各生態(tài)區(qū)域內(nèi)的種植比重及發(fā)展前景；四是目前全球氣候環(huán)境條件的變化與栽培技術(shù)的發(fā)展。根據(jù)上述原則，結(jié)合附子生態(tài)適應(yīng)性分布結(jié)果，將我國附子產(chǎn)區(qū)種植區(qū)劃劃分如下。

最適宜區(qū)：位于東經(jīng)100°～105°、北緯26°～32°區(qū)域。主要包括四川西部、云南西北部的大部分地區(qū)。該區(qū)日照充足、雨量充沛、年平均氣溫13～16 ℃，降水量900～1100 mm，無霜期220 d以上、海拔2000～3000 m，土壤為砂質(zhì)壤土。

適宜區(qū)：位于東經(jīng)98°～106°、北緯24°～33°除最適宜區(qū)外的區(qū)域。主要包括四川西部、云南西北部、西藏東部的部分地區(qū)及陜西漢中少部分區(qū)域。該區(qū)日照較為充足、雨量充沛、年平均氣溫12～17 ℃，降水量800～1200 mm，無霜期200 d以上、海拔1200～3500 m(漢中地區(qū)500 m左右)，土壤為砂質(zhì)壤土或紅壤。

次適宜區(qū)：位于東經(jīng)93°～110°、北緯23°～34.5°，除最適宜區(qū)、適宜區(qū)外的區(qū)域。包括四川西北部、云南北緯23°以北大部分區(qū)域、西藏東部、重慶東部、貴州中北部、陜西南部、河北、河南西部山區(qū)的部分地區(qū)、吉林、黑龍江東部部分區(qū)域。

不適宜區(qū)：附子基本不能生長的區(qū)域，我國除最適宜區(qū)、適宜區(qū)和次適宜外的其它地區(qū)。

附子生態(tài)適應(yīng)性理論分析表明，附子的適宜區(qū)主要分為我國西南和東北兩大區(qū)域，核心分布區(qū)位于除四川盆地外的西南部，理論適宜區(qū)與實(shí)際產(chǎn)區(qū)分布吻合。筆者通過產(chǎn)區(qū)考察，結(jié)合生態(tài)因子分析結(jié)果、文獻(xiàn)記載、自然條件、采收加工、栽培技術(shù)及文化傳承等社會(huì)經(jīng)濟(jì)條件，認(rèn)為附子的產(chǎn)區(qū)應(yīng)布局于最適宜區(qū)和適宜區(qū)，即以四川、云南、西藏三省區(qū)結(jié)合部和陜西漢中地區(qū)為最宜，包括四川江油、陜西漢中2個(gè)傳統(tǒng)產(chǎn)區(qū)和云南麗江、四川布施2個(gè)新產(chǎn)區(qū)。其它地區(qū)作為附子的潛在種植區(qū)域，在目前市場處于供需平衡的形勢下，不宜發(fā)展附子種植。

2.3 附子傳統(tǒng)產(chǎn)區(qū)與新產(chǎn)區(qū)產(chǎn)量及主要?dú)庀笠蜃颖容^分析

盡管附子在我國分布較廣，研究顯示我國適宜分布區(qū)包含四川、陜西、貴州、湖南、湖北、甘肅、云南、廣西、江西、安徽等10個(gè)省的336個(gè)縣市[22]，但我國公認(rèn)的附子道地產(chǎn)區(qū)僅有四川江油和陜西漢中，四川布施建國后也逐步成為附子產(chǎn)區(qū)。近幾年，云南逐步發(fā)展為全國最大的附子種植基地。2019年，云南麗江、大理、迪慶、昆明、曲靖等地附子種植面積已經(jīng)達(dá)到了1667 hm2以上，其中80%集中在麗江及與麗江毗鄰的賓川、鶴慶、劍川一帶，均處于相似的氣候環(huán)境條件，故本研究選擇麗江作為云南新產(chǎn)區(qū)、布施作為四川新產(chǎn)區(qū)與四川江油、陜西漢中等傳統(tǒng)道地產(chǎn)區(qū)的主要?dú)庀笠蜃舆M(jìn)行比較，為合理規(guī)劃附子種植提供更翔實(shí)的依據(jù)。從附子產(chǎn)量來看，新產(chǎn)區(qū)產(chǎn)量高于傳統(tǒng)產(chǎn)區(qū)(表3)。

附子在江油和漢中栽培的海拔高度較低(表4)，江油為450～600 m，漢中為470～1000 m，布施和麗江為2500～3000 m，屬于高海拔地區(qū)。年平均氣溫麗江最低、江油最高；幾個(gè)產(chǎn)區(qū)間降雨量沒有明顯差異。氣候環(huán)境條件的差異導(dǎo)致傳統(tǒng)產(chǎn)區(qū)與新產(chǎn)區(qū)采收期也有差異。江油產(chǎn)區(qū)一般在6—7月夏至前后必須采收，否則容易發(fā)生根腐病；漢中產(chǎn)區(qū)一般在7月至8月中下旬；布施8月底至9月初采收，麗江采收期則晚到10月底至11月初左右采收。年平均氣溫每降低1 ℃，采收期推遲1個(gè)月左右。從植物生長的角度來看，氣溫較低，有利于干物質(zhì)積累。

表3 不同產(chǎn)區(qū)附子產(chǎn)量比較

從日照條件來看，新產(chǎn)區(qū)日照時(shí)數(shù)高于傳統(tǒng)產(chǎn)區(qū)。附子屬于喜光植物，長日照有利附子產(chǎn)量形成。傳統(tǒng)產(chǎn)區(qū)四川江油日照時(shí)數(shù)比新產(chǎn)區(qū)云南麗江少1000 h以上。盡管新產(chǎn)區(qū)出苗期比傳統(tǒng)產(chǎn)區(qū)晚1～2個(gè)月，但由于采收期延遲，新產(chǎn)區(qū)附子生長期長于傳統(tǒng)產(chǎn)區(qū)。從光合作用、干物質(zhì)積累等附子生長條件來看，新產(chǎn)區(qū)屬于高海拔溫涼地區(qū)，環(huán)境條件更利于附子生長，也利于次生代謝產(chǎn)物的積累。傳統(tǒng)產(chǎn)區(qū)往往日照不足，這可能是導(dǎo)致傳統(tǒng)產(chǎn)區(qū)附子病害較多的原因之一。黨玨等[23]研究表明，江油的氣候條件最符合附子的生性，布拖和漢中的生態(tài)環(huán)境雖不及江油，亦在附子生長的適宜范圍之內(nèi)。

3 討 論

附子的栽培雖然有上千年歷史[2]，但對附子栽培的相關(guān)研究還不多，對生態(tài)適應(yīng)性及生產(chǎn)區(qū)劃研究更少。藥材的形成與生態(tài)型之間存在著密切的內(nèi)在聯(lián)系。同種植物由于生態(tài)型的不同，其植物個(gè)體在形態(tài)、生理、生化上均表現(xiàn)出差異，這些差異則是產(chǎn)生藥材產(chǎn)量高低、品質(zhì)優(yōu)劣和療效差異的實(shí)質(zhì)原因。無論從文獻(xiàn)記載還是栽培歷史來看，四川江油是公認(rèn)的附子道地產(chǎn)區(qū)，陜西漢中栽培歷史較江油晚，到現(xiàn)在約300年，質(zhì)量也不及江油[21]。本研究氣象因子分析表明，江油和漢中核心氣象要素比較接近，產(chǎn)量相差不大，且栽培、產(chǎn)地加工工藝是一脈相承，漢中也是附子道地產(chǎn)區(qū)。

表4 不同附子產(chǎn)區(qū)氣象要素比較

四川布施、云南麗江和大理等是近二十年才發(fā)展起來的附子新產(chǎn)區(qū)[22]。新產(chǎn)區(qū)的氣象特點(diǎn)為緯度較低、海拔較高，日照時(shí)數(shù)長，氣溫較低的特點(diǎn)，屬于高山溫涼氣候類型。從日照條件來看，新產(chǎn)區(qū)日照時(shí)數(shù)高于傳統(tǒng)產(chǎn)區(qū)。附子屬于喜光植物，長日照有利附子產(chǎn)量形成。傳統(tǒng)產(chǎn)區(qū)四川江油日照時(shí)數(shù)比新產(chǎn)區(qū)云南麗江少1000 h以上。盡管新產(chǎn)區(qū)出苗期比傳統(tǒng)產(chǎn)區(qū)晚1～2個(gè)月，但由于采收期延遲，新產(chǎn)區(qū)附子生長期長于傳統(tǒng)產(chǎn)區(qū)。從光合作用、干物質(zhì)積累等附子生長條件來看，新產(chǎn)區(qū)屬于高海拔溫涼地區(qū)，環(huán)境條件更利于附子生長，也利于次生代謝產(chǎn)物的積累。傳統(tǒng)產(chǎn)區(qū)往往日照不足，這可能是導(dǎo)致傳統(tǒng)產(chǎn)區(qū)附子病害較多的原因之一。黨玨等[23]研究表明，江油的氣候條件最符合附子的生性，布拖和漢中的生態(tài)環(huán)境雖不及江油，亦在附子生長的適宜范圍之內(nèi)。調(diào)查結(jié)果表明，新產(chǎn)區(qū)附子的生長期較長，產(chǎn)量也較高。從栽培技術(shù)方面看，傳統(tǒng)產(chǎn)區(qū)栽培較為規(guī)范，附子的等級規(guī)格高，產(chǎn)值也較高；新產(chǎn)區(qū)栽培比較粗放，大多不采用修根、摘尖等栽培技術(shù)，雖然總產(chǎn)量較傳統(tǒng)產(chǎn)區(qū)高，規(guī)格等級較低[25-26]。這也是普遍認(rèn)為新產(chǎn)區(qū)附子質(zhì)量不如傳統(tǒng)產(chǎn)區(qū)的主要原因。從生態(tài)角度看，新產(chǎn)區(qū)氣候環(huán)境更有利于附子的干物質(zhì)積累和次生代謝產(chǎn)物的形成，羅春梅等[27]對陜西、四川江油、布施和云南4個(gè)產(chǎn)區(qū)的化學(xué)成分分析表明，4個(gè)產(chǎn)區(qū)單酯型生物堿含量均很低，僅陜西附子、部分江油和云南附子檢測到苯甲酰新烏頭原堿，4個(gè)產(chǎn)區(qū)附子均未測出苯甲酰烏頭原堿和苯甲酰次烏頭原堿。3種雙酯型生物堿總量以云南附子最高，其次是布施、江油和陜西。模型的理論分析結(jié)果與本研究采用傳統(tǒng)的氣象因素分析結(jié)果基本一致。影響附子分布的5個(gè)環(huán)境因子中，降水量占了3個(gè)。當(dāng)然，影響附子產(chǎn)量和品質(zhì)還有土壤環(huán)境、栽培技術(shù)及產(chǎn)地加工技術(shù)等綜合因素[28]。趙潤懷等[20]認(rèn)為，附子的適宜產(chǎn)區(qū)主要是以四川盆地周圍、陜西南部、云貴高原和湘鄂地區(qū)為中心，向南最遠(yuǎn)能夠延伸到廣西北部，向東擴(kuò)展到江西和安徽西部地區(qū)，向北到甘肅南部地區(qū)。筆者分析顯示附子生態(tài)適宜區(qū)分布范圍較小，也與當(dāng)前的生產(chǎn)實(shí)際更為接近。

4 結(jié) 論

根據(jù)MaxEnt 模型模型和ArcGIS軟件分析結(jié)果，結(jié)合環(huán)境因素對附子生態(tài)適應(yīng)性的貢獻(xiàn)率和重要性研究表明，10月月降水量、海拔、植被類型、11月月降水量、5月月降水量這5個(gè)環(huán)境因子是影響附子分布的主要環(huán)境因子。結(jié)合附子生產(chǎn)種植實(shí)際，認(rèn)為附子的最適生態(tài)區(qū)為四川、云南、西藏三省區(qū)結(jié)合部以及陜西南部區(qū)域。建議附子產(chǎn)區(qū)規(guī)劃應(yīng)以四川、云南、西藏三省區(qū)結(jié)合部及陜西漢中部分區(qū)域?yàn)榈赖禺a(chǎn)區(qū)，其它地區(qū)為一般種植區(qū)域，以確保附子藥材的道地性，促進(jìn)附子產(chǎn)業(yè)的高質(zhì)量發(fā)展。