高銘 倪淑萍 沈亮

摘 要 目的：為丹參科學(xué)引種和規(guī)范栽培提供參考。方法：查閱全球生物多樣性信息網(wǎng)絡(luò)數(shù)據(jù)庫(kù)等并結(jié)合實(shí)地調(diào)研情況，收集文獻(xiàn)記載的野生分布區(qū)、栽培主產(chǎn)區(qū)及歷史分布區(qū)共697個(gè)丹參樣點(diǎn)，采用Spearman系數(shù)法對(duì)35個(gè)生態(tài)因子進(jìn)行相關(guān)性檢驗(yàn)，并結(jié)合文獻(xiàn)篩選生態(tài)因子指標(biāo)，建立MaxEnt預(yù)測(cè)模型，采用刀切法中的ROC曲線評(píng)價(jià)預(yù)測(cè)模型的準(zhǔn)確度，再根據(jù)模型計(jì)算出丹參在全球范圍內(nèi)的存在概率，進(jìn)而利用ArcGIS 10.3軟件繪制丹參全球范圍內(nèi)的潛在生態(tài)適宜產(chǎn)區(qū)圖。結(jié)果：篩選出6個(gè)環(huán)境因子指標(biāo)，分別為年均溫、最熱季均溫、最冷季均溫、年均降水量、年均日照和年均相對(duì)濕度。刀切法分析表明，氣溫是影響丹參分布的重要生態(tài)因子。丹參全球潛在生態(tài)適宜產(chǎn)區(qū)主要分布在亞洲東部和北美洲中東部等地區(qū)，其生態(tài)適宜產(chǎn)區(qū)面積約為79.16×105 km2。全球范圍內(nèi)適宜丹參種植的國(guó)家主要有中國(guó)、美國(guó)、俄羅斯、日本、韓國(guó)、朝鮮等；我國(guó)的丹參潛在生態(tài)適宜產(chǎn)區(qū)主要包括湖南、湖北、河南、江西、山東等地。結(jié)論：全球范圍內(nèi)丹參潛在生態(tài)適宜產(chǎn)區(qū)主要為亞洲東部及北美洲中東部地區(qū)，在引種栽培過(guò)程中需要注意種植基地的選擇及種植溫度的調(diào)控，以促進(jìn)高品質(zhì)丹參藥材的生產(chǎn)。

關(guān)鍵詞 丹參；MaxEnt模型；生態(tài)適宜產(chǎn)區(qū)；生態(tài)因子；全球

中圖分類號(hào) R282.2 文獻(xiàn)標(biāo)志碼 A 文章編號(hào) 1001-0408（2018）16-2243-05

ABSTRACT OBJECTIVE： To provide reference for scientific introduction and standard cultivation of Salvia miltiorrhiza. METHODS： By reviewing global biodiversity information network database， etc.， based on field research， 697 S. miltiorrhiza sample sites were collected from wild distribution area， main production area and historical distribution area recorded in literature. Spearman coefficient method was used for correlation test of 35 ecological factors. MaxEnt prediction model was established according to literature screening ecological factor index，and using the ROC curve in the Jackknife method to evaluate the accuracy of the prediction model. According to the model， existence probability of S. miltiorrhiza globally was calculated. ArcGIS 10.3 software was used to draw global potential ecological suitable areas of S. miltiorrhiza. RESULTS： Totally 6 environment factor indexes were screened， including mean annual temperature， mean temperature of the hottest season， mean temperature of the coldest season， annual precipitation， annual sunshine and annual relative humidity. Jackknife method showed temperature was an important environmental factor affecting the distribution of S. miltiorrhiza. The potential ecological suitable areas of S. miltiorrhiza were mainly in the eastern Asia and central and easten of North America， and ecological suitable areas for S. miltiorrhiza was 79.16×105 km2. At the global level， S. miltiorrhiza mainly distributed in China， the United States， Russia， Japan， South Korea， North Korea， etc. The potential ecological suitable areas in China mainly included Hunan， Hubei， Henan， Jiangxi， Shandong， etc. CONCLUSIONS： The potential ecological suitable areas for S. miltiorrhiza are mainly in the eastern Asia and central and easten of North America. In the process of introduction and cultivation， great attention should be attached to the selection of planting base and the regulation of the growth temperature of S. miltiorrhiza， so as to promote the production of high quality S. miltiorrhiza.

