高時(shí)友，何用，盧陳, 楊裕桂，吳門伍

(1.中山大學(xué) 地理科學(xué)與規(guī)劃學(xué)院，廣東 廣州 510275；2.水利部珠江河口動力學(xué)及伴生過程調(diào)控重點(diǎn)實(shí)驗(yàn)室，廣東 廣州 510611；3.珠江水利科學(xué)研究院，廣東 廣州 510611)

磨刀門口夏冬季沿岸流特征及成因分析

高時(shí)友1，2，3，何用，盧陳2，3, 楊裕桂2，3，吳門伍2，3

(1.中山大學(xué) 地理科學(xué)與規(guī)劃學(xué)院，廣東 廣州 510275；2.水利部珠江河口動力學(xué)及伴生過程調(diào)控重點(diǎn)實(shí)驗(yàn)室，廣東 廣州 510611；3.珠江水利科學(xué)研究院，廣東 廣州 510611)

根據(jù)2011—2012年磨刀門口的夏、冬季大、中、小潮定點(diǎn)觀測資料，對歐拉余流、斯托克斯余流以及拉格朗日余流進(jìn)行分析。結(jié)果表明：(1)外海測點(diǎn)的歐拉余流和拉格朗日余流，冬季在各潮型下均為一致的西南沿岸方向，夏季除東、西汊道點(diǎn)在強(qiáng)徑流下表現(xiàn)為順汊道指向外海方向，其余外海各點(diǎn)仍以西南沿岸方向?yàn)橹鳎?2)斯托克斯余流遠(yuǎn)小于歐拉余流，夏季明顯大于冬季，方向基本與歐拉余流相反；(3)冬季磨刀門口海域具有穩(wěn)定西南向沿岸流特征，南海東北季風(fēng)的驅(qū)動作用是其形成的主要原因；夏季磨刀門口各潮型下沿岸流特征各異，其影響的主要因素為徑流和風(fēng)，同時(shí)地形的影響不可忽略，特別是攔門沙形成的汊道分流作用，對強(qiáng)徑流作用下水沙輸移影響十分顯著。

沿岸流；磨刀門河口；歐拉余流；斯托克斯余流；拉格朗日余流

1 引言

珠江河口動力條件復(fù)雜，沿岸流特征顯著，形成機(jī)制復(fù)雜。珠江河口沿岸流規(guī)律及形成機(jī)制一直是眾多學(xué)者和專家關(guān)注的焦點(diǎn)，應(yīng)秩甫[1]認(rèn)為在珠江下泄的徑流和南海暖流的影響下，珠江口到湛江灣海域有一股常年存在的向西沿岸流，楊陽等[2]發(fā)現(xiàn)珠江沖淡水同時(shí)向粵西和粵東擴(kuò)展，逯玲燕等[3]認(rèn)為珠江口沿岸流冬季以粵西方向擴(kuò)展為主，而夏季則主要向粵東方向擴(kuò)展，Ou等[4]發(fā)現(xiàn)夏季珠江下泄徑流輸移方向與風(fēng)的關(guān)系極為密切，當(dāng)西南季風(fēng)盛行時(shí)珠江沖淡水向東擴(kuò)散更容易。楊士瑛等[5]發(fā)現(xiàn)夏季珠江口以西的近岸流大部分時(shí)間向西運(yùn)動。磨刀門口沿岸流既有珠江河口沿岸流的統(tǒng)一特性，同時(shí)還具有其獨(dú)特性。磨刀門作為珠江河口泄洪排沙主要通道，口門外沿岸流影響著咸淡水和泥沙的輸運(yùn)，認(rèn)識口外沿岸流特征對磨刀門口門延伸和治理意義重大。董兆英等[6]認(rèn)為磨刀門海區(qū)終年以下泄余流為主，近海區(qū)在徑流不斷擴(kuò)散和沿岸流系的影響下，余流速度在不同水域變化較大，流向以西南為主。賈良文等[7]根據(jù)實(shí)測資料計(jì)算所得的密度弗勞德數(shù)，認(rèn)為枯季磨刀門以浮力射流為主。方神光和王少波[8]根據(jù)實(shí)測水文資料，詳細(xì)分析了磨刀門漲落潮時(shí)表底層的水動力特性。可見學(xué)者對磨刀門沿岸流研究的側(cè)重點(diǎn)各不相同，系統(tǒng)地闡述磨刀門余流特征的研究并不多見。

