郭靖，杜濤

（1.中國海洋大學(xué)海洋環(huán)境學(xué)院，山東青島266100；2.中國人民解放軍91329部隊，山東威海264300）

利用SAR遙感衛(wèi)星資料研究東海內(nèi)孤立波時空分布特征

郭靖1,2，杜濤1

（1.中國海洋大學(xué)海洋環(huán)境學(xué)院，山東青島266100；2.中國人民解放軍91329部隊，山東威海264300）

通過分析歐洲空間局的地球資源衛(wèi)星（ERS-1/2）從1995—2010年的合成孔徑雷達（SAR）圖像，研究了東海內(nèi)孤立波的時間分布特征，繪制了東海內(nèi)孤立波的空間分布圖，并且反演了個別海區(qū)內(nèi)孤立波的傳播速度。研究結(jié)果表明：東海內(nèi)孤立波主要分布在兩個海域，臺灣東北部海域和長江口、杭州灣以及外舟山群島附近海域。在這兩個海域，大潮期間的強潮流有利于內(nèi)孤立波的出現(xiàn)。在長江口及附近海域，躍層強時，內(nèi)孤立波出現(xiàn)的概率相對較高。臺灣東北部內(nèi)孤立波波速約為0.6—0.8 m/s，長江口及附近海域內(nèi)孤立波波速約為0.6 m/s左右。

內(nèi)孤立波；東中國海；合成孔徑雷達

1 引言

通過合成孔徑雷達（SAR）觀測內(nèi)波始于1976年Elachi和Apel的研究[1]，而內(nèi)波的成像機理直到1985年才被Alpers揭示出來[2]。1986年Thompson和Gasparovic對內(nèi)波在L波段和X波段SAR圖像上的特征進行了仿真研究，并與該波段SAR圖像中的實際特征作了比較，驗證了Alpers的理論[3]。Brandt給出了一個描述直布羅陀海峽非線性內(nèi)波傳播的兩層數(shù)值模式，基于對150多幅內(nèi)波SAR影像的分析研究，表明所提出的模式可以較好的描述直布羅陀海峽內(nèi)波的時空演化[4]。DaSilva利用近實時的地球資源衛(wèi)星ERS-1的SAR圖像指導(dǎo)海洋內(nèi)波的調(diào)查研究[5]。

許多學(xué)者也利用SAR衛(wèi)星圖像對我國近海內(nèi)波方面作了多方面的研究。如，Liu利用ERS的SAR研究了南海內(nèi)孤立波的分布[6]。Hsu等利用衛(wèi)星圖像研究了東海部分海域的內(nèi)波特征，特別是對臺灣東北部復(fù)雜內(nèi)波生成機制進行了研究，并繪制了臺灣東北部的內(nèi)波[7]。方欣華、杜濤對黃海和整個南海海域內(nèi)孤立波的時空分布做了全面細致的分析研究[8]。

國內(nèi)外學(xué)者通過SAR對東海海域內(nèi)孤立波做了許多研究，但這些研究大多集中在臺灣東北部海域，缺乏對整個東海海域的內(nèi)孤立波分布的系統(tǒng)研究。本文通過分析從1995—2010年間的SAR圖像資料，對整個東海海域內(nèi)孤立波的時間分布特征、空間分布特征及個別海域內(nèi)孤立波的傳播速度進行了分析研究。

2 SAR及數(shù)據(jù)

SAR是一種主動式微波成像雷達。它通過適當(dāng)處理所獲得的海面后向散射測量信號，產(chǎn)生標(biāo)準(zhǔn)化后向散射截面的圖像。這些圖像的分辨率可達到幾米到幾十米量級，因而能夠詳細地顯示出海面粗糙度的空間細節(jié)變化。自20世紀70年代末開始，美國的Seasat、歐空局ERS-1/2和加拿大的Radarsat等衛(wèi)星上的合成孔徑雷達已經(jīng)獲取了大量的圖像，通過分析這些圖像中由內(nèi)波引起的海面粗糙度的變化可以獲得眾多海洋內(nèi)波的遙感觀測資料，從而為內(nèi)波研究提供豐富的信息。

