近30年海南島橡膠林時空變化分析

李廣洋，寇衛利，陳幫乾，吳志祥，張希財，云 挺，馬 俊，孫 瑞，李 瑩

(1．中國熱帶農業科學院橡膠研究所/海南儋州熱帶農業生態系統國家野外科學觀測研究站，海南 ?？?571101；2．西南林業大學大數據與智能工程學院，云南 昆明 650224；3．南京林業大學林學院，江蘇 南京 210037；4.復旦大學生物多樣性與生態工程教育部重點實驗室，復旦大學生命科學學院，上海 200438；5.中國航天科工信息技術研究院，北京 100000)

天然橡膠主要產自于巴西橡膠樹(Heveabrasiliensis)，是重要的工業原料和戰略物資[1]。為了滿足國防和經濟發展的需要，20世紀50年代開始,我國在海南島、云南省和廣東省大面積種植橡膠樹。歷經幾代人的努力，已在熱帶北緣高緯度的“次優”植膠環境取得巨大成功，極大地推動了世界植膠業的發展[2]。2002年后天然橡膠價格迅速攀升，全面促進了全球天然橡膠種植面積的飛速增長[3]，但隨后因天然橡膠庫存增加、國際經濟環境因素等影響，膠價自2011年后迅速回落并持續低迷至今。前期橡膠林面積的迅速增大加劇了熱帶地區的景觀破碎化和片段化[4]，對生態環境以及生物多樣性造成了威脅[5-7]，而近10年膠價的持續低迷直接影響膠農的種植積極性，將導致橡膠林土地利用類型的改變，從而對中國天然橡膠的供給和當地生態環境再次帶來巨大的影響。因此，準確及時地監測橡膠林的時空變化，對促進天然橡膠產業可持續發展有重要意義。

大面積橡膠林的精確識別是時空動態監測的基礎[8]。在國內，陳匯林等[9]基于250 m分辨率MODIS(中等分辨率成像光譜儀)表明落葉和抽葉期最利于橡膠林識別。但由于MODIS空間分辨率較低，識別斑塊較小的橡膠林還存在較大的局限性[10-11]，中分辨率影像逐漸成為分類的主流數據源。如Dong等[12]基于多時相Landsat影像分析了橡膠林的物候變化規律并提取了海南儋州地區的橡膠林；Xiao等[13]利用Landsat識別了西雙版納的橡膠林。由于在常年多云雨的熱帶地區難以獲得橡膠關鍵物候期的大面積無云影像[14]，該類算法在較大尺度推廣仍然十分困難，而聯合多源長時序遙感大數據監測的優勢則越來越明顯。Chen等[15]聯合PALSAR雷達和233景Landsat 光學影像實現了海南橡膠林信息的準確提取?？傮w而言，小尺度的橡膠林遙感識別已取得長足的發展[16-17]，但無云光學影像匱乏、橡膠物候差異及其他落葉林干擾[15]、景觀破碎化[18]等仍是當前將算法向國家或更大尺度推廣的主要阻礙。

在遙感識別的基礎上，國內外學者已在橡膠林時空動態監測方面開展了一些研究。比如廖諶婳等[19]對西雙版納1990—2014年的橡膠林進行了研究，表明西雙版納橡膠林面積總體呈增長趨勢，主要增長地為景洪市和勐臘縣，林分呈成熟林和老齡林態勢，海拔800 m以下大面積天然林被轉化為橡膠林?？苄l利[20]對西雙版納1990—2010年間的橡膠林時空擴張進行了監測，表明橡膠林面積受橡膠價格主導，面積總體呈增長趨勢，并有不斷向高海拔、陡坡度以及山北坡等次適宜種植地區發展的趨勢。Xiao等[21]利用多重歸一化方法對1987—2018年的西雙版納橡膠林進行了提取，表明近幾十年西雙版納的橡膠林呈現從集中到分散的擴張趨勢，并且約九成橡膠林都分布在景洪市和勐臘縣周圍。Liu等[22]探究了1980—2010年間中老緬交界區的地形因素對橡膠林種植的影響，結果表明該區域橡膠林面積擴大了近9倍，主要分布在南部和東部地區，集中分布在海拔600～1 000 m及8°～25°的斜坡上。