珠江河口余流方向具有一定指向性，與長期徑流、泥沙、營養(yǎng)物等物質(zhì)輸運(yùn)關(guān)系密切，余流分布和特征一定程度反映沿岸流的變化特征。余流是河口物質(zhì)輸運(yùn)研究的重點(diǎn)，早在20世紀(jì)60年代，Bowden[9]就將余流分解為歐拉余流和斯托克斯余流，歐拉余流是對定點(diǎn)潮流的時(shí)間序列作潮周期平均，斯托克斯余流表達(dá)的則是潮流變化項(xiàng)與水深變化項(xiàng)的潮周期相關(guān)大小，表征了水流輸移量；將兩者矢量合成可得到拉格朗日余流，表征物質(zhì)凈輸移的大小和方向。本文以2011-2012年磨刀門口的定點(diǎn)觀測資料為基礎(chǔ)，分析磨刀門口夏冬季的歐拉余流、斯托克斯余流以及拉格朗日余流的分布特征，探討磨刀門口沿岸流輸移特征及成因。

2 數(shù)據(jù)資料

2.1 資料來源

珠江水利科學(xué)研究院受珠海市橫琴島經(jīng)濟(jì)區(qū)管理委員會、珠海藍(lán)灣建設(shè)開發(fā)有限公司委托，開展了磨刀門水道橫州口外河口夏季和冬季全潮水文測驗(yàn)海洋水文觀測。如表1所示，冬季的觀測時(shí)間為2011年12月4-13日，夏季觀測時(shí)間為2012年5月7-14日，均分別監(jiān)測了當(dāng)月的大、中和小潮時(shí)間段。觀測的內(nèi)容包括海流、水溫、鹽度、懸浮泥沙以及海底表層沉積物，各測量要素每小時(shí)整點(diǎn)觀測一次，連續(xù)觀測1個(gè)潮周期。

表1 夏、冬季觀測時(shí)間表

2.2 測點(diǎn)位置

磨刀門口海域位于橫琴島南面，東鄰澳門海域，北接磨刀門水道，河口外正對一“人”字形攔門沙中心區(qū)，攔門沙中心區(qū)兩側(cè)分別為東南和西南方向的東、西汊道及東、西區(qū)攔門沙。如圖1和2所示，定點(diǎn)監(jiān)測時(shí)間分為冬季和夏季，測點(diǎn)位置以B和C來表示，各布置7個(gè)測點(diǎn)，具體坐標(biāo)見表2。

B1位于磨刀門水道出口，B2、B3、B4自東向西分布于5 m等深線位置處，其中B3位于東汊出口處，B4位于攔門沙正南方，B5位于攔門沙西側(cè)，B6位于10 m等深線并正對攔門沙，B7位于攔門沙西南側(cè)的10 m等深線處。夏季測點(diǎn)與冬季略有偏差，其中C1同樣位于磨刀門水道出口處，比冬季靠內(nèi)側(cè)，C3和C4正對于東、西汊道出口，C6點(diǎn)正對B3點(diǎn)，其他測點(diǎn)與冬季相同。

表2 測點(diǎn)位置表

圖1 冬季測點(diǎn)位置Fig.1 Location of measuring points in winter

圖2 夏季測點(diǎn)位置Fig.2 Location of measuring points in summer

2.3 分析處理方法

潮汐余流是由動量方程中的平流項(xiàng)、連續(xù)方程中的非線性項(xiàng)以及非線性底摩擦效應(yīng)種種因素引起的，使其出現(xiàn)非線性效應(yīng)，一般用拉格朗日法和歐拉法進(jìn)行研究。歐拉余流適用于空間固定地點(diǎn)，而拉格朗日余流包含歐拉余流以及斯托克斯余流，它能完整地給出一個(gè)流體元的余流，為此依據(jù)實(shí)測水深和流速觀測資料，選取夏季、冬季大、中、小潮汐周期整點(diǎn)分層數(shù)據(jù)，按以下方法進(jìn)行分解計(jì)算。

(1)

(2)

x方向瞬時(shí)流速可以分解為：

(3)

水深可以表示為：

h(t)=h0+ht.

(4)

沿x軸方向的單寬潮周期平均輸水量為：

(5)

(6)

其中，

(7)

(8)

(9)

3 余流分析

3.1 歐拉余流

3.1.1 冬季

