朱 佳 胡建宇 陳照章 孫振宇, 2 楊龍奇 彭帥星 吉皎月

大亞灣海域夏、冬季的潮汐特征及余水位與風(fēng)的相關(guān)性初步探討*

朱 佳1胡建宇1①陳照章1孫振宇1, 2楊龍奇1彭帥星1吉皎月1

(1. 廈門大學(xué)近海海洋環(huán)境科學(xué)國家重點(diǎn)實驗室 廈門 361012; 2. 南方海洋科學(xué)與工程廣東省實驗室 珠海 519082)

大亞灣及其鄰近海域冬、夏季各14個臨時水位觀測點(diǎn)1個月的實測潮位資料顯示: 各站的水位曲線均呈現(xiàn)明顯的“雙峰”現(xiàn)象, 且灣頂比灣口更為明顯。本文采用了調(diào)和分析方法, 給出M2、S2、K1、O1四個主要分潮及M4、M6、2MS6三個淺水分潮的振幅和遲角同潮圖, 分析大亞灣的主要潮汐特征, 探討了淺水分潮對雙峰結(jié)構(gòu)的貢獻(xiàn), 并采用交叉譜分析對余水位與風(fēng)的相關(guān)性進(jìn)行了討論。結(jié)果表明: (1)大亞灣海域各主要分潮振幅均由灣口向灣頂遞增; 高潮發(fā)生時間由灣口向灣頂推遲; 漲潮歷時均大于落潮歷時; 平均潮差在灣頂達(dá)到最大; (2)大亞灣內(nèi)屬于不正規(guī)半日潮, 而考洲洋及其灣外海域則屬于不正規(guī)全日潮; (3)大亞灣內(nèi)淺水效應(yīng)明顯, 從灣口至灣頂, 六分之一日分潮的振幅呈5—7倍的增長, 主導(dǎo)了大亞灣潮波系統(tǒng)的形變; (4)分潮重構(gòu)結(jié)果顯示, 四分之一日和六分之一日淺水分潮(尤其是2MS6分潮)的異常增長, 是導(dǎo)致大亞灣潮汐雙峰現(xiàn)象的主要原因; (5)冬季大亞灣內(nèi)各點(diǎn)的余水位與風(fēng)速呈現(xiàn)正相關(guān), 相關(guān)系數(shù)均在0.53以上; (6)周期為0.45—0.53 d的沿岸風(fēng)對各站余水位的影響最大。

大亞灣; 水位雙峰; 潮汐特征; 淺水分潮; 余水位與風(fēng)的相關(guān)性

大亞灣是南海沿岸的重要海灣, 北靠海岸山脈, 東、西兩側(cè)分別為平海半島與大鵬半島, 灣內(nèi)水深自北向南逐漸增加, 中部水深10余米, 灣口水深可達(dá)20多米。

鑒于其重要地理位置, 已有多位學(xué)者開展了大亞灣海域潮汐、潮流以及水交換的研究工作。李立(1998)應(yīng)用譜分析和線性模型研究了大亞灣內(nèi)冬季亞潮水位變化的特征和原因。楊國標(biāo)(2001)運(yùn)用19個站的實測數(shù)據(jù)分析了大亞灣六個主要分潮的潮流及余流特征。周巧菊(2007)采用了POM(princeton ocean model)模型研究了大亞灣溫排水交換過程, 并對潮汐進(jìn)行了模擬。王聰?shù)?2009)通過數(shù)值模擬研究了風(fēng)場和潮致余流對大亞灣水交換能力的影響。