設立國家級自然保護區對陸地生態系統的保護起到了積極作用。Yang等[23]對西雙版納納板河流域國家級自然保護區內低地森林、低地橡膠等6個主要地類的碳儲量進行了評估，結果表明橡膠園并不能補償砍伐森林所導致的碳損失，但自然保護區的建立能有助于減輕橡膠林或其他土地利用的擴張所造成的碳損失。冷靜等[24]就納板河流域國家級自然保護區內橡膠林對畜禽養殖數量的影響進行了評價，結果表明橡膠林的擴張致使區域環境退化、喪失生物多樣性，最終導致許多依賴森林生存的畜禽失去了活動場所及食物，在一定程度上破壞了生物鏈?？偠灾?，近幾十年橡膠林的快速擴張及潛在的生態問題得到了前所未有的關注，但主要集中在西雙版納及周邊的熱點區域。海南島擁有大面積的熱帶雨林，一直是印度—緬甸生物多樣性熱點區[25]，同時也是中國第二大植膠區，橡膠林總面積占海南島森林總面積的四分之一以上，是當地最大的人工林生態系統。因此，分析海南島橡膠林的時空變化特征，對實際生產和生態評估都有重要意義。本研究聯合遙感大數據和橡膠林物候等信息，以5 a為步長提取海南島1990—2020年共7期橡膠林分布圖，分析研究區近30年的橡膠林時空變化情況，以期為海南島橡膠林的科學栽培、土地合理規劃提供技術支持，并對相關國家級自然保護區的建設與監管提供一定的數據支撐。

1 材料與方法

1.1 研究區概況

海南島地處我國華南地區(108°21′～111°03′E，18°20′～20°10′N)，北以瓊州海峽為界，西臨北部灣，東瀕南海，陸地總面積3.39萬km2。海南島由山地、丘陵、平原和臺地等地貌構成，以五指山和鸚哥嶺為中心，地勢向四周逐級降低，呈明顯的梯狀結構。島內動植物資源豐富，是我國最大的熱帶寶地。海南島雨量豐沛，日照充足，受冬季寒潮影響較小，很適合橡膠樹生長，但受臺風影響較大。島內橡膠林面積僅次于云南省，是中國第二大植膠區。

研究區內有4個國家級自然保護區，分別為霸王嶺保護區、尖峰嶺保護區、五指山保護區和吊羅山保護區，分布在海南島中西部和東南部的山地區域。保護區內森林組成種類復雜，群落結構發育良好，層間植物特征豐富，是典型的熱帶雨林特征。每個保護區按功能又劃分為核心區、緩沖區和實驗區。核心區是指保存完好的天然狀態的生態系統，以及珍稀、瀕危動植物的集中分布地，且核心區內禁止任何單位和個人進入；緩沖區指核心區外圍規劃出的一定面積的區域，只允許開展科學研究觀測活動；實驗區是指緩沖區外圍的多用途區域，一般可以開展科學實驗、旅游以及馴化、繁殖珍稀與瀕危野生動植物等一系列活動。

1.2 影像數據及處理

Landsat系列表觀反射率影像來自美國地質勘探局，Sentinel-2 MSI 影像來自歐空局，已全部存于GEE(Google Earth Engine)平臺。GEE已被廣泛應用于農業、林業和自然災害等領域的監測研究[26]。由于海南島常年多云雨天氣，需要使用多年的影像才能滿足研究需求。2018年后Sentinel-2A/B星組網后數據更豐富，因此2020期橡膠林制圖使用2019年和2020年的影像數據；對于2015年及以前，使用目標年份及前后1年所有Landsat和Sentinel-2影像進行分類。覆蓋海南島需要4個WRS2條帶(124/46、124/47、123/46和123/47)，所使用的Landsat系列影像共1 710景，Sentinel-2為3 391景。

