星載SAR測量海洋流場研究進展

何宜軍 楊小波 矣娜 劉保昌

摘要多普勒信息是合成孔徑雷達（SAR）進行海洋動力遙感的一個重要參數，自從20世紀80年代末就已經開始利用該參數反演海表面動力過程的相關研究.本文重點論述了星載SAR多普勒信息反演海表面流及其應用研究，主要包括SAR干涉測量方法獲取海洋流場信息的主要技術手段及研究進展;多普勒頻移方法反演表層海洋流場的發展歷程及在上層海洋現象（包括地轉流、準地轉流以及中小尺度過程）中的應用研究;相關的國外重要數據融合計劃和研發計劃.目前國內在海洋流場的遙感方面還處于起步階段，重點關注于反演理論和方法，而將反演信息應用于上層海洋動力的研究，還十分缺乏.希望通過本文的綜述，為后續在該方向的研究提供借鑒.關鍵詞海表面流;微波遙感;多普勒中心;次級中尺度

中圖分類號P731.21;P715.7

文獻標志碼A

0引言

海氣界面是海洋與大氣之間交換物質和能量的通道，海表面的波浪場、風場、流場以及混合過程等都是海表面動力學主要的研究對象.海流是海水運動的重要方式之一，并且在海流運動過程中會對海洋上空的天氣和氣候產生重要的影響[1-2].例如，來自低緯度的海流向高緯度傳輸時，發生潛熱，海表的水汽向上輸送，就會帶來降水，因而需要增強對海洋-大氣耦合系統的了解和認識，這對海洋循環和氣候模型具有重要意義[2].海表面流也決定了海洋中生物和化學成分的輸送和分布，如海上溢油和其他污染物[3].除此之外，船舶航運路線的制定、海洋環境的監測和預報以及海上救援，都需要考慮海洋表面流的影響，因此需要對海面流場信息進行研究.

遙感方式獲取海表面的動力信息是目前一種相對高效且花費低廉的方法.微波遙感反演海表面動力場，除常用的輻射亮溫和后向散射系數外，最近30多年，多普勒信息也逐漸成為一種重要的海表動力遙感參數，尤其對海表流場和中尺度-次中尺度過程的研究，該參數具有十分重要的作用.有關海表動力場對于多普勒信息的影響以及信號仿真方面的綜述，已有專文論述[4]，本文將就星載SAR多普勒信息反演海表面流場以及研究海表面中尺度-次中尺度過程的相關文獻進行梳理和綜述，為后續研究提供借鑒.

在雷達多普勒反演海表流場方面，目前國際上主要有兩種遙感測量方式，其一是干涉測量法，最早由Goldstein等[5-6]和Shemer等[7]提出，后來以Romeiser等[8-15]為主，發展完善了順軌干涉測量（ATI）技術及M4S反演海表層流方法.此外還有英國南安普頓大學的Gommenginger等[16]提出的斜視SAR概念，馬來西亞Marghany等[17-20]提出的魯棒多普勒中心模型反演海流等.另一種是Chapron及其團隊等[21-34]先后研發出了RIM模型及其擴展模型、CDOP模型、DopRIM模型以及利用多普勒中心反演海表面流場的相關方法.

在多普勒反演海表層流的同時，還進一步研究和探討了海表存在的中尺度-次中尺度過程對海流反演的影響，進而從側面說明了多普勒信息對海表中尺度-次中尺度過程的成像能力，從而為利用多普勒信息（或多普勒與后向散射系數結合）研究中尺度-次中尺度過程提供了良好的分析基礎和借鑒案例.

除了星載SAR多普勒外，研究人員還提出了多普勒散射計反演海表層流的方法等[35-37]，這方面的內容已經超出本文的范圍，在此不再贅述.

1干涉SAR海流測量

合成孔徑雷達的干涉又稱為干涉合成孔徑雷達（InSAR）.干涉合成孔徑雷達（InSAR）是近幾十年來在射電天文學中發展起來的常規SAR技術和干涉測量技術的綜合[38].干涉雷達由于兩個接收天線與飛行軌跡的不同，又可以分為順軌干涉（ATI）和交軌干涉（XTI）.

順軌干涉法（ATI）的基本原理就是利用在順軌方向上間隔一定距離的兩個天線獲取同一場景的兩幅復圖像，這兩幅復圖像的相位差剛好反映了海流的徑向（沿雷達視線向）速度信息.順軌干涉法反演流場的處理流程包括SAR成像、配準、濾波、相位校正等幾個步驟，過程如圖1所示.

Romeiser等構建了新的海面后向散射模型[8-9]，其后正式稱為M4S模型[10].該模型與RIM模型都是在海面散射的雙尺度模型基礎上發展起來的模型，不同之處只在于對海面后向散射機制的理解和解釋.M4S模型的散射機制有鏡面散射和布拉格散射，并通過時間松弛法引入波流相互作用，來反映流場對后向散射的影響.Romeiser等[10]利用M4S模型，通過模擬結果與觀測結果之間的循環迭代逼近算法，模擬并反演了海面流場特征，反演的準確度在20%，空間分辨率50 m，主要影響因素為風和浪、水深調制、海洋動力現象的空間尺度以及反演所用的圖像數目多少，并指出高分辨率海流反演結果的驗證是最具有挑戰性的.

