李秀仲, 何宜軍, 孟俊敏, 張振華

(1. 南京信息工程大學(xué), 江蘇 南京 210044; 2. 國家海洋局第一海洋研究所, 山東 青島 266061; 3. 北京遙測技術(shù)研究所, 北京 100076)

機(jī)載波譜儀海浪譜反演方法及其驗(yàn)證

李秀仲1, 何宜軍1, 孟俊敏2, 張振華3

(1. 南京信息工程大學(xué), 江蘇 南京 210044; 2. 國家海洋局第一海洋研究所, 山東 青島 266061; 3. 北京遙測技術(shù)研究所, 北京 100076)

為了驗(yàn)證波譜儀反演二維海浪譜的功能, 根據(jù)海浪波譜儀的信號形成機(jī)制, 總結(jié)了機(jī)載波譜儀反演海浪的流程。利用機(jī)載波譜儀回波數(shù)據(jù), 通過自相關(guān)和互相關(guān)兩種功率譜估計方法, 反演了二維海浪譜。最后通過與浮標(biāo)測量的二維海浪譜進(jìn)行對比, 驗(yàn)證了該機(jī)載波譜儀探測二維海浪譜的有效性。結(jié)果表明, 無論采用自相關(guān)函數(shù)還是互相關(guān)函數(shù)進(jìn)行功率譜估計, 得到的主波波長和有效波高與實(shí)際二維海浪譜基本一致?；ハ嚓P(guān)函數(shù)法得到的交叉譜能去除180°模糊現(xiàn)象, 其在計算有效波高時相對于自相關(guān)函數(shù)會稍微偏小。在計算斜率方差時可以采用5°~12°入射角范圍的后向散射系數(shù)進(jìn)行公式擬合, 因此定標(biāo)與否并不影響最后的二維海浪譜結(jié)果, 未來星載波譜儀只有靠多波束聯(lián)合才能實(shí)現(xiàn)。

波譜儀; 功率譜估計; 二維海浪譜; 180°模糊

海況信息的獲取有著重要的業(yè)務(wù)化需求(如海洋氣象預(yù)報、導(dǎo)航、近海和沿海活動等)和科研需求(海浪動力分析、大氣和海洋邊界層的相互作用、區(qū)域和整體的海洋和大氣系統(tǒng)耦合, 電磁信號與海面相互作用等)。過去幾十年, 雖然海浪預(yù)報已經(jīng)進(jìn)步很大[1], 但是在高海況條件、極值事件(如臺風(fēng)、風(fēng)暴潮等)方面還有很大的改進(jìn)余地; 在提供準(zhǔn)確參數(shù)預(yù)報方面不僅可以預(yù)報有效波高, 還可以改進(jìn)峰值波長和波向預(yù)報。目前應(yīng)用的海浪數(shù)值模式在預(yù)測二維海浪譜方面面臨很大缺點(diǎn)[1], 比如海浪激發(fā)源和耗散項的參數(shù)化以及采用的初始場存在誤差較大等。此外, 由于數(shù)值模型離散化的限制, 在海浪能量的方向分布上, 開闊海域時會產(chǎn)生較大誤差, 這樣就會導(dǎo)致波浪在靠近海岸時的演化結(jié)果中產(chǎn)生更大的預(yù)報誤差。

目前唯一星載測量海浪的雷達(dá)系統(tǒng)是合成孔徑雷達(dá)(Synthetic Aperture Radar, SAR)。自從Seasat衛(wèi)星上搭載了SAR后[2], SAR就一直比較昂貴并且在測量海浪方向譜時并不簡單易用。為了去除SAR這種限制, 人們發(fā)展了其它類型雷達(dá): Walsh等[3]1985開發(fā)了Ka波段的海面等值線雷達(dá)(surface contour radar, SCR), Jackson等[4-5]1985年開發(fā)了Ku波段雷達(dá)海浪波譜儀(Radar Ocean Wave Spectrometer, ROWS), Hauser等[6]1992年開發(fā)了C波段雷達(dá)測波儀(Radar pour l’Etude du Spectre des Surfaces par AnalyseCirculaire, RESSAC), Hauser等[7]2003年開發(fā)了C波段極化雷達(dá)觀測系統(tǒng)(Système de Télédétection pour l’Observation par Radar de la Mer, STORM), 以及Plant等[8]2005年開發(fā)了X波段相干真實(shí)孔徑雷達(dá)(coherent real aperture radar, CORAR)。這些機(jī)載測量雷達(dá)比方向浮標(biāo)提供了更多的海浪方向性信息[9]。