KEYWORDS Salvia miltiorrhiza； MaxEnt model； Ecological suitable area； Ecological characteristics； Worldwide

丹參為唇形科鼠尾草屬多年生藥用植物丹參（Salvia miltiorrhiza Bge.）的干燥根及根莖，具有活血化瘀、養(yǎng)血安神、降低血壓等功效[1]。丹參主要含有丹參素、丹酚酸、丹參酮、原兒茶醛、咖啡酸、迷迭香酸、紫草酸等化學(xué)成分[2]。近年來(lái)，隨著對(duì)丹參抗氧化、抗凝血、抗血栓等藥理活性的深入研究，其新的藥用價(jià)值不斷被發(fā)現(xiàn)，加之心腦血管疾病患者逐年增加，其原料需求量逐年上升[3]。作為根類藥材，丹參人工種植一直存在連作障礙問(wèn)題，導(dǎo)致丹參可用栽培土地逐年減少，現(xiàn)已嚴(yán)重威脅丹參種植產(chǎn)業(yè)的健康發(fā)展[4]。為實(shí)現(xiàn)該藥材資源的可持續(xù)利用，避免盲目引種擴(kuò)種帶來(lái)的經(jīng)濟(jì)損失，利用生態(tài)位模型開(kāi)展丹參全球潛在生態(tài)適宜產(chǎn)區(qū)預(yù)測(cè)具有重要意義。

生態(tài)位模型是依據(jù)一定數(shù)學(xué)算法，利用物種已知的產(chǎn)地分布數(shù)據(jù)及與其生長(zhǎng)顯著相關(guān)的環(huán)境變量構(gòu)建的產(chǎn)地區(qū)劃模型，主要用來(lái)預(yù)測(cè)物種的潛在適宜分布區(qū)域[5]。該模型的出現(xiàn)為研究物種保護(hù)、進(jìn)化及推廣種植提供了科學(xué)高效的分析工具。目前常用的生態(tài)位預(yù)測(cè)模型主要有最大熵模型（MaxEnt）、基因遺傳算法（GARP）、藥用植物全球產(chǎn)地生態(tài)適宜性區(qū)劃信息系統(tǒng)（GMPGIS）等[6-8]，其中MaxEnt模型是基于最大熵原理提出的生態(tài)位模型，該模型具有預(yù)測(cè)精度高、操作簡(jiǎn)單、運(yùn)算速度快等優(yōu)點(diǎn)，在物種保護(hù)、種植區(qū)劃及有害生物防治等領(lǐng)域得到了廣泛應(yīng)用[9]。本研究應(yīng)用MaxEnt模型分析了全球范圍內(nèi)的丹參潛在生態(tài)適宜區(qū)域，旨在為丹參科學(xué)引種和規(guī)范栽培提供參考，并為高品質(zhì)丹參藥材生產(chǎn)提供依據(jù)。

1 材料與方法

1.1 樣點(diǎn)采集

在全球范圍內(nèi)對(duì)丹參主要分布區(qū)域進(jìn)行樣點(diǎn)選取，包括丹參野生分布區(qū)、栽培主產(chǎn)區(qū)及歷史分布區(qū)[10]。查閱全球生物多樣性信息網(wǎng)絡(luò)數(shù)據(jù)庫(kù)（GBIF，http：//www.gbif.org/）、中國(guó)數(shù)字植物標(biāo)本館（CVH，http：//www.cvh.org.cn/）、英國(guó)丘園（http：//apps.kew.org/）、中國(guó)知網(wǎng)等[11]，并結(jié)合本課題組在山東、陜西等丹參種植產(chǎn)區(qū)調(diào)查的工作記錄[12-13]，剔除產(chǎn)地有誤或重復(fù)的樣點(diǎn)；沒(méi)有經(jīng)緯度信息的標(biāo)本依據(jù)采樣地點(diǎn)采用Google Earth軟件進(jìn)行經(jīng)緯度坐標(biāo)查找，最終共收集了697個(gè)丹參樣點(diǎn)，用于丹參全球潛在生態(tài)適宜產(chǎn)區(qū)分析，詳見(jiàn)圖1[圖1經(jīng)國(guó)家測(cè)繪地理信息局審核，審圖號(hào)為：GS（2016）1766號(hào)，圖4同]。