通過以上M4S模型的構建和反演流場的實例研究，Romeiser教授開始理解海面流場在海面后向散射機制中的作用，并著手建立流場信息與干涉相位信息之間的聯系.2000年，Romeiser等[11]在基于50 m分辨率海表面流對海面波譜存在水動力調制的假設下，全面分析了通過多普勒譜計算干涉相位，并反演得到海面流場的過程，主要的結果有：1）ATI適用于數平方千米內米級（50 m分辨率）的表面流觀測，其中還包含有海面長波的軌道速度分量;2）由于構造的多普勒譜的復雜性，在計算多普勒譜前，簡化后向散射模型很重要;3）影響模擬的主要因素是入射角和R/V比率，高頻的X波段雷達更適用于海面流場的測量，VV極化更適用于大尺度的后向散射場，高風速會增強流場梯度與反演的多普勒速度之間的線性相關，低風速響應更長的去相干時間從而使得ATI圖像相位差不易測到，所以ATI測流的理想風速是5～10 m/s，風向信息也是必須的，用于校正長波和布拉格波相速度產生的軌道運動對ATI速度的影響;4）高頻和大入射角對于流場反演有利，反演的流場對流場梯度的分辨率優于0.1 m/s.2001年，他們將ATI的反演結果與ADCP和環流模型結果對比，分析ATI反演流場的可行性及局限性，結果表明ATI能高分辨地反演海岸帶區域的流場信息，并可以利用反演的流場繼續反演高分辨率的水下地形，然而反演的流場信息在近岸或島礁附近效果極差[12].Kersten等[13]利用美國海軍研究實驗室的WINSAR-ATI系統獲得的福弗爾野外觀測站附近的河流區域的數據，證明了順軌干涉SAR相比于單向SAR得到的相位誤差更小，也就是反演得到的流速誤差更小.并且，對于系統相同的功率和天線總尺寸，順軌干涉法具有更高的靈敏度[14].2009年，Romeiser等[15]首次對TerraSAR-X干涉合成孔徑雷達（InSAR，ATI）數據進行了德國易北河河口的表面流場分析，在有效空間分辨率大于1 km時，能得到0.1 m/s的測流精度，適用于沿海區域的洋流觀測.

雖然順軌干涉方法對天線方向圖沒有嚴格的要求，與任何成像模式都是兼容的，并且能達到一定的測流精度，但是該方法也有一定的局限性，比如它只能得到一維的海洋表面流場的徑向流速信息，而對于流向沒有辦法得到.因此，為了獲取流向信息，還需要增加額外的自由度.Gommenginger等[16]的研究表明，可以利用測高儀的原理解決流向的反演，也就是聯合順軌和交軌干涉，但目前存在的問題是無法將順軌的干涉相位與交軌的干涉相位進行分離.2018年，日本科學家Ouchi等[39]在順軌干涉方法的基礎上，提出了多孔徑順軌干涉合成孔徑雷達系統（Multi-Aperture Along-Track Interferometric，MA-ATI），生成四幅子圖像，兩幅為前向子圖像，另外兩幅為后向子圖像，由此得到不同方向的干涉圖，最后由方位向速度分量求出海表面流矢量方向，對機載L波段的MA-ATI系統模擬得到的方向估計與理論值較為吻合.

在順軌干涉測量的基礎上，還擴展出了雙波束干涉的概念，即利用兩個斜視波束來獲得速度的兩個分量[40].該方法彌補了順軌干涉法只能測得雷達視線向海流流速這一不足.雙波束干涉最早由Rodriguez等[41]提出，該系統的成像幾何如圖2[42]所示.飛行平臺上搭載了一對干涉儀，分別前后放置，每個干涉儀都發送前向波束和后向波束.

2005年，Buck[43]針對寬刈幅海洋高度計（WSOA）設計的不足（由于極低的入射角使得導出的表面速度矢量的主要分量在垂直方向上）提出了Wavemill方案，該方案設置了兩根天線沿軌跡放置，各產生四個波束，進行飛行方向雙邊的前向、后向的順軌干涉測量，這些波束從45°斜視前后向對表面流場正交采樣，再結合順軌干涉方法，從而得到二維表面流場.2011年10月，Wavemill的機載實驗在英國西海岸愛爾蘭海的多個地點進行[44]，Martin等[44]利用2011年10月25日到26日的Wavemill數據估算了徑向表面速度，同時也驗證了風浪對表面速度多普勒測量的強烈影響，并發現了逆風和順風下速度的不對稱性.

2017年，德國宇航中心（DLR）的Wollstadt等[45]提出了Ku波段雙波束順軌干涉SAR，并應用了混合極化以提高不同海洋表面散射機制的信息水平，相比于Wavemill的設計，他們提出的雙波束模式設計為單邊的天線，能觀測到更大范圍的海洋場景，并且采用了較大的入射角以增加對水平方向速度的敏感性和提高風速的反演性能.

雙波束干涉模式的基本原理就是利用一對干涉儀發射兩組前向波束和后向波束，就等同于從兩個方向對同一海洋表面流速度進行觀測，該速度投影到兩個方向，同時利用順軌干涉的方法得到流速，結合雙波束的兩個方向就可以得到海洋表面流的流速矢量.首先利用順軌干涉的方法計算出海流流速Usurf，用公式可表示為[44]

Usurf=Δλe2πBesin θ， ?（1）

其中，Δ是干涉相位，λe

是雷達電磁波波長，Be

是有效基線，θ

是入射角.然后通過組合不同的干涉圖得到海表面流矢量估計.

Wollstadt等[40]在2015年提出的星載雙波束順軌InSAR是建立在Ku波段的單平臺系統，刈幅為200 km.2017年他們把該方法中的多個參數在不同波段的系統進行了優化，如有效基線這一參數，在30°的入射角和不同的風速下，順軌干涉的最優基線范圍為電磁波長的1 000倍到2 000倍之間[45].如果電磁波在X波段、C波段和Ku波段，那么與之相對應的最優基線分別是50、100和22～44 m.該系統的測速精度達到5 cm/s，其二級產品分辨率達到亞中尺度水平，為4 km×4 km，即使風速低至3 m/s也能完成同等精度的海流流速測量[45].該方法最大的特點就是可以完成海流二維速度矢量的測量.該方法采用兩種極化方式，發送采用的是圓極化，接受采用的是線極化，不考慮交叉極化，因此這也成為該方法的弊端，這是因為，在很強的風速下，交叉極化占據了信號的大部分，是不能忽略的.而Martin等[46]研究并測試了雙波束斜視干涉SAR系統，獲取的海表面流速的均方根誤差優于0.1 m/s.

