基于潮汐逆模型技術對渤黃海正壓M2分潮開邊界條件的優化研究：II.潮汐特征、潮汐動力學及潮余流

宋軍，姚志剛，郭俊如，李靜，高佳，董軍興

（1.國家海洋信息中心，天津　300171；2.中國海洋大學　教育部物理海洋學重點實驗室，山東　青島　266100；3.國家海洋局海洋減災中心，北京　100194）

基于潮汐逆模型技術對渤黃海正壓M2分潮開邊界條件的優化研究：II.潮汐特征、潮汐動力學及潮余流

宋軍1，2，姚志剛2，郭俊如3，李靜1，4，高佳1，董軍興1

（1.國家海洋信息中心，天津300171；2.中國海洋大學教育部物理海洋學重點實驗室，山東青島266100；3.國家海洋局海洋減災中心，北京100194）

基于該系列文章前文建立的高分辨率數據同化模型系統，研究渤、黃海的M2分潮的潮汐特征、潮混合、潮余流、潮能及其擴散、潮動量平衡。對同潮圖結果的分析表明，渤、黃海內共存在4個無潮點，研究區域內M2分潮振幅的最大值出現在朝鮮半島西側。潮流橢圓的分析結果表明，在35°N附近區域潮流橢圓的旋轉方向有著比較突然的改變。朝鮮半島西側區域潮流的振幅較大，并且由于該區域地形復雜多變導致朝鮮半島西側海域存在著比較強烈的潮混合與潮耗散。對M2分潮動量分析的結果表明，地形等因素在近海潮汐動力學過程中發揮著比較重要的作用。在渤、黃海內大部分區域，主要的動量平衡基本上介于壓強梯度力、局地變化項以及科氏力之間，底應力和對流項都可以忽略不計。但是，線性模型對近海尤其是上述強耗散區內潮波活動的模擬還是存在著缺陷和不足。

潮汐特征；潮汐動力學；潮余流

在渤、黃海海區的眾多水文現象之中，潮汐和潮流過程以其顯著的動量、能量輸運而成為該海區最引人注目的特征之一。已有的研究已經指出渤、黃海海區的潮汐和潮流過程具有較復雜的空間結構特征（Choi，1980；Choi，2003；Fang et al，2004；Larsenetal，1985；王凱等，1999；趙保仁，1994）。此外，渤、黃海海區較淺的地理特征以及顯著的潮振幅也使得渤黃海成為世界上潮能耗散最強的海區之一，Provost等（1997）曾指出，主要太陰半日分潮M2在該海區的潮能耗散達到了180 GW，約占全球M2潮能耗散總能量（2 TW）的11.1%。因此，潮混合以及潮余流過程在該海區的動力學以及其生物、化學過程中起著重要甚至決定性作用。例如，潮混合在黃海冷水團的形成、渤海鋒面結構、鳀魚卵分布以及細菌生長等方面均起著至關重要的作用（呂新剛 等，2010；Lie，1989；Hao et al，2003；Hyun et al，1999）。因此，對潮汐潮流過程的研究成為該海區物理海洋學研究的重點之一。

1　潮汐特征

經過數據同化系統的優化調整，ROMS模式對渤、黃海區域M2分潮的模擬與觀測之間的誤差已經顯著地降低，兩者之間達到了比較不錯的符合程度。因此，本文可以在ROMS模式“最優”解的基礎上進一步的分析渤、黃海區域M2分潮的潮汐特征。從圖1的M2同潮圖可以看出，由南邊界傳入的M2潮波分成兩支向北傳播：一支折向山東半島并遭到山東半島南岸的反射；另外一支繼續向北進入北黃海后遭到黃海北岸岸線的反射，入射波與反射波的相互干擾結果，在科氏力作用下分別形成了兩個逆時針旋轉的駐波系統，無潮點的位置分別位于山東半島成山角外側和海州灣外。進入北黃海潮波的一部分繼續折向西并經過渤海海峽進入渤海，同樣受到渤海西岸的阻擋，反射波作用下分別在遼東灣西側和渤海灣南側，形成了兩個逆時針旋轉的潮波系統，無潮點的位置分別位于秦皇島和黃河口附近，渤海內的兩個無潮點的位置非常接近陸地，表現為準退化的潮波系統。

