袁凱瑞，王久娟，林佳瑜，韓 宇，王凱敏，

徐宇程3，陳勁松1,2*

(1.中國科學(xué)院深圳先進技術(shù)研究院空間信息研究中心，廣東 深圳 518055；2.深圳海洋環(huán)境大數(shù)據(jù)應(yīng)用與分析工程實驗室，廣東 深圳 518055； 3.深圳市海洋監(jiān)測預(yù)報中心，廣東 深圳 518055)

大亞灣位于大鵬半島的東部，是南海北部向陸地延伸最深的海灣。大亞灣地理區(qū)位優(yōu)勢突出，交通便利快捷，自然條件優(yōu)越。王友紹等(2004)認(rèn)為近年來大亞灣內(nèi)海岸線的改變、濕地不斷減少以及大量陸源污染物的輸入等，特別是石化區(qū)及核電站的建成，使得大亞灣海域的生態(tài)環(huán)境發(fā)生了較大的變化，海洋生態(tài)環(huán)境惡化趨勢明顯[1]。大亞灣海域水文動力的研究可為該海域的赤潮擴散、水交換能力估算、物質(zhì)運輸、環(huán)境保護等研究提供重要的數(shù)據(jù)和理論支撐。

目前對大亞灣水動力的研究主要集中在潮流與余流的特征方面。吳巖等(1998)通過潮流數(shù)值模擬發(fā)現(xiàn)，在大亞灣內(nèi)存在著兩個順時針方向的環(huán)流場，且南部環(huán)流強度大于北部環(huán)流[2]。楊國標(biāo)(2001)研究表明，大亞灣的潮流受天文潮和外海潮波的共同作用，灣內(nèi)主要為水平往復(fù)流，灣口及西部的余流較強，灣口余流流向灣外，其余基本沿岸線流動[3]。吳仁豪等(2007)利用HAMSOM模式對大亞灣的潮汐、潮流和余流進行了數(shù)值模擬，研究發(fā)現(xiàn)灣內(nèi)海域主要為南北向的往復(fù)流，流速受地形影響強烈，且灣口西側(cè)流速大于東側(cè)，夏季灣內(nèi)余流較小，灣口余流大體流向灣外，冬季余流較大，灣口余流大體流向灣內(nèi)[4]。王聰?shù)?2008)采用ECOM-si三維海洋數(shù)值模式模擬了大亞灣海流，認(rèn)為大亞灣潮汐和潮流特性主要受來自西太平洋的潮波制約即受廣東沿岸流控制，在灣口附近海域，余流從大辣甲島東部進入，西部流出[5]。Song等(2016) 、武文等(2017)、嚴(yán)聿晗等(2017)借助FVCOM海洋模式建立高時空分辨率的大亞灣三維潮汐潮流數(shù)值模型，分析了大亞灣內(nèi)潮位和潮流“雙峰”現(xiàn)象的原因，發(fā)現(xiàn)了灣內(nèi)不一致的大小潮變化以及灣外所受沿岸流的影響，揭示了大亞灣灣內(nèi)的余流多渦旋結(jié)構(gòu)[6-8]。以上研究多基于現(xiàn)場觀測資料或數(shù)值模擬結(jié)合現(xiàn)場觀測對大亞灣的潮流與余流進行分析。

本研究在大亞灣附近海域基于高頻地波雷達的長期海流觀測資料，利用Pawlowicz等(2002)發(fā)布的T_tide工具包[9]進行潮流調(diào)和分析，實現(xiàn)潮流和余流的分離，進而研究該海域余流的季節(jié)分布特征及變化規(guī)律，以期能夠為該海域的海上人員搜救、溢油擴散、污染物漂移等研究工作提供參考。

1 數(shù)據(jù)

1.1 數(shù)據(jù)簡介

1.1.1 高頻地波雷達數(shù)據(jù) 高頻地波雷達是利用雷達發(fā)射的電磁波與海浪作用時產(chǎn)生的Bragg散射和多普勒效應(yīng)，探測海表層流場、浪場和風(fēng)場的空間分布，同其他的海洋監(jiān)測設(shè)備相比，它具有監(jiān)測范圍大、探測精度高、受環(huán)境的影響較小、建設(shè)經(jīng)費少以及能夠進行連續(xù)無人值守的實時在線監(jiān)測等優(yōu)點，目前已成為了一種重要的海洋觀測手段。OSMAR-S系列便攜式高頻地波雷達是武漢大學(xué)電波傳播實驗室研制，其長周期工作的穩(wěn)定性和探測精度已被文必洋等(2009)、魏國妹等(2016)與Lai等(2017)的諸多研究認(rèn)可[10-12]。

