海南島橡膠林區(qū)域不同SIF產(chǎn)品的差異性分析及其對GPP估算的影響

安映荷，張潤卿，劉文杰,2，王 鵬，章 杰，烏 蘭，孫仲益

（1.海南大學 生態(tài)與環(huán)境學院，?？?570100; 2.海南省農(nóng)林環(huán)境過程與生態(tài)調(diào)控重點實驗室，?？?570100）

日光誘導(dǎo)葉綠素熒光（Sun-Induced Chlorophyll Fluorescence, SIF）是受太陽輻射而激發(fā)的葉綠素重新發(fā)射光子回到基態(tài)時產(chǎn)生的1種光信號，光譜范圍在650～800 nm，可在紅光和近紅外光譜區(qū)被監(jiān)測到[1]。植被光合作用分為光反應(yīng)和暗反應(yīng)，而SIF是1種受光能驅(qū)動的信號，發(fā)生在光反應(yīng)過程中，因此與光合作用聯(lián)系緊密[2]，并且與其他衛(wèi)星遙感植被指數(shù)相比，SIF信號更接近植被實際的光合作用，能更有效地反映植被生長狀況[3]。目前，關(guān)于SIF的驗證和應(yīng)用研究主要依托于模型建立[4]或高光譜傳感器測量[5-6]等工具，并已取得顯著進展。其中，基于地面實驗及高空的衛(wèi)星遙感反演方式主要證實了SIF與植被光合作用的關(guān)系[7-8]?，F(xiàn)有的大部分SIF遙感產(chǎn)品通過衛(wèi)星傳感器反演得到，這些遙感產(chǎn)品被廣泛應(yīng)用在陸地總初級生產(chǎn)力（Gross Primary Production, GPP）估算[9-11]、干旱監(jiān)測[12-14]、植被物候反演[15-16]及碳平衡觀測[17]等方面。隨著衛(wèi)星遙感技術(shù)的不斷發(fā)展，歐洲航天局（European Space Agency, ESA）首個專門為SIF探測設(shè)計的衛(wèi)星“熒光探測器”（Fluorescence Explorer, FLEX）于2022年發(fā)射，該衛(wèi)星可以提供300 m精細空間分辨率和全波段的全球植被SIF產(chǎn)品[18]，但可用的時間序列有限。SIF作為光合作用有力的探測工具，為GPP提供了1種新的估算方法[19]。因此，近年來很多學者通過日本溫室氣體觀測衛(wèi)星（Greenhouse Gases Ob-serving Satellite, GOSAT）[20]、歐洲METOP第二代臭氧業(yè)務(wù)監(jiān)測傳感器（Global Ozone Monitoring Experiment-2, GOME-2）[21]、掃描成像吸收光譜儀（SCanning Imaging Absorption SpectroMeter for Atmospheric CHartograph Y, SCIAMACHY）[22]、軌道碳觀測衛(wèi)星2號（Orbiting Carbon Observatory-2,OCO-2）[23]及中國首顆二氧化碳探測衛(wèi)星TanSat[24]進行SIF的反演，并構(gòu)建了一系列全球SIF衛(wèi)星遙感產(chǎn)品，為SIF遙感產(chǎn)品在估算GPP等方面提供有效數(shù)據(jù)來源。其中以CSIF（Continuous Solar Induced Fluorescence）[25]、GOSIF（Global OCO-2 of SIF）[26]、SIFOCO2-1696（High Resolution Global Contiguous SIF of OCO-2）[27]、SIF005（Harmonized long-term SIF）[28]、SIFLUE（Light Use Efficiency SIF）[29]5種產(chǎn)品最具有代表性。