圖1　同潮圖

應當指出的是，目前關于渤海灣內無潮點的位置尚存在著爭議。Fang等（2004）的結果認為渤海灣內無潮點的位置應當在陸地之上，表現為退化潮波系統。而王永剛等（2004）基于趨近法的結果與本文的模式結果非常類似，即渤海灣內無潮點的位置非常接近陸地。造成這一分歧的一個主要原因是渤海灣內的無潮點非常接近于黃河口三角洲附近區域，由于黃河的巨大輸沙量導致該區域的地形和岸線處在不斷的變化之中，使得對于黃河三角洲附近區域潮汐活動的觀測和模擬都較為困難。在本文的研究區域內，M2分潮的最大振幅出現在朝鮮半島西側的江華灣灣頂，最大振幅超過了2.6 m，整個朝鮮半島西側海域是世界上潮差最為顯著的區域之一。繼續向西至靠近中國大陸，M2分潮的振幅顯著減小，海州灣頂M2分潮的最大振幅僅有1.4 m左右。北部的渤海海區由于水深較淺，從而底摩擦相對比較強導致渤海區域內M2的潮振幅相對黃海偏弱。

圖2　深度平均的M2分潮潮流橢圓（每10個格點一個潮汐橢圓，藍色（紅色）代表反時針（順時針）旋轉）

同時，還給出了對渤、黃海區域內對深度平均的M2分潮潮流橢圓的分析結果，如圖2所示。在渤海，M2分潮流的最大流速在20～60 cm/s之間，渤海灣、遼東灣流速相對較強，萊州灣流速較弱，其大小僅有10～20 cm/s。渤海海峽附近區域M2分潮流最大流速可以達到40 cm/s。對潮流橢圓扁率的分析表明，渤海中部區域的潮流更多的表現為旋轉流，而遼東灣和渤海灣內的潮流更多表現為往復流的形式。在黃海，M2分潮流一個比較明顯的特征是流速東西方向的不對稱，即東岸流速大于西岸流速。西岸流速在20～40 cm/s之間，東岸流速大小介于60～80 cm/s，尤其是江華灣內，最大流速可以達到100 cm/s，為黃海區域內的最強潮流區。黃海中部區域潮流較弱，大小僅有20 cm/s左右。對黃海內潮流橢圓扁率的分析表明，朝鮮半島西側海域M2分潮流更多的表現為往復流的形式，而在靠近中國大陸一側，M2分潮流則更多的表現為旋轉流的特征。此外，M2潮流橢圓一個比較顯著的特征是，在35每10個格點一個潮汐橢圓，藍色（紅色）代表反時針（順時針）旋轉N以南區域M2分潮流的旋轉方向基本上為順指針，而在35°N以北，除山東半島南端的一小部分區域之外，M2分潮流基本變為逆時針旋轉，這一發現與Lie等（2002）之前的觀測以及Kang等（1998）的模式結果基本一致。山東半島南端一小部分區域內M2分潮流的順時針旋轉結果也和之前的研究（Lie et al，2002；Kang et al，1998；Lee et al，1999）相一致。在渤海，M2分潮流順時針旋轉的區域主要出現在渤海中部，在其他區域基本上為逆時針旋轉。注意到M2分潮流旋轉方向的之一空間分布特征與之前Prandle（1982）的關于分潮流旋轉方向的理論分析相符，即在分潮等振幅線間隔比較均勻的區域，分潮流的旋轉方向為逆時針，本文中北黃海以及山東半島南端M2等振幅線和分潮流旋轉方向之間的關系基本與Prandle的結果相一致。

2　潮汐動力學及潮余流

