基于MaxEnt模型的胖大海全球生態適宜區預測及生態特征研究

王梅 謝彩香

摘要 利用35個環境因子和53個胖大海（Sterculia lychnophora Semen）樣本點數據，采用Maxent模型預測胖大海全球潛在生態適宜區。用受試者工作特性曲線（ROC，receive operating characteristic）對模型進行評價。模型結果顯示，胖大海全球生態適宜區主要分布于赤道附近的熱帶雨林氣候國家，遍布于亞洲東南部、美洲中南部和非洲中西部；此外，在印度半島南部和非洲中東部也有部分存在。模型評價結果顯示，測試集和訓練集的AUC值分別為0.988和0.996，表明該模型的預測效果良好。結果顯示，影響胖大海分布的主要環境因子包括年溫度變化范圍（bio7）、最干月份降水量（bio14）、年平均降水量（bio12）、土壤電導率（T-ECE）、晝夜溫差與年溫差比值（bio3）及土壤鹽基飽和度（T-BS）。該研究結果為胖大海人工引種栽培和提高其品質質量提供參考。

關鍵詞 胖大海；最大信息熵模型；生態適宜區；生態特性

中圖分類號 S-058? 文獻標識碼 A

文章編號 0517-6611（2023）14-0218-07

作者簡介

王梅（1971—），女，山東武城人，副教授，博士，從事地理信息系統及其應用研究。*通信作者，副研究員，從事藥用植物地理分布及產地適宜性研究。

胖大海（Sterculia lychnophora Semen）為梧桐科（Sterculiaceae）蘋婆屬（ Sterculia? L.），胖大海的成熟種子又稱大海子、通大海、大洞果等，是傳統常用中藥，味甘性寒，具有清熱潤肺，利咽開音，潤腸通便的功能［1］。胖大海傳入中國的時間已不可考，可追溯的最早記載是清代趙學敏所著《本草綱目拾遺》［3-4］。現代藥理研究表明，胖大海在抑制病毒、緩瀉、鎮痛、抑菌、抗炎、增強胸腺和脾臟功能等方面都取得的一定進展［5-10］。

胖大海原產于越南、緬甸、柬埔寨、老撾、馬來西亞、泰國和印度尼西亞等東南亞和東亞國家的熱帶森林生態系統［2］，海南、云南、廣東、廣西、福建等省區也有少量栽培。性喜熱，要求年平均溫度為21.0～24.9? ℃，成齡樹較耐旱，適宜生長的月降雨量為100～300 mm。胖大海是喜陽植物，因莖稈細高，木質部松脆，根系不發達，抗風能力極差。

近年來，對胖大海樹木的不可再生性砍伐和對胖大海果實的過度利用造成了胖大海資源匱乏，甚至導致了個別地區的生態不平衡 ［11-12］。盡管胖大海在云南、海南、廣東等省已經有引種栽培的歷史，但從生態地理角度判斷這些地區是否真正適合胖大海的生產，在其他地區是否還有新的適宜區可以發展胖大海種植，這些問題都需要用科學數據來回答［13-14］。因此，擴大胖大海的適宜栽培面積，對其進行科學引種和栽培區劃、生產高質高產的胖大海是胖大海產業發展的重要方向。氣候是決定物種分布的主要因素之一，植物潛在的分布區域與環境因子密切相關［15］。在不考慮生態環境的情況下盲目引種會造成不必要的經濟損失。為避免此類情況的發生，提高胖大海的產量和質量，該研究結合GIS技術與Maxent最大信息熵模型，對胖大海的全球適宜性區劃及其生態特征進行研究，分析胖大海生長的最適宜區域和影響其生長的主要生態因子，為胖大海的科學引種、栽培、區劃提供科學依據。

1 材料

1.1 胖大海樣本分布數據

根據全球生物多樣性信息網絡（GBIF，http：//www.gbif.org/）、中國植物圖像庫（PPBC，http：//ppbc.iplant.cn/）、英國邱園（http：//www.kew.org/）等物種資源信息共享平臺以及文獻檢索，為保證樣本點信息的準確，確保符合物種實際生境，對重復的樣本點進行去除，同時通過“谷歌地球”對文獻及網絡共享資源中記載有詳細地點但丟失經緯度的樣本點進行精準定位，從而得到該樣本點的經緯度信息［16］。經過篩選后共得到53個胖大海全球分布樣本點，主要分布在亞洲赤道附近（圖1）。

1.2 環境數據