圖1 內(nèi)波、表面波和SAR圖像灰度的關(guān)系

圖2 SAR影像觀測到的內(nèi)孤立波

SAR對內(nèi)波成像主要包含以下3個物理過程[9]：

（1）內(nèi)波在傳播過程中引起表層流場發(fā)生輻聚輻散，見圖1a；

（2）表層流場的輻聚輻散對海表微尺度波的調(diào)制，改變了海表的粗糙度，見圖1b；

（3）粗糙度的改變影響了雷達波后向散射的強度，使SAR圖像中的灰度值發(fā)生變化；粗糙海表面的后向散射強度大，在SAR圖像上呈現(xiàn)亮條紋，反之則呈現(xiàn)暗條紋，見圖1c。

當(dāng)衛(wèi)星SAR圖片的明暗條紋是由內(nèi)波的影響產(chǎn)生時，可以通過觀察這些明暗條紋的位置、條紋數(shù)量、間距、輻射方向以及相鄰兩條紋組的相對距離和發(fā)生時間等，研究內(nèi)波出現(xiàn)的空間位置、傳播方向和傳播速度等內(nèi)波要素。圖2是ERS-2衛(wèi)星的合成孔徑雷達拍攝到的臺灣東北部附近的圖像，從圖像中可以明顯地看到眾多明暗相間的條紋。利用這些條紋，可以提取內(nèi)波的波長和波向等要素信息，再結(jié)合其他海洋要素就能夠反演內(nèi)波波速等要素。

