海南島森林植被NEP 季節(jié)性時(shí)空變化規(guī)律及氣候驅(qū)動(dòng)因素分析

劉應(yīng)帥，余 瑞，鄭彬彬，劉嘉慧，宋 奇，陳榮昊，嚴(yán) 哲

（1. 海南大學(xué) 生態(tài)與環(huán)境學(xué)院，海口 570228; 2. 海南集思勘測規(guī)劃設(shè)計(jì)有限公司，海口 570203）

定量描述陸地生態(tài)系統(tǒng)碳動(dòng)態(tài)變化及影響因素對(duì)區(qū)域生態(tài)系統(tǒng)碳循環(huán)的研究很重要[1?2]。而森林凈生態(tài)系統(tǒng)生產(chǎn)力（Net Ecosystem Productivity，NEP）的大小能夠表征生態(tài)系統(tǒng)的固碳能力，它最初由WOODWELL 等[3]在分析陸地生物圈源、匯問題時(shí)提出，其生態(tài)學(xué)的含義表示為生態(tài)系統(tǒng)總的生產(chǎn)力與生態(tài)系統(tǒng)呼吸之差。通過NEP 來指示生態(tài)系統(tǒng)的固碳狀態(tài)，NEP 為正值，說明生態(tài)系統(tǒng)呈現(xiàn)碳匯狀態(tài)，反之則為碳源。陸地生態(tài)系統(tǒng)一直扮演著強(qiáng)大的碳匯角色[4]，而陸地碳匯其中的很大一部分位于熱帶森林區(qū)域。熱帶森林覆蓋面積約占全球表面積的10%，其森林碳匯約占陸地總碳庫的25%，并且占全球植被碳儲(chǔ)存的50%[5?6]。由于熱區(qū)異常氣候頻發(fā)，熱帶森林生態(tài)系統(tǒng)碳匯的估算對(duì)于氣候的敏感性的分析研究存在較大的不確定性[7]。IPCC 的報(bào)告指出，氣候變化會(huì)對(duì)全球大部分區(qū)域的水文循環(huán)造成影響，進(jìn)而導(dǎo)致區(qū)域干旱頻發(fā)，這在亞馬遜地區(qū)已經(jīng)得到證實(shí)[8?10]。因此，對(duì)關(guān)于亞馬遜熱帶森林的干旱敏感性問題產(chǎn)生了一些爭論。一些學(xué)者認(rèn)為，亞馬遜熱帶森林旱季呈碳匯狀態(tài)，雨季則為碳源[11]；反駁者認(rèn)為，熱帶森林長期的維持碳匯功能，但在遭遇干旱事件時(shí)會(huì)逆轉(zhuǎn)為碳源[8,12]。從這些爭論中可以看出，氣候因素會(huì)對(duì)熱帶森林植被碳匯大小產(chǎn)生影響。因此，對(duì)熱帶地區(qū)森林植被碳匯進(jìn)行氣候驅(qū)動(dòng)因素分析極為重要。溫度和降雨是影響生態(tài)系統(tǒng)凈

生產(chǎn)力的最主要的2 個(gè)氣候因素[13?15]，它能通過影響植物的光合作用和呼吸作用[16?17]，進(jìn)而影響到NEP。當(dāng)然還有其他的因素在調(diào)控，如土地利用變化、土壤水分、CO2施肥效應(yīng)以及氮沉降效應(yīng)等[18?21]。例如CO2升高對(duì)生產(chǎn)力的潛在影響可能受到全球變暖和降水模式改變的影響，顯示出了水的可用性和溫度在全球植被光合作用和呼吸作用中承擔(dān)著強(qiáng)大的驅(qū)動(dòng)因素的角色[22]。一方面，熱帶林木長期處于熱穩(wěn)定的環(huán)境下，未來氣溫的持續(xù)升高，很有可能導(dǎo)致熱帶樹木的熱不育狀態(tài)[15]；另一方面，LIU 等[23]的研究表明了降水閾值對(duì)于生態(tài)系統(tǒng)生產(chǎn)和呼吸的調(diào)節(jié)作用。因此，了解熱帶森林生態(tài)系統(tǒng)凈生產(chǎn)力對(duì)于氣溫和降雨的響應(yīng)至關(guān)重要。當(dāng)然，也有研究表明，地形因素會(huì)影響到氣候因素，同一塊區(qū)域較高海拔的地方溫度和降雨會(huì)與低海拔地方呈現(xiàn)不一樣的模式，植被類型也會(huì)有很大的區(qū)別[24?25]。