我國學者對SAR模擬和反演海表層流的工作開展得相對較晚.2015年，何宜軍等[1]論述了反演海表層流的常用方法，常亮等[2]詳細論述了星載SAR反演海表層流的兩種方法，楊小波等[4]分析了微波遙感反演海流中的若干主要問題.在干涉SAR海流反演方面，任永政[47]利用ATI技術結合迭代算法以及STRM數據，反演了西歐海域流場，并與實測結果及Romeiser等[15]的結果進行了對比驗證，2017年又采用MCC方法對星座SAR圖像上海表層流進行了反演和驗證[48];此外，還有文獻[35-37]建立的多普勒散射計海表層流仿真模擬.在借鑒國外研究經驗的基礎上，文獻[49-50]也從基礎理論出發，嘗試構建新型干涉測量SAR系統（雙波束和雙通道雙波束干涉SAR）來模擬流場環境下的原始SAR數據，以及抑制方位向信號模糊造成的InSAR測流的影響.Pan等[51]基于機載圓形掃描SAR系統的實測數據，提取不同角度下的多普勒中心偏移，并考慮了雷達波束的方位指向和入射角引起的補償誤差，使得海表面流速估計結果更為接近實際海流流速.

2SAR多普勒中心偏移法

在研究海面多普勒信息之前，以法國國家海洋資源開發研究院（IFREMER）的Charpon等為核心的歐洲團隊利用后向散射系數已經逐步建立起了分析海流存在狀況下的海面后向散射模型[21-29].該團隊的Kudryavtsev[21]首先建立了綜合考慮布拉格散射和非布拉格散射（鏡面散射和波破碎）時的后向散射系數模型，并引入了風速的影響;2005年Kudryavtsev[22]建立了二維的雷達成像模型RIM，Johannessen[23]在波浪譜中加入了背景流場，利用RIM模型分析了背景流場對后向散射的影響.至此，該團隊已建立了以后向散射系數為中心，系統分析海表面流的研究方案.其后，該團隊又將RIM模型的理念進一步擴展，2012年發展了雙共極化SAR的海面成像模型，將極化差異和極化比引入后向散射系數模型中[24-25].2014—2016年又發展了全極化SAR的海面成像模型[26-28]，分別用于研究海面流特征、C波段和L波段共同成像的流特征以及浮油中的波破碎，另外還用RIM模型分析了存在潮流、徑流和水深變化情況下的長江口微波成像[29].然而值得注意的是，RIM模型引入背景流場的方式是通過波浪譜，但是西邊界流或入海徑流可以在多普勒圖像上觀測到，明顯不是因為流作用于局部的風浪譜，所以從RIM模型或者后向散射系數分析的海表面流場，只能是流與浪相互作用的部分，并不是全部的海表面流場信息.從RIM模型的分析可以知道，后向散射系數模型比較適宜分析中小尺度運動中背景流場的輻散輻合作用，對于大尺度的海表面流場研究能力有限，需要引入多普勒信息.

在多普勒信息的建模和分析方面，Charpon[30]首先從理論上分析了多普勒信號反演海流的可行性，分析了多普勒信號與后向散射系數之間的關系，給出了簡化的多普勒中心模型，并分析了可能造成多普勒信息偏差的因素.Johannessen等[31]在Charpon分析的基礎上建立了RIM模型與多普勒信息的關系，構建了DopRIM模型，在RIM模型的理論框架內引入了多普勒信息，用于分析多普勒或速度成像的海面狀況下風、浪、流的相互影響.Johannessen 等[32-33]將DopRIM模型用于分析南非南部海域阿古勒斯流的多普勒成像.Hansen等[34]用DopRIM模型模擬了強潮流區的波流相互作用，進一步說明了強表面流對海面粗糙度和斜距多普勒信號的調制作用，以及通過極化性質和入射角的差異，可以從DopRIM模型中分析和分離出非布拉格散射的部分.從DopRIM模型的相關論文可以看出，在構建該模型的過程中，Charpon及其團隊也在逐漸理解多普勒信息的成像機制、成像能力以及影響因素，雖然在此過程中，有些問題暫時沒有明晰，但是對于研究多普勒信息反演海表動力過程，起到了很大的推動作用.

在研制模型的同時，Charpon團隊也在嘗試著用多普勒信息反演大洋環流及其中尺度-次中尺度特征[52-62].Danilo等[52]首先比較了強潮流區HF雷達與SAR測量的徑向速度的差異，發現HF雷達與SAR分別響應不同尺度的波浪場，SAR反演結果主要受風場影響，源于風浪與表層流的相互作用，且風浪會限制流反演的精度.Collard等[53] 開始采用新的反演流的方法，利用神經網絡方法構建了CDOP模型計算風生多普勒信號，再從星載SAR測量的多普勒信號中減去風影響的部分，反演得到徑向流場.Rouault等[54]繼續采用移除風場多普勒信息的方法，詳細地反演了阿古勒斯流系.Backeberg等[55]進一步利用HYCOM和ROMS 模式，對比分析該流系的地形剪切作用和上層流強度.Hansen等[56]繼續采用從多普勒中心異常中移除風場影響的方法，反演了進入挪威海的大西洋水的兩條主要分支挪威大西洋坡流（NwASC）和挪威大西洋鋒面流（NwAFC），另外還有挪威沿岸流（NCC）.Hansen等[57]繼續采用校正風影響的方法，反演多普勒流，并依據地轉流時間尺度較長的特征，利用升降軌數據的時間平均分別反演地轉流的緯向和經向分量.Dagestad等[58]依據地轉流時間尺度較長的特點，利用時間平均的升降軌數據，反演得到了二維的多普勒流場，以及相應的緯向和經向地轉流分量，并對比了高度計反演的地轉流，在流的連續性和兩者的相關性方面均令人滿意，從而實現了多普勒信息反演二維地轉流的重要突破.Johannessen等[59]對海洋流場反演和應用進行了一次較全面的總結，認為對于次中尺度渦旋（直徑小于50 km）的觀測，SAR顯示出獨有的能力.Johannessen等[60]還反演了灣流區的徑向多普勒速度，并與GOCE測量的MDT反演的地轉流做比較，進一步驗證了多普勒在強流區的測量效果.Johannessen等[61]基于對阿古勒斯流的反演，總結認為將浮標數據、多普勒流場、高度計反演的地轉流和SST等光學數據進行融合，多數據綜合研究海面，不僅可以研究動力過程，還可以計算多種海洋現象的體積和熱量輸運，包括地轉流的強度和輸運量、渦度和輻散輻合場的非地轉和垂向運動及其輸運量.Husson等[62]利用多源數據，驗證了Sentinel-1衛星在海態反演（包括表面風場、方向波數譜、有效波高以及海面流場）方面的能力.