雙峰結(jié)構(gòu)是大亞灣潮波的主要特征, 即在半日潮汐周期中, 水位停止上升, 略有下降, 然后再次上升, 而后又迅速下降至低潮的狀態(tài)(Godin, 1993)。這種由于高頻振蕩疊加半日潮產(chǎn)生的雙峰現(xiàn)象出現(xiàn)在不少港灣中, 如美國的Buzzards灣(Redfield, 1980), 英國的Solent海峽、荷蘭的Den Helde港(Pugh, 1987)和英格蘭南海岸的Swanage灣(Godin, 1993)都觀測到了該現(xiàn)象。Dyer(1973)指出, 當(dāng)潮波傳播至淺水時, 受到漏斗效應(yīng)和淺水成分的影響, 越靠近上游潮差越大。Boon III (1975)認(rèn)為M6和M8這類高頻潮可導(dǎo)致漲潮和落潮歷時的差異。Blanton等(2002)的研究證明在美國佛羅里達(dá)州的Charles Creek, M6分潮對潮流形變的影響甚至大于M4分潮。Kim等(2001)推測Narragansett灣的M4分潮振幅可能因共振增大了2—3倍。雙峰現(xiàn)象目前在國內(nèi)的報道和研究并不多, 且多集中于大亞灣。多位學(xué)者根據(jù)2個潮位站1個月的水位數(shù)據(jù), 對潮汐做了系列研究, 指出M6分潮的異常增長是引起大亞灣潮位雙峰的主要原因(Song et al, 2016; 刁希梁等, 2017); 用三維潮汐、潮流數(shù)值模型(武文等, 2017)和數(shù)學(xué)解析方法(嚴(yán)聿晗等, 2017)揭示了水位雙峰、漲潮雙峰以及高頻分潮的強(qiáng)耗散等現(xiàn)象。由以上研究可見, 淺水分潮的振幅雖不大, 但在港灣內(nèi)的作用卻不容忽視, 其在港灣內(nèi)的淺水效應(yīng)及與地形呼應(yīng)導(dǎo)致的潮汐共振, 是引起雙峰結(jié)構(gòu)、潮波變形以及振幅劇增的重要原因。

近海海域的海平面及水位變化, 除受地形、淺水成分等的影響之外, 對風(fēng)的響應(yīng)也較明顯。Chuang等(1983)調(diào)查了墨西哥灣北部陸架的沿岸風(fēng)和海陸風(fēng)(包括向岸風(fēng)和離岸風(fēng))對海平面變化的影響, 發(fā)現(xiàn)不同的響應(yīng)主要是由水深變化引起的。Craig (1989a, 1989b)建立了海洋對固定頻率的周期性風(fēng)強(qiáng)迫的響應(yīng)理論。根據(jù)觀測和數(shù)值模擬結(jié)果, 多位學(xué)者探討了海平面和潮流在不同海域里對海陸風(fēng)的響應(yīng)(Clancy, 1979; Vesecky, 1997; Gallop, 2012; Lin, 2016), 特別是在臨界緯度附近(約30°N和30°S),由于晝夜強(qiáng)迫的頻率與當(dāng)?shù)貞T性頻率重合而出現(xiàn)近共振, 還經(jīng)常能觀察到增強(qiáng)的海洋響應(yīng)(Simpson, 2002; Zhang, 2009; Nam, 2013)。

目前對大亞灣水文動力方面的研究雖包含了實測、模型及數(shù)學(xué)解析等方法, 但仍需更多的實測數(shù)據(jù)進(jìn)行補(bǔ)充和佐證。因此, 我們在大亞灣及其鄰近海域布設(shè)了冬、夏季各14個臨時水位觀測點(diǎn), 進(jìn)行了一個月的同步觀測, 旨在分析冬季和夏季的潮汐特征、季節(jié)差異及成因, 探討淺水分潮及共振對潮汐畸變的貢獻(xiàn)和影響, 并初步討論余水位與風(fēng)速、風(fēng)向及周期的相關(guān)性。

1 數(shù)據(jù)及方法說明

1.1 水位及余水位數(shù)據(jù)

選用大亞灣內(nèi)冬季、夏季各14個臨時水位觀測點(diǎn)的連續(xù)水位數(shù)據(jù), 水位觀測點(diǎn)分布見圖1, 其基準(zhǔn)面是以平均海平面為準(zhǔn), 所有站位都統(tǒng)一。雖然部分站點(diǎn)的實際觀測時間遠(yuǎn)大于1個月, 為求一致, 我們均截取一個月的長度。冬季的水位觀測時間為2018年1月8日—2月8日, 夏季的觀測時間為2018年5月17日—6月17日。