1.2 氣侯變量的選擇

本研究氣候變量數(shù)據(jù)來(lái)自CliMond數(shù)據(jù)庫(kù)（http：//www.climond.org/）1950-2000年的監(jiān)測(cè)數(shù)據(jù)。該數(shù)據(jù)庫(kù)共有35個(gè)生態(tài)因子，坐標(biāo)系為WGS84，主要包括年均、季均及月均溫/濕度、年均日照、年均降水量等指標(biāo)，氣候柵格數(shù)據(jù)的分辨率為10′ [14]。基礎(chǔ)地理信息數(shù)據(jù)來(lái)源于國(guó)家基礎(chǔ)地理信息中心（http：//www.geodata.cn/Portal/index.jsp）。為了避免各生態(tài)因子間共線性造成模型預(yù)測(cè)結(jié)果的過(guò)度擬合[15]，本研究采用Spearman系數(shù)法對(duì)35個(gè)生態(tài)因子進(jìn)行了相關(guān)性檢驗(yàn)：若兩個(gè)因子的相關(guān)系數(shù)≥0.75，則留下其中生態(tài)意義較為重要的因子，相關(guān)系數(shù)<0.75的因子則全部保留[16]；同時(shí)，依據(jù)影響藥用植物生長(zhǎng)的關(guān)鍵環(huán)境因子，結(jié)合相關(guān)文獻(xiàn)報(bào)道，本研究最終選取了6個(gè)生態(tài)因子用于后續(xù)分析，即年均溫、最熱季均溫、最冷季均溫、年均降水量、年均日照和年均相對(duì)濕度[17-18]。

1.3 模型的建立與評(píng)價(jià)

由于MaxEnt模型可以根據(jù)物種分布的樣點(diǎn)信息和環(huán)境變量數(shù)值分析該物種的生態(tài)位需求，并預(yù)測(cè)其潛在分布地區(qū)[9]，因此本研究擬采用MaxEnt模型進(jìn)行分析。先通過(guò)丹參實(shí)際分布樣點(diǎn)的環(huán)境變量采用MaxEnt 3.3.3軟件建立MaxEnt預(yù)測(cè)模型，再根據(jù)模型計(jì)算出丹參在全球范圍內(nèi)的存在概率，利用ArcGIS 10.3軟件處理得到丹參全球范圍內(nèi)的潛在生態(tài)適宜產(chǎn)區(qū)。MaxEnt 3.3.3軟件由其網(wǎng)站下載（http：//www.cs.princeton.edu/～schapire/MaxEnt/）；ArcGIS 10.3軟件由美國(guó)ESRI公司提供。具體操作方法為：首先，將丹參分布樣點(diǎn)數(shù)據(jù)和環(huán)境因子數(shù)據(jù)進(jìn)行格式轉(zhuǎn)換，然后導(dǎo)入MaxEnt 3.3.3軟件中，運(yùn)行參數(shù)設(shè)置為隨機(jī)選取25%的分布點(diǎn)作為測(cè)試集（Test data），剩余的樣點(diǎn)數(shù)據(jù)作為訓(xùn)練集（Training data）進(jìn)行模型準(zhǔn)確性分析，最大迭代次數(shù)為500[6]。

刀切法（Jackknife）是一種評(píng)估模型可靠性的分析方法，因此本研究采用刀切法中的ROC曲線（Receiver operating characteristic curve）來(lái)評(píng)價(jià)預(yù)測(cè)模型的準(zhǔn)確度[19]。ROC曲線是目前物種潛在分布預(yù)測(cè)模型評(píng)價(jià)中應(yīng)用較為廣泛的模型[20]。ROC曲線下面積（Area under curve，AUC）數(shù)值范圍一般為0.5～1.0，其值越大表示預(yù)測(cè)越精確，通常數(shù)值在0.8～0.9為好，0.9～1.0為非常好[21]。

通常進(jìn)行生態(tài)適宜產(chǎn)區(qū)劃分的方式有很多種，本研究根據(jù)丹參生長(zhǎng)的存在概率，以劃分效果較好的人工分級(jí)方法為基礎(chǔ)，將丹參生態(tài)適宜產(chǎn)區(qū)劃分為3級(jí)，分別為高適生區(qū)（存在概率為0.6～1）、中適生區(qū)（存在概率為0.4～<0.6）及不宜適生區(qū)（存在概率為0～<0.4）[22]。