該高頻雷達觀測系統(tǒng)是由兩臺武漢德威斯公司生產(chǎn)的OSMAR-S50(主要技術(shù)指標(biāo)如表1 所示)組成，分別位于深圳大鵬半島的鹿嘴山莊和惠東平海鎮(zhèn)的海龜島(圖1)。本研究采用的表層海流數(shù)據(jù)的起止時間為2015年1月1日至2015 年12月31日，數(shù)據(jù)由深圳市海洋監(jiān)測預(yù)報中心提供。該雷達系統(tǒng)的波束覆蓋扇角為180°，即每臺雷達的探測覆蓋范圍為法向方向±90°的半圓區(qū)域。矢量海流為兩臺雷達的單站探測徑向流合成，空間分辨率為0. 01°×0. 01°，數(shù)據(jù)采樣間隔為20 min。圖1給出了2015年該海域雷達系統(tǒng)數(shù)據(jù)的空間分布，在兩個雷達連線附近公共覆蓋區(qū)域的空白區(qū)數(shù)據(jù)系手動剔除，高精度區(qū)主要分布在灣外海域，因此本研究主要關(guān)注灣外海域。

表1 便攜式高頻地波雷達(OSMAR-S50)的主要技術(shù)指標(biāo)Tab.1 Main technical indexes of the portable high frequency ground wave radar (OSMAR-S50)

圖1 大亞灣附近海域雷達站位及數(shù)據(jù)有效率分布Fig.1 Distribution of the radar station and data availability in adjacent waters of Daya Bay