為了解渤、黃海海區M2分潮對水動力學環境的影響，本節對與M2分潮相關的潮混合、潮余流、潮能通量以及潮能耗散等進行了分析研究，最后還對渤、黃海海區內M2分潮的動量平衡進行了診斷分析。

2.1潮混合

在對渤、黃海潮汐潮流動力學作用的研究中，一個比較重要的方面就是研究海區內強潮流所導致的潮致湍混合對鋒面的形成、結構及其分布的影響。Simpson等（1974）對愛爾蘭海域的潮混合鋒面進行研究時提出可以用U3/H作為確定潮混合鋒面位置的依據，其中U為深度平均的M2分潮潮流流速，H代表局地水深。基于Simpson等提出的方法，計算了渤、黃海區域內潮混合參數（Tidal Mixing Parameter，以下簡稱TMP）：

TMP=log10（h/U3）

潮混合參數TMP的空間分布結果如圖3所示，可以看到存在4個比較明顯的潮混合參數TMP的低值區域，分別位于：西朝鮮灣、江華灣、朝鮮半島西南部附近海域以及34°N-35°N之間的江蘇外海。Lie等（2002）和Naimie等（2001）的研究指出潮混合參數TMP=2.0處對應為潮混合鋒面發生時的位置（即圖3中的白色等值線），從圖3中可以看到，在本文的研究中渤、黃海內存在著4個比較明顯的潮混合鋒面，分別位于：西朝鮮灣外、江華灣外、朝鮮半島西南部以及江蘇外海。這一發現與之前的衛星觀測的結果（Lie，1989；Hickox et al，2000）基本一致。除強潮流所導致的潮致混合，淡水注入以及海水的層結等浮力通量的變化也對潮混合鋒面的形成、強度以及位置有著重要的影響。例如Lie等（1989）曾經指出長江的巨大淡水通量強化了江蘇外海附近的海水層結，從而抑制和影響了江蘇外海潮混合鋒面的形成。因此，受到以上因素的影響，潮混合鋒面的位置、形狀、結構等隨時間存在較明顯的變化。

圖3　M2分潮的潮混合參數的空間分布

2.2潮余流

由于摩擦、海底地形變化、邊界形狀等原因，使得潮流發生非線性現象，從而水質點在一個潮周期之后并不返回到原先的位置 （Zimmerman，1978），稱這種由于潮流非線性相互作用所導致的流動為潮余流。潮余流在近海以及河口區域對懸浮物、泥沙以及其他物質的長期輸運有著重要影響，因此在本節對渤、黃海海區潮余流的空間分布特征進行分析。由于本文的研究中僅包括M2一個天文分潮，因此通過對ROMS模式輸出在12.42 h（M2天文分潮的周期）上進行平均，可以得到M2分潮的潮余流，一般稱這種在固定空間點上將潮流在一個潮周期T內時間平均所獲得的余流為歐拉余流。如圖4所示，給出了渤、黃海海區歐拉余流的空間分布，從圖中可以看出，在渤海和黃海內部，潮余流的強度比較弱，最大為2 cm/s左右，而且分布比較有規律；而近岸附近的潮余流量值較大，最大量值超過了5 cm/s并且分布比較凌亂。在渤海，潮余流較強的區域出現在黃河口、渤海海峽北部的老鐵山水道以及遼東灣內，這與Fang等（1985）的結果相一致；黃海的余流較渤海略大，強余流出現的地點分別為西朝鮮灣、江華灣等朝鮮半島西側海區以及江蘇外海，Lee等（1999）的結果也給出了類似的潮流分布特征。在遼東半島最南端我們可以看到存在著一對氣旋和反氣旋環流，這與之前黃大吉等的研究結果基本一致，導致這一現象出現的原因可能是該海域漲、落潮過程的不對稱（Huang et al，1999）。此外，在萊州灣內M2潮余流存在一個反氣旋的渦旋結構，但強度較弱。繼續向南，我們看到黃海內部區域余流的方向基本偏南，余流的大小偏弱；在南黃海靠近中國大陸一側，余流沿山東半島南下，之后在江蘇外海折向東從而形成一個氣旋結構的渦旋。