本研究結(jié)合多種遙感數(shù)據(jù)產(chǎn)品，首先利用時(shí)間序列分析和一元線性回歸的方法，通過季節(jié)性分解得到了NEP 的年際趨勢、年內(nèi)趨勢以及季節(jié)性趨勢。其次基于降水的年內(nèi)變化劃分了干季（1?3 月），濕季（7?9 月），并從時(shí)空層面分析了不同季節(jié)的NEP 變化，以掌握NEP 的時(shí)空動(dòng)態(tài)特征，并采用增強(qiáng)回歸樹的方法，探索了氣候因子（氣溫、降水）、地形因子（海拔、坡度）對(duì)于NEP 的貢獻(xiàn)程度，在此基礎(chǔ)上，分海拔討論不同季節(jié)NEP 趨勢和降水趨勢、溫度趨勢的相關(guān)性，進(jìn)一步確定海南島森林NEP 對(duì)氣候因素的響應(yīng)，旨在為理解全球氣候變化背景下區(qū)域碳循環(huán)研究提供資料。

1 材料與方法

1.1 研究區(qū)概況海南島位于中國的南部，熱帶北緣，地理位置介于108°03′～111°03′E 和18°10′～20°10′N 之間。海南島的氣候類型屬熱帶海島季風(fēng)氣候，長夏無冬，1 年分干濕2 季，全年平均溫度22.5～25.6 ℃，年均降雨量900～2 500 mm[26?27]。海南島擁有較高的森林覆蓋度，其中，熱帶山地雨林是海南島熱帶森林植被中面積最大、分布集中的垂直自然地帶性的植被類型[28]，主要分布在吊羅山、五指山、霸王嶺、尖峰嶺以及鸚哥嶺等600 m以上海拔的山地。而低地森林主要分布在600 m以下中低海拔區(qū)域，包括紅樹林、橡膠林等[29]。

1.2 土地利用數(shù)據(jù)集及海南島森林提取土地利用數(shù)據(jù)集源自中國研制的30 m 空間分辨率全球地表覆蓋數(shù)據(jù)GlobalLand30[30]，共包含2000、2010、2020 年3 期的土地利用數(shù)據(jù)（http://www.globallandcover.com/）?；诤Ｄ蠉u的邊界圖裁出海南島的土地利用數(shù)據(jù)集，然后統(tǒng)計(jì)近20 a 海南島的林地區(qū)域，繪制出穩(wěn)定的林地區(qū)域圖像，即在以上3 組土地利用數(shù)據(jù)中均為林地的區(qū)域會(huì)被提取出來，形成穩(wěn)定的林地分布圖（圖1），再將其分辨率重采樣至0.01°以匹配NEP 數(shù)據(jù)的空間分辨率。

圖 1 研究區(qū)概況圖(a)2020 年海南島土地利用圖；(b)基于3 期的土地利用提取的林地區(qū)域圖；(c)海南島高程圖。

1.3 NEP 數(shù)據(jù)集本研究使用的NEP 是來自全球環(huán)境研究中心（Center for Global Environmental Research，CGER），該產(chǎn)品提供了1999—2018 年空間分辨率為0.01°，時(shí)間分辨率為10 d 的全球凈生態(tài)系統(tǒng)交換（Net Ecosystem Exchange，NEE）數(shù)據(jù)[31]（https://db.cger.nies.go.jp/DL/10.17595/20200227.001.html.en）。初始NEE 數(shù)據(jù)格式被轉(zhuǎn)換為柵格，取相反數(shù)即為本研究所使用的NEP 數(shù)據(jù)。 全年共分為36 幅柵格圖像，再分別將每10 d 的NEP均值數(shù)據(jù)合成月尺度數(shù)據(jù)和年尺度數(shù)據(jù)。