圖1 大亞灣內(nèi)冬季、夏季的臨時水位觀測點(diǎn)分布圖

注:表示風(fēng)向角度; T01—T14、HY001、HY002、HY004、HY005表示站位名稱

水位觀測使用的儀器為RBR TDR-2050潮位儀及KELLER DCX-22水位記錄儀。所用儀器在施測之前均通過調(diào)試檢驗, 且比測結(jié)果一致性良好。數(shù)據(jù)采樣間隔分為兩種, RBR TDR-2050采樣間隔為1 min, KELLER DCX-22采樣間隔為5 min, 兩者精度均為0.01 m。

使用t_tide Matlab軟件包(Pawlowicz, 2002)進(jìn)行潮汐分析。與經(jīng)典的諧波分析相似, t_tide預(yù)測潮汐信號系與天文參數(shù)相關(guān)的特定頻率的有限組正弦信號的和(Pawlowicz, 2002)。余水位由原始水位減去潮汐信號得到。在t_tide諧波分析中使用了35個潮汐成分, 頻率為每小時0.0015 (MM潮汐成分)至0.3220 (M8潮汐成分)周期。

1.2 風(fēng)場數(shù)據(jù)

采用對應(yīng)冬季、夏季水位觀測期間的歐洲氣象中心(European Centre for Medium-Range Weather Forecasts, ECMWF)風(fēng)場數(shù)據(jù), 其空間分辨率為0.125°×0.125°, 時間分辨率為3 h。為減小全球模式近岸區(qū)域個別網(wǎng)格點(diǎn)的誤差, 取調(diào)查海域114°—115°E, 22.0°—23.0°N區(qū)域內(nèi)的平均風(fēng)場進(jìn)行計算。

文中對余水位和風(fēng)應(yīng)用了相關(guān)分析和交叉譜分析方法, 當(dāng)兩類數(shù)據(jù)時間分辨率不同時, 選較低時間分辨率做插值計算。例如, 水位數(shù)據(jù)和風(fēng)場數(shù)據(jù)的時間分辨率分別為1 h和3 h, 兩者進(jìn)行相關(guān)分析時取時間分辨率為3 h。

2 潮汐特征

2.1 潮位特征

選擇水位變化特征較有代表性的灣頂T06站為例, 繪制冬季和夏季的水位時序曲線。從圖2中可看出水位過程曲線在每個太陰日經(jīng)歷了2次高潮和2次低潮,且每次高潮或低潮的水位高度均不相同, 故存在潮汐的日不等現(xiàn)象。此外, 無論是高潮還是低潮, 水位曲線存在明顯的雙峰現(xiàn)象。雙峰之間的時間間隔約為3.77—4.32 h不等; 雙峰值間的波峰和波谷的平均水位差約為40 cm, 最大水位差可達(dá)93 cm。

圖2 大亞灣T06冬季(a)和夏季(b)水位過程曲線

2.2 潮汐特征

在冬季和夏季, 平均潮差由灣口向灣頂逐漸增大, 在灣頂達(dá)到最大。各站漲潮歷時均遠(yuǎn)大于落潮歷時。整體上看, 灣頂?shù)母叱睍r要早于灣口, 灣頂?shù)牡统睍r要晚于灣口(圖3)。

冬季, 灣內(nèi)平均潮差最大值出現(xiàn)在T06站, 為115 cm; 而T09、T11和T14站的平均潮差則最小, 均為80 cm。各站平均漲潮歷時大于平均落潮歷時, 平均漲潮歷時在6.38—7.65 h范圍內(nèi), 平均落潮歷時在4.78—7.11 h范圍內(nèi)。

