潘明婕，楊 芳，荊 立，羅照陽(yáng)，王 青，孔 俊

(1.海岸災(zāi)害及防護(hù)教育部重點(diǎn)實(shí)驗(yàn)室(河海大學(xué)), 江蘇 南京 210098; 2.南京昊控軟件技術(shù)有限公司, 江蘇 南京 211100; 3.珠江水利委員會(huì)珠江水利科學(xué)研究院，廣東 廣州 510611; 4.生態(tài)環(huán)境部南京環(huán)境科學(xué)研究所, 江蘇 南京 210042)

磨刀門水道是中山、珠海、澳門特別行政區(qū)的主要飲用水水源地，其水質(zhì)安全對(duì)周邊城市的經(jīng)濟(jì)、社會(huì)發(fā)展有重要的影響。2000年以來(lái)，受河道疏浚、圍墾工程、人工采砂等人類活動(dòng)的影響，磨刀門鹽水入侵問(wèn)題日益嚴(yán)峻，極大地威脅到周圍城市的供水安全[1]。磨刀門水道地處亞熱帶季風(fēng)氣候區(qū)，臨近南海，容易受到熱帶氣旋的襲擊。根據(jù)1949—2008年在珠江三角洲登陸的熱帶氣旋資料統(tǒng)計(jì)，平均每年有5.7個(gè)熱帶氣旋在此登陸。受臺(tái)風(fēng)天氣的影響，磨刀門水道鹽水入侵的形勢(shì)變得更為復(fù)雜。

針對(duì)河口地區(qū)的咸潮上溯問(wèn)題，近年來(lái)國(guó)內(nèi)外眾多學(xué)者開展了深入的研究[2-3]。Li等[4]揭示了強(qiáng)臺(tái)風(fēng)影響下，半封閉海灣經(jīng)歷了強(qiáng)鹽水上涌、鹽淡水分層破壞以及臺(tái)風(fēng)過(guò)后受重力調(diào)整的再次層化過(guò)程。Gong等[5]基于大量實(shí)測(cè)數(shù)據(jù)，利用EFDC建立了磨刀門水道三維水流鹽度數(shù)值模型，分析了枯季鹽度輸運(yùn)動(dòng)力特征，結(jié)果顯示磨刀門鹽水入侵主要受潮汐和徑流的相互作用，大的徑流量能有效抑制鹽水入侵。劉吉等[6]通過(guò)分析2008年強(qiáng)臺(tái)風(fēng)“黑格比”登陸期間磨刀門水道縱向鹽度變化發(fā)現(xiàn)，風(fēng)暴潮期間磨刀門水體鹽度具有突發(fā)、單峰特征，且鹽度變化過(guò)程與增水過(guò)程基本一致。Burchard等[7]指出，在較弱鹽度梯度力作用下，潮汐應(yīng)變控制交換余環(huán)流(residual exchange circulation)形成縱向河口環(huán)流，且河口環(huán)流強(qiáng)度隨鹽度梯度的增大而明顯增大。目前，關(guān)于臺(tái)風(fēng)對(duì)咸潮的影響，大多集中于定性和現(xiàn)象的描述，對(duì)其動(dòng)力作用機(jī)制缺乏深入的探討。

本文選取具有代表性的臺(tái)風(fēng)“納沙”，基于SCHISM模型，建立了臺(tái)風(fēng)天氣下磨刀門水道三維潮流鹽度數(shù)值模型，并結(jié)合辛普森數(shù)(Simpson number)分析了臺(tái)風(fēng)過(guò)境期間水道河口環(huán)流形成的動(dòng)力機(jī)制，同時(shí)采用鹽通量機(jī)制分析法解釋臺(tái)風(fēng)期鹽度分布“雙峰”特征形成的動(dòng)力機(jī)制。

1 臺(tái)風(fēng)“納沙”過(guò)境期間特殊動(dòng)力條件與現(xiàn)象

