何淑婷　 白碧玉　 但佳惠等

摘要 [目的]對(duì)南丹參在中國(guó)的潛在分布區(qū)進(jìn)行預(yù)測(cè)，并對(duì)適生性進(jìn)行分析。[方法]通過(guò)收集中草藥南丹參的分布信息，利用最大熵模型（MaxEnt）與地理信息系統(tǒng)（ArcGIS），綜合相關(guān)68項(xiàng)環(huán)境因子，分析南丹參在中國(guó)的分布信息。[結(jié)果]全國(guó)適合于南丹參分布的地區(qū)主要集中在湖南、江西、浙江和福建等地，其高適生區(qū)面積依次為湖南省（319.44 km2 ）、江西?。?51.39 km2）、浙江?。?22.92 km2）、福建?。?2.5 km2）；刀切法（Jackknife）分析顯示，南丹參分布主要受到降水的影響。[結(jié)論]該研究結(jié)果為南丹參的野生資源利用及其他南丹參相關(guān)研究提供了一定技術(shù)手段和理論依據(jù)。

關(guān)鍵詞 南丹參（Salvia bowleyana Dunn.）；潛在適生區(qū)；MaxEnt；ArcGIS；生態(tài)位

中圖分類(lèi)號(hào) S567.5+3 文獻(xiàn)標(biāo)識(shí)碼

A 文章編號(hào) 0517-6611（2014）08-02311-04

Prediction of Potential Distribution Areas of Salvia bowleyana Dunn. in China Based on MaxEnt and Suitability Analysis

HE Shuting， JING Pengfei et al （School of Life Science， Shaanxi Normal University， Xian， Shaanxi 710062）

Abstract [Objective] To predict potential distribution area of Salvia bowleyana Dunn. in China and analyze its suitability. [Method] The distribution information of Salvia bowleyana Dunn. and 68 related environmental factors were collected. The potential suitable distribution range of Salvia bowleyana Dunn. in China was analyzed by the software of MaxEnt and ArcGIS. [Result] The results showed that Salvia bowleyana Dunns suitable distribution areas are mainly concentrated in Hunan， Jiangxi， Zhejiang， Fujian. The high suitable area of province above are as follows： Hunan Province（319.44 km2）， Jiangxi Province （251.39 km2）， Zhejiang Province （122.92 km2）， Fujian Province （62.5 km2）. Jackknife analysis showed that precipitation has obvious influence on Salvia bowleyana Dunns distributions. [Conclusion] The results can provide reference for Salvia bowleyana Dunns wild resource utilization and study.

Key words Salvia bowleyana Dunn.； Potential distribution； MaxEnt； ArcGIS； niche

南丹參（Salvia bowleyana Dunn.）系唇形科（Labiaceae）鼠尾草屬（Salvia）多年生草本植物，根肥厚入藥，功效同丹參，可代中藥丹參（S.salviae）入藥，具有祛瘀、活血、調(diào)經(jīng)和抗氧化等功效；主產(chǎn)于浙江、湖南、江西、福建、廣東和廣西等南部地方；生于山地、山谷、路旁、林下或水邊，海拔30～960 m[1-3]，在浙江和江西等一些南方地區(qū)南丹參作丹參藥用[4-5]。施天慧等對(duì)云南產(chǎn)的丹參資源植物化學(xué)成分進(jìn)行測(cè)定，結(jié)果顯示南丹參中總酚酸的含量為4.53%，丹酚酸B為41.33 μg/g，與丹參中含的量接近，具有較高的利用價(jià)值[7]。根據(jù)《中國(guó)藥典》2005版中用丹酚酸B作為評(píng)價(jià)丹參藥材中的酚酸類(lèi)成分的指標(biāo)，并規(guī)定丹參中丹酚酸B的含量不低于30 mg/g[8]。雖然南丹參還沒(méi)有被藥典收錄，但有研究表明南丹參（浙江產(chǎn)）中丹酚酸B含量（82.52 mg/g）較高，不僅符合藥典的規(guī)定，而且比丹參（云南產(chǎn)）中丹酚酸B含量（45.38 mg/g）高，是丹參藥材中的優(yōu)良品種[9]。南丹參在常用8種丹參品系中生長(zhǎng)冠幅最大，單株產(chǎn)量及有效成分含量（丹參素）最高，根的優(yōu)質(zhì)品最高，抗根結(jié)線蟲(chóng)病能力為中度感病，屬3級(jí)，是最為適宜推廣種植的種質(zhì)（品種）[10]。

