溫靜涵, 何宜軍, 2, 陳忠彪

基于X波段海洋雷達(dá)的南海西北部?jī)?nèi)波參數(shù)特征研究

溫靜涵1, 何宜軍1, 2, 陳忠彪1

(1. 南京信息工程大學(xué) 海洋科學(xué)學(xué)院, 江蘇 南京 210044; 2. 國(guó)家衛(wèi)星海洋應(yīng)用中心 空間海洋遙感與應(yīng)用研究重點(diǎn)實(shí)驗(yàn)室, 北京 100081)

中國(guó)南海是內(nèi)波頻發(fā)海域, 衛(wèi)星遙感在內(nèi)波參數(shù)特征的統(tǒng)計(jì)分析中得到了廣泛應(yīng)用, 但是衛(wèi)星軌道的重復(fù)訪問(wèn)時(shí)間長(zhǎng), 不能連續(xù)觀測(cè)內(nèi)波參數(shù)的變化特征。X波段海洋雷達(dá)具有高時(shí)間和空間分辨率, 可以長(zhǎng)期連續(xù)觀測(cè)內(nèi)波的變化, 本文提出了利用X波段海洋雷達(dá)圖像提取內(nèi)波參數(shù)的方法, 并利用連續(xù)觀測(cè)的數(shù)據(jù)研究了各參數(shù)的分布特征。首先, 對(duì)雷達(dá)圖像進(jìn)行預(yù)處理, 包括平均處理及斜坡校正兩步; 利用二維快速傅里葉變換確定內(nèi)波的傳播方向, 進(jìn)而根據(jù)該方向上的徑向廓線確定內(nèi)波的相速度大小、波長(zhǎng)及周期。最后, 利用在南海石油平臺(tái)上觀測(cè)的X波段雷達(dá)圖像對(duì)內(nèi)波各參數(shù)進(jìn)行提取和統(tǒng)計(jì)分析, 結(jié)果表明, 研究區(qū)域內(nèi)波的傳播方向多為西北向及西向傳播, 相速度多為0.6～0.8 m/s, 波長(zhǎng)一般為400～600 m, 周期大部分不超過(guò)1000 s。

海洋內(nèi)波; 南海西北部; X波段海洋雷達(dá); 參數(shù)統(tǒng)計(jì)

中國(guó)南海是海洋內(nèi)波的多發(fā)地, 內(nèi)波在海洋內(nèi)部能量的傳遞、海洋生物的生長(zhǎng)、海上的生產(chǎn)活動(dòng)和海洋工程等方面都有重要影響。隨著遙感圖像分辨率的提高, 衛(wèi)星遙感技術(shù)可以用于獲取大范圍的海洋內(nèi)波圖像。目前用于海洋內(nèi)波觀測(cè)的傳感器主要有合成孔徑雷達(dá)(Synthetic Aperture Radar, SAR)、可見光光譜儀和X波段海洋雷達(dá)等。SAR具有全天時(shí)、全天候觀測(cè)成像等優(yōu)點(diǎn), 可以獲得高分辨率的海面圖像[1]。利用SAR進(jìn)行內(nèi)波觀測(cè)最早始于SAR圖像中內(nèi)波調(diào)制條紋的發(fā)現(xiàn)[2], Alpers于1985年針對(duì)SAR圖像提出了內(nèi)波的雷達(dá)成像理論為后人利用SAR觀測(cè)內(nèi)波打下理論基礎(chǔ)[3]。Li等[4]通過(guò)測(cè)量?jī)?nèi)波波包之間的距離確定內(nèi)波的群速度, 并結(jié)合兩層有限深度模型導(dǎo)出混合層深度。Zheng等[5]于2001年提出了曲線擬合法和peak–peak法來(lái)確定內(nèi)孤立波的半振幅寬度, 進(jìn)而利用SAR圖像計(jì)算內(nèi)孤立波的振幅。Shen等[6]通過(guò)收集分析SAR圖像, 統(tǒng)計(jì)了南海西北部海洋內(nèi)波的分布并提出混合層深度的反演方法。Ning等[7]通過(guò)收集44幅哨兵一號(hào)及高分3號(hào)的SAR圖像, 統(tǒng)計(jì)了馬六甲海峽處內(nèi)波的振幅及內(nèi)波波包傳播的群速度的分布范圍。Jackson[8]于2007年將中分辨率成像光譜儀(MODIS)圖像中的太陽(yáng)耀斑用于分析全球范圍內(nèi)的高頻非線性內(nèi)孤立波, 同時(shí)分析了全球各個(gè)海域的內(nèi)波發(fā)生頻率。Wang等[9]于2012年分析了2500多幅光學(xué)衛(wèi)星和SAR圖像后, 對(duì)整個(gè)南海的內(nèi)波傳播進(jìn)行了系統(tǒng)的研究, 針對(duì)南海的不同區(qū)域分析內(nèi)波的傳播機(jī)制。Sun等[10]利用MODIS圖像對(duì)安達(dá)曼海中內(nèi)波的時(shí)空分布進(jìn)行統(tǒng)計(jì)并推測(cè)該地內(nèi)波發(fā)生源。但是, 由于衛(wèi)星軌道的限制, 星載SAR和MODIS的重復(fù)訪問(wèn)周期長(zhǎng)、時(shí)間分辨率不足, 不便研究?jī)?nèi)波連續(xù)的發(fā)展變化過(guò)程, 導(dǎo)致內(nèi)波參數(shù)統(tǒng)計(jì)不完全; 光學(xué)傳感器容易受云、霧、觀測(cè)角度等影響。