中法海洋衛(wèi)星(China–France Oceanography Satellite, CFOSAT)是一個由中國和法國空間機(jī)構(gòu)正在準(zhǔn)備的新衛(wèi)星。它是為了同步海浪譜和海面風(fēng)場觀測而設(shè)計的。該衛(wèi)星上將搭載一個Ku波段的散射計測量風(fēng)速, Ku波段的波譜儀(Surface Waves Investigation and Monitoring, SWIM)來測量海面浪場。CFOSAT是一個極軌衛(wèi)星, 飛行高度為520 km。SWIM雷達(dá)能向海面發(fā)射6個不同入射角的波束: 0°、2°、4°、6°、8°和10°, 每個波束天線孔徑都是2°,天線旋轉(zhuǎn)速度是5.6 rad/min。SWIM能提供如下的地球物理參數(shù): 0°~10°入射角的標(biāo)準(zhǔn)化后向散射系數(shù)(normalized radar cross-section, NRCS); 星下點(diǎn)的有效波高和風(fēng)速; 6°、8°、10°波束能提供海浪方向譜。

在CFOSAT的準(zhǔn)備階段, 北京遙測技術(shù)研究所開發(fā)了機(jī)載波譜儀, 并在山東半島以南的黃海區(qū)域進(jìn)行了校飛試驗(yàn), 以檢驗(yàn)SWIM測量海浪方向譜的能力。本文針對機(jī)載波譜儀70圈的回波數(shù)據(jù), 通過總結(jié)國內(nèi)外反演算法, 反演了二維海浪方向譜, 并與浮標(biāo)實(shí)測的二維海浪譜對比, 檢驗(yàn)了波譜儀測量海浪二維譜的有效性, 為星載波譜儀的反演算法的開發(fā)奠定了基礎(chǔ)。

圖1 處理的波譜儀校飛試驗(yàn)數(shù)據(jù)所在軌跡Fig. 1 Tracks of the processed data obtained using an airborne sea-wave spectrometer

表1 機(jī)載波譜儀主要參數(shù)Tab. 1 Main parameters of the airborne sea-wave spectrometer

1 機(jī)載波譜儀校飛介紹

機(jī)載波譜儀的校飛試驗(yàn)位置選在黃海北部, 圖1中三角形和六角形分別為2014年6月12日10: 21和2014年6月21日11: 07浮標(biāo)所在位置(分別是122.84°E、36.62°N和122.61°E、36.62°N), 浮標(biāo)每半小時給出一個海浪譜結(jié)果。校飛時無風(fēng)速計實(shí)測風(fēng)速, 但根據(jù)浮標(biāo)所測海浪波高較小可知, 兩個架次的風(fēng)速在2~5 m/s變化。三角形附近的粗曲線為飛機(jī)2014年6月12日11: 6: 20~11: 11: 20共301 s的飛行軌跡, 也是本文處理的30圈回波數(shù)據(jù)的軌跡, 始末位置相差15.62 km, 起始位置與浮標(biāo)相差7.17 km;細(xì)直線為飛機(jī)2014年6月21日10: 48: 58~10: 55: 37共400 s的飛行軌跡, 是本文處理的40圈的回波數(shù)據(jù)軌跡, 始末位置相差28.54 km, 與六角形所代表的浮標(biāo)最近距離約為11.58 km。

對于機(jī)載波譜儀, 由于飛行高度不可能達(dá)到衛(wèi)星的高度, 因此采用一個大波束, 以增加足印面積。該機(jī)載波譜儀的主要參數(shù)在表1中給出。

