穆珊珊 李海艷 吳明柏, 2

基于Sentinel-1A的全球有效波高的反演研究*

穆珊珊1李海艷1①吳明柏1, 2

(1. 中國(guó)科學(xué)院大學(xué) 北京 100049; 2. 中國(guó)科學(xué)院地理科學(xué)與資源研究所 北京 100101)

本文利用神經(jīng)網(wǎng)絡(luò)的技術(shù)手段, 針對(duì)Sentinel-1A二級(jí)波模式數(shù)據(jù)提出一種用于海浪有效波高(Hs)反演的模型——N_N模型。該模型在基于ERS2 SAR波模數(shù)據(jù)開(kāi)發(fā)的雙參數(shù)模型的基礎(chǔ)上, 加入經(jīng)度、緯度、方位向截?cái)嗖ㄩL(zhǎng)(c)、圖像偏斜(skewness, skew)、圖像峰度(kurtosis, kurt)、衛(wèi)星平臺(tái)距目標(biāo)物的距離與衛(wèi)星飛行速度之比()等其他參數(shù)信息, 根據(jù)不同輸入?yún)?shù)的組合, 建立了14個(gè)模型用于Hs反演, 旨在分析各參數(shù)對(duì)有效波高反演的影響。通過(guò)分析表明, 14個(gè)N_N模型相關(guān)系數(shù)都在0.8以上。隨著c、參數(shù)的加入, N_N模型性能均大幅上升, 且c參數(shù)對(duì)模型性能的改善作用更加明顯, 相關(guān)系數(shù)提升0.06左右, 均方根誤差(Root Mean Squared Error, RMSE)下降0.12m左右。另外, skew與kurt的加入也使N_N模型性能有所改善, RMSE下降0.03m左右, 相關(guān)系數(shù)提升0.01左右。其中, N_N10模型效果最佳且性能最穩(wěn)定, 與歐洲中程天氣預(yù)測(cè)中心(the European Centre for Medium-Range Weather Forecasts, ECMWF)數(shù)據(jù)對(duì)比, 相關(guān)系數(shù)(CORR)達(dá)到0.905, 散射指數(shù)(Scattering Index, SI)與RMSE最低, 分別為18.74%、0.502m, 與獨(dú)立測(cè)量的浮標(biāo)數(shù)據(jù)的相關(guān)系數(shù)達(dá)到了0.894。

神經(jīng)網(wǎng)絡(luò); 有效波高; 方位向截?cái)嗖ㄩL(zhǎng); 歸一化雷達(dá)后向散射系數(shù)

海浪是海洋最明顯的表面特征, 波長(zhǎng)范圍從幾厘米到數(shù)百米, 波高范圍從海洋表面的微小擾動(dòng)到數(shù)十米, 對(duì)海洋工程、船舶設(shè)計(jì)、海上運(yùn)輸以及海洋污染的消散等都有很大的影響。因此, 通過(guò)有效的技術(shù)手段了解海浪的統(tǒng)計(jì)特性, 如有效波高、平均波周期等尤為重要。

自1978年發(fā)射Seasat1搭載合成孔徑雷達(dá)(Synthetic Aperture Radar, SAR)開(kāi)始, ERS-1/2、RADARSAT-1和ENVISAT、Sentinel-1持續(xù)運(yùn)行, 星載SAR不斷提供全球范圍內(nèi)實(shí)時(shí)動(dòng)態(tài)的海面波浪信息, 憑借其間接的、大范圍的測(cè)量方式和全天候、全天時(shí)、高分辨率的測(cè)量能力, 成為海浪信息獲取的重要途徑之一(Jackson, 2004)。利用SAR進(jìn)行海浪統(tǒng)計(jì)參數(shù)的反演有兩種常見(jiàn)的方式, 一種是通過(guò)SAR圖像譜獲取二維海浪方向波譜, 繼而獲取海浪統(tǒng)計(jì)參數(shù)。例如, 有效波高(Hs)可以表示為方向波譜的積分:

第一種通過(guò)SAR圖像譜獲取海浪譜, 反演海浪參數(shù)的方式在海浪遙感應(yīng)用中是很常見(jiàn)的。但海洋的SAR圖像常常在方位方向顯得十分模糊且導(dǎo)出的波譜容易發(fā)生畸變(Grieco, 2016; Stopa, 2017)。這是由于SAR特殊的成像機(jī)制——速度聚束調(diào)制引起的。受這一調(diào)制作用的影響, 海浪的隨機(jī)運(yùn)動(dòng)會(huì)引起SAR后向散射回波信號(hào)的多普勒頻移, 造成SAR不能對(duì)低于某一波長(zhǎng)(即方位向截?cái)嗖ㄩL(zhǎng))的波浪進(jìn)行成像。在這種情況下, 海浪的成像過(guò)程中常表現(xiàn)出很強(qiáng)的非線性, 對(duì)SAR圖像影響很大(Hasselmann, 1985; Jackson, 2004)。因此, 通過(guò)SAR圖像推導(dǎo)二維海浪譜絕非一項(xiàng)簡(jiǎn)單的任務(wù), 必須通過(guò)引入其他數(shù)值模式結(jié)果作為先驗(yàn)信息才有可能估計(jì)完整的二維波譜(Hasselmann, 1991; Hasselmann, 1996; Schulz-Stellenfleth, 2005)。或者利用SAR圖像交叉譜計(jì)算得到海浪譜, 此計(jì)算雖不需要初猜信息, 且可以消除海浪傳播方向180°模糊的問(wèn)題, 但產(chǎn)生的海浪譜無(wú)法完全解析高頻波, 即丟失短波信息, 多數(shù)情況下是涌浪部分的譜分量(Engen, 1995)。由此可見(jiàn), 通過(guò)SAR圖像譜獲取海浪譜, 進(jìn)一步獲得Hs等海浪統(tǒng)計(jì)參數(shù)這一途徑, 具有一定的困難性和局限性。

因此, 很多學(xué)者針對(duì)第二種方法, 即在不獲取二維海浪譜的情況下反演海浪參數(shù)的方法進(jìn)行了研究。Schulz-Stellenfleth等(2007)基于ERS2數(shù)據(jù), 提出了一個(gè)經(jīng)驗(yàn)二次模型CWAVE-ERS, 在不需要先驗(yàn)信息的條件下, 通過(guò)SAR圖像直接獲取海浪積分參數(shù)。Li等(2011)與Stopa等(2017)在CWAVE-ERS模型的基礎(chǔ)上, 分別提出了針對(duì)ENVISAT ASAR數(shù)據(jù)的CWAVE- ENV模型以及針對(duì)S1A數(shù)據(jù)的CWAVE-S1A模型, 用于計(jì)算有效波高、平均周期等海浪參數(shù)。但由于上述CWAVE算法, 僅針對(duì)特定的SAR數(shù)據(jù), 因此對(duì)于其他衛(wèi)星數(shù)據(jù)的普適性不強(qiáng)(Stopa, 2017)。