X波段雷達(dá)具有高時(shí)間(幾秒至幾分鐘)和空間分辨率(幾米)的優(yōu)點(diǎn), 可以實(shí)時(shí)獲得海面連續(xù)觀測(cè)的圖像, 近年來(lái)被廣泛應(yīng)用于海浪和海流等參數(shù)的觀測(cè)。隨著X波段雷達(dá)技術(shù)的發(fā)展, 研究人員也開始探索利用X波段雷達(dá)反演海洋內(nèi)波的方法。盡管與SAR相比, X波段雷達(dá)的探測(cè)范圍較小, 但它可以對(duì)海面進(jìn)行實(shí)時(shí)、連續(xù)的觀測(cè), 在海洋內(nèi)孤立波的動(dòng)態(tài)信息研究中有其獨(dú)到的優(yōu)勢(shì)。1990年, Watson和Robinson使用X波段海洋雷達(dá)數(shù)據(jù)研究直布羅陀海峽內(nèi)波的演化過(guò)程[11]。2009年, Ramos等[12]利用X波段雷達(dá)的圖像, 提出了一種基于Radon變換的方法, 計(jì)算內(nèi)波的傳播方向、非線性速度、孤立波之間的距離和每個(gè)波包中孤立波的數(shù)量等參數(shù)。在此基礎(chǔ)上, Lu等[13]基于X波段雷達(dá)、利用Radon變換的方法計(jì)算出南海東北部?jī)?nèi)孤立波的傳播方向和傳播速度等參數(shù), 并應(yīng)用海洋的兩層有限深度理論計(jì)算混合層深度。Plant等[14]通過(guò)對(duì)2005年和2007年部署在中國(guó)南海的艦載雙極化相干X波段雷達(dá)的圖像進(jìn)行分析, 發(fā)現(xiàn)南海至少存在兩種不同類型的內(nèi)波, 并指出兩類內(nèi)波各自的特征。Lund等[15]提出了一種從雷達(dá)圖像中全自動(dòng)檢索內(nèi)波的方法, 并且針對(duì)振幅超過(guò)20 m的大振幅內(nèi)波進(jìn)行了單獨(dú)討論。2016年, Badiey等[16]通過(guò)分析由錨定溫度傳感器組成的陣列和X波段雷達(dá)所得到的數(shù)據(jù), 統(tǒng)計(jì)了一系列內(nèi)波參數(shù), 包括內(nèi)波的傳播速度、傳播方向和振幅等。

