陳博，曹后康，呂蓉，王莉，成希，晉玲，2

1.甘肅中醫(yī)藥大學(xué)藥學(xué)院，甘肅 蘭州 730000；2.中（藏）藥資源研究所，甘肅 蘭州 730000

車(chē)前子為車(chē)前科植物車(chē)前Plantago asiaticaL.或平車(chē)前Plantago depressaWilld.的干燥成熟種子。具有清熱、利尿通淋、滲濕止瀉、明目、祛痰等功效，用于治療熱淋澀痛、水腫脹滿(mǎn)、暑濕泄瀉、目赤腫痛、痰熱咳嗽等癥[1]。由于市場(chǎng)需求較大，為保證車(chē)前子的品質(zhì)，選擇適宜地區(qū)進(jìn)行人工栽培是目前行之有效的方法。

最大熵（MaxEnt）模型可根據(jù)物種已知的分布數(shù)據(jù)和環(huán)境變量判斷其生態(tài)需求，并以此進(jìn)行該物種的分布預(yù)測(cè)，結(jié)合地理信息系統(tǒng)ArcGIS分析技術(shù)，可對(duì)物種的潛在分布進(jìn)行預(yù)測(cè)，已廣泛應(yīng)用于中藥相關(guān)領(lǐng)域[2]。本研究通過(guò)查閱資料及實(shí)地采樣，收集車(chē)前和平車(chē)前樣點(diǎn)信息，結(jié)合環(huán)境因子數(shù)據(jù)，采用MaxEnt模型預(yù)測(cè)車(chē)前和平車(chē)前生態(tài)適宜分布區(qū)，并利用ArcGIS軟件將其展示出來(lái)，以期為優(yōu)質(zhì)中藥車(chē)前子人工栽培選地提供依據(jù)。

1 資料與方法

1.1 車(chē)前和平車(chē)前樣點(diǎn)分布數(shù)據(jù)

通過(guò)自采及查閱中國(guó)數(shù)字植物標(biāo)本館（http://www.cvh.ac.cn），收集車(chē)前和平車(chē)前樣點(diǎn)信息，其中車(chē)前325個(gè)、平車(chē)前255個(gè)，記錄其經(jīng)緯度、海拔等分布信息，樣點(diǎn)分布詳見(jiàn)圖1。

圖1 車(chē)前和平車(chē)前樣點(diǎn)分布

1.2 生態(tài)因子數(shù)據(jù)

本研究使用來(lái)源于“中藥資源空間信息網(wǎng)格數(shù)據(jù)庫(kù)”（http://www.tcm-resources.com/）的55個(gè)生態(tài)因子數(shù)據(jù)（降水、溫度、植被類(lèi)型、土壤亞類(lèi)），選取適應(yīng)車(chē)前和平車(chē)前生長(zhǎng)的所有生態(tài)因子數(shù)據(jù)。

1.3 最大熵模型參數(shù)設(shè)置與運(yùn)算

將車(chē)前和平車(chē)前樣點(diǎn)經(jīng)緯度信息導(dǎo)入MaxEnt3.4.0軟件，設(shè)置建模參數(shù)：隨機(jī)測(cè)試比例25%，最大迭代次數(shù)106，收斂門(mén)限0.000 5，其余值默認(rèn)。進(jìn)行迭代運(yùn)算，舍棄貢獻(xiàn)率為0的生態(tài)因子，得到各生態(tài)因子對(duì)車(chē)前和平車(chē)前生長(zhǎng)貢獻(xiàn)率，依據(jù)響應(yīng)曲線(xiàn)得到主生態(tài)因子適宜范圍[3]。用受試者工作特征（ROC）曲線(xiàn)下面積（AUC）評(píng)價(jià)預(yù)測(cè)結(jié)果的精度，取值范圍0～1，若AUC＞0.9，表明此模型精確度較高。

1.4 生態(tài)適宜區(qū)劃分

將MaxEnt運(yùn)算結(jié)果導(dǎo)入ArcGIS軟件，依據(jù)生態(tài)適宜度最小值及正態(tài)分布曲線(xiàn)的均值（μ）和均方差（δ），將車(chē)前和平車(chē)前生長(zhǎng)區(qū)劃分為不適宜區(qū)、次適宜區(qū)和適宜區(qū)3個(gè)等級(jí)，并繪制車(chē)前和平車(chē)前生態(tài)適宜性區(qū)劃分布圖。

2 結(jié)果與分析

2.1 模型預(yù)測(cè)精度