圖1為研究區(qū)域和水位、鹽度測(cè)站位置示意圖(圖中S—S為水道沿程縱斷面，斷面TR為橫斷面，1、2、3、4為所選分析點(diǎn))。“納沙”為2011年登陸珠江的年度最強(qiáng)臺(tái)風(fēng)。結(jié)合表1可知，移動(dòng)過(guò)程中“納沙”先后于9月26日23時(shí)和29日7時(shí)兩次加強(qiáng)到強(qiáng)臺(tái)風(fēng)強(qiáng)度且長(zhǎng)期維持在臺(tái)風(fēng)強(qiáng)度。此次臺(tái)風(fēng)影響范圍廣、作用強(qiáng)度大，對(duì)我國(guó)海南、廣東和廣西以及菲律賓多地造成重大影響。

圖1 研究區(qū)域和水位、鹽度測(cè)站位置

“納沙”過(guò)境期間，磨刀門水道出現(xiàn)嚴(yán)重的鹽水入侵，平崗泵測(cè)站鹽度過(guò)程呈現(xiàn)出特殊的“雙峰”現(xiàn)象(圖2)。該現(xiàn)象的形成還受到上游徑流量影響，由于2011年珠江流域遭受干旱災(zāi)害，9月流量長(zhǎng)期維持在2 000 m3/s左右，低徑流量為咸潮上溯創(chuàng)造了有利條件。

圖2 “納沙”過(guò)境期間磨刀門水道表層鹽度和流量過(guò)程

2 研究方法

2.1 數(shù)值模型

2.1.1模型建立

磨刀門水道三維水流鹽度數(shù)值模型采用Zhang等[8]開發(fā)的、廣泛用于河口海洋水動(dòng)力問(wèn)題研究的SCHISM跨尺度湖泊-河流-河口-海洋水動(dòng)力模型。模型水平方向采用無(wú)結(jié)構(gòu)網(wǎng)格(圖3)，垂向采用SZ混合坐標(biāo)，輸運(yùn)方程采用TVD2(高階隱式平流格式)求解。上游流量邊界和外海水位、鹽度邊界從珠江大范圍潮流數(shù)學(xué)模型提取，二維模型范圍及網(wǎng)格如圖4所示。大范圍二維模型上游流量邊界為實(shí)測(cè)流量，外海水位由潮汐預(yù)報(bào)值和實(shí)測(cè)水位值綜合率定調(diào)整得到，鹽度邊界為定值33 psu；風(fēng)速條件采用構(gòu)建的臺(tái)風(fēng)場(chǎng)計(jì)算得出，臺(tái)風(fēng)場(chǎng)由CCMP/NECP背景風(fēng)場(chǎng)和臺(tái)風(fēng)經(jīng)驗(yàn)?zāi)Ｐ惋L(fēng)場(chǎng)合成得到。詳細(xì)構(gòu)建方法參見文獻(xiàn)[9]。模型運(yùn)行102 d達(dá)到穩(wěn)定狀態(tài)。

圖3 磨刀門水道三維模型網(wǎng)格

表1 “納沙”臺(tái)風(fēng)中心位置、氣壓、風(fēng)速特征

圖4 珠江二維模型網(wǎng)格

2.1.2模型驗(yàn)證

采用“納沙”過(guò)境期間實(shí)測(cè)水位和表層鹽度數(shù)據(jù)對(duì)磨刀門水道三維模型進(jìn)行驗(yàn)證，水位和鹽度測(cè)站位置見圖1。模型驗(yàn)證結(jié)果如圖5和圖6所示，由于臺(tái)風(fēng)期外海風(fēng)浪較大，大橫琴測(cè)站部分鹽度數(shù)據(jù)缺測(cè)。驗(yàn)證結(jié)果顯示模型模擬結(jié)果較好，可用于后期計(jì)算分析。

2.2 辛普森數(shù)分析法

潮汐河口受水平水體密度梯度影響，潮周期過(guò)程水道中主要呈現(xiàn)出3種基本鹽淡水狀態(tài)：完全混合狀態(tài)、應(yīng)變致周期性層化(strain-induced periodic stratification，SPIS)狀態(tài)以及持續(xù)層化狀態(tài)[10]。采用辛普森數(shù)Si可以刻畫潮汐應(yīng)變的動(dòng)態(tài)變化，其定義為應(yīng)變引起的勢(shì)能變化與湍動(dòng)能產(chǎn)生率之比[7,11-12]：