最大熵理論（the principle of maximum entropy）是一種基于有限的已知信息對(duì)未知分布進(jìn)行無(wú)偏推斷的預(yù)測(cè)方法。對(duì)其預(yù)測(cè)結(jié)果進(jìn)行AUC（ areas under curve）分析后顯示，其結(jié)果要優(yōu)于同類(lèi)預(yù)測(cè)模型如GARP、Climex和BIOCLIM，特別是在物種分布數(shù)據(jù)不全的情況下，MaxEnt仍然能得到較為滿意的結(jié)果[17-21]。目前，該模型已成功應(yīng)用于預(yù)測(cè)保護(hù)物種、入侵種、以及有害生物和藥用植物的潛在生境分布[21-24]。

隨著近年來(lái)心腦血管疾病患者的增加及南丹參新的藥理學(xué)活性的逐步發(fā)現(xiàn)，其利用價(jià)值和需求量也明顯提高[12]。為了充分利用南丹參的野生資源，適應(yīng)南丹參的收集和生產(chǎn)開(kāi)發(fā)需要，筆者對(duì)南丹參的潛在分布區(qū)進(jìn)行預(yù)測(cè)，并對(duì)其適生性進(jìn)行分析，以期為南丹參的野生資源利用及其他南丹參相關(guān)研究提供理論依據(jù)。

1 數(shù)據(jù)收集

1.1 南丹參分布數(shù)據(jù)收集和處理 試驗(yàn)需要的所有南丹參分布數(shù)據(jù)來(lái)源于中國(guó)數(shù)字標(biāo)本館（CVH）（http：//www.cv h.or g.cn/cms/）、中國(guó)自然保護(hù)區(qū)資源平臺(tái)（http：//www.papc.cn/）、全球生物多樣性信息平臺(tái)（GBIF）（http：//www.gbif.org/）、教學(xué)標(biāo)本標(biāo)準(zhǔn)化整理整合與資源共享平臺(tái)（http：//mnh.scu.edu.cn/）和中國(guó)國(guó)內(nèi)21所高校及研究所標(biāo)本館（華中師范大學(xué)CCNU、中科院成都生物所CDBI、湖北中醫(yī)藥大學(xué)ECM等）。經(jīng)過(guò)對(duì)以上平臺(tái)南丹參數(shù)據(jù)的收集，將處理所得的所有數(shù)據(jù)轉(zhuǎn)換為經(jīng)緯度坐標(biāo)。對(duì)于沒(méi)有記錄經(jīng)緯度坐標(biāo)的分布點(diǎn)，均以最小可查的行政單位為準(zhǔn)。根據(jù)MaxEnt 軟件要求，將實(shí)際分布點(diǎn)按物種名、分布點(diǎn)經(jīng)度和緯度順序儲(chǔ)存成后綴名為.csv格式的文件，其中東經(jīng)和北緯為正，西經(jīng)和南緯為負(fù)。