現(xiàn)有的SIF遙感產(chǎn)品在全球尺度上應(yīng)用與驗證效果雖好，但是在區(qū)域的普適性上還有待證明[20,30-31]。中國位于亞洲東部和太平洋西岸，生態(tài)系統(tǒng)類型多樣，大部分地區(qū)處于溫帶大陸，少部分地處熱帶，楊鳳珠等[32]以中國生態(tài)系統(tǒng)光譜觀測網(wǎng)絡(luò)的4個農(nóng)田站SIF觀測數(shù)據(jù)為參考，對不同SIF遙感產(chǎn)品進行驗證，發(fā)現(xiàn)各產(chǎn)品在空間格局上保持高度一致性，在年均值變化趨勢上面，CSIF和GOSIF的表現(xiàn)最好；但對于中國熱帶地區(qū)SIF遙感產(chǎn)品適用情況的表述并不明確。目前，在中國乃至全世界范圍內(nèi)，鮮有針對常年高云覆蓋、蒸發(fā)量強及植被群落結(jié)構(gòu)復(fù)雜的熱帶地區(qū)的SIF遙感產(chǎn)品進行適用性分析。熱帶地區(qū)是重要的碳匯增長點，熱帶森林作為覆蓋率最高、群落結(jié)構(gòu)最復(fù)雜、碳匯能力最強的生態(tài)系統(tǒng)，對生物地球循環(huán)和大氣功能都有著重大影響[33]，但也正因如此，對熱帶森林開展研究相對也存在更多的困難。橡膠樹（Hevea brasiliensis）是全球熱帶地區(qū)重要的經(jīng)濟林木和人工種植生態(tài)系統(tǒng)之一[34]，其物種多樣性在一定尺度上低于熱帶雨林[35-36]，生態(tài)系統(tǒng)的結(jié)構(gòu)和功能與熱帶雨林相比較也較簡單[37]，是研究熱帶森林的理想突破口。橡膠林具有水土保持、固定CO2和維持土壤肥料等生態(tài)功能[38]，同時，作為熱帶地區(qū)的經(jīng)濟林木，又具有一定的實際生產(chǎn)意義[39]。此外，海南島橡膠林種植區(qū)域占全島森林總面積的1/4，是我國第二大天然橡膠種植區(qū)。綜上所述，研究在海南島橡膠林GPP估算中的不同SIF遙感產(chǎn)品的適用性對比，對于重要農(nóng)產(chǎn)品保護、實際生產(chǎn)和生態(tài)環(huán)境保護都有重要實踐意義。同時，為SIF遙感產(chǎn)品在熱帶地區(qū)GPP估算中提供理論參考和借鑒意義。本研究選擇海南島橡膠林為研究對象，以2001—2019年海南島橡膠林生態(tài)系統(tǒng)GPP數(shù)據(jù)為參考，從不同的時間、空間角度通過相對系數(shù)和絕對系數(shù)等指標，對比了CSIF、GOSIF、SIFOCO2-1696、SIF005、SIFLUE5種產(chǎn)品的差異性及與海南島橡膠林生態(tài)系統(tǒng)GPP的時空一致性，旨在分析5種SIF遙感產(chǎn)品在海南島橡膠林的時空格局差異，并通過空間分布一致性、季節(jié)變化一致性和與GPP的相關(guān)性等相關(guān)指標的計算，探討不同SIF遙感產(chǎn)品對海南島橡膠林GPP估算中的可能帶來的影響。

1 材料與方法

1.1 研究區(qū)概括海南島位于（108°03′—111°03′）E和（18°10′—20°10′）N之間，西臨北部灣，東瀕南海，是我國第二大島，全島面積約為3.54萬km2。海南島屬于熱帶季風氣候，年均溫25℃，年降水量在1 600 mm以上，降水量充沛，日照充足，受冷空氣影響小。海南島的森林覆蓋率為53.51%，主要分布有熱帶濕潤雨林、常綠季雨林、山地雨林、落葉季雨林和熱帶稀樹草原，其余為人工種植的經(jīng)濟林木，如橡膠林[40]。島內(nèi)的橡膠林種植面積僅次于云南省，是我國第二大橡膠種植區(qū)。