除以上兩種方法之外, 隨著海浪參數(shù)反演研究的不斷進(jìn)步, 許多學(xué)者認(rèn)識(shí)到方位向截?cái)嗖ㄩL(zhǎng)與海況條件之間具有強(qiáng)相關(guān)性, 可用于海浪統(tǒng)計(jì)參數(shù)反演的研究(Jackson, 1985; Kerbaol, 1998; Ren, 2016; Wang, 2012)。方位向截?cái)嗖ㄩL(zhǎng)為SAR數(shù)據(jù)所特有的海況參數(shù), 可作為對(duì)SAR方位向分辨率的一種度量(Greico, 2016)。在假設(shè)線性波的情況下, 方位向截?cái)嗖ㄩL(zhǎng)可表示為(Lyzenga, 1986; Kerbaol, 1998):

因此, 本文擬在相關(guān)有效波高反演研究的基礎(chǔ)上進(jìn)行拓展, 利用S1A二級(jí)波模式數(shù)據(jù), 采用神經(jīng)網(wǎng)絡(luò)的技術(shù)手段建立有效波高與0、c等SAR參數(shù)之間的聯(lián)系, 旨在分析各參數(shù)對(duì)有效波高反演的影響。本次實(shí)驗(yàn)建立的神經(jīng)網(wǎng)絡(luò)模型, 除用于模型訓(xùn)練的歐洲中程天氣預(yù)測(cè)中心(the European Centre for Medium-Range Weather Forecasts, ECMWF)模式數(shù)據(jù)外, 不需要額外模式數(shù)據(jù)作為先驗(yàn)信息, 且模型輸入的SAR參數(shù)可在S1A二級(jí)WV數(shù)據(jù)中直接獲取, 大大降低了數(shù)據(jù)存儲(chǔ)內(nèi)存, 減化了數(shù)據(jù)處理過(guò)程, 對(duì)于硬件條件受限的用戶尤為重要。

1 數(shù)據(jù)源介紹

本文采用了三種類型的數(shù)據(jù), 分別為S1A衛(wèi)星二級(jí)波模式數(shù)據(jù)、ECMWF再分析數(shù)據(jù)和美國(guó)國(guó)家數(shù)據(jù)浮標(biāo)中心(National Data Buoy Center, NDBC)浮標(biāo)數(shù)據(jù)。其中, S1A數(shù)據(jù)與ECMWF再分析數(shù)據(jù)進(jìn)行時(shí)空匹配, 用于神經(jīng)網(wǎng)絡(luò)模型的訓(xùn)練, 獨(dú)立測(cè)量的浮標(biāo)數(shù)據(jù)用于神經(jīng)網(wǎng)絡(luò)訓(xùn)練結(jié)果的驗(yàn)證。

1.1 Sentinel-1A 二級(jí)波模式數(shù)據(jù)

Sentinel-1A衛(wèi)星于2014年4月發(fā)射成功, 在太陽(yáng)同步軌道運(yùn)行, 高度()約為693km, 飛行速度約為7570m/s, 傾角98.18°, 重復(fù)周期為12d, 搭載了基于C波段的雷達(dá)成像系統(tǒng)。S1A衛(wèi)星具有四種成像模式, 其中波模式(Wave Mode, WV)是S1A在大洋上的操作模式, 主要應(yīng)用于海洋參數(shù)的獲取。在此模式下, S1A沿衛(wèi)星軌道每100km進(jìn)行一次采集, 以23°(WV1)和36°(WV2)兩個(gè)入射角交替工作, 在近距處獲取一個(gè)圖斑, 在遠(yuǎn)距處獲取下一個(gè)圖斑, 具有相同入射角的圖斑間隔為200km, 空間分辨率為5m, 圖斑大小約為20km×20km。