通過以上的研究可以看出，DOPRIM主要的作用是研究風、浪和流對平均的海面粗糙度的影響，以及造成的多普勒異常，但若用于反演海表面流則實際過程會非常復雜.Collard等[53]、Johannessen等[59-61]和Hansen等[62]，對Charpon等提出的反演方法進行簡化，發展CDOP模型，采用移除風多普勒影響的方法來反演海表層流，從結果及驗證結果來看，均較理想.而Hansen等[63]繼續對Charpon等的方法進行完善，提出了一個較全面的反演方法.從反演方法上看，該反演體系非常有助于理解多普勒中心，然而反演過程相當復雜，其中充分考慮了單個多普勒中心像元中后向散射系數的偏移，并依據陸地靜止不動的特征來校正衛星與海洋的相對運動，考慮因素非常全面，然而從反演的結果上來看效果不佳，甚至沒有給出反演的徑向流場圖像.分析其原因，除影響因素較多外，對多普勒中心的過度校正，也是很重要的一方面，而Collard等[53]和Johannessen等[59-61]提供的通過CDOP校正風影響的方法，非常符合實際的動力過程.Rouault等[54]分析入射角小于30°情況下多普勒速度的明顯異常，更充分說明了風與海面流場之間密切的互相關關系.基于這些研究，以及Mouche 等[64]發展的利用后向散射系數及CDOP模型聯合進行風場反演，結果都表明，采用風校正的方法反演海表層流，是目前最好的方法.正如文中結論所指出的，SAR成像海表面，目前實現了多普勒信息和高分辨率的后向散射系數的結合，使得風、浪和流的聯合反演成為可能，從而使得SAR遙感海表面動力具備了其他遙感手段無法比擬的優勢.

國內也開展了海表層流反演的相關工作，Fu 等[65]用SAR后向散射系數圖像，通過迭代算法，反演了臺灣海峽某淺灘的海表層流;Wang等[66]用Hansen等[63]提出的方法，反演了長江口附近的表面流場分布;Yu等[29]用RIM模型模擬和分析長江口的遙感成像;楊小波[67]用均值濾波的方法從多普勒中心中移除衛星運動造成的多普勒頻移，反演得到江浙沿岸流的流場分布.對比國內和國外的研究進展可以看出，我國在該研究領域目前尚處在起始發展階段，還需要投入大量的研究力量來發展和提高我國的研究水平.

3海流遙感的研發計劃

通過梳理文獻發現，國外在海洋流場的遙感方面先后提出和發展了幾個重要的研發和數據融合計劃，分別是GlobCurrent計劃、Wavemill計劃和SKIM計劃.這些計劃的實施，對于研究海洋流場的遙感具有極大的促進作用，在此簡要進行介紹.

GlobCurrent計劃啟動較早，最早開始于2012年的GlobCurrent論壇[59]，2014年進行了相關的概念解釋，主要的目的是利用衛星數據協同融合來提供海洋表面流場的定量估計，改善和增強衛星數據遙感上層海洋的應用，并展示以用戶為主導的科學應用、業務化應用和商業應用的重要影響[68].2015年，該計劃融合了歐空局Sentinel-1、Sentinel-3衛星與其他輔助支撐數據[68];2016年，分析了GlobCurrent數據的區域性特征[69];2017年在澳大利亞開展試驗，對GlobCurrent計劃的產品數據進行了充分的檢驗和驗證，并在2018年將Sentinel-2衛星數據也融入到了該數據產品中[70]，并對GlobCurrent數據的誤差校正進行了估算[71].GlobCurrent數據產品的開發和分布，對于反演和分析包括海表層流在內的上層海洋動力現象，提供了極大的便利.

Wavemill計劃采用順軌干涉測量來反演流場信息，最早的設計理念源于2001年的雙波束干涉測量海表面流，2004年發展成為雙邊順軌干涉測量海表流場，并于2005年正式提出Wavemill的概念.2011年，在英國中部的愛爾蘭海進行了機載驗證[16]，其后由英國南安普頓大學的Gommenginger及其研究團隊相繼開展了一系列研究，對風浪影響和風流聯合反演等工作進行了分析和仿真[44，72-75].該計劃是由歐空局資助，英國對地觀測儀器和空間技術中心（CEOI-ST）支持的研究項目，截至目前還沒有相關的衛星發射計劃發布.

SKIM是英文Sea Surface KInematics Multiscalee monitoring的縮寫，主要任務是發射雷達衛星觀測海面流場和海面波浪場，并由海面流場和波浪場擴展開來，綜合研究上層海洋現象和動力過程.SKIM計劃于2017年正式發布[76-77]，主要技術特征是結合了高度計和波譜儀的優點.其中雷達發射天線中，采用了Ka波段的8個波束進行干涉測量，其中6°和12°入射角的所有波數均為旋轉波數，與波譜儀采用旋轉多波束的特征類似，可以進行不同入射角的干涉測量，由此能夠提供精確的海面流場和浪場測量[77].SKIM主要的測量參數為多普勒中心和后向散射系數，根據此測量特征，歐空局稱其為多普勒海洋雷達，但遺憾的是未能列入未來衛星發射計劃.