已有的研究主要是利用SAR及MODIS圖像進(jìn)行內(nèi)波時(shí)間及空間分布的統(tǒng)計(jì), 由于衛(wèi)星的重復(fù)訪問(wèn)周期較長(zhǎng), 使得統(tǒng)計(jì)時(shí)間不連續(xù), 不能體現(xiàn)內(nèi)波一段時(shí)間內(nèi)連續(xù)的變化特征。為解決傳統(tǒng)方法的不足, 本文提出一種基于X波段雷達(dá)圖像提取海洋內(nèi)波參數(shù)方法, 實(shí)驗(yàn)時(shí)間內(nèi)X波段雷達(dá)固定在石油平臺(tái)上, 基于X波段雷達(dá)高時(shí)間和空間分辨率、實(shí)時(shí)連續(xù)觀測(cè)的特點(diǎn), 對(duì)海洋內(nèi)波參數(shù)進(jìn)行定量提取, 并對(duì)一段時(shí)間內(nèi)南海西北部的內(nèi)波參數(shù)進(jìn)行統(tǒng)計(jì)分析。本文結(jié)構(gòu)如下。第二節(jié)介紹所用的實(shí)驗(yàn)數(shù)據(jù), 第三節(jié)介紹本文提出的利用X波段雷達(dá)圖像提取內(nèi)波參數(shù)及參數(shù)統(tǒng)計(jì)的方法, 第四節(jié)利用實(shí)驗(yàn)數(shù)據(jù)對(duì)方法進(jìn)行驗(yàn)證, 第五節(jié)對(duì)研究結(jié)果進(jìn)行討論與總結(jié)。

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

本實(shí)驗(yàn)開展于2017年4—12月, 地點(diǎn)位于中國(guó)南海西北部的石油平臺(tái)上(圖1)。觀測(cè)海區(qū)位于海南島西南部約100 km處, 該地水深約為90 m。觀測(cè)海區(qū)的地形變化為由東南向西北水深逐漸減小, 潮汐以不規(guī)則全日潮為主, 前人已證實(shí)該地常有內(nèi)波出現(xiàn)[17]。實(shí)驗(yàn)所用的X波段海洋雷達(dá)基于一套標(biāo)準(zhǔn)的船載導(dǎo)航雷達(dá)開發(fā), 工作頻率為9.4 GHz, 工作在掠入射角, 極化方式為HH極化, 加裝了40 MHz的數(shù)據(jù)采集卡和數(shù)據(jù)采集軟件, 用于將雷達(dá)接收的海面后向散射回波轉(zhuǎn)換為8 bit的灰度圖像。雷達(dá)的主要參數(shù)如表1所示, 其中, 雷達(dá)天線的旋轉(zhuǎn)速度為42 r/min, 故雷達(dá)的時(shí)間分辨率為1.43 s, 實(shí)驗(yàn)中一組雷達(dá)序列采集32至128幅雷達(dá)圖像, 采集相鄰兩個(gè)圖像序列的時(shí)間間隔為3～ 10 min。雷達(dá)的徑向距離分辨率為3.75 m, 方位分辨率為1°, 最大探測(cè)范圍約為6 km。圖2為本實(shí)驗(yàn)中內(nèi)波通過(guò)觀測(cè)海區(qū)時(shí)X波段雷達(dá)采集的一幅圖像, 其中明暗變化的粗條紋是內(nèi)波造成的海面粗糙度變化引起的雷達(dá)回波, 從圖中目視可以清楚地判斷出該組內(nèi)波是向西北方向傳播。

表1 X波段海洋雷達(dá)系統(tǒng)的參數(shù)

由于實(shí)驗(yàn)條件的限制, 本實(shí)驗(yàn)中沒有內(nèi)波的現(xiàn)場(chǎng)實(shí)測(cè)數(shù)據(jù)。目前, 衛(wèi)星遙感圖像已經(jīng)被廣泛用于海洋內(nèi)波的觀測(cè), 本文選取與X波段海洋雷達(dá)同步觀測(cè)的MODIS圖像來(lái)估計(jì)內(nèi)波到達(dá)觀測(cè)海區(qū)的時(shí)間。MODIS搭載于NASA的極軌衛(wèi)星Terra和Aqua上, 其中Terra衛(wèi)星每天10: 30前后(當(dāng)?shù)貢r(shí)間)過(guò)境、Aqua衛(wèi)星每天13: 30前后(當(dāng)?shù)貢r(shí)間)過(guò)境。MODIS圖像的刈幅為2 330 km, 空間分辨率為250 m, 可以大范圍觀測(cè)海洋內(nèi)波。由于光學(xué)衛(wèi)星圖像容易受云層遮擋、實(shí)驗(yàn)中雷達(dá)調(diào)試等原因, MODIS與X波段海洋雷達(dá)同步觀測(cè)的數(shù)據(jù)較少, 為了驗(yàn)證利用X波段海洋雷達(dá)觀測(cè)內(nèi)波方法的可靠性, 選取雷達(dá)觀測(cè)區(qū)域內(nèi)出現(xiàn)內(nèi)波的鄰近時(shí)間內(nèi)較清晰的MODIS圖像。例如, 圖3即為本文獲取到的與X波段海洋雷達(dá)圖像(圖2)匹配的MODIS衛(wèi)星圖像, 圖中紅色五角星所在位置即為實(shí)驗(yàn)中所用的X波段海洋雷達(dá)所處位置, 黑色虛線框中有清晰的內(nèi)波條紋, 根據(jù)內(nèi)波的傳播方向和速度可以估計(jì)該內(nèi)波達(dá)到雷達(dá)觀測(cè)海區(qū)的時(shí)間(圖2), 表明利用X波段海洋雷達(dá)可以有效觀測(cè)海洋內(nèi)波。

