巫明焱　何蘭　陳佳麗　董光　程武學(xué)

[摘要]研究川牛膝在全國范圍的潛在分布及種植空缺，有利于為其實(shí)地保護(hù)及種植推廣工作提供依據(jù)。該文結(jié)合63個(gè)川牛膝地理分布數(shù)據(jù)和49個(gè)相關(guān)生態(tài)因子，利用3S技術(shù)及MaxEnt模型對川牛膝在全國的潛在分布進(jìn)行預(yù)測及種植空缺分析。研究結(jié)果表明：①川牛膝在全國潛在分布面積約634 38580 km2，主要分布于四川省的東北部、東南部，云南省的北部、東南部，貴州省的西部、西南部，重慶市的西南部、東北部，陜西省的西南部，甘肅省的東南部，湖北省的西部，廣西省的西部，西藏自治區(qū)的東南部等區(qū)域；②影響川牛膝生長的主要生態(tài)因子是降水量、海拔、最冷月極低溫度、土壤類型、月均溫等；③全國川牛膝種植空缺區(qū)域主要位于四川省廣元市、綿陽市、雅安市、樂山市、涼山州、攀枝花市，陜西省漢中市，云南省大理自治州、怒江自治州、楚雄自治州、保山市、曲靖市、文山自治州，貴州省黔西南自治州等。該研究結(jié)果對掌握川牛膝的生長環(huán)境、預(yù)測在全國的潛在分布及其栽培推廣具有重要意義。

[關(guān)鍵詞]川牛膝； MaxEnt模型； 生態(tài)因子； 分布； 適宜性； 種植空缺

[Abstract]Research on predictive distribution and planting GAP of Cyathula officinalis in China is helpful to provide scientific basis for its protection and planting popularization According to the data in 63 distribution sites and 49 ecological variables， using MaxEnt ecological niche model and 3S technology， we performed a quantitative analysis of suitable distribution and planting GAP of C officinalis in China Our results show that： ① the area of suitable distribution of C officinalis is about 634 38580 km2 in total， and mainly in Northeastern and Southeastern Sichuan， Northern and Southeastern Yunnan， Western and Southwestern Guizhou， Southwestern and Northeastern Chongqing， Southwestern Shaanxi， Southeastern Gansu， Western Guangxi， Southeastern Tibet ② The main ecological factors determining the potential distribution are precipitation， altitude， minimum temperature of coldest month， soil type， monthly mean temperature ③ The planting GAP region are mainly in Guangyuan， Mianyang， Ya′an， Leshan， Liangshan， Panzhihua of Sichuan province， Hanzhong of Shaanxi province， Dali， Nujiang， Chuxiong， Baoshan， Qujing， Wenshan of Yunnan province， southwestern autonomous prefecture in Guizhou province The results are of great significance for realizing the growth environment， predicting the potential distribution and promoting planting popularization for C officinalis

[Key words]Cyathula officinalis； MaxEnt modeling； ecological factors； distribution； suitability； planting GAP

川牛膝又稱白牛膝、拐牛膝，是莧科植物川牛膝Cyathula officinalis Kuan的干燥根[1]。其味甘、性平、微苦，主要化學(xué)成分為植物留酮、多糖及阿魏酸等，具有逐瘀通經(jīng)，通利關(guān)節(jié)，利尿通淋的效用[2]。川牛膝是著名的川產(chǎn)道地藥材之一[3]，主要分布在四川、貴州、云南、湖北、陜西、湖南等地。近年來隨著對中藥材需求的增加，川牛膝的人工栽培范圍也在逐步擴(kuò)大，但針對其適種條件、適種區(qū)域、種植管理方式等還沒有較深入的評價(jià)指標(biāo)來權(quán)衡種植的適宜性，而且目前人工栽培的方式較粗放，川牛膝品種混雜混交，故市場上川牛膝的品質(zhì)也參差不齊。因此，采用科學(xué)的方法，準(zhǔn)確合理的掌握影響川牛膝生長及栽培的生態(tài)因子，并開展對川牛膝潛在分布區(qū)的預(yù)測與評價(jià)工作，對掌握川牛膝的實(shí)地生長狀況、科學(xué)推廣人工栽培、提升川牛膝質(zhì)量水平具有重要意義。