圖4　深度平均的M2潮余流（每5個格點一個潮流矢量）

2.3潮能及其耗散

之前對于渤、黃海海區的研究曾經指出，全球的M2潮汐能量的大約8%耗散在黃海區域（Provost et al，1997；Kantha et al，1995；Munk，1997），因此，對渤、黃海正壓潮能通量及其耗散的研究室該區域潮汐研究的重要方向之一，本文基于ROMS模式的結果，給出了渤黃海潮汐通量及其耗散的研究結果。根據Greenberg（1979），平均潮能通量輸運，包括潮動能和潮勢能，可以定義為：

上式中，ρ0為海水密度，取為1 026 kg/m3；h為局地水深，單位為m；η代表由于潮汐過程導致的海表面高度的漲落；u和v分別是深度平均的潮流流速；g為重力加速度，一般取為9.8 m/s2。其中，海表面高度η、深度平均的流速u、v分別來自于優化后的ROMS模式輸出。T表示一個M2分潮的潮周期，為12.42 h。因此，最終得到的潮能通量F是在一個潮周期上的平均值，潮能通量F的空間分布如圖5所示。從圖5可以看出，進入渤海和黃海的潮能有兩個主要來源：大部分潮能來自于東海，經模式的南邊界進入黃海，還有一小部分能量來自于東側的對馬海峽。

圖5　M2分潮的潮能通量（每5個格點一個矢量）

渤、黃海潮能通量空間分布最顯著的特征是在黃海潮能通量的東西不對稱，東岸靠近中國大陸的一側潮能通量比較小，而臨近朝鮮半島的西側海區潮能通量比較強。這與前面提到的潮流橢圓的空間分布特征是基本一致的。對潮流能量的輸運路徑分析表明，進入南黃海的M2潮流通量主要沿朝鮮半島一側向北輸運，在36°N附近海區分為兩支：其中一小部分向西折向山東半島南岸，之后沿山東半島及江蘇沿岸南下進入蘇北外海，形成一個巨大的渦旋結構；大部分潮能通量則繼續北上，進入北黃海的潮能通量在到達黃海北岸后折向西沿岸線流動，并與山東半島北岸的東向潮能輸運一起形成氣旋性的渦旋。進入北黃海潮能通量的一小部分繼續向西經渤海海峽（主要沿老鐵山水道）而進入渤海，之后一部分折向北進入遼東灣，其他的則繼續向西進入渤海中部。從渤、黃海潮能通量輸運的路徑以及強度變化可以看出，潮能通量的大部分沿朝鮮半島西側輸運的過程中分別在朝鮮半島西南部海域、江華灣以及西朝鮮灣附近匯聚并最終在局地耗散掉，因此這3個區域是渤黃區域海潮能能量最主要的“匯”。

根據動力學的基本定理，在數學上功率可以表達為力與速度的乘積。因此由底摩擦所導致的潮能耗散率可以定義為（Greenberg，1979；Taylor，1919）：

其中，Cd為底摩擦系數，在本文的研究中取為0.003；ρ0為海水密度，取為1 026 kg/m3；u和v分別是深度平均的潮流流速；T表示一個M2分潮的潮周期，為12.42 h。這里對潮能耗散率的計算基于二次的底應力公式，這與本文數據同化系統中ROMS模型底應力的計算公式是一致的。渤、黃海整潮周期平均潮能耗散的空間分布如圖6所示，最大的潮能耗散出現在朝鮮半島西側海域，這與前面對潮能通量輸運的分析相一致。在以上這些區域，水深相對較淺，地形變化比較復雜，從而導致這些區域的底粗糙度較大，此外，沿朝鮮半島西側海域的M2潮流流速振幅也比較顯著，因此導致了在這些區域M2分潮的潮能耗散非常驚人，這也與之前對該區域潮能耗散的研究相吻合（Lefevre et al，2000；Kang et al，2002）。以log10D=0.3 W/m2作為判定強潮能耗散的臨界值（Naimie et al，2001），在本文的研究區域內可以發現3個潮能耗散非常強的區域，分別位于：西朝鮮灣、江華灣以及朝鮮半島西南部海域，其中，潮能耗散最強的區域出現在江華灣內部，潮能耗散的最大值高達10 W/m2。此外，在江蘇外海的蘇北淺灘附近，雖然潮能耗散的強度較朝鮮半島西側海區偏弱，但也屬于潮能耗散較顯著的區域。以上4個強潮耗散區所具有的一個共同的特征是其均具有非常強的潮流和較粗糙的海底地形，從而導致了M2分潮在這些區域內較強的耗散。在渤、黃海的內部區域，潮能耗散非常微弱，與以上幾個區域相比基本可以忽略。