1.4 氣候因子數(shù)據(jù)集氣溫和降水?dāng)?shù)據(jù)來源于國家科技基礎(chǔ)條件平臺(tái)—國家地球系統(tǒng)科學(xué)數(shù)據(jù)共享服務(wù)平臺(tái)?黃土高原科學(xué)數(shù)據(jù)中心(http://loess.geodata.cn)，該數(shù)據(jù)集是根據(jù)CRU 發(fā)布的全球0.5°氣候數(shù)據(jù)集以及WorldClim 發(fā)布的全球高分辨率氣候數(shù)據(jù)集，通過Delta 空間降尺度方案在中國地區(qū)降尺度生成的，并使用496 個(gè)獨(dú)立氣象觀測點(diǎn)數(shù)據(jù)進(jìn)行了數(shù)據(jù)驗(yàn)證[32]。

1.5 Digital Elevation Model（DEM）數(shù)據(jù)集

DEM 數(shù)據(jù)集來自中科院資源環(huán)境科學(xué)與數(shù)據(jù)中心（https://www.resdc.cn/），該數(shù)據(jù)集為基于最新的SRTM V4.1 數(shù)據(jù)經(jīng)整理拼接生成的90 m 的分省數(shù)據(jù)。數(shù)據(jù)采用WGS84 橢球投影。隨后將其空間分辨率采樣至0.01°，并基于1.2 中提取的林地區(qū)域進(jìn)一步提取林地所處的海拔高度。筆者使用ArcGIS 中的空間分析工具根據(jù)海拔計(jì)算出坡度和坡向。

1.6 NEP 和降水的時(shí)間序列分解針對(duì)存在季節(jié)性變化的月度NEP，采用季節(jié)性分解，將其分解為趨勢因子（Trend Component，TC）、季節(jié)性因子（Seasonal Component，SC）和 誤 差 因 子（Error Component，EC）。TC 能準(zhǔn)確把握數(shù)據(jù)的長期變化；SC 能捕捉到數(shù)據(jù)1 年內(nèi)的周期性變化；而EC能反映那些不能被趨勢或季節(jié)效應(yīng)解釋的變化[33]。通過相加模型可以表示為：

式中： NEPt指 每月的NEP，T Ct、S Ct、 E Ct分別指趨

勢、季節(jié)性以及隨機(jī)誤差。

線性濾波器是估計(jì)時(shí)間序列趨勢常用的方法，最常見的線性濾波器之一是滑動(dòng)平均。針對(duì)NEP 月度數(shù)據(jù)，選擇了12 點(diǎn)移動(dòng)平均法，線性濾波器如下：

式中： TCt指去除季節(jié)性效應(yīng)的NEP 趨勢值；1999—2018 共20 年，每年12 個(gè)月，合計(jì)240 個(gè)月，t 從第7 月開始取值。

在獲取趨勢因子之后，通過下式可得到季節(jié)性因子：

式中：S Ct為第t月時(shí)算得的季節(jié)性因子。

式中：Slope 是NEP 逐像元線性回歸方程的斜率；代表年數(shù)，n 為時(shí)間跨度。當(dāng)Slope>0 時(shí)，NEP呈增加趨勢；當(dāng)Slope=0 時(shí)，NEP 基本穩(wěn)定，無明顯變化；當(dāng)Slope<0 時(shí)，NEP 呈減少趨勢。

這個(gè)季節(jié)性因子的估計(jì)也包含了每個(gè)t 時(shí)間的隨機(jī)誤差因子 ECt，通過計(jì)算每個(gè)月季節(jié)性效應(yīng)估計(jì)的平均值并在所有年份重復(fù)此序列來獲取整體的季節(jié)性效應(yīng)。

式中： i 為年份，j 為月份， j 取1~12， N EPij為第i 年第 j 月的平均NEP，N EPj為20 年的j 月平均NEP。

1.7 基于降雨效應(yīng)的干濕季劃分及趨勢分析

基于像元尺度的趨勢分析法能模擬研究區(qū)中每個(gè)柵格單元的變化趨勢，從而反映NEP 變化的方向和速率。計(jì)算公式為：

1.8 基于BRT 的相對(duì)貢獻(xiàn)率分析使用BRT 來評(píng)估1999—2018 年氣候因子與地形因子對(duì)于海南島森林NEP 變化的相對(duì)影響程度。BRT 分析具有容納任何數(shù)據(jù)分布的能力，因此，在分析過程中無需進(jìn)行數(shù)據(jù)的轉(zhuǎn)換。在進(jìn)行BRT 分析之前，對(duì)森林NEP 與氣候因子以及地形因子進(jìn)行皮爾遜相關(guān)分析和顯著性分析。使用R 語言中的GBM包進(jìn)行BRT 分析，將1999—2018 年平均NEP 的逐像元數(shù)據(jù)作為響應(yīng)變量，同一時(shí)期的氣候因子和地形因子作為解釋變量。BRT 的本質(zhì)就是1 個(gè)加法模型，每次建立模型是在之前建立模型損失函數(shù)的梯度下降方向。

1.9 海拔劃分及趨勢相關(guān)性分析基于海南島森林的分布情況[29]，以600 m 海拔為界限，將海南島森林分為低海拔森林和高海拔森林分別研究。筆者計(jì)算了近20 年NEP 趨勢與氣候因子趨勢的相關(guān)性，如下所示：

式中： xi為1999—2018 年溫度或降雨的逐像元趨勢， yi為1999—2018 年NEP 的逐像元趨勢，為像元數(shù)。隨后對(duì)NEP 趨勢、溫度趨勢以及降雨趨勢構(gòu)建了多元線性回歸模型，并以回歸系數(shù)來估計(jì)各因子的貢獻(xiàn)度，進(jìn)而確定不同條件下對(duì)NEP 影響的主導(dǎo)因素。

2 結(jié)果與分析

2.1 NEP 和降雨的時(shí)間序列分解對(duì)1999—2018 年月度NEP 數(shù)據(jù)進(jìn)行時(shí)間序列分析，發(fā)現(xiàn)海南島森林NEP 隨時(shí)間的變化趨勢及季節(jié)性波動(dòng)。通過使用滑動(dòng)平均濾波器，得到了近20 年去除季節(jié)性影響的較為平穩(wěn)的NEP 變化趨勢（圖2）。總體來看，NEP 在20 年間變化起伏波動(dòng)，在2010 年前后出現(xiàn)最大趨勢，而在2016 年出現(xiàn)最低趨勢。筆者進(jìn)一步分析了NEP 年內(nèi)的季節(jié)性效應(yīng)，發(fā)現(xiàn)年內(nèi)NEP 呈現(xiàn)先增后減的季節(jié)影響，1?5 月，NEP 逐漸增大，在5 月達(dá)到最大值，隨后的6?12 月，NEP 逐漸減小。同理，可以得到降雨的年內(nèi)季節(jié)性效應(yīng)。1?9 月，降雨隨月份逐漸增加，在9 月達(dá)到最大值，隨后的10?12 月，降雨逐漸減少?；诖丝梢詫⒑Ｄ蠉u1?3 月劃分為干季，7?9 月劃分為濕季。

2.2 海南島森林NEP 的年度及季節(jié)變化通過計(jì)算每年海南島森林NEP 的均值，得到 1999—2018 年逐年NEP 的年際變化趨勢圖（圖3-a），整體來看NEP 隨時(shí)間呈現(xiàn)不顯著的下降趨勢，有機(jī)碳的變化率為?0.57 g·m?2·a（P>0.05），分段來看，前10 年NEP 呈現(xiàn)不顯著的增長趨勢，有機(jī)碳的變化率為3.3 g·m?2·a（P>0.05）；后10 年NEP 呈現(xiàn)顯著下降趨勢，有機(jī)碳的變化率為?8.4 g·m?2·a（P<0.05）。近20 年NEP 的均值(有機(jī)碳)為483.23 g·m?2·a，其中，有9 年的NEP 高于均值。在2010 年，NEP 達(dá)到最大值，為534.68 g·m?2·a。在2016 年，NEP 有最小值，為439.47 g·m?2·a。

同樣，在1999—2018 年逐年干季和濕季的NEP 趨勢變化圖中（圖3-b）可以看到濕季NEP 有不顯著的下降趨勢，而干季NEP 有不顯著的上升趨勢。濕季NEP 的均值為136.76 g·m?2(有機(jī)碳)，其中有11 年NEP 超過均值，最大值出現(xiàn)在2005年，為156.21 g·m?2，最小值出現(xiàn)在2015 年，為117.10 g·m?2。干季NEP 的均值為102.09 g·m?2，僅有9 年NEP 超過均值，最大值出現(xiàn)在2009 年，為129.54 g·m?2，最小值出現(xiàn)在2005 年，為81.24 g·m?2。干、濕2 季對(duì)比，除2009 年干季NEP 超過濕季NEP，其余年份濕季NEP 的值均高于干季。

圖 2 時(shí)間序列趨勢圖(a)1999—2018 年NEP 月趨勢圖；(b)經(jīng)過時(shí)間序列分解去除季節(jié)性影響后的趨勢圖；(c)NEP 的季節(jié)性變化圖；(d)隨機(jī)誤差；(e)季節(jié)性影響下NEP 的年內(nèi)變化；(f)季節(jié)性影響下降雨的年內(nèi)變化。

圖 3 不同年份和季節(jié)下NEP 的變化(a)NEP 的年際變化；(b)干季和濕季NEP 的變化。

2.3 海南島森林NEP 的空間分布及趨勢變化

通過繪制全年、干季以及濕季的NEP 空間分布圖（圖4-a、b、c）來觀察NEP 在空間方位上的分布。年際的NEP 最高值為825.30 g·m?2·a(有機(jī)碳)，最低值為?54.64 g·m?2·a，其中NEP 的高值集中在海南島中部和東南部，而在海南島北部以及海岸線附近多為NEP 低值聚集。此外，大部分地區(qū)的NEP 均大于0。干季的NEP 最高值為236.03 g·m?2，最低值為?69.89 g·m?2，與年際NEP 相似，干季NEP 的高值多分布在海南島中部五指山、霸王嶺以及尖峰嶺一帶，NEP 低值則聚集在海南島東北部。濕季NEP 的分布較為破碎化，高值和低值相間，其中最高值為211.84 g·m?2，最低值為?19.74 g·m?2，濕季NEP 最高值較干季減少了24.19 g·m?2，但最低值較干季增加了50.15 g·m?2。

圖 4 1999—2018 年NEP 的空間分布及趨勢變化(a)、(b)、(c)分別為年際、干季、濕季NEP 的空間分布；(d)、(e)、(f)分別為年際、干季、濕季NEP 的空間趨勢變化。

應(yīng)用一元線性回歸對(duì)海南島NEP 近20 年的年際趨勢變化進(jìn)行分析，進(jìn)而繪制了NEP 的趨勢變化空間分布圖（圖4-d, 圖4-e, 圖4-f）。年際的趨勢變化如圖4 所示，NEP 呈現(xiàn)顯著下降趨勢的區(qū)域占比為18.49%（P<0.05），NEP 呈現(xiàn)顯著增長趨勢的區(qū)域占比為17.99%（P<0.05），NEP 無顯著變化的區(qū)域占比為63.52%。從分布來看，NEP 顯著增長區(qū)域主要集中在海南島東北部，而顯著下降區(qū)域主要集中在海南島中部。干季的趨勢變化：NEP 呈現(xiàn)顯著下降趨勢的區(qū)域占比為6.39%（P<0.05），NEP 呈現(xiàn)顯著增長趨勢的區(qū)域占比為23.05%（P<0.05），NEP 無顯著變化的區(qū)域占比為70.56%。在空間分布上，干季的NEP 趨勢變化與年際的趨勢變化相似，但是呈顯著下降趨勢的區(qū)域明顯減少。濕季的趨勢變化：NEP 呈現(xiàn)顯著下降趨勢的區(qū)域占比為23.66%（P<0.05），NEP呈現(xiàn)顯著增長趨勢的區(qū)域占比為12.65（P<0.05），NEP 無顯著變化的區(qū)域占比為63.69%。與年際趨勢相比，濕季的NEP 顯著增長趨勢占比降低了5.43%，顯著下降趨勢占比提高了5.17%；與干季趨勢相比，濕季的NEP 顯著增長趨勢占比降低了10.40%，顯著下降趨勢占比提高了17.27%。

2.4 基于BRT 的各因素貢獻(xiàn)度分析通過對(duì)不同季節(jié)的NEP 與氣候因子、地形因子進(jìn)行皮爾遜相關(guān)分析以及顯著性檢驗(yàn)，發(fā)現(xiàn)年際NEP 與氣溫、降雨、海拔以及坡度均有極顯著的相關(guān)性（P<0.01），相關(guān)系數(shù)分別為?0.53、?0.03、0.57、0.50，而與坡向無顯著相關(guān)性（圖5-a）；干季NEP 與氣溫、降雨、海拔以及坡度均有極顯著的相關(guān)性（P<0.01），相關(guān)系數(shù)分別為?0.54、?0.07、0.73、0.66，而與坡向無顯著相關(guān)性（圖5-b）；濕季NEP 與氣溫、降雨、海拔、坡度以及坡向均有極顯著的相關(guān)性（P<0.01），相關(guān)系數(shù)分別為0.07、?0.04、?0.09、?0.10、?0.03（圖5-c）?；谝陨戏治觯x擇氣溫、降雨、海拔和坡度作為自變量因子用于BRT 分析，以便計(jì)算各因素在不同季節(jié)對(duì)NEP 的貢獻(xiàn)程度。在進(jìn)行參數(shù)調(diào)優(yōu)后，在R 語言中使用GBM 包進(jìn)行BRT 分析，得到不同季節(jié)不同因子對(duì)NEP 的相對(duì)貢獻(xiàn)率。如圖6-a 所示，對(duì)年際NEP 影響程度最大的是海拔，其相對(duì)貢獻(xiàn)率為45.46%，其次分別為降雨、氣溫和坡度，相對(duì)貢獻(xiàn)率依次為25.78%、16.17%、12.59%；如圖6-b 所示，對(duì)干季NEP 影響程度最大的是海拔，其相對(duì)貢獻(xiàn)率為40.58%，其次為坡度、降雨和溫度，相對(duì)貢獻(xiàn)率依次為28.13%、16.12%、15.17%；如圖6-c 所示，對(duì)濕季NEP 影響程度最大的是降雨，其相對(duì)貢獻(xiàn)率為30.75%，其次為溫度，海拔和坡度，相對(duì)貢獻(xiàn)率依次為26.77%、21.88%、20.60%。除濕季外，海拔在年際NEP 和干季NEP 中扮演著重要的角色。

圖 5 相關(guān)性熱力圖pre、tem、dem、aspect、slope 分別表示降雨、溫度、海拔、坡向以及坡度；(a)、(b)、(c)分別表示年際、干季和濕季3 個(gè)時(shí)期。

圖 6 基于BRT 的貢獻(xiàn)度分析pre、tem、dem、slope 分別表示降雨、溫度、海拔以及坡度；(a)、(b)、(c)分別表示年際、干季和濕季3 個(gè)時(shí)期。

2.5 不同海拔區(qū)域NEP 的驅(qū)動(dòng)因子分析為進(jìn)一步說明海南島森林NEP 的主要影響因素，分別對(duì)高海拔區(qū)域和低海拔區(qū)域進(jìn)行了NEP 趨勢和氣候因子（氣溫、降雨）趨勢的相關(guān)性分析。結(jié)果顯示，在海拔600 m 以上的森林生態(tài)系統(tǒng)中，全年的NEP 趨勢和濕季的NEP 趨勢與降雨趨勢均有顯著的負(fù)相關(guān)性，相關(guān)系數(shù)分別為?0.11 和?0.19（圖7-d, 圖7-f），而干季的NEP 趨勢則與降雨趨勢沒有顯著的相關(guān)性（圖7-e）。不論是年際的NEP趨勢還是干季、濕季的NEP 趨勢均與溫度的趨勢沒有顯著的相關(guān)性（圖7-a, 圖7-b, 圖7-c）。而在

海拔600 m 以下的森林生態(tài)系統(tǒng)中，則呈現(xiàn)出十分不同的模式。年際的NEP 趨勢對(duì)于溫度趨勢有微弱的顯著負(fù)相關(guān)性，相關(guān)系數(shù)為?0.02（圖8-a），對(duì)于降雨趨勢有極顯著的正相關(guān)性，相關(guān)系數(shù)為0.46（圖8-d）；干季的NEP 趨勢與溫度趨勢沒有顯著的相關(guān)性（圖8-b），與降雨趨勢有極顯著的負(fù)相關(guān)性，相關(guān)系數(shù)為?0.06（圖8-e）；濕季的NEP 趨勢與溫度趨勢、降雨趨勢均有極顯著的正相關(guān)關(guān)系，相關(guān)系數(shù)分別為0.11 和0.30（圖8-c,f）。如圖9 所示，多元線性回歸結(jié)果表明，在干季600 m 以下的森林生態(tài)系統(tǒng)中，降水趨勢對(duì)NEP 趨勢有顯著的負(fù)貢獻(xiàn)，貢獻(xiàn)率為?53%。在濕季600 m 以下的森林生態(tài)系統(tǒng)中，溫度趨勢對(duì)NEP 趨勢有顯著的正貢獻(xiàn)，貢獻(xiàn)率為90%。

圖 7 海拔600 m 以上森林NEP 趨勢與溫度、降雨趨勢的相關(guān)性散點(diǎn)圖(a)、(b)、(c)分別表示年際、干季、濕季與溫度趨勢的相關(guān)性；(d)、(e)、(f)分別表示年際、干季、濕季與降雨趨勢的相關(guān)性。

圖 8 海拔600 m 以下森林NEP 趨勢與溫度、降雨趨勢的相關(guān)性散點(diǎn)圖(a)、(b)、(c)分別表示年際、干季、濕季與溫度趨勢的相關(guān)性；(d)、(e)、(f)分別表示年際、干季、濕季與降雨趨勢的相關(guān)性。

圖 9 基于多元線性回歸的不同時(shí)期不同海拔下氣候因子的貢獻(xiàn)度分析圖

紅色五角星表示方程具有顯著性，黑色星號(hào)表明回歸系數(shù)或殘差通過顯著性P=0.05 水平檢驗(yàn)。

3 討 論

3.1 海南島森林穩(wěn)定的碳匯功能海南島森林NEP 表現(xiàn)出了季節(jié)性變化，但是與亞馬遜地區(qū)不同的是，亞馬遜地區(qū)在季節(jié)性變換時(shí)會(huì)存在碳動(dòng)態(tài)的轉(zhuǎn)變，如碳源和碳匯的相互轉(zhuǎn)換[8,11?12,34?35]，而海南島熱帶森林不論是從年際角度出發(fā)，還是從干季、濕季角度出發(fā)，均呈現(xiàn)穩(wěn)定碳匯狀態(tài)。所不同的是，海南島森林NEP 在濕季時(shí)普遍要比干季高。地理位置差異所帶來的氣候類型差別可能是海南島熱帶森林和亞馬遜熱帶雨林呈現(xiàn)不同碳動(dòng)態(tài)格局的原因。亞馬遜地區(qū)屬熱帶雨林氣候，長年高溫多雨，對(duì)干旱的敏感性更高[36]；而海南島則屬熱帶海島季風(fēng)氣候，受季風(fēng)影響分明顯的干季和濕季。相對(duì)固定的季節(jié)模式，使得海南島森林固碳動(dòng)態(tài)雖有季節(jié)性變化，但總體趨勢保持穩(wěn)定。更為重要的一點(diǎn)是，不只是氣候變化的影響，農(nóng)業(yè)擴(kuò)張所帶來的森林砍伐，火災(zāi)和干旱之間的相互作用也是造成亞馬遜森林碳損失的原因[36?37]。相比而言，海南島的森林保存的較為完整。此外，本研究采用基于多期土地利用遙感數(shù)據(jù)提取穩(wěn)定林地的方法，一定程度上規(guī)避了森林變化或者損失而帶來的凈生態(tài)系統(tǒng)生產(chǎn)力的損失。

3.2 不同海拔不同季節(jié)的NEP 與氣候因素的相關(guān)性分析從趨勢的相關(guān)性上來看，相比氣溫，降雨與NEP 有著更為顯著的相關(guān)性。但是在不同海拔下，降雨與NEP 的相關(guān)性表現(xiàn)出了相反的模式。在海拔600 m 以下，年際、濕季降雨趨勢分別與年際、濕季NEP 趨勢呈現(xiàn)出顯著的正相關(guān)，但是在海拔600 m 以上，又呈現(xiàn)出截然相反的顯著負(fù)相關(guān)性結(jié)果。群落結(jié)構(gòu)和生態(tài)系統(tǒng)過程往往隨海拔梯度變化[24]，一般而言對(duì)海拔的響應(yīng)都是由于溫度變化所驅(qū)動(dòng)，但也不絕對(duì)，降雨等因素也會(huì)隨海拔梯度變化[24,38]。海拔梯度下，海南島降雨趨勢的變化不是十分明顯，所以這種結(jié)果可能是因?yàn)槿郝浣Y(jié)構(gòu)發(fā)生變化。海南島海拔600 m 以上保存較為原始的熱帶山地雨林，在年際和濕季雨量較為充沛的階段，隨著降雨趨勢的增加，附生植物旺盛生長會(huì)導(dǎo)致NEP 下降。而在低海拔區(qū)域，由于人為干擾，林型更為復(fù)雜，人工林的固碳能力相較高海拔原始林更為強(qiáng)大。

