基于生態(tài)位模型的南、北五味子潛在適生區(qū)預(yù)測(cè)及其影響因子研究

張迪 薛明珂 劉明 尚曉峰 張旭彤 劉金娜

摘要：基于184條南五味子（Schisandra sphenanthera）地理分布記錄、129條北五味子（Schisandra chinensis）地理分布記錄以及23個(gè)環(huán)境因子數(shù)據(jù)，采用MaxEnt生態(tài)位模型預(yù)測(cè)南、北五味子潛在適宜分布區(qū)，評(píng)估了制約其地理分布的主要環(huán)境因子。結(jié)果表明，北五味子和南五味子MaxEnt模型訓(xùn)練集和驗(yàn)證集曲線下面積（AUC）均大于0.90，建立的模型預(yù)測(cè)結(jié)果準(zhǔn)確。北五味子最佳生境主要分布在遼寧省和吉林省的長(zhǎng)白山區(qū)域以及北京市；南五味子的最佳生境主要分布在秦巴山區(qū)內(nèi)，以及湖南省武陵山脈，郴州市、株洲市和衡陽(yáng)市，浙江省金華市等地區(qū)。限制北五味子生長(zhǎng)的重要環(huán)境因子變量分別是最濕月份降水量、最熱月份最高溫度、土壤酸堿度、最暖季平均降水量、最冷季平均降水量和晝夜溫差月均值；限制南五味子分布的環(huán)境因子變量分別是最冷月份最低溫溫度、年平均降水量、溫度變化方差和降水量變化方差。南、北五味子的潛在適生區(qū)呈區(qū)域化分布。

關(guān)鍵詞：南五味子（Schisandra sphenanthera）； 北五味子（Schisandra chinensis）；MaxEnt模型； 潛在適生區(qū)；環(huán)境因子

中圖分類號(hào)：S567.5? ? ? ? ?文獻(xiàn)標(biāo)識(shí)碼：A

文章編號(hào)：0439-8114（2024）01-0116-06

DOI：10.14088/j.cnki.issn0439-8114.2024.01.021 開(kāi)放科學(xué)（資源服務(wù)）標(biāo)識(shí)碼（OSID）：

Prediction of potential distribution of Schisandra chinensis and Schisandra sphenanthera based on MaxEnt model and analysis of its important influencing factors

ZHANG Di1，2， XUE Ming-ke1， LIU Ming1， SHANG Xiao-feng1， ZHANG Xu-tong1， LIU Jin-na1

（1.Yangling Vocational & Technical College， Yangling? 712100， Shannxi，China；2.Institute of Environment and Sustainable Development in Agriculture， Chinese Academy of Agricultural Sciences， Beijing? 100081，China）

Abstract： Based on the geographical distribution records of 184 Schisandra sphenanthera， 129 Schisandra chinensis and 23 environmental factors， the maximum entropy （MaxEnt） ecological niche model was used to predict the potential suitable distribution areas of Schisandra sphenanthera and Schisandra chinensis. The main environmental factors restricting their geographical distribution were evaluated. The results showed that the area under the curve （AUC） of the training set and verification set of MaxEnt model of Schisandra chinensis and Schisandra sphenanthera were both greater than 0.90， and the prediction results of established models were accurate. The best habitats of Schisandra chinensis were mainly distributed in Changbai Mountain area of Liaoning Province and Jilin Province and Beijing City. The best habitats of Schisandra sphenanthera mainly distributed in Qinba Mountain area， Wuling Mountain Range of Hunan Province， Chenzhou City， Zhuzhou City and Hengyang City， Jinhua City of Zhejiang Province and other areas. The important environmental factors limiting the growth of Schisandra chinensis were precipitation in the wettest month， maximum temperature in the hottest month， soil pH， average precipitation in the warmest season， average precipitation in the coldest season and the monthly mean temperature difference between day and night. The environmental factors limiting the distribution of Schisandra sphenanthera were the lowest temperature in the coldest month， the average annual precipitation， the variance of temperature change and the variance of precipitation change. The potentially suitable areas of Schisandra sphenanthera and Schisandra chinensis showed regional distribution.

