劉芙蓉，羅建勛，楊馬進(jìn)

(1.四川省林業(yè)科學(xué)研究院林業(yè)研究所，四川 成都 610081； 2.成都理工大學(xué)地球科學(xué)學(xué)院，四川 成都 610059)

山桐子的地理分布及其潛在適宜栽培區(qū)區(qū)劃

劉芙蓉1，2，羅建勛1*，楊馬進(jìn)1

(1.四川省林業(yè)科學(xué)研究院林業(yè)研究所，四川 成都 610081； 2.成都理工大學(xué)地球科學(xué)學(xué)院，四川 成都 610059)

山桐子；Maxent；地理分布區(qū)；適宜栽培區(qū)

山桐子(IdesiapolycarpaMaxim.)為大風(fēng)子科(Flacourtiaceae)山桐子屬(Idesia)落葉喬木，樹干通直，枝條輪生，樹形美觀[1]。山桐子為亞熱帶陽(yáng)性速生樹種，其果實(shí)產(chǎn)量高，呈串狀，入秋后紅艷奪目，被喻為“美麗的樹上油庫(kù)”，是我國(guó)重要的木本油料樹種。果實(shí)和種子含油率高達(dá)20%以上[2]，油脂中不飽和脂肪酸的含量達(dá)70%左右[3]，對(duì)高血壓、高血脂等心腦血管疾病有很好的預(yù)防作用和輔助療效，具有很高的價(jià)值。作為重要的經(jīng)濟(jì)林樹種之一，山桐子具有較大的產(chǎn)業(yè)推廣前景，目前相關(guān)研究主要集中在良種選育、苗木繁殖、栽培技術(shù)和油脂成分等方面[3-5]，而其地理分布區(qū)及氣候特征等方面的研究還較少，對(duì)該樹種的適宜栽培區(qū)進(jìn)行區(qū)劃有利于推動(dòng)山桐子產(chǎn)業(yè)化發(fā)展。

物種分布模型充分利用物種現(xiàn)有的分布資料和環(huán)境數(shù)據(jù)，分析環(huán)境因素與物種分布區(qū)的關(guān)系，并將該關(guān)系特性用于預(yù)測(cè)物種的潛在分布區(qū)。目前，大量基于GIS技術(shù)和數(shù)學(xué)算法發(fā)展而來(lái)的物種分布模型被用于瀕危植物保護(hù)[6-7]、外來(lái)物種入侵風(fēng)險(xiǎn)評(píng)估[8-10]，預(yù)測(cè)氣候變化對(duì)物種分布的影響[11]、野生植物引種栽培區(qū)區(qū)劃[12]等多個(gè)領(lǐng)域，如生態(tài)位因子模型(ENFA)、廣義線性模型(GLM)、最大熵值模型(Maxent)、神經(jīng)網(wǎng)各模型(ANN)、距離模型(Domain)和隨機(jī)森林模型(MARS)均是重要的物種分布模型[13-14]，這些生態(tài)學(xué)模型為物種資源的科學(xué)管理提供了重要依據(jù)。許多研究者對(duì)不同模型的模擬效果進(jìn)行了對(duì)比分析，曹向鋒等[15]利用Maxent、GARP、ENFA、Bioclim和Domain 5種生態(tài)位模型預(yù)測(cè)了黃頂菊在中國(guó)范圍內(nèi)的潛在適生分布區(qū)，并利用受試者工作特征曲線(ROC)比較不同模型的模擬精度，認(rèn)為Maxent的預(yù)測(cè)結(jié)果對(duì)黃頂菊的擬合度較高，預(yù)測(cè)的效果較好。殷曉潔等[16]利用Maxent模擬的蒙古櫟潛在分布區(qū)域覆蓋了蒙古櫟實(shí)際地理分布點(diǎn)的98%，表明模擬的我國(guó)蒙古櫟地理分布與實(shí)際分布非常符合，并利用我國(guó)蒙古櫟的地理分布概率與主導(dǎo)氣候因子的關(guān)系，得出了主導(dǎo)氣候因子的閾值。Maxent模型是基于最大熵理論發(fā)展起來(lái)的生態(tài)位模型，在物種潛在分布區(qū)的預(yù)測(cè)中效果較好，被廣泛用于經(jīng)濟(jì)作物的栽培區(qū)適宜性評(píng)價(jià)。采用物種分布模型對(duì)其地理分布范圍和主要環(huán)境因子進(jìn)行分析，對(duì)于指導(dǎo)山桐子的栽培推廣和產(chǎn)業(yè)布局具有重要意義。因此，本研究以山桐子的野生分布數(shù)據(jù)為基礎(chǔ)，采用Maxent模型對(duì)其適宜栽培區(qū)進(jìn)行預(yù)測(cè)，以期為山桐子推廣栽培的生產(chǎn)實(shí)踐工作提供指導(dǎo)依據(jù)。

