榮成外海海流特征及影響因素

魏詩晏, 楊 偉, 趙 亮

榮成外海海流特征及影響因素

魏詩晏1, 楊 偉2, 趙 亮1

(1. 天津科技大學 海洋與環境學院, 天津 300457; 2. 天津大學 海洋科學與技術學院, 天津 300072)

榮成是我國重要的海產養殖城市, 為掌握榮成外海的海流特征, 基于2016年1—10月于榮成外海楮島岬角南側布放的浮標觀測數據, 通過譜分析、調和分析等方法, 分析了實測流速時間序列剖面及季節變化, 分離出該觀測點的潮流和余流并進行定性定量分析, 探討了該海域余流的影響因素。結果表明: 1) 榮成外海浮標觀測點流速東西分量總體大于南北分量, 實測流速均值為0.20 m/s, 最大值為0.87 m/s, 實測流速存在明顯的季節差異, 秋季流速最大, 夏季次之, 冬春季較小且冬春季節差異小。實測流存在漲落潮不對稱現象。2) M2分潮絕對占優, 觀測海域潮流性質為規則半日潮流, 潮流以往復流形式運動, 各層潮流均為順時針旋轉。淺水分潮十分顯著, 該點觀測得到的最大潮流流速為0.38 m/s。3) 余流相對較小, 日平均余流速度均值為冬季<春季<夏季<秋季; 流向上E和ENE向占優, 余流總體向東。弱風條件下, 余流流向主要集中于NNE—ESE, 余流流速大于0.15 m/s的概率約為25%; 強風條件下風速和余流流速的互相關系數可達0.95, 夏秋季較強的風速可能是流速增快的原因。

實測流; 調和分析; 潮流; 余流; 榮成外海

淺海水動力過程是海流、波浪、風和地形等多因素綜合作用的結果。狹義的海流定義為海流水平運動分量[1], 包含周期性潮流、由風和地形等因素引起的余流和隨機高頻擾動信號, 潮流影響物質短期輸運, 而余流則影響物質長期輸運過程[2]。海流的影響因素包括風、徑流、海區地形和大小等, 同時養殖活動也是淺海水動力變化的主要原因之一[3], 大量養殖物會消耗海域營養物質, 同時可能造成水交換不及時, 從而影響海域水質環境和養殖效益[4]。因此, 掌握養殖海域自然狀態的海流特征, 對海區生態環境保護、沿海人民生產生活和海洋工程建設等具有重要意義。

榮成市三面環海, 海岸線總長約500 km[5], 有榮成灣、桑溝灣、黑泥灣等多個海灣。榮成是我國重要的海產養殖城市, 其外海水產養殖發達, 貝類、藻類等多種海產品養殖密度高。除了發達的養殖業, 近年來榮成海島核電工程、臨海風力發電項目和濱海旅游業也有所發展[6]。

為加深對榮成外海海域水文環境的認識, 近年來學者們對該海域開展了諸多研究工作。水動力方面, 以風浪為主[7], 冬半年波浪大, 夏季減弱[8]。風向控制的風浪和潮流的相干作用決定了黑泥灣海域懸浮體分布[9]。Fan等[10]在楮島南側海域研究發現往復式潮流導致漲落潮周期內懸浮物質量濃度的不對稱, 大量沉積物輸運出研究區, 有侵蝕事件發生。此外, 山東多處海岸處于侵蝕狀態, 其中, 榮成沿岸屬于嚴重侵蝕型海灘。莊振業等研究表明, 楮島連島沙壩形成是多向泥沙縱向運移的結果[11-12]。黑泥灣海區在風浪較大時, 尤其10 m水深以淺的海域泥沙易起動并隨海流運移。而10 m等深線直逼海岸線, 侵蝕主要發生在岬角處, 侵蝕物堆積在海底或在海灣成為袋狀沙灘[9]。宮立新等[13]指出楮島—鏌铘島砂質海岸的侵蝕速率為每年2～3 m。