近年來，隨著3S技術(shù)的不斷發(fā)展，利用遙感數(shù)據(jù)快速準(zhǔn)確的分析中藥材的適生條件、適生區(qū)域越來越廣泛，特別是3S技術(shù)在物種保護(hù)、物種適宜栽培區(qū)域及物種儲(chǔ)量估算等方面發(fā)揮著重要的作用[46]。生態(tài)位模型是利用物種存在數(shù)據(jù)和影響因子數(shù)據(jù)，按照特定的方法估算物種的生態(tài)位，并投影到待預(yù)測區(qū)域中，以概率的形式呈現(xiàn)物種生境特征[7]。目前將MaxEnt生態(tài)位模型與3S技術(shù)結(jié)合應(yīng)用于中藥材的潛在分布預(yù)測及評價(jià)已較為廣泛，為中藥材的栽培管理提供了準(zhǔn)確的數(shù)據(jù)支撐。張東方等[8]利用最大信息熵模型和GIS技術(shù)對當(dāng)歸進(jìn)行全球適宜性區(qū)劃，為當(dāng)歸的人工引種栽培及精細(xì)化種植管理提供了一定的參考；郭杰等[9]利用MaxEnt模型預(yù)測黨參的全球生態(tài)適宜區(qū)，并得出影響其地理分布的生態(tài)因子主要是溫度因子、降水因子；唐策[10]等利用MaxEnt和GIS技術(shù)對翼首草進(jìn)行適宜性分析，為其野生撫育及GAP種植選址提供了科學(xué)的依據(jù)。endprint

縱觀針對川牛膝的現(xiàn)有研究，大多數(shù)都是資源調(diào)查概述、品質(zhì)評價(jià)、藥材化學(xué)研究及部分區(qū)域的適宜性分布研究[1114]，對全國范圍川牛膝適生區(qū)的劃分較少，對影響川牛膝栽培條件的定量研究鮮見，對全國川牛膝適宜栽培地的空缺分析則未曾有研究。本文擬結(jié)合川牛膝實(shí)際分布數(shù)據(jù)，利用3S技術(shù)及MaxEnt生態(tài)位模型，對川牛膝在全國范圍的潛在分布進(jìn)行預(yù)測，探索影響川牛膝生長的生態(tài)因子，并對全國范圍內(nèi)川牛膝適宜栽培區(qū)的空缺分析，以期為全國范圍川牛膝的物種保護(hù)及人工栽培推廣提供理論支持及數(shù)據(jù)支撐。

1材料與方法

11分布數(shù)據(jù)

川牛膝分布數(shù)據(jù)主要來源于全球生物多樣性信息網(wǎng)絡(luò)（GBIF，http：//www.gbif.org/）、中國數(shù)字植物標(biāo)本館（CVH，http：//wwwcvhaccn/）、文獻(xiàn)檢索及野外手持GPS測量獲取。剔除經(jīng)緯度重復(fù)及實(shí)際不符合的存在點(diǎn)，最終得到63個(gè)存在點(diǎn)信息，并轉(zhuǎn)換為模型支持的CSV格式文件。其中模型運(yùn)算的47個(gè)川牛膝分布點(diǎn)數(shù)據(jù)見表1。

12環(huán)境數(shù)據(jù)