1.2 數(shù)據(jù)來源

1.2.1 SIF衛(wèi)星遙感產(chǎn)品使用的5款SIF遙感產(chǎn)品見表1。（1） CSIF。該產(chǎn)品利用中分辨率成像光譜儀（Moderate-resolution Imaging Spectroradiometer, MODIS）的地表反射率和OCO-2傳感器的SIF數(shù)據(jù)，生成了2款全球空間連續(xù)SIF產(chǎn)品，分別用于晴空瞬時條件與全天條件。其中，全天條件的SIF數(shù)據(jù)集具有較強的空間、季節(jié)和年際動態(tài)變化[25]。為便于與其他產(chǎn)品進行對比，本研究選擇晴空瞬時條件（CSIFclear-inst）的數(shù)據(jù)產(chǎn)品作為研究對象。（2） GOSIF。該產(chǎn)品是Li等[26]基于數(shù)據(jù)驅(qū)動的方法開發(fā)的1種新的全球SIF數(shù)據(jù)集。GOSIF的預(yù)測模型基于OCO-2探測器、MODIS遙感數(shù)據(jù)和氣象再分析數(shù)據(jù)，具有更精細的空間分辨率。本研究選擇了時間分辨率為逐月的、有效時間范圍為2000—2020年的數(shù)據(jù)產(chǎn)品作為研究對象。（3） SIFOCO2-1696。該產(chǎn)品是利用機器學習方法在充分考慮到植物在空間和時間上的不同生理特征基礎(chǔ)上，對生物群落和時間進行分層訓(xùn)練，得到的全球空間連續(xù)SIF產(chǎn)品。該產(chǎn)品準確地保存了OCO-2軌道的時空變化，并與獨立機載SIF測量結(jié)果相比較，具有很強的一致性[27]。（4） SIF005。該產(chǎn)品為德國地球科學研究中心（GFZ）開發(fā)，綜合SCIAMACHY與GOME-2的SIF產(chǎn)品，基于時間和空間分層訓(xùn)練神經(jīng)網(wǎng)絡(luò)的機器學習算法，解釋不同的生物生理和物理因素與SIF在時空變化上的關(guān)系，將原始數(shù)據(jù)空間分辨率細化至0.05°，并利用累積分布函數(shù)匹配技術(shù)進行校正，得到最終的SIF005產(chǎn)品。（5） SIFLUE。該產(chǎn)品基于2版GOME-2衛(wèi)星傳感器的遙感數(shù)據(jù)（分別用JJ和PK表示）[22]，利用構(gòu)建光能利用效率（LUE）模型的方法獲取到的SIF數(shù)據(jù)集[29]。該數(shù)據(jù)集根據(jù)建立在SIF與高空間分辨率遙感解釋變量間的半經(jīng)驗關(guān)系，利用來自O(shè)CO-2的SIF觀測數(shù)據(jù)進行獨立驗證，選擇了最優(yōu)解釋變量集，完成SIF產(chǎn)品生產(chǎn)。

表1 SIF衛(wèi)星遙感數(shù)據(jù)集信息

1.2.2 海南島橡膠林GPP數(shù)據(jù)以農(nóng)業(yè)農(nóng)村部儋州熱帶作物科學觀測實驗站（又稱儋州橡膠林通量站；該站點位于橡膠林區(qū)域內(nèi)，裝有開路渦度相關(guān)系數(shù)（OPEC）和植物冠層垂直微氣候剖面系統(tǒng)，已連續(xù)觀測10余年[41]）觀測數(shù)據(jù)為基礎(chǔ)，本地化的改進EC-LUE（eddy covariance light use efficiency）模型并作為估算工具[42]，獲得2001—2020年逐月海南島橡膠林生態(tài)系統(tǒng)GPP數(shù)據(jù)，其空間分辨率為1 000 m。原始EC-LUE模型[43]參數(shù)包括歸一化差值植被指數(shù)（NDVI）、光合作用輻射（PAR）、氣溫以及顯熱通量與潛熱通量的比率，本研究使用的EC-LUE模型綜合了飽和水汽壓差（VPD）對GPP的限制影響[44]。