Landsat系列影像采用自帶掩碼波段對云和陰影以及ETM+的缺失條帶等無效觀測進行掩膜[14]，Sentinel-2 MSI使用SentinelHub提供的云概率圖層進行掩膜(閾值為50%)。Landsat與Sentinel-2傳感器歸一化的參數來自文獻[27]。云掩膜結束后，計算歸一化植被指數(NDVI，normalized difference vegetation index，式中記為INDVI)和地表水分指數(land surface water index，LSWI，式中記為ILSWI)。

(1)

(2)

式中：ρRed、ρNIR、ρSWIR分別為紅光、近紅外和短波紅外波段的反射率。植被指數計算完畢后，分別對目標年份1—3月[橡膠快速落葉和抽葉期，年積日(DOY為1～90 d]獲取的影像集的LSWI進行最小值和標準差合成，分別命名為LSWImin和LSWIstd，對于4—12月(DOY為90～365 d)獲取的影像集的NDVI進行最大值合成，記為NDVImax。

1.3 地面數據源及處理

2010、2015和2020年的地面數據來源于實地調研和谷歌高清影像目視解譯。野外調研始于2011年11月，幾乎每年都會在海南島開展橡膠林等不同典型植被調查，動態標記不同植被樣本和拍攝帶GPS信息的景觀照片，先后獲得2010、2015和2020年800多個橡膠林及500余個非橡膠林樣本[15,28]?；?015年的橡膠林分布圖和1987年以來的所有Landsat系列影像，獲得了年尺度的橡膠林建立時間分布圖(等效于林齡圖)[28]。橡膠林的經濟壽命周期為30 a左右，因此2005年及以前的地面樣本從該林齡圖中抽樣生成，主要步驟為：①在林齡圖中生成種植時間比目標年份早8 a(幼苗期)的橡膠林掩碼圖層；②在掩碼圖層上利用GEE抽樣函數隨機抽樣200個樣點，以點為中心生成半徑為30 m的緩沖區；③結合歷年Landsat影像植被指數時序和快視圖對抽樣結果進行二次確認，刪除部分不合格的抽樣點。檢視結束后隨機抽取地面樣本的60%用于確定算法閾值，40%用于分類精度評估。

1.4 其他數據源及處理

數字高程模型(digital elevation model,DEM)產品為SRTM(Shuttle Radar Topography Mission)的V3版本，來自美國宇航局，空間分辨率為1弧秒(約30 m)[29]。海拔信息直接從DEM中獲取，坡度通過GEE提供的地形函數計算獲得。

我認為，在目前的發展階段，部分產業談共享經濟也許是過早的。但我并不認為有些產業注定沒有辦法走網絡協同的道路，很可能未來它們會在一個更大的網絡協同或者說是協同網絡中去發展。像租車，本身并不適合共享經濟，但是智能交通可能會出現一個巨大的網絡協同平臺，而租車只是其中的一個服務。所以，共享平臺或者說共享經濟、網絡協同也許能在一個更大的產業范圍內實現。這就需要更長時間的積累。

1.5 橡膠林提取方法

橡膠樹是高大的熱帶喬木，每年2—3月會經歷一個快速集中的落葉和抽葉過程，此時與熱區廣泛分布的常綠森林有明顯的可分離性[12]。另外，由于人工定期施肥管理，在生長季其冠層也比許多其他林分要茂盛。前期研究表明，落葉期橡膠林的LSWI值很低(LSWImin)，快速落葉與抽葉期因LSWI變化劇烈從而產生較大的標準差(LSWIstd)，而生長季茂盛的冠層催生了較大的NDVI值(NDVImax)，據此建立了基于物候的橡膠林分類算法[15]。2010和2015年的橡膠林分布圖來源于前期研究[15,28]，流程是先從PALSAR/PALSAR-2雷達數據獲得高精度的森林本底圖，再使用LSWImin、LSWIstd和NDVImax從森林本底圖中分離橡膠林。對于缺乏雷達數據的年份，結合當年的地面樣本和合成的LSWImin、LSWIstd和NDVImax圖層，繪制直方圖并提取2.5%和97.5%的分位數為分類閾值(適當圓整)，將對應指標位于分位數區間的像元的交集定義為橡膠林。由于橡膠樹抽葉期的光學影像非常利于目視解譯，結合單景無云影像可進一步提升分類精度。因此，在目標年份前后1年內查找橡膠樹抽葉期的1～2景關鍵影像，再次結合地面樣本繪制LSWI和NDVI等直方圖，提取閾值對未受云影像區域的合成影像分類結果進行優化。分類完成后，計算斑塊連接度，去掉連接度小于4的小斑塊。用混淆矩陣進行精度評估表明，1990—2020年共7期橡膠林分類圖(最終為橡膠林和非橡膠林兩大類)的總體精度分別為90.21%、91.57%、92.31%、91.87%、95.24%、96.77%和94.67%，能夠滿足對林分時空變化分析的需要。

1.6 橡膠林時空分析方法

結合省級和市/縣級行政邊界矢量圖、DEM等基礎數據，在GEE云平臺中對1990—2020年海南島7期橡膠林分布圖進行面積統計、地形分布及變化分析，利用ArcGIS 10.4.1和Origin 2017進行制圖。主要步驟如下：

1)上傳海南島省級和市/縣尺度的行政邊界矢量圖至GEE平臺，分別裁剪1990、1995、2000、2005、2010、2015和2020年的海南島橡膠林分布圖，計算全省和各市縣區的橡膠林面積，分析各市縣區內橡膠林面積增減情況。