1.1.2 風(fēng)場數(shù)據(jù) 風(fēng)場數(shù)據(jù)為美國國家環(huán)境預(yù)測中心和國家大氣研究中心的ERDDAP服務(wù)器中提供的All Metop ASCAT中的月平均數(shù)據(jù)(https://coastwatch.pfeg.noaa.gov/erddap/griddap/erdQMwindmday.html), 空間分辨率為0.25°×0.25°。為了更直觀地分析觀測海域風(fēng)速風(fēng)向的整體變化過程，本研究還下載了單點A(22.25°N，114.75°E)的All Metop ASCAT的日平均數(shù)據(jù)，主要用于單點風(fēng)場時間序列分析。

1.1.3 地形數(shù)據(jù) 地形數(shù)據(jù)為Etopo1的水深數(shù)據(jù)，該數(shù)據(jù)是美國的國家海洋和大氣局(NOAA)的國家地球物理數(shù)據(jù)中心(NGDC)公布的數(shù)據(jù)(http://www.ngdc.noaa.gov/mgg/global/)，分辨率為1′。

1.2 數(shù)據(jù)質(zhì)量控制

高頻地波雷達在長時間的持續(xù)工作中，不僅會受到天氣、海況的影響，同時雷達系統(tǒng)本身穩(wěn)定性也會對數(shù)據(jù)質(zhì)量有影響，導(dǎo)致數(shù)據(jù)缺測、數(shù)據(jù)覆蓋率不穩(wěn)定、虛假流速等問題的存在。因此，在對雷達數(shù)據(jù)進行分析前需要對雷達數(shù)據(jù)進行嚴(yán)格的數(shù)據(jù)質(zhì)量控制，這樣才能得到較高質(zhì)量的雷達數(shù)據(jù)，為后期的數(shù)據(jù)應(yīng)用提供可靠保障。本研究依據(jù)朱大勇等(2008)和徐粲等(2014)的經(jīng)驗[13-14]對雷達數(shù)據(jù)做如下的質(zhì)量控制：①剔除異常值，即剔除觀測中東分量和北分量的流速超過2 m/s的數(shù)據(jù)和由于數(shù)據(jù)傳輸不穩(wěn)定導(dǎo)致的連續(xù)幾組完全相同的明顯異常數(shù)據(jù)；②剔除跳躍值，即與平均值之差絕對值大于2倍標(biāo)準(zhǔn)差的數(shù)據(jù)和前后兩次觀測絕對變化值大于2倍標(biāo)準(zhǔn)差的數(shù)據(jù)；③剔除孤立值，即剔除前后均無有效采樣點的數(shù)據(jù)；④利用時間與空間插值將缺失時刻和缺失的雷達網(wǎng)格點數(shù)據(jù)補充完整，并進行時間序列平滑與空間平滑，最后獲得時空連續(xù)的流場序列。通過以上質(zhì)控手段，認(rèn)為最后獲得時空連續(xù)的雷達數(shù)據(jù)是可靠的。

2 結(jié)果與討論

余流是實際海流濾掉潮流后剩下的部分，包括潮余流、風(fēng)海流和密度流等，余流對海水中物質(zhì)的輸運有重要作用，余流分析是海域泥沙輸運，污染物、赤潮等擴散的重要研究手段。本研究將時空連續(xù)的雷達表層流數(shù)據(jù)減去調(diào)和分析得到的潮流即得到了余流的時空序列，并對各月余流序列進行月平均，得到月均余流場。

2.1 1—3月余流特征

如圖2～4所示，1—3月，大亞灣海域整體風(fēng)速為6～8 m/s的東北風(fēng)，風(fēng)向逐漸由NE向轉(zhuǎn)向ENE方向(圖2、3)。大亞灣灣口余流總體表現(xiàn)為西進，且西部流速大于東部流速；灣外東南部外海有一股約28～33 cm/s的強余流，流向為SW向，基本沿著等深線方向流動，流速梯度與水深梯度近乎一致，其中2月余流強度最強，3月開始減弱，這股流在越過114.80°E后隨著等深線的轉(zhuǎn)向也逐漸向陸地側(cè)偏轉(zhuǎn)，同時隨著水深的減少，流速也逐漸減弱，流速方向變化趨勢與等深線相同，即余流的整體梯度變化與水深變化一致(圖4)。

2.2 4—6月余流特征

如圖5～7所示，4月，季風(fēng)逐漸開始轉(zhuǎn)向，大亞灣海域整體為約5.0 m/s的東風(fēng)(圖5、6)，灣外強余流消失，整體余流流速減小到一年之中的最小值(小于5 cm/s)，流向指向灣內(nèi)，灣口外東側(cè)余流近乎消失；至5月，海域內(nèi)的風(fēng)轉(zhuǎn)為3.5 m/s的南風(fēng)，灣口外的強余流逐漸加強，灣口外的流向完全反轉(zhuǎn)為NE向，可能是由于風(fēng)的轉(zhuǎn)向?qū)е?；到?月，海域內(nèi)的風(fēng)轉(zhuǎn)為約3.8 m/s的西南風(fēng)，灣外的強余流加強達到20 cm/s以上，流向為NE向，等流速線與等深線幾乎平行，流速梯度與水深梯度近乎一致，且此時風(fēng)向與流向夾角大約為45°(圖7)，與Ekman流的特點一致。此外，在4月流速變化期，余流整體較小，灣外的強余流近乎消失，流速與雷達本身的測流誤差相當(dāng)。因此，4月的流向、流速可能存在較大誤差，如A點的風(fēng)速矢量顯示4月為風(fēng)速轉(zhuǎn)向期，特別是4月下旬風(fēng)速呈不穩(wěn)定的轉(zhuǎn)向，整體風(fēng)向大體為E向，但余流流向為SW向，與風(fēng)向不一致。

2.3 7—9月余流特征

如圖8～10所示，7—8月，海域整體為1.5～2.6 m/s的南風(fēng)和西南風(fēng)(圖8、9)，余流分布與6月一致，流向為NW向，8月流速達到最大，灣外的強余流最大流速可達23 cm/s，此時灣口外的等流速線與等深線近乎平行；9月海域盛行東風(fēng)，灣口東側(cè)有一股約10 cm/s的SW向流，且在越過114.80°E后轉(zhuǎn)向灣內(nèi)(圖10)。同8月相比，余流流向再一次發(fā)生翻轉(zhuǎn)，且僅在一個月之內(nèi)完成，比4、5月的翻轉(zhuǎn)時間縮短近一個月，說明該海域的海水對東北季風(fēng)發(fā)展的響應(yīng)迅捷。

2.4 10—12月余流特征

如圖11～13所示，10—12月該海域又再次變?yōu)闁|北風(fēng)，風(fēng)速為6.5～8.8 m/s(圖11、12)，余流整體變化趨勢與1—3月變化一致，流向為SW向，流速可達30 cm/s以上，且灣外強余流逐漸增強，流速梯度逐漸變大，等流速線逐漸與等深線平行，流速梯度趨向與水深梯度一致，灣外強余流對灣口余流的影響趨于穩(wěn)定(圖13)。

圖2 2015年1—3月大亞灣附近海域ASCAT遙感風(fēng)場Fig. 2 ASCAT remote sensing wind field in adjacent waters of Daya Bay from January to March, 2015圖中“”是A點(22.25°N，114.75°E)的位置，為All Metop ASCAT日平均數(shù)據(jù)點；箭頭指向表示風(fēng)向；下同。