關于榮成外海筏式養殖對水動力的影響等方面也有諸多研究。樊星等[14]通過準調和分析法分析桑溝灣潮流流速垂向特征, 并建立了雙阻力模型, 發現養殖區底應力比自然海區高一個量級。嚴立文等[9]通過SWAN模式分析了該海域海帶養殖區域的沉積環境效應, 結果顯示海帶架群消浪系數在0.2～0.4。此外, 榮成外海存在自北向南的沿岸流, 由于該海域有大面積養殖區, 沿岸流主要影響養殖區外圍[15]。綜上所述, 該區域水動力是海流、風、波浪、地形、養殖活動等多重因素綜合作用的結果。基于前人的研究, 對于該海域的水文背景環境等有了一定認識, 然而至今為止的研究主要基于較短的時間觀測序列或模型模擬, 長期觀測資料的匱乏導致目前該海域海流季節變化特征和影響因素等尚未得到分析和解釋。

長時間高質量的實測海流數據能夠在垂向分布和不同時間尺度上多角度反映該海域水動力條件。利用2016年1—10月在榮成外海楮島岬角南側布放的浮標, 基于該長期海洋綜合觀測資料, 統計分析榮成外海海流特征, 分析潮流和余流特征, 并探究海流變化特征與影響因素。

1 數據與預處理

1.1 數據來源

本文采用的數據源自2016年在榮成外海楮島岬角南側(122°33′38″E, 37°01′56″N)布放的1套綜合觀測浮標。圖1(a)中五角星標注處即為觀測站點位置。該浮標站點在黑泥灣養殖區的西北側邊緣外, 與北側楮島岬角直線距離為0.8 km, 距離西側岸線約2.5 km。該觀測點水深在7 m左右, 海底沉積物為砂質。為保證數據質量, 進行了3次浮標的布放和回收, 以完成對儀器的檢修和校準工作。本文數據的3個觀測階段分別為1月22日—2月13日、3月17日—6月20日、7月12日—10月17日。

圖1 榮成外海(a)浮標位置和(b)浮標結構示意圖

如圖1(b)所示, 浮標觀測系統上部為總高度3.5 m的不銹鋼支架, 在高度3.2 m處裝載了Young公司的風向風速儀, 其測量范圍為0～100 m/s, 風速精度為讀數1%, 方向精度為±3°。風數據實時回傳至CR1000數據采集器中并傳回地面基站。在浮標中部, 裝載了2款Nortek公司的聲學多普勒海流剖面儀: 觀測階段1為采樣頻率2 MHz的“闊龍”流速剖面儀(Aquapro Profiler)和觀測階段2、3為采樣頻率600 kHz的“浪龍”流速剖面儀(AWAC), 2款測流儀器的采樣分辨率均為1 mm/s, 精確度均是海流流速測量值的1%。兩個儀器均采用下視工作方式, 換能器位置分別距水面1.5 m和1.2 m。除上述測量儀器, 浮標配置了太陽能供電板用于實時供電, 支架上放置了GPS16X-HVS用于對浮標位置進行實時監控, 浮標系統底部有坐底重石以固定浮標位置。具體參數如表1所示。

表1 2016年觀測起止時間和聲學多普勒海流剖面儀參數設置

1.2 數據預處理

觀測使用的測流儀器在輸出實測數據文件時會進行移動校準, 儀器的俯仰角(pitch)和橫滾角(roll)均在±15°以內。以同一站點同一浮標相同測流儀器觀測到的底層流速數據, 與楮島岬角南側相鄰海床基觀測點S6站(122°33′54″E, 37°01′59″N)同時期使用ADV觀測的流速數據進行對比和相關性分析, 計算得到浮標數據的底層流速與海床基測得的流速東西分量之間的相關系數值為0.91, 證實了浮標系統與數據的可靠性。

觀測使用的聲學多普勒海流剖面儀采取下視工作方式, 可以通過回聲強度特征實現對流速剖面數據的判斷和預處理[16]。由于2016年1—2月數據未全部測量到海底, 因此以回聲強度小于50為判斷閾值; 2016年3月之后剖面數據則通過對6 m以深的相鄰層回聲強度差進行統計, 最終對每月流速剖面數據以相應的第+1層回聲強度減去第層回聲強度大于25或30來對數據進行篩選。

由于觀測過程中可能存在的不穩定性和信號噪聲擾動, 在計算分析前, 基于統計理論的Grubbs準則[17], 對海流、風速等實測數據進行了質量控制, 并對異常值進行去除和線性插值。本文將浮標原位觀測的風速換算至10 m高度風速并進行6 h平均, 隨后與ERA-Interim 10 m風速再分析數據(ERA, 即ECMWF Re-Analysis, 來源于European Center for Medium-Range Weather Forecasts[18])離浮標點經緯度最近的網格進行對比驗證, 經計算, 二者10 m風速相關系數為0.73。浮標數據采集器未能采集到的第三觀測階段9月、10月部分風速數據, 通過ERA- Interim風速再分析數據進行插值補齊。