車(chē)前、平車(chē)前的ROC曲線(xiàn)訓(xùn)練數(shù)據(jù)集AUC分別為0.968、0.977，測(cè)試數(shù)據(jù)集AUC分別為0.927、0.940，均大于0.9，表明MaxEnt模型對(duì)預(yù)測(cè)車(chē)前和平車(chē)前的潛在分布準(zhǔn)確度和可信度都很高。車(chē)前、平車(chē)前潛在分布MaxEnt模型ROC曲線(xiàn)見(jiàn)圖2、圖3。

圖2 車(chē)前分布MaxEnt模型ROC曲線(xiàn)

圖3 平車(chē)前分布MaxEnt模型ROC曲線(xiàn)

2.2 主生態(tài)因子選擇

55個(gè)生態(tài)因子中，42個(gè)對(duì)車(chē)前生態(tài)適宜性有貢獻(xiàn)，18個(gè)為影響車(chē)前適宜性的主要生態(tài)因子，貢獻(xiàn)率總和為94.09%，車(chē)前主生態(tài)因子適宜值見(jiàn)表1；47個(gè)對(duì)平車(chē)前生態(tài)適宜性有貢獻(xiàn)，19個(gè)為影響平車(chē)前適宜性的主要生態(tài)因子，貢獻(xiàn)率總和為92.95%，平車(chē)前主生態(tài)因子適宜值見(jiàn)表2。降水量、坡度、溫度、海拔、植被類(lèi)型、土壤亞類(lèi)為影響車(chē)前和平車(chē)前生態(tài)適宜性的主要生態(tài)因子。

表1 車(chē)前主生態(tài)因子貢獻(xiàn)率及適宜值

表2 平車(chē)前主生態(tài)因子貢獻(xiàn)率及適宜值

2.3 生態(tài)適宜性區(qū)劃

車(chē)前生態(tài)適宜性區(qū)劃見(jiàn)表3、圖4。適宜區(qū)包括甘肅省南部邊緣、陜西南部、安徽南部、四川東部、重慶、湖北中西部、浙江、福建、貴州、湖南中南部、江西、廣西北部、廣東南部，次適宜區(qū)包括遼寧南部、河北北部和西部、北京、天津大部、山東、山西、河南南部、寧夏、江蘇中南部、安徽北部、甘肅中部、四川中部、西藏南部邊緣、廣西南部、廣東南部、臺(tái)灣中部。

圖4 車(chē)前生態(tài)適宜性區(qū)劃

表3 車(chē)前生態(tài)適宜性分布區(qū)

平車(chē)前主要生態(tài)適宜性區(qū)劃見(jiàn)圖5。適宜區(qū)包括吉林東部邊緣、遼寧東部、青海東北角、甘肅南部、寧夏南部、陜西大部、山西、河北北部、北京、山東南部、河南西部、四川中部、云南北部、安徽南部，次適宜區(qū)包括黑龍江中部、吉林中部、遼寧西部、內(nèi)蒙古中南部、西藏東南部、青海西北角、甘肅北部、寧夏北部、河北南部、天津北部、山東北部、河南東部、安徽大部、江蘇中西部、浙江中北部、福建中部、江西大部、湖北大部、重慶大部、云南中部、貴州北部、湖南中北部、廣西北部邊緣、廣東北部邊緣、臺(tái)灣西南邊緣。見(jiàn)表4。

圖5 平車(chē)前生態(tài)適宜性區(qū)劃

表4 平車(chē)前生態(tài)適宜性分布區(qū)

續(xù)表4

3 小結(jié)

本研究對(duì)580個(gè)車(chē)前子樣點(diǎn)進(jìn)行生態(tài)適宜性區(qū)劃研究，預(yù)測(cè)模型AUC均大于0.9，表明預(yù)測(cè)結(jié)果良好、可信度高，可為合理選擇種植區(qū)提供依據(jù)。

車(chē)前和平車(chē)前為耐寒、耐旱型植物，對(duì)土壤要求不嚴(yán)格，但更適宜溫暖、潮濕、向陽(yáng)、沙質(zhì)沃土環(huán)境。車(chē)前生長(zhǎng)海拔為0～3 687.5 m，平車(chē)前生長(zhǎng)海拔為388～4 000 m，超出該環(huán)境范圍會(huì)影響植物營(yíng)養(yǎng)吸收；3、4月份為發(fā)芽期，9、10月份為開(kāi)花結(jié)果期，適宜的溫度、充足的降水量能保證植物獲取足夠的營(yíng)養(yǎng)。有報(bào)道，車(chē)前種子重量與年均溫度及年降水量之間關(guān)系顯著，種子重量隨著年均溫度及年降水量的增加而增加[4]，與本研究結(jié)果一致。