圖3 2015年1—3月大亞灣附近海域A點日均風(fēng)矢量變化Fig. 3 Change of the average daily wind vector at point A in adjacent waters of Daya Bay from January to March， 2015

圖4 2015年1—3月大亞灣附近海域地波雷達表層余流場Fig. 4 Surface residual current field observed by ground wave radar in adjacent waters of Daya Bay from January to March, 2015圖中藍色等值線為水深，白色箭頭指向為余流方向，下同。

圖5 2015年4—6月大亞灣附近海域ASCAT遙感風(fēng)場Fig. 5 ASCAT remote sensing wind field in adjacent waters of Daya Bay from April to June, 2015

圖6 2015年 4—6月大亞灣附近海域A點日均風(fēng)矢量變化 Fig. 6 Change of the average daily wind vector at point A in adjacent waters of Daya Bay from April to June， 2015

圖7 2015年4—6月大亞灣附近海域地波雷達表層余流場Fig. 7 Surface residual current field observed by ground wave radar in adjacent waters of Daya Bay from April to June, 2015

圖8 2015年7—9月大亞灣附近海域ASCAT遙感風(fēng)場Fig. 8 ASCAT remote sensing wind field in adjacent waters of Daya Bay from July to September, 2015

圖9 2015年7—9月大亞灣附近海域A點日均風(fēng)矢量變化Fig. 9 Change of the average daily wind vector at point A in adjacent waters of Daya Bay from July to September, 2015

圖10 2015年7—9月大亞灣附近海域地波雷達表層余流場Fig. 10 Surface residual current field observed by ground wave radar in adjacent waters of Daya Bay from July to September, 2015

2.5 討論

2.5.1 大亞灣附近海域余流季節(jié)變化特征 本研究通過對2015年全年的高頻地波雷達連續(xù)觀測數(shù)據(jù)分析發(fā)現(xiàn)，大亞灣海域表層余流的變化主要呈冬、夏兩種形態(tài)，其中，冬季形態(tài)持續(xù)半年以上(從10月到次年3 月)，流向為SW向，流速可達30 cm/s以上；夏季形態(tài)僅在5—8 月間出現(xiàn)，流向為NE向，流速可達23 cm/s。本研究的余流變化特征與朱大勇等[13]的觀測結(jié)論類似，而且這種變化過程與季風(fēng)的轉(zhuǎn)變基本同步。

2.5.2 大亞灣附近海域余流驅(qū)動機制討論 本研究的余流空間分布數(shù)據(jù)顯示：10月到次年3月，風(fēng)向與強余流的方向基本一致，此時強余流特征與粵東沿岸流基本一致，且余流總在越過114.80°E后逐漸轉(zhuǎn)為向岸流，這與水深等值線在這里發(fā)生向岸轉(zhuǎn)向有關(guān)；6—8月，風(fēng)速較小，但風(fēng)向與流向的夾角小于90°，體現(xiàn)了風(fēng)對余流的影響，同時也發(fā)現(xiàn)風(fēng)只能影響表層的流速大小而不能改變其主要流向，這與陳希榮等(2018)在珠江口海域的觀測結(jié)論[15]相同。此外，除4月與9月，余流等值線與水深等值線近乎平行，流速梯度均為NW—SE向，這與研究區(qū)域地形梯度的分布一致，又體現(xiàn)了地形對余流的影響。綜合全年的余流變化特征，不難看出大亞灣附近海域的余流變化主要受東亞季風(fēng)影響，其變化特征與每年10月至次年3月盛行的東北季風(fēng)，6—8月盛行的西南季風(fēng)相一致。此外，地形對余流分布也存在重要影響。這與其他學(xué)者的研究相一致，如徐粲等在南黃海輻射沙脊群的研究發(fā)現(xiàn)，近岸海域余流分布主要受潮汐與水下地形共同作用的影響，而靠外海域余流格局主要受風(fēng)場的控制[14]；史軍強(2014)在長江口海域的研究也同樣證明了風(fēng)場和地形對余流分布的影響[16]。