(1)

其中

(a)三灶測(cè)站

(b)大橫琴測(cè)站

(c)竹銀測(cè)站

(a)大橫琴測(cè)站

(b)掛定角測(cè)站

(c)平崗泵測(cè)站

η為潮位幅值。當(dāng)辛普森數(shù)較小時(shí)，落潮時(shí)湍動(dòng)能作用超過(guò)潮汐應(yīng)變的層化穩(wěn)定作用，導(dǎo)致水體垂向完全混合；當(dāng)辛普森數(shù)處于中間大小時(shí)，水體呈現(xiàn)SPIS狀態(tài)，即落潮期水體層化、漲潮期水體混合；當(dāng)辛普森數(shù)較大時(shí)，河道長(zhǎng)期處于層化狀態(tài)[11]。針對(duì)不同河口，辛普森數(shù)的臨界值有所不同，Simpson等[12]針對(duì)Liverpool Bay劃定：完全混合和SPIS狀態(tài)的臨界值為8.8×10-2，SPIS和持續(xù)層化狀態(tài)的臨界值為8.4×10-1；而Stacey等[13]取0.2作為落潮期持續(xù)層化和混合的臨界值，即此時(shí)垂向混合和SPIS作用相平衡。此外，辛普森數(shù)的臨界值大小還受風(fēng)應(yīng)力[14]、地球自轉(zhuǎn)和相對(duì)潮汐頻率[15]的影響。

2.3 鹽通量機(jī)制分析法

為分析臺(tái)風(fēng)期鹽度輸運(yùn)的動(dòng)力機(jī)制，采用Lerczak等[16]提出的鹽通量機(jī)制分析法，其中總鹽通量Fs表示為

Fs=〈?usdA〉

(2)

式中：u為橫斷面法向流速；s為鹽度；A為橫斷面面積；角括號(hào)表示33 h的低通濾波，角括號(hào)內(nèi)橫斷面積分則得到瞬時(shí)的鹽通量值。Fs可被分解為

Fs=〈?(u0+uE+uT)(S0+SE+ST)dA〉≈

〈?(u0S0+uESE+uTST)dA〉=

QfS0+FE+FT

(3)

式中：u0、uE、uT分別為潮平均斷面平均的流速、潮平均斷面變化的流速和潮汐變化斷面變化的流速、S0、SE、ST分別為潮平均斷面平均的鹽度、潮平均斷面變化的鹽度和潮汐變化斷面變化的鹽度。機(jī)制分解結(jié)果顯示，鹽通量輸運(yùn)取決于潮平均斷面平均的平流輸運(yùn)QfS0、潮平均剪切擴(kuò)散FE和潮汐震蕩FT三者動(dòng)力輸運(yùn)的平衡。此外，Lerczak等[16]特別指出，F(xiàn)E和FT項(xiàng)主要驅(qū)動(dòng)鹽分向陸輸運(yùn)，而QfS0項(xiàng)主要受上游徑流量影響，驅(qū)動(dòng)鹽分向海輸運(yùn)。

3 結(jié)果分析與討論

針對(duì)縱斷面沿程的流速、鹽度分布，分析臺(tái)風(fēng)過(guò)境前、中、后期的變化差異，研究臺(tái)風(fēng)對(duì)磨刀門水道流速、鹽度分布的影響；并采用辛普森數(shù)來(lái)定量分析磨刀門水道潮周期鹽淡水混合狀態(tài)，探究河口縱向環(huán)流的形成機(jī)制。此外，針對(duì)“納沙”過(guò)境期間鹽度分布的“雙峰”性，設(shè)置數(shù)值試驗(yàn)并結(jié)合鹽通量機(jī)制分析法，探究其形成的動(dòng)力機(jī)制。

3.1 縱向水流鹽度時(shí)空分布特征

