SST融合產(chǎn)品在渤黃海的質(zhì)量評(píng)估與修正研究

曹凱翔, 靳熙芳, 孫偉富, 王 進(jìn)

?

SST融合產(chǎn)品在渤黃海的質(zhì)量評(píng)估與修正研究

曹凱翔1, 靳熙芳2, 孫偉富1, 王 進(jìn)3

(1. 國(guó)家海洋局第一海洋研究所, 山東 青島 266061; 2. 國(guó)家海洋局北海預(yù)報(bào)中心, 山東 青島 266061; 3. 青島大學(xué), 山東 青島 266003)

為獲取高空間分辨率與空間覆蓋率的海表面溫度(Sea Surface Temperature, SST)產(chǎn)品, 基于最優(yōu)插值方法, 對(duì)微波輻射計(jì)WindSat、AMSR-E、ASMR2、HY-2 RM和紅外輻射計(jì)MODIS、AVHRR的SST觀測(cè)數(shù)據(jù)進(jìn)行融合, 生成了一種0.1°空間分辨率的每日SST融合產(chǎn)品, 利用浮標(biāo)數(shù)據(jù)在渤黃海區(qū)域進(jìn)行了精度評(píng)估和修正, 并分析了該區(qū)域的SST時(shí)空分布特征。結(jié)果表明: SST融合產(chǎn)品的在中國(guó)近海的精度為1.1℃, 利用浮標(biāo)數(shù)據(jù)修正后的精度略有改善; 利用修正的SST產(chǎn)品對(duì)渤黃海區(qū)域SST分布特征進(jìn)行了分析, 分析結(jié)果顯示, 渤黃海海域冬季海溫最為均勻, 春季在海水升溫過(guò)程中海溫不均勻性明顯。

海表溫度; 數(shù)據(jù)融合; 渤黃海; 質(zhì)量評(píng)估; 時(shí)空分布特征

海表面溫度(Sea Surface Temperature, SST)是海洋研究的基本參數(shù)之一, 其直接影響著大氣和海洋間的熱量、動(dòng)量和水汽交換, 是影響全球水循環(huán)的重要參數(shù), 影響著全球表面的能量收支平衡[1]。基于衛(wèi)星平臺(tái)的被動(dòng)遙感觀測(cè)技術(shù)是全球海表溫度監(jiān)測(cè)的主要手段。目前可進(jìn)行SST觀測(cè)的星載傳感器主要包括微波輻射計(jì)和紅外輻射計(jì)兩類。微波輻射計(jì)觀測(cè)可以穿透云霧, 受氣候影響小, 每天可覆蓋全球近90%的海域, 但是其空間分辨率較低, 為幾十千米。紅外輻射計(jì)工作頻段遠(yuǎn)高于微波輻射計(jì), 因此其空間分辨率可達(dá)千米量級(jí), 但紅外輻射計(jì)受天氣條件影響大, 不能穿透云霧, 單個(gè)紅外輻射計(jì)的每日有效觀測(cè)數(shù)據(jù)僅能覆蓋全球海洋的10%左右。采用一定的數(shù)據(jù)融合方法, 對(duì)星載微波輻射計(jì)和紅外輻射計(jì)SST觀測(cè)結(jié)果進(jìn)行融合, 可以有效地改善微波輻射計(jì)空間分辨率低、紅外輻射計(jì)空間覆蓋率差的問(wèn)題, 從而生成可覆蓋全球海洋的高分辨率SST多源融合產(chǎn)品。近年來(lái)相關(guān)研究受到各國(guó)研究機(jī)構(gòu)的普遍關(guān)注, 如美國(guó)國(guó)家海洋大氣管理局NOAA(National Oceanic and Atmospheric Administration)、NASA下屬噴氣推進(jìn)實(shí)驗(yàn)室JPL(Jet Propulsion Laboratory)、日本氣象廳JMA(Japan Meteorological Agency)、Met Office、遙感系統(tǒng)RSS(Remote Sensing System)均開(kāi)發(fā)了SST融合數(shù)據(jù)產(chǎn)品。這些產(chǎn)品空間分辨率為0.05°~0.25°, 時(shí)間分辨率一般為日均, 并普遍采用浮標(biāo)等實(shí)測(cè)數(shù)據(jù)進(jìn)行了精度驗(yàn)證[2]。