如圖3所示，冬季歐拉余流值整體表現(xiàn)為磨刀門出口位置小，外海各點(diǎn)位置大的特征。其中B1的大、中、小潮歐拉余流值分別為11cm/s、6cm/s、2cm/s。大、中潮的方向以落潮流方向?yàn)橹?，即順磨刀門口指向外海，約為145°；小潮由于河道下泄徑流量小，在外海較強(qiáng)的西南向輸運(yùn)動力下，其歐拉余流方向發(fā)生偏轉(zhuǎn)，變?yōu)?28°，即指向西南方向。外海B2～B7歐拉余流方向基本為240°，表現(xiàn)出明顯的沿著橫琴島岸線，向西南方向輸移的整體趨勢。

3.1.2 夏季

如圖4各潮型下的夏季歐拉余流值，表現(xiàn)為磨刀門出口大于外海各點(diǎn)的特征。在夏季強(qiáng)泄徑流的影響下，C1歐拉余流值均較大，大、中、小潮分別為31cm/s、22cm/s、34cm/s，方向均為順磨刀門口指向外海。夏季大潮外海各點(diǎn)除C3外，歐拉余流基本以沿岸向西南方向?yàn)橹鳌3由于位于東汊出口位置，在磨刀門水道徑流量足夠大的情況下，沖淡水沿東汊道下泄，因而形成了C3歐拉余流方向?yàn)檠貣|汊道方向，即133°左右。夏季中潮的外海各點(diǎn)的歐拉余流表現(xiàn)為一致的西南沿岸流特征，夏季小潮則表現(xiàn)為C2、C3、C4歐拉余流向西南方向，C6和C7向東方向。

圖3 冬季歐拉余流Fig.3 Euler residual current in winter

圖4 夏季歐拉余流Fig.4 Euler residual current in summer

3.2 斯托克斯余流

3.2.1 冬季

冬季各潮型斯托克斯余流均較小，最大值均出現(xiàn)在B1，大、中、小潮分別為4cm/s、1cm/s、1cm/s，方向基本與漲潮流方向一致，為順著河道指向上游。受河口的束窄作用，潮流的斯托克斯輸移作用明顯加強(qiáng)，如圖5。外海B2～B7斯托克斯余流均小于1cm/s，方向以東北方向?yàn)橹?，總體上各點(diǎn)的斯托克斯余流從大到小依次為大潮、中潮、小潮，斯托克斯余流大小值與潮汐呈正相關(guān)關(guān)系。

3.2.2 夏季

如圖6，夏季斯托克斯余流與冬季類似，同樣的是最大值出現(xiàn)在C1，方向?yàn)轫樦焙拥乐赶蛏嫌?，外海C2～C7斯托克斯余流均較小，除大、中潮的C3、C4以外，方向均以55°的東北方向?yàn)橹?。而大、中潮的C3由于受東汊漲潮流影響，其斯托克斯余流向?yàn)?35°的西北方向，C4則由于受西汊漲潮流影響，其斯托克斯余流向?yàn)?5°的東北方向。各潮型靠外海側(cè)的C6和C7的斯托克斯余流值，小于靠口門的C3、C4、C5。分析其原因是由于漲潮過程中靠口門的水深逐漸變淺，使得漲潮流變化量與潮汐變化量正相關(guān)性更強(qiáng)。

圖5 冬季斯托克斯余流Fig.5 Stokes residual current in winter

圖6 夏季斯托克斯余流Fig.6 Stokes residual current in summer

3.3 拉格朗日余流

3.3.1 冬季

如圖7，冬季拉格朗日余流均與歐拉余流相近，其中B1的歐拉余流為斯托克斯余流的3～5倍，方向正好相反，為此矢量合成的拉格朗日余流與歐拉余流方向一致。其余的外海點(diǎn)歐拉余流均為斯托克斯余流的幾十倍甚至上百倍，其大小和方向主要由歐拉余流決定。B1拉格朗日余流方向?yàn)轫樅拥老蛲夂５?45°為主，外海除B5以外，均為230°～250°之間，表現(xiàn)出一致的沿岸向西南方向流動的特征。

3.3.2 夏季

如圖8，夏季拉格朗日余流同樣表現(xiàn)為與歐拉余流相近的趨勢，河道內(nèi)的C1在強(qiáng)徑流作用下，沿口門向外海的拉格朗日余流特征顯著，外海各點(diǎn)的拉格朗日余流在不同潮型作用下則表現(xiàn)方向有所差異，總體表現(xiàn)為西南方向?yàn)橹鞯奶卣鳌?/p>