此外，粵東沿岸流與珠江沖淡水對大亞灣灣口余流空間分布具有調(diào)控作用。孫湘平(2006)、張志欣(2014)和舒業(yè)強等(2018)的研究結(jié)果認(rèn)為廣東沿岸存在著一股沿岸流即粵東沿岸流，它主要是隨季風(fēng)風(fēng)向變化，冬季為SW向流，流幅窄；夏季為NE向流，流幅寬，且珠江沖淡水在夏季可影響至大亞灣外海域[17-19]。10月至次年3月研究海域主要受到東北風(fēng)影響下的SW向強沿岸流影響，粵東沿岸流在越過114.80°E后，由于受地形的影響，部分流發(fā)生轉(zhuǎn)向作為補償流逐漸轉(zhuǎn)向成為垂直岸線的余流；4月東北風(fēng)逐漸減弱，到5月初已表現(xiàn)為較弱的西南風(fēng)，但整體上表現(xiàn)為東北風(fēng)，由于東北季風(fēng)的減弱，Ekman效應(yīng)制約沖淡水向外海擴散的能力明顯減弱，沖淡水漸漸影響大亞灣灣口以外的海域； 6—8月隨著珠江徑流量的增加與西南季風(fēng)的建立，在西南季風(fēng)的驅(qū)使及粵東沿岸流的誘導(dǎo)下，表層沖淡水向東擴散至大亞灣附近海域，這與前人的研究結(jié)果[20]是一致的。歐素英(2005)、曾淦寧等(2005)、Ou等(2007)、楊陽等(2014)與Bai等(2015)的研究表明，本研究區(qū)域，大亞灣附近海域的余流受到了珠江沖淡水和粵東沿岸流的共同作用[20-24]。

2.5.3 高頻地波雷達在沿岸流觀測中的作用 對于粵東沿岸流的觀測大多基于一次或多次的航次觀測和定點觀測數(shù)據(jù)，盡管衛(wèi)星遙感觀測和數(shù)值模擬方法很大程度上豐富了我們對粵東沿岸流的觀測手段，但缺少大范圍連續(xù)觀測數(shù)據(jù)，特別是長時間序列的連續(xù)海流觀測數(shù)據(jù)，導(dǎo)致我們無法進行準(zhǔn)確直觀的觀測。本研究利用的長周期高頻地波雷達觀測分析余流變化特征與Fang等(1998)對粵東海域沿岸流的研究結(jié)果相似，即沿岸流主要受季節(jié)性反轉(zhuǎn)的風(fēng)場強迫，夏季呈現(xiàn)NE向，冬季為W向流等特征[25]，但又略有差異，主要體現(xiàn)在本研究的觀測結(jié)果顯示粵東沿岸流在大亞灣沿岸冬季為SW向，這可能與雷達觀測結(jié)果比較精細化有關(guān)。本研究的觀測結(jié)果在一定程度上說明高頻地波雷達對粵東沿岸流、珠江沖淡水等大尺度流系觀測中的突出作用，對比以往只能通過現(xiàn)場觀測結(jié)合數(shù)值模擬的方式，雷達有著顯著優(yōu)勢。由于雷達海流探測精度的限制，在余流較小的4月、5月和9月，余流的分布特征可能存在較大誤差，分析結(jié)論僅供參考。

圖11 2015年10—12月大亞灣附近海域ASCAT遙感風(fēng)場Fig. 11 ASCAT remote sensing wind field in adjacent waters of Daya Bay from October to December, 2015

圖12 2015年10—12月大亞灣附近海域A點日均風(fēng)矢量變化 Fig. 12 Change of the average daily wind vector at point A in adjacent waters of Daya Bay from October to December， 2015

圖13 2015年10—12月大亞灣附近海域地波雷達表層余流場Fig. 13 Surface residual current field observed by ground wave radar in adjacent waters of Daya Bay from October to December, 2015

3 結(jié)論

本研究通過利用T_tide工具包對高頻地波雷達觀測數(shù)據(jù)進行潮流調(diào)和分析，獲得如下結(jié)論：

(1)大亞灣海域表層余流的變化特征主要呈冬、夏兩種形態(tài)，其變化過程與季風(fēng)的轉(zhuǎn)變基本同步。其中冬季形態(tài)持續(xù)半年以上(從10月到次年3月)，流向為SW向，流速可達30 cm/s以上; 夏季形態(tài)僅在5—8月間出現(xiàn)，流向為NE向，流速可達23 cm/s，這與該區(qū)域受到東亞季風(fēng)控制，每年10月至次年3月盛行東北季風(fēng)，6—8月盛行西南季風(fēng)相一致。

(2)大亞灣附近海域的余流分布特征不僅與海表面風(fēng)場變化相關(guān)，而且與海底地形也存在顯著相關(guān)關(guān)系。此外，粵東沿岸流與珠江沖淡水對大亞灣灣口余流空間分布具有調(diào)控作用。

(3)高頻地波雷達可在粵東沿岸流、珠江沖淡水等大尺度流系的觀測中發(fā)揮重要作用。