1.3 研究方法

1.3.1 技術(shù)路線技術(shù)路線如圖1所示。（1）對來源不同的5種SIF產(chǎn)品進行統(tǒng)一時間分辨率的預(yù)處理，為消除隨機誤差的影響，將CSIF和SIFLUE2款產(chǎn)品分別按每4 d和8 d的最大值統(tǒng)一為逐月的數(shù)據(jù)集，其余數(shù)據(jù)集按每月均值統(tǒng)一為相同時間分辨率的數(shù)據(jù)集；（2）對5種SIF產(chǎn)品和海南島橡膠林GPP原始數(shù)據(jù)進行掩膜提取，保留包含65%以上面積為橡膠林的像元，確定產(chǎn)品選擇的空間范圍；（3）計算5種SIF產(chǎn)品有效時間范圍內(nèi)的多年均值、月均值和趨勢值；各產(chǎn)品的月均值進行歸一化處理，以消除傳感器不同波段的影響，繼而計算季節(jié)性強度指數(shù)和分析季節(jié)性變化；（4）逐像元分別計算5種SIF產(chǎn)品與海南島橡膠林GPP的解釋系數(shù)，確定空間相關(guān)關(guān)系；（5）評估5種SIF產(chǎn)品的時空格局差異性以及與海南島橡膠林GPP的一致性，最終完成不同SIF產(chǎn)品對GPP估算影響分析。具體評估方法與指標見表2。

圖1 技術(shù)路線圖

表2 SIF遙感數(shù)據(jù)集差異性評估指標

1.3.2 空間一致性檢驗方法以各產(chǎn)品在有效時間段內(nèi)的年均值進行空間分布特征對比分析；同時，計算各產(chǎn)品逐像元的變化趨勢，分析各產(chǎn)品空間上是否具有一致性。

趨勢計算公式如下：

式中，Slope為SIF年均值變化趨勢線的斜率，即SIF的變化趨勢；n代表每種數(shù)據(jù)集的有效時間段；i代表年份（i=1，2，3，……，n）；SIFi表示第i年的年均值。當Slope大于0時，則呈現(xiàn)上升趨勢，反之則呈下降趨勢。

最后，通過計算各產(chǎn)品與海南島橡膠林GPP之間的決定系數(shù)，對比分析5種SIF產(chǎn)品在空間上與GPP之間的擬合程度，以此來判斷各產(chǎn)品在GPP估算中的適用性。

1.3.3 時間一致性檢驗方法很多研究理論與實踐證明，GPP與SIF之間在時間尺度上呈強線性相關(guān)[45]。本研究旨在探討SIF產(chǎn)品在海南島橡膠林GPP估算適用性差異，為能更好地反映SIF值與橡膠樹生長的季節(jié)變化一致性情況，因此只選擇分析SIF數(shù)據(jù)集時間序列的季節(jié)性變化。利用標準化處理后的月均值，從趨勢值、多年每月均值、與GPP的相關(guān)系數(shù)及季節(jié)性指數(shù)4個方面評估不同SIF數(shù)據(jù)集時間序列在季節(jié)上的變化特征，以此作為時間一致性的分析方法。

其中，季節(jié)性指數(shù)法是以時間序列含有季節(jié)性周期變動的特征，計算描述該變動的季節(jié)變動指數(shù)，計算公式如下[46]：

式中，Si為第i月的季節(jié)性指數(shù)；n為總的年數(shù)；SIF（i,j）為第j年第i月的SIF值；12表示1年當中的12個月。每個月SIF均值的平均季節(jié)性指數(shù)為100%，當某個月的Si超過100%時，表示該月的SIF均值在季節(jié)尺度上達到高峰期，Si最大則說明對應(yīng)月份的SIF均值為峰值，反之，則為最低值。