3.3 降水與溫度的主導(dǎo)性分析僅從氣候因素考慮，在旱季時(shí)，降雨對(duì)低海拔森林NEP 變化有顯著的負(fù)向作用，而在濕季時(shí)，溫度對(duì)低海拔森林NEP 變化有積極作用。氣溫和降雨是影響生態(tài)系統(tǒng)生產(chǎn)力的主要的因素[39]，但是隨著全球變化加速，陸地生態(tài)系統(tǒng)的碳交換對(duì)于氣候變化的響應(yīng)和反饋仍舊存在不確定性[40?41]。尤其在熱帶森林生態(tài)系統(tǒng)，降雨量在年內(nèi)的波動(dòng)導(dǎo)致季節(jié)性干旱的產(chǎn)生，干旱事件對(duì)熱帶森林生態(tài)系統(tǒng)樹木的生長以及生態(tài)系統(tǒng)的功能有著嚴(yán)重的影響[42]。低海拔區(qū)域更多為人工林，以橡膠林為例，干旱會(huì)導(dǎo)致其落葉甚至死亡，光合作用受到抑制[43]。另一方面，旱季的降雨會(huì)增加土壤水分[44]，進(jìn)一步影響了土壤呼吸[45]。雨季時(shí)，溫度呈現(xiàn)出顯著的正向作用。降水的充沛使得水分并不是NEP 的主要限制因子，而此時(shí)的溫度卻有著積極的調(diào)節(jié)能力，與NEP 變化呈現(xiàn)出一致性。研究表明，溫度的升高是利于生態(tài)系統(tǒng)呼吸[16]、生態(tài)系統(tǒng)總產(chǎn)量[46]以及凈碳吸收[47]。可能在雨季時(shí)，由于不再受到水分條件的限制，海南島森林植被得到較大恢復(fù)，處于生長階段的旺盛期，此時(shí)溫度對(duì)于整個(gè)生態(tài)系統(tǒng)的光合作用影響遠(yuǎn)遠(yuǎn)大于呼吸作用。

4 結(jié) 論

近20 年海南島森林NEP 的年際變化波動(dòng)較大，存在明顯的季節(jié)性。濕季的NEP 高于干季的NEP?？臻g上，NEP 呈顯著增長趨勢主要分布在海南島東北部，顯著下降趨勢主要分布在海南島中部。不同季節(jié)NEP 與氣候因子的關(guān)系有所差異，在干季低海拔區(qū)域，降雨對(duì)NEP 有顯著的負(fù)向貢獻(xiàn)，貢獻(xiàn)度為?53%（P<0.05）;在濕季低海拔區(qū)域，溫度對(duì)NEP 有顯著的正向貢獻(xiàn)，貢獻(xiàn)度為90%（P<0.05）。