4討論與總結

4.1存在的問題

星載SAR遙感海面流場方面，目前存在主要的問題主要有兩個：1）遙感反演的誤差分析模型不完善;2）反演結果的驗證方法不完善.在第一個問題方面，干涉測量方法目前已經發展了一些分析反演誤差的模型[35-37]，而多普勒中心偏移方法目前僅有一些對反演誤差的定性分析結果，還沒有發展出有效的誤差分析模型.對于第二個問題，Chapron團隊通過與大量的高度計及光學圖像的對比，對多普勒中心偏移法測量的徑向流場進行了大量的結果驗證，證明了SAR無論對于大尺度的地轉流還是中小尺度過程，均具有較好的成像能力，反演的結果真實準確.而干涉測量方法方面，結果的驗證則相對較少，且僅限于觀測數據驗證，不同遙感數據的交叉驗證以及對反演結果的分析尚缺乏完善的研究.

雖然經過了十幾年的研究，以上兩個問題依然存在，且未得到較好的解決，主要原因是SAR遙感海洋流場過程中，受到了衛星運動、復雜的海面散射以及復雜的風浪運動影響，使得解決這兩個問題比較困難.就誤差分析模型而言，反演結果的好壞受到諸多因素的影響，比如成像分辨率、海面的后向散射機制、反演的精度和反演對象的空間尺度等.這就使得在分析反演誤差時，需要針對不同尺度的海洋現象在不同分辨率下分別做具體的分析，測量對象的時空尺度決定測量誤差.因此每一種海洋現象都需要獨立分析其反演誤差，從而造成了反演結果驗證和誤差分析的復雜性.

然而，即便存在上述兩個不確定的問題，合成孔徑雷達依然是目前最理想的遙感手段.相比散射計、高度計等測量手段，SAR或InSAR遙感海洋流場信息將測量后向散射系數、多普勒頻率和干涉相位三個參數，它們分別具有不同的成像分辨率.對于時變海面的遙感而言，成像的時空分辨率決定了圖像所能探測的上層海洋現象時空尺度的范圍.后向散射系數圖像和干涉SAR的相位圖像分辨率較高，達到米級;多普勒頻率的分辨率較低，為千米級.SAR或InSAR具有的這種多分辨率優勢，將極大擴展衛星所能探測的上層海洋動力的時空尺度，這是其他遙感手段所不具備的，因而它是最適合的海面流場遙感觀測手段.

4.2研究展望

隨著遙感探測時空分辨率的逐步提高，海洋遙感所能探測的對象也日益豐富.星載SAR作為微波遙感海面動力的一個重要組成部分，可以預見在中尺度、次中尺度乃至小尺度海洋動力過程的研究中，將會逐漸起到主導性的作用.將不同時空尺度海洋動力過程的散射模型和反演方法建立起來，將會是未來很長一段時間內SAR動力遙感的核心研究方向.而隨著上層海洋觀測數據的日益完善，在傳統觀測和反演大尺度海洋環流和中尺度渦的基礎上，海洋遙感開始重點關注海表面存在的次中尺度和小尺度混合過程，為全面理解上層海洋的物質能量演變乃至氣候變化等重要課題奠定數據基礎.

目前，國內在相關的研究方向上尚處在初期發展階段，研究的成果還不是很豐富，研究的方向也相對單一，重點關注在遙感技術手段獲得反演流場信息方面.遙感信息在海洋學的應用研究，目前還十分缺乏，希望通過本文的討論，引起國內同行在該研究方向的關注和重視.

參考文獻

References

[1]何宜軍，劉保昌，張彪，等.海面流場衛星遙感方法綜述[J].廣西科學，2015，22（3）：294-300

HE Yijun，LIU Baochang，ZHANG Biao，et al.Overview on satellite remote-sensing methods for seasurface-current measurement[J].Guangxi Sciences，2015，22（3）：294-300

[2]常亮，高郭平，郭立新.星載SAR海洋表層流場反演綜述[J].海洋科學進展，2015，33（1）：107-117

CHANG Liang，GAO Guoping，GUO Lixin.Review on ocean surface current field measurement by space-borne SAR[J].Advances in Marine Science，2015，33（1）：107-117

[3]Klemas V.Remote sensing of coastal and ocean currents：an overview[J].Journal of Coastal Research，2012，28（3）：576-586

[4]楊小波，劉保昌，何宜軍，等.海表面流場微波遙感研究進展及其關鍵技術問題分析[J].廣西科學，2016，23（6）：499-506

YANG Xiaobo，LIU Baochang，HE Yijun，et al.Research progress in microwave remote sensing of seasurface currents and analysis of the key problems[J].Guangxi Sciences，2016，23（6）：499-506

[5]Goldstein R M，Zebker H A，Barnett T P.Remote sensing of ocean currents[J].Science，1989，246（4935）：1282-1285

[6]Goldstein R M，Zebker H A.Interferometric radar measurement of ocean surface currents[J].Nature，1987，328（6132）：707-709

[7]Shemer L，Marom M，Markman D.Estimates of currents in the nearshore ocean region using interferometric synthetic aperture radar[J].Journal of Geophysical Research：Oceans，1993，98（C4）：7001-7010

[8]Romeiser R，Alpers W，Wismann V.An improved composite surface model for the radar backscattering cross section of the ocean surface：1.theory of the model and optimization/validation by scatterometer data[J].Journal of Geophysical Research：Oceans，1997，102（C11）：25237-25250

[9]Romeiser R，Alpers W.An improved composite surface model for the radar backscattering cross section of the ocean surface：2.model response to surface roughness variations and the radar imaging of underwater bottom topography[J].Journal of Geophysical Research：Oceans，1997，102（C11）：25251-25267