夏季, 平均潮差最大值出現(xiàn)在T06站, 為121 cm; 而T10站的平均潮差最小, 為85 cm。平均漲潮歷時大于平均落潮歷時, 平均漲潮歷時在6.42—7.43 h范圍內(nèi), 平均落潮歷時在4.98—5.96 h范圍內(nèi)。

海洋圖集編委會(2004)曾指出: 在南海海域, 幾乎沒有正規(guī)半日潮出現(xiàn), 而是以不正規(guī)全日潮類型占優(yōu)勢, 并有明顯的日不等現(xiàn)象; 而考洲洋的地理位置恰好位于南海北部不正規(guī)半日潮和不正規(guī)全日潮的分界點(diǎn)上。這就是大亞灣灣內(nèi)和灣口為不正規(guī)半日潮, 而考洲洋灣內(nèi)、灣口有時是不正規(guī)全日潮的原因。

圖3 冬季(a—c)和夏季(d—f)各站潮汐分布特征圖

圖4 冬季(a)和夏季(b)大亞灣及其鄰近海域潮汐特征值分布

2.3 同潮圖

為了解大亞灣及其鄰近海域的潮汐特征, 繪制M2、S2、K1、O1、P1、Q1、K2、N2、M4、M6、MS4、2MS6十二個分潮的同潮圖, 但限于篇幅, 本文僅列出冬季、夏季M2、S2、K1、O1四個主要分潮及M4、M6、2MS6三個淺水分潮的同潮圖(圖5)。由結(jié)果可得: 在14個水位測站中, 振幅最大的分潮為M2和K1, 其次是O1分潮, 而后依次是S2分潮, P1和M4分潮, N2、Q1、MS4、M6、K2分潮, 最后是2MS6分潮; 大多數(shù)分潮的振幅均由灣口至灣頂逐漸增大。

冬季, M2分潮振幅由大亞灣灣口T09站的30 cm到灣頂T05站的35 cm, 增長了15%, K1、O1、S2分潮振幅分別增長了6%、4%和27%, 而振幅增幅較大的是淺水分潮M4、MS4(圖略)、M6和2MS6分潮, 增幅分別為81%、97%、525%和616%。夏季, 從灣口到灣頂, M2分潮振幅增幅為14%, K1、O1、S2分潮振幅增幅為3%、2%和16%, 而淺水分潮M4、MS4(圖略)、M6和2MS6的增幅則分別達(dá)到60%、86%、547%和707%。

遲角方面, M2、S2、K1、O1和P1、Q1、K2、N2(圖略)這八個主要分潮的等遲角線在大亞灣內(nèi)分布稀疏, 從灣口到灣頂遲角變化不大, 尤其是全日分潮, 灣頂站位的位相比灣口站位僅滯后2°—4°, 說明全日分潮的潮波在大亞灣內(nèi)傳播最快, 半日分潮次之; 而M4、M6、MS4(圖略)和2MS6四個淺水分潮的遲角分布呈現(xiàn)從考洲洋灣口到大亞灣灣頂迅速增大的趨勢。尤其是M6和2MS6這兩個六分之一日淺水分潮, 其等值線在T11—T10站之間分布密集, 表明在這些海域, 六分之一日分潮的傳播速度最慢。

2.4 潮位雙峰現(xiàn)象的成因分析

2.4.1 分潮共振現(xiàn)象 由上述分析可見, 在潮波傳入大亞灣之后, 六分之一日分潮的振幅呈5—7倍的增長, 主導(dǎo)了大亞灣潮波系統(tǒng)的形變。這是由于分潮共振現(xiàn)象引起的, 即當(dāng)某分潮的頻率與半封閉海灣的固有(本征)頻率相近時, 激發(fā)出潮汐共振(董禮先等, 1999a, 1999b; Zhong, 2008; 刁希梁等, 2017; 武文等, 2017; 嚴(yán)聿晗等, 2017)。

根據(jù)四分之一波長理論公式(Proudman, 1952; Defant, 1961; Song, 2016)計算大亞灣的固有共振周期。

