許玉壯,何賢強(qiáng),白雁,朱乾坤,龔芳
1.上海交通大學(xué) 海洋學(xué)院, 上海 200240;
2.自然資源部第二海洋研究所 衛(wèi)星環(huán)境海洋環(huán)境動力學(xué)國家重點實驗室, 杭州 310012
中國于2018年9月發(fā)射了第3顆海洋水色衛(wèi)星HY-1C,其上搭載海洋水色水溫掃描儀COCTS(Chinese Ocean Color and Temperature Scanner)。COCTS 具有8 個可見光(VIS)與近紅外(NIR)波段以及2 個熱紅外波段,星下點空間分辨率為1.1 km(沈亞峰 等,2020;Liu 等,2021)。COCTS實現(xiàn)每天全球覆蓋,探測海洋水色要素(如葉綠素、懸浮泥沙和可溶有機(jī)物等)和海表溫度等,為海洋資源調(diào)查、生態(tài)環(huán)境監(jiān)測、應(yīng)急減災(zāi)等提供觀測數(shù)據(jù)支持(沈亞峰 等,2020;范文龍 等,2022)。對水色遙感產(chǎn)品的真實性檢驗是衛(wèi)星產(chǎn)品定量化應(yīng)用的前提。Bailey 和 Werdell(2006)利用全球?qū)崪y數(shù)據(jù)確認(rèn)了SeaWiFS水色遙感產(chǎn)品在全球大部分開闊大洋水體不確定性均小于5%,推動了海洋水色遙感的發(fā)展;Hlaing等(2013)利用兩個AERONET-OC(WaveCIS 和LISCO)站位評價了VIIRS業(yè)務(wù)化產(chǎn)品在近岸水體的表現(xiàn),在兩個站點的不確定性分別為23.4% 和26.3%,證明了VIIRS 在近岸海洋水色的觀測能力;Liu 等(2022)利用全球AERONET-OC 實測數(shù)據(jù)對采用POLYMER、C2RCC、BAC、SeaDAS 等7 種不同大氣校正算法的OLCI 水色遙感產(chǎn)品進(jìn)行了評價,結(jié)果表明POLYMER算法表現(xiàn)最好(412 nm、443 nm、490 nm、560 nm 和665 nm 的不確定性分別為46.23%、34.89%、15.79%、11.78%和21.63%),為OLCI 遙感器在水環(huán)境監(jiān)測應(yīng)用中選擇合適的大氣校正算法提供了參考。SeaWiFS、MODIS、VIIRS、OLCI 等水色衛(wèi)星遙感器的產(chǎn)品真實性檢驗工作(Hlaing 等,2013;Ladner 等,2014;Wang 等,2014;Zibordi 等,2018;Barnes 等,2019;Kyryliuk和Kratzer,2019)極大地推動了各自遙感產(chǎn)品的廣泛應(yīng)用(Wang 等,2018;Balasubramanian等,2020;Hernández 等,2020)。因此,COCTS遙感產(chǎn)品真實性檢驗對發(fā)揮其應(yīng)用價值也十分重要。
在HY-1C 衛(wèi)星COCTS 定標(biāo)檢驗方面,Song 等(2019)利用MODIS 氣溶膠遙感產(chǎn)品和MOBY 實測光譜數(shù)據(jù)建立了COCTS 的替代定標(biāo)系數(shù)估計模型,并在短期數(shù)據(jù)驗證中取得了較好的結(jié)果;Chen 等(2021)從多星水色產(chǎn)品數(shù)據(jù)融合的角度,利用SeaDAS 7.5.3 軟件對中國海洋衛(wèi)星數(shù)據(jù)服務(wù)系統(tǒng)發(fā)布HY1C/COCTS 的Level 1 級產(chǎn)品進(jìn)行大氣校正獲得遙感反射率(Rrs)產(chǎn)品,研究了COCTS 的信噪比和水體Rrs 產(chǎn)品在南太平洋環(huán)流區(qū)域(穩(wěn)定清潔的一類水體)的不確定性,發(fā)現(xiàn)在開闊大洋水體,HY-1C/COCTS 的藍(lán)光波段Rrs 不確定性小于5%。Rrs 是進(jìn)一步反演葉綠素、懸浮泥沙等水色要素產(chǎn)品的基礎(chǔ)產(chǎn)品,其精度直接決定了水色要素產(chǎn)品的精確性(Xu 等,2021;仲蘇珂 等,2022);但是其在近海水體精度有待進(jìn)一步認(rèn)識。隨著COCTS 在軌運行積累了大量觀測數(shù)據(jù),而目前對COCTS 近海業(yè)務(wù)化產(chǎn)品精度表現(xiàn)研究工作較少,因此對COCTS全球近海業(yè)務(wù)化Rrs產(chǎn)品進(jìn)行真實性檢驗具有必要性。
本文利用AERONET-OC 站點實測水體光譜數(shù)據(jù),對COCTS的近海業(yè)務(wù)化Rrs產(chǎn)品進(jìn)行系統(tǒng)性評估。AERONET-OC 站點分布在全球近海不同水域,覆蓋了不同類型的水體,可較全面地評估COCTS的Rrs產(chǎn)品精度。
本文使用從2018年9至12月和2019年11月至2020年12月內(nèi)與AERONET-OC 實測資料匹配的HY-1C/COCTS Level-2A 的各波段Rrs 數(shù)據(jù)產(chǎn)品,直接下載自中國海洋衛(wèi)星數(shù)據(jù)服務(wù)系統(tǒng)(https://osdds.nsoas.org.cn[2021-01-25]),系統(tǒng)數(shù)據(jù)文件版本號10。因進(jìn)行本文工作時中國海洋衛(wèi)星數(shù)據(jù)服務(wù)系統(tǒng)正在重新處理2019年1—10月的COCTS 數(shù)據(jù),故本文剔除了該段時間數(shù)據(jù),共獲得205 景有效影像。COCTS 波段信息如表1 所示(表格信息國家衛(wèi)星海洋應(yīng)用中心官方網(wǎng)站,http://www.nsoas.org.cn/news/content/2018-11/23/44_696.html[2022-01-11]。
表1 COCTS可見光及近紅外波段Table 1 Technical parameters for COCTS at VIS and NIR band
AERONET-OC 站點使用改裝的CE-318 太陽光度計來自動測量水面上行輻亮度(LT)、天空光下行輻亮度(Li)和向下輻照度(ES)。最新的太陽光度系列(CE318-T)可以提供12 個波段的數(shù)據(jù),其離水輻亮度LW可由下式得到:式中,ρ為海面反射率,λ、θ(θ')、?、θ0、W分別為波長、觀測天頂角、相對方位角、太陽天頂角和風(fēng)速。在對LW進(jìn)行觀測幾何方位校正和二向性校正后,進(jìn)一步將其轉(zhuǎn)換為歸一化離水輻亮度Lwn(λ):
式中,D為日地距離,td為大氣漫射透過率,可由向下輻照度得到??捎上率綄wn(λ)轉(zhuǎn)換為本文中使用遙感反射率Rrs(λ):
式中,F(xiàn)0為平均日地距離的大氣層頂太陽輻照度(Zibordi 等,2009;Liu 等,2022)。本文中使用的F0值來自Thuillier等(2003)。
AERONET-OC 站點從當(dāng)?shù)貢r間8:00 到16:00每30 min 進(jìn)行一次光譜測量,在陰雨天氣自動停止觀測(Zibordi 等,2009)。AERONET-OC 通過在所有站點使用統(tǒng)一的測量系統(tǒng)和數(shù)據(jù)處理方法,保證了實測數(shù)據(jù)的質(zhì)量,其在668 nm 的不確定性約為7.8%,在短波長的不確定性約為5%(Zibordi等,2009和2010)。在2018年9月到2020年12月期間,與COCTS 有匹配樣本的AERONET-OC 站點共有12個,分布在近岸和內(nèi)陸湖(Zibordi 等,2010)。本文從AERONET官方網(wǎng)站(https://aeronet.gsfc.nasa.gov/[2021-01-25])獲取了2018年9月到2020年12月期間的12個站點數(shù)據(jù)(表2)。
表2 本文使用的AERONET-OC站點信息Table 2 The AERONET-OC sites information used in this paper
本文將實測Rrs 光譜歸一化,獲得歸一化遙感反射率(nRrs)譜線,
式中,N為COCTS 可見光范圍的特征波段數(shù),本文取值為6(對應(yīng)COCTS 的6 個可見光波段),λi分別對應(yīng)412 nm、 443 nm、 490 nm、 520 nm、565 nm 和670 nm(Wei 等,2016)。nRrs光譜和原始Rrs光譜形狀一致、取值范圍為0—1?;跉w一化的AERONET-OC 實測Rrs 光譜,本文結(jié)合Wei等(2016)提出的水體光學(xué)分類方法和中國近海水體光學(xué)特性調(diào)查的光譜分類方法(李銅基,2012),將AERONET-OC 實測光譜數(shù)據(jù)分為A(清潔)、B(比較清潔)、C(輕微渾濁)、D(比較渾濁)共4種。圖1(a)和圖1(b)分別給出了4 種水體類型的nRrs譜線和平均nRrs特征譜線,4種水體的基本特點(李銅基,2012)如表3。光譜數(shù)據(jù)質(zhì)量評價通過將目標(biāo)Rrs 光譜與特征光譜進(jìn)行比較,對目標(biāo)光譜從0 到1 進(jìn)行質(zhì)量評分,分值越高,光譜質(zhì)量越高(Wei 等,2016),本文剔除了質(zhì)量評分小于1的實測光譜數(shù)據(jù)。本文利用分類和質(zhì)量控制后AERONET-OC 光譜數(shù)據(jù)對COCTS 的Rrs 產(chǎn)品在不同光學(xué)水體類型的精度進(jìn)行評價,同時也對各波段Rrs 產(chǎn)品在全部類型水體的精度進(jìn)行評價。
表3 4種水體類型的Rrs光譜基本特點Table 3 The basic characteristics of water and Rrs spectra of 4 water types
圖1 4種水體類型的歸一化Rrs光譜和平均歸一化Rrs光譜Fig.1 Normalized Rrs spectra and Mean normalized Rrs spectra characteristics of four water types
在波段匹配上,AERONET-OC 站點只提供幾個特定波長的數(shù)據(jù),且不同站點提供的波段組合不完全相同,本文中使用AERONET-OC 站點提供數(shù)據(jù)波段信息如表4 所示。因此,AERONET-OC數(shù)據(jù)的波段和COCTS 波段并不完全匹配,特別是COCTS 的520 nm 和565 nm 波段AERONET-OC 未提供實測數(shù)據(jù)。因此本文對AERONET-OC 數(shù)據(jù)根據(jù)COCTS 的波段進(jìn)行了波長插值處理。AERONET-OC 站點未提供750 nm 波段數(shù)據(jù),因此本文僅對COCTS可見光波段的業(yè)務(wù)化Rrs產(chǎn)品進(jìn)行了評價。由于AERONET-OC 波段和COCTS 存在差異,可能會導(dǎo)致評估誤差存在高估。
表4 AERONET-OC數(shù)據(jù)波段信息Table 4 AERONET-OC site data band information
空間匹配上,本文選取了以AERONET-OC 站點為中心的COCTS影像的5 × 5 像元空間窗口,并進(jìn)行質(zhì)量控制:(1)舍棄無效像元數(shù)大于像元空間窗口總像元數(shù)一半的數(shù)據(jù);(2)使用空間變異系數(shù)(空間變異系數(shù)=標(biāo)準(zhǔn)差/均值)來衡量像元空間窗口內(nèi)的空間異質(zhì)性,并剔除了空間變異系數(shù)大于0.4 的數(shù)據(jù)(Harding 等,2005;Barnes 等,2019)。在同一景衛(wèi)星影像中,COCTS 數(shù)據(jù)存在只有部分波段為有效數(shù)據(jù)的情況,并且565 nm 波段數(shù)據(jù)失效率顯著高于其他波段,為保證有足夠的統(tǒng)計數(shù)據(jù),本文保留了565 nm 無效但其他波段均有效的Rrs匹配數(shù)據(jù)。
在時間匹配上,為減少水體和大氣隨時間快速變化的影響,要求COCTS 與AERONET-OC 觀測時間窗口差異小于1 h。
本文利用平均相對偏差PD(average Percent Difference)、平均絕對百分比誤差A(yù)PD(average Absolute Percent Difference)作為精度評價指標(biāo)(Hlaing 等,2013)。PD 為衛(wèi)星數(shù)據(jù)和實測數(shù)據(jù)的偏差,正值代表衛(wèi)星Rrs數(shù)據(jù)相較于實測Rrs高估,負(fù)值代表低估;APD 用來衡量對比量之間的不確定性。PD 值為匹配的每一組數(shù)據(jù)的PDi均值,計算公式為
相似的,APD 的值為匹配的每一組數(shù)據(jù)的APDi的均值:
式中,i為一組匹配數(shù)據(jù)的計數(shù)下標(biāo),xi和yi分別代表匹配的每一組數(shù)據(jù)中的實測數(shù)據(jù)和衛(wèi)星反演數(shù)據(jù)(Hlaing 等,2013)。
COCTS 和AERONET-OC 匹配的Rrs 光譜曲線對比見圖2,A、B、C、D這4種水體類型中分別獲得了28對、40對、48對和89對匹配樣本。在A 類水體中,COCTS 和AERONET-OC 光譜曲線隨波長增加變化趨勢一致,但COCTS 值高于AERONETOC,兩者差異隨波長呈現(xiàn)先增加后減小的變化趨勢,在520 nm處差異最大;在B類水體中,COCTS的Rrs 光譜均值在< 565 nm 的波段高于實測光譜均值曲線,且COCTS 光譜標(biāo)準(zhǔn)差大于實測光譜,總體上衛(wèi)星光譜均值曲線和實測光譜均值曲線變化趨勢一致;在C類水體中,COCTS光譜均值曲線和實測光譜均值曲線具有一致的變化趨勢,且兩者的均值接近,在412 nm 波段,COCTS 的標(biāo)準(zhǔn)差大于實測數(shù)據(jù)標(biāo)準(zhǔn)差;在D類光學(xué)中,COCTS光譜和實測光譜具有相似的變化趨勢,COCTS 光譜均值曲線分布顯著低于實測光譜,在412 nm 波段,存在均值減去標(biāo)準(zhǔn)差為負(fù)值的情況,這表明COCTS 光譜在412 nm波段存在負(fù)值。
圖2 每類水體匹配的光譜對比(實線為均值,顏色填充范圍為均值±標(biāo)準(zhǔn)差)Fig.2 Spectral comparison for each water type(solid lines indicate mean and the coloring range indicate mean±standard deviation)
本文對4 種光學(xué)水體所有匹配的COCTS 和AERONET-OC 的Rrs 光譜進(jìn)行了散點對比分析(圖3)。在A、B、C、D 這4 種水體散點數(shù)分別為148 個、255 個、265 個和495 個,衛(wèi)星Rrs 數(shù)據(jù)和實測Rrs 數(shù)據(jù)相關(guān)系數(shù)(R)分別為0.879、0.680、0.780 和0.690。在A 類水體(PD=6.79%,APD=38.79%)中,散點緊密不均勻的分布在1∶1 線兩側(cè),1∶1 線上側(cè)散點多,下側(cè)散點少;在B 類水體(PD=18.37%,APD=44.44%)中,1∶1 線上側(cè)存在較多散點偏離1∶1 線較遠(yuǎn),其余散點相對均勻緊密分布在1∶1 線兩側(cè);在A 類水體和B 類水體中均呈現(xiàn)COCTS 數(shù)據(jù)呈現(xiàn)輕微高估的現(xiàn)象;在C 類水體(PD=-14.38%,APD=40.85%)中COCTS數(shù)據(jù)呈現(xiàn)輕微低估,多數(shù)散點相對均勻分布在1∶1 線兩側(cè),部分散點的縱軸為負(fù)值,在1∶1 線下側(cè)有少量散點偏離較遠(yuǎn);在D 類水體(PD=-32.35%,APD=47.14%)中,散點分布相對分散,存在較多散點的縱軸為負(fù)值的現(xiàn)象,且多數(shù)散點分布在1∶1線下側(cè),COCTS數(shù)據(jù)顯著低估。
圖3 4種水體類別中衛(wèi)星與實測數(shù)據(jù)遙感反射率散點圖Fig.3 Comparison of COCTS and AERONET-OC Rrs in each water type
本文對4 種水體類型中COCTS 和AERONETOC 匹配的Rrs 光譜進(jìn)行了單波段散點對比分析(圖4),A、B、C、D 這4 種水體類型在565 nm 波段匹配散點數(shù)分別為8 個、25 個、25 個和50 個,在其他波段匹配的散點數(shù)分別為28 個、40 個、48 個和89 個。在412 nm(PD=11.21%、APD=23.13%)和443 nm(PD=26.38%、APD=27.32%)波段,A 類水體衛(wèi)星數(shù)據(jù)和實測數(shù)據(jù)一致性均較好。在490 nm 波段,A、B、C 這3 類水體的APD值均小于30%;520 nm 波段,在C 類水體(PD=15.74%、APD=25.37%)、D 類水體(PD=1.73%、APD=29.76%)衛(wèi)星數(shù)據(jù)和實測數(shù)據(jù)的一致性相對較好,且比在A、B 類水體(APD 值均大于30%)的表現(xiàn)更好;565 nm 波段,在A 類水體(PD=32.36%、APD=30.68%)衛(wèi)星數(shù)據(jù)表現(xiàn)為高估,在B、C、D這3種水體類型中APD值均小于23%,衛(wèi)星數(shù)據(jù)和實測數(shù)據(jù)一致性較好;在670 nm 波段,4 種水體類型APD 值均在45%—75%,PD 值均小于-45%,衛(wèi)星數(shù)據(jù)在4 種水體類型中均表現(xiàn)為顯著低估。COCTS 不同波段Rrs 產(chǎn)品在4 類水體的精度亦存在差異。
圖4 4種水體類別中衛(wèi)星與實測數(shù)據(jù)單波段遙感反射率散點圖Fig.4 Comparison of COCTS and in-situ Rrs in each band of each water type
本文對匹配的AERONET-OC 數(shù)據(jù)和COCTS 數(shù)據(jù)分波段進(jìn)行了時間序列對比分析(圖5)。在整體上,衛(wèi)星數(shù)據(jù)和實測數(shù)據(jù)的散點交錯分布、數(shù)據(jù)值相近,并且隨時間序列具有相似的變化趨勢,COCTS 數(shù)據(jù)和實測數(shù)據(jù)具有較好一致性。但在2020年9月15日之后,COCTS 數(shù)據(jù)相較于AERONET-OC 數(shù)據(jù)始終表現(xiàn)為低估且COCTS 數(shù)據(jù)在412 nm 和443 nm 多為負(fù)值。經(jīng)檢查原始數(shù)據(jù)發(fā)現(xiàn)A、B、C這3種水體類型在該時間段內(nèi)無有效匹配數(shù)據(jù),該段時間全部為D 類水體(圖6)。在剔除了該段數(shù)據(jù)后,D 類水體各波段散點圖如圖7,與圖4中D類水體統(tǒng)計數(shù)據(jù)相比,各波段APD值減小均不超過5%,衛(wèi)星數(shù)據(jù)和實測數(shù)據(jù)間R 值和PD 值變化也不大;同時,在D 類水體獲得有效數(shù)據(jù)的其他時間段中,在412 nm 和443 nm 波段的COCTS 數(shù)據(jù)也有較多負(fù)值出現(xiàn),在490 nm、520 nm、565 nm、670 nm 波段COCTS 的數(shù)據(jù)低估現(xiàn)象也普遍存在。因此,本文認(rèn)為2020年9月15日之后時間段的數(shù)據(jù)同D 類水體的其它時間段相比無明顯異常。
圖5 不區(qū)分水體類別的衛(wèi)星和實測數(shù)據(jù)各波段Rrs時間序列對比圖Fig.5 Comparison of Rrs time series of COCTS and AERONET-OC for each band without distinguishing water types
圖6 D類水體的衛(wèi)星和實測數(shù)據(jù)各波段Rrs時間序列對比圖Fig.6 Comparison of Rrs time series of COCTS and AERONET-OC for each band in class D water type
圖7 衛(wèi)星與實測數(shù)據(jù)單波段遙感反射率散點圖(D類光學(xué)水體)Fig.7 Comparison of single-band Rrs between COCTS and AERONET-OC of optical water in Class D
Chen 等(2021)研究發(fā)現(xiàn)在水體極為清潔的南太平洋海盆區(qū)域,COCTS 的Rrs 產(chǎn)品在藍(lán)光波段的不確定性小于5%,衛(wèi)星Rrs 數(shù)據(jù)和實測數(shù)據(jù)具有較好的一致性。而本文所用的AERONETOC 站點主要分布在近海區(qū)域,因近海的水體懸浮泥沙濃度相對較高,導(dǎo)致近紅外波段的離水輻亮度顯著大于0,從而導(dǎo)致大氣校正過程中氣溶膠散射貢獻(xiàn)被高估,使得反演Rrs 數(shù)據(jù)出現(xiàn)低估問題。Zibordi等(2006)和Hlaing等(2013)利用AERONET-OC 實測數(shù)據(jù)分別對MODIS、VIIRS 的歸一化離水輻亮度(Lwn)產(chǎn)品進(jìn)行了評價,發(fā)現(xiàn)在渾濁水體,兩顆遙感器在412 nm、443 nm 附近波段均出現(xiàn)負(fù)值,存在低估,在670 nm 附近波段表現(xiàn)為低估,在490—565 nm 的波段,兩顆遙感器的反演數(shù)據(jù)和實測數(shù)據(jù)具有較好的一致性,這與本文對COCTS的檢驗結(jié)果表現(xiàn)一致。
本文利用全球12個AERONET-OC 站點實測數(shù)據(jù)對HY-1C/COCTS 業(yè)務(wù)化Rrs 產(chǎn)品進(jìn)行了精度評價。本文將實測光譜水體分為A(清潔)、B(比較清潔)、C(輕微渾濁)、D(比較渾濁)共4 種水體類型。在4 種水體類型中,A、C 類水體R值分別為0.879 和0.780,B、D 類水體R值分別為0.680 和0.690,衛(wèi)星Rrs 數(shù)據(jù)和實測Rrs 數(shù)據(jù)相關(guān)性較好。在A 類水體衛(wèi)星和實測數(shù)據(jù)PD 值為6.79%,APD值38.79%,衛(wèi)星數(shù)據(jù)相較于實測數(shù)據(jù)輕微高估。在B類水體衛(wèi)星Rrs數(shù)據(jù)和實測數(shù)據(jù)PD值為18.73%,APD 值小于45%,衛(wèi)星數(shù)據(jù)略微高估。在C 類水體,412 nm、443 nm 波段出現(xiàn)遙感Rrs 數(shù)據(jù)為負(fù)值的問題,PD 值為-14.38%,APD 值為40.85%,衛(wèi)星數(shù)據(jù)存在輕微低估。在D 類水體,412 nm、443 nm 波段存在較多Rrs 為負(fù)值的散點,衛(wèi)星Rrs數(shù)據(jù)和實測數(shù)據(jù)PD 值為-32.35%,APD 值47.14%,衛(wèi)星數(shù)據(jù)顯著低估。
COCTS各個波段在4種水體類型的表現(xiàn)也存在差異。在412 nm 和443 nm 波段,A 類水體中衛(wèi)星Rrs 數(shù)據(jù)和實測數(shù)據(jù)一致性較好,而在B、C、D 這3 類水體衛(wèi)星數(shù)據(jù)和實測數(shù)據(jù)一致性相對較差。在490 nm波段,A、B、C這3種水體類型中衛(wèi)星和實測數(shù)據(jù)APD 均小于30%,一致性相對較好;在520 nm 波段,COCTS 在C、D 類水體中衛(wèi)星數(shù)據(jù)和實測數(shù)據(jù)一致性相對較好,且一致性表現(xiàn)表現(xiàn)好于在A、B 類水體;在560 nm 波段,除A 類水體外,在其他3類水體APD值均小于23%,衛(wèi)星數(shù)據(jù)和實測數(shù)據(jù)一致性較好;然而,在670 nm 波段四種水體類型的PD 值均小于-45%、APD 均大于45%,衛(wèi)星Rrs 數(shù)據(jù)和實測一致性較差,遙感Rrs 數(shù)據(jù)嚴(yán)重低估。由于AERONET-OC波段和COCTS存在差異,可能會導(dǎo)致評估誤差存在高估,未來需要使用連續(xù)光譜的實測數(shù)據(jù)進(jìn)一步評價COCTS產(chǎn)品精度。
志 謝感謝國家衛(wèi)星海洋應(yīng)用中心提供HY-1C/COCTS Level-2 數(shù)據(jù)產(chǎn)品,感謝Joji Ishizaka,Giuseppe Zibordi,Marco Talone,Young-Je Park,Burton Jones,Matthew Ragan,Bill Gibson,Sherwin Ladner 提供的ARIAKE_TOWER,Casablanca Platform,Galata Platform,Helsinki Lighthouse,Irbe Lighthouse,Gloria,Gustav Dalen Tower,Section-7 Platform,Venise,Socheongcho,USC SEAPRISM,WaveCIS Site CSI 6 站點的AERONET-OC 實測水體光譜數(shù)據(jù)。感謝衛(wèi)星海洋環(huán)境動力學(xué)國家重點實驗室衛(wèi)星地面站和海洋遙感數(shù)據(jù)在線分析平臺(SatCO2)提供數(shù)據(jù)收集及處理幫助。