通過(guò)分析臺(tái)風(fēng)過(guò)境前、中、后期縱向流速和鹽度分布特點(diǎn)，詳細(xì)對(duì)比不同時(shí)期縱斷面差異性的動(dòng)力特征，突出臺(tái)風(fēng)對(duì)縱向水流、鹽度結(jié)構(gòu)的影響。圖7為三灶測(cè)站“納沙”過(guò)境期間水位過(guò)程線(圖中紅色為潮周期歷時(shí))，統(tǒng)一選取落憩漲初時(shí)刻(圖7中虛線)作為研究時(shí)刻，得到圖1中S—S斷面(從口門外攔門沙至竹銀上游)流速、鹽度分布如圖8所示。

a. 臺(tái)風(fēng)來(lái)臨前，磨刀門水道正處于小潮轉(zhuǎn)大潮時(shí)期(圖8(a))。21日小潮期垂向等鹽度線分布最為密集，表底層鹽度差較大，鹽淡水分層最為明顯。由于落憩時(shí)刻受鹽度梯度影響，水道中普遍存在著表層向海、底層向陸的縱向環(huán)流，鹽度從底層不斷上溯。洪灣水道處鹽度沿程分布不連續(xù)，鹽度明顯高于兩側(cè)，可見洪灣水道向主水道輸送較多鹽分，有助于鹽分向上游擴(kuò)散。此外，由于受倒坡地形影響，鹽水被凸起的地形阻擋，小潮動(dòng)力較弱，難以將鹽水排出，從而囤積在上游水道中。

圖7 “納沙”過(guò)境期間三灶測(cè)站水位過(guò)程線

(a)“納沙”過(guò)境前小潮期(2011-09-21T13：00)

(b)“納沙”過(guò)境前中潮期(2011-09-24T16：00)

(c)“納沙”過(guò)境前大潮期(2011-09-27T17：00)

(d)“納沙”過(guò)境期間(2011-09-29T19：00)

(e)“納沙”過(guò)境后(2011-10-01T9：00)

b. 24日中潮期(圖8(b))，表底層鹽度差減小，鹽淡水分層相對(duì)小潮期減弱，水道中縱向環(huán)流強(qiáng)度也隨著鹽度梯度減小而有所減弱。但由于潮動(dòng)力的增強(qiáng)，小潮期在水道中囤積的鹽分得以向上游不斷擴(kuò)散，導(dǎo)致鹽水上溯距離進(jìn)一步增大。

c. 隨著27日大潮期(圖8(c))的到來(lái)，表底層鹽度差銳減，鹽淡水分層減弱，混合加強(qiáng)，鹽水近乎垂直地退出水道，縱向環(huán)流基本消失。

d. 臺(tái)風(fēng)過(guò)境期間，磨刀門水道處于大潮轉(zhuǎn)中潮時(shí)期(圖8(d))。咸潮上溯本應(yīng)較大潮期更弱，但受臺(tái)風(fēng)引發(fā)的外海增水影響，口門外高濃度鹽水短時(shí)間大量涌入水道，導(dǎo)致水道中、下游鹽度整體偏高且分層明顯，進(jìn)而在水道中部形成了較大范圍的縱向環(huán)流。

e. 臺(tái)風(fēng)過(guò)境后，磨刀門水道大致處于中潮期(圖8(e))，臺(tái)風(fēng)期涌入水道的鹽水不斷向外海排出。10月1日9時(shí)，臺(tái)風(fēng)影響基本消失，水道恢復(fù)正常天氣下大潮后中潮期的流速、鹽度分布狀態(tài)，即水道混合較強(qiáng)，落潮時(shí)等鹽度線垂直地退出水道。

3.2 河口環(huán)流形成的動(dòng)力特征