本研究使用的衛(wèi)星數(shù)據(jù)為S1A二級(jí)WV海洋產(chǎn)品(L2 Ocean Product, OCN)中的海洋涌浪譜(Ocean Swell spectra, OSW)數(shù)據(jù), 可通過(guò)ESA網(wǎng)站(https:// scihub.copernicus.eu/dhus/#/home)免費(fèi)下載獲取。該數(shù)據(jù)可提供0、NV、、skew、kurt、對(duì)數(shù)極坐標(biāo)系下的海浪方向譜、c、Hs等參數(shù)。但受到速度聚束調(diào)制作用的影響, 其測(cè)量的海浪Hs被嚴(yán)重低估, 且受方位向截?cái)嗖ㄩL(zhǎng)的影響較大, 僅可用于較低海況(Li, 2011)。如圖1所示, 我們將S1A二級(jí)數(shù)據(jù)提供的Hs以及根據(jù)公式(1)將S1A提供的海浪譜數(shù)據(jù)進(jìn)行積分獲得的Hs與ECMWF再分析數(shù)據(jù)提供的Hs進(jìn)行了比較, 可以看出, 圖1a、b中散點(diǎn)基本分布在=直線下方, 表明S1A二級(jí)數(shù)據(jù)提供的以及通過(guò)海浪譜計(jì)算的Hs低于ECMWF再分析數(shù)據(jù)中的Hs。

圖1 S1A二級(jí)波模式數(shù)據(jù)提供的Hs與ECMWF再分析數(shù)據(jù)提供的Hs的對(duì)比散點(diǎn)圖

注：a: S1A二級(jí)波模式數(shù)據(jù)提供的Hs與ECMWF再分析數(shù)據(jù)提供的Hs的比較散點(diǎn)圖; b: 通過(guò)S1A二級(jí)波模式數(shù)據(jù)提供的海浪譜積分獲得的Hs與ECMWF再分析數(shù)據(jù)提供的Hs的比較散點(diǎn)圖; Hs: 有效波高

為保證本研究實(shí)驗(yàn)結(jié)果不受船舶、海冰、溢油等影響, 選擇雷達(dá)截面的歸一化方差在1—2之間, 緯度在±60°, 時(shí)間為2018年1—12月的S1A數(shù)據(jù)作為實(shí)驗(yàn)數(shù)據(jù)。

1.2 ECMWF再分析數(shù)據(jù)

由于本研究選用的S1A波模式數(shù)據(jù)通常為全球公海數(shù)據(jù), 而在公海僅有少量浮標(biāo)測(cè)量數(shù)據(jù)可用。考慮到神經(jīng)網(wǎng)絡(luò)對(duì)訓(xùn)練數(shù)據(jù)高準(zhǔn)確度、廣覆蓋度的要求, 因此, 本研究同時(shí)選擇ECMWF再分析數(shù)據(jù)ERA5與S1A數(shù)據(jù)進(jìn)行時(shí)空匹配, 組成神經(jīng)網(wǎng)絡(luò)訓(xùn)練數(shù)據(jù)集。

ERA5是第五代ECMWF再分析數(shù)據(jù)集, 是ECMWF使用其預(yù)測(cè)模型和數(shù)據(jù)同化系統(tǒng)“重新分析”存檔的觀測(cè)結(jié)果, 用于創(chuàng)建描述大氣、陸地表面和海洋最近歷史的全球數(shù)據(jù)集, 可由ECMWF網(wǎng)站下載獲得(http://www. ecmwf.int/)。本研究主要采用ERA5海浪數(shù)據(jù)中的有效波高數(shù)據(jù), 其格網(wǎng)大小為0.5°×0.5°, 時(shí)間間隔為1h。本研究按照小于0.25°的空間窗口、小于30min的時(shí)間窗口對(duì)二者進(jìn)行時(shí)空匹配, 共獲得982351對(duì)匹配數(shù)據(jù)對(duì), 其中328081對(duì)用于模型訓(xùn)練, 其余數(shù)據(jù)用于模型應(yīng)用。

1.3 NDBC浮標(biāo)數(shù)據(jù)

為了對(duì)神經(jīng)網(wǎng)絡(luò)模型輸出結(jié)果進(jìn)行對(duì)比分析, 本文采用獨(dú)立測(cè)量的具有海浪譜信息的浮標(biāo)數(shù)據(jù)對(duì)其進(jìn)行驗(yàn)證, 其數(shù)據(jù)可由NDBC網(wǎng)站下載獲得(https://www.ndbc.noaa.gov/)。浮標(biāo)多布放于北半球, 數(shù)據(jù)獲取頻率約為1次/h, 根據(jù)其頻率測(cè)量范圍分為兩類: 0.03—0.4Hz, 離散為38個(gè)頻率; 0.02—0.485Hz, 離散為47個(gè)頻率。可根據(jù)公式(1), 將海浪譜進(jìn)行積分獲取有效波高。本研究將浮標(biāo)數(shù)據(jù)和S1A數(shù)據(jù)按照小于0.5°的空間窗口和小于0.5h的時(shí)間窗口進(jìn)行時(shí)空匹配, 獲得400對(duì)數(shù)據(jù)用于模型驗(yàn)證。

2 研究方法

本研究擬在Schulz-Stellenfleth等(2007)開(kāi)發(fā)的利用ERS2 SAR數(shù)據(jù)0、NV及其乘積項(xiàng)(02、NV2、0×NV)共5個(gè)參數(shù)用于海浪參數(shù)反演的雙參數(shù)模型基礎(chǔ)上, 加入經(jīng)緯度、c等信息進(jìn)行反演研究。同時(shí), 為增加模型的計(jì)算效率, 本文將利用神經(jīng)網(wǎng)絡(luò)技術(shù)進(jìn)行多元非線性回歸, 用于描述SAR不同輸入?yún)?shù)與有效波高之間的關(guān)系。首先將S1A數(shù)據(jù)與時(shí)空匹配后的ECWMF數(shù)據(jù)用于模型的訓(xùn)練, 然后利用獨(dú)立的浮標(biāo)數(shù)據(jù)對(duì)模型性能進(jìn)行驗(yàn)證, 同時(shí), 選取均值偏差(Bias)、均方根誤差、散射指數(shù)、相關(guān)系數(shù)四個(gè)統(tǒng)計(jì)參數(shù)對(duì)模型性能進(jìn)行定量評(píng)估。

2.1 神經(jīng)網(wǎng)絡(luò)的選取

神經(jīng)網(wǎng)絡(luò)是一種方便而通用的預(yù)測(cè)結(jié)果的方法, 近年來(lái)被多次應(yīng)用于遙感研究中。本研究基于MATLAB進(jìn)行神經(jīng)網(wǎng)絡(luò)訓(xùn)練, 選用的BP(back propagation)神經(jīng)網(wǎng)絡(luò), 是一種多層前饋神經(jīng)網(wǎng)絡(luò), 包括輸入層, 輸出層和隱藏層。其主要特點(diǎn)是信號(hào)前向傳遞, 誤差反向傳播。當(dāng)神經(jīng)網(wǎng)絡(luò)的輸入和輸出數(shù)據(jù)確定后, 可調(diào)整的只有隱藏層的神經(jīng)元數(shù)目, 而優(yōu)化隱藏層結(jié)點(diǎn)數(shù), 成為神經(jīng)網(wǎng)絡(luò)訓(xùn)練的首要任務(wù)。實(shí)驗(yàn)表明, 如隱藏層結(jié)點(diǎn)數(shù)過(guò)少, 網(wǎng)絡(luò)則不能具有必要的學(xué)習(xí)和信息處理能力。反之, 若過(guò)多, 不僅會(huì)增加網(wǎng)絡(luò)結(jié)構(gòu)的復(fù)雜性, 還會(huì)降低網(wǎng)絡(luò)學(xué)習(xí)的速度。因此, 在神經(jīng)網(wǎng)絡(luò)模型測(cè)試過(guò)程, 可對(duì)隱藏層結(jié)點(diǎn)數(shù)進(jìn)行調(diào)整, 最終取最優(yōu)結(jié)果來(lái)確定神經(jīng)元數(shù)目(王小川等, 2013; Stopa, 2017)。本文在進(jìn)行神經(jīng)網(wǎng)絡(luò)訓(xùn)練的實(shí)驗(yàn)中, 根據(jù)不同N_N模型輸入?yún)?shù)的個(gè)數(shù), 將隱藏層神經(jīng)元個(gè)數(shù)設(shè)定為4—30之間, 發(fā)現(xiàn)其性能均在隱藏層神經(jīng)元個(gè)數(shù)大于輸入?yún)?shù)個(gè)數(shù)時(shí)有所提升, 且在大于20時(shí)趨于穩(wěn)定, 但隨著神經(jīng)元數(shù)目增多, 訓(xùn)練時(shí)間也逐漸增加, 因此, 考慮到不同模型之間的對(duì)比, 本文將N_N模型的隱藏層神經(jīng)元數(shù)目均定為20。

基于輸入與輸出參數(shù)之間非線性關(guān)系的考慮, 神經(jīng)網(wǎng)絡(luò)采用正切S型傳遞函數(shù)(tansig)隱藏神經(jīng)元和線性傳遞函數(shù)(purelin)輸出神經(jīng)元。另外, 基于Levenberg-Marquardt反向傳播算法(trainlm)是現(xiàn)階段應(yīng)用最廣泛的優(yōu)化算法, 可為非線性最小二乘問(wèn)題提供解決方案, 本研究采用trainlm進(jìn)行神經(jīng)網(wǎng)絡(luò)訓(xùn)練。最后, 使用均方誤差(MSE)和回歸分析評(píng)估其性能(王小川等, 2013)。

2.2 N_N模型建立

為避免神經(jīng)網(wǎng)絡(luò)訓(xùn)練效果受到進(jìn)行不同范圍的Hs數(shù)據(jù)分層篩選等人為因素影響, 同時(shí)考慮到對(duì)全球范圍內(nèi)的有效波高反演的目的以及有效波高隨季節(jié)明顯變化的特性, 本文擬在2018年所有S1A數(shù)據(jù)以及時(shí)空匹配后的ECMWF的Hs數(shù)據(jù)中每個(gè)季節(jié)選擇一個(gè)月份的數(shù)據(jù)(本文選擇1、4、7、10月份的數(shù)據(jù))作為訓(xùn)練數(shù)據(jù), 同時(shí)擬采用常規(guī)神經(jīng)網(wǎng)絡(luò)模型訓(xùn)練方法, 在模型訓(xùn)練過(guò)程中將訓(xùn)練數(shù)據(jù)隨機(jī)分為60%、20%、20%, 分別用于模型建立過(guò)程中的訓(xùn)練、測(cè)試和驗(yàn)證。