2 方法

本節(jié)介紹利用X波段海洋雷達(dá)圖像提取內(nèi)波參數(shù)的方法, 包括內(nèi)波的傳播方向、內(nèi)波的相速度、內(nèi)波的波長(zhǎng)和周期等參數(shù)。

2.1 數(shù)據(jù)預(yù)處理

2.1.1 平均處理

2.1.2 斜坡校正

由于雷達(dá)收到的電磁能量隨距離增大而衰減, X波段雷達(dá)的后向散射強(qiáng)度隨距離的增加急劇減小。為了有效提取內(nèi)波的變化特征, 需要去除雷達(dá)圖像的灰度值隨距離的變化趨勢(shì), 即對(duì)雷達(dá)圖像進(jìn)行斜坡校正[15]。

其中,與為待定系數(shù), 可以采用最小二乘法進(jìn)行擬合確定。

然后, 用平均的灰度值減去擬合值:

完成斜坡校正過(guò)程后, 再將雷達(dá)圖像的灰度值調(diào)整至0～255

2.2 內(nèi)波參數(shù)的提取

其中, Δ為前導(dǎo)波波峰在兩個(gè)徑向廓線中的距離差, Δ為相鄰兩組雷達(dá)圖像觀測(cè)的時(shí)間差。根據(jù)上述其中一個(gè)徑向廓線中相鄰波峰的位置差值, 可以獲得內(nèi)波波包中相鄰內(nèi)孤立波之間的距離, 本文中將該距離定義為內(nèi)波波長(zhǎng), 但是并不是所有內(nèi)波都以波包的形式存在, 研究區(qū)域中也存在許多內(nèi)孤立波[19], 對(duì)于內(nèi)孤立波則不能以上述方法計(jì)算內(nèi)波波長(zhǎng), 故本文利用以下公式計(jì)算內(nèi)波波長(zhǎng):

= 2, (7)

其中,為傳播方向徑向廓線上前導(dǎo)波點(diǎn)(最亮點(diǎn))與相鄰最暗點(diǎn)之間的距離。已知內(nèi)波相速度及內(nèi)波波長(zhǎng)后可以利用下式求得內(nèi)波周期:

2.3 內(nèi)波參數(shù)的統(tǒng)計(jì)

對(duì)一段時(shí)間內(nèi)的研究區(qū)域中的X波段海洋雷達(dá)圖像進(jìn)行篩選, 標(biāo)記每一組有內(nèi)波出現(xiàn)的雷達(dá)圖像, 利用上述方法提取內(nèi)波相速度、內(nèi)波傳播方向、內(nèi)波波長(zhǎng)及周期等參數(shù), 統(tǒng)計(jì)觀測(cè)時(shí)間內(nèi)各參數(shù)的分布規(guī)律。

由于X波段雷達(dá)連續(xù)觀測(cè)海面的變化, 而內(nèi)波通過(guò)觀測(cè)海區(qū)需要一定的時(shí)間, 為了避免同一個(gè)內(nèi)波被重復(fù)統(tǒng)計(jì), 本文根據(jù)以下原則統(tǒng)計(jì)內(nèi)波參數(shù):

(1) 相鄰3 h內(nèi)只統(tǒng)計(jì)一次內(nèi)波, 即, 將3 h內(nèi)雷達(dá)圖像中出現(xiàn)且傳播方向相近的內(nèi)波作為同一個(gè)內(nèi)波;

(2) 當(dāng)一個(gè)內(nèi)波在相鄰的多組雷達(dá)圖像內(nèi)連續(xù)出現(xiàn)時(shí), 將雷達(dá)回波最強(qiáng)的圖像用于統(tǒng)計(jì)內(nèi)波的參數(shù)。

圖4為本文方法的流程圖。

3 結(jié)果