采用辛普森數(shù)分析法，針對(duì)臺(tái)風(fēng)過(guò)境期間的特征潮周期歷時(shí)(圖7中紅線)，分析水道沿程1、2、3、4點(diǎn)(圖1)位置處的縱向環(huán)流形成的動(dòng)力機(jī)制。4個(gè)沿程點(diǎn)分別位于磨刀門水道口門、下游水道中部、洪灣水道與主水道交匯處以及水道中游彎曲處。由于不同河口地形、動(dòng)力條件的差異，辛普森數(shù)的臨界值有所不同(表2)。本文結(jié)合磨刀門水道鹽淡水混合狀態(tài)的變化情況，參考Simpson等[12]以及Stacey等[13]研究成果，初步擬定辛普森數(shù)在0.08 ～ 0.12之間為SPIS狀態(tài)，高于0.12為持續(xù)層化狀態(tài)，低于0.08為完全混合狀態(tài)。

表2 不同時(shí)期4個(gè)沿程點(diǎn)的辛普森數(shù)

a. 小潮期2、4點(diǎn)的辛普森數(shù)較大(>0.1)，其中2點(diǎn)的值超過(guò)0.15，可見該點(diǎn)附近產(chǎn)生持續(xù)層化，下游水道中部的河口環(huán)流不僅是由潮汐應(yīng)變引起的，更主要的是強(qiáng)鹽水體密度梯度力驅(qū)動(dòng)下的重力環(huán)流；4點(diǎn)的值略小于0.15，水道中游鹽度值相對(duì)下游整體偏低，水體密度梯度較小，局部縱向環(huán)流的主要驅(qū)動(dòng)力還是潮汐應(yīng)變，但是由于所處位置河道較為彎曲易于產(chǎn)生側(cè)向環(huán)流，進(jìn)而可能對(duì)縱向環(huán)流產(chǎn)生影響；1點(diǎn)的辛普森數(shù)最小，表示水體混合強(qiáng)烈，這主要與其所處位置密切相關(guān)，即口門處受潮汐混合作用最為強(qiáng)烈；3點(diǎn)的辛普森數(shù)為0.083 4，受洪灣水道潮周期性輸水輸鹽影響，加強(qiáng)了水體混合，打破了主水道的高度層化狀態(tài)，具有SPIS效應(yīng)。

b. 中潮期仍是2點(diǎn)和4點(diǎn)的辛普森數(shù)較大(在0.09左右)，3點(diǎn)次之，1點(diǎn)最小。除口門1點(diǎn)變化不大，其余點(diǎn)位均較小潮期有所減小，表明潮汐混合動(dòng)力增強(qiáng)，水道整體呈現(xiàn)SPIS狀態(tài)。由于水道中水體密度梯度較小，重力環(huán)流較弱，河口環(huán)流主要由潮汐應(yīng)變驅(qū)動(dòng)且?guī)缀踔辉诼涑逼诔霈F(xiàn)。

c. 大潮期潮汐混合作用進(jìn)一步加大，各點(diǎn)位的辛普森數(shù)進(jìn)一步減小，水道以完全混合為主，僅落潮期出現(xiàn)弱層化狀態(tài)，河口環(huán)流基本消失。小、中、大潮期辛普森數(shù)的變化與臺(tái)風(fēng)過(guò)境前縱斷面流速、鹽度分布變化(圖8(a)(b)(c))相一致，解釋了各時(shí)期環(huán)流形成的主要原因。

d. 臺(tái)風(fēng)期辛普森數(shù)的分布發(fā)生較大變化，1、2、3點(diǎn)較大潮期進(jìn)一步減小，而4點(diǎn)卻大幅增大，表明下游水道受臺(tái)風(fēng)影響混合更為強(qiáng)烈，而中游處隨著鹽水囤積增多，形成一定的水體密度梯度，導(dǎo)致水體層化。水道中縱向環(huán)流的形成，主要受潮汐應(yīng)變作用，在落潮期出現(xiàn)，與圖8(d)中環(huán)流出現(xiàn)時(shí)間及位置一致。

e. 大潮后中潮期辛普森數(shù)除1點(diǎn)外，均較大潮期有所減小，經(jīng)過(guò)大潮期以及臺(tái)風(fēng)期強(qiáng)烈的混合作用，水道處于完全混合狀態(tài)，與小潮后中潮的SIPS狀態(tài)差別較大。