采用刀切法在MaxEnt 3.3.3軟件中生成氣候因子響應(yīng)圖，本研究設(shè)定存在概率大于0.5的結(jié)果數(shù)值為物種最適宜的生態(tài)因子范圍[23]。

2 結(jié)果

2.1 影響丹參生長(zhǎng)的生態(tài)因子及ROC曲線

基于丹參在全球分布區(qū)的697個(gè)樣點(diǎn)，經(jīng)過(guò)MaxEnt模型分析得到丹參各生態(tài)因子統(tǒng)計(jì)參數(shù)，詳見(jiàn)表1。由表1可知，丹參分布地區(qū)的最冷季均溫?cái)?shù)值變異范圍較大，其次是年均降水量，而最熱季均溫和年均日照數(shù)值變異范圍較小。本研究得到的ROC訓(xùn)練集及測(cè)試集的AUC值均為0.979（如圖2所示），說(shuō)明模型的預(yù)測(cè)效果達(dá)到了極好的測(cè)試水平，因此通過(guò)該模型預(yù)測(cè)得到的丹參生態(tài)適宜產(chǎn)區(qū)具有較高的準(zhǔn)確性和可信度。刀切法分析表明，最熱季均溫的貢獻(xiàn)度最大（為34.47%），其次是最冷季均溫（為31.14%），另外年均降水量及年均日照的累積貢獻(xiàn)度約為31.33%，以這4個(gè)指標(biāo)繪制生態(tài)因子響應(yīng)圖，詳見(jiàn)圖3。由圖3可知，存在概率大于0.5的適宜指標(biāo)結(jié)果為：最適宜丹參生長(zhǎng)的地區(qū)最熱季均溫范圍為22.5～27.5 ℃，最冷季均溫范圍為-12～8 ℃，年均降水量范圍為700～1 800 mm，年均日照范圍為130～160 W/m2。

2.2 丹參全球潛在生態(tài)適宜產(chǎn)區(qū)

由上文分析可得丹參全球潛在生態(tài)適宜產(chǎn)區(qū)，詳見(jiàn)圖4。由圖4可知，丹參全球潛在生態(tài)適宜產(chǎn)區(qū)主要包括亞洲東部、歐洲南部、北美洲中東部以及非洲東北部地區(qū)，其全球潛在生態(tài)適宜產(chǎn)區(qū)面積約為79.16×105 km2，其中高適生區(qū)主要分布在中國(guó)東部、美國(guó)中東部地區(qū)。適宜丹參種植的國(guó)家主要有中國(guó)、美國(guó)、俄羅斯、日本、韓國(guó)、朝鮮、哈薩克斯坦等，其中美國(guó)適宜產(chǎn)區(qū)面積為30.53×105 km2、中國(guó)適宜產(chǎn)區(qū)面積為27.11×105 km2、俄羅斯適宜產(chǎn)區(qū)面積為4.45×105 km2。

2.3 丹參在我國(guó)的潛在生態(tài)適宜產(chǎn)區(qū)

利用ArcGIS 10.3軟件將MaxEnt 3.3.3軟件運(yùn)行的結(jié)果與我國(guó)省級(jí)行政分區(qū)圖疊加，得到丹參在我國(guó)的潛在生態(tài)適宜產(chǎn)區(qū)及各省適宜種植面積。結(jié)果，丹參潛在生態(tài)適宜產(chǎn)區(qū)主要分布于華北平原、長(zhǎng)江中下游平原以及新疆、遼寧的部分地區(qū)；其中，湖南、湖北、河南、江西、山東、四川、安徽、陜西、河北等省的適宜種植面積均超過(guò)1.50×105 km2。

3 討論