2 結(jié)果與分析

2.1 SIF產(chǎn)品的差異性分析

2.1.1 多年均值空間分布特征五種SIF遙感數(shù)據(jù)集的多年均值空間分布圖如圖2所示。5種SIF產(chǎn)品表現(xiàn)出相似的空間分布格局，北部多為低值，而南部多為高值，這與海南島橡膠林植被覆蓋率由北向南逐漸升高的分布相一致。其中，CISF（圖2-a）、GOSIF（圖2-b）和SIF005（圖2-c）3款產(chǎn)品從數(shù)據(jù)分布上看，具有較高的完整性；SIFLUE（圖2-e）產(chǎn)品在東南部表現(xiàn)出數(shù)據(jù)缺失，但高低值的分布特征與前三者相同。

圖2 海南島橡膠林5種SIF產(chǎn)品多年均值空間分布圖

2.1.2 趨勢分布特征5款SIF遙感數(shù)據(jù)集多年均值的趨勢空間分布特征如圖3所示。僅有CSIF（圖3-a）和SIFLUE（圖3-e）2款產(chǎn)品在島南及西南部的零星橡膠林區(qū)域表現(xiàn)出降低趨勢，其余3款產(chǎn)品均呈現(xiàn)出全島范圍的增加趨勢。在趨勢的高低分布特征上，5款產(chǎn)品表現(xiàn)一致，高值普遍分布在海南島北部，其中SIFLUE（圖3-e）數(shù)據(jù)集在島中部也有高值分布。經(jīng)統(tǒng)計發(fā)現(xiàn)，SIFOCO2-1696（圖3-d）數(shù)據(jù)集趨勢值均值最高；GOSIF（圖3-b）的趨勢值均值最低。

圖3 海南島橡膠林5種SIF產(chǎn)品多年均值趨勢值空間分布圖

2.1.3 季節(jié)性特征5款SIF遙感數(shù)據(jù)集年內(nèi)月均值季節(jié)性變化如圖4所示，所有數(shù)據(jù)均經(jīng)各自標準化處理后。其中，5種產(chǎn)品的月均值均表現(xiàn)出先升高再降低的變化趨勢，說明產(chǎn)品在季節(jié)性變化特征上相一致，峰值出現(xiàn)在5—8月，低值出現(xiàn)在1月和2月，這與海南島橡膠林的生長狀況一致。此外，如表3所示，各產(chǎn)品在4—10月期間的季節(jié)性指數(shù)普遍大于100%，表示在這期間的SIF值處于高峰期。5款產(chǎn)品分別在6月、8月、8月、6月和5月的季節(jié)性指數(shù)最高，說明各產(chǎn)品在5—8月中達到峰值。

圖4 海南島橡膠林5種SIF產(chǎn)品每年月均值季節(jié)性變化圖

表3 5種SIF遙感產(chǎn)品每月季節(jié)性指數(shù)

2.1.4 季節(jié)性趨勢特征5款SIF遙感數(shù)據(jù)集的各月趨勢變化如圖5所示。CSIF、GOSIF和SIF0053款產(chǎn)品各季節(jié)均呈現(xiàn)增加的趨勢（圖5右側(cè)坐標軸），而SIFOCO2-1696與SIFLUE則表現(xiàn)出不規(guī)則的降低趨勢（圖5左側(cè)坐標軸）。SIFOCO2-1696數(shù)據(jù)集在12月的負值最大，說明該數(shù)據(jù)集在12月的SIF值明顯降低；SIF005數(shù)據(jù)集在3月有明顯的增加趨勢，此時的趨勢值為0.026；CSIF、GOSIF和SIFLUE3款數(shù)據(jù)集的趨勢變化較為穩(wěn)定，而其余2款數(shù)據(jù)的趨勢變化的突變性較明顯。