2.4.2 分潮重構(gòu) 從調(diào)查海域冬季、夏季各14個水位站所繪制水位過程曲線圖上看, 無論是高高潮還是低高潮, 無論是大潮期還是小潮期, 都存在明顯的水位雙峰現(xiàn)象, 而且這一現(xiàn)象在灣頂要明顯強(qiáng)于灣口。為探究到底是哪些分潮導(dǎo)致雙峰現(xiàn)象的出現(xiàn), 我們利用調(diào)和常數(shù)對各站潮位進(jìn)行了重構(gòu)。在重構(gòu)過程中逐步剔除某個分潮, 而后用剩余分潮回構(gòu)潮位。經(jīng)過多次嘗試, 發(fā)現(xiàn)四分之一日和六分之一日的淺水分潮(尤其是2MS6分潮)的異常增長是導(dǎo)致大亞灣及其鄰近海域潮汐雙峰現(xiàn)象的主要原因。

以水位雙峰現(xiàn)象最明顯的灣頂站位T06站為例。首先剔除振幅增幅最大的淺水分潮2MS6, 回構(gòu)曲線中的雙峰程度已大幅削減, 僅在個別波峰處存在小于10 cm水位差(圖略); 而后剔除M6和2MS6分潮, 在回構(gòu)曲線(圖6b)中可見, 雙峰現(xiàn)象幾乎消失, 僅有略微殘存的雙峰且水位差小于2 cm; 而在剔除M4、MS4、M6、2MS6這四個淺水分潮之后, 圖6c中的回構(gòu)潮位曲線中已完全不見雙峰現(xiàn)象。由此可見, 大亞灣內(nèi)潮位的雙峰現(xiàn)象是由M4、MS4、M6、2MS6產(chǎn)生的, 尤其是M6及2MS6這兩個六分之一日分潮占主導(dǎo)的,說明六分之一日分潮主導(dǎo)了大亞灣及其鄰近海域潮波系統(tǒng)的形變, 該結(jié)論進(jìn)一步佐證了以往的研究結(jié)果 (Song et al, 2016; 武文等, 2017)。

注: a: 冬季T06站的實測潮位曲線; b: 剔除M6和2MS6分潮的潮位回構(gòu)曲線; c: 剔除M6、2MS6、M4、MS4四個淺水分潮的潮位回構(gòu)曲線

3 討論

3.1 余水位與風(fēng)速的相關(guān)性

為了探討余水位與風(fēng)速的關(guān)系, 取冬季調(diào)查期間的風(fēng)場為例。由圖7可見(選取灣口、中部及灣頂4個測點(diǎn)展示), 冬季大亞灣內(nèi)水位測點(diǎn)的余水位與風(fēng)速曲線的變化趨勢較吻合, 風(fēng)速大時余水位也相應(yīng)升高, 兩者相關(guān)系數(shù)均在0.53以上, 呈現(xiàn)正相關(guān)關(guān)系。

3.2 余水位與風(fēng)向及周期的相關(guān)性

為進(jìn)一步探究余水位與風(fēng)的關(guān)系, 我們對兩者進(jìn)行交叉譜分析, 目的是找出與余水位相關(guān)性最高的風(fēng)向與頻率(周期)。如圖1所示, 風(fēng)向以西風(fēng)為0°(指向正東為0°), 逆時針旋轉(zhuǎn)增大, 東風(fēng)為180°(指向正西為180°), 例如,=30°即表示東偏北30°。由圖8可見, 風(fēng)與余水位的相關(guān)性最高值約為0.75, 其對應(yīng)的風(fēng)頻率為1.9—2.2 d-1(平均2.1 d-1), 風(fēng)轉(zhuǎn)角為20°—79°(平均=36°)。也就是說, 對各站余水位影響最大的風(fēng)周期約為0.45—0.53 d, 對各站余水位影響最大的平均風(fēng)向為東偏北36° (或36°+180°=216°), 該風(fēng)向基本與大亞灣岸線平行(標(biāo)注于圖1), 因此可以認(rèn)為, 沿岸風(fēng)對各站余水位的影響最大。此外, 圖8中還可見相關(guān)系數(shù)大于0.6的等值線分布在轉(zhuǎn)角10°—120°、頻率1.7—3.5 d-1的范圍里, 可能與陸海風(fēng)與局地天氣變化有關(guān)?？紤]到北半球的大氣變量在冬季最為活躍, 我們將在接下來的工作里, 用小波分析進(jìn)一步研究余水位與其他氣象因子以及陸海風(fēng)之間的關(guān)系。