丹參作為根莖類藥材，存在較為嚴(yán)重的連作障礙問(wèn)題，重茬種植丹參易導(dǎo)致其抗逆性（抗旱、澇、蟲(chóng)?。┫陆?、產(chǎn)量減少、外觀粗糙、須根增多等問(wèn)題；隨著種植年限增加，容易導(dǎo)致死苗、根腐爛、產(chǎn)生大量寄生線蟲(chóng)等問(wèn)題[12]。開(kāi)展丹參生態(tài)適宜產(chǎn)區(qū)預(yù)測(cè)，可有效指導(dǎo)丹參科學(xué)選地及生產(chǎn)高品質(zhì)藥材。通過(guò)物種生態(tài)預(yù)測(cè)模型可有效預(yù)測(cè)物種的生態(tài)適宜產(chǎn)區(qū)，其中MaxEnt模型因其設(shè)置的環(huán)境約束條件，使其預(yù)測(cè)結(jié)果相對(duì)保守，但準(zhǔn)確性較高，預(yù)測(cè)面積適中，應(yīng)用較為廣泛[24]，因此本研究采用MaxEnt模型對(duì)丹參生態(tài)適宜產(chǎn)區(qū)進(jìn)行預(yù)測(cè)。結(jié)果表明，該模型訓(xùn)練集和測(cè)試集的AUC值均大于0.9，說(shuō)明模型預(yù)測(cè)結(jié)果準(zhǔn)確度較高，適用于丹參藥材生態(tài)適宜產(chǎn)區(qū)的預(yù)測(cè)分析。另外，該模型預(yù)測(cè)的適宜產(chǎn)區(qū)基本包括了目前我國(guó)已有的丹參種植產(chǎn)區(qū)，如陜西商洛、河南方城、四川中江、遼寧凌源、河北撫寧、山東平邑、江蘇南京等地[13，25-26]，進(jìn)一步驗(yàn)證了本研究MaxEnt模型預(yù)測(cè)的準(zhǔn)確性。

本研究結(jié)果表明，丹參潛在生態(tài)適宜產(chǎn)區(qū)主要分布在亞洲東部和北美洲中東部地區(qū)，全球范圍內(nèi)丹參生態(tài)相似度最高的地區(qū)主要集中在中國(guó)、美國(guó)、俄羅斯、日本、韓國(guó)、朝鮮等地。我國(guó)的丹參潛在生態(tài)適宜產(chǎn)區(qū)主要包括湖南、湖北、河南、江西、山東、四川、安徽、陜西、河北等地，另外浙江、貴州等地也適宜引種丹參進(jìn)行栽培。楊先富等[27]研究表明，在貴州地區(qū)引種的丹參藥材質(zhì)量在10～11月份品質(zhì)較好，各丹參酮成分含量較高，均高于2015年版《中國(guó)藥典》規(guī)定，這也說(shuō)明了本研究預(yù)測(cè)產(chǎn)區(qū)的準(zhǔn)確性。魏志華等[28]采用中藥材產(chǎn)地適宜性分析地理信息系統(tǒng)（TCMGIS）對(duì)裕丹參在我國(guó)的生態(tài)適宜產(chǎn)區(qū)進(jìn)行了區(qū)劃，結(jié)果表明，裕丹參除河南外，在四川、陜西、湖北等地均可生長(zhǎng)。何淑婷等[25]采用MaxEnt模型對(duì)南丹參在我國(guó)的潛在分布區(qū)進(jìn)行了預(yù)測(cè)，結(jié)果表明南丹參適宜生長(zhǎng)的地區(qū)包括湖南、江西、浙江和福建等地。景鵬飛[26]采用MaxEnt模型研究表明，我國(guó)丹參的生態(tài)適宜分布區(qū)主要在華北地區(qū)。本研究預(yù)測(cè)得到的我國(guó)丹參生態(tài)適宜產(chǎn)區(qū)基本包含了前人的預(yù)測(cè)地區(qū)，且預(yù)測(cè)所得全球生態(tài)適宜產(chǎn)區(qū)面積更大、覆蓋范圍更廣，對(duì)未來(lái)丹參引種及科學(xué)種植具有一定參考價(jià)值。