Key words： Schisandra chinensis； Schisandra sphenanthera； MaxEnt model； potential distribution； environmental factor

收稿日期：2022-08-05

基金項(xiàng)目：陜西省教育廳2019年度專項(xiàng)科學(xué)研究計(jì)劃項(xiàng)目（19JK0984）；楊凌職業(yè)技術(shù)學(xué)院“博士+高職生”工作室建設(shè)項(xiàng)目（BG202005； BG2022003）；楊凌職業(yè)技術(shù)學(xué)院自然科學(xué)研究基金項(xiàng)目（ZK20-73）；楊凌職業(yè)技術(shù)學(xué)院基地建設(shè)項(xiàng)目（JD20-02）

作者簡(jiǎn)介：張 迪（1988-），男，河北辛集人，講師，博士，主要從事農(nóng)業(yè)生態(tài)系統(tǒng)模擬研究，（電話）18731229918（電子信箱）zhangdi0101@sina.com；通信作者，劉金娜，女，講師，博士，主要從事中藥材栽培與質(zhì)量評(píng)價(jià)研究，（電子信箱）sasaas_lt@163.com。

五味子的藥用歷史悠久，始載于《神農(nóng)本草經(jīng)》，具有收斂固澀、益氣生津、補(bǔ)腎寧心的功效［1］。五味子有南北之分，長(zhǎng)期以來(lái)兩者均同作為五味子入藥，但是南、北五味子在傳統(tǒng)用藥上存在差異［2］。北五味子功效偏向補(bǔ)益心腎，南五味子偏于斂肺止咳［3］，2000年《中國(guó)藥典》已將南北五味子分開(kāi)收錄。北五味子為木蘭科植物五味子［Schisandra chinensis （Turcz.） Baill.］的干燥成熟果實(shí)，南五味子為木蘭科植物華中五味子［Schisandra sphenanthera Rehd. et Wils.］干燥成熟果實(shí)［1］。南、北五味子均被納入《國(guó)家重點(diǎn)保護(hù)野生藥材物種名錄》，屬于三級(jí)保護(hù)野生藥材物種。五味子的藥用價(jià)值日益受到重視，市場(chǎng)需求量激增，由于過(guò)度采摘，導(dǎo)致野生資源遭到嚴(yán)重破壞，產(chǎn)量和質(zhì)量逐年下降。因此，保護(hù)五味子野生資源，實(shí)行規(guī)范化種植成為亟待解決的問(wèn)題［4］。

南、北五味子均是重要的護(hù)肝中藥，其制劑已被廣泛應(yīng)用于臨床治療，如以南五味子為原料的五酯膠囊、五酯片、五酯軟膠囊等，以北五味子為原料的制劑五味子顆粒、五味子糖漿等［5］。雖然南、北五味子均有護(hù)肝作用，但是不同產(chǎn)地來(lái)源的五味子藥效差異較大［6，7］。此外，南、北五味子為同科屬植物，其性質(zhì)、顯微特征等相似，導(dǎo)致南、北五味子在藥品、食品、保健品等使用方面存在嚴(yán)重混用［6］。解決南、北五味子制劑在臨床上混用的問(wèn)題已經(jīng)引起廣大醫(yī)學(xué)家和研究者的共鳴［5］。正確劃分五味子區(qū)域成為解決南、北五味子混用問(wèn)題的前提。環(huán)境條件是影響植物地理空間分布以及化學(xué)成分種類和含量差異的重要因素［8］?？茖W(xué)預(yù)測(cè)南、北五味子的生態(tài)適宜區(qū)，在此基礎(chǔ)上建立野生資源保護(hù)區(qū)和人工栽培基地，這對(duì)不同五味子生境保護(hù)、合理開(kāi)發(fā)和規(guī)范使用具有重要意義。