4 結(jié)論

本文根據(jù)大亞灣及其鄰近海域冬季、夏季各14個臨時水位觀測點(diǎn)的同步觀測數(shù)據(jù), 分析其潮汐特征及成因, 結(jié)論總結(jié)如下:

(1)冬季和夏季的水位曲線存在明顯雙峰現(xiàn)象, 且灣頂明顯強(qiáng)于灣口。

圖7 余水位與風(fēng)速的相關(guān)性

注:表示相關(guān)系數(shù)

圖8 T01、T03、T05、T07站的余水位與風(fēng)向及頻率的相關(guān)性

注: 橫坐標(biāo)的轉(zhuǎn)角表示風(fēng)向, 風(fēng)向按逆時針旋轉(zhuǎn), 0°為西風(fēng)(指向正東為0°)

(2)各站漲潮歷時均遠(yuǎn)大于落潮歷時。平均潮差由灣口向灣頂逐漸增大, 在灣頂達(dá)到最大。

(3)大亞灣內(nèi)各站潮型屬于不正規(guī)半日潮, 而考洲洋灣內(nèi)和灣口海域則屬于不正規(guī)全日潮。

(4)大多數(shù)分潮振幅均呈現(xiàn)從灣口向灣頂遞增的趨勢; 各站高潮發(fā)生時間由灣口向灣頂推遲。

(5)灣內(nèi)淺水效應(yīng)明顯, 淺水分潮的振幅遞增幅度遠(yuǎn)大于八個主要分潮。從灣口至灣頂, M6和2MS6分潮冬季增幅分別可達(dá)525%和616%, 夏季分別可達(dá)547%和707%。

(6)由分潮重構(gòu)和潮汐共振分析結(jié)果來看, 六分之一日淺水分潮(尤其是2MS6分潮)的異常增長是導(dǎo)致大亞灣潮汐雙峰現(xiàn)象的主要原因。

(7)冬季大亞灣內(nèi)的余水位與風(fēng)速呈現(xiàn)正相關(guān), 兩者相關(guān)系數(shù)r均在0.53以上;

(8)風(fēng)向為東偏北36°(或216°)且周期為0.45—0.53 d的沿岸風(fēng)對余水位影響最大, 相關(guān)系數(shù)值為0.75。此外, 陸海風(fēng)與局地天氣變化對余水位的影響也不可忽視。

刁希梁, 丁 揚(yáng), 鮑獻(xiàn)文, 2017. 大亞灣海域潮位“雙峰”現(xiàn)象生成機(jī)制研究. 中國海洋大學(xué)學(xué)報, 47(9): 1—10

王 聰, 林 軍, 陳丕茂等, 2009. 年平均風(fēng)場作用下大亞灣水交換的數(shù)值模擬. 上海海洋大學(xué)學(xué)報, 18(3): 351—358

嚴(yán)聿晗, 武 文, 宋德海等, 2017. 大亞灣的潮汐動力學(xué)研究——Ⅱ. 潮位和潮流雙峰現(xiàn)象的產(chǎn)生機(jī)制. 熱帶海洋學(xué)報, 36(3): 46—54

李 立, 1998. 大亞灣冬季水位的亞潮變化及其與南海的耦合. 臺灣海峽, 17(4): 383—390

楊國標(biāo), 2001. 大亞灣海區(qū)潮流運(yùn)動特征. 人民珠江, (1): 30—32