SST是目前海洋學(xué)界的研究熱點(diǎn), 它直接影響著大氣和海洋間的熱量、動(dòng)量和水汽交換過(guò)程, 對(duì)全球表面的能量收支起重要作用。本文基于最優(yōu)插值方法, 對(duì)微波輻射計(jì)WindSat、AMSR-E、ASMR2、HY-2 RM和紅外輻射計(jì)MODIS、AVHRR數(shù)據(jù)進(jìn)行融合, 制作了0.1°空間分辨率的SST融合產(chǎn)品, 并基于實(shí)測(cè)數(shù)據(jù)在中國(guó)渤海和黃海北部開(kāi)展了SST產(chǎn)品的精度驗(yàn)證和修正研究。這為開(kāi)展海氣相互作用的研究提供了重要的科學(xué)依據(jù), 具有良好的應(yīng)用前景。

1 SST融合產(chǎn)品介紹

1.1 SST多源遙感數(shù)據(jù)

SST融合產(chǎn)品所使用的遙感數(shù)據(jù)包括: WindSat、ASMR-E、AMSR2、HY-2 RM、MODIS(Aqua、Terra)和AVHRR等7個(gè)微波、紅外輻射計(jì)數(shù)據(jù)。

微波輻射計(jì)AMSR-E、AMSR2和WindSat的每日SST數(shù)據(jù)由遙感系統(tǒng)RSS提供, 其空間分辨率為25 km, 數(shù)據(jù)格式為二進(jìn)制, 數(shù)據(jù)產(chǎn)品中包括每日升軌/降軌的SST、風(fēng)速、水汽和云液水產(chǎn)品[3]。ASMR-E的在軌時(shí)間為2002年3月至2011年10月, 目前已停止工作。AMSR2是AMSR-E的后繼星, 于2012年5月發(fā)射入軌, 至今正常在軌運(yùn)行[4]。WindSat是歷史上第一個(gè)星載全極化微波輻射計(jì), 于2003年1月發(fā)射, 至今在軌正常運(yùn)行。本文提取了該3顆星的SST數(shù)據(jù), 剔除了SST默認(rèn)值和異常值。WindSat和AMSR-E的數(shù)據(jù)樣例如圖1所示。

圖1 WindSat(a)與AMSR-E(b)SST產(chǎn)品樣例

HY-2 RM微波輻射計(jì)數(shù)據(jù)由國(guó)家衛(wèi)星海洋應(yīng)用中心提供, 為L(zhǎng)2級(jí)沿軌數(shù)據(jù), 格式為HDF, 其中包含HY-2微波輻射計(jì)反演的SST、水汽含量、云液水含量等數(shù)據(jù)[5]。本文提取了HY-2的原始分辨率SST數(shù)據(jù)并進(jìn)行了數(shù)據(jù)有效性檢查。HY-2 SST數(shù)據(jù)樣例及掃描時(shí)間如圖2所示。

圖2 HY-2 SST網(wǎng)格數(shù)據(jù)

MODIS Aqua/Terra數(shù)據(jù)由NASA提供, 空間分辨率為4 km, 數(shù)據(jù)格式為HDF[6]。本文讀取了SST數(shù)據(jù)和相應(yīng)的質(zhì)量標(biāo)記位, 并進(jìn)行了質(zhì)量控制。MODIS/Terra于1999年12月入軌, MODIS/Aqua于2002年3月入軌, 目前兩星皆在軌正常運(yùn)行。MODIS SST數(shù)據(jù)樣例如圖3所示。

圖3 MODIS SST數(shù)據(jù)樣例

AVHRR Pathfinder每日SST數(shù)據(jù)由NOAA下屬國(guó)家海洋數(shù)據(jù)中心NODC(National Ocean Data Center)提供, 空間分辨率為4 km, 數(shù)據(jù)格式為NetCDF[6]。本文讀取了該星的SST數(shù)據(jù)和相應(yīng)的質(zhì)量標(biāo)記位, 并進(jìn)行了質(zhì)量控制。AVHRR首星于1981年6月發(fā)射入軌, 目前在軌的是其后繼星。AVHRR SST數(shù)據(jù)樣例如圖4所示。

圖4 AVHRR SST數(shù)據(jù)樣例

1.2 SST融合方法

圖5為融合SST產(chǎn)品的數(shù)據(jù)樣例。

圖5 SST融合個(gè)例

2 SST產(chǎn)品在渤黃海的質(zhì)量評(píng)估與修正

2.1 數(shù)據(jù)時(shí)空匹配