共選取環境因子35個（表1），其中氣候因子19個，土壤因子15個，地形因子1個。氣候數據來源于世界氣候信息網worldclim（http：//www.worldclim.org/），選取圖層的分辨率為30”（1 km×1 km），是1950—2000年世界氣候監測數據的平均值，主要包括19個綜合氣候因子和月均氣候因子（月均氣溫、月均降水），土壤數據來源于世界土壤數據庫（http：//www.fao.org/soils-portal/soils-survey/soilmaps-and-databases? /harmonized-World-soil-database-V12-en），坐標系為WGS84，柵格大小約為1 km2。地圖數據來自于國家基礎地理信息系統（http：//nfgis.nsdi.gov.cn/）。圖層中的溫度數值（ ℃）為實際數值×10。

1.3 最大信息熵模型及相關軟件來源

最大信息熵模型（MaxEnt）是基于最大熵原理所提出的一種生態位模型，可利用物種分布點的環境因子作為約束條件，找尋最大熵的概率分布作為物種的可能分布，從而推測出物種的適宜區［17］。該研究使用3.3.3 k版本的最大信息熵模型MaxEnt軟件來源于官網（http：//www.cs.princeton.edu/～schapire/maxent/），Arcgis10.0軟件來源于Ersi中國（http：//support.esrichina.com.cn/support/download/），Graphpad軟件同樣源于官網（https：//www.graphpad.com/）。

2 方法

2.1 MaxEnt模型參數設置

在MaxEnt軟件中設置參數，將胖大海樣本分布點數據（csv格式）與氣候因子數據集（.asc）導入并運行。參數具體設置為：測試集為樣本分布數據隨機抽取的25%，訓練集為樣本分布點數據的75%，采用自然分類法檢測權重，設置最大迭代次數為500，文件輸出格式選為.asc。MaxEnt自定義設置受試者工作特性曲線（ROC，receive operating characteristic）對評測模型預測結果精度，其余參數選項采用軟件默認設置。

2.2 胖大海生態適宜性區域預測

將MaxEnt模型分析計算結果以ASC II類型文件輸出，導入Arcgis 10.0軟件中轉化為Raster格式，利用Arcgis 10.0軟件進行疊加分析和地圖制作，繪制出基于氣候因子的胖大海生態適宜性區劃圖，采用空間分析工具（Spatial Analyst Tools）中的自然間斷點分級法（Jenksnatural breaks）對胖大海進行生態適宜性等級分類。

3 結果與分析

3.1 預測區劃結果評價

ROC曲線與橫坐標所圍成的面積叫做AUC值，其大小一般作為預測模型準確度的衡量指標。AUC值的數值范圍為0～1，數值越大，預測模型模擬的精確度則越高，預測效果越好。一般認為AUC值大于0.9時預測模擬效果的價值較高；0.7～0.9時的預測模擬效果的價值中等；0.5～0.7時的預測模擬效果的價值較低。胖大海分布預測的ROC曲線如圖2所示，訓練集AUC值為0.996，測試集AUC值為0.988，說明預測模擬效果達到的水平較高，由此模型運算得出的胖大海生態適宜性區劃具有較高的準確度和可信度。

3.2 氣候因子貢獻率 由表3可知，Maxent模型運算得出的35個環境因子貢獻率中，列出了貢獻率大于0.1的18個環境因子及其貢獻率。選擇前6個因子作為影響胖大海分布的主要因子，6個環境因子總貢獻率超過90%（累計總貢獻率達到90.75%），共包括2個土壤因子和4個氣候因子，其中土壤因子的總貢獻率為15.82%，氣候因子的總貢獻率74.93%。

由表3可知，對胖大海資源分布影響最大的2個環境因子分為年溫度變化范圍（bio7）和最干月份降水量（bio14），貢獻率分別為41.02%和12.08%。主要環境因子貢獻度從大到小依次為年溫度變化范圍（bio7）、最干月份降水量（bio14）、年平均降水量（bio12）、土壤電導率（T-ECE）、晝夜溫差與年溫差比值（bio3）及土壤鹽基飽和度（T-BS）。

3.3 胖大海主要氣候因子適宜值范圍