綜合考慮川牛膝栽培及適生地的環(huán)境條件[15]，本文共選取了生物氣候因子、氣溫降水因子、地形因子、土壤因子、土地利用類型因子及植被覆蓋度度因子共6類49個(gè)生態(tài)因子用于研究，見表2。生物氣候因子及氣溫降水因子源于世界氣候數(shù)據(jù)庫（http：//wwwworldclimorg/），包含19個(gè)生物年均氣候因子、12個(gè)月降水量因子及12個(gè)月均溫因子，該氣候環(huán)境數(shù)據(jù)是利用克里金插值法將全球1950—2000年全球有記錄的氣象站點(diǎn)觀測數(shù)據(jù)生成的柵格數(shù)據(jù)；地形因子包含高程、坡度、坡向共3個(gè)因子，高程數(shù)據(jù)源于USGS網(wǎng)站（http：//glovis.usgs.gov/），其中坡度、坡向由DEM數(shù)據(jù)衍生獲取；土壤因子源于世界土壤數(shù)據(jù)庫（http：//webarchive.iiasa.ac.at/），含1個(gè)土壤類型因子；土地利用類型因子及植被覆蓋度度因子源于世界地圖數(shù)據(jù)庫（https：//globalmaps.github.io/），含1個(gè)土地利用類型因子和1個(gè)植被覆蓋度因子。

13模型與方法

生態(tài)位模型是利用獲取的物種分布數(shù)據(jù)和相關(guān)影響因子，依據(jù)一定算法來推算影響物種適宜分布的因素，再以算出的數(shù)據(jù)投影到其他的時(shí)空來預(yù)測物種潛在適生區(qū)域[1617]。本文采用3S技術(shù)結(jié)合生態(tài)位模型MaxEnt （http：//wwwcsprincetonedu/～schapire/MaxEnt）對川牛膝在中國的適宜生境及影響其生長的主要生態(tài)因子進(jìn)行探究。首先，利用ENVI 54對所有柵格數(shù)據(jù)進(jìn)行預(yù)處理，并裁剪出全國范圍的柵格數(shù)據(jù)，再利用ArcGIS 102提取研究的所有生態(tài)因子，并統(tǒng)一坐標(biāo)系和統(tǒng)一像元大小，并轉(zhuǎn)換為模型所支持ASCII文件格式；然后，隨機(jī)選取75%的川牛膝分布點(diǎn)用于構(gòu)建模型，剩余的25%用于模型的驗(yàn)證與測試，重復(fù)迭代運(yùn)算500次，利用MaxEnt模型初步運(yùn)算川牛膝的潛在分布；最后篩選主要貢獻(xiàn)率的生態(tài)因子進(jìn)行二次MaxEnt模型運(yùn)算，再利用SPSS 200 評價(jià)結(jié)果的精確性和穩(wěn)定性，分析影響川牛膝生長的主要生態(tài)因子，并使用ArcGIS 102生成川牛膝潛在分布區(qū)及栽培種植空缺區(qū)域。本文研究的技術(shù)路線見圖1。

2結(jié)果與分析

21模型精度評價(jià)

ROC曲線是根據(jù)一系列不同的二分類方式，以真陽性率（靈敏度）為縱坐標(biāo)，假陽性率（1特異度）為橫坐標(biāo)繪制的曲線[18]。本文利用ROC曲線與橫坐標(biāo)圍成的面積來表征預(yù)測結(jié)果的準(zhǔn)確性，值范圍為[0，1]，其值越高，預(yù)測結(jié)果越準(zhǔn)確。本文研究結(jié)果中運(yùn)算模型數(shù)據(jù)的AUC達(dá)0965，驗(yàn)證數(shù)據(jù)的AUC達(dá)0946，均達(dá)較高的預(yù)測精度，表明運(yùn)算結(jié)果的準(zhǔn)確性較高，可靠性較強(qiáng)，見圖2。

22主要生態(tài)因子貢獻(xiàn)率

生態(tài)因子是影響川牛膝生長的主要因素，不同因子對川牛膝適生區(qū)及栽培的影響程度不同。通過刀切法（Jackknife）定量研究了各生態(tài)因子對川牛膝適生區(qū)的貢獻(xiàn)率，結(jié)果表明10月降水量對川牛膝生長的貢獻(xiàn)率最大，達(dá)到398%，其他對川牛膝生長具有影響的生態(tài)因子主要是海拔、最冷月極低溫度、土壤類型、5，7，8，9，11月降水量、土地利用類型、4，8月平均溫度、植被覆蓋度、坡度等，累積貢獻(xiàn)率達(dá)984%。由此可知，充足的降雨量是川牛膝生長的必要條件，高程值限定了川牛膝生長的區(qū)域，土壤影響了川牛膝生長的品質(zhì)，月均溫、坡度一定程度的影響川牛膝的生長環(huán)境。各主要生態(tài)因子的貢獻(xiàn)率見表3。