本文利用渤海和黃海北部(35°~41°N, 117°~ 127°E)6個(gè)浮標(biāo)的實(shí)測(cè)數(shù)據(jù), 對(duì)SST融合產(chǎn)品在該區(qū)域的精度進(jìn)行評(píng)估, 浮標(biāo)數(shù)據(jù)的時(shí)間范圍是2011—2014年。遙感-現(xiàn)場(chǎng)數(shù)據(jù)時(shí)空匹配過(guò)程中, 首先對(duì)浮標(biāo)數(shù)據(jù)進(jìn)行預(yù)處理, 剔除了浮標(biāo)觀測(cè)數(shù)據(jù)中的異常值; 由于浮標(biāo)對(duì)SST的觀測(cè)頻次很高, 因此對(duì)同一浮標(biāo)的觀測(cè)結(jié)果進(jìn)行了日平均; 進(jìn)而采用50 km、1 d的時(shí)空匹配窗口對(duì)浮標(biāo)和SST融合產(chǎn)品進(jìn)行了時(shí)空匹配, 共獲得匹配數(shù)據(jù)4 850組。

2.2 SST產(chǎn)品精度評(píng)估

基于以上時(shí)空匹配數(shù)據(jù), 本文計(jì)算了SST融合產(chǎn)品相對(duì)浮標(biāo)數(shù)據(jù)的平均偏差、均方根(root mean square, RMS)誤差和相關(guān)系數(shù), 作為SST產(chǎn)品精度評(píng)估的統(tǒng)計(jì)指標(biāo)。結(jié)果表明, SST融合產(chǎn)品的平均偏差為0.1℃, RMS誤差為1.1℃,相關(guān)系數(shù)為0.99。圖6是SST反演誤差的散點(diǎn)圖和誤差分布直方圖。

圖6 SST產(chǎn)品與浮標(biāo)數(shù)據(jù)的比較(a)和誤差分布直方圖(b)

一般而言, 星載微波輻射計(jì)和紅外輻射計(jì)在大洋區(qū)域的SST反演精度優(yōu)于1.0℃, HY-2微波輻射計(jì)的SST反演精度稍差, 約1.7℃[8]。本文生產(chǎn)的SST融合產(chǎn)品在中國(guó)近海的誤差略大于大洋區(qū)域, 其原因可能在于輻射計(jì)在近海區(qū)域的觀測(cè)數(shù)據(jù)容易受到陸地污染和無(wú)線電射頻干擾等因素的影響, 導(dǎo)致SST反演誤差增大[9]。

2.3 SST融合產(chǎn)品修正

本節(jié)利用SST融合產(chǎn)品和浮標(biāo)匹配數(shù)據(jù), 對(duì)SST融合產(chǎn)品進(jìn)行修正, 以改善其在中國(guó)近海的精度。首先將4 850組匹配數(shù)據(jù)隨機(jī)平均分為建模數(shù)據(jù)集和檢驗(yàn)數(shù)據(jù)集, 利用建模數(shù)據(jù)集建立SST融合產(chǎn)品與現(xiàn)場(chǎng)數(shù)據(jù)的經(jīng)驗(yàn)關(guān)系, 對(duì)融合產(chǎn)品進(jìn)行修正; 利用檢驗(yàn)數(shù)據(jù)集對(duì)修正后的SST產(chǎn)品進(jìn)行檢驗(yàn), 并與修正前的結(jié)果進(jìn)行比較。

由圖6可知, SST融合數(shù)據(jù)和浮標(biāo)觀測(cè)數(shù)據(jù)在散點(diǎn)圖中基本都分布在1︰1線附近, 因此本文選擇線性方程對(duì)SST融合數(shù)據(jù)和浮標(biāo)數(shù)據(jù)進(jìn)行擬合, 得到線性擬合方程為:

其中,Buoy為浮標(biāo)觀測(cè)數(shù)據(jù),Fusion為融合數(shù)據(jù)。

利用以上回歸方程對(duì)檢驗(yàn)數(shù)據(jù)集中的SST融合結(jié)果進(jìn)行了修正, 并與浮標(biāo)數(shù)據(jù)進(jìn)行比較。結(jié)果表明, 修正后的SST融合產(chǎn)品的平均偏差為0.01℃, 明顯優(yōu)于修正前的結(jié)果, RMS誤差為1.0℃, 略優(yōu)于修正前的結(jié)果。下圖為修正后的SST產(chǎn)品與浮標(biāo)實(shí)測(cè)數(shù)據(jù)的比較。

圖7 修正后的SST產(chǎn)品與浮標(biāo)的比較

3 渤黃海區(qū)域SST分布特征分析