國(guó)內(nèi)外對(duì)特定物種生境適宜性評(píng)價(jià)及潛在分布區(qū)預(yù)測(cè)的研究較多。楊航等［9］應(yīng)用主成分分析和相關(guān)分析，探索了吉林省靖宇縣北五味子單株產(chǎn)量與土壤養(yǎng)分和地理因子的關(guān)系。閆伯前等［8］分析了東北17個(gè)野生五味子種群與環(huán)境因子的關(guān)系。實(shí)地調(diào)查研究具有調(diào)研周期長(zhǎng)、成本高以及難以全面統(tǒng)計(jì)物種全部分布區(qū)的缺點(diǎn)［10］。物種分布模型由于其工作效率高、預(yù)測(cè)分布區(qū)全面、省時(shí)省力等優(yōu)點(diǎn)，已被廣泛應(yīng)用于植物生境適宜性預(yù)測(cè)研究中［10-13］。郭彥龍等［14］利用模糊數(shù)學(xué)模型預(yù)測(cè)了秦嶺地區(qū)野生華中五味子的潛在空間分布。胡理樂(lè)等［15］利用最大熵模型（MaxEnt）預(yù)測(cè)了A2和A1B未來(lái)氣候變化背景下北五味子的分布范圍。MaxEnt模型是眾多生態(tài)模型中應(yīng)用較為廣泛的模型，被廣泛應(yīng)用于入侵生物學(xué)、保護(hù)生物學(xué)、全球氣候變化對(duì)物種分布的影響等領(lǐng)域［16］。盡管關(guān)于北五味子和南五味子種植分布已有研究，但是在全國(guó)范圍內(nèi)預(yù)測(cè)南、北五味子的潛在空間分布還需進(jìn)一步研究。為此，本研究利用MaxEnt模型分別對(duì)南、北五味子生境適宜性進(jìn)行定性、定量和定位評(píng)價(jià)，明確其生境適宜性特征與空間分布差異，可為生境保護(hù)與管理提供參考，同時(shí)也為野生撫育GAP基地建設(shè)提供借鑒。

1 材料與方法

1.1 物種分布數(shù)據(jù)的獲取

南五味子和北五味子地理分布數(shù)據(jù)收集自國(guó)家標(biāo)本資源共享平臺(tái)（http：//nsii.org.cn/2017/home.php）和中國(guó)數(shù)字植物標(biāo)本館（https：//www.cvh.ac.cn/）。初步獲得南五味子標(biāo)本數(shù)據(jù)2 165條，北五味子標(biāo)本數(shù)據(jù)1 509條，為保證其樣本點(diǎn)信息的準(zhǔn)確性，利用高德開(kāi)發(fā)平臺(tái)（https：//lbs.amap.com/console/show/picker）對(duì)所選取樣本點(diǎn)核實(shí)經(jīng)緯度信息，去除采集地信息不全、重復(fù)記錄的標(biāo)本數(shù)據(jù)，最終得到南五味子184個(gè)分布點(diǎn)，北五味子129個(gè)分布點(diǎn)（圖1）。依照MaxEnt 3.4.1模型文件輸入要求，將數(shù)據(jù)整理為csv格式備用。

1.2 地理和環(huán)境數(shù)據(jù)

從WorldClim數(shù)據(jù)庫(kù)（https：//www.worldclimorg/）下載1970—2000年世界氣候環(huán)境圖層數(shù)據(jù)（Version 2.1），包括11個(gè)溫度數(shù)據(jù)和8個(gè)降雨數(shù)據(jù)（表1）。圖層分辨率為30″（約1.0 km）。從國(guó)家冰川凍土沙漠科學(xué)數(shù)據(jù)中心（http：//www.crensed.ac.cn/portal/）下載基于世界土壤數(shù)據(jù)庫(kù)（HWSD）的中國(guó)土壤數(shù)據(jù)集（V1.1）和中國(guó)（1 km）分辨率數(shù)字高程模型數(shù)據(jù)集，選取pH、有機(jī)碳含量、土壤質(zhì)地類型、海拔高度4個(gè)變量因子。