23川牛膝適生條件定量分析

231氣溫降水因子分析降水量是影響川牛膝生長的最主要因子，氣溫對其生長環(huán)境能造成一定的影響。川牛膝人工栽培一般在3—4月，采收期在11月，由研究可知川牛膝生長周期內(nèi)需要充足的降雨量，特別是采收期前1個(gè)月（10月）對降水量的需求較大。研究表明川牛膝生長全周期適宜的月降水量為100～800 mm，但月降雨量在100～200 mm最合適，見圖3。適量并充足的水分能促進(jìn)其根系的生長發(fā)育，但不同月份的降水量不能超過其承受的極限；氣溫對川牛膝生長的影響主要是生長初期和中期，研究表明月平均溫度為15～23 ℃是川牛膝生長的最適宜條件，特別是研究結(jié)果中4月和8月平均溫度對川牛膝的生長影響較明顯，見圖3。在人工栽培時(shí)應(yīng)合理控制上述研究中川牛膝生長環(huán)境土壤水分和溫度以提供良好的生長外界環(huán)境。

232地形因子分析海拔是影響川牛膝生長的次重要生態(tài)因子，其生長區(qū)的海拔具有一定的限制，坡度也會(huì)輕微的影響川牛膝的分布見圖3。研究表明川牛膝主要生長在1 000～2 800 m的高程，其中1 500～2 500 m是其生長的最佳高度范圍；坡度也能一定程度影響川牛膝的生長，研究表明平緩坡更有利于川牛膝的生長，坡度在35°的范圍內(nèi)隨坡度的增大，適宜性逐漸增強(qiáng)。坡度在35°范圍以上的區(qū)域川牛膝分布的可能性相對較小。endprint

233土壤因子分析土壤能為川牛膝的生長提供營養(yǎng)，能影響植株的整體發(fā)育。研究表明，川牛膝主要生長在土層深厚、富含有機(jī)質(zhì)、排水性強(qiáng)的山地棕壤、黃壤和紫色土，見圖3。所以在進(jìn)行人工栽培時(shí)應(yīng)當(dāng)充分考慮栽培區(qū)的土壤類型以提供良好的土壤生長環(huán)境。

234生物氣候因子分析生物氣候因子的定量數(shù)據(jù)能揭示適宜川牛膝生長的生物氣候變量幅度范圍[19]。研究結(jié)果表明所有生物氣候因子中最冷月極低溫度對川牛膝的生長影響最大，-1～5 ℃是最適宜川牛膝生長的最冷月極低溫度范圍，見圖3。

此外，除上述因子對川牛膝的生長有較大影響外，研究表明，土地利用類型、植被覆蓋度對川牛膝的生長也存在輕微影響。川牛膝主要生長在林地邊緣及草叢地中，且因其喜陰涼，故生長地的植被覆蓋度為20%～60%最佳。

24川牛膝潛在分布區(qū)劃

將模型的運(yùn)行結(jié)果導(dǎo)入到ARCGIS 102中，依據(jù)采樣點(diǎn)適生指數(shù)P的平均值、標(biāo)準(zhǔn)差及專家經(jīng)驗(yàn)法將其潛在分布區(qū)劃分為3個(gè)等級：P<0349 27不適宜區(qū)，0349 27≤P≤0652 35適宜區(qū)，P>0652 35最適宜區(qū)。并以此生成川牛膝在全國的潛在分布范圍，見圖4。