本文的ERS-1/2衛(wèi)星SAR圖像來自“中國科學(xué)院對地觀測與數(shù)字地球科學(xué)中心(http://cs.rsgs.ac. cn/cs_cn/cshome.asp)”。按需求選擇區(qū)域、衛(wèi)星、數(shù)據(jù)獲取時間、圖像云量，下載衛(wèi)星影像資料。本文選擇的是ERS-1/2衛(wèi)星，從1995—2010年衛(wèi)星圖像。衛(wèi)星圖片的像素為512×512，100 km×100 km。本文選取圖像覆蓋范圍為（21°—32°N，119°—125°E）

3 統(tǒng)計方法

3.1 區(qū)域劃分

由于東海海域東側(cè)一些區(qū)域存在大量資料缺失，且該區(qū)域能觀測到內(nèi)孤立波的圖像極少，所以將東側(cè)一些區(qū)域忽略。根據(jù)ERS-1/2衛(wèi)星的運行軌跡，將所研究區(qū)域劃成許多個大小約為100 km× 100 km的地理區(qū)塊，它們與SAR影像圖觀測的范圍一致。整個區(qū)域劃分后如棋盤一般（見圖3）。然后，分別以字母和數(shù)字為橫、縱坐標(biāo)將每個小區(qū)塊進行編號。例如，圖3中左下角的小區(qū)塊編號為C1。

3.2 統(tǒng)計方法

在劃分的每一個小區(qū)塊中，分別統(tǒng)計覆蓋該區(qū)塊的遙感圖像總數(shù)和含有內(nèi)孤立波的圖像數(shù)量。根據(jù)圖像中內(nèi)孤立波出現(xiàn)的位置、形狀，將該內(nèi)孤立波繪制到實際區(qū)域的相應(yīng)位置上，可以得到東海內(nèi)孤立波的總體分布圖（見圖4）。對于內(nèi)孤立波出現(xiàn)次數(shù)較多的區(qū)域，需要進一步統(tǒng)計內(nèi)孤立波出現(xiàn)次數(shù)的月分布和農(nóng)歷日分布情況。據(jù)此，可以了解內(nèi)孤立波在全年各月份中的出現(xiàn)情況是否一致，內(nèi)孤立波的出現(xiàn)和潮流的變化關(guān)系是否密切等。

當(dāng)內(nèi)孤立波是由潮流和地形作用產(chǎn)生或者SAR圖像中的內(nèi)孤立波是從內(nèi)潮波中演變產(chǎn)生時，參照楊勁松等的方法[10-11]，內(nèi)孤立波傳播速度可以這樣近似獲得：首先，從SAR圖像上沿著內(nèi)孤立波的弧線輻射方向、測量兩相鄰的獨立內(nèi)孤立波群的頭波之間的距離Λ，依據(jù)該內(nèi)孤立波群是由半日潮還是全日潮產(chǎn)生的確定這兩相鄰內(nèi)孤立波群之間的時間間隔（半日潮或全日潮周期）T，內(nèi)孤立波的傳播速度C可這樣計算：

圖3 東海內(nèi)孤立波區(qū)域劃分

圖4 東海西部內(nèi)孤立波空間分布圖

4 統(tǒng)計結(jié)果

4.1 空間分布

將所下載的衛(wèi)星SAR圖片按照劃分后的小區(qū)塊一一定位并編號。找出其中含有內(nèi)孤立波的圖片，依據(jù)其定位，將圖片中的內(nèi)孤立波按照位置和形狀繪制到代表實際海域的圖中，最終得到東海的內(nèi)孤立波的空間分布圖（見圖4）。依據(jù)整個東海內(nèi)孤立波的空間分布圖可知，內(nèi)孤立波分布較為集中的區(qū)域一個是臺灣東北部海域，另一個是長江口、杭州灣外以及舟山群島附近的海域。剩余的海域內(nèi)孤立波也有出現(xiàn)，但出現(xiàn)的數(shù)量相對較少。

進一步分析這兩個海域中內(nèi)孤立波在遙感圖片上的特征，不難發(fā)現(xiàn)出現(xiàn)在兩個海域中的內(nèi)孤立波各具特色。出現(xiàn)在臺灣東北部海域的內(nèi)孤立波具有如下特征：（1）每一組內(nèi)孤立波的波峰線都比較短；（2）各種傳播方向的內(nèi)孤立波相互交錯。

而出現(xiàn)在長江口、杭州灣外和舟山群島附近海域的內(nèi)孤立波則具有這樣的特征：（1）與一般陸架淺海區(qū)內(nèi)孤立波不同，這里既有向岸傳播的內(nèi)孤立波，也有離岸傳播的內(nèi)孤立波；（2）波峰線既有與等深線平行趨勢的，也有不具備這一特征的。

4.2 時間分布

4.2.1 月變化分布

為了解東海內(nèi)孤立波出現(xiàn)較多的兩個海域（長江口、杭州灣外舟山群島附近和臺灣東北部）中、內(nèi)孤立波在一年內(nèi)各月份分布情況的變化規(guī)律，需要分別計算這兩個海域內(nèi)孤立波出現(xiàn)次數(shù)和出現(xiàn)概率隨月份的分布情況。內(nèi)孤立波在某月份出現(xiàn)的概率=該月觀測到的內(nèi)孤立波數(shù)（含有內(nèi)孤立波的SAR圖像數(shù)）/該月SAR圖像總數(shù)。圖5和圖6分別是長江口附近海域內(nèi)孤立波出現(xiàn)次數(shù)和出現(xiàn)概率的月分布情況。在該海域，雖然內(nèi)孤立波在每個月出現(xiàn)的次數(shù)不同，但在一年中每個月都會出現(xiàn)；從概率上看，以3—8月內(nèi)孤立波的出現(xiàn)概率比較大，1月內(nèi)孤立波的出現(xiàn)概率最低。

在東海的陸架海域，季節(jié)性躍層在每年的3月份后進入成長期，至5月份已相當(dāng)強盛，9月份以后又開始減明，從11月到次年2月基本消失。所以，圖5和圖6所示的長江口、杭州灣外已及舟山群島附近海域的內(nèi)孤立波出現(xiàn)次數(shù)和出現(xiàn)概率隨月份變化是與東海陸架海域海洋躍層隨月份變化有著基本一致變化規(guī)律。躍層強的月份，內(nèi)孤立波出現(xiàn)的次數(shù)多或概率大?；蛘呖梢哉f，兩者的強弱變化規(guī)律是基本一致的。而在沒有季節(jié)躍層的月份仍觀測到有內(nèi)孤立波出現(xiàn)，說明這些內(nèi)孤立波的產(chǎn)生與季節(jié)性躍層關(guān)系不大，它們很可能產(chǎn)生于長江口附近的陸地淡水和海水的對流混合[8]。

圖5 長江口附近海域內(nèi)孤立波出現(xiàn)次數(shù)的月分布規(guī)律

圖6 長江口附近海域內(nèi)孤立波出現(xiàn)概率的月分布規(guī)律

圖7 臺灣東北部海域內(nèi)孤立波出現(xiàn)次數(shù)的月分布規(guī)律

圖8 臺灣東北部海域內(nèi)孤立波出現(xiàn)概率的月分布規(guī)律

圖7 和圖8是臺灣東北部海域內(nèi)孤立波出現(xiàn)次數(shù)和出現(xiàn)概率的月分布情況?？梢钥闯觯号_灣東北部海域的內(nèi)孤立波也是每個月都出現(xiàn)，但每個月出現(xiàn)的次數(shù)不同；內(nèi)孤立波出現(xiàn)概率的月分布情況比較復(fù)雜，6、7、9、12月出現(xiàn)概率較高，其它月份出現(xiàn)概率相對較低。這與東海陸架海域海洋躍層隨月份變化的規(guī)律明顯不完全一致。

比較兩個海域內(nèi)孤立波出現(xiàn)次數(shù)和出現(xiàn)概率月分布情況，共同點是：兩個海域的內(nèi)孤立波在每個月都會出現(xiàn)。不同點是：在長江口、杭州灣外已及舟山群島附近海域，內(nèi)孤立波出現(xiàn)概率大小的變化和東海陸架海域海洋躍層強弱的變化有著較為一致的密切關(guān)系。即，躍層強時內(nèi)孤立波的出現(xiàn)概率就大；躍層弱時內(nèi)孤立波的出現(xiàn)概率就小。在臺灣東北部海域，夏季躍層較強時，內(nèi)孤立波的出現(xiàn)概率也較大；但在其他季節(jié)，兩者的強弱變化和大小變化步調(diào)不是特別一致。

4.2.2 農(nóng)歷日變化分布

將東海內(nèi)孤立波SAR影像的時間數(shù)據(jù)由公歷日期轉(zhuǎn)換為農(nóng)歷日期，然后分別統(tǒng)計計算長江口、杭州灣外已及舟山群島附近海域和臺灣東北部海域內(nèi)孤立波出現(xiàn)次數(shù)和概率的農(nóng)歷日分布情況，繪制出相應(yīng)的日分布圖。（圖中日期均為農(nóng)歷日期，1至30分別代表初一至三十）

圖9和圖10分別是長江口、杭州灣以及外舟山群島附近海域從農(nóng)歷初一至三十內(nèi)孤立波的出現(xiàn)次數(shù)和出現(xiàn)概率的日分布情況。可以看出，在農(nóng)歷十六至十九和農(nóng)歷初二，內(nèi)孤立波出現(xiàn)概率有小的峰值，說明強的潮流對內(nèi)孤立波的出現(xiàn)可能會有促進作用，但是整體來看，內(nèi)孤立波出現(xiàn)概率的大小與潮流強弱的變化規(guī)律并不完全一致。因此，在長江口、杭州灣以及外舟山群島附近海域，當(dāng)有強潮流出現(xiàn)時內(nèi)孤立波可能在該海域出現(xiàn)，也可能不出現(xiàn)。這一點與呂宋海峽的內(nèi)孤立波情況[12-13]不一樣。在呂宋海峽，強潮流雖然也不是導(dǎo)致內(nèi)孤立波產(chǎn)生的唯一因素，但只要強潮流出現(xiàn)，內(nèi)孤立波就可以產(chǎn)生[14]。

圖9 長江口附近海域內(nèi)孤立波出現(xiàn)次數(shù)的農(nóng)歷日分布規(guī)律

圖10 長江口附近海域內(nèi)孤立波出現(xiàn)概率的農(nóng)歷日分布規(guī)律

圖11 臺灣東北部海域內(nèi)孤立波出現(xiàn)次數(shù)的農(nóng)歷日分布規(guī)律

圖12 臺灣東北部海域內(nèi)孤立波出現(xiàn)概率的農(nóng)歷日分布規(guī)律

圖11 和圖12分別是臺灣東北部海域內(nèi)孤立波出現(xiàn)次數(shù)和概率的農(nóng)歷日分布情況。內(nèi)孤立波出現(xiàn)概率在農(nóng)歷初一、初四和十四、十七的大潮期間呈現(xiàn)較高的峰值，表明大潮前后較強的潮流可能有利于內(nèi)孤立波的產(chǎn)生或出現(xiàn)。但就整個月份來看，這個作用并不是決定性的。因為大潮期間內(nèi)孤立波的出現(xiàn)概率并不總是處在高峰，且小潮期間，內(nèi)孤立波的出現(xiàn)概率也并非一定特別低。

限于每天觀測的SAR圖像的總數(shù)不一樣，并且觀測樣本數(shù)目較小，對內(nèi)孤立波出現(xiàn)概率的農(nóng)歷日分布規(guī)律的分析結(jié)果，存在誤差是必然的，這里的結(jié)果僅作為參考。

圖13 長江口附件和臺灣東北部海域內(nèi)孤立波波速統(tǒng)計（單位：m/s）

4.3 內(nèi)孤立波波速

依據(jù)本文3.2所述從遙感圖像中推算內(nèi)孤立波傳播速度的方法，選取一些典型的圖像，分別對長江口、杭州灣外以及舟山群島附近海域和臺灣東北部海域內(nèi)孤立波的傳播速度進行了計算，結(jié)果如圖13。在臺灣東北部海域，內(nèi)孤立波傳播波速約為0.6—0.8 m/s。在長江口、杭州灣外和舟山群島附近海域，內(nèi)孤立波傳播波速約為0.6 m/s左右。臺灣東北部內(nèi)孤立波波速計算結(jié)果與楊勁松等的結(jié)果一致[10]。

5 結(jié)論

利用從1995—2010年間的ERS-1和ERS-2 SAR衛(wèi)星數(shù)據(jù)，研究了東海內(nèi)孤立波的時、空分布，通過分析研究得到以下結(jié)論:

（1）從內(nèi)孤立波的空間分布特征發(fā)現(xiàn)，東海內(nèi)孤立波集中出現(xiàn)在兩個海域：臺灣東北部海域和長江口、杭州灣外以及舟山群島附近海域；

（2）東海內(nèi)孤立波的時間分布特征：在長江口、杭州灣外以及舟山群島附近海域，內(nèi)孤立波在全年各月都有可能發(fā)生，3—8月內(nèi)孤立波的出現(xiàn)概率較高，1月內(nèi)孤立波的出現(xiàn)概率最低；內(nèi)孤立波出現(xiàn)概率的高低與東海陸架海域躍層的出現(xiàn)及強弱變化規(guī)律有一致的地方；即，躍層強時，內(nèi)孤立波出現(xiàn)的概率變大；躍層弱時，內(nèi)孤立波出現(xiàn)的概率變小。同樣，農(nóng)歷大潮期間潮流較強時，內(nèi)孤立波的出現(xiàn)概率也可能會較高。在臺灣東北部海域，內(nèi)孤立波也是在全年各月都可能出現(xiàn)，在6、7、9、12月出現(xiàn)概率較高，其它月份出現(xiàn)概率相對較低；內(nèi)孤立波出現(xiàn)概率和東海陸架海域的躍層變化沒有密切關(guān)系；農(nóng)歷大潮期間潮流較強時，內(nèi)孤立波的出現(xiàn)概率也可能提高；

（3）內(nèi)孤立波的傳播速度：臺灣東北部海域內(nèi)孤立波傳播波速約為0.6—0.8 m/s，長江口、杭州灣外以及舟山群島附近海域內(nèi)孤立波傳播波速約為0.6 m/s；

（4）較強的潮流會增加內(nèi)孤立波在臺灣東北部海域和長江口、杭州灣外以及舟山群島附近海域的出現(xiàn)或產(chǎn)生的概率或可能性。

[1]Elachi,Apel J R.Internal wave observations made with an airborne synthetic aperture imaging radar[J].Geophys Res,1976,3: 647-650.

[2]Alpers.Theory of Radar Imaging of Internal Waves[J].Nature, 1985,314:245-247.

[3]Thompson,Gasparovic.Intensity modulation in SAR images of internal waves[J].Nature,1986,320(27):345-348.

[4]Brandt,Werner Alpers and Jan O.Backhaus.Study of the Generation and Propagation of Internal Wave in the Strait ofGibraltar Using a Numerical Mode and Synthetic Aperture Radar Images of the European ERS-1 Satellite[J].Geophys Res,1996, 101(C6):14237-14252.

[5]DaSilva J C,Jeans D R G,Robinson I S,et al.The application of near real-time ERS-1 SAR data to the prediction the location of the internal waves at sea[J].Remote Sensing,1997,18(16):3507-3517.

[6]Liu A K,Hsu M K.Internal wave study in the South China Sea using Synthetic Aperture Radar(SAR)[J].Remote Sensing,2004,4 (25):1261-1264.

[7]Hsu M K,Liu A K,Liu C.A Study of Internal Waves in the China Seas and Yellow Sea Using SAR[J].Continental Shelf Res,2000, 20:389-410.

[8]方欣華,杜濤.海洋內(nèi)波基礎(chǔ)和中國海內(nèi)波[M].青島:中國海洋大學(xué)出版社，2005:291-296.

[9]李海艷,杜濤.海洋內(nèi)波SAR影像成像理論[J].科學(xué)視野,2005, 29(10):85-89.

[10]楊勁松,周長寶,黃韋艮,等.合成孔徑雷達圖像內(nèi)波參數(shù)提取方法研究[J].遙感技術(shù)與應(yīng)用,2000,15(1):6-9.

[11]甘錫林,黃韋良,楊勁松,等.利用多源遙感衛(wèi)星數(shù)據(jù)研究南海內(nèi)波的時空分布特征[J].遙感技術(shù)與應(yīng)用,2007,22(2): 242-245.

[12]Xu Z H,Yin B S,Hou Y J,et al.A study of internal solitary waves observed on the continental shelf in the northwestern South China Sea[J].Acta Oceanol Sin,2010,29(3):18-25.

[13]Xu Z H,Yin B S,Hou Y J,et al.Variability of internal tides and near-inertial waves on the continental slope of the northwestern South China Sea[J].Geophys Res Oceans,2013,118,doi:10.1029/ 2012JC008212.

[14]Du T,Yu H T,et al.Impacts of tidal currents and Kuroshio intrusion on the generation of nonlinear internal waves in Luzon Strait[J].GeophysRes,2008,113,C08015,doi:10.1029/2007 JC004294.

Spatial and temporal characteristics of internal solitary waves in the East China Sea from SAR images

GUO Jing1,2，DU Tao1
（1.College of Physical and Environmental Oceanography,Ocean University of China,Qingdao 266003 China; 2.PLA Unit 91329,Weihai 264300 China）

By analyzing the SAR images of ERS-1/2 from 1995 to 2010,the spatial and temporal distribution of internal solitary waves in the whole East China Sea are studied.The speeds of internal solitary waves in two special areas are also calculated.The research shows that internal solitary waves in the East China Sea occur mainly in two areas:the northeast sea area of Taiwan and the sea area around Yangtze Estuary,Hangzhou Bay, and the Zhoushan Islands.Tidal current is not the only factor that can affect the occurrence of internal solitary waves in the two areas.When the pycnocline is strong in the sea area around Yangtze Estuary,Hangzhou Bay, and the Zhoushan Islands,internal waved occur more possibly.The calculated internal wave speed is about 0.6 m/s in the northeast sea area to Taiwan while it is about 0.6—0.8 m/s in the other area.

internal solitary waves；East China Sea；SAR

P731.24

：A

：1003-0239（2014）05-0001-07

10.11737/j.issn.1003-0239.2014.05.001

2013-10-30

國家“973”計劃（2007CB416605）

郭靖(1986-)，男,碩士研究生，主要從事海洋內(nèi)波研究。E-mail:ggguojinggg@gmail.com