中藥材產(chǎn)量及有效成分與其生長(zhǎng)環(huán)境生態(tài)因子密切相關(guān)[29]。謝彩香等[30]研究表明，適當(dāng)?shù)蜏靥幚碛欣谌藚⒃碥粘煞值目焖俜e累。本研究結(jié)果表明，影響丹參分布的主要生態(tài)因子有最熱季均溫、最冷季均溫、年均降水量和年均日照。其中，最熱季均溫的貢獻(xiàn)率高達(dá)34.47%，其次是最冷季均溫（31.14%），由此表明這兩個(gè)因子是影響丹參生長(zhǎng)最重要的生態(tài)因子。由我國(guó)的丹參潛在生態(tài)適宜產(chǎn)區(qū)圖可以看出，丹參適宜生長(zhǎng)在溫帶地區(qū)，在我國(guó)最北端的適宜產(chǎn)區(qū)以北京和河北周邊地區(qū)為界限，這說(shuō)明氣溫對(duì)丹參的分布影響較大。張辰露等[31]研究表明，在溫暖濕潤(rùn)、晝夜溫差相對(duì)小的地區(qū)，丹參中丹參酮類成分的含量較高，而隨著緯度及海拔升高，其含量逐漸降低。這與本研究得出的影響丹參的主要生態(tài)因子為氣溫的結(jié)論相符。因此，高海拔及高緯度地區(qū)不適宜丹參栽培，在引種過(guò)程中需要注意種植基地的選擇及種植溫度的調(diào)控。

綜上所述，本研究預(yù)測(cè)的丹參潛在生態(tài)適宜產(chǎn)區(qū)可為未來(lái)丹參藥材引種和栽培提供科學(xué)參考。但由于影響丹參生長(zhǎng)的環(huán)境因素眾多，本研究預(yù)測(cè)的部分產(chǎn)區(qū)實(shí)際可能不適宜丹參生長(zhǎng)，因此需要先進(jìn)行小范圍栽培試驗(yàn)，成功后方可大規(guī)模推廣種植。

參考文獻(xiàn)

[ 1 ] 國(guó)家藥典委員會(huì). 中華人民共和國(guó)藥典：一部[S]. 2015年版. 北京：中國(guó)醫(yī)藥科技出版社，2015：76-77.

[ 2 ] 柴瑞震. 丹參的藥理研究近況[J]. 中國(guó)中醫(yī)藥科技，2003，10（6）：390-391.

[ 3 ] 李耿，孟繁蘊(yùn)，楊洪軍，等. UPLC 法同時(shí)測(cè)定丹參中 11 種成分的含量[J]. 中國(guó)藥房，2014，25（19）：1766-1768.

[ 4 ] 張辰露，孫群，葉青. 連作對(duì)丹參生長(zhǎng)的障礙效應(yīng)[J]. 西北植物學(xué)報(bào)，2005，25（5）：1029-1034.

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

[ 6 ] ZHANG Q，WEN J，CHANG ZQ，et al. Evaluation and prediction of ecological suitability of medicinal plant American Ginseng （Panax quinquefolius）[J]. Chin Herb Med，2018，10（1）：80-85.

[ 7 ] RAY D，BEHERA MD，JACOB J. Evaluating ecological niche models：a comparison between Maxent and GARP for predicting distribution of Hevea brasiliensis in India [J]. Proc Natl Acad Sci，India Sect B，2017，14（3）：1-7.

[ 8 ] DU Z，WU J，MENG X，et al. Predicting the global potential distribution of four endangered Panax species in middle-and low-latitude regions of China by the Geographic Information System for Global Medicinal Plants （GMPGIS）[J]. Molecules，2017，22（10）：e1630.

[ 9 ] PHILLIPS SJ，ANDERSON RP，SCHAPIRE RE. Maximum entropy modeling of species geographic distributions[J]. Ecol Model，2006，190（3/4）：231-259.

[10] 陳士林. 中國(guó)藥材產(chǎn)地生態(tài)適宜性區(qū)劃[M].2版. 北京：科學(xué)出版社，2017：14.

[11] 沈亮，吳杰，李西文. 北細(xì)辛全球產(chǎn)地生態(tài)適宜性分析及品質(zhì)生態(tài)學(xué)研究[J]. 世界科學(xué)技術(shù)（中醫(yī)藥現(xiàn)代化），2016，18（8）：1295-1302.

[12] 趙魁，郭曉恒，宋杰，等. 全國(guó)丹參生產(chǎn)現(xiàn)狀的調(diào)查和分析[J]. 時(shí)珍國(guó)醫(yī)國(guó)藥，2010，21（9）：2307-2310.