2 海流特征分析

2.1 海流特征

觀測期間內, 該浮標觀測點流速東西分量變化范圍為–0.72～0.79 m/s(東向為正, 西向為負, 下文同), 流速南北分量的變化范圍為–0.62～0.64 m/s(北向為正, 南向為負, 下文同), 平均流速大小為0.20 m/s, 最大值出現在10月13日, 最大值為0.87 m/s。從圖2(a)和圖2(b)可以看出, 流速東西分量總體大于流速南北分量, 從圖2(c)中可以看出流速大小存在明顯季節變化, 夏秋季節的流速較大。

圖2 榮成外海流速時間序列剖面圖

以1、4、7、10月為代表月(具體時間以1月22日—2月13日代表冬季, 以4月1日—4月30日代表春季, 以7月17日—8月15日代表夏季, 以9月18日—10月17日代表秋季), 對四季實測流速raw進行分級分向統計。結果如圖3所示, 可以看出, 該浮標觀測點流向以E向、W向和WNW方向為主導, 流速主要集中在0～0.3 m/s, 流速流向存在明顯季節變化, 冬季以東西向往復流為主; 春季流速偏小, 流向較分散, 夏秋季漲落潮不對稱現象比較顯著, 東向流速增大。Fan等[10]通過觀測流速和懸浮顆粒物濃度等數據分析研究得到, 榮成近海潮流呈東西向往復運動, 且漲落潮不對稱, 近岸泥沙、海底沉積物被帶出近海區域, 從而導致侵蝕事件的發生, 該觀測點漲落潮不對稱特性與Fan等[10]在鄰近站位的研究結果一致。統計得到, 冬季、春季、夏季、秋季的流速的平均值分別為0.18 m/s、0.18 m/s、0.22 m/s、0.24 m/s, 實測流速均值冬春季<夏季<秋季。

圖3 榮成外海實測流速玫瑰圖

2.2 頻譜分析

為分析海流時間序列的周期特征, 對三個觀測階段實測流速的東西分量、南北分量分別作垂向平均, 其后對三個觀測階段垂向平均流分別作功率譜分析, 譜分析結果如圖4所示, 縱坐標頻譜幅值log10C/(m2?s–2?d) 指示著海流在不同的頻率上的能量。圖4(a)中的譜線相對于圖4(b)、圖4(c)較稀疏, 這是由于觀測階段1僅有23 d, 而觀測階段2和3時間長度近100 d, 因此階段一的譜圖在頻域上的分辨率較低。從圖4可以看出, 流速存在多個顯著周期。浮標觀測站的垂向平均流能量在主要半日分潮M2對應的頻率處最大, 具有較強的半日潮信號, 圖4(b)、圖4(c)中主要半日分潮S2存在譜峰。全日潮信號相對較弱, 全日分潮K1存在譜峰, 全日分潮O1相對不顯著。淺水分潮M4能量值相比全日分潮高2個量級, 淺水分潮MS4和M6能量也十分顯著。浮標觀測點的近慣性振蕩頻率為1.204 5 d–1, 所對應能量沒有出現明顯譜峰, 說明該點近慣性振蕩不顯著。總的來說, 三個觀測時間段的譜線趨勢較為一致, 該觀測點以半日潮為主, 淺水分潮也有顯著貢獻。

3 潮流、余流特征分析

海流是多時間尺度信號共同作用的結果, 通過調和分析、濾波等方法可以將海流分解成潮流和余流等部分。潮流是潮波現象在水平方向的表現形式, 主要影響海域物質在潮周期內的運輸, 長時間序列流速數據能夠得到更準確的回報潮流值。余流是指實測海流中去掉周期性天文分潮后的部分, 余流主要由海面風場、外海環流、河流徑流等驅動。

圖4 實測流垂向平均流速頻譜分析

3.1 潮流特征