2)疊加DEM數據進行地形分布及變化分析。根據GB/T 15772—1995《水土保持綜合治理規劃通則》的坡度分級體系，將研究區坡度分為6個等級，分別為微坡(0～5°)、較緩坡(>5°～8°)、緩坡(>8°～15°)、較陡坡(>15°～25°)、陡坡(>25°～35°)和急陡坡(>35°)。海拔按0、50、100、200、300、400和600 m為界分為7個等級。分別統計各個年份不同坡度和海拔的橡膠林面積，并分析其1990—2020年的時空變化特征。

3)上傳海南島國家級自然保護區矢量邊界數據至GEE平臺，加載并裁剪保護區內的橡膠林，計算保護區內的實驗區、緩沖區和核心區3個區域內橡膠林面積并分析變化情況。

2 結果與分析

2.1 海南島橡膠林空間分布及變化

1990—2020年共7期橡膠林分布圖如圖1所示。從圖中可以看出，在2000年以前，海南島的橡膠林分布相對穩定，空間變化不明顯(圖1A—1C)。但是，2000年后橡膠林面積迅速增長，尤其在海南島中西部及中北部地區，空間變化明顯，種植面積顯著增加(圖1D—1G)?？傮w上看，近30年間海南島橡膠林由中部山區向西北和北部地區擴張顯著，南部地區有所減少，總體呈現出“北增南減”的變化趨勢(圖1H)。

圖1 近30年海南島橡膠林分布圖Fig.1 The spatial distribution maps of rubber plantations on Hainan Island over the past 30 years

1990—2000年面積變化不明顯，2000—2015年迅速增加，到2020年略有下降(圖1I)。1990年海南島橡膠林面積約為34.44 萬hm2，2000年以后迅速增加，至2015年達到59.77萬hm2，2020年降至58.58 萬hm2。近30年橡膠林面積凈增長70.11%。1995—2000年間橡膠林面積減少了1.68 萬hm2，約減少4.34%；2010—2015年間增長最多，約增長16.63%，達到了8.52萬hm2，2015—2020年減少了約1.99%，為1.19萬hm2?？傮w看，海南橡膠林的面積保持明顯的增長趨勢(圖1I)。

2.2 海南島各市縣橡膠林分布及變化