圖5 海南島橡膠林5種SIF產(chǎn)品逐月趨勢變化

2.2 SIF 對 GPP 估算的影響分析

2.2.1 空間影響分析5款SIF遙感產(chǎn)品數(shù)據(jù)集與海南島橡膠林GPP的多年均值的決定系數(shù)空間分布特征如圖6所示。通過分析R2的空間分布情況得知，CSIF（圖6-a）、GOSIF（圖6-b）和SIF005（圖6-c）3款產(chǎn)品的分布特征表現(xiàn)出一致性，均呈現(xiàn)出北高南低的分布特征，能夠較好地捕捉到海南島橡膠林生態(tài)系統(tǒng)GPP的空間分布格局，可用于估算GPP空間分布特征。而在此結(jié)果中，SIFOCO2-1696因可用數(shù)據(jù)集太少，無法計算與GPP的決定系數(shù)；SIFLUE（圖6-d）數(shù)據(jù)集的計算結(jié)果分布特征的表現(xiàn)次于CSIF和GOSIF。

圖6 海南島橡膠林GPP與5種SIF產(chǎn)品多年均值的決定系數(shù)空間分布圖

2.2.2 時間影響分析5款SIF遙感數(shù)據(jù)集與海南島橡膠林GPP的相關(guān)系數(shù)季節(jié)性變化特征如圖7所示。通過分析可知，CSIF和GOSIF 2款產(chǎn)品的多年月均值和GPP的多年月均值，在季節(jié)變化上呈現(xiàn)高度一致性，在4月GPP與SIF之間的相關(guān)性最高（R=0.77（P＜0.01）和R=0.74（P＜0.01））；SIF005和SIFLUE在前4個月的變化與前兩者一致，5—12月則呈現(xiàn)出不規(guī)則變化；SIFOCO2-1696數(shù)據(jù)集與GPP的季節(jié)相關(guān)性較差，且季節(jié)變化不穩(wěn)定。

圖7 海南島橡膠林GPP與5種SIF產(chǎn)品多年均值的相關(guān)性系數(shù)季節(jié)變化圖

2.3 SIF產(chǎn)品的GPP估算適用性5款SIF產(chǎn)品在空間分析、長期趨勢分析和對GPP估算的影響3個方面的總結(jié)結(jié)果如表4所示，5款產(chǎn)品中只有GOSIF無論在空間分析還是長期趨勢分析中的結(jié)果都優(yōu)于其他4款產(chǎn)品（圖2，圖5），其次是CSIF產(chǎn)品。在海南島橡膠林區(qū)域適用情況最差的是SIFLUE產(chǎn)品，這是由于該產(chǎn)品在該地區(qū)的數(shù)據(jù)完整性較差。而針對海南島橡膠林區(qū)域GPP估算的適用性可以看出，GOSIF和CCSIF 2款產(chǎn)品與橡膠林GPP無論在空間還是長期趨勢上都保持高度一致，適用情況較好。

表4 5種SIF遙感產(chǎn)品適用情況

3 討 論