采用MATLAB的T_TIDE工具包[19]對2016年1—10月2～6 m水深中各層流速進行調和分析。結果顯示, 通過信噪比檢驗的分潮共39個。統計得到整個觀測期間的潮流東西分量范圍為–0.38～0.35 m/s(以東向為正), 潮流南北分量的范圍為–0.06～0.08 m/s(以北向為正), 潮流速度平均值為0.15 m/s, 潮流速度最大值為0.38 m/s。將四大分潮M2、S2、O1、K1和淺水分潮M4、MS4的潮流橢圓垂向分布繪制于圖5中, 圖中右上角的圖例正圓代表振幅單位大小0.01 m/s, 繪圖程序中預設黑色代表潮流橢圓順時針旋轉, 紅色代表橢圓逆時針旋轉, 由圖中可見, 均為順時針旋轉。從圖5中可以看出, M2分潮絕對占優, 對潮流起主導作用。S2分潮為次顯著分潮, 淺水分潮M4、MS4比全日分潮O1、K1更為顯著。橢圓長軸從海水上層到下層逐漸減小, 潮流會出現越來越扁的情況。潮流在垂向上越接近海底, 通常長半軸越小; 上層橢圓為右旋時, 下層橢圓將逐漸變扁[20], 本文結果與前人觀測到的潮流垂向變化的一般結果相符。

各深度的分潮的長軸、短軸分別列于表2中。如表中所示, M2分潮的最大流速的最大值出現在水深2 m處, 大小為0.220 m/s, 隨深度增加M2分潮最大流速逐漸減小, 垂向上M2分潮最大流速變化量為0.02 m/s。S2分潮最大流速隨深度增加而單調減小, 最大流速約為0.05 m/s。淺水分潮M4最大流速出現在中層, 但整體上差異微小。淺水分潮MS4和全日分潮K1、O1的最大流速相對較小。各分潮短軸垂向上整體差異較小。M2分潮短軸垂向差異小, 計算得到M2分潮橢圓率在–0.064至–0.05。各主要分潮橢圓短軸均為負值, 因此該海域各主要分潮的旋轉方向一致, 均為順時針右旋。方國洪[20]研究認為, 分潮流最大流流向的垂向變化取決于分潮角頻率和科式參量,>時, 變動較小,<時, 隨著接近海底而向左偏轉。M2分潮角頻率大于該觀測點科式參量, 因此分潮流流向的垂向變化較小。近底層水深5.5 m處半日分潮和淺水分潮的調和常數傾角均大于上層水深5 m處和下層水深6 m處的調和常數傾角, 這可能是由于下邊界層底摩擦影響產生。但總體來說, 主要分潮傾角垂向變動很小。

圖5 分潮潮流橢圓垂向分布圖

表2 分潮橢圓長軸、短軸

從圖5可以看出, 主要分潮特征垂向變化較小, 因此可以根據垂向平均流調和分析結果對潮流特征進行定性定量描述[21]。潮流性質的判別可根據式(1)確定, 式中值為各分潮的潮流橢圓長半軸長度。當≤0.5時, 海域潮流性質為規則半日潮流。

根據得到的潮流調和常數計算得到本海域值為0.06, 因此觀測海域潮流性質為規則半日潮流。這一結果與海洋圖集結果相同[22]。根據運動形式, 潮流可以分為往復流和旋轉流, 一般以橢圓率值判斷, 具體公式如式(2)所示, 其中為潮流橢圓短半軸長度,為潮流橢圓長半軸長度。的正負, 即橢圓短軸的正負,值為負時表示潮流橢圓順時針旋轉,值為正時表示潮流橢圓逆時針旋轉[20]。的絕對值越大, 旋轉流性質越顯著, 若的絕對值在0～0.2時, 則判定為往復流。

通過計算得到M2分潮值為–0.06, 即該海域表現為往復流特征。原因可能是觀測站點離海岸和北側楮島岬角近, 受岸線和地形約束。由于觀測站點海域水深較淺, 因此淺水分潮的效應不可忽視。一般以值[公式(3)]來判斷淺水分潮是否顯著, 當值>0.04時, 表明淺水效應顯著。

通過計算得到值為0.15, 大于判別標準近4倍, 即該海域淺水分潮效應顯著。

3.2 余流特征