1990—2020年海南島各市縣的橡膠林面積如圖2所示。近30年來儋州市橡膠林平均面積穩居第一，瓊中縣、白沙縣及澄邁縣緊居其次，這4個市縣在近30年一直是海南島的橡膠主產區域，其次是樂東縣以及屯昌縣。橡膠林面積較小的市縣主要是五指山市、?？谑?、文昌市等，其中陵水縣面積最小。從面積變化看，儋州市、白沙縣和澄邁縣等橡膠林面積增長趨勢非常顯著，而東方市、保亭縣以及三亞市的橡膠林面積呈負增長趨勢，2020年橡膠林面積較1990年均有所下降。

從具體數值看(圖2)，儋州市、白沙縣、澄邁縣以及瓊中縣4個市縣在2020年橡膠林面積共計約30.22 萬hm2，占全島橡膠林面積的51.58%。其中，儋州市2020年橡膠林面積約為12.65萬hm2，占全島橡膠林總面積的21.60%，較1990年增長了約7.83 萬hm2。2020年三亞市、東方市、文昌市以及陵水縣4個橡膠林面積最小，共計約3.06 萬hm2，僅占全島橡膠林面積的5.22%。海南島北部的儋州市、白沙縣、澄邁縣等9個市縣橡膠林面積近30年共計增長了27.09 萬hm2，占全島橡膠林面積的46.24%；而南部的東方市、保亭縣、三亞市、樂東縣以及陵水縣共計減少了3.54 萬hm2，為海南島橡膠林主要減少區域。

圖2 1990—2020年海南島各市/縣橡膠林面積動態Fig.2 The area dynamics and statistics of rubber plantations in different cities/counties of Hainan Island from 1990 to 2020

2.3 海南島橡膠林的地形分布及變化

2020年海南島海拔0～50 m區域橡膠林面積為5.27萬hm2(圖3A)，較1990年增長了約2.77倍;海拔>50～100 m區域2020年橡膠林面積為14.95 萬hm2，增長2.91倍；>100～200 m橡膠林面積增長1.78倍，2020年達到24.18 萬hm2；>200～300 m區域增長了1.23倍，達到了9.05 萬hm2；>300～400 m區域2020年為3.50 萬hm2，增長0.95倍；海拔600 m以上區域橡膠林面積極少，近30年雖略有增長，但也僅占橡膠林總面積的0.08%。海南島橡膠林主要種植在海拔0～300 m的區域，截至2020年，300 m以下區域橡膠林面積約53.45萬hm2，占海南島橡膠林總面積的91.24%。

圖3 1990—2020年海南島不同海拔(A)和坡度(B)區域橡膠林面積的變化Fig.3 The area changes of rubber plantations in different altitudes(A) and slopes(B) in Hainan Island from 1990 to 2020

海南島0°～5°的微坡區域橡膠林面積在2020年達到25.89 萬hm2，較1990年增長約2.25倍(圖3B)；>5°～8°的較緩坡區域在2020年達到11.84 萬hm2，增長約1.71倍；>8°～15°的緩坡區域近30年增長1.38倍，達到了12.74 萬hm2；>15°～25°的較陡坡區域增長了1.18倍，2020年達到了6.79 萬hm2；>25°～35°的陡坡區域增長1.18倍，達到了1.22 萬hm2，35°以上的急陡坡區域橡膠林雖然也有小幅增長，但面積極小，近30年增長了0.08 萬hm2。海南島橡膠林面積在0°～5°微坡區上升幅度最大，主要種植坡度在0°～25°。截至2020年，25°以下坡地的橡膠林面積已經達到57.26 萬hm2，占橡膠林總面積的97.75%，而35°為橡膠林種植的上限坡度。

2.4 海南國家級自然保護區內橡膠林分布及變化

海南四大自然保護區內均有橡膠林分布，在2020年，霸王嶺保護區內橡膠林面積最大，其次是吊羅山保護區、尖峰嶺保護區和五指山保護區。霸王嶺保護區的橡膠林主要分布在實驗區內，吊羅山和五指山保護區則主要分布在緩沖區和實驗區內，尖峰嶺保護區其主要集中在緩沖區。近30年來，四大自然保護區的橡膠林面積均呈持續減少趨勢。