本節(jié)首先結(jié)合衛(wèi)星圖像驗(yàn)證提出的方法的有效性, 然后根據(jù)X波段雷達(dá)長(zhǎng)期的觀測(cè)數(shù)據(jù)統(tǒng)計(jì)觀測(cè)海區(qū)中內(nèi)波的變化特征。

3.1 內(nèi)波參數(shù)提取

為了進(jìn)一步說(shuō)明提出的利用X波段雷達(dá)圖像提取內(nèi)波參數(shù)的方法, 并檢驗(yàn)方法的有效性, 選取圖2中X波段雷達(dá)觀測(cè)的內(nèi)波為例, 提取內(nèi)波參數(shù)并與MODIS觀測(cè)的內(nèi)波做比較。

首先對(duì)X波段雷達(dá)圖像進(jìn)行平均處理, 如圖5所示, 可看出平均處理可以一定程度上平滑掉圖像中海表面的風(fēng)浪和涌浪, 使得內(nèi)波條紋更容易被識(shí)別, 雷達(dá)圖像中內(nèi)波波包的傳播方向與MODIS圖像(圖3)中內(nèi)波的傳播方向一致。

接下來(lái)對(duì)雷達(dá)圖像進(jìn)行斜坡校正, 即去除雷達(dá)圖像的灰度值隨距離的變化趨勢(shì), 圖6展示的為該組圖像灰度值隨距離擬合結(jié)果的示意圖(式(1)), 該擬合結(jié)果所得R–square值(確定系數(shù))為0.95, 接近于1, 說(shuō)明所選的冪函數(shù)模型適用。

3.2 內(nèi)波參數(shù)統(tǒng)計(jì)

X波段雷達(dá)具有可以連續(xù)、長(zhǎng)期觀測(cè)的優(yōu)點(diǎn), 可以有效獲取內(nèi)波參數(shù)的分布特征。下面利用第2節(jié)的X波段雷達(dá)數(shù)據(jù)統(tǒng)計(jì)觀測(cè)海區(qū)中內(nèi)波的變化特征。實(shí)驗(yàn)期間各月份中內(nèi)波的出現(xiàn)次數(shù)如圖9所示, 可以看出2017年7月至9月內(nèi)波出現(xiàn)次數(shù)較多, 4—6月出現(xiàn)內(nèi)波的次數(shù)呈遞增趨勢(shì), 10—12月出現(xiàn)內(nèi)波的次數(shù)遞減, 均少于7—9月出現(xiàn)內(nèi)波的次數(shù), 即這3個(gè)月為研究區(qū)域內(nèi)(南海西北部)內(nèi)波頻發(fā)的月份, 因此本文主要統(tǒng)計(jì)7—9月中內(nèi)波各參數(shù)的分布規(guī)律。

1) 內(nèi)波的傳播方向

首先對(duì)內(nèi)波的傳播方向進(jìn)行統(tǒng)計(jì)。圖10展示了7月、8月、9月及7—9月內(nèi)波傳播方向的統(tǒng)計(jì)圖, 在7月和8月, 內(nèi)波傳播方向最多的是西北偏西向, 在9月內(nèi)波多為西向傳播; 圖10(d)展示的是7—9月方向統(tǒng)計(jì)集合, 可以看出研究時(shí)間內(nèi), 內(nèi)波多為西北偏西向、西向傳播, 也存在少數(shù)內(nèi)波向南向、東南向及西南向傳播, 這與前人利用遙感圖像統(tǒng)計(jì)結(jié)果一致[20]。

2) 內(nèi)波的相速度

利用式(6)計(jì)算每組內(nèi)波相速度, 統(tǒng)計(jì)7—9月各月內(nèi)波相速度大小的分布情況, 如圖11所示。圖11(a)可以看出, 7月研究區(qū)域中內(nèi)波相速度集中分布于0.2～ 0.8 m/s, 也有少部分分布在0～0.2 m/s及0.8～1.8 m/s; 8月內(nèi)波速度分布較為均勻, 集中分布在0.2～0.8 m/s, 少部分分布在0.8～1.6 m/s中, 極少部分內(nèi)波相速度小于0.2 m/s或超過(guò)1.8 m/s; 9月內(nèi)波速度大部分處于0.6～1 m/s, 觀察7—9月內(nèi)波相速度大小分布圖, 內(nèi)波相速度最多落在0.6～0.8 m/s的區(qū)間內(nèi), 這個(gè)結(jié)論與前人利用SAR圖像總結(jié)出的經(jīng)驗(yàn)公式結(jié)果相吻合[21]。