Maxent軟件模型計算得出各主要生態因子響應曲線存在概率與生態因子閾值范圍關系，如圖3所示（圖中溫度數值為實際數值×10）。一般認為，植物的環境因子適宜值范圍，對應環境因子響應曲線存在概率（probability of presence）應大于0.5［18］。根據胖大海的響應曲線可得出各主要環境因子響應曲線存在概率大于0.5的范圍，即各主要環境因子的適宜值范圍為：年溫度變化范圍（貢獻率41.02%）適宜值范圍為0～10.1? ℃，存在概率最大值在7.8? ℃左右，且大于7.8 ℃時存在概率明顯下降；最干月份降水量（貢獻率12.08%）適宜值范圍為103～2 56 mm，存在概率最大值存在于156 mm左右；年平均降水量（貢獻率10.98%）適宜值范圍為2 360～5 440 mm，存在概率最大值位于3 070 mm左右；晝夜溫差與年溫差比值（貢獻率10.85%）適宜值范圍為0.83～0.92，在比值為0.88時存在概率最大，并且晝夜溫差與年溫差比值越小，存在概率越小。

3.4 胖大海潛在生態適宜區預測

Maxent模型預測結果說明，胖大海潛在生態適宜區主要分布在赤道周圍18°S～22°N的熱帶雨林區域（圖4A），主要集中在東南亞地區、印度半島南部、大洋洲北部、非洲幾內亞灣沿岸及剛果雨林、南美洲北部及北美洲加勒比海沿岸。

亞洲（圖4B）的最適宜區主要分布在印度、越南、馬來西亞、新加坡、印度尼西亞、菲律賓等國；較適宜區主要分布在柬埔寨、緬甸，泰國等；次適宜區主要分布在中國、孟加拉國、老撾等。大洋洲（圖4B）的最適宜區主要分布在巴布亞新幾內亞；較適宜區主要分布在澳大利亞北部；次適宜區主要分布在斐濟、瓦努阿圖等。分布在非洲（圖4C）的最適宜區主要在塞拉利昂、尼日利亞、喀麥隆等國家；較適宜區主要在加納、剛果、幾內亞、剛果共和國等國家；次適宜區主要在安哥拉、坦桑尼亞等國家。分布在美洲（圖4D）的最適宜區主要在巴拿馬、哥倫比亞、委內瑞拉、巴西等國家；較適宜區主要在危地馬拉、厄瓜多爾、秘魯、牙買加等國家；次適宜區主要在墨西哥、巴哈馬等國家。

根據胖大海各個適宜等級在不同地區的分布情況（圖5）可知，胖大海最適宜區涉及的國家數目最多的是亞洲，其次是美洲、非洲和大洋洲；較適宜區覆蓋的國家數最多的是美洲，其次是非洲、亞洲和大洋洲；次適宜區覆蓋國家數最多的是美洲、其次是非洲，亞洲和大洋洲。

4 不同適宜區差異性比較

4.1 不同生態適宜區主要氣候因子分析

從表3可以看出，對胖大海資源分布影響最大的5類氣候因子利用ArcGIS軟件，在所得各個適宜區中進行取點，保證點位的隨機性和全面性，盡可能囊括各個適宜區中所有可能不同的生態特征，并提取所選擇點位的主要環境因子進行分析。各適宜區所提取點位的環境因子范圍及平均值見表4。

4.2 不同適宜區主要氣候因子差異分析 從表4可以看出，最適宜區的年溫度變化范圍相較于較適宜區和次適宜區明顯偏小，即1年中溫度變化幅度較小，年溫度保持相對穩定，而次適宜區溫度變化幅度明顯較大，較不適宜胖大海的生長。最干月份降水量方面，最適宜區雖也存在降水偏少的情況，但整體平均值較適宜區和次適宜區明顯高，且不存在無有效降水的月份，即不存在最干月份降水量為0的情況。各適宜區年平均降水量雖差異不明顯，但最適宜區和次適宜區的年平均降水量較為充沛，而最適宜區不存在降水量過大的情況，說明過澇或年均降水不足均不適宜胖大海的生長。最適宜區和較適宜區在晝夜溫差與年溫差比值方面相差較小，而次適宜區與兩者相比則具有明顯差異，次適宜區的晝夜溫差與年溫差比值明顯偏低，即說明了胖大海適宜在1年中溫度變化差距不大且晝夜溫差也較小的地區生長。

5 討論

5.1 胖大海生態特征

胖大海原產于熱帶，喜高溫多雨的氣候，忌低溫冷害，成齡樹木較耐旱。該研究結果表明，在影響胖大海地理分布的6個主要環境因子中，4個氣候因子的貢獻率總和大于2個土壤因子的貢獻率總和。年溫度變化范圍和最干月份降水量是對胖大海適宜區分布影響最大的2個因子，由上述圖表可知，適宜區的年溫度變化范圍適宜值范圍為0～10.1 ℃，最干月份降水量為103～256 mm。