圖2顯示了此次機(jī)載波譜儀校飛的幾何示意圖,根據(jù)3 dB波束寬度, 每個脈沖在海面上形成的等天線增益值線可以簡化成一個橢圓形狀。根據(jù)傾斜調(diào)制的原理, 雷達(dá)接收圖中橢圓內(nèi)圓環(huán)里的海面反射能量, 當(dāng)波浪局地傾斜導(dǎo)致的海面法線對著雷達(dá)時海面回波強(qiáng), 反之當(dāng)波浪局地傾斜導(dǎo)致的海面法線偏離雷達(dá)時, 雷達(dá)接收的回波弱。因此雷達(dá)信號變化的幅值與海浪局地傾角有關(guān)。根據(jù)Jackson等[4-5]、Hauser等[10], 海面斜率與后向散射系數(shù)的關(guān)系為:

圖2 機(jī)載波譜儀觀測幾何Fig. 2 Geometry of observation of the airborne sea-wave spectrometer

式中, p為海面斜率概率密度函數(shù), a為定標(biāo)參數(shù),x為波面高度, q為入射角。因此獲取海面斜率的雷達(dá)調(diào)制信號, 求取調(diào)制信號功率譜, 可得某一方向的海浪斜率譜通過雷達(dá)天線360°水平面旋轉(zhuǎn), 即可得到二維方向譜, 即公式(3)中的

式中，Ly為波束方位向3 dB寬度,為調(diào)制譜, k為波數(shù), j為方位角。

圖3模擬了機(jī)載波譜儀校飛1 min形成的足跡。波譜儀在海面形成的足跡與以下幾個參數(shù)有關(guān): 飛機(jī)飛行高度、飛機(jī)飛行速度、天線旋轉(zhuǎn)速度以及脈沖重復(fù)頻率。為了顯示清晰, 該圖只顯示了5 Hz的脈沖重復(fù)頻率結(jié)果, 飛行高度為3 000 m、飛機(jī)飛行速度為60 m/s、天線轉(zhuǎn)速是6 rad/min?？紤]到準(zhǔn)鏡面散射理論適用范圍以及其水平分辨率因素, 圖中展示了入射角5°~12°時的波束足印, 水平徑向距離為375 m, 短半軸長則為214 m。

圖3 機(jī)載波譜儀探測足跡Fig. 3 Footprint detected by the sea-wave spectrometer

2 機(jī)載海浪譜反演方法

本文根據(jù)波譜儀的工作原理, 在國內(nèi)外波譜儀研究基礎(chǔ)上, 采用以下流程進(jìn)行了二維海浪譜的反演。

圖4中, 回波數(shù)據(jù)是雷達(dá)接收的不同時刻電磁波能量, 天線增益是由實(shí)驗(yàn)室實(shí)測后的二維天線增益經(jīng)過飛機(jī)姿態(tài)角校正后得到的增益[11], r代表斜距, m( r)是沿斜距的調(diào)制信號, 由本文公式(4)中的組成, x是水平距離,是沿水平方向的調(diào)制信號。此處海浪斜率譜為由調(diào)制譜根據(jù)公式(3)轉(zhuǎn)換得到斜率譜為標(biāo)準(zhǔn)化過程,j為方位角。

2.1 計算0σ

波譜儀獲取回波能量之后, 為了獲取海浪的無量綱的標(biāo)準(zhǔn)化后向散射系數(shù)0s, 采用雷達(dá)方程為

式中,rP和tP是接收和發(fā)射的能量, R是斜距,tG和rG是天線發(fā)射和接收的增益, l是單位為米的載頻電磁波長。雷達(dá)接收的截面(Radar Cross-Section, RCS)s用公式(5)表達(dá)

圖4 機(jī)載波譜儀海浪譜反演方法Fig. 4 Retrieving method of sea-wave spectrum using an airborne spectrometer

S是電磁波波束照射的面積, 這個面積能夠用一個矩形的寬度dx和長度dy來表示, 并且是方位向的波束寬度, 所以。將此公式代入公式(4), 并令可得到公式(6)

2.2 計算斜率方差和調(diào)制系數(shù)

基于小入射角下NRCS, 有一個簡便的方法能用來計算斜率方差(mean square slope, MSS)。Barrick[12]和Valenzuela[13]導(dǎo)出了有限粗糙面的電磁散射公式, 總結(jié)出小入射角范圍是0°~15°時準(zhǔn)鏡面散射是海面散射的主要部分。準(zhǔn)鏡面散射公式是此處天底點(diǎn)的菲尼爾反射系數(shù),是二維海面斜率概率密度分布函數(shù), zx表示距離向海面斜率, zy表示方位向海面斜率。正比于滿足條件的微小面元的數(shù)量, 所以也可以寫成。如果是高斯分布, 那么等式(8)變成