1.3 模型預(yù)測(cè)及評(píng)估

首先將南、北五味子地理分布點(diǎn)和環(huán)境變量數(shù)據(jù)轉(zhuǎn)換成ASCII格式，然后導(dǎo)入MaxEnt 3.4.1軟件中進(jìn)行建模運(yùn)算。隨機(jī)選取25%的分布點(diǎn)數(shù)據(jù)作為測(cè)試集，剩余的75%分布點(diǎn)數(shù)據(jù)作為訓(xùn)練集，設(shè)置10次重復(fù)，其他參數(shù)默認(rèn)，進(jìn)行模型準(zhǔn)確性分析［17］。利用受試者工作特征曲線（Receiver operating characteristic curve，ROC）對(duì)MaxEnt模型運(yùn)算結(jié)果進(jìn)行檢驗(yàn)。ROC曲線下面積（Area under curve，AUC）數(shù)值為0～1，當(dāng)AUC為0～0.50時(shí)，表明模型預(yù)測(cè)效果差；當(dāng)AUC為0.70～0.80時(shí)，認(rèn)為模型預(yù)測(cè)結(jié)果較準(zhǔn)確；當(dāng)AUC為0.80～0.90時(shí)，認(rèn)為模型預(yù)測(cè)結(jié)果很準(zhǔn)確；當(dāng)AUC大于0.90時(shí)，認(rèn)為模型預(yù)測(cè)結(jié)果極準(zhǔn)確，AUC越接近1，表明模型預(yù)測(cè)結(jié)果越好［18］。

1.4 適生等級(jí)劃分和主導(dǎo)環(huán)境因子分析

將MaxEnt軟件分析得到的結(jié)果導(dǎo)入Arcgis10.2軟件中，利用Arcgis 10.2軟件空間分析工具中重分類工具進(jìn)行適生區(qū)劃分。參考Yuan等［19］的方法，將五味子適宜性等級(jí)分為5級(jí)：不適宜生境（0～0.2）、低適宜生境（0.2～0.4）、中適宜生境（0.4～0.6）、高適宜生境（0.6～0.7）、最高適宜生境（0.7～1.0）。

2 結(jié)果與分析

2.1 MaxEnt模型的準(zhǔn)確性評(píng)價(jià)

根據(jù)MaxEnt模型模擬的預(yù)測(cè)結(jié)果，北五味子和南五味子ROC曲線訓(xùn)練集的AUC分別為0.913和0.952，測(cè)試集的AUC分別為0.907和0.944，均大于0.90。表明本研究建立的MaxEnt模型能夠精確模擬北五味子和南五味子的地理分布（圖2）。

2.2 北五味子和南五味子在中國(guó)的生態(tài)適宜性區(qū)劃

南、北五味子生態(tài)適宜性區(qū)劃見(jiàn)圖3。從空間格局來(lái)看，北五味子最佳生境主要分布在遼寧省和吉林省的長(zhǎng)白山區(qū)域以及北京市，中高適宜區(qū)主要分布在陜西省、山西省、湖北省和湖南省等。南五味子的最佳生境主要分布在秦巴山區(qū)，包括河南省、陜西省和湖北省交界處，甘肅省、四川省和陜西省交界處，重慶市北部地區(qū)。湖南省武陵山脈以及郴州市、株洲市和衡陽(yáng)市，浙江省金華市的南五味子生態(tài)適宜性也較高。

2.3 影響北五味子和南五味子分布的關(guān)鍵環(huán)境因子