[10]Romeiser R，Ufermann S，Alpers W.Remote sensing of oceanic current features by synthetic aperture radar：achievements and perspectives[J].Annals of Telecommunications，2001，56（11/12）：661-671

[11]Romeiser R，Thompson D R.Numerical study on the along-track interferometric radar imaging mechanism of oceanic surface currents[J].IEEE Transactions on Geoscience and Remote Sensing，2000，38（1）：446-458

[12]Romeiser R，Hirsch O.Possibilities and limitations of current measurements by airborne and spaceborne along-track interferometric SAR[C]∥IEEE International Geoscience & Remote Sensing Symposium，2001，DOI：10.1109/IGARSS.2001.976227

[13]Kersten P R，Toporkov J V，Ainsworth T L，et al.Estimating surface water speeds with a single-phase center SAR versus an along-track interferometric SAR[J].IEEE Transactions on Geoscience and Remote Sensing，2010，48（10）：3638-3646

[14]Frasier S J，Camps A J.Dual-beam interferometry for ocean surface current vector mapping[J].IEEE Transactions on Geoscience and Remote Sensing，2001，39（2）：401-414

[15]Romeiser R，Suchandt S，Runge H，et al.First analysis of TerraSAR-X along-track InSAR-derived current fields[J].IEEE Transactions on Geoscience and Remote Sensing，2009，48（2）：820-829

[16]Gommenginger C，Chapron B，Marquez J，et al.Wavemill：a new mission for high-resolution mapping of total ocean surface current vectors[C]∥EUSAR 2014：10th European Conference on Synthetic Aperture Radar，2014：1-4

[17]Marghany M，Hashim M，Cracknell A P.Hopfield neural network for sea surface current tracking from tiungsat-1 data[C]∥International Conference on Computational Science and Its Applications，2008：950-958

[18]Marghany M.Simulation sea surface current from RADARSAT-2 SAR data using Hopfield neural network[C]∥IEEE 5th Asia-Pacific Conference on Synthetic Aperture Radar （APSAR），2015，DOI：10.1109/APSAR.2015.7306326

[19]Marghany M.Hopfield neural network and Pareto optimal algorithms for retrieving sea surface current from TanDEM-X data[J].IOP Conference Series：Earth and Environmental Science，2018，169：012023

[20]Marghany M.Robust model for retrieval sea surface current from different RADARSAT-1 SAR mode data[C]∥IEEE International Conference on Signal and Image Processing Applications，2009，DOI：10.1109/ICSIPA.2009.5478705

[21]Kudryavtsev V.A semiempirical model of the normalized radar cross-section of the sea surface 1.background model[J].Journal of Geophysical Research Atmospheres，2003，108（C3）：8054

[22]Kudryavtsev V.On radar imaging of current features：1.model and comparison with observations[J].Journal of Geophysical Research Atmospheres，2005，110（C7）：C07016

[23]Johannessen J A.On radar imaging of current features：2.mesoscale eddy and current front detection[J].Journal of Geophysical Research Atmospheres，2005，110（C7）：C07017

[24]Kudryavtsev V N，Chapron B，Myasoedov A G，et al.On dual co-polarized SAR measurements of the ocean surface[J].IEEE Geoscience and Remote Sensing Letters，2013，10（4）：761-765

[25]Myasoedov A，Kudryavtsev V，Chapron B.Dual co-polarized SAR imaging of the ocean surface phenomena[C]∥Conference Proceedings of 2013 Asia-Pacific Conference on Synthetic Aperture Radar （APSAR），2014，DOI：10.1088/1742-6596/478/1/012001

[26]Kudryavtsev V，Kozlov I，Chapron B，et al.Quad-polarization SAR features of ocean currents[J].Journal of Geophysical Research：Oceans，2014，119（9）：6046-6065

[27]Kudryavtsev V，Kozlov I，Chapron B，et al.Quad-polarized SAR measurements of ocean currents in C-and L-bands[C]∥IEEE International Geoscience and Remote Sensing Symposium （IGARSS），2015，DOI：10.1109/IGARSS.2015.7326755

[28]Hansen M W，Kudryavtsev V，Chapron B，et al.Wave breaking in slicks：impacts on C-band quad-polarized SAR measurements[J].IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing，2016，9（11）：4929-4940

[29]Yu P，Johannessen J A，Kudryavtsev V，et al.Radar imaging of shallow water bathymetry：a case study in the Yangtze Estuary[J].Journal of Geophysical Research：Oceans，2016，DOI：10.1002/2016jc011973

[30]Chapron B.Direct measurements of ocean surface velocity from space：interpretation and validation[J].Journal of Geophysical Research Atmospheres，2005，110（C7）：C07008

[31]Johannessen J A，Kudryavtsev V，Chapron B，et al.Backscatter and Doppler signals of surface current in SAR images：a step towards inverse modelling[C]∥Advances in SAR Oceanography from Envisat and ERS Missions，2006

[32]Johannessen J A，Collard F，Chapron B，et al.Combined backscatter and Doppler shifts from SAR images：observations，synergies and simulations[C]∥SEASAR 2008-Advances in SAR Oceanography from ENVISAT and ERS Missions，2008

[33]Johannessen J A，Chapron B，Collard F，et al.Direct ocean surface velocity measurements from space：improved quantitative interpretation of Envisat ASAR observations[J].Geophysical Research Letters，2008，35（22）：L22608

[34]Hansen M W，Kudryavtsev V，Chapron B，et al.Simulation of radar backscatter and Doppler shifts of wave—current interaction in the presence of strong tidal current[J].Remote Sensing of Environment，2012，120：113-122

[35]鮑青柳.多普勒需達散射計的系統設計與仿真[D].北京：中國科學院空間科學與應用研究中心，2015

BAO Qingliu.System design and simulation of Doppler scatterometer：wide swath ocean surface current measurement[D].Beijing：Center for Space Science and Applied Research，2015