武 文, 嚴(yán)聿晗, 宋德海, 2017. 大亞灣的潮汐動力學(xué)研究——I. 潮波系統(tǒng)的觀測分析與數(shù)值模擬. 熱帶海洋學(xué)報, 36(3): 34—45

周巧菊, 2007. 大亞灣海域溫排水三維數(shù)值模擬. 海洋湖沼通報, (4): 37—46

海洋圖集編委會, 2004. 南海海洋圖集——水文分冊. 北京: 海洋出版社

董禮先, 蘇紀(jì)蘭, 1999a. 象山港潮波響應(yīng)和變形研究I. 觀測和分析. 海洋學(xué)報, 21(1): 1—10

董禮先, 蘇紀(jì)蘭, 1999b. 象山港潮波響應(yīng)和變形研究II. 象山港潮波數(shù)值研究. 海洋學(xué)報, 21(2): 1—8

Blanton J O, Lin G Q, Elston S A, 2002. Tidal current asymmetry in shallow estuaries and tidal creeks. Continental Shelf Research, 22(11—13): 1731—1743, doi: 10.1016/S0278- 4343(02)00035-3

Boon III J D, 1975. Tidal discharge asymmetry in a salt marsh drainage system. Limnology and Oceanography, 20(1): 71—80, doi: 10.4319/lo.1975.20.1.0071

Chuang W S, Wiseman Jr W J, 1983. Coastal sea level response to frontal passages on the Louisiana-Texas shelf. Journal of Geophysical Research: Oceans, 88(C4): 2615—2620

Clancy R M, Thompson J D, Hurlburt H E, 1979. A model of mesoscale air-sea interaction in a sea breeze-coastal upwelling regime. Monthly Weather Review, 107(11): 1476—1505

Craig P D, 1989a. Constant-eddy-viscosity models of vertical structure forced by periodic winds. Continental Shelf Research, 9(4): 343—358

Craig P D, 1989b. A model of diurnally forced vertical current structure near 30° latitude. Continental Shelf Research, 9(11): 965—980

Defant A, 1961. Physical Oceanography. New York, USA: Pergamon Press, 1—598

Dyer K R, 1973. Estuaries: A Physical Introduction. Chichester, UK: Wiley-Interscience, 1—140

Gallop S L, Verspecht F, Pattiaratchi C B, 2012. Sea breezes drive currents on the inner continental shelf off southwest Western Australia. Ocean Dynamics, 62(4): 569—583

Godin G, 1993. An investigation of the phenomenon of double high water or double low water at some harbours. Deutsche Hydrografische Zeitschrift, 45(2): 87—106, doi: 10.1007/ BF02226316

Kim H S, Swanson J C, 2001. Modeling of double-flood currents in the Sakonnet River. In: Proceedings of the Seventh International Conference on Estuarine and Coastal Modeling. St. Petersburg, Florida: American Society of Civil Engineers, 418—433

Lin H Y, Hu J Y, Zhu J, 2017. Tide- and wind-driven variability of water level in Sansha Bay, Fujian, China. Frontiers of Earth Science, 11(2): 332—346, doi: 10.1007/s11707-016-0588-x

Nam S, Send U, 2013. Resonant diurnal oscillations and mean alongshore flows driven by sea/land breeze forcing in the coastal southern California Bight. Journal of Physical Oceanography, 43(3): 616—630

Pawlowicz R, Beardsley B, Lentz S, 2002. Classical tidal harmonic analysis including error estimates in MATLAB using T_TIDE. Computers & Geosciences, 28(8): 929—937

Proudman J, 1952. Dynamical Oceanography. London, UK: Methuen, 1—409

Pugh D T, 1987. Tides, Surges and Mean Sea-Level. New York, USA: John Wiley, 1—472

Redfield A C, 1980. The Tides of the Waters of New England and New York. Woods Hole, USA: Woods Hole Oceanographic Institution, 1—108, doi: 10.1575/1912/1136