總體看來(lái)，辛普森數(shù)沿程分布存在差異，口門處1點(diǎn)辛普森數(shù)持續(xù)較小，水體長(zhǎng)期處于完全混合狀態(tài)，2、4點(diǎn)位置處在正常潮周期內(nèi)，即從小潮至大潮后中潮期辛普森數(shù)不斷變小，水道大體經(jīng)歷了“層化—SPIS—完全混合”的過(guò)程，3點(diǎn)受洪灣水道影響，比2或4點(diǎn)混合強(qiáng)烈，但是弱于1點(diǎn)。在洪枯季之交，上游徑流量(馬口站)為2 000 m3/s時(shí)，當(dāng)水道處于高度層化狀態(tài)時(shí)(小潮期)，縱向環(huán)流主要是受水體密度梯度驅(qū)動(dòng)下的重力環(huán)流；水道處于SPIS狀態(tài)時(shí)(小潮后中潮期)，縱向環(huán)流的形成主要受潮汐應(yīng)變影響；水道處于完全混合狀態(tài)時(shí)(大潮期及大潮后中潮期)，縱向環(huán)流幾乎不可見。此外，辛普森數(shù)對(duì)臺(tái)風(fēng)的響應(yīng)明顯，改變了大潮后期水體完全混合的狀態(tài)，水道中游出現(xiàn)特殊的SPIS狀態(tài)，受潮汐應(yīng)變影響落潮時(shí)形成縱向環(huán)流。

3.3 鹽度分布“雙峰”特征形成的動(dòng)力機(jī)制

結(jié)合前文實(shí)測(cè)數(shù)據(jù)(圖2)分析可知，臺(tái)風(fēng)“納沙”過(guò)境期間，磨刀門水道鹽度測(cè)站(平崗泵和掛定角)鹽度變化過(guò)程呈現(xiàn)出“雙峰”特征。針對(duì)這一特殊現(xiàn)象，設(shè)計(jì)去除“納沙”影響的數(shù)值試驗(yàn)，并與經(jīng)過(guò)實(shí)測(cè)數(shù)據(jù)驗(yàn)證的“納沙”臺(tái)風(fēng)模型模擬結(jié)果進(jìn)行對(duì)比，探究“雙峰”特征形成的動(dòng)力機(jī)制，結(jié)果見圖9。

圖9 掛定角測(cè)站鹽度過(guò)程線模擬

從掛定角測(cè)站處的表層鹽度變化過(guò)程可以直觀地看出，鹽度分布的“前峰”變化較為緩慢，而“后峰”的形成及衰落都十分迅速，且峰值高于“前峰”。去除臺(tái)風(fēng)影響后，鹽度分布的“后峰”隨之消失，表明“后峰”的形成與此次臺(tái)風(fēng)密切相關(guān)，而“前峰”的鹽度幾乎不變，應(yīng)與此次臺(tái)風(fēng)無(wú)關(guān)。從潮周期來(lái)看，9月20日為小潮期，27日為大潮期，鹽度“前峰”出現(xiàn)在25日，正是小潮后的中潮期。胡溪等[17]研究發(fā)現(xiàn)，磨刀門水道咸潮上溯過(guò)程具有明顯的半月周期變化，鹽水入侵在小潮期加強(qiáng)，在大潮期較弱，咸界上溯最遠(yuǎn)出現(xiàn)在小潮之后的中潮期，可見鹽度分布的“前峰”很可能是周期性的鹽水入侵所致。

通過(guò)計(jì)算“納沙”臺(tái)風(fēng)模型和去臺(tái)風(fēng)影響模型三灶測(cè)站位置處的水位差，得到臺(tái)風(fēng)“納沙”過(guò)境期間磨刀門水道大致的增水情況，如圖10(圖中時(shí)間以9月28日0時(shí)為起點(diǎn)“0”)所示?？梢?，此次臺(tái)風(fēng)增水約持續(xù)了36 h，最大達(dá)到1.0 m左右，使得磨刀門水道的外海動(dòng)力加強(qiáng)，推動(dòng)鹽水向河道內(nèi)涌入。