1.2 環(huán)境變量及相關(guān)地理信息數(shù)據(jù) 查詢(xún)世界氣候數(shù)據(jù)庫(kù)（http：//www.worldclim.org），精度2.5弧分 （2.5 km），得到研究所需的氣候環(huán)境圖層數(shù)據(jù)，共6類(lèi)，即地形因子：alt（海拔）；生物氣候因子：bio1（年均溫）、 bio2（晝夜溫差月均值）、 bio3（等溫性[（Bio2 / Bio7） ×100]）、 bio4（溫度季節(jié)性變化的標(biāo)準(zhǔn)差）、bio5（最暖月最高溫）、bio6（最冷月最低溫）、bio7（年均溫變化范圍）、bio8（最濕季度平均溫度）、bio9（最干季度平均溫度）、bio10（最暖季度平均溫度）、bio11（最冷季度平均溫度）、bio12（年均降水量）、bio13（最濕月降水量）、bio14（最干月降水量）、bio15（降水量變異系數(shù)）、bio16（最濕季度降水量）、bio17（最干季度降水量）、bio18（最暖季度降水量）、bio19（最冷季度降水量）；單月因子：tmax1-tmax12（1～12月最高溫）； tmean1～tmean12（1～12月平均氣溫）；tmin1～tmin12（1～12月最低氣溫）； prec1～prec12（1～12月降雨量）。在國(guó)家地理信息系統(tǒng)網(wǎng)站（http：//nfgis.nsdi.gov.cn/）下載研究所采用的1：400萬(wàn)中國(guó)地圖及中國(guó)行政區(qū)劃圖。

2 研究方法

2.1 數(shù)據(jù)處理 將獲得的南丹參分布數(shù)據(jù)和相應(yīng)的環(huán)境數(shù)據(jù)按相應(yīng)的格式轉(zhuǎn)換后導(dǎo)入MaxEnt，之后隨機(jī)選取25%的分布點(diǎn)作為測(cè)試集（test data），剩余的作為訓(xùn)練集（training data），其他參數(shù)均為軟件默認(rèn)值。

2.2 模型預(yù)測(cè)評(píng)價(jià) 設(shè)定環(huán)境圖層導(dǎo)入的每個(gè)環(huán)境因子選擇離散（categorical）變量模式，用刀切法將68項(xiàng)環(huán)境因子在模型中按訓(xùn)練得分排列后，選出得分較高的前15項(xiàng)進(jìn)行進(jìn)一步分析，評(píng)價(jià)影響南丹參分布的最關(guān)鍵因素。

采用MaxEnt軟件中受試者工作曲線方法進(jìn)行模型精度檢驗(yàn)，Hanley（1882）闡釋了ROC曲線的意義及使用方法，曲線下的面積（AUC值）代表著模型的測(cè)試能力[19]。AUC值在0.5～0.6，表明模型模擬效果為失?。?.6～0.7表明效果較差； 0.7～0.8表明模擬效果一般；0.8～0.9表明模擬效果好；0.9～1.0表明模擬效果非常好[20]。

2.3 環(huán)境因子權(quán)重分析 用刀切法分析各個(gè)環(huán)境因子對(duì)預(yù)測(cè)結(jié)果的影響程度，來(lái)對(duì)影響物種分布的環(huán)境因素進(jìn)行權(quán)重排序。該方法用所有的變量建立模型（深灰色條帶表示）；單獨(dú)用每一變量建立模型（黑色條帶表示）；每次排除一個(gè)變量，用剩余的變量建立模型（淺灰色條帶表示），觀察這個(gè)模型（淺灰色條帶）與所有變量都存在的模型（深灰色條帶）之間的變化，如果變化大，說(shuō)明排除的那個(gè)變量含有很重要的信息，反之，有用信息少[24]。

3 結(jié)果與分析