Simpson J H, Hyder P, Rippeth T P, 2002. Forced oscillations near the critical latitude for diurnal-inertial resonance. Journal of Physical Oceanography, 32(1): 177—187

Song D H, Yan Y H, Wu W, 2016. Tidal distortion caused by the resonance of sexta-diurnal tides in a micromesotidal embayment. Journal of Geophysical Research: Oceans, 121(10): 7599—7618

Vesecky J F, Teague C C, Onstott R G, 1997. Surface current response to land-sea breeze circulation in Monterey Bay, California as observed by a new multifrequency HF radar. In: Proceedings of Oceans '97. MTS/IEEE Conference Proceedings. Halifax, NS, Canada: IEEE, 1019—1024

Zhang X Q, DiMarco S F, Smith IV D C, 2009. Near-resonant ocean response to sea breeze on a stratified continental shelf. Journal of Physical Oceanography, 39(9): 2137—2155

Zhong L J, Li M, Foreman M G G, 2008. Resonance and sea level variability in Chesapeake Bay. Continental Shelf Research, 28(18): 2565—2573

A PRELIMINARY STUDY ON THE TIDAL CHARACTERISTICS IN SUMMER AND WINTER AND THE CORRELATION BETWEEN RESIDUAL WATER LEVEL AND WIND IN DAYA BAY

ZHU Jia1, HU Jian-Yu1, CHEN Zhao-Zhang1, SUN Zhen-Yu1, 2, YANG Long-Qi1, PENG Shuai-Xing1, JI Jiao-Yue1

(1. State Key Laboratory of Marine Environmental Science(Xiamen University), Xiamen 361012, China;2. Southern Marine Science and Engineering Guangdong Laboratory (Zhuhai), Zhuhai 519082, China)

One-month tidal data were measured at 14 temporary observation stations in winter and summer in Daya Bay, and the dada show that the water level at each station presented obvious bimodal phenomenon, which is more obvious at the top than at the mouth of the bay. By using harmonic analysis method, the cotidal charts of the amplitude and the phase lag of four main tidal components M2, S2, K1and O1, and three shallow water tidal components M4, M6, and 2MS6were plotted; and the major tidal characteristics were analyzed. In addition, the contribution of shallow water tide to the bimodal structure was deduced through the tidal constituent reconstruction, and the correlation between residual water level and wind was discussed. The results are as followed. First, the amplitudes of each tidal constituent in the Daya Bay increased from the mouth to the endpoint of the bay, the occurrence time of high tide was postponed from the mouth to the endpoint of the bay, and the flood duration lasts much longer than the ebb duration, the mean tidal range reached its maximum at the endpoint. Secondly, the tidal type in the Daya Bay was irregular semidiurnal tide, while those in the Kaozhouyang and the areas outside the bay were irregular diurnal tide. Thirdly, the shallow water effect dominated the deformation of the tidal wave system in the Daya Bay, and the amplitude of the quarter-diurnal and sexta-diurnal tides increased 5-7 times from the mouth to the endpoint of the bay. Fourthly, the abnormal growth of the quarter-diurnal and sexta-diurnal tides (especially 2MS6) was the main cause of the bimodal phenomenon in the bay. Fifthly, in winter, the residual water level in Daya Bay was correlated positively with the wind speed, and the correlation coefficient was above 0.53. Sixthly, the alongshore wind with a period of 0.45—0.53 d had the greatest influence on the residual water level of each station.

Daya Bay; water level bimodality; tidal characteristics; shallow water constituents; correlation between residual water level and wind

* 國家自然科學(xué)基金項目, 41776027號; 惠州市2017年市級統(tǒng)籌油補(bǔ)資金項目, F2017-01-1號。朱 佳, 工程師, E-mail: zhujia@xmu.edu.cn

胡建宇, 教授, E-mail: hujy@xmu.edu.cn

2020-08-06,

2020-11-12

P717

10.11693/hyhz20200800230