本文首先對(duì)2018年4個(gè)月訓(xùn)練數(shù)據(jù)的ECMWF的Hs數(shù)據(jù)的分布及其占2018年總Hs數(shù)據(jù)的比例分布進(jìn)行了分析。1、4、7、10月份數(shù)據(jù)包括了全球范圍內(nèi)不同的海況, 且分別代表冬、春、夏、秋四個(gè)季節(jié), 因此最終選擇1、4、7、10月份的數(shù)據(jù)進(jìn)行訓(xùn)練, 并將2018年其他月份數(shù)據(jù)作為獨(dú)立的數(shù)據(jù)(應(yīng)用數(shù)據(jù))用于神經(jīng)網(wǎng)絡(luò)的應(yīng)用, 不參與模型的訓(xùn)練。訓(xùn)練數(shù)據(jù)的Hs分布情況以及在2018年所有數(shù)據(jù)中的比例分布情況的柱狀圖, 分別如圖2a和b所示, 可以看出, 其主要Hs范圍在1—4m, 且在2018年總數(shù)據(jù)中分布較均勻, 比例均在33%左右。

圖2 訓(xùn)練數(shù)據(jù)的Hs分布情況

注：a: 訓(xùn)練數(shù)據(jù)的Hs在每個(gè)區(qū)間內(nèi)的數(shù)據(jù)量分布情況; b: 訓(xùn)練數(shù)據(jù)的Hs與總數(shù)據(jù)在每個(gè)區(qū)間內(nèi)的比例, 其中, 由于較低海況以及高海況的數(shù)據(jù)較少, 因此將Hs按照以下方式進(jìn)行劃分: 橫軸代表不同的Hs區(qū)間, 共17個(gè), 1代表Hs<1m; 2—15代表19m

本文選用0、NV、02、NV2、0×NV、c、skew、kurt、lon、lat、作為神經(jīng)網(wǎng)絡(luò)訓(xùn)練模型的輸入?yún)?shù), 通過(guò)改變輸入?yún)?shù)的個(gè)數(shù), 進(jìn)行神經(jīng)網(wǎng)絡(luò)模型調(diào)整, 共建立14個(gè)N_N模型, 如表1所示。

表1 不同有效波高反演模型一覽表

Tab.1 List of different Hs inversion models

N_N模型輸入?yún)?shù)中,0通常與風(fēng)速和風(fēng)向有關(guān), 常用于風(fēng)場(chǎng)反演算法中, 因此可表示海浪的短波信息(Horstmann, 2003; Li, 2011)。NV提供關(guān)于圖像均勻性和海況的信息, 反映了由于SAR圖像中的波浪而引起的回波強(qiáng)度調(diào)制(Stopa, 2017)?？紤]到有效波高具有明顯的緯向分布的特征, 因此在進(jìn)行神經(jīng)網(wǎng)絡(luò)訓(xùn)練時(shí), 加入經(jīng)緯度信息, 增加輸入與輸出參數(shù)之間的空間對(duì)應(yīng)性。c反映SAR非線性成像機(jī)制, 與Hs具有很強(qiáng)的相關(guān)性(Beal, 1983; Shao, 2016)。另外, 有研究表明, 海表面的隨機(jī)運(yùn)動(dòng)對(duì)SAR成像的影響取決于SAR衛(wèi)星系統(tǒng)的表征參數(shù),值越大, SAR成像海洋表面波的能力就越差, 因此,信息可用于提供海表面運(yùn)動(dòng)的信息(Beal, 1983)。最后利用skew、kurt來(lái)考慮SAR圖像重要的高階特征信息(Stopa, 2017)。

2.3 統(tǒng)計(jì)分析

本研究對(duì)上述N_N模型反演得到的海浪有效波高, 進(jìn)行對(duì)比分析, 選擇Bias、RMSE、SI、CORR四種統(tǒng)計(jì)參數(shù)進(jìn)行評(píng)價(jià), 定義如下:

3 結(jié)果分析

為了更全面的對(duì)不同N_N模型的反演結(jié)果進(jìn)行分析, 本研究分別從模式數(shù)據(jù)、獨(dú)立測(cè)量的浮標(biāo)數(shù)據(jù)以及不同環(huán)境條件下的模型性能分析三個(gè)方面對(duì)海浪有效波高反演結(jié)果進(jìn)行了對(duì)比分析。

3.1 N_N模型反演結(jié)果與模式數(shù)據(jù)對(duì)比

為評(píng)估上述14個(gè)神經(jīng)網(wǎng)絡(luò)模型的性能, 本研究對(duì)應(yīng)用數(shù)據(jù)進(jìn)行了有效波高的反演, 對(duì)其結(jié)果進(jìn)行了統(tǒng)計(jì)分析。本節(jié)主要分析輸出結(jié)果與進(jìn)行時(shí)空配后的ECMWF數(shù)據(jù)之間的對(duì)比, 統(tǒng)計(jì)參數(shù)如表2所示。

表2 不同神經(jīng)網(wǎng)絡(luò)模型輸出的Hs與ECMWF數(shù)據(jù)提供的Hs對(duì)比統(tǒng)計(jì)參數(shù)一覽表

Tab.2 Comparison statistical parameters between Hs output from different N_N models and Hs provided by ECMWF

注：Bias: 均值偏差; RMSE: 均方根誤差; CORR: 相關(guān)系數(shù); SI: 散射指數(shù)

根據(jù)表2不同模型間的Bias、RMSE、CORR、SI四個(gè)統(tǒng)計(jì)參數(shù)對(duì)比可知:

1) N_N 2與N_N 1相比, 方位向截?cái)嗖ㄩL(zhǎng)的加入使神經(jīng)網(wǎng)絡(luò)模型的性能顯著提高, 散射指數(shù)明顯下降為19.93%, Bias與RMSE分別由-0.017m、0.654m降為-0.009m、0.534m, 相關(guān)系數(shù)由0.832提升至0.892, 這與Stopa等(2017)中的結(jié)果相近, 證明了方位向截?cái)嗖ㄩL(zhǎng)與有效波高之間的強(qiáng)相關(guān)性。

2) N_N3模型在N_N1模型的基礎(chǔ)上, 加入02、NV2、0×NV三個(gè)參量, 使得模型的統(tǒng)計(jì)參數(shù)略有提升, 均方根誤差的由0.654m降為0.643m, 以下模型均在N_N3基礎(chǔ)上進(jìn)行改進(jìn)。

3) N_N4、N_N5模型分別加入、c參數(shù), N_N6模型則同時(shí)加入、c參數(shù), N_N5與N_N6的統(tǒng)計(jì)參數(shù)相差較小, N_N4與N_N5相比, RMSE由0.557m降為0.532m, 這可能是由于c與成比例, 且其包含除信息外的更多信息的關(guān)系(Beal, 1983)。但N_N4與N_N3相比, 性能有明顯提升, 相關(guān)系數(shù)由0.838升為0.882, 因此在方位向截?cái)嗖ㄩL(zhǎng)信息缺失或者不準(zhǔn)確時(shí), 可以選擇信息作為有效波高反演模型的輸入?yún)?shù)。

4) N_N7、N_N8、N_N9與N_N10、N_N11、N_N12這兩組模型分別在N_N4和N_N5的基礎(chǔ)上依次加入skew、kurt、skew和kurt。通過(guò)分別對(duì)這兩組模型進(jìn)行分析, 可以看出, skew加入比kurt的加入, 對(duì)模型性能的改善更有優(yōu)勢(shì), 而二者的同時(shí)加入與單獨(dú)加入skew對(duì)模型的影響相近, 統(tǒng)計(jì)參數(shù)的差別均在千分之一數(shù)量級(jí)。另外, 在模型訓(xùn)練過(guò)程中, skew對(duì)模型的穩(wěn)定性有所提升, 而kurt的單獨(dú)加入會(huì)造成模型結(jié)果出現(xiàn)與ECMWF數(shù)據(jù)相差5m以上的奇異點(diǎn)。