3.1 南丹參在中國(guó)潛在分布和分析 利用GIS的空間分析工具，采用人工（Manual）分級(jí)方法分析得到南丹參在中國(guó)的分布圖（圖1），以南丹參在中國(guó)的潛在分布圖（ESRI ASCII格式）為底圖，掩膜抽提各?。ㄊ些p自治區(qū)）的分布圖，采用Manual分級(jí)方法將各個(gè)省份的潛在分布區(qū)面積進(jìn)行詳細(xì)統(tǒng)計(jì)（表2）。圖1表明，南丹參在我國(guó)的適生范圍主要集中在華南、華中和東南地區(qū)，高適生范圍分布在地處24°～31°N，109°～120°E的湖南東南部，江西西北部，浙江北部和福建東北部，同時(shí)西藏、陜西、新疆和甘肅等西北地區(qū)也有少量分布。由表2可知，高適生區(qū)依次為湖南?。?19.44 km2 ）、江西省（251.39 km2）、浙江省（122.92 km2）、福建省（62.5 km2）和廣西?。?6.53 km2）；較適宜生長(zhǎng)區(qū)包括湖南?。? 421.25 km2）、江西省（1 272.22 km2）、浙江?。?91.67 km2）、安徽?。?45.83 km2）和福建省（268.75 km2）；內(nèi)蒙古、吉林、天津、遼寧、青海和黑龍江等省區(qū)則基本沒(méi)有其適生區(qū)分布。其中，湖南省的高適生區(qū)和較適宜生長(zhǎng)區(qū)主要分布在長(zhǎng)沙市長(zhǎng)沙縣、瀏陽(yáng)縣和岳陽(yáng)市岳陽(yáng)縣境內(nèi)；江西省主要分布在東鄉(xiāng)縣、九江縣以及宜豐縣境內(nèi)；浙江省主要分布在蕭山市、紹山市和富陽(yáng)縣、建德縣境內(nèi)；福建省主要分布在福州市、崇安縣、建寧縣境內(nèi)。

3.2 預(yù)測(cè)結(jié)果評(píng)價(jià) 經(jīng)受試者工作曲線檢驗(yàn)，MaxEnt模型的ROC曲線下面積訓(xùn)練集和測(cè)試集的AUC值依次為0.975和0.968，接近于1，表明MaxEnt模型對(duì)南丹參中國(guó)潛在分布區(qū)的預(yù)測(cè)效果非常好（圖3）。

3.3 影響南丹參分布的不同環(huán)境因子的權(quán)重評(píng)價(jià)和相關(guān)性分析 對(duì)最初的68項(xiàng)環(huán)境因子進(jìn)行刀割法分析，其中測(cè)試得分最高的前15項(xiàng)得分由高到低依次為：溫差季節(jié)性變化標(biāo)準(zhǔn)差（bio4）、年均降水量（bio12）、最濕季節(jié)降水量（bio16）、海拔（alt）、最冷季節(jié)降水量（bio19）、4月平均降水量（prec4）、5月平均降雨量（prec5）、3月平均降雨量（prec3）、6月平均降雨量（prec6）、最干月降水量（bio14）、2月平均降雨量（prec2）、最干季度降水量（bio17）、9月份最低溫（tmin9）、最暖季度降水量（bio18）、最干季度平均溫度（bio9）（圖4）。這15項(xiàng)環(huán)境因子對(duì)預(yù)測(cè)結(jié)果有較大的影響。結(jié)果表明，與降雨量相關(guān)的環(huán)境變量對(duì)南丹參的地理分布有重要影響，比如溫差季節(jié)性變化標(biāo)準(zhǔn)差、年均降水量、最濕季節(jié)降水量等。

鑒于環(huán)境因子之間存在相關(guān)性，試驗(yàn)利用ArcGIS中的Band Collection Statistics工具對(duì)得分較高的幾項(xiàng)環(huán)境因子進(jìn)行了相關(guān)性分析。由于年均降水量（bio12）、最濕季節(jié)降水量（bio16）、最冷季節(jié)降水量（bio19）、4月平均降水量（prec4）這幾項(xiàng)環(huán)境因子的相關(guān)性較高，所以選擇相關(guān)性較低而權(quán)重較高的幾項(xiàng)環(huán)境因子進(jìn)行分析，得出對(duì)南丹參適宜生長(zhǎng)區(qū)域生態(tài)位影響較大的環(huán)境因子是：溫度季節(jié)性變化的標(biāo)準(zhǔn)差（bio4）700～850；最干季度降水量（bio17）15～20 mm；海拔（alt）0～200 m；9月份最低溫（tmin9）20.5～22.0 ℃；最干季度平均溫度（bio9）10～14 ℃。

圖3 刀切法的環(huán)境變量重要性分析