3) 內(nèi)波的波長(zhǎng)

內(nèi)波波長(zhǎng)的統(tǒng)計(jì)結(jié)果如表2所示, 大部分內(nèi)波波長(zhǎng)為200～800 m(出現(xiàn)次數(shù)為74次), 尤其在400～ 600 m內(nèi)居多(出現(xiàn)31次), 部分為800～1000 m(出現(xiàn)14次), 但也有小部分處于1 000 m以上及200 m以下(共出現(xiàn)6次), 不同月份的波長(zhǎng)分布略有差別。

表2 內(nèi)波波長(zhǎng)的分布

注: 表中的數(shù)值為內(nèi)波的出現(xiàn)次數(shù)。

4) 內(nèi)波的周期

利用式(8)計(jì)算內(nèi)波每組內(nèi)波的周期, 統(tǒng)計(jì)于表3。周期集中分布在0～2 000 s(出現(xiàn)次數(shù)為82次), 其中大部分內(nèi)波周期小于1 000 s(出現(xiàn)62次), 2 000～3 000 s內(nèi)也有分布(出現(xiàn)8次), 極少數(shù)周期大于3 000 s(共出現(xiàn)4次), 該結(jié)果在研究的3個(gè)月內(nèi)基本一致。

表3 內(nèi)波周期的分布

注: 表中的數(shù)值為內(nèi)波的出現(xiàn)次數(shù)。

此外, 前人的研究表明內(nèi)孤立波的周期(持續(xù)時(shí)間)一般在10～30 min(600～1 800 s)[22], 而超出這個(gè)范圍的值可能是涌浪或潮波。由于本文在提取內(nèi)波參數(shù)時(shí)對(duì)雷達(dá)圖像做了平均處理(第3.1節(jié)), 每次平均時(shí)所用的雷達(dá)圖像一般大于32幅, 而涌浪的周期一般小于20 s, 所以平均后的圖像中一般不包括涌浪; 潮波造成的海面粗糙度變化較小, 在雷達(dá)圖像不易觀察。本文分析了周期大于1 800 s的X波段雷達(dá)圖像, 結(jié)果表明圖像中的條紋的特征與內(nèi)波相同, 所以本文將3 000 s左右的條紋仍然認(rèn)為是內(nèi)波造成的。

4 討論與總結(jié)

南海是內(nèi)波頻發(fā)海域, 前人的研究大多集中于南海東北部, 且大多利用遙感影像對(duì)內(nèi)波參數(shù)進(jìn)行統(tǒng)計(jì)分析。由于衛(wèi)星軌道的重復(fù)訪問(wèn)時(shí)間長(zhǎng)不能連續(xù)觀測(cè), 統(tǒng)計(jì)時(shí)間跨度較大, 不能獲得連續(xù)時(shí)間內(nèi)波參數(shù)的分布特征。本文采用的X波段雷達(dá)具有高時(shí)間分辨率(幾秒至幾分鐘)和空間分辨率(幾米), 可以對(duì)實(shí)驗(yàn)地點(diǎn)進(jìn)行全天時(shí)的實(shí)時(shí)連續(xù)觀測(cè), 從而可以獲得觀測(cè)海區(qū)中內(nèi)波的大范圍、連續(xù)變化, 有利于全面分析內(nèi)波的變化特征?；赬波段海洋雷達(dá)圖像, 本文提出了內(nèi)波參數(shù)的提取方法, 并利用長(zhǎng)期的觀測(cè)數(shù)據(jù)統(tǒng)計(jì)了內(nèi)波參數(shù)的分布。觀測(cè)時(shí)間中內(nèi)波的主要分布特征如下:

(1) 南海西北部7—9月內(nèi)波出現(xiàn)次數(shù)最多, 5、6月次之, 11、12月最少;

(2) 南海西北部?jī)?nèi)波多為西北向及西向傳播, 也有少部分內(nèi)波為西南及東南向傳播;

(3) 內(nèi)波相速度大小最多集中于0.6～0.8 m/s之間, 0.4～0.6 m/s也存在多數(shù), 該海域內(nèi)波相速度最大值超過(guò)2 m/s;

(4) 內(nèi)波波長(zhǎng)多數(shù)為400～600 m, 其次為600～ 800 m、200～400 m, 也有極少數(shù)內(nèi)波波長(zhǎng)超過(guò)1 000 m;