5) N_N13是參照Schulz-Stellenfleth等(2007)的雙參數(shù)模型, 在N_N5基礎(chǔ)上加入方位向截?cái)嗖ㄩL(zhǎng)的乘積參數(shù)(c2、0×c、NV×c)。研究表明, N_N13和N_N5統(tǒng)計(jì)參數(shù)之間的差別均在千分之一的數(shù)量級(jí), RMSE下降了0.001m, 相關(guān)系數(shù)提升了0.001, 對(duì)模型的結(jié)果影響很小, 但輸入?yún)?shù)卻增加了8個(gè)。

6) N_N14在N_N10基礎(chǔ)上去掉了經(jīng)緯度信息, 旨在分析經(jīng)緯度信息對(duì)有效波高反演的重要性。研究表明, N_N14模型較N_N10模型, RMSE提升0.01m, SI增大0.36%, N_N10模型的性能優(yōu)于N_N14模型, 表明經(jīng)緯度信息對(duì)有效波高的反演具有一定的影響, 但影響很小。

綜上所述, 14個(gè)神經(jīng)網(wǎng)絡(luò)模型中, N_N10為最佳模型, 其性能最穩(wěn)定, 相關(guān)系數(shù)最高達(dá)到0.905, SI與RMSE最低, 分別為18.74%、0.502m。

為了更加直觀的對(duì)模型的性能進(jìn)行分析, 本研究將N_N模型的輸出結(jié)果與ECMWF數(shù)據(jù)的輸出結(jié)果的對(duì)比散點(diǎn)圖及二者之間的殘差全球分布圖展示如下, 由于篇幅限制, 本文僅對(duì)N_N5、N_N10模型的輸出結(jié)果與ECMWF數(shù)據(jù)的對(duì)比散點(diǎn)圖由圖3表示, 將N_N10模型與ECMWF數(shù)據(jù)的對(duì)比殘差全球分布圖由圖4表示。

通過(guò)與ECMWF數(shù)據(jù)相比, 由圖3a可知, N_N10模型訓(xùn)練結(jié)果較為理想, 相關(guān)系數(shù)接近0.91, SI為19.07%。圖3cN_N5模型與N_N10的散點(diǎn)圖分布相似度極高, 且相關(guān)系數(shù)也接近0.9, 但其Bias、RMSE均低于N_N10模型。圖3b的應(yīng)用效果也較好, Bias接近于零, RMSE為0.502m, SI低于20%。其數(shù)據(jù)點(diǎn)集中在1—4m之間, 且基本圍繞在=左右, 與Stopa等(2017)中的結(jié)果類似。當(dāng)Hs大于4m時(shí), 隨著訓(xùn)練數(shù)據(jù)點(diǎn)的減少, 訓(xùn)練效果也逐漸變差, 尤其當(dāng)Hs大于10m時(shí), 數(shù)據(jù)點(diǎn)基本在=以下, 這是由于在極端海況下, 數(shù)據(jù)點(diǎn)較少, 神經(jīng)網(wǎng)絡(luò)模型不能達(dá)到較好的學(xué)習(xí)效果造成的。

注：a: N_N10訓(xùn)練數(shù)據(jù)的輸出結(jié)果; b: N_N10應(yīng)用數(shù)據(jù)的輸出結(jié)果; c: N_N5訓(xùn)練數(shù)據(jù)的輸出結(jié)果; d: N_N5應(yīng)用數(shù)據(jù)的輸出結(jié)果; 藍(lán)色實(shí)線表示=

圖4 Hs的殘差及S1A應(yīng)用數(shù)據(jù)密度的全球分布圖

注：a: S1A數(shù)據(jù)通過(guò)N_N10模型得到的Hs與ECMWF再分析數(shù)據(jù)提供的Hs的殘差全球分布圖; b: 2018年S1A應(yīng)用數(shù)據(jù)全球分布密度圖

雖然圖3dN_N5模型的結(jié)果較圖3bN_N10模型稍差, 但二者的統(tǒng)計(jì)參數(shù)差別在百分之一的數(shù)量級(jí), 均值偏差接近, 相關(guān)系數(shù)僅相差0.013。也就是說(shuō), 即使無(wú)法獲取代表雷達(dá)截面高階特征的skew、kurt時(shí), 神經(jīng)網(wǎng)絡(luò)模型反演的Hs仍然可接受。

圖4a展示了S1A二級(jí)波模式數(shù)據(jù)通過(guò)N_N10模型反演得到的Hs與ECMWF模式數(shù)據(jù)提供的Hs之差的全球分布圖?？梢钥闯? 二者之間的殘差整體在±1m, 南半球殘差正值居多, 而北半球殘差多為負(fù)值, 且中高緯度地區(qū)殘差水平高于低緯度地區(qū), 整體色調(diào)為白色, 沿海地區(qū)有少量的紅色和藍(lán)色。數(shù)據(jù)分析得到, 二者之間的殘差小于二者的RMSE(0.502m)的數(shù)據(jù)占整體數(shù)據(jù)的90%左右, 而殘差水平大于3倍的RMSE (HsN_N10-HsECMWF>3*RMSE)的數(shù)據(jù)異常值僅為整體數(shù)據(jù)的0.1%, 異常值在全球分布且沒(méi)有空間相關(guān)性, 這與Stopa等(2017)的結(jié)果類似。表明整體殘差控制在理想的范圍之內(nèi)。由圖4b可知, 在太平洋海域, S1A數(shù)據(jù)最多, 沿海地區(qū)相對(duì)稀疏, 這是由于在一些沿海地區(qū), SAR天線經(jīng)常以其他模式(例如, 干涉寬幅模式)獲取數(shù)據(jù), 無(wú)法獲得波模式數(shù)據(jù)造成的(Li, 2011)。通過(guò)對(duì)比圖4a、b還可以觀察到, 殘差值的大小與應(yīng)用數(shù)據(jù)點(diǎn)的密度成正相關(guān)關(guān)系, 數(shù)據(jù)點(diǎn)密集區(qū)域, 殘差為正值, 數(shù)據(jù)點(diǎn)稀疏的區(qū)域, 殘差為負(fù)值, 這是由于訓(xùn)練數(shù)據(jù)點(diǎn)過(guò)少導(dǎo)致的神經(jīng)網(wǎng)絡(luò)模型不能達(dá)到較好的學(xué)習(xí)效果, 使得訓(xùn)練結(jié)果限制在一定范圍內(nèi), 相反, 數(shù)據(jù)點(diǎn)過(guò)于密集的區(qū)域, 可能存在過(guò)擬合現(xiàn)象。由圖4a可以看出, 殘差的高值分布區(qū)域, 一般為數(shù)據(jù)點(diǎn)比較少的區(qū)域, 證明了神經(jīng)網(wǎng)絡(luò)模型對(duì)訓(xùn)練數(shù)據(jù)的高要求。

3.2 N_N模型反演結(jié)果與浮標(biāo)數(shù)據(jù)對(duì)比

本節(jié)擬將N_N模型的反演結(jié)果與獨(dú)立測(cè)量的浮標(biāo)數(shù)據(jù)結(jié)果進(jìn)行比較, 由于篇幅限制, 本文只將N_N10模型的與浮標(biāo)數(shù)據(jù)的對(duì)比分析結(jié)果進(jìn)行展示, 如圖5所示。