4 結(jié)論與討論

試驗(yàn)利用MaxEnt模型與地理信息系統(tǒng)對(duì)南丹參在我國(guó)的適生區(qū)進(jìn)行預(yù)測(cè)，直觀且定量地獲得了南丹參在我國(guó)的潛在分布區(qū)，其預(yù)測(cè)結(jié)果經(jīng)ROC曲線分析法驗(yàn)證，得到模型的AUC值為0.968，預(yù)測(cè)效果相當(dāng)好。經(jīng)Jacknife對(duì)環(huán)境因子分析后得出，與降雨量相關(guān)的環(huán)境變量對(duì)南丹參的地理分布有關(guān)鍵性的影響。預(yù)測(cè)結(jié)果得出，南丹參的高適生區(qū)主要分布于湖南、江西、浙江和福建等地；較適宜生長(zhǎng)區(qū)主要分布在華東和華南地區(qū)（浙江、安徽、廣西、廣東和福建等）；四川省、重慶市以及陜西省等西北和西南地區(qū)也有分布。預(yù)測(cè)結(jié)果與南丹參主要產(chǎn)地（浙江、湖南、江西、福建和廣東）相吻合，表明這些地區(qū)的氣候適合南丹參生長(zhǎng)繁殖，南丹參在這些地區(qū)的野生資源也可能比較豐富。

試驗(yàn)得出南丹參在浙江省適生區(qū)面積為91 624.98 km2，高適生面積為122.92 km2 ，云南省只有低適生區(qū)，面積為341 688 km2。而施天慧等對(duì)云南產(chǎn)的南丹參化學(xué)成分測(cè)定，結(jié)果顯示其丹酚酸B為41.33 μg/g[7]。李旻輝等人對(duì)浙江產(chǎn)南丹參中測(cè)定丹酚酸B含量為82.52 mg/g[9]，得出浙江產(chǎn)南丹參的丹酚酸B含量比云南產(chǎn)的南丹參要高。故推測(cè)南丹參在潛在適生區(qū)生長(zhǎng)，其個(gè)體發(fā)育狀況會(huì)比其非適生區(qū)好，且有效活性物質(zhì)的積累量會(huì)更高。因此建議：采收和使用南丹參時(shí)可以參考研究中得出的南丹參野生資源理論分布點(diǎn)；積極推廣、擴(kuò)大、保護(hù)南丹參野生藥用資源。研究得到南丹參藥材的適生區(qū)可以對(duì)其進(jìn)行野生資源保護(hù)，生態(tài)環(huán)境保護(hù)，合理采購(gòu)、并且結(jié)合其生長(zhǎng)生理特性進(jìn)行可持續(xù)資源開(kāi)發(fā)利用。

綜上所述，研究采集了相當(dāng)數(shù)量的南丹參分布數(shù)據(jù)，基本考慮到絕大多數(shù)南丹參的實(shí)際生態(tài)位，從理論上明確了南丹參的具體適生區(qū)，可以為南丹參野生資源的收集利用，保護(hù)資源活動(dòng)和科學(xué)研究提供參考。研究中的不足之處是僅考慮了氣候和地區(qū)的影響，而沒(méi)有考慮土壤類(lèi)型和人類(lèi)活動(dòng)范圍以及社會(huì)的城鎮(zhèn)化發(fā)展對(duì)其分布的影響。

參考文獻(xiàn)

[1]

林來(lái)官.福建植物志.第4卷[M].福州：福建科學(xué)技術(shù)出版社，1989：560.

[2] 沈建芳，王強(qiáng)，汪紅，等.南丹參化學(xué)成分研究[J].中國(guó)野生植物資源，2006，25（2）：55-58.

[3] 楊慶華，周翔宇，魏宇昆，等.華東區(qū)系鼠尾草屬藥用植物資源[J].中國(guó)野生植物資源，2011，30（2）：21-26.