3.1 不確定性分析CSIF、GOSIF、SIF005、SIFOCO2-1696和SIFLUE5款產(chǎn)品從空間和時間兩個角度對相對系數(shù)、決定系數(shù)、均值、趨勢值等指標進行評估，不確定性主要來源于：（1） 5款產(chǎn)品的時間序列長度不同。5款產(chǎn)品來自不同的傳感器，估算方法不同，每款數(shù)據(jù)的有效時間范圍也不同，其中，GOSIF擁有最長的時間序列長度，為21 a；而SIFOCO2-1696的時間序列長度最短，只有4 a，因此，統(tǒng)一有效時間序列長度對本研究沒有實際意義，而這種不一致性會導(dǎo)致在季節(jié)性變化分析上的不確定性，因為SIFOCO2-1696可用年份較少，使用均值計算與其他長時間序列的數(shù)據(jù)相比，缺乏代表性。（2）空間覆蓋度的完整性。從空間覆蓋度的角度來說，5款數(shù)據(jù)也缺少一致性，其中，SIFLUE的數(shù)據(jù)缺失最為嚴重，空間覆蓋率較低，這是因為SIFLUE的原數(shù)據(jù)集較為稀疏[29]，并且缺少解釋變量的數(shù)據(jù)，由此產(chǎn)生的細化產(chǎn)品雖然空間分辨率提高了，但仍存在數(shù)據(jù)稀疏的差距，在島嶼或半島上就會出現(xiàn)空白，即鄰近網(wǎng)格單元數(shù)量不足，無法細化。因此，會對SIFLUE數(shù)據(jù)集在空間分布特征和季節(jié)變化特征上的分析帶來不確定性。其次，SIF產(chǎn)品有效柵格選取過程中以橡膠林占面積超過65%作為依據(jù)，地物的混合效應(yīng)也是引入不確定性的一個原因；海南島橡膠林的種植面積在近30 a內(nèi)呈現(xiàn)出“北增南減”的空間趨勢[47]，北部橡膠林連片程度更高，SIF數(shù)據(jù)能更好的反映GPP的變化情況。（3）SIF數(shù)據(jù)在熱帶地區(qū)的局限性。IF數(shù)據(jù)與其他光學遙感植被指數(shù)相比，受到云、雪或冰覆蓋的影響小[48]，但熱帶地區(qū)的氣候條件復(fù)雜，云層較厚、海鹽氣溶膠濃度高，雖沒有研究證明SIF產(chǎn)品是否與其他植被指數(shù)一樣，在熱帶地區(qū)的適用性與普遍得出的理論存在偏差[49]，但熱帶地區(qū)的復(fù)雜環(huán)境條件對衛(wèi)星遙感產(chǎn)品反演帶來的不確定性無法忽視。（4）參考值的不確定性。本研究的GPP數(shù)據(jù)來源于改進型EC-LUE的GPP估算產(chǎn)品，站點尺度的模擬結(jié)果與觀測值相比RMSE達到2.00 gC·m-2·d-1[44]，并且該數(shù)據(jù)集在空間尺度升級過程中驅(qū)動數(shù)據(jù)的不確定性會進一步增加全島尺度的GPP準確性。

3.2 SIF產(chǎn)品的差異性分析本研究選擇的5款SIF遙感數(shù)據(jù)集的差異性來源主要是衛(wèi)星傳感器和算法不同。OCO-2衛(wèi)星傳感器發(fā)射于2014年7月2日，可用波段為757～775 nm[23]，由此傳感器反演出的3款數(shù)據(jù)CSIF、GOSIF和SIFOCO2-1696的波長為757 nm，CSIF產(chǎn)品相對于原始數(shù)據(jù)，缺乏生理和物理意義[25]，GOSIF則充分考慮了植被、氣候和土地覆蓋條件3種類型的SIF數(shù)據(jù)解釋變量[26]，SIFOCO2-1696與CSIF的傳感器來源相同，但存在是否包含生物和物理意義的差異，前者因算法的不同，所以具有合理的生理和物理意義，因此具有比CSIF更強的季節(jié)性。SIF005和SIFLUE分別來自SCIAMACHY/GOME-2和GOME-2。SCIAMACHY和GOME-2兩者的載荷設(shè)置類似，光譜分辨率為0.48 nm。SIF005數(shù)據(jù)集是利用SCIAMACHY和GOME-2月值SIF產(chǎn)品通過降低空間分辨率得到的產(chǎn)品，實現(xiàn)了2種傳感器的融合；SIFLUE是GOME-2在空間細化尺度的產(chǎn)品，基于光能利用效率模型，選擇植被、水分和溫度參數(shù)作為SIF的解釋變量，提高了該數(shù)據(jù)的降尺度產(chǎn)品的精度[50]。