[13] 王志芬，劉喜民. 山東省丹參藥材產(chǎn)業(yè)現(xiàn)狀與發(fā)展策略[J]. 山東農(nóng)業(yè)科學(xué)，2012，44（8）：131-132、136.

[14] KRITICOS DJ，WEBBER BL，LERICHE A，et al. CliMond：global high resolution historical and future scenario climate surfaces for bioclimatic modelling[J]. Methods Ecol Evol，2012，3（1）：53-64.

[15] GRAHAM MH. Confronting multicollinearity in ecological multiple regression[J]. Ecology，2003，84（11）：2809- 2815.

[16] ZHANG MG，SLIK JWF，MA KP. Using species distribution modeling to delineate the botanical richness patterns and phytogeographical regions of China[J]. Sci Rep，2016. DOI：org.10.1038/srep22400.

[17] 沈亮，吳杰，李西文，等. 人參全球產(chǎn)地生態(tài)適宜性分析及農(nóng)田栽培選地規(guī)范[J]. 中國(guó)中藥雜志，2016，41（18）：3314-3322.

[18] 孟祥霄，黃林芳，董林林，等. 三七全球產(chǎn)地生態(tài)適宜性及品質(zhì)生態(tài)學(xué)研究[J]. 藥學(xué)學(xué)報(bào)，2016，51（9）：1483- 1493.

[19] DORFMAN DD，BERBAUM KS，METZ CE. Receiver operating characteristic rating analysis：generalization to the population of readers and patients with the Jackknife method[J]. Invest Radiol，1992，27（9）：723-731.

[20] PETERSON AT，MONICA P，EATON M. Transferability and model evaluation in ecological niche modeling：a comparison of GARP and MaxEnt[J]. Ecography，2007，30（4）：550-560.

[21] SWETS JA. Measuring the accuracy of diagnostic systems[J]. Science，1988，240（4857）：1285-1293.

[22] YANG XQ，KUSHWAHA SPS，SARAN S，et al. Maxent modeling for predicting the potential distribution of medicinal plant，Justicia adhatoda L. in Lesser Himalayan foothills[J]. Ecol Eng，2013，51（1）：83-87.

[23] HEMERY LG，MARION SR，ROMSOS CG，et al. Henkel ecological niche and species distribution modelling of sea stars along the Pacific Northwest continental shelf[J]. Divers Distrib，2016，22（12）：1314-1327.

[24] 張琴，宋經(jīng)元，邵飛，等. 防風(fēng)固沙優(yōu)良樹(shù)種歐李的潛在適生區(qū)及生態(tài)特征[J]. 北京林業(yè)大學(xué)學(xué)報(bào)，2018，40（3）：66-74.

[25] 何淑婷，白碧玉，但佳惠，等. 基于MaxEnt的南丹參在中國(guó)的潛在分布區(qū)預(yù)測(cè)及適生性分析[J]. 安徽農(nóng)業(yè)科學(xué)，2014，42（8）：2311-2314.

[26] 景鵬飛. 丹參組藥用植物資源調(diào)查和多樣性評(píng)價(jià)研究[D]. 西安：陜西師范大學(xué)，2013.

[27] 楊先富，徐文靜，簡(jiǎn)衛(wèi)光，等. 貴州引種丹參質(zhì)量動(dòng)態(tài)變化研究[J]. 中國(guó)現(xiàn)代中藥，2018，20（3）：301-304.

[28] 魏志華，王新民，喬卿梅，等. 裕丹參藥材適宜產(chǎn)區(qū)區(qū)劃研究[J]. 廣東農(nóng)業(yè)科學(xué)，2011，38（17）：24-25.

[29] 王歡，李慧，曾凡琳，等. 黃花蒿空間分布及全球潛在氣候適宜區(qū)[J]. 中藥材，2015，38（3）：460-466.

[30] 謝彩香，索風(fēng)梅，賈光林，等. 人參皂苷與生態(tài)因子的相關(guān)性[J]. 生態(tài)學(xué)報(bào)，2011，31（24）：7551-7563.

[31] 張辰露，梁宗鎖，郭宏波，等. 不同氣候區(qū)丹參生物量、有效成分變化與氣象因子的相關(guān)性研究[J]. 中國(guó)中藥雜志，2015，40（4）：607-613.

（收稿日期：2018-02-24 修回日期：2018-07-02）

（編輯：孫 冰）