此處us和cs分別是海面斜率在逆風(fēng)向和側(cè)風(fēng)向的標(biāo)準(zhǔn)偏差。MSS是方位角j處的斜率方差。如果海面斜率概率分布函數(shù)(Probability Density Function, PDF)是高斯的, 那么我們可以得到[14-15]

所以, 當(dāng)0s通過式(7)獲取后, 利用方位角j處的回波對上式擬合, 可以獲取MSS。此處選擇5°入射角作為擬合上限是考慮了地面分辨率隨著入射角減小而變差這個因素。Freilich等[16]假設(shè)入射角0°~18°時降雨雷達(dá)散射能量在海面風(fēng)速0~20 m/s范圍內(nèi)主要是由準(zhǔn)鏡面散射機(jī)制產(chǎn)生的, 其中較大入射角時布拉格散射機(jī)制變得相對重要了。此處選擇12°入射角作為上限并假設(shè)此時準(zhǔn)鏡面散射占主要部分是合適的。需要注意的是, 我們獲取的是系數(shù)的負(fù)倒數(shù), 所以0s是否進(jìn)行了定標(biāo)并沒有關(guān)系。

根據(jù)海面斜率分布近似高斯函數(shù)形式的結(jié)論,則通過公式(2)調(diào)制系數(shù)計算公式變?yōu)?/p>

由于海浪譜是由信號調(diào)制譜得來的, 因此計算海浪譜之前, 需要通過調(diào)制信號計算調(diào)制譜。采用兩種方法估計調(diào)制譜, 第一種是自相關(guān)函數(shù)功率譜估計方法

該方法是采用同一信號進(jìn)行自相關(guān)得到的功率譜, 在進(jìn)行譜估計之前, 需要對調(diào)制信號進(jìn)行加窗處理, 以消除由于信號截斷造成的頻譜泄露問題。第二種是交叉譜估計方法

j￠為同圈內(nèi)與j有一定時間間隔的方位角, 交叉譜計算時, 主要針對相鄰的兩個回波, 以達(dá)到去除噪聲的目的。

這兩種方法求取的海浪譜都存在180°模糊的問題, 本文借鑒SAR消除180°模糊的方法[17], 根據(jù)公式(13)中求取的互相關(guān)函數(shù), 則觀測到的海浪傳播的徑向軌道速度為

2.4 海浪參數(shù)的確定

有效波高(significant wave height, SWH)HSW代表了海浪能量的大小, 其計算公式為

主波波數(shù)kp為F(k,j)最大值所對應(yīng)的波數(shù), 主波波長則為主波波向?yàn)榉较蜃V取最大值處對應(yīng)的方向j。在求取浮標(biāo)所測海浪譜主波波長和主波波向時, 對于浮標(biāo)提供的頻率譜, 需要轉(zhuǎn)換成波數(shù)譜。根據(jù)

w是角頻率, f為波浪頻率, h和g分別為水深和重力加速度, 水深h在30 m左右, 根據(jù)上式可以得到每個頻率對應(yīng)的波數(shù)。

3 二維海浪譜及其參數(shù)的驗(yàn)證

浮標(biāo)測量的二維海浪譜如圖5、圖6所示, 其中HSW分別為0.47, 0.6 m, 主波波長是83, 39 m, 波浪傳播方向是南偏東20°, 12°, 若正北為0°, 逆時針為正, 則浪向?yàn)?00°, 192°。

圖5 2014年6月12日10: 21波浪騎士所測二維海浪譜Fig. 5 Two-dimensional spectrum detected from a sea wave rider at 10: 21 on June 12, 2014

圖6 2014年6月21日11: 07波浪騎士所測二維海浪譜Fig. 6 Two-dimensional spectrum detected from a sea wave rider at 11: 07 on June 21, 2014