利用MaxEnt模型［20］運(yùn)行得出的貢獻(xiàn)率和置換重要值對(duì)北五味子和南五味子在中國(guó)潛在地理分布的環(huán)境變量重要性進(jìn)行評(píng)估（表2）。用于MaxEnt模型預(yù)測(cè)的23個(gè)環(huán)境變量中，對(duì)北五味子分布貢獻(xiàn)率排名前三位的環(huán)境因子變量分別為最濕月份降水量（39.2%）、最熱月份最高溫度（12.7%）和土壤酸堿度（9.8%）；對(duì)南五味子分布貢獻(xiàn)率排名前三位的環(huán)境因子分別為最冷月份最低溫度（40.7%）、年平均降水量（28.3%）和溫度變化方差（6.8%）。

置換重要值是隨機(jī)置換每個(gè)環(huán)境變量在訓(xùn)練存在和背景數(shù)據(jù)上的數(shù)值，其數(shù)值越大表明模型對(duì)該變量的依賴性越強(qiáng)［21］。北五味子置換重要值排名前三位的環(huán)境因子變量分別為最暖季平均降水量（14.4%）、最冷季平均降水量（10.5%）和晝夜溫差月均值（9.2%）；南五味子分別為溫度變化方差（18.0%）、降水量變化方差（17.0%）和最冷月份最低溫（13.6%）。

利用環(huán)境因子變量響應(yīng)曲線判斷五味子的存在概率與環(huán)境因子之間的關(guān)系，當(dāng)五味子的存在概率大于0.5時(shí)，表明其對(duì)應(yīng)的環(huán)境變量值有利于五味子的生長(zhǎng)。由圖4可以看出，適宜北五味子生長(zhǎng)的最熱月份最高溫度的范圍為25～31 ℃，在溫度達(dá)27.5 ℃時(shí)北五味子的生存概率達(dá)最大值。適宜北五味子生長(zhǎng)的最濕月份降水量的范圍為135～280 mm，當(dāng)降水量達(dá)200 mm時(shí)北五味子的生存概率達(dá)最大值（0.56）。適宜北五味子生長(zhǎng)的最暖季平均降水量和最冷季平均降水量的范圍分別為331～688 mm和9～116 mm。

南五味子存在概率對(duì)主要?dú)夂蛞蜃拥捻憫?yīng)曲線見(jiàn)圖5。由圖5可以看出，年平均降水量在765～? ? ? 1 690 mm適宜南五味子的分布及生長(zhǎng)；最冷月份最低溫度在-3.1～3.7 ℃適宜南五味子的分布及生長(zhǎng)；溫度變化方差和降水變化方差分別在706～890 ℃和44～72 mm時(shí)適宜南五味子的分布及生長(zhǎng)。

3 討論

MaxEnt是以氣候相似性為原理，利用已有的物種分布數(shù)據(jù)和環(huán)境數(shù)據(jù)找出物種概率分布的最大熵，從而對(duì)物種生境適宜性進(jìn)行評(píng)價(jià)［20］。建立的生態(tài)位模型需要通過(guò)驗(yàn)證后才能應(yīng)用于預(yù)測(cè)物種潛在分布的研究中。本研究隨機(jī)選取75%分布點(diǎn)數(shù)據(jù)作為訓(xùn)練集，建立MaxEnt模型參數(shù)，剩余25%的分布點(diǎn)數(shù)據(jù)作為測(cè)試集，對(duì)已經(jīng)建立的模型進(jìn)行驗(yàn)證，利用ROC曲線評(píng)估MaxEnt模型模擬的準(zhǔn)確性。結(jié)果顯示，南、北五味子的訓(xùn)練集和驗(yàn)證集AUC均在0.90以上，表明模型模擬效果好，可信度高，可用于五味子在中國(guó)地理分布與環(huán)境關(guān)系的研究。