本文利用以上制作的SST融合數(shù)據(jù), 以中國(guó)渤海和黃海北部(35°~41°N, 117°~127°E)為研究區(qū)域, 對(duì)該海域SST的時(shí)空分布特征進(jìn)行了分析, 計(jì)算了該海域2004—2014年各年4季的海表溫度平均值和標(biāo)準(zhǔn)差, 其中春季取3—5月, 夏季取6—8月, 秋季取9—11月, 冬季取12月—次年2月, 結(jié)果如圖8所示。

圖8 2004~2014年海表溫度平均值和標(biāo)準(zhǔn)差

圖中可見(jiàn), 渤黃海區(qū)域夏季海溫一般在22℃以上, 其中2011年夏季的平均海溫為20.9℃, 較其他年份同期海溫低1℃以上; 秋季海溫一般在18.5℃以上, 其中2011年秋季平均海溫18.1℃較往年略有偏低; 該海域春季平均海溫在2010—2011年也出現(xiàn)一個(gè)低值, 約7.1℃; 冬季海溫比較穩(wěn)定, 為3~3.5℃。海溫的標(biāo)準(zhǔn)差反映了海水溫度的均勻程度, 由圖8可知, 該海域春季海溫標(biāo)準(zhǔn)差最大, 為1.2℃以上, 表明海水在升溫過(guò)程中溫度梯度比較明顯; 冬季海溫標(biāo)準(zhǔn)差最小, 一般在0.6℃以下, 表明海溫此時(shí)最為均勻; 2011年夏季海溫的標(biāo)準(zhǔn)差明顯低于2010年和2012年, 表明上述2011年夏季平均海溫低于常年的現(xiàn)象是由于海溫整體偏低導(dǎo)致。

4 結(jié)論

本文利用國(guó)內(nèi)外多個(gè)微波、紅外輻射計(jì)的SST產(chǎn)品, 基于最優(yōu)插值方法制作了一種0.1°空間分辨率的全球每日SST產(chǎn)品; 并基于浮標(biāo)實(shí)測(cè)數(shù)據(jù)在我國(guó)渤黃海區(qū)域?qū)ST產(chǎn)品進(jìn)行了精度評(píng)估與修正, 最后對(duì)該海域11年的海溫分布特征進(jìn)行了分析, 得到結(jié)論如下: (1)相對(duì)于浮標(biāo)實(shí)測(cè)數(shù)據(jù), SST融合產(chǎn)品在渤黃海區(qū)域的偏差為0.1℃, RMS誤差為1.1℃; 由于輻射計(jì)在近海區(qū)域的數(shù)據(jù)容易受到陸地污染和無(wú)線電射頻干擾, 融合產(chǎn)品在近海的精度低于大洋區(qū)域。(2)本文利用浮標(biāo)實(shí)測(cè)SST數(shù)據(jù)對(duì)SST融合產(chǎn)品進(jìn)行了修正, 修正后SST融合產(chǎn)品的偏差為0.01℃, RMS誤差為1.0℃, 較修正前有所改善。(3)利用SST融合產(chǎn)品對(duì)渤黃海海域的SST分布特征進(jìn)行了分析, 結(jié)果表明, 渤黃海海域夏季水溫一般高于22℃, 冬季水溫一般為3~3.5℃; 該海域冬季海溫均勻, 標(biāo)準(zhǔn)差小于0.6℃; 春季海水升溫過(guò)程中海溫最不均勻, 海溫標(biāo)準(zhǔn)差大于1.2℃。

[1] Dash P, Ignatov A, Martin M, et al. Group for high resolution sea surface temperature (GHRSST) analysis field sinter-comparisons—Part 2: Near real time web-based level 4 SST Quality Monitor (L4-SQUAM)[J]. Deep-Sea Research II, 2012, 80: 31-43.

[2] Hosoda K, Qin Huiling. Algorithm for estimating sea surface temperatures based on aqua/MODIS global ocean data. 1. Developmentand validation of the algorithm [J]. Journal of Oceanography, 2011, 67 (1): 135-145.

[3] H?yer J L, Karagali I, Dybkj?r G, et al. Multi sensor validation and error characteristics of arctic satellite sea surface temperature observations [J]. Remote Sensing of Environment, 2012, 121: 335-346.

[4] Kurihara Y, Sakurai T, Kuragano T. Global daily sea surface temperature analysis using data from satellite microwave radiometer, satellite infrared radiometer and in-situ observations [J]. Weather Bull, 2006, 73: 1-18.

[5] Martin M, Dash P, Ignatov A, et al. Group for high resolution sea surface temperature (GHRSST) analysis field sinter-comparisons. Part 1: A GHRSST multi-product ensemble (GMPE) [J]. Deep-Sea Research II, 2012, 80: 21-30.