由于本次校飛時海況較低, 因此海浪信息并不夠規(guī)則, 此處將這兩塊數(shù)據(jù)所得結(jié)果各自進(jìn)行了平均。利用前面兩種功率譜計算方法得到兩種海浪譜反演結(jié)果如下圖, 對于第一塊數(shù)據(jù), 自相關(guān)函數(shù)方法得到的結(jié)果顯示, 主波波長93 m, 波浪傳播方向?yàn)?0°, 與浮標(biāo)結(jié)果相反, 有效波高為0.46 m。而交叉譜方法得到的這3個參數(shù)分別為93 m、210°、0.41 m, 海浪傳播方向基本與浮標(biāo)所測結(jié)果一致。

第2塊數(shù)據(jù)中, 自相關(guān)函數(shù)方法得到的主波波長、波浪傳播方向、有效波高分別為128 m, 270°和 0.57 m, 交叉譜方法得到的這3個參數(shù)分別為37 m、170°、0.55 m。很顯然, 該結(jié)果中自相關(guān)函數(shù)法由于噪聲的影響無法測得正確的結(jié)果。

圖7 2014年6月12日自相關(guān)函數(shù)法所得海浪譜Fig. 7 Sea-wave spectrum obtained from the auto-correlation method on June 12, 2014

圖8 2014年6月12日互相關(guān)函數(shù)法所得海浪譜Fig. 8 Sea-wave spectrum obtained from the cross-correlation method on June 12, 2014

圖9 2014年6月21日自相關(guān)函數(shù)法所得海浪譜Fig. 9 Sea-wave spectrum obtained from the auto-correlation method on June 21, 2014

圖10 2014年6月21日互相關(guān)函數(shù)法所得海浪譜Fig. 10 Sea-wave spectrum obtained from the cross-correlation method on June 21, 2014

這兩種功率譜估計方法所測主波波長和有效波高與浮標(biāo)一致性都較好, 誤差在可接受范圍之內(nèi)。根據(jù)交叉譜方法中相位差的方法, 基本可以去除180°模糊。圖7顯示方位向30°時能量大于210°, 這應(yīng)該是由于在210°方位向時異常回波較多, 因此會有較多的回波被刪除掉, 導(dǎo)致結(jié)果顯示波浪傳播在30°方位向。圖8中顯示有幾圈180°模糊并沒有去除掉, 因此會在30°方位向時產(chǎn)生一個較小的峰值。交叉譜方法得到的有效波高較小, 這應(yīng)該與互相關(guān)函數(shù)功率譜估計的過程中用到了相鄰兩束回波進(jìn)行互相關(guān),相關(guān)性強(qiáng)度必定會比自相關(guān)函數(shù)減弱, 但這樣做的優(yōu)點(diǎn)是可以去除兩束回波中不相關(guān)的信息, 對去噪有一定幫助。

若都采用互相關(guān)函數(shù)法, 6月12日結(jié)果得到的海浪傳播方向10°的誤差比21日22°的誤差小, 但主波波長10 m的誤差大于21日2 m的誤差。12日和21日校飛時海面有效波高分別是0.47 m和0.6 m, 所以海況均很小, 若單獨(dú)處理出某幾圈的結(jié)果, 誤差會較大, 因此本文兩塊數(shù)據(jù)分別進(jìn)行了30圈數(shù)據(jù)的平均和40圈數(shù)據(jù)的平均。通過本次結(jié)果處理可知, 交叉譜方法得到的海浪譜在去除180°模糊方面有較大優(yōu)勢, 并且在去噪方面有一定優(yōu)勢, 因此交叉譜比自相關(guān)函數(shù)譜更有優(yōu)勢。由于此次校飛試驗(yàn)所在的海況環(huán)境較低, 對于更高海況和更大的飛行高度尚沒有試驗(yàn), 以及風(fēng)浪和涌浪混合時海浪反演的準(zhǔn)確性無法得到驗(yàn)證。

4 結(jié)論