圖8 夏季拉格朗日余流Fig.8 Lagrange residual current in summer

4 沿岸流成因分析

由前面的余流分析可知，磨刀門口“2011.12”測次表現(xiàn)為顯著的西南沿岸流特征，“2012.5”測次的沿岸流方向則表現(xiàn)為不確定性，而表征物質(zhì)凈輸移大小和方向的拉格朗日余流與歐拉余流相近，說明沿岸流主要由形成歐拉余流的季風(fēng)、徑流、地形等因素決定。

(1)枯季盛行東北風(fēng)(9月中下旬至翌年4月)；洪季盛行西南風(fēng)(5月至9月上旬)；9月為西南季風(fēng)轉(zhuǎn)東北季風(fēng)時(shí)段；5月為東北季風(fēng)轉(zhuǎn)西南季風(fēng)時(shí)段。

南海是北太平洋西部最大的一個(gè)邊緣海，是連接太平洋和印度洋的一條重要通道。海域上空，枯季盛行東北風(fēng)，通常從9月中下旬到翌年的4月，洪季盛行西南風(fēng)，通常始于5月，終于9月上旬。

(2)枯季東北季風(fēng)是形成磨刀門口枯水期西南沿岸流的主要因素。

如圖9所示，測量期間B2、B6測站表明，測區(qū)風(fēng)向主要為N、NE、ENE向風(fēng)，枯季風(fēng)向穩(wěn)定，平均風(fēng)速接近6.0m/s，風(fēng)力較為強(qiáng)勁，磨刀門口海域在南海東北季風(fēng)的驅(qū)動下形成穩(wěn)定西南向沿岸流，各潮次的風(fēng)速從大到小依次為大潮、小潮、中潮，與枯季沿岸流強(qiáng)度有較好的一致相關(guān)性，為此，在枯季磨刀門口下泄徑流量偏小的條件下，枯季季風(fēng)是形成磨刀門口西南沿岸流的主要因素。

圖9 枯季逐時(shí)風(fēng)速、風(fēng)向Fig.9 Hourly wind speed and wind direction in the dry season

(3)“2012.5”測次的沿岸流方向不確定性，主要是由于5月轉(zhuǎn)季風(fēng)引起的，徑流和風(fēng)是影響沿岸流的主要動力，同時(shí)潮汐和地形作用也不可忽略。

洪季南海盛行西南季風(fēng)，而5月份則正逢東北季風(fēng)轉(zhuǎn)西南季風(fēng)階段[10]。如圖10所示的C2、C6測站風(fēng)況數(shù)據(jù)表明，風(fēng)向以東南向、南向?yàn)橹?，其中大潮以東南向?yàn)橹鳎谐眲t逐漸由南向轉(zhuǎn)為東北向，小潮以西南和東南為主。

圖10 洪季逐時(shí)風(fēng)速、風(fēng)向Fig.10 Hourly wind speed and wind direction in the flood season

由前面的數(shù)據(jù)統(tǒng)計(jì)結(jié)果可知，C1洪季大、中、小潮歐拉余流值分別為31cm/s、22cm/s、34cm/s，方向均為順直磨刀門水道指向外海。根據(jù)珠江委水文局發(fā)布的流量數(shù)據(jù)，西江上游馬口站5月日均下泄徑流量在5 000～15 000m3/s之間，觀測期大、中、小潮的馬口實(shí)際流量分別為6 900m3/s、5 520m3/s、9 190m3/s。結(jié)合對應(yīng)的馬口流量分析可知，馬口流量為從小到大依次為中潮、大潮、小潮，歐拉余流值也正好表現(xiàn)為小潮最大，中潮最小，即磨刀門河道出口的歐拉余流值的大小主要由徑流量決定。