圖6　M2平均潮能耗散速率（單位：W/m2，以log10為底）

2.4動量分析

為了對渤、黃海區域潮汐動力學過程有更加清晰的認識，在本節對渤、黃海區域的潮汐過程進行動量分析和診斷。基于ROMS優化之后的模擬結果，分別對深度積分的動量方程中各項的量值進行了計算，ROMS模型分別輸出并保存了動量方程中各項的u分量和v分量，并在一個M2分潮的潮周期（12.42 h）上分別計算了各項的均方差。最后基于其中5個比較重要的項給出了對動量的結果分析,這五項分別是：壓強梯度力 （pressure gradient）、科氏力（corriolis）、底應力 （bottom friction）、對流項（advection）和局地變化項（local acceleration），結果如圖7所示。從對動量方程各項均方差的對比可以看出，渤、黃海區域最顯著的兩項分別是壓強梯度力和局地變化項，科氏力的均方差相對前兩項偏弱，位居第3。科氏力的大小雖然較前兩項偏弱，但是其具有和壓強梯度力相同的量級，此外科氏力也在渤、黃海旋轉潮波系統的形成中起著重要作用。水平對流項和底應力項在絕大部分研究區域內都非常弱，較壓強梯度力等前三項至少小一個量級，這說明在渤、黃海的大部分區域，潮波系統的非線性作用非常微弱，潮波基本表現為線性過程。但是在一些沿岸地區，例如西朝鮮灣，江華灣等，可以看到底應力等非線性項的量級顯著增加到和壓強梯度力相同的大小，甚至基本的動量平衡變為壓強梯度力和底應力之間的平衡。

圖7　一個M2潮周期內，基于深度積分動量方程的U、V分量，計算的壓強梯度力、科氏力、底應力、對流項和局地變化項的均方差

基于動量方程中各項的大小以及相對權重的變化，可以把渤、黃海研究區域劃分成3個區域：弱潮耗散區、強潮耗散區以及中間過渡區。為對動量方程各項之間各平衡關系及時間演變有一個更直觀的了解，在以上3個區域中分別給出了一個典型站位的動量分析時間序列，3個典型站位的位置如圖8所示。其中，站位A位于山東半島以南，接近南黃海的中部40 m等深線位置附近；站位B位于渤海灣中部，該站位的水深大約為10 m左右；站位C位于朝鮮半島西側的江華灣內，該站位水深較淺僅有6 m左右。結果如圖9所示，給出了以上3個站位共25 h的動量方程各項的時間序列圖（分別包括U分量和V分量）。

圖8　典型站位的位置

圖9　基于深度積分動量方程的U、V分量A，B，C三站M2分潮的壓強梯度力、科氏力、底應力、對流項和局地變化項對比。

在南黃海內的A站位，可以看到其動量平衡關系與渤黃海絕大部分地點類似。在東西方向上，主要的動量平衡介于壓強梯度力、局地變化項以及科氏力之間，其中壓強梯度力和局地變化項的大小基本相同，位相相反，科氏力較局地變化項明顯偏弱。底應力和水平對流項在A站位的東西方上基本上可以忽略不計。在該站位的南北向分量情況基本與東西向分量類似，但南北向的科氏力大小較之東西向有比較顯著的增加，其大小和位相基本上與局地變化項相同，科氏力和局地變化項共同作用來平衡掉壓強梯度力。這說明在A站位，潮流主要為為東西向流動，這與前面對該區域潮流橢圓的分析相符合。在南黃海內的A站位，東西向和南北向動量分析的都表明底應力和水平對流項都可以忽略不計，在A站位的潮流運動主要表現為線性。