本研究表明，北五味子最適宜種植區(qū)域集中在遼寧省、吉林省和北京市，在陜西省、山西省、湖北省和湖南省等也有較高適宜種植區(qū)。本研究得出的北五味子分布與《中國(guó)高等植物圖鑒》的北五味子空間分布一致，與胡理樂(lè)等［15］的研究結(jié)果略有差異，可能是適宜性等級(jí)劃分標(biāo)準(zhǔn)不一致造成的。南五味子的最適宜種植區(qū)域集中在甘肅省、陜西省、河南省、湖北省、四川省、重慶市，湖南省和浙江省等也有較高適宜種植區(qū)，本研究得出的南五味子分布與《中國(guó)高等植物圖鑒》的南五味子空間分布一致。最大熵模型潛在生境分析方法是利用環(huán)境數(shù)據(jù)預(yù)測(cè)植物潛在最佳生長(zhǎng)區(qū)域，而并未將植物有效成分含量加入到模型中。本研究得出的南五味子分布與郭彥龍等［14］基于五味子甲素和乙素的野生南五味子潛在空間分布一致。一般認(rèn)為，植物生長(zhǎng)在最適宜的地區(qū)，其體內(nèi)的主要內(nèi)含物含量也較高［22-25］。如何利用最大熵模型（MaxEnt）適宜分析方法預(yù)測(cè)植物內(nèi)含物最佳含量潛在分布，還有待進(jìn)一步研究。

楊航等［9］研究發(fā)現(xiàn)北五味子的基因型與土壤養(yǎng)分呈顯著正相關(guān)，與海拔高度呈負(fù)相關(guān)，而與經(jīng)緯度不相關(guān)。閆伯前等［8］報(bào)道，年均氣溫是影響北五味子果實(shí)質(zhì)量的最主要因子，然后依次是土壤全磷、全氮，有機(jī)質(zhì)含量和年降水量。胡理樂(lè)等［15］研究認(rèn)為溫度季節(jié)性變化、最濕月份降水量、最熱季降水量、年均溫、最冷月最低溫度和最熱月最高溫度是限制北五味子分布的主要因子。毛亞娟等［26］指出最熱月的最高氣溫、最冷月的最低氣溫、最熱季降水量、年均溫、年降水量是影響北五味子果實(shí)中總木質(zhì)素含量的重要因子。本研究根據(jù)各環(huán)境變量對(duì)建模的貢獻(xiàn)率和置換重要性，限制北五味子生長(zhǎng)的重要環(huán)境因子變量分別是最濕月份降水量、最熱月份最高溫度、土壤酸堿度、最暖季平均降水量、最冷季平均降水量和晝夜溫差月均值。北五味子生長(zhǎng)的最熱月份最高溫度適宜范圍為25～31 ℃，北五味子生長(zhǎng)的最暖季平均降水量和最冷季平均降水量的范圍分別為331～688 mm和9～116 mm，與毛亞娟等［26］的研究結(jié)果一致（最熱月份最高溫度為25.89～28.51 ℃，最暖季平均降水量為368～654 mm），說(shuō)明北五味子適宜生長(zhǎng)在濕潤(rùn)涼爽的氣候環(huán)境中。限制南五味子分布的環(huán)境因子變量分別是最暖季平均降水量、最冷季平均降水量、晝夜溫差月均值、溫度變化方差、降水量變化方差和最冷月份最低溫。適宜南五味子生長(zhǎng)的年平均降水量在765～1 690 mm；最冷月份最低溫度為-3.1～3.7 ℃。說(shuō)明南五味子不耐寒，適宜生長(zhǎng)在潮濕溫暖的氣候環(huán)境中。