從分布區(qū)域來看，川牛膝在全國適宜生長的潛在區(qū)域主要集中在中國西南部，即包含四川省的東北部、東南部，云南省的北部、東南部，貴州省的西部、西南部，重慶市的西南部、東北部，陜西省的西南部，甘肅省的東南部，湖北省的西部，廣西省的西部，西藏自治區(qū)的東南部等區(qū)域。

從最適宜分布區(qū)域來看，川牛膝生長的最適宜區(qū)域主要位于四川省東北—西南走向的沿線區(qū)域，云南省北部的東西走向沿線區(qū)域，貴州省與云南省接壤的南北走向沿線區(qū)域。

從分布區(qū)面積，川牛膝在全國潛在適宜生長區(qū)域的面積為634 38580 km2，最適宜區(qū)生長區(qū)面積為128 71491 km2。全國范圍川牛膝潛在分布面積最大的是云南省，適宜生長面積達(dá)236 63197 km2；其次是四川省，適宜生長面積達(dá)184 31699 km2；另貴州省、陜西省、西藏自治區(qū)、廣西壯族自治區(qū)、甘肅省、湖北省適生面積依次減少，見表4。

25川牛膝適宜栽培區(qū)空缺分析

依據(jù)目前已有調(diào)查研究[2021]及查閱相關(guān)資料，本文掌握了川牛膝目前在全國范圍內(nèi)的主要栽培種植區(qū)有四川省的天全縣、寶興縣、滎經(jīng)縣、漢源縣、名山縣、樂山市金口河區(qū)、鹽邊縣、會(huì)東縣、安縣；重慶市的奉節(jié)縣、巫溪縣；湖北省的恩施市；湖南省的龍山縣、隆回縣；云南省的騰沖縣；貴州省的畢節(jié)市等區(qū)域。將MaxEnt模型預(yù)測的川牛膝在全國最適宜生長區(qū)與全國已有人工栽種培育區(qū)域進(jìn)行對比分析，利用ArcGIS 102進(jìn)一步對川牛膝適宜種植區(qū)進(jìn)行空缺分析。研究結(jié)果表明川牛膝還有較大面積的適宜栽培的空缺區(qū)域，主要分布在四川省廣元市的西北部、綿陽的東北部、雅安市的東北部和西南部、樂山市的西南部、涼山州的中部、攀枝花市的東部；陜西省漢中市的西南部；云南省大理自治州、怒江自治州、楚雄自治州、保山市的東部、曲靖市的南部、文山自治州；貴州省的黔西南自治州西部。較典型的縣域有寧強(qiáng)縣、青川縣、平武縣、大邑縣、蘆山縣、石棉縣、馬邊縣、峨邊縣、賓川縣、呂寧縣、永平縣、南華縣、羅平縣、廣南縣、硯山縣、興義市等區(qū)域。因此可依據(jù)川牛膝適宜栽培的空缺區(qū)域，有依據(jù)的合理擴(kuò)大川牛膝的種植范圍，但還需考慮當(dāng)?shù)貙?shí)際種植環(huán)境，先進(jìn)行小面積范圍的試栽培，再依據(jù)試驗(yàn)結(jié)果進(jìn)行推廣種植，見圖5。

3討論

川牛膝是川產(chǎn)道地常用中藥材，全國每年對該藥材的需求量巨大，定量分析影響其生長的生態(tài)因子、研究其在全國范圍的潛在分布并進(jìn)行種植空缺分析對掌握川牛膝的適宜生境、促進(jìn)川牛膝種植空缺區(qū)域的推廣種植、提升藥材品質(zhì)具有重要意義。

31川牛膝潛在分布結(jié)果分析

本文利用MaxEnt模型及3S技術(shù)研究了川牛膝在全國范圍的潛在分布及種植空缺，利用ROC曲線對研究結(jié)果精度評價(jià)，其AUC達(dá)0965，精度較高，研究結(jié)果較準(zhǔn)確。但受采樣區(qū)川牛膝與其它品種混交、樣點(diǎn)范圍具有一定的區(qū)域局限性等影響，研究結(jié)果的精度會(huì)受到影響。