(5) 內(nèi)波的周期絕大部分小于1 000 s, 少部分為1 000～2 000 s, 很少部分內(nèi)波的周期大于2 000 s, 最大的周期大于3 000 s。

由于實(shí)驗(yàn)中沒有內(nèi)波的實(shí)測(cè)數(shù)據(jù), 下一步將結(jié)合現(xiàn)場(chǎng)實(shí)測(cè)及數(shù)值模擬結(jié)果, 對(duì)觀測(cè)海區(qū)內(nèi)波的源地及成因進(jìn)行研究。

[1] PAN J, JAY D A, ORTON P M. Analyses of internal solitary waves generated at the Columbia River plume front using SAR imagery[J]. Journal of Geophysical Research, 2007, 112(C7): C07014.

[2] ELACHI C, APEL J R. Internal wave observations made with an airborne synthetic aperture imaging radar[J]. Geophysical Research Letters, 1976, 3(11): 647-650.

[3] ALPERS W. Theory of radar imaging of internal waves[J]. Nature, 1985, 314(6008): 245-247.

[4] LI X, CLEMENTE-COLON P, FRIEDMAN K S. Estimating oceanic mixed-layer depth from internal wave evolution observed from Radarsat-1 SAR[J]. Johns Hopkins Apl Technical Digest, 2000, 21(1): 130-135.

[5] ZHENG Q, YUAN Y, KLEMAS V, et al. Theoretical expression for an ocean internal soliton synthetic aperture radar image and determination of the soliton characteristic half width[J]. Journal of Geophysical Research Oceans, 2001, 106(C12): 31415-31423.

[6] SHEN H, HE Y. SAR imaging simulation of horizontal fully two-dimensional internal waves[J]. Journal of Hy-drodynamics, 2006, 18(3): 294-302.

[7] NING J, SUN L, CUI H, et al. Study on characteristics of internal solitary waves in the Malacca Strait based on Sentinel-1 and GF-3 satellite SAR data[J]. Acta Oceanologica Sinica, 2020, 39(5): 151-156.

[8] JACKSON C. Internal wave detection using the Moderate Resolution Imaging Spectroradiometer (MODIS)[J]. Journal of Geophysical Research, 2007, 112(C11): C11012.

[9] WANG J, HUANG W, YANG J, et al. Study of the propagation direction of the internal waves in the South China Sea using satellite images[J]. Acta Oceanologica Sinica, 2013, 32(5): 42-50.

[10] SUN L, ZHANG J, MENG J. A study of the spatial- temporal distribution and propagation characteristics of internal waves in the Andaman Sea using MODIS[J]. Acta Oceanologica Sinica, 2019, 38(7): 121-128.

[11] WATSON G, ROBINSON I S. A Study of internal wave propagation in the Strait of Gibraltar using shore-based marine radar images[J]. Journal of Physical Oceanography, 1990, 20(3): 374-395.

[12] RAMOS R J, LUND B, GRABER H C. Determination of internal wave properties from X-Band radar observations[J]. Ocean Engineering, 2009, 36(14): 1039-1047.

[13] LV H, HE Y, SHEN H, et al. A new method for the estimation of oceanic mixed-layer depth using shipboard X-band radar images[J]. Chinese Journal of Oceanology and Limnology, 2010, 28(5): 962-967.

[14] PLANT W J, KELLER W C, HAYES K, et al. Characteristics of Internal Waves in the South China Sea Observed by a Shipboard Coherent Radar[J]. IEEE Journal of Oceanic Engineering, 2011, 36(3): 441-446.

[15] LUND B, GRABER H C, XUE J, et al. Analysis of internal wave signatures in Marine Radar Data[J]. IEEE Transactions on Geoscience & Remote Sensing, 2013, 51(9): 4840-4852.

[16] BADIEY M, WAN L, LYNCH J F. Statistics of Nonlinear Internal Waves during the Shallow Water 2006 Experiment[J]. Journal of Atmospheric & Oceanic Technology, 2016, 33(4): 839-846.

[17] 孫麗娜, 張杰, 孟俊敏. 基于遙感與現(xiàn)場(chǎng)觀測(cè)數(shù)據(jù)的南海北部?jī)?nèi)波傳播速度[J]. 海洋與湖沼, 2018, 49(3): 471-480.