在余流的計算方面, 前人主要通過實測流速減去調和分析回報潮流, 或進行不同截止頻率的低通濾波。本文實測流速數據為300 s高頻采樣, 在調和分析后去潮余流結果中, 可能存在對風瞬時響應的信號, 而低通濾波則會過濾掉高頻組分。為下一步分析余流的影響因素和對大氣強迫的響應, 選擇實測流經調和分析去掉潮流后得到的余流res結果進行分析。統計得到整個觀測時間段內, 榮成外海浮標觀測點天平均余流東西分量變化范圍為–0.05～0.16 m/s,天平均余流南北分量變化范圍為–0.06～0.09 m/s, 最大值為0.16 m/s, 均值為0.05 m/s, 平均速度對應方向為69°, 呈現往東運移的特征對泥沙物質的長期向外輸運作用不容忽視。

進一步分級分向統計, 得到了該觀測站點的日平均余流季節變化特征, 四季余流玫瑰圖如6所示, 四季均使用同一個區間分級分向。從圖中可以看出, 余流流速流向季節變化均十分明顯, 四季余流流向均以東向或東北向為主, 較大的速度值也出現在落潮流方向。計算得到余流冬季、春季、夏季、秋季的季節平均值分別為0.03 m/s、0.04 m/s、0.07 m/s、0.10 m/s。

圖6 榮成外海余流玫瑰圖

3.3 風對余流的影響

余流特征往往是風場、陸源輸入、地形、人類養殖活動、海岸工程建設等多種因素綜合作用的結果。在淺海區域, 局地風對海流影響顯著[23]。風場的非線性作用對余流有一定影響, 如在潮流相對較弱時, 余流對風響應較為顯著, 余流和風的相關性較好[24]。本文涉及的研究站位水深約7 m, 余流對風的響應更快, 且風的影響可以直至海底。為具體分析風w對余流的影響大小, 首先統計了與流速對應的四季風速風向概率, 結果如圖7所示。從圖7中可以發現, 榮成外海實測風速風向季風特征顯著。通過統計得到, 冬季盛行W、WNW、NW、NNW、N向風, 春季風向分散于16個分方向, 夏季集中于WSW、W、WNW3個方向內, 秋季主要集中于NNE、SSE、S、SSW、NNW方向內。結合風向與余流流向(圖6), 可以發現, 冬、夏季風向和余流流向的主要方向有較好的對應關系。秋季風向和余流流向關系不明顯, 可能受多種因素影響。并且由于秋季原位觀測的風數據缺測而使用ERA數據代替, 盡管ERA數據與原位風速變化趨勢相符, 但仍然不能排除ERA數據在近岸的不確定性。將四季風速分級的概率統計列于表3中, 從表中可以看出, 冬季各級風速概率大致相當, 春、夏季主要集中在0～6 m/s, 秋季大于8 m/s的風速概率顯著增高。

統計榮成外海2016年1—10月浮標觀測站的所有風速風向情況, 以風速大于等于13.9 m/s(風級大于等于7)為強風天氣情況, 風速小于等于3.4 m/s(風級小于等于2)為弱風天氣情況, 得到強風和弱風天氣下10 min平均的風和余流的對比, 如圖8所示。弱風條件下, 余流信號以潮致余流為主, 余流的特征變化較小。共統計弱風天氣累計時長為97 d。弱風風向較為分散, 概率高于10%的風向有WSW、W、WNW。從圖8(c)中可以看出, 在弱風條件下, 余流流向主要集中于NNE—ESE, 余流流速大于0.15 m/s的概率為25%。風速大于13.9 m/s屬于少發天氣事件, 從圖8(d)中可以看出, 余流流向在7級以上的強風事件中, 余流流向發生明顯的變化。

圖7 榮成外海風玫瑰圖

表3 榮成外海四季風速分級概率統計對比表

風速實測數據中, 2月13日強風持續時間最長。如圖9所示, 在2月9日到2月10日風速較小時, 余流基本穩定地向北流動, 在2月10日至2月11日北風增強的過程中, 北向余流有相應的小幅增大。值得注意的是, 2月13日風往南吹且風速逐漸增大, 在2月13日18時之后, 持續保持在風力7級以上, 此時的余流由之前的穩定向北流動轉而向南流動, 但由于和弱風條件下的背景余流方向相反, 因此2月13日的余流速度值并未明顯增大。進一步使用N-CCF函數(the normalised cross-correlation function)[25]計算得到2月9日至2月13日的風速和余流流速的互相關系數, 圖9(c)的軸代表風速和余流流速的超前滯后時間, 當軸值為正時, 代表風速變化超前于余流流速變化。如圖9(c)所示, 兩者的互相關系數最大值為0.68。而僅計算2月13日18時之后的強風事件時間段, 風速和余流流速的互相關系數可達到0.95, 風速和實測流速的相關系數達–0.49, 和余流流速的相關系數為–0.35。從中可以發現, 強風事件對該海域的余流有著顯著的影響, 甚至使余流轉向。