通過研究，降水量是影響川牛膝生長的關(guān)鍵因素，說明川牛膝生長周期內(nèi)對水分的需求相對較大，但是水分必須保證在合適的范圍內(nèi)，才能有利于其根部的生長發(fā)育。其次是海拔限制了川牛膝潛在分布的區(qū)域，最佳的海拔高度為1 500～2 500 m。另外最冷月極低溫度、土壤類型等對川牛膝的生長也具有一定程度的影響。本文利用模型對生態(tài)因子對川牛膝生長貢獻(xiàn)率的研究結(jié)果與利用化學(xué)成分與川牛膝產(chǎn)區(qū)生態(tài)環(huán)境因素的灰色關(guān)聯(lián)度[22]的結(jié)果一致，表明本文對影響川牛膝生長的生態(tài)因子進(jìn)行定量分析的結(jié)果較準(zhǔn)確。

本文中分析了川牛膝在全國范圍的潛在分布，主要集中在四川省、云南省、貴州省、重慶市、陜西省、甘肅省、湖北省、廣西省、西藏自治區(qū)等地。該研究結(jié)果與前人對川牛膝分布區(qū)的調(diào)查結(jié)果[21]相吻合。但是模型預(yù)測結(jié)果中云南省的面積相對較大，略微超過了四川省的潛在分布面積，這與參與運(yùn)算的川牛膝樣點(diǎn)數(shù)據(jù)有一定關(guān)系，但從研究結(jié)果來看，云南省降水充足、土壤肥沃的條件的確有利于川牛膝的生長。另研究中西藏自治區(qū)中低海拔的東南區(qū)域也存在川牛膝潛在分布的概率，依據(jù)中國數(shù)字植物標(biāo)本館記錄，西藏自治區(qū)的確有少量野生川牛膝標(biāo)本采集的記錄。

32川牛膝種植空缺及栽培推廣

川牛膝種植空缺的分析對其種植推廣具有重要意義。本文基于已有川牛膝栽培區(qū)域的調(diào)查數(shù)據(jù)進(jìn)行了全國川牛膝的種植空缺分析。該結(jié)果是依據(jù)川牛膝最佳適宜生長區(qū)域而進(jìn)行的分析，忽略了模型預(yù)測結(jié)果中川牛膝適宜度次于最適宜分布的區(qū)域，但該結(jié)果能更準(zhǔn)確的為全國范圍內(nèi)川牛膝的栽培推廣提供依據(jù)。

川牛膝的品質(zhì)與其生長的環(huán)境有重要的聯(lián)系。在進(jìn)行種植推廣時(shí)，為保證川牛膝的品質(zhì)，務(wù)必依據(jù)其生長條件及環(huán)境，精細(xì)化的管理栽培。首先，保證其生長周期內(nèi)水分的足而不過；其次，溫度的控制，在其生長初期和中期需保證15～23 ℃的最佳適宜溫度范圍；再次，栽培育種地的選擇，需考慮1 000～2 800 m合適的高程及生長地的土壤需考慮排水性強(qiáng)、土層深厚、富含有機(jī)質(zhì)的棕壤型土；還有，依據(jù)最冷月極度溫度對川牛膝生長影響，過低的溫度會(huì)直接凍壞該物種，故應(yīng)隨時(shí)關(guān)注其生長地夜晚的最低溫度，但研究表明一定的低溫有利于川牛膝生長發(fā)育，這與其葉片夜晚呼吸作用的能耗有關(guān)，故在采收前期可適度的增大晝夜溫差，減少夜晚呼吸能量的損耗，有利于其根部能量的儲(chǔ)存，促進(jìn)其產(chǎn)量的增收；最后，任何植物在進(jìn)行栽培推廣前，必須結(jié)合推廣地生長環(huán)境，先小范圍的試種，再依據(jù)種植成效進(jìn)行種植推廣。endprint

[參考文獻(xiàn)]