本研究在全國(guó)尺度上分別分析了南、北五味子的潛在適生區(qū)，研究結(jié)果可為南、北五味子的規(guī)范種植和野生資源保護(hù)提供一定的科學(xué)依據(jù)。本研究由于數(shù)據(jù)有限，僅考慮了溫度、降水、海拔高度和部分土壤理化指標(biāo)對(duì)其適生分布的影響，并未考慮植物遺傳特性、復(fù)雜地理環(huán)境下的小氣候條件、土壤養(yǎng)分、種群特征、人類活動(dòng)、社會(huì)經(jīng)濟(jì)結(jié)構(gòu)、栽培管理水平等因素［27-29］。在未來(lái)的工作中，應(yīng)盡可能獲取全面而準(zhǔn)確的分布數(shù)據(jù)、氣象數(shù)據(jù)、地理環(huán)境數(shù)據(jù)、地理經(jīng)濟(jì)數(shù)據(jù)等，將這些數(shù)據(jù)納入到模型中以進(jìn)一步提高模型預(yù)測(cè)的準(zhǔn)確性。

參考文獻(xiàn)：

［1］ 國(guó)家藥典委員會(huì).中華人民共和國(guó)藥典（2020年版）［M］.北京：中國(guó)醫(yī)藥科技出版社，2020.

［2］ 焦美鈺，佟 月，吳浩善，等.南、北五味子“辨色論質(zhì)”機(jī)制及其定性識(shí)別模型研究［J］.中國(guó)藥房，2020，31（24）： 3007-3012.

［3］ 劉宇靈，付 賽，樊麗姣，等.南北五味子化學(xué)成分、藥理作用等方面差異的研究進(jìn)展［J］.中國(guó)實(shí)驗(yàn)方劑學(xué)雜志，2017，23（12）：228-234.

［4］ 包京姍，楊世海，傅金泉.北五味子種質(zhì)資源及栽培技術(shù)研究進(jìn)展［J］.人參研究，2017，29（4）：26-30.

［5］ 馮 格，翟健秀，陳萬(wàn)生，等.南北五味子及其制劑與其他藥物相互作用的研究進(jìn)展［J］.中國(guó)醫(yī)院藥學(xué)雜志， 2017，37（21）：2206-2209，2215.

［6］ 牟倩倩，賀敬霞，盧 一，等.電子鼻結(jié)合化學(xué)計(jì)量學(xué)方法在南五味子和北五味子鑒別上的應(yīng)用研究［J］.時(shí)珍國(guó)醫(yī)國(guó)藥，2017，? 28（6）：1356-1359.

［7］ 左曉彬，孫小茗，王春宇，等.南、北五味子醇提物對(duì)對(duì)乙酰氨基酚誘導(dǎo)的藥物性肝損傷保護(hù)作用及藥效差異研究［J］.中國(guó)中藥雜志，2019，44（6）：1238-1245.

［8］ 閆伯前，羅遵蘭，胡理樂(lè).東北野生五味子種群特征及其與環(huán)境因子的關(guān)系［J］.環(huán)境科學(xué)研究，2013，26（6）： 618-623.

［9］ 楊 航，張麗鵬，孫 偉，等.北五味子品種適宜栽培地及其與環(huán)境因子的關(guān)系［J］.西北農(nóng)業(yè)學(xué)報(bào)，2014，23（4）： 92-98.

［10］ 張 童，黃治昊，彭?xiàng)罹福?基于Maxent模型的軟棗獼猴桃在中國(guó)潛在適生區(qū)預(yù)測(cè)［J］.生態(tài)學(xué)報(bào)，2020，40（14）：4921-4928.

［11］ 孫李勇，蔣 政，劉晨妮，等.基于DIVA-GIS的紫玉蘭地理分布及適生性分析［J］.植物科學(xué)學(xué)報(bào)，2018，36（6）： 804-811.

［12］ 段義忠，王 馳，王海濤，等.不同氣候條件下沙冬青屬植物在我國(guó)的潛在分布——基于生態(tài)位模型預(yù)測(cè) ［J］.生態(tài)學(xué)報(bào)，2020，40（21）：7668-7680.