1 材料與方法

1.1 山桐子分布數(shù)據(jù)的收集與整理

查詢中國(guó)數(shù)字植物標(biāo)本館(http://www.cvh.org.cn/)和國(guó)家科技部教學(xué)標(biāo)本資源共享平臺(tái)(http://mnh.scu.edu.cn/)等數(shù)據(jù)庫(kù)，收集全國(guó)范圍內(nèi)山桐子標(biāo)本信息，選取野生分布的位置，并查閱相關(guān)文獻(xiàn)記錄。結(jié)合山桐子實(shí)地調(diào)查，采用GPS對(duì)山桐子野生資源分布位置進(jìn)行定位，整理分布樣點(diǎn)，剔除重復(fù)的數(shù)據(jù)分布點(diǎn)，共確定了220個(gè)野生山桐子分布點(diǎn)，利用ArcGIS軟件整理成shp格式文件，用于繪制山桐子地理分布圖。

1.2 環(huán)境數(shù)據(jù)集來(lái)源

本文選擇了影響植物生長(zhǎng)和分布的環(huán)境因子數(shù)據(jù)集，主要包括：氣候、海拔和土壤條件，共24個(gè)因子作為環(huán)境數(shù)據(jù)集。氣候和海拔數(shù)據(jù)來(lái)源于世界氣候數(shù)據(jù)庫(kù)(http://www.worldclim.org)，氣候數(shù)據(jù)是采用各地氣象站的氣象信息進(jìn)行插值生成的全球氣候數(shù)據(jù)，包括與溫度和降水量相關(guān)的19個(gè)生物氣候變量[17]，包括：年平均氣溫(bio1)、平均月氣溫變化范圍(bio2)、等溫性(bio3)、氣溫季節(jié)變化方差(bio4)、極端最高氣溫(bio5)、極端最低氣溫(bio6)、年氣溫幅度(bio7)、雨季平均氣溫(bio8)、干季平均氣溫(bio9)、最暖季平均氣溫(bio10)、最冷季平均氣溫(bio11)、年降水量(bio12)、最濕月降水量(bio13)、最干月降水量(bio14)、季降水量變異系數(shù)(bio15)、最濕季降水量(bio16)、最干季降水量(bio17)、最暖季降水量(bio18)、最冷季降水量(bio19)，數(shù)據(jù)空間分辨率為30″。

1.3 潛在分布區(qū)模擬及檢驗(yàn)

2 結(jié)果與分析

2.1 地理分布

2.2 氣候特征

利用最大熵模型評(píng)價(jià)25個(gè)環(huán)境因子對(duì)野生山桐子的貢獻(xiàn)率，以篩選并確定影響山桐子分布的主導(dǎo)因子。本文選擇累積貢獻(xiàn)率達(dá)85%的因子作為影響山桐子分布的主要因子，其貢獻(xiàn)率見(jiàn)表1，按順序大小依次為：最冷季平均氣溫(bio11)>極端最低氣溫(bio6)>平均月氣溫變化范圍(bio2)>氣溫季節(jié)變化方差(bio4)>雨季平均氣溫(bio8)>年降水水(bio12)>海拔>極端最高氣溫(bio5)>年均氣溫(bio1)。貢獻(xiàn)率的排序表明:低溫對(duì)山桐子分布區(qū)的影響極其重要，限制著植物是否能安全越冬；月平均氣溫變化范圍和氣溫季節(jié)變化方差則反映了物種對(duì)溫差變化的需求是否得到滿足，影響植物的正?；ㄑ糠只伴_(kāi)花結(jié)實(shí)等；雨季平均氣溫反映了山桐子在生長(zhǎng)季對(duì)水熱同步的需求；年降水量反映了植物對(duì)水分條件的需求和耐受能力；海拔反映了山桐子在垂直方向上的分布范圍；極端最高氣溫反映了山桐子對(duì)高溫的耐受范圍；年均氣溫反映了植物對(duì)全年總熱量的需求范圍。前9個(gè)累積貢獻(xiàn)率達(dá)86.6%，這表明山桐子對(duì)溫度的要求很高，而土壤因子的貢獻(xiàn)率相對(duì)較小，這可能是由于山桐子根系發(fā)達(dá)，耐貧瘠，對(duì)土壤養(yǎng)分條件的要求不太嚴(yán)格，因此，山桐子在重慶等地區(qū)的石漠化治理造林中應(yīng)用較廣泛[20]。

表1 影響山桐子分布的環(huán)境因子的貢獻(xiàn)率

2.3 潛在適宜栽培區(qū)

Maxent模型預(yù)測(cè)結(jié)果的AUC值為0.956，表明所構(gòu)建的模型預(yù)測(cè)精度達(dá)到了非常好的水平。山桐子潛在適宜栽培區(qū)的預(yù)測(cè)結(jié)果與實(shí)際分布區(qū)有很好的一致性，其潛在適宜栽植區(qū)主要集中在我國(guó)南方地區(qū)。山桐子的潛在最適宜和較高適宜栽培區(qū)主要分布在四川盆周山區(qū)、大巴山、武陵山、苗嶺、雪峰山、武功山、武夷山、九連山、戴云山、大別山、伏牛山、云貴高原以東、大瑤山和云開(kāi)大山以北地區(qū)；而四川盆中丘陵地區(qū)、兩湖平原和鄱陽(yáng)湖平原地區(qū)為一般適宜區(qū)，經(jīng)人工栽培試點(diǎn)的驗(yàn)證，山桐子在該區(qū)域可進(jìn)行正常的生長(zhǎng)和結(jié)實(shí)；在一般適宜區(qū)與不適宜區(qū)的過(guò)渡地帶，為較低適宜區(qū)；而當(dāng)適宜性P<0.05的地區(qū)，不宜栽植山桐子。

由表2可知：山桐子在我國(guó)的最適宜區(qū)面積為431 717 km2，較適宜生長(zhǎng)區(qū)面積為725 871 km2，一般適宜區(qū)為677 229 km2，較低適宜區(qū)為337 631 km2。適宜栽培區(qū)總面積較大的主要包括：四川、重慶、貴州、陜西、湖北、湖南、江西、安徽、浙江、河南、廣東、廣西、福建和臺(tái)灣地區(qū)；而云南、甘肅、江蘇、山東等省的適宜栽培區(qū)面積相對(duì)較小，統(tǒng)一歸并在其它地區(qū)。

表2 山桐子主要分布行政區(qū)內(nèi)不同等級(jí)的適宜栽培區(qū)面積

從最適宜栽培區(qū)來(lái)看，江西的面積最大，為65 744 km2，其次為湖南、湖北、貴州、四川、廣東、廣西、陜西和福建，而重慶、浙江、安徽、河南和臺(tái)灣的最適宜栽培區(qū)面積相對(duì)較小。從較高適宜栽培區(qū)來(lái)看，湖南的面積最大，為113 756 km2，其次為貴州、江西、福建、廣西、四川、湖北、浙江、廣東和重慶，而河南、陜西、臺(tái)灣的較高適宜栽培區(qū)面積相對(duì)較小。從一般適宜栽培區(qū)來(lái)看，廣西的面積最大，為98 603 km2，其次為四川、河南、湖北、湖南、廣東、浙江、安徽、貴州、陜西、重慶和福建，而江西、臺(tái)灣的一般適宜栽培區(qū)面積相對(duì)較小。從較低適宜栽培區(qū)來(lái)看，河南的面積最大，為34 439 km2，其次為廣西、廣東、四川、安徽、陜西、湖北、貴州地區(qū)，而浙江、臺(tái)灣、湖南、福建、重慶、江西的較低適宜栽培區(qū)面積相對(duì)較小。

3 討論

本文收集野生山桐子分布點(diǎn)220個(gè)，利用Maxent生態(tài)學(xué)模型預(yù)測(cè)了山桐子在中國(guó)的適生范圍，為山桐子適宜栽培區(qū)區(qū)劃提供了依據(jù)。Maxent模型在物種的潛在分布區(qū)研究中具有很高的價(jià)值[22]，本文預(yù)測(cè)的適宜栽培區(qū)與野生山桐子的實(shí)際分布區(qū)具有較高的一致性，精確度高，這與大多數(shù)研究結(jié)論一致[23]，可為該經(jīng)濟(jì)林樹種的產(chǎn)業(yè)化栽培范圍提供指導(dǎo)依據(jù)，有效避免盲目推廣發(fā)展。

4 結(jié)論

通過(guò)Maxent最大熵值模型對(duì)我國(guó)范圍內(nèi)的山桐子適宜栽培區(qū)進(jìn)行預(yù)測(cè)和區(qū)劃，結(jié)果表明：該樹種的栽培區(qū)與野生山桐子分布區(qū)范圍基本一致，主要集中在我國(guó)南方地區(qū)。潛在最適宜栽培區(qū)面積為431 717 km2，較適宜區(qū)面積為725 871 km2，一般適宜區(qū)為677 229 km2，較低適宜區(qū)為337 631 km2。山桐子的潛在最適宜和較高適宜栽培區(qū)主要分布在四川盆周、大巴山、武陵山、苗嶺、雪峰山、武功山、武夷山、九連山、戴云山、大別山、伏牛山、云貴高原以東、大瑤山和云開(kāi)大山以北地區(qū)；而四川盆中丘陵地區(qū)、兩湖平原和鄱陽(yáng)湖平原地區(qū)為一般適宜區(qū)。

綜上所述，本文對(duì)我國(guó)范圍內(nèi)野生山桐子的分布點(diǎn)、氣候和土壤因子進(jìn)行了分析，基本考慮到了山桐子分布的實(shí)際生態(tài)位，給出了山桐子的天然分布范圍，以及適生區(qū)的年均氣溫、降水量等氣候因子，區(qū)劃了山桐子的栽培適宜區(qū)，包括山桐子的潛在最適宜栽培區(qū)、較高適宜栽培區(qū)、一般適宜栽培區(qū)和較低適宜栽培區(qū)。本文在生物氣候因子的基礎(chǔ)上，加入了土壤因子，優(yōu)化分析了山桐子栽培區(qū)范圍和適宜性等級(jí)，為我國(guó)山桐子產(chǎn)業(yè)的推廣栽培提供參考。在今后的研究中，應(yīng)加強(qiáng)區(qū)域栽培試驗(yàn)工作，更好地對(duì)適宜栽培區(qū)預(yù)測(cè)工作進(jìn)行驗(yàn)證，以提高經(jīng)濟(jì)林樹種在產(chǎn)業(yè)化發(fā)展進(jìn)程中的推廣栽培效果。

[1] 江錫兵，龔榜初，李大偉，等. 山桐子自然群體表型性狀變異分析[J].林業(yè)科學(xué)研究，2013,26(1):113-117.

[2] 龔榜初，李大偉，江錫兵，等.不同產(chǎn)地山桐子果實(shí)含油率及其理化指標(biāo)變異分析[J].西北植物學(xué)報(bào)，2012a,32(8):1680-1685.

[3] 王艷梅，王東洪，劉 震，等. 不同種源山桐子種子含油率和脂肪酸GC/MS分析[J].經(jīng)濟(jì)林研究，2011,29(2):14-21.

[4] 歲立云，劉曉敏，李周岐，等. 山桐子果實(shí)性狀的自然變異及類型劃分[J].西北農(nóng)林科技大學(xué)學(xué)報(bào):自然科學(xué)版，2009,37(8):115-120.

[5] 龔榜初，李大偉，江錫兵，等.不同種源山桐子果實(shí)脂肪酸組成變異分析[J].植物生理學(xué)報(bào)，2012b,48(5):505-510.

[6] Adhikari D, Barik S K, Upadhaya K. Habitat distribution modelling for reintroduction ofIlexkhasianaPurk., a critically endangered tree species of northeastern India[J]. Ecological Engineering, 2012, 40(3):37-43.

[7] 姜建福，樊秀彩，張 穎，等.中國(guó)三種瀕危葡萄屬(VitisL.)植物的地理分布模擬[J].生態(tài)學(xué)雜志，2014,33(6):1615-1622.

[8] Giljohann K M, Hauser C E, Williams N S G,etal. Optimizing invasive species control across space: willow invasion management in the Australian Alps[J]. Journal of Applied Ecology, 2011, 48(5):1286-1294.

[9] Richter R, Dullinger S, Essl F,etal. How to account for habitat suitability in weed management programmes[J]. Biological Invasions, 2013, 15(3):657-669.

[10] Chapman D S, Makra L, Albertini R,etal. Modelling the introduction and spread of non-native species: international trade and climate change drive ragweed invasion[J]. Global Change Biology, 2016, 22(9):3067-3079.

[11] Remya K, Ramachandran A, Jayakumar S. Predicting the current and future suitable habitat distribution ofMyristicadactyloidesGaertn. using MaxEnt model in the Eastern Ghats, India[J]. Ecological Engineering, 2015, 82(9):184-188.

[12] 胡 秀，郭 微，吳福川，等．MaxEnt生態(tài)學(xué)模型在野生植物近自然林引種區(qū)劃中的應(yīng)用—以紅姜花為例[J].廣西植物，2015， 35(3):325-330.

[13] 李國(guó)慶，劉長(zhǎng)成，劉玉國(guó)，等.物種分布模型理論研究進(jìn)展[J].生態(tài)學(xué)報(bào)，2013,33(16):4827-4835.

[14] Yuan H S, Wei Y L, Wang X G. Maxent modeling for predicting the potential distribution of Sanghuang, an important group of medicinal fungi in China[J]. Fungal Ecology, 2015, 17:140-145.

[15] 曹向鋒，錢國(guó)良，胡白石，等.采用生態(tài)位模型預(yù)測(cè)黃頂菊在中國(guó)的潛在適生區(qū)[J].應(yīng)用生態(tài)學(xué)報(bào)，2010,21(12):3063-3069.

[16] 殷曉潔，周廣勝，隋興華，等.蒙古櫟地理分布的主導(dǎo)氣候因子及其閾值[J].生態(tài)學(xué)報(bào)，2013,33(1):103-109.

[17] 張興旺，李 垚，方炎明.麻櫟在中國(guó)的地理分布及潛在分布區(qū)預(yù)測(cè)[J].西北植物學(xué)報(bào)，2014,34(8):1685-1692.

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

[19] 劉少軍，周廣勝，房世波.中國(guó)橡膠種植北界[J].生態(tài)學(xué)報(bào)，2016,36(5):1272-1280.

[20] 馮大蘭，黃小輝，劉 蕓,等.4種木本植物在石漠化地區(qū)的生長(zhǎng)狀況及光合特性[J].北京林業(yè)大學(xué)學(xué)報(bào)，2015，37(5):62-69.

[21] 李 垚，張興旺，方炎明.氣候變暖對(duì)中國(guó)栓皮櫟地理分布格局影響的預(yù)測(cè)[J].應(yīng)用生態(tài)學(xué)報(bào)，2014,25(12):3381-3389.

[22] Wang Y L. Predicting the potential geographic distribution of crofton weed (Ageratinaadenophora) around the world using Maxent modeling[J]. International Journal of Plant Research, 2012, 25(2):324-335.

[23] 王 偉，田榮榮，那立妍，等.基于MaxEnt生態(tài)軟件劃分澳洲堅(jiān)果的潛在地理適生區(qū)[J].林業(yè)科學(xué)研究，2017,30(3):444-449.

[24] Wong M H G, Li R, Xu M,etal. An integrative approach to assessing the potential impacts of climate change on the Yunnan snub-nosed monkey[J]. Biological Conservation, 2013, 158(2):401-409.

[25] 王東洪.不同種源山桐子冬芽休眠的溫度特性研究[D].鄭州河南農(nóng)業(yè)大學(xué)， 2012.

[26] 段居琦，周廣勝.中國(guó)雙季稻種植區(qū)的氣候適宜性研究[J].中國(guó)農(nóng)業(yè)科學(xué)，2012，45(2):218-227.

[27] 胡 秀，吳福川，郭 微.基于MaxEnt生態(tài)學(xué)模型的毛姜花潛在園林引種區(qū)預(yù)測(cè)[J].中國(guó)城市林業(yè)，2013,11(4):28-31.

GeographicDistributionandRegionalizationofPotentialSuitableCultivatedAreaforIdesiapolycarpa

LIUFu-rong1,2,LUOJian-xun1,YANGMa-jin1

(1. Research Institute of Forestry, Sichuan Academy of Forestry, Chengdu 610081, Sichuan, China ; 2. College of Earth Sciences, Chengdu University of Technology, Chengdu 610059, Sichuan, China)

Idesiapolycarpa; Maxent; geographic distribution; suitable cultivated area

10.13275/j.cnki.lykxyj.2017.06.021

2017-05-05

山桐子高產(chǎn)、高含油良種選育與栽培示范(2015BAD15B02)；四川省科技計(jì)劃項(xiàng)目“木本油料突破性新品種選育、育種材料與方法創(chuàng)新”(2016NYZ0035-06)；四川省植物資源共享平臺(tái)建設(shè)項(xiàng)目(TJPT20160021)

劉芙蓉(1987—)，女，四川成都人，在讀博士.主要研究方向：林業(yè)資源遙感.E-mail:liu729103540@163.com

* 通訊作者：羅建勛(1964—)，男，四川成都人，博士，研究員.主要研究方向：林木種質(zhì)資源和林木遺傳育種.E-mail:jianxunl@163.com

S794

A

1001-1498(2017)06-1028-06

詹春梅)