SUN Lina, ZHANG Jie, MENG Junmin. On propagation velocity of internal solitary waves in the northern South China Sea with remote sensing and in-situ observation data[J]. Oceanologia et Limnologia Sinica, 2018, 49(3): 471-480.

[18] CHEN Z, ZHANG B, KUDRYAVTSEV V, et al. Estimation of sea surface current from X-Band marine radar images by cross-spectrum analysis[J]. Remote Sensing, 2019, 11(9): 1031-1045.

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

FANG Xinhua, DU Tao. Fundamentals of oceanic internal waves and internal waves in the China Seas[M]. Qingdao: China Ocean University Press, 2005.

[20] 王雋. 基于衛(wèi)星遙感觀測(cè)的南海內(nèi)波發(fā)生源與傳播路徑分析[D]. 青島: 中國(guó)海洋大學(xué), 2012.

WANG Juan. The sources and propagation analysis of inter-nal waves in the South China Sea based on satellite remote sensing[D]. Qingdao: Ocean University of China, 2012.

[21] 湯琦璇. 基于SAR圖像的海洋內(nèi)波識(shí)別與參數(shù)提取技術(shù)研究[D]. 哈爾濱: 哈爾濱工程大學(xué), 2019.

Tang Qixuan. Research on oceanic internal wave detection and parameter extraction technology based on SAR image[D]. Harbin: Harbin Engineering University, 2019.

[22] RAMP S R, YANG Y, BAHR F L. Characterizing the nonlinear internal wave climate in the northeastern South China Sea[J]. Nonlinear Processes in Geophysics, 2010, 17(5): 481-498.

Study of the characteristics of internal wave parameters in the northwestern South China Sea based on X-band marine radar

WEN Jing-han1, HE Yi-jun1, 2, CHEN Zhong-biao1

(1. College of Marine Sciences, Nanjing University of Information Science & Technology, Nanjing 210044, China; 2. Key Laboratory of Space Ocean Remote Sensing and Applications, National Satellite Ocean Application Service, Beijing 100081, China)

Internal waves (IWs) frequently occur in the South China Sea (SCS), and satellite remote sensing images have been widely used to analyze the characteristics of IW parameters. However, the repeated access time of the satellite orbit is long; thus, the characteristics of IW parameters cannot be continuously observed. X-band marine radar has high temporal and spatial resolutions, which can be used to observe the change of IWs continuously. In this study, a method for determining IW parameters from X-band marine radar images is proposed, and the characteristics of each parameter are investigated using continuous observation data. Firstly, X-band marine radar images are preprocessed by averaging and ramp correction. Secondly, the propagation directions of IWs are determined by the two-dimensional fast Fourier transform, and the parameters of IWs are determined by the radial profile in this direction, such as the phase velocity, wavelength, and period. Finally, the X-band radar images observed on the SCS oil platform are used to extract and statistically analyze the parameters of IWs. The results show that the propagation directions of IWs in the study area are mostly northwest and west, the phase velocities are mostly 0.6–0.8 m/s, the wavelengths are generally 400–600 m, and the periods are mostly less than 1, 000 s.

internal wave; the northwestern South China Sea; X-band marine radar; parametric statistics

Oct. 9, 2021

P733

A

1000-3096(2022)09-0055-09

10.11759/hykx2021109002

2021-10-09;

2022-04-12

國(guó)家自然科學(xué)基金國(guó)際合作項(xiàng)目(41620104003); 國(guó)家自然科學(xué)基金重大儀器設(shè)備研制項(xiàng)目(42027805); 江蘇省海洋科技創(chuàng)新專項(xiàng)項(xiàng)目(HY2018-12); 國(guó)家自然科學(xué)基金青年基金(41506199); 江蘇省自然科學(xué)基金青年基金(BK20150905)

[The International Cooperation Project of the National Natural Science Foundation of China, No. 41620104003; The National Natural Science Foundation of China Major Instrument and Equipment Development Project, No. 42027805; Jiangsu marine science and technology innovation project, No. HY2018-12; The National Natural Science Youth Foundation of China, No. 41506199; The Natural Science Youth Foundation of Jiangsu Province, No. BK20150905]

溫靜涵(1997—), 女, 山東煙臺(tái)人, 碩士研究生, 主要從事海洋內(nèi)波研究, E-mail: jinghanwen@nuist.edu.cn; 陳忠彪(1987—),通信作者, 副研究員, 研究方向: 海洋微波遙感, E-mail: chenzhongbiao@nuist.edu.cn

(本文編輯: 康亦兼)