4 結論

基于2016年榮成外海的浮標觀測數據, 本文分析了榮成外海海流特征和影響因素, 得到結論如下:

1) 榮成外海浮標觀測點流速東西分量總體大于南北分量, 實測流速平均流速大小為0.20 m/s, 最大值為0.87 m/s, 實測流速存在明顯的季節差異, 流速均值冬春季<夏季<秋季。該海域存在漲落潮不對稱現象, 指征侵蝕事件。

圖8 (a)弱風, (b)強風, (c)弱風時余流和(d)強風時余流對比玫瑰圖(風向為來向, 流向為去向)

圖9 榮成外海2月9日—13日(a)風、(b)余流及(c)風速和余流速度的互相關系數

注: (a)(b)中黑色加粗線為速度值, 黑色箭頭為速度矢量; (c)中軸代表風速和余流流速的超前滯后時間, 值為正時, 代表風速變化超前于余流流速變化。風矢量和流矢量均為去向

2) M2分潮絕對占優, 觀測海域潮流性質為規則半日潮流, 潮流以往復流形式運動, 各層潮流均為順時針旋轉, 淺水分潮十分顯著。該點觀測得到的最大潮流流速為0.38 m/s。

3) 余流顯著弱于潮流, 天平均余流均值為冬季<春季<夏季<秋季; 流向上E和ENE向占優。天平均余流均值為0.05 m/s, 相對潮流流速較小; 余流平均速度對應方向為69°, 由此看來, 該海域總體呈現往東運移的特征。

4) 弱風條件下, 余流流向主要集中于NNE—ESE, 余流流速大于0.15 m/s的概率為25%。強風條件下風速和流速的互相關系數可達0.95。

[1] 馮士筰, 李鳳歧, 李少菁. 海洋科學導論[M]. 北京: 高等教育出版社, 1999.

FENG Shizuo, LI Fengqi, LI Shaojing. Introduction to ocean science[M]. Beijing: Higher Education Press, 1999.

[2] 曹慧慧, 陳元, 趙騫, 等. 遼東灣東部秋季海流特征的同步觀測研究[J]. 海洋環境科學, 2019, 38(4): 615- 620.

CAO Huihui, CHEN Yuan, ZHAO Qian, et al. Study on the current characteristics in the eastern area of Liao-dong bay in autumn based on synchronous observations[J]. Marine Environmental Science, 2019, 38(4): 615-620.

[3] JACKSON G A, WINANT C D. Effect of a kelp forest on coastal currents[J]. Continental Shelf Research, 1983, 2(1): 75-80.

[4] 蔡繼晗, 李凱, 鄭向勇, 等. 水產養殖對環境的影響及其防治對策分析[J]. 水產養殖, 2010, 31(5): 32-38.

CAI Jihan, LI Kai, ZHENG Xiangyong, et al. The influences of aquaculture on environment and the preven-tion strategy analysis[J]. Journal of Aquaculture, 2010, 31(5): 32-38.

[5] 徐宗軍, 張緒良, 張朝暉. 山東半島和黃河三角洲的海岸侵蝕與防治對策[J]. 科技導報, 2010, 28(10): 90-95.

XU Zongjun, ZHANG Xuliang, ZHANG Zhaohui, et al. Coastal erosion in Shandong Peninsula and Yellow River Delta and related countermeasures[J]. Science & Technology Review, 2010, 28(10): 90-95.

[6] 馬立杰, 崔迎春, 王海榮. 運用因子分析評價山東黑泥灣沉積物重金屬污染來源[J]. 海洋地質前沿, 2012, 28(4): 57-61.

MA Lijie, CUI Yingchun, WANG Hairong. Source of heavy metal contamination in Heini Bay sediments: A case study using factor analysis[J]. Marine Geology Frontiers, 2012, 28(4): 57-61.