向北進入渤海灣內，可以看到在B站位情況與A站位情況略有不同，底應力的大小及作用在B站位附近海區顯著增加。在東西方向上，動量方程中最基本的平衡與A站位相同，依然介于壓強梯度力和局地變化項之間；不同的是在B站位，底應力的大小已經超過了科氏力，位居第3位。科氏力的作用變得非常微弱，與水平對流項具有相同的量級，其作用基本可以忽略。在南北方向，B站位的動量平衡介于壓強梯度力和科氏力之間，即基本的地轉平衡。而局地變化項的作用在該方向上變得比較弱。底應力與水平對流項在B站位的南北向上依然可以忽略不計。這與前面對該區域潮流橢圓的分析結果一致，在渤海灣內，潮流的主要流動方向為東西向，導致南北向的底應力基本為零，因此在南北方向上基本的平衡為地轉平衡。在渤海灣內的B站位，東西向的底應力作用的顯著增加表明，在該站位的潮流運動已經表現出相當程度的非線性。

江華灣內的動量平衡與以上兩個站位相比截然不同。從前文對潮能通量等的分析可知，在朝鮮半島西側的潮汐能量非常巨大，因此在江華灣內的C站位動量分析各因子的量級較前面的A和B兩個站位顯著增加，增大到B站位的4倍大小左右。C站位動量平衡關系最為顯著的特征即是底應力和水平對流項在該站附近增加到和其他幾項相同的量級。在東西向，動量方位潮汐-地形間非線性相互作用導致了淺水分潮M4的發生。此外，還注意到底應力與壓強梯度力基本反位相，并且兩者之間存在一個小時的延遲，這也說明了在該站位潮汐非線性作用的重要性。考慮到江華灣內島嶼眾多，地形復雜以及自南海進入該海域的巨大潮汐能量，潮汐非線性作用在該海域的顯著增加是比較自然的事情。在C站位，其南北向的動量分析結果與東西方向類似。在C站位，淺水分潮的形成也說明了在該區域非線性作用的顯著增強，這與我們前面在潮混合和潮能耗散的分析結果相一致。

3　結論

本系列文章上一篇中建立了一套高精度的數據同化系統用來對正壓潮汐開邊界條件進行優化。該數據同化系統主要包括一個正模型ROMS和逆模塊TRUXTON，基于ROMS模式結果與觀測數據之間的誤差，驅動TRUXTON模型產生對潮汐開邊界條件的訂正項，從而優化原潮汐開邊界條件。

ROMS模式調整前和調整后的同潮圖結果對比表明，優化后的M2分潮同潮圖更加接近于Fang等（2004）所給出的“真實值”，其中最大的改善出現在渤海、渤海海峽以及北黃海附近區域。對同潮圖結果的分析表明，渤、黃海內共存在4個無潮點，分別位于山東半島成山角外側、海州灣、秦皇島外海以及舊黃河口附近；研究區域內M2分潮振幅的最大值出現在朝鮮半島西側的江華灣內，最大值達到了2.6 m。潮流橢圓的分析結果表明，在35°N附近區域潮流橢圓的旋轉方向有著比較突然的改變。朝鮮半島西側區域潮流的振幅較大，并且由于該區域地形復雜多變導致朝鮮半島西側海域存在著比較強烈的潮混合與潮耗散，尤其是在江華灣、西朝鮮灣以及朝鮮半島西南部海域。對M2分潮動量分析的結果表明，地形等因素在近海潮汐動力學過程中發揮著比較重要的作用。在渤、黃海內大部分區域，主要的動量平衡基本上介于壓強梯度力、局地變化項以及科氏力之間，底應力和對流項在大部分區域內都可以忽略不計，潮流基本上表現為線性。