圖5 浮標(biāo)數(shù)據(jù)提供的Hs與ECMWF再分析數(shù)據(jù)提供的Hs及S1A數(shù)據(jù)通過(guò)N_N10模型得到的Hs的對(duì)比散點(diǎn)圖

注：a: 浮標(biāo)數(shù)據(jù)提供的Hs與ECMWF再分析數(shù)據(jù)提供的Hs比較結(jié)果; b: 浮標(biāo)數(shù)據(jù)提供的Hs與S1A數(shù)據(jù)通過(guò)N_N10模型得到的Hs的比較結(jié)果; 藍(lán)色實(shí)線表示=

為了更好的分析N_N10模型輸出結(jié)果與浮標(biāo)數(shù)據(jù)的差別, 本研究首先對(duì)浮標(biāo)數(shù)據(jù)與時(shí)空匹配后的ECMWF模式數(shù)據(jù)進(jìn)行對(duì)比, 見(jiàn)圖5a。由于N_N10模型是以ECMWF數(shù)據(jù)為基準(zhǔn)進(jìn)行訓(xùn)練的神經(jīng)網(wǎng)絡(luò), 因此本文判斷其結(jié)果應(yīng)與ECMWF數(shù)據(jù)具有相似的趨勢(shì)。

由圖5a可知, 浮標(biāo)和ECMWF的有效波高整體上圍繞著=直線, 相關(guān)系數(shù)高達(dá)0.980, 散射系數(shù)為14.96%, 這表明了ECMWF再分析數(shù)據(jù)的可靠性及其用于提供可靠海況信息的能力。當(dāng)Hs<2m時(shí), 分布于=直線上方的數(shù)據(jù)點(diǎn)居多, 表明二者之間存在一個(gè)正偏差(HsECMWF-Hsbuoy), 數(shù)據(jù)點(diǎn)比較集中; 當(dāng)Hs>2m時(shí), 數(shù)據(jù)點(diǎn)分布在=直線左右, 數(shù)據(jù)點(diǎn)相對(duì)較為離散。

圖5b中散點(diǎn)圖分布趨勢(shì)與圖5a類似, 當(dāng)Hs<2m時(shí), 數(shù)據(jù)點(diǎn)基本分布于=直線上方, 當(dāng)Hs>2m時(shí), 數(shù)據(jù)點(diǎn)則離散的分布于=直線兩側(cè)。與圖5a相比, 圖5b的數(shù)據(jù)點(diǎn)則更加離散, SI值更大, 為32.27%, 圖5a、b的Bias分別為0.043m和0.263m, 均為正值, 這是由Hs<2m部分的數(shù)據(jù)所造成的。圖5b相關(guān)系數(shù)為0.894, 雖然沒(méi)有圖5a相關(guān)性高, 但也達(dá)到了較好的結(jié)果, 且二者之間相似的分布趨勢(shì)也進(jìn)一步證明N_N10模型提供了可靠的Hs估計(jì)。

3.3 不同環(huán)境條件下N_N10模型的性能分析

本節(jié)擬對(duì)不同SAR觀測(cè)模式(WV1/2), 不同的c以及不同海況下的N_N10模型的Hs反演情況進(jìn)行分析, 以有效波高的大小表示海況的變化。不同條件下的統(tǒng)計(jì)參數(shù)變化如表3所示。

表3 不同環(huán)境條件下的N_N10模型性能分析統(tǒng)計(jì)參數(shù)一覽表(N_N10 VS ECMWF)

Tab.3 List of statistical parameters for performance analysis of N_N10 Model under different environment conditions (N_N10 VS ECMWF)

由表3對(duì)不同入射角、Hs、c條件下的統(tǒng)計(jì)參數(shù)比較可知:

1) S1A衛(wèi)星在WV模式下以23°(WV1)和36°(WV2)兩個(gè)入射角交替工作, WV2模式的數(shù)據(jù)點(diǎn)較WV1的數(shù)據(jù)點(diǎn)略多, 但N_N10模型在WV2模式的反演結(jié)果比WV1稍差, RMSE相差0.029m, WV2的SI為19.22%, WV1的SI為18.19%, WV2的低性能是由于它的傾斜調(diào)制較弱, 信噪比較低, 在0和c中產(chǎn)生較大噪聲值的比率造成的(Stopa, 2017)。整體上N_N10模型在兩種SAR觀測(cè)模式下均取得了較好的反演結(jié)果。

2) 對(duì)于不同的Hs, N_N10模型性能具有明顯差異。其中, 平穩(wěn)海況(18m), N_N10模型性能大幅降低, 相關(guān)系數(shù)均在0.15以下。分析可知, 針對(duì)低海況的情況, 一方面由于用于進(jìn)行模型訓(xùn)練的低海況數(shù)據(jù)點(diǎn)較少, 神經(jīng)網(wǎng)絡(luò)的學(xué)習(xí)要求未達(dá)標(biāo), 另一方面, 當(dāng)Hs<1m時(shí), 存在一些N_N10模型輸出結(jié)果與ECMWF相差很大的數(shù)據(jù)點(diǎn)(由圖3可知)。而高海況下, 數(shù)據(jù)點(diǎn)更少, N_N10與N_N5模型均處于低估的狀態(tài)。

3) 方位向截?cái)嗖ㄩL(zhǎng)表示為SAR圖像譜在方位方向上受約束的程度, 當(dāng)方位向截?cái)嗖ㄩL(zhǎng)較大時(shí)(如,c>500m), SAR圖像譜將變得不清楚, 多數(shù)波譜成分將不可靠(Kerboal, 1998), 但N_N10模型在高c條件下, 卻取得了較好的結(jié)果, 相關(guān)系數(shù)高達(dá)0.944, Bias僅為0.102m, 而SI為35.94%表明數(shù)據(jù)較為離散。對(duì)于c小于100m的情況, 其Bias接近于0, 但相關(guān)系數(shù)很低, 為0.418, 分析可知, 由于c較小的情況下, 海面大多處于較為平穩(wěn)的狀態(tài), 訓(xùn)練數(shù)據(jù)較少, 訓(xùn)練效果不佳導(dǎo)致相關(guān)系數(shù)降低。

圖6進(jìn)一步對(duì)S1A的c、N_N10模型輸出的Hs、以及N_N10模型與ECMWF對(duì)比的殘差全球分布進(jìn)行展示, 從全球分布的角度, 分析三者之間的關(guān)系, 本研究選擇春秋兩個(gè)季節(jié)的數(shù)據(jù)進(jìn)行分析, 以1°的經(jīng)緯度窗口進(jìn)行季度平均計(jì)算。

從圖6a、b、d、e可知,c與Hs均具有緯向分布的特征, 且二者具有明顯正相關(guān)關(guān)系。春、秋兩個(gè)季節(jié), 南半球平均海浪有效波高均高于北半球, 尤其在30°—50°S, 平均Hs達(dá)到5m以上, 與c分布一致, 秋季則更為明顯。秋季時(shí)期北大西洋與南印度洋均出現(xiàn)了Hs高值區(qū), Hs在5m左右, 此處c也處于高值區(qū), 而春季時(shí)期, 南印度洋、南太平洋與南大西洋Hs較秋季時(shí)期減小, 但也為高值區(qū)與相應(yīng)的c的全球分布情況一致, 這也表明c可用于表示Hs的季節(jié)變化。從兩個(gè)季節(jié)的平均殘差分布圖可知, N_N10模型整體性能良好, 其殘差水平在沿岸地區(qū)高于大洋內(nèi)部, 且隨著c與Hs的提高, 殘差有增大的趨勢(shì), 且中高緯度地區(qū)的殘差水平高于低緯度地區(qū)。