從具體數值(表1)可以看出，霸王嶺保護區近30年橡膠林面積約減少756.01 hm2，吊羅山保護區減少了約167.03 hm2，尖峰嶺保護區減少約116.89 hm2，橡膠林面積最小的五指山保護區約減少248.14 hm2，四大國家級自然保護區橡膠林總面積近30年間約減少1 288.08 hm2。霸王嶺保護區在2020年的橡膠林面積占四大保護區內橡膠林總面積的85.77%，種植面積居首；吊羅山保護區和尖峰嶺保護區共計占四大保護區內總面積的13.77%，而種植面積最小的五指山保護區橡膠林面積僅占0.46%。霸王嶺保護區在近30年間受橡膠林的擴張影響最為嚴重，五指山保護區受影響最小。

表1 1990—2020年海南島不同國家級自然保護區內橡膠林面積統計Table 1 The statistics of rubber plantation area in different national protection zones in Hainan Island from 1990 to 2020

3 討 論

3.1 海南島橡膠林面積變化及驅動因素

1990—2000年間海南島橡膠林面積呈緩慢增長，2000年后迅速增加，2015年達到峰值，2020年略有下降?？傮w而言，近30年橡膠種植面積增長非常顯著，2000—2015年間增長最顯著；近30年凈增24.14萬hm2，增長了70.11%，空間分布上呈現明顯的“北增南減”特征。2000年前橡膠種植面積增長緩慢的原因主要有兩方面：一是20世紀90年代早期，國內天然橡膠市場還未完全開放，民營橡膠銷售渠道較窄，抑制農戶種植意愿，同時農墾大部分土地已開發種植；二是20世紀90年代中期基本完成市場化改革后，受亞洲金融危機沖擊，天然橡膠價格下跌，農戶對橡膠種植的觀望情緒較濃。2003年后天然橡膠價格開始迅速攀升，2002年均價約8 900 元/t，2011年均價漲至3.4 萬元/t(最高價約4.3萬元/t)，極大地刺激了天然橡膠林種植面積的增長，尤其是民營橡膠部分，在2009年國內總面積首次超過農墾[30]。2011年后天然橡膠價格迅速回落并持續低迷至今，最低曾跌破1萬元/t，植膠生產積極性遭到破壞，部分膠園被棄割、棄種，這也是2015年后橡膠林面積略有下降的原因[31]。

近30年海南島西北和北部地區橡膠面積增長顯著，特別是儋州市、澄邁縣、臨高縣、瓊海市和定安縣，2000年后增長最明顯，以儋州市為中心的周邊市縣構成海南天然橡膠核心種植區；南部地區的橡膠面積縮減明顯，尤其是東方市和南部的三亞市以及保亭縣，總面積逐年降低。近30年海南島橡膠林總體呈現“北增南減”的變化趨勢，與之前的研究結果相似[32-33]。導致這一變化的主要原因有兩個：其一，東南沿海是臺風破壞重災區，在1952—2005年間登陸海南島大于12級的臺風(風速≥32.7 m/s)有35次，其中33次均從東南沿海登陸，超強臺風如2005年的“達維”和2011年“納沙”給天然橡膠產業造成了巨大損失[34-35];因此，即使在價格行情較好的21世紀前10年，膠農也更愿意將橡膠種植在西北和北部的儋州市、白沙縣和澄邁縣等風災相對較小的市縣。其二，東南部地區地勢相對平緩，城市化的過程也會加快土地利用方式的轉變[36]，南部三亞等市縣還具有天然的熱帶氣候優勢，種植熱帶水果可以提前上市，加快了橡膠林向其他果園轉換的進程[37]。此外，海南島的雨量雖然非常充沛，但空間分布差異很大，整體呈東濕西干、中部山區多雨的格局，比如東方市是典型的雨影區，年均降雨量多在800～1 000 mm，而橡膠樹是喜高溫多雨的多年生喬木，年均降雨量要求1 500 mm以上，這也是導致海南西部地區橡膠林面積逐漸下降的原因之一。