[6] Meissner T, Wentz F. Ocean Retrievals for WindSat[J]. IEEE Microrad, 2006, 7: 119-124.

[7] Reynolds R, Smith T. Improved global sea surface temperature analyses using optimum interpolation[J]. Journal of Climate, 1994, 7(6): 929-948.

[8] Reynolds W, Smith T, Liu C, et al. Daily high-resolu-tion-blended analyses for sea surface temperature[J]. Journal of Climate, 2007, 20(11): 5473- 5496.

[9] Shaw A, Vennell R. A front-following algorithm for AVHRR SST imagery[J]. Remote Sensing of Environment, 2000, 72(3): 317-327.

[10] Zabolotskikh E, Mitnik L, Reul N, et al. New possibilities for geophysical parameter retrieval opened by GCOM-W1 AMSR2[J]. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing, 2014, 8(9): 4248-4261.

[11] Zhao Yili, Zhu Jianhua, Lin Mingsen, et al. Assessment of the initial sea surface temperature product of the scanning microwave radiometer aboard on HY-2 Satellite[J]. Acta Oceanologica Sinica, 2014, 33(1): 109- 113.

[12] 蔣興偉, 林明森, 宋清濤. 海洋二號(hào)衛(wèi)星主被動(dòng)微波遙感探測(cè)技術(shù)研究[J]. 中國(guó)工程科學(xué), 2013, 15(7): 4-11.Jiang Xingwei, Lin Mingsen, Song Qingtao. Active and passive microwave remote sensing technology of the HY-2A ocean satellite mission[J]. Engineering Sciences, 2013, 15(7): 4-11.

[13] 孫廣輪, 關(guān)道明, 趙冬至, 等. 星載微波遙感觀測(cè)海表溫度的研究進(jìn)展[J]. 遙感技術(shù)與應(yīng)用, 2013, 28(4): 721-730. Sun Guanglun, Guan Daoming, Zhao Dongzhi, et al. Research on observing sea surface temperature(SST) based on microwave remote sensing by satellite[J]. Remote Sensing Technology and Application, 2013, 28(4): 721-730.

(本文編輯: 劉珊珊)

Quality assessment and correction of SST fusion product in the Bohai Sea and the Huanghai Seas

CAO Kai-xiang1, JIN Xi-fang2, SUN Wei-fu1, WANG Jin3

(1. The First Institute of Oceanography, State Oceanic Administration, Qingdao 266061, China; 2. North China Sea Marine Forecasting Center of State Oceanic Administration, Qingdao 266061, China; 3. Qingdao University, Qingdao 266003, China)

Sea surface temperature (SST) is a key parameter used in physical ocean research and has an important influence on sea-air interaction and global climate change. Space-borne microwave and infrared radiometers are the main methods used in SST remote sensing. However, although microwave radiometers provide all-weather technology their spatial resolution is sparse, and although infrared radiometers have a good resolution they are affected by cloud. In this study, a 0.1° daily SST fusion product is analyzed using the optimum interpolation method. Data resources include WindSat, AMSR-E, ASMR2, HY-2 RM, MODIS, and AVHRR. Buoy data are used to validate and correct the fusion SST product and its SST characteristics in the Bohai and the Huanghai Seas are also studied. Results show that the RMS error of this SST product is 1.1℃, but this decreases slightly after correction. SST characteristics show the lowest standard deviation of SST is in winter, and SST tends to be inhomogeneous in spring when sea temperatures rise.

SST(sea surface temperature); data fusion; the Bohai Sea and the Huanghai Sea; quality assessment; distribution characteristics

Aug. 31, 2016

曹凱翔(1994-), 男, 山東諸城人, 碩士, 研究方向?yàn)槲⒉ㄝ椛溆?jì)信息提取, 電話: 0532-58750699, E-mail: 676751961@ qq.com; 王進(jìn),通信作者, 博士, 主要從事被動(dòng)微波海洋遙感方向研究, 電話: 0532-85950680, E-mail: wangjin@qdu.edu.cn

P731.11

A

1000-3096(2017)09-0050-06

10.11759/hykx20160831001

2016-08-31;

2016-09-20

國(guó)家重點(diǎn)研發(fā)計(jì)劃項(xiàng)目(2016YFA0600102); 國(guó)家高技術(shù)研究發(fā)展計(jì)劃(863計(jì)劃)(2013AA122803); 海洋公益性行業(yè)科研專項(xiàng)(2013418032)

[The National Key Research and Development Program of China, No. 2016YFA0600102; National High-tech R&D Program of China (863 Program), No. 2013AA122803; Public Science and Technology Research Projects of Ocean, No. 2013418032]