［13］ 王國(guó)崢，耿其芳，肖孟陽(yáng)，等.基于4種生態(tài)位模型的金錢松潛在適生區(qū)預(yù)測(cè)［J］.生態(tài)學(xué)報(bào)，2020，40（17）：6096-6104.

［14］ 郭彥龍，顧 蔚，路春燕，等.基于五味子甲素與乙素的秦嶺地區(qū)野生華中五味子的潛在空間分布［J］.植物學(xué)報(bào)，2013，48（4）：411-422.

［15］ 胡理樂(lè)，張海英，秦 嶺，等.中國(guó)五味子分布范圍及氣候變化影響預(yù)測(cè)［J］.應(yīng)用生態(tài)學(xué)報(bào)，2012，23（9）： 2445-2450.

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

［17］ WU B C，ZHOU L J，QI S，et al. Effect of habitat factors on the understory plant diversity of Platycladus orientalis plantations in Beijing mountainous areas based on MaxEnt model［J］.Ecological indicators，2021， 129：107917.

［18］ WANG G，WANG C，GUO Z R，et al. Integrating Maxent model and landscape ecology theory for studying spatiotemporal dynamics of habitat： Suggestions for conservation of endangered Red-crowned crane［J］.Ecological indicators，2020，116：106472.

［19］ YUAN H-S，YULIAN W，WANG X. Maxent modeling for predicting the potential distribution of Sanghuang， an important group of medicinal fungi in China［J］.Fungal ecology，2015，17：140-145.

［20］ 張 華，趙浩翔，王 浩.基于Maxent模型的未來(lái)氣候變化情景下胡楊在中國(guó)的潛在地理分布［J］.生態(tài)學(xué)報(bào)，2020，40（18）：6552-6563.

［21］ FENG L，WANG H，MA X，et al. Modeling the current land suitability and future dynamics of global soybean cultivation under climate change scenarios［J］.Field crops research，2021，263：108069.

［22］ 賴長(zhǎng)江生，魏旭雅，邱子棟，等.道地藥材質(zhì)量標(biāo)準(zhǔn)設(shè)計(jì)與應(yīng)用［J］.中國(guó)中藥雜志，2020，45（24）：6072-6080.

［23］ 何冬梅，王 海，陳金龍，等.中藥微生態(tài)與中藥道地性［J］.中國(guó)中藥雜志，2020，45（2）：290-302.

［24］ 雷 捷，徐建國(guó)，張曉佳，等.銀柴胡的道地性評(píng)價(jià)初探［J］.時(shí)珍國(guó)醫(yī)國(guó)藥，2020，31（8）：1973-1975.

［25］ 孟祥才，沈 瑩，杜虹韋.道地藥材概念及其使用規(guī)范的探討［J］.中草藥，2019，50（24）：6135-6141.

［26］ 毛亞娟，衛(wèi)海燕，尚忠慧，等.基于GIS與模糊物元模型的東北地區(qū)五味子生境適宜性評(píng)價(jià)［J］.應(yīng)用與環(huán)境生物學(xué)報(bào)，2016，22（2）：241-248.

［27］ 劉麗娟，秦佳梅，叢 毓，等.栽培北五味子土壤營(yíng)養(yǎng)因子變化規(guī)律研究［J］.廣東農(nóng)業(yè)科學(xué)，2013，40（22）： 79-81.

［28］ 曹 倩，高慶波，郭萬(wàn)軍，等.基于MaxEnt模擬人類活動(dòng)與環(huán)境因子對(duì)青藏高原特有植物祁連獐牙菜潛在分布的影響［J］.植物科學(xué)學(xué)報(bào)，2021，39（1）：22-31.

［29］ 張正海，李愛(ài)民，張 悅，等.修剪強(qiáng)度對(duì)北五味子新梢萌發(fā)和生長(zhǎng)的影響［J］.江蘇農(nóng)業(yè)科學(xué)，2014，42（9）： 220-222.