洪季的外海各點(diǎn)拉格朗日余流在不同潮型下有所差異，大潮階段在徑流和東南向風(fēng)的共同影響下，東西兩泄洪汊道表現(xiàn)為順汊道的輸移方向，外海余流主要為西南向沿岸輸移特征。中潮階段由于東南向風(fēng)的存在，也使得余流為一致向西南向輸移的沿岸流特征，小潮階段由于逐漸轉(zhuǎn)向西南風(fēng)向，則出現(xiàn)了近岸點(diǎn)西南向飄移，遠(yuǎn)岸點(diǎn)東北向飄移的特征。即磨刀門口外海C2～C7測站洪季沿岸流以受徑流和風(fēng)的影響為主，主要偏向西南方向，同時(shí)潮汐和地形的影響也不可忽略，特別是攔門沙形成的汊道分流在強(qiáng)徑流作用下輸移特征十分明顯。

5 結(jié)論

本文基于2011-2012年磨刀門口夏冬季的大、中、小潮7個(gè)同步海流及風(fēng)況測點(diǎn)資料，統(tǒng)計(jì)分析了歐拉余流、斯托克斯余流以及拉格朗日余流的分布狀況，討論了磨刀門口沿岸流的成因，主要得到以下結(jié)論：

(1)歐拉余流和拉格朗日余流的大小和方向基本一致，磨刀門水道夏、冬季均以順河道指向口外為主，夏季明顯大于冬季。外海各點(diǎn)的歐拉余流和拉格朗日余流，冬季在各潮型下均為一致的西南沿岸方向，夏季的東、西汊道點(diǎn)在強(qiáng)徑流下表現(xiàn)為順汊道指向外海方向，其余外海各點(diǎn)仍以西南沿岸方向?yàn)橹鳌?/p>

(2)斯托克斯余流遠(yuǎn)小于歐拉余流，夏季明顯大于冬季，方向基本與歐拉余流相反。

(3)不同潮型的歐拉余流、斯托克斯余流和拉格朗日余流的大小和方向有所差異，其中斯托克斯余流與潮差有較好的正相關(guān)性。

(4)冬季磨刀門口海域具有穩(wěn)定西南向沿岸流特征，南海東北季風(fēng)的驅(qū)動作用是其形成的主要原因；夏季磨刀門口各潮型下沿岸流特征各異，其影響的主要因素為徑流和風(fēng)，同時(shí)地形的影響不可忽略，特別是攔門沙形成的汊道分流作用，對強(qiáng)徑流作用下水沙輸移影響十分顯著。

[1]應(yīng)秩甫. 粵西沿岸流及其沿岸沉積[J]. 中山大學(xué)學(xué)報(bào)(自然科學(xué)版), 1999, 38(3): 85-89.

YingZhifu.OnthecoastcurrentanditsdepositalongwesterncoastinGuangdong[J].ActaScientiarumNaturaliumUniversitatisSun-YatSeni, 1999, 38(3): 85-89.

[2] 楊陽, 孟強(qiáng), 夏華永, 等. 2006年夏季珠江沖淡水?dāng)U展及響應(yīng)[J]. 熱帶海洋學(xué)報(bào), 2010, 29(6): 15-21.

YangYang,MengQiang,XiaHuayong,etal.ExpansionofthePearlRiverdilutedwaterin2006anditsecologicalresponse[J].JournalofTropicalOceanography, 2010, 29(6): 15-21.

[3] 逯玲燕, 詹杰民, 耿兵緒. 基于通量收支分析的珠江口沖淡水?dāng)U展規(guī)律研究[J]. 水動力學(xué)研究與進(jìn)展, 2013, 28(3): 252-259.

LuLingyan,ZhanJieming,GengBingxu.StudyofthePearlRiverplumedispersionbasedonfluxbudgetanalysis[J].ChineseJournalofHydrodynamics, 2013, 28(3): 252-259.

[4]OuSuying,ZhangHong,WangDongxiao.HorizontalcharacteristicsofbuoyantplumeoffthePearlRiverEstuaryduringsummer[J].JournalofCoastalResearch, 2007, 50: 652-657.

[5] 楊士瑛, 鮑獻(xiàn)文, 陳長勝, 等. 夏季粵西沿岸流特征及其產(chǎn)生機(jī)制[J]. 海洋學(xué)報(bào), 2003, 25(6): 1-7.

YangShiying，BaoXianwen，ChenChangsheng,etal.AnalysisoncharacteristicsandmechanismofcurrentsysteminwestcoastofGuangdongProvinceinthesummer[J].HaiyangXuebao, 2003, 25(6): 1-7.

[6] 董兆英, 李素瓊, 盛晨. 珠江口余流特征初步分析[J]. 熱帶地理, 1985, 5(3): 177-185.