但在西朝鮮灣，江華灣等區域，潮流流速的顯著增加以及海底地形的復雜變化使得在這些區域內底應力的大小和重要性顯著增加，從而在上述區域基本的動量平衡介于壓強梯度力和底應力之間，對江華灣內C站位動量的時間序列分析可以看出，潮波-地形間非線性相互作用導致了淺水分潮M4的形成，因此線性模型對近海尤其是上述強耗散區內潮波活動的模擬還是存在著缺陷和不足。

Choi B H,1980.A tidal model of the Yellow Sea and the Eastern China Sea,Korea Ocean Research and Development Institute,Seoul.

Choi B H,2003.A synchronously coupled tide-wave-surge model of the Yellow Sea,Coast Eng,47(4):381-398.

Fang G,Yang J,1985.A two-dimensional numerical model of the tidal motions in the Bohai Sea.Chinese Journal of Oceanology and Limnology,3(2):135-152.

Fang G,2004.Empirical cotidal charts of the Bohai,Yellow,and East China Seas from 10 years of TOPEX/Poseidonaltimetry.J GEOPHYS RES-OCEANS,109(C11):C11006.

Greenberg D A,1979.A numerical model investigation of tidal phenomena in the Bay of Fundy and Gulf of Maine.MAR GEODESY,2(2):161-187.

Hao W,2003.Tidal front and the convergence of anchovy(Engraulis japonicus)eggs in the Yellow Sea.Fish Oceanogr,12(45):434-442.

Hickox R,2000.Climatology and seasonal variability of ocean fronts in the East China,Yellow and Bohai Seas from satellite SST data. Geophys Res Lett,27(18):2945-2948.

Huang D,1999.Modelling the seasonal thermal stratification and baroclinic circulation in the Bohai Sea.Cont Shelf Res,19(11): 1485-1505.

Hyun J H,1999.Tidally induced changes in bacterial growth and viability in the macrotidal Han River estuary.Yellow Sea,Estuarine Coastal Shelf Sci,48(2):143-153.

Kang S K,1998.Fine grid tidal modeling of the Yellow and East China Seas.Cont Shelf Res,18(7):739-772.

Kang S K,2002.Two-layer tidal modeling of the Yellow and East China Seas with application to seasonal variability of the M2tide.J GEOPHYS RES,107(C3):3020.

Kantha L H,1995.Barotropic tides in the global oceans from a nonlinear tidalmodelassimilatingaltimetrictides:2.Altimetricand geophysical implications J Geophys Res,100(25):309-325.

Larsen L H,1985.East China Sea tide currents.Cont Shelf Res,4(1-2):77-103.

Lee J C,K T Jung,1999.Application of eddy viscosity closure models for the M2 tide and tidal currents in the Yellow Sea and the East China Sea.Cont Shelf Res,19(4):445-475.

Lee S H,R C Beardsley,1999.Influence of stratification on residual tidal currents in the Yellow Sea.J GEOPHYS RES-OCEANS,104 (C7):15:615-679,701.

Lefevre F,2000.How can we improve a global ocean tide model at a regional scale?A test on the Yellow Sea and the East China Sea.J GEOPHYS RES,105(C4):8707-8725.

Lie H J,1989.Tidal fronts in the southeastern Hwanghae(Yellow Sea). Cont Shelf Res,9(6):527-546.

Lie H J,2002.Computation methods of major tidal currents from satellite-tracked drifter positions,with application to the Yellow and East China Seas.J GEOPHYS RES-OCEANS,107(C1): 3003.

Munk W,1997.Once again:once again-tidal friction.Prog Oceanogr,40 (1-4):7-35.

Naimie C E,2001.Seasonal Mean Circulation in the Yellow Sea-A Model-Generated Climatology.Cont Shelf Res,21(6-7):667-695.

Naimie C E,2001.Seasonal Mean Circulation in the Yellow Sea-A Model-Generated Climatology.Cont Shelf Res,21(6-7):667-695.