由此可知，胖大海是適宜生長在常年溫度在20 ℃以上，溫差變化較小，干旱月份仍有有效降水的地區，年溫度變化在15 ℃以上時胖大海的存在概率幾乎為0，這些都說明胖大海適宜在熱帶地區生長，但不適宜在旱季較長的地區種植。

除溫度變化方差和最干月份降水量2個最主要環境因子外，年平均降水量及晝夜溫差與年溫差比值等次要因子也說明，胖大海的生態適宜區應在溫度較高、晝夜溫差較小且最暖季降水充沛的地區，這些都與胖大海喜高溫多雨忌冷害的生物學特性相符。

大量文獻表明，溫度是影響植物開花結果的主要因素［19］。巫金華等［20］發現，在胖大海樹花芽發育及果實成熟等多個階段對溫度都有較高的要求，3月份胖大海開花結果時，要求氣溫必須高于22 ℃，升溫較早的年份里，開花結果時間也較早。這些研究更說明了在胖大海生長發育的各個過程中溫度都起著十分重要的作用。除溫度與降水外，由于胖大海莖稈細高、木質部松脆、根系不發達，導致其樹木抗風能力極差，在臺風頻繁地區種植時極易造成斷稈或倒株。同時，胖大海屬喜陽植物，張巧巧等［21］研究表明，過強或過弱的光照都會對胖大海的生長起抑制作用。

綜上所述，對于氣候因子而言，胖大海適宜生長在高溫、多雨、向陽、少臺風的地區，風力較大時會導致胖大海樹木折斷，溫度稍低則會導致胖大海停止生長。而在胖大海的實際栽培和引種過程中，除了受氣候因子的影響較大，還有很多其他因素也影響著胖大海的生長發育過程［22-25］。

5.2 胖大海全球潛在生態適宜區預測

模型預測結果顯示，胖大海的全球潛在生態適宜區主要位于熱帶雨林氣候帶。熱帶雨林氣候分布于赤道兩側南北緯10°之間。具體說來，主要分布在南美洲的亞馬孫河流域，中美東部，非洲的剛果河流域、幾內亞灣北岸，馬達加斯加島東部，馬來群島，菲律賓群島南部，印度尼西亞，大洋洲的蘇門答臘島至新幾內亞島一帶、澳大利亞的東北部等地。熱帶雨林氣候全年高溫多雨、氣溫年變化較小，符合胖大海的生長需求。

根據模型預測，胖大海在中南美洲的潛在最適宜區主要位于亞馬孫河流域及加勒比海沿岸地區，亞洲地區的馬來群島，菲律賓群島南部，印度尼西亞等東南亞國家以及位于幾內亞灣北部地區的尼日利亞，喀麥隆等非洲西部國家。這些地區都是極為顯著的熱帶雨林氣候地區，氣候條件適合胖大海的生長。而在中南美洲，所得潛在生態適宜區與張婕等［14］所得結果存在較大差異，分析原因可能在于其所選胖大海空間分布樣本點多為國內中藥材植物野生分布點，氣候等因子與胖大海原產地馬來西亞等地存在一定差異，從而導致預測結果有一定差別。

同時，預測結果顯示，印度南部地區、越南及菲律賓境內存在部分胖大海生態最適宜區，這些地區屬于熱帶季風氣候，旱季雨季分明且旱季較長，與胖大海的生態特征有一定差異，但也具有高溫多雨的特點，較適宜胖大海的生長。同樣屬于熱帶季風氣候的地區還有中國海南、泰國和老撾等地。而墨西哥南部、坦桑尼亞及馬達加斯加島等熱帶草原氣候地區有生態次適宜區分布，但文獻及標本館中未出現相應的樣本記載，因此可考慮作為胖大海引種栽培的試驗區域。

5.3 胖大海在中國的潛在生態適宜區

由胖大海潛在生態適宜區的模型預測結果可知，在中國臺灣北部、廣東南部、廣西南部、海南島、云南西雙版納及福建南部等地存在一定面積的生態次適宜區。這些地區均屬于熱帶季風氣候，全年高溫，降雨量大，但旱季也較為明顯。對比這些地區的氣候條件與胖大海原產地的氣候條件可知，海南、廣東等地的年平均溫度較低且最干月份有效降水不足，對胖大海的生長發育存在一定影響。由于海南等地常年遭受臺風襲擊，極易造成胖大海斷稈或倒株，所以可能導致在此地區大面積種植的難度較大。因此，在分析結果的基礎上，各地需結合當地自然條件及社會經濟情況等方面，進行適當引種栽培。