DongZhaoying,LiShuqiong,ShengChen.PreliminaryanalysisoftheresidualcurrentintheZhujiangestuary[J].TropicalGegoraphy, 1985, 5(3): 177-185.

[7] 賈良文, 吳超羽, 任杰, 等. 珠江口磨刀門枯季水文特征及河口動力過程[J]. 水科學(xué)進(jìn)展, 2006, 17(1): 82-88.

JiaLiangwen，WuChaoyu，RenJie,etal.HydrologiccharacteristicsandestuarinedynamicprocessduringthedryseasoninModaomenEstuaryofthePearlRiver[J].AdvancesinWaterScience, 2006, 17(1): 82-88.

[8] 方神光, 王少波. 磨刀門水道枯季水動力特性分析[J]. 水文, 2013, 33(5): 70-74.

FangShenguang,WangShaobo.AnalysisofwaterdynamiccharacteristicsofModaomenwaterwayduringdryseason[J].JournalofChinaHydrology, 2013, 33(5): 70-74.

[9]BowdenKF.Themixingprocessesintidalestuary[J].InternationalJournalofAirandWaterPollution, 1963, 7: 343-356.

[10]楊士瑛, 鮑獻(xiàn)文, 陳長勝, 等. 夏季粵西沿岸流特征及其產(chǎn)生機(jī)制[J]. 海洋學(xué)報(bào), 2003, 25(6):1-8.

YangShiying,BaoXianwen,ChenChangsheng,etal.AnalysisoncharacteristicsandmechanismofcurrentsysteminwestcoastofGuangdongProvinceinthesummer[J].HaiyangXuebao, 2003, 25(6):1-8.

The characteristics of alongshore current in flood and the dry season in Modaomen Estuary of the Pearl River and the reason analysis

Gao Shiyou1,2,3, He Yong2,3, Lu Chen2,3,Yang Yugui2,3,Wu Menwu2,3

(1.SchoolofGeographyandPlanningofSunYat-SenUniversity,Guangzhou510275,China;2.KeyLaboratoryofthePearlRiverEstuarineDynamicsandAssociatedProcessRegulation,MinistryoftheWaterResources,Guangzhou510611,China;3.PearlRiverHydraulicResearchInstitute,Guangzhou510611,China)

According to the observation data in spring tide, middle tide, neap tide in flood and dry season at Modaomen Estuary from 2011 to 2012, this paper analyzes Euler residual current, Stokes residual current and Lagrange residual flow. The results show that: (1) the main direction of Euler residual current and Lagrange residual flow are southwest in winter and summer but the current in east and west branch show offshore under the strong runoff; (2) the Stokes residual current is much less than Euler residual current and much greater in summer than in winter, the direction of the Stokes residual current is in the opposite of Euler residual current; (3) there is a stable southwest alongshore current in winter in Modaomen Estuary, the driving effect of northeast monsoon in the South China Sea is the main reason, there is no consistent characteristics of alongshore current in summer, the main influencing factors is runoff and the wind, and the influence of the terrain especially the shunt effect caused by the mouth bar should not be neglected, it will influence the transport of water and sediment which under the strong runoff significantly.

alongshore current; Modaomen Estuary; Euler residual current; Stokes residual current; Lagrange residual flow

10.3969/j.issn.0253-4193.2017.05.001

2016-04-20；

2016-11-09。

水利部公益性行業(yè)科研專項(xiàng)(201501010)；國家自然科學(xué)基金資助項(xiàng)目(51409286)；珠科院開放基金([2013]KJ05)。

高時(shí)友(1978—)，男，安徽省廬江縣人，主要從事河道與港航工程。E-mail：pearlgsy@qq.com

*通信作者：何用(1977—)，男，湖北省黃岡市人，工學(xué)博士，主要從事河口海岸動力學(xué)研究。E-mail：8459585@qq.com

P731.2

A

0253-4193(2017)05-0001-09

高時(shí)友，何用，盧陳，等. 磨刀門口夏冬季沿岸流特征及成因分析[J]. 海洋學(xué)報(bào), 2017, 39(5): 1-9,

Gao Shiyou, He Yong, Lu Chen, et al. The characteristics of alongshore current in flood and the dry season in Modaomen Estuary of the Pearl River and the reason analysis[J]. Haiyang Xuebao, 2017, 39(5): 1-9, doi：10.3969/j.issn.0253-4193.2017.05.001