圖10 三灶測(cè)站臺(tái)風(fēng)“納沙”過(guò)境期間增水過(guò)程線

進(jìn)一步計(jì)算水道斷面TR(圖1)的鹽通量(圖11)，從鹽水入侵動(dòng)力機(jī)制上對(duì)比分析鹽度“雙峰”特征形成的原因。由圖11可見，9月20日小潮期總鹽通量向陸輸送，鹽分不斷向水道上游累積，其中平流通量和潮汐震蕩通量變化不大，而剪切通量向陸輸送明顯增大，成為總鹽通量向陸輸送增大的主導(dǎo)因素。剪切通量與橫斷面流速分布有關(guān)，剪切通量增大表明表底層流速差大，鹽淡水分層明顯，重力環(huán)流不斷加強(qiáng)，高濃度鹽水從底層上溯，鹽度分布不斷上漲。23日之后水道進(jìn)入中潮期，總鹽通量轉(zhuǎn)為向海輸出，但量級(jí)較小，剪切通量有所減小，鹽水從底層上溯變緩，鹽度分布達(dá)到峰值。隨著27日大潮的到來(lái)，剪切通量降到零點(diǎn)，總鹽通量主要受平流通量控制，以向海輸出為主，鹽度分布逐漸降低。因此，在小潮轉(zhuǎn)大潮的過(guò)程中，形成了鹽度分布的“前峰”。此外，磨刀門鹽水入侵長(zhǎng)度受徑流量影響很大，此次“前峰”產(chǎn)生的重要前提條件是，水道上游徑流量長(zhǎng)時(shí)間維持在2 000 m3/s左右。

(a)“納沙”臺(tái)風(fēng)模型

(b)去臺(tái)風(fēng)影響模型

而在臺(tái)風(fēng)影響下，28—30日總鹽通量向陸輸送大幅增加，相比于去臺(tái)風(fēng)影響后，總鹽通量較小且向海輸出，可見此次臺(tái)風(fēng)對(duì)鹽度輸運(yùn)產(chǎn)生了重要影響。臺(tái)風(fēng)“納沙”過(guò)境期間，強(qiáng)烈的氣旋風(fēng)應(yīng)力導(dǎo)致外海增水，使得平流通量轉(zhuǎn)而向陸輸送，大量高濃度鹽水從口門向上游輸送，導(dǎo)致水道鹽度的突然升高，由此產(chǎn)生鹽度分布的“后峰”。

4 結(jié) 論

a. 臺(tái)風(fēng)促進(jìn)外海水體混合，不同于小潮期鹽水逐漸在水道內(nèi)囤積，而是高濃度鹽水在漲潮期近乎垂直地向水道內(nèi)推進(jìn)，導(dǎo)致鹽水上溯距離在短時(shí)間內(nèi)大幅增加。此外，高濃度鹽水受地形等因素影響難以及時(shí)排出，在水道內(nèi)形成明顯的鹽淡水分層，進(jìn)而引發(fā)較強(qiáng)的縱向環(huán)流。

b. 結(jié)合辛普森數(shù)定量分析發(fā)現(xiàn)，與臺(tái)風(fēng)前、后期相比，臺(tái)風(fēng)期辛普森數(shù)變化較大，水道下游仍處于強(qiáng)烈混合狀態(tài)，而水道中游隨著鹽水囤積，出現(xiàn)特殊的SPIS狀態(tài)。水道中縱向環(huán)流的形成，主要是受潮汐應(yīng)變的驅(qū)動(dòng)。

c. 臺(tái)風(fēng)“納沙”過(guò)境期間咸潮上溯導(dǎo)致鹽度分布呈現(xiàn)“雙峰”特征。鹽通量機(jī)制分析結(jié)果顯示，“前峰”是持續(xù)低流量條件下，受剪切通量控制，小潮轉(zhuǎn)大潮時(shí)期鹽水隨潮波周期性上溯；“后峰”則是臺(tái)風(fēng)“納沙”的增水效應(yīng)導(dǎo)致平流通量轉(zhuǎn)而向陸輸送，導(dǎo)致鹽度突然增大。