本文分析了二維海浪譜的近實(shí)時探測在業(yè)務(wù)化和科研方面的需求, 介紹了未來中法星星載波譜儀探測海浪譜的幾何特性, 根據(jù)機(jī)載波譜儀校飛試驗(yàn)和機(jī)載波譜儀的信號形成機(jī)制, 總結(jié)了波譜儀反演海浪的流程, 并利用該流程處理了機(jī)載校飛的回波數(shù)據(jù), 通過自相關(guān)函數(shù)和互相關(guān)函數(shù)兩種功率譜估計方法, 反演了二維海浪譜。最后根據(jù)浮標(biāo)得到的二維海浪譜對機(jī)載波譜儀探測二維海浪譜的有效性進(jìn)行了驗(yàn)證, 并對比了兩種海浪譜計算方法。主要結(jié)論如下:

1) 工作在Ku 波段的機(jī)載海洋波譜儀是一種真實(shí)孔徑雷達(dá)系統(tǒng), 波譜儀波束基于長波對海面微尺度波的傾斜調(diào)制, 通過方位向360°掃描測量海浪譜,通過本文的驗(yàn)證, 這種方法探測二維海浪譜以及波向、波長和有效波高是有效的。

2) 根據(jù)波譜儀信號形成機(jī)制, 采用互相關(guān)函數(shù)進(jìn)行功率譜估計在獲取主波波長、波向和有效波高時, 能與浮標(biāo)所測結(jié)果較為一致且精度較高。交叉譜方法相對于自相關(guān)函數(shù)方法能去除海浪傳播的180°模糊, 因此交叉譜方法更有優(yōu)勢。

3) 互相關(guān)函數(shù)法得到的交叉譜, 在計算有效波高時會偏小, 這與互相關(guān)函數(shù)功率譜估計的方法有關(guān)。由于互相關(guān)法進(jìn)行時, 兩個回波信號必須有差別,因此相關(guān)強(qiáng)度會變?nèi)? 導(dǎo)致結(jié)果總能量偏小。海浪信息越規(guī)則, 這兩種方法結(jié)論就會越一致。

4) 在計算斜率方差時采用了5°~12°入射角范圍進(jìn)行公式擬合, 根據(jù)本文結(jié)論可知, 低海況時定標(biāo)與否對最后反演的二維海浪譜結(jié)果不會產(chǎn)生決定性影響, 但該方法對于更高海況適用性尚未可知。未來星載波譜儀只有靠多波束聯(lián)合才能實(shí)現(xiàn)。

致謝: 非常感謝北京遙測技術(shù)研究所提供的機(jī)載波譜儀數(shù)據(jù)和國家海洋局第一海洋研究所提供的浮標(biāo)數(shù)據(jù)。

[1] Cavaleri L. Wave modeling: Where to go in the future[J]. Bulletin of the American Meteorological Society, 2006, 87(2): 207-214.

[2] Gonzalez F I, Beal R C, Brown W E, et al. Seasat synthetic aperture radar: Ocean wave detection capabilities[J]. Science, 1979, 204(4400): 1418-1421.

[3] Walsh E J, Hancock III D W, Hines D E, et al. Directional wave spectra measured with the surface contour radar[J]. Journal of Physical Oceanography, 1985, 15(5): 566-592.

[4] Jackson F C, Walton W T, Peng C Y. A comparison of in situ and airborne radar observations of ocean wave directionality[J]. Journal of Geophysical Research: Oceans (1978–2012), 1985, 90(C1): 1005-1018.

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

[6] Hauser D, Caudal G C, Rijckenberg G J, et al. RESSAC: A new airborne FM/CW radar ocean wave spectrometer[J]. Geoscience and Remote Sensing, 1992, 30(5): 981-995.

[7] Danièle H, Podvin T, Dechambre M, et al. Polarimetric measurements over the sea-surface with the airborne STORM radar in the context of the geophysical validation of the ENVISAT ASAR[J/OL]. [2015-09-17]. http: //earth.esa.int/workshops/polinsar2003/participants/ hauser5/POLINSAR-VALPARESO.pdf, 2003.

[8] Plant W J, Keller W C, Hayes K. Simultaneous measurement of ocean winds and waves with an airborne coherent real aperture radar[J]. Journal of Atmospheric and Oceanic Technology, 2005, 22(7): 832-846.

[9] Pettersson H, Graber H C, Hauser D, et al. Directional wave measurements from three wave sensors during the FETCH experiment[J]. Journal of Geophysical Research:Oceans (1978–2012), 2003, 108(C3): 209.