[7] 陳吉余. 中國海岸帶地貌[M]. 北京: 海洋出版社, 1996.

CHEN Jiyu. China’s coastal landscape[M]. Beijing: China Ocean Press, 1996.

[8] 劉艷霞, 黃海軍, 樊建勇, 等. 黑泥灣近岸懸沙分布與擴散變化監測[J]. 海洋科學, 2009, 33(1): 25-29.

LIU Yanxia, HUANG Haijun, FAN Jianyong, et al. De-te-ction of distribution change and diffusion of suspended sediment concentration around Heini Bay[J]. Marine Sciences, 2009, 33(1): 25-29.

[9] 嚴立文. 淺海區海帶養殖的沉積環境效應及動力機制[D]. 青島: 中國科學院海洋研究所, 2008.

YAN Liwen. Sedimentary environment evolution in representative kelp (laminaria japonica)-cultured region (Harny Bay) and inner hydrodynamic mechanism[D]. Qingdao: Institute of oceanology, Chinese academy of science, 2008.

[10] FAN R F, WEI H, ZHAO L, et al. Identify the impacts of waves and tides to coastal suspended sediment concentration based on high-frequency acoustic observations[J]. Marine Geology, 2019, 408: 154-164.

[11] 莊振業, 劉冬雁, 劉承德, 等. 海岸帶地貌調查與制圖[J]. 海洋地質動態, 2008, 24(9): 25-32.

ZHUANG Zhenye, LIU Dongyan, LIU Chengde, et al. Geo--morphological survey and mapping of coastal zone[J]. Marine Geology Letters, 2008, 24(9): 25-32.

[12] 莊振業, 曹立華, 李兵, 等. 我國海灘養護現狀[J]. 海洋地質與第四紀地質, 2011, 31(3): 133-139.

ZHUANG Zhenye, CAO Lihua, LI Bing, et al. The cur-rent state of beach maintenance in our country[J]. Marine Geology & Quaternary Geology, 2011, 31(3): 133-139.

[13] 宮立新. 山東半島東部海灘侵蝕現狀與保護研究[D]. 青島: 中國海洋大學, 2014.

GONG Lixin. erosion situation and protection of beach in eastern part of shangdong peninsula[D]. Qingdao: Ocean University of China, 2014.

[14] 樊星. 典型養殖海區潮動力結構特征的初步研究——觀測與數值模擬[D]. 青島: 中國海洋大學, 2008.

FAN Xing. Preliminary studies on the features of tidal- dynamic structure in a typically high density mariculture coastal bay-observation and simulations[D]. Qing-dao: Ocean University of China, 2008.

[15] 朱述龍, 張占睦. 遙感圖像獲取與分析[M]. 北京: 科學出版社, 2002.

ZHU Shulong, ZHANG Zhanmu. Remote sensing image acquisition and analysis[M]. Beijing: Science Press, 2002.

[16] 任少華, 周沖, 葉小凡. 海流測量數據的后處理探究[J]. 海洋技術, 2013, 32(4): 63-66.

REN Shaohua, ZHOU Chong, YE Xiaofan. Research on the post-processing of current velocity data[J]. Journal of Ocean Technology, 2013, 32(4): 63-66.

[17] 劉首華, 陳滿春, 董明媚, 等. 一種實用海洋浮標數據異常值質控方法[J]. 海洋通報, 2016, 35(3): 264- 270.

LIU Shouhua, CHEN Manchun, DONG Mingmei, et al. Research on the Post-Processing of Current Velocity Data[J]. Marine Science Bulletin, 2016, 35(3): 264- 270.

[18] SONG L N, LIU Z L, WANG F. Comparison of wind data from ERA-Interim and buoys in the Yellow and East China Seas[J]. Chinese Journal of Oceanology and Limnology, 2015, 33(1): 282-288.

[19] PAWLOWICZ R, BEARDSLEY B, LENTZ S. Classical tidal harmonic analysis including error estimates in MATLAB using T_TIDE[J]. Computers & Geosciences, 2002, 28(8): 929-937.

[20] 方國洪. 潮流垂直結構的基本特征——理論和觀測的比較[J]. 海洋科學, 1984, 8(3): 1-11.

FANG Guohong. Basic characteristics of vertical stru-c-ture of tidal current-comparison between theory and observation[J]. Marine Sciences, 1984, 8(3): 1-11.