3.3 SIF產(chǎn)品的GPP估算適用性從年均值的時空角度分析，CSIF和GOSIF在空間分布特征分析中比其余3種產(chǎn)品更適合應(yīng)用于海南島橡膠林GPP估算，但是CSIF有效時間序列短，不適合長期趨勢分析。由圖5可知，CSIF在3月的趨勢值最低，增加趨勢不明顯，且季節(jié)性變化沒有GOSIF穩(wěn)定；這是因為本研究中選擇的是晴空瞬時條件下的CSIF產(chǎn)品，相比較另一款CSIF產(chǎn)品來說，在季節(jié)變化上缺少與其他SIF產(chǎn)品的一致性[51]。另外，就CSIF產(chǎn)品而言，該產(chǎn)品的反演過程中缺少植被生理和環(huán)境物理因素，因此無法有效體現(xiàn)光合生理狀態(tài)[25]，從而季節(jié)性變化不明顯，而降水、溫度、橡膠生長等因素具有很明顯的季節(jié)性變化，因此，CSIF較適合應(yīng)用于GPP的空間分布格局的估算應(yīng)用。

5種產(chǎn)品分別與海南島GPP進行影響分析，結(jié)果表明，CSIF和GOSIF比其余3種產(chǎn)品更能捕捉到GPP的變化，它們的決定系數(shù)空間分布特征最完整且最能反映橡膠林的GPP狀況。而SIFOCO2-1696和SIFLUE在海南島的可用數(shù)據(jù)缺失嚴重，在GPP估算中不如CISF、GOSIF和SIF005產(chǎn)品表現(xiàn)。

4 結(jié) 論

SIF遙感數(shù)據(jù)作為光合作用的“探測工具”，能有效估算生態(tài)系統(tǒng)GPP，因此，一系列SIF遙感產(chǎn)品應(yīng)運而生。這些產(chǎn)品來源不同，估算方式不同，個別產(chǎn)品在時空分辨率和有限時間范圍等方面在特定區(qū)域上的應(yīng)用尚不明確。本研究以海南島橡膠林為切入點，選取CSIF、GOSIF、SIFOCO2-1696、SIF005、SIFLUE5款產(chǎn)品，利用海南島橡膠林GPP為參考，評估不同SIF數(shù)據(jù)產(chǎn)品間的時空差異，及其對海南島橡膠林GPP估算的影響，結(jié)果表明：CSIF與GOSIF在時空變化的分析中具有更好的適用性，但CSIF不適用于長期趨勢分析，二者能夠良好地對海南島橡膠林生態(tài)系統(tǒng)GPP的時空分布特征進行指示；SIF005與GPP一致性分析中的表現(xiàn)僅次于前兩者；而SIFOCO2-1696在海南島近20 a內(nèi)的數(shù)據(jù)缺失較為嚴重，SIFLUE則是缺失部分海南島南部區(qū)域的數(shù)據(jù)，適用性較差。

致謝：海南島橡膠林GPP數(shù)據(jù)來源于Wei Cui等人的研究成果（https://doi.org/10.3390/f13060893）；CSIF數(shù)據(jù)（https://osf.io/8xqy6/）來源于Zhang Y等人研究成果；GOSIF數(shù)據(jù)（http://data.globalecology.unh.edu/data/GOSIF_v2/Monthly/）來源于Li和Xiao研究成果；SIF005（https://cornell.app.box.com/s/gkp4moy4grvqsus1q5oz7u5lc 30i7o41/folder/128864739508）由Wen J提供；SIFOCO2-1696（https://daac.ornl.gov/）由Yu等人提供；SIFLUE數(shù)據(jù)（http://jeodpp.jrc.ec.europa.eu/ftp/jrc-opendata/ECOCLIM/Downscaled-GOME2-SIF/v2.0/）由Gentine等人提供。特此感謝！