[10] Hauser D, Soussi E, Thouvenot E, et al. SWIMSAT: A real-aperture radar to measure directional spectra of ocean waves from space-main characteristics and performance simulation[J]. Journal of Atmospheric and Oceanic Technology, 2001, 18(3): 421-437.

[11] Li Xiuzhong, He Yijun, Zhang Biao, et al. The construction of a three-dimensional antenna gain matrix and its impact on retrieving sea surface mean square slope based on aircraft wave spectrometer data[J]. Journal of Atmospheric and Oceanic Technology, 2016, 33(4): 847-856.

[12] Barrick D E. Rough surface scattering based on the specular point theory[J]. Antennas and Propagation, 1968, 16(4): 449-454.

[13] Valenzuela G R. Theories for the interaction of electromagnetic and oceanic waves-A review[J]. Boundary-Layer Meteorology, 1978, 13(1-4): 61-85.

[14] Jackson F C, Walton W T, Hines D E, et al. Sea surface mean square slope from Ku-band backscatter data[J]. Journal of Geophysical Research: Oceans (1978–2012), 1992, 97(C7): 11411-11427.

[15] Hesany V, Plant W J, Keller W C. The normalized radar cross section of the sea at 10 incidence[J]. Geoscience and Remote Sensing, 2000, 38(1): 64-72.

[16] Freilich M H, Vanhoff B A. The relationship between winds, surface roughness, radar backscatter at low incidence angles from TRMM Precipitation Radar measurements[J]. J Atmos Ocean Technol, 2003, 20(4): 549-562.

[17] Engen G, Johnsen H. SAR-ocean wave inversion using image cross spectra[J]. Geoscience and Remote Sensing, 1995, 33(4): 1047-1056.

Received:Dec. 24, 2015

Retrieval method of an ocean wave spectrum using an airborne spectrometer and performing the validation

LI Xiu-zhong1, HE Yi-jun1, MENG Jun-min2, ZHANG Zhen-hua3
(1. Nanjing University of Science Information and Technology, Nanjing 210044, China; 2. First Institute of Oceanography, State Oceanic Administration, Qingdao 266061, China; 3. Beijing Research Institute of Telemetry, Beijing 100076, China)

spectrometer; power spectrum estimation; two-dimensional ocean wave spectra; 180° ambiguity

The ocean wave retrieval method is designed on the basis of the signal formation principle of an ocean wave spectrometer. By using the spectra estimation methods via the auto-correlation and cross-correlation functions, the two-dimensional ocean wave spectrum is obtained. Finally, after comparing the spectrum received from the buoy and that retrieved from the spectrometer, the effectiveness of detecting a two-dimensional spectrum from an airborne spectrometer is evaluated. We observe that in the environment of flight, the results of the methods using the auto-correlation and cross-correlation functions for retrieving ocean wave spectrum are consistent with that obtained from the buoy. From the cross spectrum, the ambiguity of 180° is excluded, although the significant wave height is smaller than that from auto-correlation method. When the sea-slope variance is calculated, the radar backscattering coefficients of the incidence angles at 5°–12° are fitted. Therefore, calibration of the radar backscattering coefficients is not required. Moreover, the future spaceborne spectrometer will be able to attain calibration of the radar backscattering coefficients using multibeam joints.

P731.22; TP732.1

A

1000-3096(2016)12-0123-08

10.11759/hykx20151224004

(本文編輯: 李曉燕)

2015-12-24;

2016-03-19

國家重點(diǎn)研發(fā)計劃子課題(2016YFC1401005); 江蘇省自然科學(xué)基金項目(BK2011008); 國家自然科學(xué)基金(41476158); 江蘇省高等教育優(yōu)勢學(xué)科項目

[Foundation: Research and Development of the National Key Program Corpus, No.2016YFC1401005; Natural Science Foundation of Jiangsu Province, No.BK2011008; National Science Foundation of China, No.41476158; The Preponderant Discipline Project of High Education in Jiangsu Province]

李秀仲(1985-), 山東濟(jì)寧人, 博士研究生, 主要從事海洋微波遙感研究, E-mail: qdlixiuzhong@163.com; 何宜軍, 通信作者,教授, 博士生導(dǎo)師, E-mail: yjhe@nuist.edu.cn