Prandle D,1982.The vertical structure of tidal currents.Geophysical& Astrophysical Fluid Dynamics,22(1-2):29-49.

Provost C L,F Lyard,1997.Energetics of the M2 barotropic ocean tides: an estimate of bottom friction dissipation from a hydrodynamic model.Prog Oceanogr,40(1-4):37-52.

Simpson J H,J R Hunter,1974.Fronts in the Irish Sea.Nature,250:404-406.

Taylor G I,1919.Tidal friction in the Irish Sea.Proceedings of the Royal Society of London.Series A,96(678):330-330.

Zimmerman J,1978.Topographic generation of residual circulation by oscillatory(tidal）currents.Geophysical&Astrophysical Fluid Dynamics,11(1):35-47.

呂新剛，趙昌，夏長水，等，2010.膠州灣水交換及灣口潮余流特征的數值研究.海洋學報，32（02）:20-30.

王凱，方國洪，馮士，1999.渤海、黃海、東海M_2潮汐潮流的三維數值模擬.海洋學報，21（4）:1-3.

王永剛，方國洪，曹德明，等，2004.渤、黃、東海潮汐的一種驗潮站資料同化數值模式.海洋科學進展，15（3）:253-274.

趙保仁，方國洪，曹德明，1994.渤、黃、東海潮汐潮流的數值模擬.海洋學報，13（5）:1-10.

（本文編輯：袁澤軼）

Refining barotropic tide simulations for the Bohai and Yellow Seas using a hybrid data assimilation approach

SONG Jun1,2,YAO Zhi-gang2,GUO Jun-ru3,LI Jing1,4,GAO Jia1,DONG Jun-Xing1
（1.National Marine DataandInformation Service,Tianjin 300171,China;2.The Ministryof Education Key Laboratoryof Physical Oceanography,Ocean University of China,Qingdao266100,China;3.National Marine Hazard Mitigation Service,SOA,Beijing 100194, China;4.College of Marine Science,Shanghai Ocean University,Shanghai 201306,China）

A high-resolution hybrid data assimilative(DA)modeling system is adopted to study on tidal characteristics, tidal mixing,residual currents,energy flux,dissipation and momentum term balance dynamics in the Bohai and Yellow Seas. Our solutions show together there are four amphidromic points in the Bohai and Yellow Seas.Significantly larger tidal ranges are found along the west coast of Korea Peninsula.Tidal current ellipse maps reveal an abrupt change in tidal rotation directions in the Yellow Sea around 35°N,and currents are the strongest along the west coast of Korea Peninsula,which are combined with rough and irregular bottom bathymetry,producing significant tidal mixing and energy dissipation.Momentum balance analyses on M2further confirm the complexity of tidal dynamics in this area.For most regions in our study area,the main dynamical balance is among the pressure gradient,local acceleration terms,and Coriolis force,and bottom friction and horizontal advection are smaller.But a simple linear,analytical model will not adequately reproduce tidal currents and elevations.

tidal characteristics;tidal dynamics;tidal residual circulation

P731.23

A

1001-6932（2016）05-0535-09

10.11840/j.issn.1001-6392.2016.05.008

2015-10-08；

2015-11-23

國家自然科學基金（41206013；41376014；41430963；41206004）；教育部物理海洋重點實驗室開放基金；2011年度高等學校博士學科點專項科研基金（20110132130001）；海洋公益性行業科研專項（201205018；201005019）；國家科技支撐計劃項目（2014BAB12B02）；天津市科技支撐計劃項目（14ZCZDSF00012）；國家海洋局青年科學基金重點項目（2012202；2013203；2012223）；國家建設高水平大學公派研究生項目（留金出 [2008]3019； [2012]3013）；國家海洋局空間遙感與應用研究重點實驗室開放基金重點課題（201601003）。

宋軍（1983-），男，博士，主要從事近海動力學、業務化海洋學方面的研究。

姚志剛，電子郵箱：874623647@qq.com。