[36]鮑青柳，董曉龍，朱迪，等.基于筆形波束掃描雷達散射計的海洋表面流測量[J].電子學報，2015，43（6）：1200-1204

BAO Qingliu，DONG Xiaolong，ZHU Di，et al.Ocean surface current measurement using rotating pencil-beam scatterometer[J].Acta Electronica Sinica，2015，43（6）：1200-1204

[37]Miao Y J，Dong X L，Bao Q L，et al.Perspective of a Ku-Ka dual-frequency scatterometer for simultaneous wide-swath ocean surface wind and current measurement[J].Remote Sensing，2018，10（7）：1042

[38]Rosen P A，Hensley S，Joughin I R，et al.Synthetic aperture radar interferometry[J].Proceedings of the IEEE，2000，88（3）：333-382

[39]Ouchi K，Yoshida T，Yang C S.Multi-aperture along-track interferometric SAR for estimating velocity vector of ocean currents[C]∥IEEE International Geoscience and Remote Sensing Symposium，2018，DOI：10.1109/IGARSS.2018.8517425

[40]Wollstadt S，Lopez-Dekker P，de Zan F，et al.A Ku-band SAR mission concept for ocean surface current measurement using dual beam ATI and hybrid polarization[C]∥IEEE International Geoscience and Remote Sensing Symposium （IGARSS），2015

[41]Rodriguez E，Imel D，Houshmand B.Two-dimensional surface currents using vector along-track interferometry[C]∥Proceedings PIERS'95，1995：763

[42]Buck C，Donlon C，Gebert N.A status update of investigations into the Wavemill concept[C]∥EUSAR 2014：10th European Conference on Synthetic Aperture Radar，2014：1-4

[43]Buck C.An extension to the wide swath ocean altimeter concept[C]∥Proceedings of 2005 IEEE International Geoscience and Remote Sensing Symposium，2005，8：5436-5439

[44]Martin A，Gommenginger C，Chapron B，et al.Dual beam along-track interferometic SAR to MAP total ocean surface current vectors with the airborne wavemill proof-of-concept instrument：impact of wind-waves[C]∥IEEE International Geoscience and Remote Sensing Symposium （IGARSS），2015，DOI：10.1109/IGARSS.2015.7326719

[45]Wollstadt S，Lopez-Dekker P，de Zan F，et al.Design principles and considerations for spaceborne ATI SAR-based observations of ocean surface velocity vectors[J].IEEE Transactions on Geoscience and Remote Sensing，2017，55（8）：4500-4519

[46]Martin A C H，Gommenginger C P，Quilfen Y.Simultaneous ocean surface current and wind vectors retrieval with squinted SAR interferometry：geophysical inversion and performance assessment[J].Remote Sensing of Environment，2018，216：798-808

[47]任永政.從衛星TerraSAR-X圖像反演海面風場和海表流場方法研究[D].青島：中國海洋大學，2009

REN Yongzheng.Study on retrieval algorithms of sea surface wind fields and sea surface current fields from TerraSAR-X images[D].Qingdao：Ocean University of China，2009

[48]Ren Y Z，Li X M，Gao G P，et al.Derivation of sea surface tidal current from spaceborne SAR constellation data[J].IEEE Transactions on Geoscience and Remote Sensing，2017，55（6）：3236-3247

[49]Liu B C，He Y J，Li Y K，et al.A new azimuth ambiguity suppression algorithm for surface current measurement in coastal waters and rivers with along-track InSAR[J].IEEE Transactions on Geoscience and Remote Sensing，2019，57（6）：3148-3165

[50]端和陽，劉保昌，何宜軍，等.一種改進的基于逆Omega-K算法的海面場景SAR原始數據仿真方法[J].海洋與湖沼，2019，50（2）：291-300

DUAN Heyang，LIU Baochang，HE Yijun，et al.An improved raw data simulation method for measuring ocean surface current[J].Oceanologia et Limnologia Sinica，2019，50（2）：291-300

[51]Pan X L，Liao G S，Yang Z W，et al.Sea surface current estimation using airborne circular scanning SAR with a medium grazing angle[J].Remote Sensing，2018，10（2）：178

[52]Danilo C，Chapron B，Mouche A，et al.Comparisons between HF radar and SAR current measurements in the Iroise Sea[C]∥OCEANS 2007，2007，DOI：10.1109/OCEANSE.2007.4302424

[53]Collard F，Mouche A，Chapron B，et al.Routine high resolution observation of selected major surface currents from space[C]∥Proceedings of SEASAR，2008：21-25

[54]Rouault M J，Mouche A，Collard F，et al.Mapping the Agulhas current from space：an assessment of ASAR surface current velocities[J].Journal of Geophysical Research Atmospheres，2010，115（C10）：C10026

[55]Backeberg B C，Johannessen J，Rouault M，et al.Preliminary inter-model comparison of the Agulhas current with direct range Doppler velocity estimates from Envisat's Advanced Synthetic Aperture Radar （ASAR）[C]∥ESA Living Planet Symposium，2010

[56]Hansen M W，Johannessen J A，Dagestad K F，et al.Monitoring the surface inflow of Atlantic Water to the Norwegian Sea using Envisat ASAR[J].Journal of Geophysical Research Atmospheres，2011，116（C12）：C12008

[57]Hansen M W，Johannessen J A，Raj R.Mapping the Nordic Seas surface velocity using Envisat ASAR[C]∥Pan Ocean Remote Sensing Conference （PORSEC），2013，709：13

[58]Dagestad K F，Hansen M W，Johannessen J A，et al.INCUSAR：a method to retrieve temporal averages of 2D ocean surface currents from synthetic aperture radar Doppler shift[M]∥Fernndez-Prieto D，Sabia R.Remote Sensing Advances for Earth System Science.Berlin，Heidelberg：Springer Berlin Heidelberg，2013：59-67