[4] 嚴(yán)玉平，由會(huì)玲，徐樹(shù)楠，等.丹參種源分布及道地性研究[J].時(shí)珍國(guó)醫(yī)國(guó)藥，2007，18（8）：1882-1883.

[5] 汪紅，王強(qiáng).丹參及鼠尾草屬植物的rDNAITS序列分析[J].中草藥，2005，36（9）：1381-1385.

[6] 施天慧，王強(qiáng)，袁理春，等.云南丹參資源植物的成分測(cè)定與質(zhì)量評(píng)價(jià)[J].植物資源與環(huán)境學(xué)報(bào)，2005，15（2）：73-74.

[7] 國(guó)家藥典委員會(huì).《中華人民共和國(guó)藥典》2005版，一部[S].北京：化學(xué)工業(yè)出版社，2005：52-53.

[8] 李旻輝，陳建民，彭勇，等.中國(guó)鼠尾草屬酚酸類(lèi)活性成分的分布規(guī)律研究[J].世界科學(xué)技術(shù)：中醫(yī)藥現(xiàn)代化，2008，10（5）：46-52.

[9] 舒志明，梁宗鎖，孫群，等.不同丹參種質(zhì)生物學(xué)性狀比較與評(píng)價(jià)[J].西安文理學(xué)院學(xué)報(bào)：自然科學(xué)版，2007，10（2）：24-29.

[10] 劉曉龍，王曉明.安徽鼠尾草屬藥用植物資源及開(kāi)發(fā)利用[J].中國(guó)中藥雜志，1996，21（2）：73-75.

[11] 段英姿，客紹英，曹靜，等.秋水仙堿誘導(dǎo)南丹參多倍體的研究[J].中國(guó)中藥雜志，2006，31（6）：445-448.

[12] 段英姿，牛應(yīng)澤，彭朝忠，等.南丹參離體快速繁殖與多倍體誘導(dǎo)[J].植物生理學(xué)通訊，2003，39（3）：201-205.

[13] 嚴(yán)玉平，由會(huì)玲，徐樹(shù)楠，等.丹參種源分布及道地性研究[J].時(shí)珍國(guó)醫(yī)國(guó)藥，2007，18（8）：1882-1883.

[14] PHILLIPS S J，ANDERSON R P，SCHAPIRE R E.Maximumentropymodeling of species geographic distributions[ J].Ecological Modelling，2006，190：231-259.

[15] 王運(yùn)生，謝丙炎，萬(wàn)方浩，等.ROC曲線分析在評(píng)價(jià)入侵物種分布模型中的應(yīng)用[J].生物多樣性，2007，15（4）：365-372.

[16] WISZ M S，HIJMANS R J，LI J，et al.Effects of sample size on the performance of species distribution models[J].Diversity and Distributions，2008，14（5）：763-773.

[17] SWETS J A.Measuring the accuracy of diagnostic systems[J].Science，1988，240：1285-1293.

[18] 陳新美，雷淵才，賈宏炎，等.樣本量對(duì)模型預(yù)測(cè)物種分布精度和穩(wěn)定性的影響[J].林業(yè)科學(xué)，2012，48（1）：53-58.

[19] 陳林，程登發(fā)，田喆，等.多種生物潛在分布預(yù)測(cè)分析軟件功能對(duì)比分析[C]//第一屆全國(guó)生物入侵學(xué)術(shù)研討會(huì)論文摘要集.福州，2012：53-59.

[20] 徐衛(wèi)華，羅翀.MaxEnt模型在秦嶺川金絲猴生境評(píng)價(jià)中的應(yīng)用[J].森林工程，2010，26（2）：1-3.

[21] 張海娟，陳勇，黃烈健，等.基于生態(tài)位模型的薇甘菊在中國(guó)適生區(qū)的預(yù)測(cè)[J].農(nóng)業(yè)工程學(xué)報(bào)，2011，27（5）：413-419.

[22] 常志隆，周益林，趙遵田，等.基于MaxEnt的小麥印度腥黑穗病在中國(guó)的適生性分析[J].植物保護(hù)，2010，36（3）：110-112.