參考文獻

［1］國家藥典委員會.中華人民共和國藥典：2015 年版一部［S］.北京：中國醫藥科技出版社，2015：261.

［2］ LE M T，LE T H.Vietnam plants forest［M］.Hanoi：Hanoi Agriculture Press，2000：97-103.

［3］ 徐祥浩，徐頌軍.中藥胖大海名實考［J］.中藥材，1987（2）：52-53.

［4］ 徐祥浩.中藥胖大海的引進歷史和拉丁學名考證［J］.華南農學院學報，1982，3（3）：100-104.

［5］ 安忠蘭，王尚德，任保成.胖大海治療腹瀉 560 例臨床觀察［J］.中醫藥研究，1994（5）：12.

［6］ 啟明.胖大海治便秘［J］.家庭中醫藥，1994（4）：45.

［7］ 周長堅，林鍵，許凌枚，等.復方胖大海的藥理研究［J］.福建中醫學院學報，1994，4（ 3） ：30-33.

［8］ 杜力軍，孫紹美，於蘭，等.國產與進口胖大海對小鼠抗炎和小腸推進作用比較［J］.中藥材，1995，18（8）：409-411.

［9］ OPPONG M B，LI Y，BANAHENE P O，et al.Ethnopharmacology，phytochemistry，and pharmacology of Sterculia lychnophora Hance （Pangdahai）［J］.Chinese journal of natural medicines，2018，16（10）：721-731.

［10］ 李娜，高昂，鞏江，等.胖大海藥學研究概況［J］.安徽農業科學，2011，39（16）：9609-9610.

［11］ THAI V T.The rainforest ecosystem of Vietnam［M］.Hanoi：Science Technology Press，1998：321-345.

［12］ THAI V M.Technical and economic feasibility to develop of reserve phase 2015［M］.Dongnai：Dong Nai Natural and Culture Reserve，2010.

［13］ 王世敏，吳德勛，張仲偉，等.胖大海在海南島生長狀況觀察［J］.中藥材，1998，21（6）：271-274.

［14］ 張婕，李西文，許翔翔，等.胖大海全球產地適宜性分析［J］.世界科學技術-中醫藥現代化，2016，18（8）：1280-1288.

［15］ 邵璞，曾曉東.氣候年際變率對全球植被平均分布的影響［J］.生態學報，2011，31（6）：1494-1505.

［16］ WILKIE P.A revision of Scaphium（Sterculioideae，Malvaceae/Sterculiaceae）［J］.Edinburgh journal of botany，2009，66（2）：283-328.

［17］ 郭杰，劉小平，張琴，等.基于Maxent模型的黨參全球潛在分布區預測［J］.應用生態學報，2017，28（3）：992-1000.

［18］ 張東方，張琴，郭杰，等.基于MaxEnt模型的當歸全球生態適宜區和生態特征研究［J］.生態學報，2017，37（15）：5111-5120.

［19］ 朱云美，張玲，高用順，等.溫度對植物開花時間調控的研究進展［J］.亞熱帶農業研究，2016，12（2）：130-135.

［20］ 巫金華，李榮英，里二.胖大海樹開花結果習性觀察［J］.中藥材，1998，21（10）：487-489.

［21］ 張巧巧，VAN HUONG P，陳昌雄，等.遮蔭處理對胖大海苗期生長特點的影響［J］.江西農業大學學報，2016，38（5）：927-934.

［22］ 張巧巧，VAN HUONG P，陳昌雄，等.母樹徑級與貯藏對胖大海種子發芽率的影響［J］.福建農林大學學報（自然科學版），2017，46（2）：159-165.

［23］ VAN HUONG P，陳昌雄，張巧巧，等.灌草對胖大海種群更新幼樹幼苗密度的影響［J］.西南林業大學學報，2016，36（4）：1-8.

［24］ VAN HUONG P，陳昌雄，洪偉，等.表層土壤濕度和pH值對胖大海幼樹幼苗更新的影響［J］.福建農林大學學報（自然科學版），2016，45（2）：173-178.

［24］ VAN HUONG P，陳昌雄，洪偉，等.不同林分下胖大海幼樹幼苗更新分布格局［J］.森林與環境學報，2015，35（4）：377-384.