[59]Johannessen J A，Holt B，Chapron B，et al.Ocean current retrievals and applications[C]∥SeaSAR 2012，2013，DOI：10.1002/（SICI）1521-3773（19981204）37：223.0.CO;2-#

[60]Johannessen J A，Chapron B，Collard F，et al.Range Doppler centroid shift-altimetry comparison[C]∥EUSAR 2014：10th European Conference on Synthetic Aperture Radar，2014

[61]Johannessen J A，Chapron B，Collard F，et al.Use of SAR data to monitor the Greater Agulhas current[M]∥Barale V，Gade M.Remote Sensing of the African Seas.Dordrecht：Springer Netherlands，2014：251-262

[62]Husson R，Mouche A，Chapron B，et al.Taking advantage of Sentinel-1 acquisition modes to improve ocean sea state retrieval[C]∥IEEE International Geoscience and Remote Sensing Symposium，2016，DOI：10.1109/IGARSS.2016.7730009

[63]Hansen M W，Collard F，Dagestad K F，et al.Retrieval of sea surface range velocities from envisat ASAR Doppler centroid measurements[J].IEEE Transactions on Geoscience and Remote Sensing，2011，49（10）：3582-3592

[64]Mouche A A，Collard F，Chapron B，et al.On the use of Doppler shift for sea surface wind retrieval from SAR[J].IEEE Transactions on Geoscience and Remote Sensing，2012，50（7）：2901-2909

[65]Fu B，Huang W G，Fan K G，et al.Method for ocean surface currents measurement by SAR[J].Proceedings of SPIE-the International Society for Optical Engineering，2008，DOI：10.1117/12.804897

[66]Wang L H，Zhou Y X，Ge J Z，et al.Mapping sea surface velocities in the Changjiang coastal zone with advanced synthetic aperture radar[J].Acta Oceanologica Sinica，2014，33（11）：141-149

[67]楊小波.基于ASAR的時變海表面流場反演研究[D].上海：上海海洋大學，2016

YANG Xiaobo.Sea surface current retrieval based on ASAR data[D].Shanghai：Shanghai Ocean University，2016

[68]Johannessen J A，Chapron B，Collard F，et al.GlobCurrent：sentinel-3 synergy in action[C]∥Proceedings of Sentinel-3 for Science Workshop，2015：SP-734

[69]Danielson R E，Korosov A，Johannessen J A，et al.A regional characterization of the GlobCurrent ocean surface current analysis[C]∥Proceedings of Living Planet Symposium，2016

[70]Johannessen J A，Chapron B，Collard F，et al.Advances in surface current observations from space：the GlobCurrent case[C]∥IEEE International Geoscience and Remote Sensing Symposium，2018，DOI：10.1109/IGARSS.2018.8519451

[71]Danielson R E，Johannessen J A，Quartly G D，et al.Exploitation of error correlation in a large analysis validation：GlobCurrent case study[J].Remote Sensing of Environment，2018，217：476-490

[72]Doody S，Marquez J，Dobke B，et al.A mission for measuring ocean surface current vectors[C]∥Proceedings of EUSAR 2016：11th European Conference on Synthetic Aperture Radar，2016：1-4

[73]Martin A C H，Gommenginger C，Marquez J，et al.Wind-wave-induced velocity in ATI SAR ocean surface currents：first experimental evidence from an airborne campaign[J].Journal of Geophysical Research：Oceans，2016，121（3）：1640-1653

[74]Martin A C H，Gommenginger C.Towards wide-swath high-resolution mapping of total ocean surface current vectors from space：airborne proof-of-concept and validation[J].Remote Sensing of Environment，2017，197：58-71

[75]Martin A C H，Gommenginger C P，Quilfen Y.Simultaneous ocean surface current and wind vectors retrieval with squinted SAR interferometry：geophysical inversion and performance assessment[J].Remote Sensing of Environment，2018，216：798-808

[76]Nouguier F，Chapron B，Collard F，et al.Sea surface kinematics from near-nadir radar measurements[J].IEEE Transactions on Geoscience and Remote Sensing，2018，56（10）：6169-6179

[77]Ardhuin F，Aksenov Y，Benetazzo A，et al.Measuring currents，ice drift，and waves from space：the Sea surface KInematics Multiscale monitoring （SKIM） concept[J].Ocean Science，2018，14（3）：337-354

Progress in sea surface current retrieval from

spaceborne SAR measurements

HE Yijun1YANG Xiaobo1YI Na1LIU Baochang1

1School of Marine Sciences，Nanjing University of Information Science & Technology，Nanjing210044

AbstractDoppler signal is an important synthetic aperture radar （SAR） parameter for ocean dynamics remote sensing.It was first used for the retrieval of ocean surface dynamic processes in the late 1980s.In this paper，we review the researches on retrievals of sea surface dynamics from spaceborne SAR Doppler measurements including：（1） SAR interferometry techniques and related research progress for ocean surface currents retrieval;（2） SAR Doppler shift and its application in retrieval of upper ocean dynamics （i.e.geostrophic currents，quasi-geostrophic currents，and meso-?and small-scale processes）;（3） Relevant data fusion and important R&D programs around the world.However，China's remote sensing of sea currents is still in the initial stage，and highly focused on retrieval theory and techniques.The application of the retrieved information in study of upper ocean dynamics is still limited，so this review can serve as a reference for future studies in this direction.

Key wordssea surface currents;microwave remote sensing;Doppler centroid;submeso-scale

收稿日期2019-10-19

資助項目國家重點研發計劃項目（2016YFC1401002）;國家自然科學基金國際合作重點項目（41620104003）;國家自然科學基金青年基金（41606201）

作者簡介何宜軍，男，博士，教授，博士生導師，長期從事海洋微波遙感機理、海上臺風遙感、海流遙感技術以及海氣相互作用等工作.yjhe@nuist.edu.cn

1南京信息工程大學海洋科學學院，南京，210044