圖6 S1A數(shù)據(jù)的λc、N_N10模型輸出的Hs及Hs殘差的全球分布圖

注：a、d表示S1A數(shù)據(jù)的c全球分布圖; b、e表示S1A數(shù)據(jù)通過(guò)N_N10模型得到的Hs分布圖; c、f表示N_N10模型輸出的Hs以及ECMWF再分析數(shù)據(jù)提供的Hs對(duì)比的殘差(HsN_N10-HsECMWF)分布情況; 以上數(shù)據(jù)均為季度平均的結(jié)果, 其中圖a、b、c為春季, 時(shí)間為2018年三月、五月, 圖d、e、f為秋季, 時(shí)間為2018年九月、十一月

4 結(jié)論

本文在Schulz-Stellenfleth等(2007)提出的雙參數(shù)模型的基礎(chǔ)上, 利用神經(jīng)網(wǎng)絡(luò)技術(shù)對(duì)2018年的S1A二級(jí)波模式數(shù)據(jù)進(jìn)行了海浪有效波高的反演, 對(duì)具有不同輸入?yún)?shù)的神經(jīng)網(wǎng)絡(luò)模型的性能進(jìn)行了分析。為了更準(zhǔn)確的進(jìn)行全球范圍的Hs反演, 本文在N_N模型中加入經(jīng)緯度的信息。對(duì)模型反演結(jié)果, 除利用ECMWF再分析數(shù)據(jù)進(jìn)行對(duì)比外, 還利用獨(dú)立測(cè)量的浮標(biāo)數(shù)據(jù)進(jìn)行了驗(yàn)證, 初步得到以下結(jié)論:

1) 總體上, 14個(gè)N_N模型的反演結(jié)果與ECMWF數(shù)據(jù)相比, 均達(dá)到了可接受的訓(xùn)練結(jié)果, 相關(guān)系數(shù)均在0.8以上, 且大部分模型RMS控制在0.6m以內(nèi), 其中N_N10的性能最佳, 相關(guān)系數(shù)達(dá)到0.9以上, 散射指數(shù)在19%以內(nèi)。

2) 通過(guò)表2對(duì)14個(gè)模型性能的分析, 隨著c的加入, N_N模型的性能大幅提升, 相關(guān)系數(shù)達(dá)到0.88以上。其次,參數(shù)的加入也會(huì)小幅度提高模型反演的準(zhǔn)確度, 因此, 在c缺失或者不準(zhǔn)確時(shí), 可選擇信息作為Hs反演模型的輸入?yún)?shù)。

3) 浮標(biāo)測(cè)量數(shù)據(jù)用于驗(yàn)證N_N模型算法, N_N10與浮標(biāo)對(duì)比的RMS為0.607m, Bias為0.263m, 且相關(guān)系數(shù)達(dá)到了0.894, 表明將N_N10模型用于Hs的反演, 整體上達(dá)到了較好的結(jié)果。同時(shí), 通過(guò)散點(diǎn)圖可知, 當(dāng)Hs<2m時(shí), N_N10模型、ECMWF均與浮標(biāo)觀測(cè)結(jié)果之間存在較小的正偏差; 當(dāng)Hs>2.5m時(shí), 數(shù)據(jù)點(diǎn)的偏差趨于0, 但數(shù)據(jù)點(diǎn)相對(duì)較為離散, N_N10與浮標(biāo)的散點(diǎn)圖分布情況與ECMWF和浮標(biāo)的對(duì)比散點(diǎn)圖具有相似趨勢(shì), 也進(jìn)一步表明N_N10模型提供了可靠的Hs估計(jì)。由于本研究用于評(píng)估的浮標(biāo)數(shù)據(jù)僅有400對(duì), 數(shù)量較少, 因此, 為了更好的確認(rèn)N_N性能, 需要進(jìn)行更密集的現(xiàn)場(chǎng)浮標(biāo)測(cè)量。

4) 通過(guò)對(duì)不同入射角下N_N10模型的性能分析, 在入射角為36°(WV2)模式的反演結(jié)果比23°(WV1)較差, RMSE相差0.029m, WV1的SI控制在19%以內(nèi)。但總體上N_N10模型在兩種SAR觀測(cè)模式下均取得了不錯(cuò)的反演結(jié)果。

5) 對(duì)于不同海況條件, N_N10模型性能具有明顯差異。其中, 平穩(wěn)海況(18m)的模型性能大幅降低, 相關(guān)系數(shù)均在0.15以下, 這是由于用于模型訓(xùn)練的低海況及高海況的數(shù)據(jù)點(diǎn)較少, 神經(jīng)網(wǎng)絡(luò)的學(xué)習(xí)要求未達(dá)標(biāo)所造成, 但如果通過(guò)加入更大的反映所有海況下的數(shù)據(jù)集, 該模型性能將會(huì)有所改進(jìn)。

6) 歷史研究表明, 較大的c, 會(huì)引起SAR圖像譜, 在方位方向上變得十分模糊, 很難用于海浪參數(shù)的提取(Kerbaol, 1998)。但本研究通過(guò)N_N10模型進(jìn)行Hs反演, 在高c條件下, 卻取得了較好的結(jié)果, 相關(guān)系數(shù)高達(dá)0.944。同時(shí), 通過(guò)對(duì)S1A數(shù)據(jù)的c以及N_N10模型反演的Hs進(jìn)行春、秋兩個(gè)季節(jié)的全球分布分析可知,c與Hs的全球分布具有大致相同的緯向分布的特征, 表明c還可用來(lái)表示Hs的季節(jié)變化。

7) 通過(guò)對(duì)N_N10與ECMWF的Hs殘差全球分布圖可知, N_N10模型在全球大部分海域達(dá)到了較好的反演效果, 但在沿岸地區(qū)的有效波高反演結(jié)果仍不理想, 遠(yuǎn)不及遠(yuǎn)海地區(qū), 該問(wèn)題與數(shù)據(jù)點(diǎn)的分布密度有關(guān)。因此, 通過(guò)加入包含更多沿岸數(shù)據(jù)的訓(xùn)練數(shù)據(jù)用于神經(jīng)網(wǎng)絡(luò)的訓(xùn)練, 可以完善N_N模型在全球范圍內(nèi)Hs反演的性能。

綜上分析, 本研究建立的用于Hs反演的神經(jīng)網(wǎng)絡(luò)模型, 取得了較好的反演效果, 在全球范圍內(nèi)提供了可靠的海浪有效波高估計(jì)。盡管本研究?jī)H限于S1A衛(wèi)星, 但該方法還可應(yīng)用于當(dāng)前任何具有二級(jí)波模式數(shù)據(jù)的其他合成孔徑雷達(dá)數(shù)據(jù), 如ENVISAT、Sentinel 1B等。此外, 若在模型訓(xùn)練過(guò)程中考慮海浪的傳播方向以及主波波長(zhǎng)等信息, 模型性能將會(huì)有所改進(jìn)(Stopa, 2017)。