[1]吳向莉 川牛膝與麻牛膝的鑒別[J] 貴陽中醫(yī)學(xué)院學(xué)報(bào)，2000，22（1）：62

[2]趙華杰，舒光明，周先健，等 我國川牛膝資源分布及生產(chǎn)狀況調(diào)查[J] 資源開發(fā)與市場，2012，28（5）：414

[3]葉品良，彭娟，劉娟 川牛膝研究概況[J] 中醫(yī)藥學(xué)報(bào)，2007（2）：51

[4]劉金欣，潘敏，李耿，等 3S技術(shù)在藥用植物資源調(diào)查研究中的應(yīng)用[J] 中草藥，2016，47（4）：695

[5]劉金欣，李耿，夏艷華，等 基于3S技術(shù)內(nèi)蒙古地區(qū)野生黃芩儲(chǔ)量估算研究[J] 中草藥，2016，47（6）：997

[6]郭蘭萍，黃璐琦，蔣有緒 “3S”技術(shù)在中藥資源可持續(xù)利用中的應(yīng)用[J] 中國中藥雜志，2005，30（18）：1397

[7]Guisa A，Zimmermann N E Predictive habitat distribution models in ecology[J] Ecol Model，2000，135（2）：147

[8]張東方，張琴，郭杰，等 基于MaxEnt模型的當(dāng)歸全球生態(tài)適宜區(qū)和生態(tài)特征研究[J] 生態(tài)學(xué)報(bào)，2017，37（15）：1

[9]郭杰，劉小平，張琴，等 基于Maxent模型的黨參全球潛在分布區(qū)預(yù)測[J] 應(yīng)用生態(tài)學(xué)報(bào)，2017，28（3）：992

[10]唐策，文檢，張雯，等 基于最大熵模型的藏藥翼首草生態(tài)適宜性研究[J] 中國中藥雜志，2017，42（10）：1871

[11]范巧佳，劉帆，高雨榮，等 天全川牛膝生產(chǎn)現(xiàn)狀調(diào)查[J] 中藥研究與信息，2001，3（11）：16

[12]童凱 川牛膝品質(zhì)評價(jià)研究[D] 雅安：四川農(nóng)業(yè)大學(xué)，2013

[13]孫傳鑫，郭晶，王秋紅，等 川牛膝化學(xué)拆分組分的制備及其免疫調(diào)節(jié)作用[J] 中國實(shí)驗(yàn)方劑學(xué)雜志，2016，22（2）：70

[14]尚雪，董麗君，文路軍，等 基于遙感與GIS的四川省川牛膝資源適宜性分布研究[J] 中草藥，2016，47（24）：4445

[15]羅緒和 川牛膝栽培方法有講究[N] 中國中醫(yī)藥報(bào)，20060316（007）

[16]朱耿平，劉國卿，卜文俊，等 生態(tài)位模型的基本原理及其在生物多樣性保護(hù)中的應(yīng)用[J] 生物多樣性，2013，21（1）：90

[17]Phillips S J，Anderson R P，Schapire R E Maximum entropy modeling of species geographic distributions[J] Ecol Model，2006，190（3/4）：231

[18]王漢卿，馬玲，王慶，等 甘草藥材生產(chǎn)區(qū)劃研究[J] 中國中藥雜志，2016，41（17）：3122

[19]Waltari E，Schroeder R，Mcdonald K，et al Bioclimatic variables derived from remote sensing： assessment and application for species distribution modelling[J] MEE，2014，5（10）：1033

[20]趙華杰 川牛膝資源調(diào)查及品質(zhì)評價(jià)研究[D] 瀘州：瀘州醫(yī)學(xué)院，2012

[21]劉維，裴瑾，楊梅，等 川產(chǎn)道地藥材川牛膝產(chǎn)地變遷探討[J] 中草藥，2016，47（9）：1625

[22]劉維，張祎楠，裴瑾，等 川牛膝品種與品質(zhì)的灰色關(guān)聯(lián)度分析研究[J] 中國藥學(xué)雜志，2014，49（20）：1796

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