[21] 中交第一航務工程勘察設計院有限公司. JTS 145— 2015, 港口與航道水文規范[S]. 北京: 人民交通出版社, 2015: 107-112.

CCC First Harboar Consultants Co Ltd (FDINE). JTS 145—2015, Code of Hydrology for Harbour and Waterway[S]. Beijing: China Communications Press, 2015: 107-112.

[22] 海洋水文圖集編輯委員會. 渤海黃海東海海洋圖集: 水文[M]. 北京: 海洋出版社, 1993.

Editorial Board for Marine Atlas of Bohai Sea Yellow Sea East China Sea. Marine Atlas of Bohai Sea Yellow Sea East China Sea: Hydrology[M]. Beijing: China Ocean Press, 1993.

[23] 趙保仁. 局地風對黃海和東海近岸淺海海流影響的研究[J]. 海洋與湖沼, 1982(6): 479-490.

ZHAO Baoren. A study of the effects of local wind on currents in the off-store areas of the east China sea and huanghai Sea[J]. Oceanologia et Limnologia Sinica, 1982(6): 479-490.

[24] 胡鋒濤, 陶宇. 浙中沿岸流及其對風和潮汛的響應[J]. 海洋科學, 2019, 43(12): 50-58.

HU Fengtao, TAO Yu. Coastal current and its response to wind and spring—Neap cycles in the middle Zhejiang coastal area[J]. Marine Sciences, 2019, 43(12): 50-58.

[25] TSAI D M, LIN C T, CHEN J F. The evaluation of nor-ma-lized cross correlations for defect detection[J]. Pattern Recognition Letters, 2003, 24(15): 2525-2535.

Characteristics and influencing factors of the currents in the Rongcheng offshore region

WEI Shi-yan1, YANG Wei2, ZHAO Liang1

(1. College of Marine and Environmental Science, Tianjin University of Science and Technology, Tianjin 300457, China; 2. School of Marine Science and Technology, Tianjin University, Tianjin 300072, China)

Rongcheng is an important mariculture city in China. To determine the characteristics of ocean currents in the Rongcheng offshore region, we investigated the features of the measured currents using spectral and harmonic analysis methods based on the in situ observations obtained from a moored buoy south of the Chudao Island from January to October 2016. The tidal and residual currents are analyzed separately. The influencing factors of the residual currents are highlighted. The main conclusions are as follows. (1) The zonal component of the velocity generally has larger amplitudes than the meridional component. The average and maximum values of the measured velocity are 0.20 and 0.87 m/s, respectively. In addition, the velocity of the currents shows strong seasonal variation; the velocity is the strongest in autumn, followed by summer and weakest in winter and spring. Significant asymmetries between flood and ebb flow are observed. (2) The observed currents are dominated by the M2tidal component and the shallow water tide is significant. Further analysis shows that the dominant tide belongs to a type of regular semidiurnal tidal current that rotates clockwise in all layers. The observed tidal current velocity is 0.38 m/s. (3) The residual current is relatively small compared to the tidal component. The daily averaged residual flow is the largest and smallest in autumn and winter, respectively. E and ENE are dominant, and the residual current generally exhibits the characteristics of eastward migration. Under the condition of weak wind, the residual current is primarily in the NNE–ESE direction. Approximately 25% of residual flow velocity is greater than 0.15 m/s. The correlation coefficient between the wind speed and velocity under strong wind conditions can reach 0.95. The strong wind speeds in summer and autumn may account for the increase of the residual velocity during that time.

observed current; tidal harmonic analysis; tidal current; residual current; Rongcheng offshore region

Apr. 12, 2021

P731

A

1000-3096(2022)04-0055-12

10.11759/hykx20210412003

2021-04-12;

2021-05-14

國家自然科學基金項目(41876018); 天津市自然科學基金項目(19JCZDJC40600)

[Project of National Natural Science Foundation of China, No. 41876018; Tianjin Natural Science Foundation project, No. 19JCZDJC40600]

魏詩晏 (1995—), 女, 江蘇鹽城人, 碩士, 主要從事淺海動力學研究, E-mail: 1254510293@qq.com; 趙亮(1975—),通信作者, 教授, 主要從事淺海動力學、海洋生態動力學研究, E-mail: zhaoliang@tust.edu.cn

(本文編輯: 叢培秀)