致謝 本文使用的Sentinel 1A二級(jí)波模式數(shù)據(jù)由歐洲太空局(ESA)網(wǎng)站提供, ECMWF再分析數(shù)據(jù)集ERA5由ECMWF網(wǎng)站提供, 浮標(biāo)數(shù)據(jù)由美國(guó)國(guó)家數(shù)據(jù)浮標(biāo)中心(NDBC)網(wǎng)站下載, 在此表示感謝。

王小川, 史 峰, 郁 磊等, 2013. MATLAB神經(jīng)網(wǎng)絡(luò)43個(gè)案例分析. 北京: 北京航空航天大學(xué)出版社, 1—32

Beal R C, Tilley D G, Monaldo F M, 1983. Large-and small-scale spatial evolution of digitally processed ocean wave spectra from SEASAT synthetic aperture radar. Journal of Geophysical Research: Oceans, 88(C3): 1761—1778

Engen G, Johnsen H, 1995. SAR-ocean wave inversion using image cross spectra. IEEE Transactions on Geoscience and Remote Sensing, 33(4): 1047—1056

Grieco G, Lin W, Migliaccio M, 2016. Dependency of the Sentinel-1 azimuth wavelength cut-off on significant wave height and wind speed. International Journal of Remote Sensing, 37(21): 5086—5104

Hasselmann K, Hasselmann S, 1991. On the nonlinear mapping of an ocean wave spectrum into a synthetic aperture radar image spectrum and its inversion. Journal of Geophysical Research: Oceans, 96(C6): 10713—10729

Hasselmann K, Raney R K, Plant W J, 1985. Theory of synthetic aperture radar ocean imaging: A MARSEN view. Journal of Geophysical Research: Oceans, 90(3): 4659—4686

Hasselmann S, Brüning C, Hasselmann K, 1996. An improved algorithm for the retrieval of ocean wave spectra from synthetic aperture radar image spectra. Journal of Geophysical Research: Oceans, 101(C7): 16615—16629

Horstmann J, Schiller H, Schulz-Stellenfleth J, 2003. Global wind speed retrieval from sar. IEEE Transactions on Geoscience and Remote Sensing, 41(10): 2277—2286

Jackson C R, Apel J R, 2004. Synthetic aperture radar marine user’s manual. U.S. Department of Commerce National Oceanic and Atmospheric Administration. http:// sarusersmanual.com

Jackson F C, Walton W T, Baker P L, 1985. Aircraft and satellite measurement of ocean wave directional spectra using scanning-beam microwave radars. Journal of Geophysical Research: Oceans, 90(C1): 987—1004, doi: 10.1029/ JC090iC01p00987

Kerbaol V, Chapron B, Vachon P W, 1998. Analysis of ERS-1/2 synthetic aperture radar wave mode imagettes. Journal of Geophysical Research: Oceans, 103(C4): 7833—7846

Li X M, Lehner S, Bruns T, 2011. Ocean wave integral parameter measurements using Envisat ASAR wave mode data. IEEE Transactions on Geoscience and Remote Sensing, 49(1): 155—174

Lyzenga D R, 1986. Numerical simulation of synthetic aperture radar image spectra for ocean waves. IEEE Transactions on Geoscience and Remote Sensing, GE-24(6): 863—872

Marghany M, Ibrahim Z, Van Genderen J, 2002. Azimuth cut-off model for significant wave height investigation along coastal water of Kuala Terengganu, Malaysia. International Journal of Applied Earth Observation and Geoinformation, 4(2): 147—160

Ren L, Yang J S, Zheng G, 2015. Significant wave height estimation using azimuth cutoff of C-band RADARSAT-2 single-polarization SAR images. Acta Oceanologica Sinica, 34(12): 93—101

Schulz-Stellenfleth J, K?nig T, Lehner S, 2007. An empirical approach for the retrieval of integral ocean wave parameters from synthetic aperture radar data. Journal of Geophysical Research: Oceans, 112(C3): C03019

Schulz-Stellenfleth J, Lehner S, Hoja D, 2005. A parametric scheme for the retrieval of two‐dimensional ocean wave spectra from synthetic aperture radar look cross spectra. Journal of Geophysical Research: Oceans, 110(C5): C05004, doi: 10.1029/2004JC002822

Shao W Z, Zhang Z, Li X F, 2016. Ocean wave parameters retrieval from Sentinel-1 SAR imagery. Remote Sensing, 8(9): 707, doi: 10.3390/rs8090707

Stopa J E, Ardhuin F, Chapron B, 2015. Estimating wave orbital velocity through the azimuth cutoff from space-borne satellites. Journal of Geophysical Research: Oceans, 120(11): 7616—7634

Stopa J E, Mouche A, 2017. Significant wave heights from Sentinel-1 SAR: validation and applications. Journal of Geophysical Research: Oceans, 122(3): 1827—1848

Wang H, Zhu J H, Yang J S, 2012. A semiempirical algorithm for SAR wave height retrieval and its validation using Envisat ASAR wave mode data. Acta Oceanologica Sinica, 31(3): 59—66

INVERSION OF GLOBAL SIGNIFICANT WAVE HEIGHT BASED ON SENTINEL-1A

MU Shan-Shan1, LI Hai-Yan1, WU Ming-Bo1, 2

(1. University of Chinese Academy of Sciences, Beijing 100049, China; 2. Institute of Geographic Sciences and Natural Resources Research, Chinese Academy of Science, Beijing 100101, China)

By using the technique of neural network, a new neural network model a new neural network model (N_N model) was proposed for the inversion of wave effective wave height (HS) from Sentinel-1A level-2 wave model data. Based on the two-parameter model developed by ERS2 SAR wave mode data, the model was added with other parameters including longitude, latitude, the azimuth cutoff (c), skewness, kurtosis, and the ratio () of the distance between satellite platform and the target to the satellite flight speed. The influence of each parameter on the inversion of significant wave height was analyzed using different combinations of input parameters, based on 14 models which were established for HSinversion. Results show that the correlation coefficients of all the 14 models were above 0.8. With the addition ofcand β parameters, the performance of the N_N model increased significantly, and the improvement effect ofcon the model performance was more obvious. The correlation coefficient increased by about 0.06, and RMSE decreased by about 0.12m. In addition, the addition of skewness and kurtosis also improved the performance of the N_N model as the RMSE decreased by about 0.03m, and the correlation coefficient increased by about 0.01. Among them, the N_N10 model had the best effect and the most stable performance. Compared with the ECMWF (European Centre for medium range weather forecasts), the correlation coefficient (CORR) was 0.905, and the scattering index (SI) and RMSE were the lowest, being 18.74% and 0.502m, respectively. The correlation coefficient with the independently measured buoy data reached 0.894.

neural network; significant wave height; azimuth cutoff; normalized radar cross-section

2019-12-07

P7; TP3

李海艷, 碩士生導(dǎo)師, 副教授, E-mail: lihaiyan@ucas.edu.cn

2019-09-21,

* 國(guó)家自然科學(xué)基金項(xiàng)目, 41776197號(hào); 中國(guó)科學(xué)院大學(xué)優(yōu)秀青年教師科研能力提升項(xiàng)目, Y95401N號(hào)。穆珊珊, 碩士研究生, E-mail: mushanjiangnan@163.com

10.11693/hyhz20190900177