3.2 海南島橡膠林地形變化及驅動因素

在地形分布上，近30年橡膠林的增長集中在海拔低于200 m和坡度<15°的地區，橡膠林的增長未對位于中高海拔地區的熱帶雨林造成影響。海南島約91.24%的橡膠都種植在海拔300 m以下地區，97.75%種植在25°坡度以下，主要以平原和丘陵地帶為種植區。即使在保護區外，近30年來也沒有向高海拔和陡峭山坡大面積擴張，總的來說增長模式良好，這與云南西雙版納和老撾北部地區的橡膠增長模式是不同的?，F有的研究表明，西雙版納橡膠林主要分布在海拔600～1 200 m和坡度在5°～25°的范圍內[20]，老撾北部5省的橡膠林也呈相似分布特征[38]。造成海南島橡膠與西雙版納及老撾北部兩地種植地形差異的主要原因有以下幾點：①海南省出臺一系列的森林保護政策，如2005年的“森林生態補償機制”等，加強了熱帶森林的保護，避免大面積森林轉換為橡膠林[39]；②海南島是獨特的島嶼性氣候，沒有橫斷山脈阻擋寒潮，高海拔地區的橡膠極易受寒害和風害等自然災害影響，而西雙版納和老撾北部冬季有高原和橫斷山脈阻擋寒潮，夏季的水系引導又與橫斷山脈走向一致，季風可帶來暖濕氣團，這一獨特的逆溫現象利于橡膠生存[40-41]；③2000年后，橡膠種植中心向西北和北部地區遷移的過程中，因該地區多為低海拔的平原、丘陵或臺地，有大量耕地轉換為橡膠林[28]，從而降低了平均種植海拔和坡度。下一步可分析氣候等因素對海南島橡膠林種植的影響，為橡膠林種植規劃和適宜性評價提供科學依據及決策支持。

3.3 自然保護區內橡膠林面積變化及生態影響

熱帶雨林能夠促進水循環，對氣候調節和生物多樣性的保護起著至關重要的作用[42]。為保護海南島生物多樣性[43]并促進社會經濟的可持續發展，我國從20世紀60年代起先后在海南島成立了尖峰嶺(1960年)、霸王嶺(1980年)、五指山(1985年)和吊羅山(1994年)自然保護區，并分別于2002、1988、2003年和2008年上升為國家級自然保護區。由于歷史遺留或保護力度不夠的原因，在四大保護區內均有少量的橡膠林分布。但是，建立自然保護區和等級提升的保護效果是非常明顯的，四大保護區內橡膠林面積下降趨勢明顯。近30年間共計減少了1 288.08 hm2，較1990年減少了約68.02%，且四大保護區的核心區內幾乎沒有橡膠種植，說明保護區的建立對該地區的生態恢復起到了至關重要的作用。

相比之下，四大國家級自然保護區中霸王嶺保護區橡膠林面積最大，主要因為它是位于中西部輕風害的山區，水熱條件優越，非常適合橡膠樹種植；其次，它接近海南島的橡膠核心種植區域，周邊大面積種植橡膠形成輻射效應，刺激了保護區中實驗區和緩沖區橡膠的非法種植。五指山自然保護區位于海南島的東南部，周邊橡膠種植規模較小，且平均海拔和坡度最高，因此保護區橡膠分布最小。從國家級自然保護區內的橡膠林面積變化情況來看，種植橡膠并沒有對位于中高海拔地區的熱帶雨林造成破壞，也沒有進一步對海南島的森林和生物多樣性造成明顯威脅，這與其他研究結果是一致的[44-45]。保護區內橡膠林面積的不斷減少，逐步被恢復植被所取代，說明政府建立保護區取得了積極的成效。盡管如此，各保護區內仍有少量橡膠林，建議當地政府出臺一系列森林生態補償機制，并加強對保護區的建設與監管，盡早用其他樹種逐步替代保護區內的橡膠林，在保護生態環境的同時實現可持續發展。