珠江口雞啼門近海海洋對臺風(fēng)的響應(yīng)*

陳杏文 邱春華, 3 張 恒 劉 東 石 慧 馬永貴 王東曉

珠江口雞啼門近海海洋對臺風(fēng)的響應(yīng)*

陳杏文1, 2邱春華1, 2, 3張 恒1, 2劉 東4石 慧1, 5馬永貴1, 2王東曉1, 2①

(1. 中山大學(xué)海洋科學(xué)學(xué)院 廣東珠海 519082; 2. 珠江口海洋生態(tài)環(huán)境教育部野外科學(xué)觀測研究站 廣東珠海 519082; 3. 自然資源部海洋環(huán)境探測技術(shù)與應(yīng)用重點(diǎn)實(shí)驗(yàn)室 廣東廣州 510310; 4. 廣州睿海海洋科技有限公司 廣東廣州 511400; 5. 中水珠江規(guī)劃勘測設(shè)計有限公司 廣東廣州 510610)

海洋對臺風(fēng)的響應(yīng)對臺風(fēng)預(yù)報具有重要意義, 但是由于動力結(jié)構(gòu)比較復(fù)雜, 近岸水體對臺風(fēng)的響應(yīng)仍然不明確?；谥榻诤Ｑ笊鷳B(tài)環(huán)境教育部野外科學(xué)觀測研究站的1號浮標(biāo)觀測資料(2021年7～8月), 分析了珠江口雞啼門海域在臺風(fēng)“查帕卡”和“盧碧”期間海洋動力學(xué)和熱力學(xué)響應(yīng)特征。臺風(fēng)過境引起水體的緯向流速顯著增強(qiáng), “查帕卡”激發(fā)了強(qiáng)烈的順時針近慣性振蕩, “盧碧”則激發(fā)了接近慣性頻率的逆時針海水運(yùn)動。受當(dāng)?shù)販\水深的影響, 近慣性能量在一個慣性周期(32 h)后快速耗散。兩個臺風(fēng)引起的海表溫度(sea surface temperature, SST)降溫幅度均大于2.2 °C, 這種降溫受潛熱通量影響較大; 同時, 臺風(fēng)過后晚上的SST均大于白天的SST, 這也是受到潛熱通量的影響。

臺風(fēng); 近慣性振蕩; 海表溫度; 熱通量; 雞啼門; 珠江口

臺風(fēng)是世界上最具破壞力的自然災(zāi)害之一。登陸時, 臺風(fēng)帶來的強(qiáng)風(fēng)暴雨、巨浪以及風(fēng)暴潮會對沿海地區(qū)造成破壞。1983～2006年, 臺風(fēng)造成中國平均每年死亡人數(shù)達(dá)472人, 直接經(jīng)濟(jì)損失達(dá)287億元, 占國內(nèi)生產(chǎn)總值的0.38% (Zhang, 2009)。分析海洋對臺風(fēng)路徑響應(yīng)過程及其相應(yīng)機(jī)理, 可以提高臺風(fēng)強(qiáng)度和路徑的預(yù)測精度, 這對減少地區(qū)經(jīng)濟(jì)損失, 保障人民生命和財產(chǎn)安全具有重要的社會意義。

臺風(fēng)過境時會引起上層海洋物理過程的相應(yīng)調(diào)整。在開放海域, 上層海洋對臺風(fēng)的響應(yīng)可分為兩個階段(管守德, 2014)。臺風(fēng)過境前的12～24 h被稱為上層海洋對臺風(fēng)響應(yīng)的“強(qiáng)迫階段”。在此階段, 臺風(fēng)快速變化的風(fēng)場會引起海洋上層產(chǎn)生強(qiáng)烈剪切不穩(wěn)定進(jìn)而引起強(qiáng)烈的垂向混合(卷挾)和混合層加深(Shay,1998; Sanford, 2011); 同時, 臺風(fēng)的正風(fēng)應(yīng)力旋度會引起海表水輻散, 引發(fā)局地水體上升運(yùn)動。臺風(fēng)過境后, 通常會在海洋中產(chǎn)生強(qiáng)烈的近慣性運(yùn)動, 臺風(fēng)輸入局地混合層的近慣性能量以近慣性內(nèi)波的形式向海洋溫躍層甚至海洋深層傳播(Chen, 2017; Mackinnon, 2017), 近慣性振蕩衰退時間一般一周至10 d以上(Qi, 1995; Teague, 2007), 此階段為上層海洋對臺風(fēng)響應(yīng)的另一個階段, 被稱為“松弛階段”(Brooks, 1983)。相對于開闊大洋, 近岸淺水海域受潮汐、徑流、地形等影響, 海洋的近慣性振蕩過程對臺風(fēng)的響應(yīng)更為復(fù)雜(Glenn, 2016; Seroka, 2016)。Millot等(1981)發(fā)現(xiàn)地中海西北部的獅子灣海域在陣風(fēng)持續(xù)作用后近岸上升流會同時伴隨近慣性振蕩的產(chǎn)生。De Freitas等(2019)在沿巴西陸架上的觀測中發(fā)現(xiàn), 近岸的近慣性振蕩向低緯趨于減小, 并且擁有較高的垂向剪切, 達(dá)到了10–3/s。Li等(2021)在南海北部發(fā)現(xiàn)臺風(fēng)引起的近慣性振蕩在跨陸坡分布上存在斜坡強(qiáng)化的現(xiàn)象。Igeta等(2011)發(fā)現(xiàn)在季風(fēng)和臺風(fēng)的共同作用下, 日本能登半島海域的沿岸捕獲波和近慣性振蕩疊加形成強(qiáng)沿岸流。理論上近岸區(qū)域潮下流伴隨的負(fù)渦度類似于鋒面的作用, 能夠捕獲近慣性內(nèi)波(Kunze, 1984), 但鮮有直接觀測的報道。

海表溫度對臺風(fēng)響應(yīng)的兩個階段也對應(yīng)著兩次降溫過程(Price, 1981)。由于臺風(fēng)路徑右側(cè)風(fēng)場和近慣性運(yùn)動共同作用, 臺風(fēng)路徑右側(cè)(南半球?yàn)樽髠?cè))的海洋流動與垂向混合強(qiáng)于左側(cè), 并引起更強(qiáng)的海洋上層降溫(Oey, 2008; Huang, 2015; Sun, 2015)。臺風(fēng)在開闊海洋中引起的海表溫度(sea surface temperature, SST)降溫幅度通常在1～6 °C, 最大可達(dá)9 °C (Lin, 2003; Black, 2008)。在近岸海區(qū), 受局地環(huán)境條件影響, 如水深急劇變淺和岸線阻擋等, SST對臺風(fēng)的響應(yīng)具有與開闊大洋不一樣的特征(賴巧珍等, 2013); 臺風(fēng)在近岸海域引起的最大SST降溫可能出現(xiàn)在臺風(fēng)路徑左側(cè)(Bingham, 2007; 梁曉紅等, 2014), 甚至臺風(fēng)過后出現(xiàn)海表升溫的現(xiàn)象(謝玲玲等, 2017), 謝玲玲等(2017)16認(rèn)為這可能是臺風(fēng)在外海引起的非線性海面孤立波向近岸的熱輸送造成。

臺風(fēng)過境會引起強(qiáng)烈的風(fēng)暴潮過程, 不同臺風(fēng)登陸所產(chǎn)生的風(fēng)暴潮具有很大差異。鄧國通等(2021)利用ROMS (regional ocean model system)模型研究了深圳近海臺風(fēng)所致風(fēng)暴潮的影響因素, 發(fā)現(xiàn)臺風(fēng)移動速度總體上對風(fēng)暴潮影響不大, 但不同登陸地點(diǎn)存在明顯差異。在深圳西邊登陸的臺風(fēng)比在深圳東邊登陸的臺風(fēng)產(chǎn)生的最大增水高1.5 m左右; 由東往西移動并登陸深圳的臺風(fēng), 比由南向北移動的臺風(fēng)產(chǎn)生的最大增水高 1.0 m左右。趙長進(jìn)等(2015)對影響長江口及其鄰近海區(qū)的風(fēng)暴潮影響因素進(jìn)行了研究, 發(fā)現(xiàn)向岸風(fēng)力的大小與持續(xù)時間對風(fēng)暴潮增水有著重要的作用, 并且風(fēng)應(yīng)力對增水的作用相較于氣壓更為顯著。

珠江口位于南海北部, 每年有發(fā)源于南海及西北太平洋地區(qū)的臺風(fēng)頻繁影響該區(qū)域, 或直接在珠江口登陸。臺風(fēng)影響珠江口海域時, 不僅會引起強(qiáng)烈的風(fēng)暴潮過程, 還會產(chǎn)生極端海流和波浪過程, 引起近海環(huán)流過程的急劇變化, 這會對泥沙、污染物等物質(zhì)在河口及近海的輸運(yùn)造成顯著影響, 并最終改變河口地形和引起近海地區(qū)浮游植物生長的短暫暴發(fā)。前人在風(fēng)暴潮、泥沙輸運(yùn)及初級生產(chǎn)力變化對臺風(fēng)響應(yīng)方面有相當(dāng)多的研究, 但上述研究主要是利用水位觀測資料、遙感數(shù)據(jù)和數(shù)值模型對海洋災(zāi)害及物質(zhì)輸運(yùn)和生態(tài)響應(yīng)過程進(jìn)行研究, 目前尚未見利用現(xiàn)場長期定點(diǎn)觀測浮標(biāo)數(shù)據(jù)對臺風(fēng)影響期間的海溫快速變化和淺海的近慣性振蕩過程開展研究。本文將基于浮標(biāo)觀測數(shù)據(jù), 研究珠江口雞啼門海域的海洋動力學(xué)和熱力學(xué)對臺風(fēng)的響應(yīng)。

1 數(shù)據(jù)和方法

1.1 浮標(biāo)觀測數(shù)據(jù)

依托珠江口海洋生態(tài)環(huán)境教育部野外科學(xué)觀測研究站, 中山大學(xué)海洋科學(xué)學(xué)院于2021年7月5日在珠江口雞啼門海域(水深10 m)投放了“1號”海洋環(huán)境觀測浮標(biāo)(圖1)。珠江口雞啼門位于珠海市金灣區(qū), 磨刀門以西約20 km, 最大漲潮差2.90 m, 最大落潮差2.71 m, 平均漲潮差為1.01 m。浮標(biāo)位于113°20.04′E, 21°55.34′N, 離高欄島4 km。浮標(biāo)上部的觀測平臺安裝有測風(fēng)儀、氣象儀、紅外等觀測氣象要素的傳感器, 下部裝載海流計、波浪儀等傳感器用于觀測海表波浪、水溫和流速。氣象觀測儀主要是對海表面高度2 m進(jìn)行采樣, 所有傳感器的采樣周期均設(shè)定為30 min, 獲取的氣象要素主要包括風(fēng)速風(fēng)向、空氣溫度、濕度、氣壓以及海表面溫度。

浮標(biāo)在7～8月期間共經(jīng)歷了臺風(fēng)“查帕卡”(Cempaka, 2021年7月18～24日, 7月20日登陸陽江市)和“盧碧”(Lupit, 2021年8月2～9日, 8月5日登陸汕頭市南澳島), 兩個臺風(fēng)的路徑見圖1。

1.2 再分析數(shù)據(jù)

1.2.1 海表面溫度數(shù)據(jù) 本文利用Remote Sensing System (RSS)提供的SST再分析資料研究了臺風(fēng)過境前后的SST變化。該產(chǎn)品時間分辨率是1 ｄ, 空間分辨率為0.25°×0.25°, 是由多個微波和紅外衛(wèi)星傳感器通過最優(yōu)插值融合得到。文中選取了2021年7月15～24日、7月30日到8月8日兩個臺風(fēng)經(jīng)過前后的數(shù)據(jù)。

1.2.2 臺風(fēng)路徑數(shù)據(jù) 臺風(fēng)路徑數(shù)據(jù)取自中央氣象臺臺風(fēng)網(wǎng)(typhoon.nmc.cn)發(fā)布的臺風(fēng)信息系統(tǒng), 數(shù)據(jù)包括臺風(fēng)的強(qiáng)度、時間、臺風(fēng)中心位置、中心氣壓、風(fēng)速等。臺風(fēng)移動速度是根據(jù)６ h內(nèi)臺風(fēng)移動的距離計算得到。

1.3 數(shù)據(jù)分析方法

基于浮標(biāo)觀測的風(fēng)速、空氣溫度、相對濕度、氣壓和水溫等數(shù)據(jù), 利用COARE3.0塊體算法對海氣熱通量進(jìn)行了計算; 再采用巴特沃茲濾波器對流速分量進(jìn)行帶通濾波, 提取海表近慣性頻段(0.85～1.2,為局地科里奧利頻率)的水平流速分布(u和v), 該頻段可以在包含近慣性振蕩的同時盡可能避免受到如全日潮等其他因素的影響。

圖1 浮標(biāo)位置及臺風(fēng)移動路徑

注: 圖中每個點(diǎn)代表時隔6 h臺風(fēng)的位置, 虛曲線代表水深(單位: m); S1、S2分別代表臺風(fēng)“查帕卡”和“盧碧”生成的地方

上層海洋近慣性動能(NIKE)可用式(1)進(jìn)行估計(張騫等, 2019):

2 結(jié)果與分析

2.1 海表響應(yīng)過程

2.1.1 海表流速與波浪的響應(yīng) 圖2展示了2021年7月5日至8月6日海表流速與氣象參數(shù)的變化。臺風(fēng)“查帕卡”加強(qiáng)階段(7月20日), 相對濕度增大, 氣壓降低, 這時浮標(biāo)位于臺風(fēng)的右側(cè), 風(fēng)向發(fā)生了順時針方向變化, 并一直持續(xù)到臺風(fēng)消失前一天。而海表流速幾乎保持向南流動, 當(dāng)海流方向與風(fēng)向一致時, 緯向流速快速增強(qiáng)到0.8 m/s (圖2a、2b的風(fēng)場、波浪為原始數(shù)據(jù), 其余數(shù)據(jù)經(jīng)過1 h低通濾波處理)。19日20點(diǎn), 由于風(fēng)向變化, 流速大小逐漸降低。臺風(fēng)“盧碧”期間, 浮標(biāo)在臺風(fēng)左側(cè), 相對濕度和氣壓變化不明顯, 測得風(fēng)場為逆時針變化, 風(fēng)向與海流方向一致, 最大緯向流達(dá)到0.5 m/s (圖2d)。

“查帕卡”經(jīng)過浮標(biāo)附近時, 有效波高從0.5 m突增至1.5 m, 值得注意的是, 在7月7日(圖2c), 有效波高突然增加, 而“盧碧”經(jīng)過浮標(biāo)附近時, 有效波高變化較小。這可能是因?yàn)樵诖箫L(fēng)過程中風(fēng)向持續(xù)保持北向, 南北走向岸線的沿岸海域的風(fēng)浪能持續(xù)通過風(fēng)應(yīng)力獲取能量累加而更加充分地成長; 而在“查帕卡”和“盧碧”過程中風(fēng)向有明顯的變向過程, 當(dāng)風(fēng)從陸域(特別是東側(cè))方向吹向海面時, 較短的風(fēng)區(qū)造成局部海域風(fēng)浪成長不充分, 此過程中, 波浪受到地形折射影響, 波向則仍以偏南向?yàn)橹?。因? 該海區(qū)強(qiáng)冷空氣過程對波浪影響可能大于臺風(fēng)過程。

2.1.2 SST對臺風(fēng)的響應(yīng) 兩次臺風(fēng)過程均引起浮標(biāo)所在位置SST顯著降溫。“查帕卡”期間SST最大降溫幅度為2.6 °C (與臺風(fēng)生成前3 d對比, 下同); 而“盧碧”生成的第一天和第二天SST先增溫了0.6 °C, 隨后逐漸降溫, 最大降溫幅度為2.2 °C。在臺風(fēng)過境時, 海表面空氣溫度(surface air temperature, SAT)和SST均呈下降趨勢, 且下降幅度大于SST(圖2e), 表明臺風(fēng)加強(qiáng)階段, 海氣之間的溫差在增大, 但差值取決于不同的臺風(fēng)特性, 比如“查帕卡”期間SST與SAT最大差值達(dá)到4.2 °C, 而“盧碧”期間差值為2.5 °C。

表1總結(jié)了臺風(fēng)過境期間與正常天氣情況的晝夜平均SST、潛熱通量、風(fēng)速以及濕度的對比。正常天氣是指臺風(fēng)生成前3天的天氣情況, 晝夜時間分別為6:00～20:00和20:00～6:00。由表1可知, 兩次臺風(fēng)過境期間夜間SST均高于白天(兩者之差為0.23和0.09 °C)。潛熱通量顯示, 白天海洋向大氣輸送的熱量高于晚上, 即白天海洋失熱更多, 這可能是白天低溫的主要原因。而在臺風(fēng)“盧碧”生成前, 白天海洋失熱更多, 但白天SST卻高于晚上, 表明此時海溫并非主要受海氣熱通量的影響, 在近岸淺水區(qū)域, 垂向混合、平流等因素也可能主導(dǎo)著海表溫度的變化(徐海明等, 2018; 趙一丁等, 2019; 易香妤等, 2021)。

圖2 浮標(biāo)觀測到的海氣要素時間序列

注: 灰色陰影表示臺風(fēng)事件, 黑色虛線方框表示一次大風(fēng)過程; a: 風(fēng)向(藍(lán)色箭頭表示方向, 向上為北); b: 流向(紅色箭頭表示方向, 向上為北); c: 有效波高(紅線)和波周期(藍(lán)線); d: 緯向流速(, 紅線)和經(jīng)向流速(, 黑線); e: 海表溫度SST(紅線)和海表面空氣溫度SAT(藍(lán)線); f: 相對濕度(紅線)和氣壓(藍(lán)線)

表1 臺風(fēng)過境期間與正常天氣情況的海氣要素日變化

Tab.1 Diurnal variation of ocean-atmosphere factors during the passage of typhoon and the normal weather conditions

另外, 潛熱通量與風(fēng)速、溫度、濕度等因素有關(guān), 前人的研究也表明: 潛熱通量隨風(fēng)速增加而增大, 潛熱通量與風(fēng)速的相關(guān)性大于與海溫的相關(guān)性(孫啟振等, 2010)。從圖3的潛熱通量與風(fēng)速分布的相關(guān)關(guān)系可知, 潛熱通量與風(fēng)速大小之間存在較明顯的線性趨勢。從表1計算得到的風(fēng)速和濕度變化來看, 臺風(fēng)過境期間白天風(fēng)速更強(qiáng), 夜間濕度更低, 這表明白天和夜間的潛熱通量差異主要與風(fēng)速大小有關(guān)。

基于再分析數(shù)據(jù), 圖4和圖5展示了臺風(fēng)期間整個海區(qū)SST的降溫情況?！安榕量ā鄙僧?dāng)天, 浮標(biāo)附近海域SST降溫幅度為1.35 °C (與臺風(fēng)生成前3天對比, 下同), 第二天降溫幅度為0.95 °C, 第三天降溫幅度為0.37 °C, 此時SST達(dá)到了最低值, 臺風(fēng)致使浮標(biāo)附近海域SST降溫了2.67 °C, 這與浮標(biāo)監(jiān)測數(shù)據(jù)(2.6 °C)幾乎一致?！氨R碧”生成的第一天和第二天(圖5a和5b), 浮標(biāo)附近的海域先增溫了1 °C左右, SST的增溫現(xiàn)象與之前的研究結(jié)果一致(謝玲玲等, 2017)12。

圖3 臺風(fēng)過境期間潛熱通量與風(fēng)速分布的相關(guān)關(guān)系

注: a為臺風(fēng)“查帕卡”期間, b為“盧碧”期間;為相關(guān)系數(shù)

圖4 “查帕卡”登陸過程海表面溫度變化分布圖

注: 圖中變化是指臺風(fēng)當(dāng)天與臺風(fēng)生成前3天對比; 實(shí)線框?yàn)榕_風(fēng)中心當(dāng)日位置

再分析數(shù)據(jù)表明“查帕卡”引起海表溫度的冷卻集中在路徑的右側(cè)(圖4), 這與以往的研究一致(Sun, 2015), 但“盧碧”期間, 浮標(biāo)的右側(cè)并不比左側(cè)降溫明顯。臺風(fēng)登陸當(dāng)天(圖5e), 左側(cè)海溫降低幅度與右側(cè)相當(dāng)。原因可能是受左側(cè)淺水深的影響, 海氣通量變化過程和風(fēng)浪攪拌作用使得海水垂向混合更加充分, 加上地形、徑流等因素的影響, 左側(cè)能量迅速耗散, 使表層海溫更容易降低。

2.2 臺風(fēng)引起的熱量及能量變化

2.2.1 海氣界面熱通量 為探究海表溫度的響應(yīng), 本部分探究了海氣界面熱通量的變化。2021年7月15日至8月6日, 珠江口雞啼門海域海氣界面的感熱通量變化范圍是-7.2～75 W/m2, 正值表示海洋向大氣輸送熱量。圖6表明, 雞啼門海域在觀測期間以海洋向大氣釋放熱量為主。兩個臺風(fēng)生成初期, 感熱通量與潛熱通量迅速增加至極大值, 這與海氣溫差、大氣濕度的變化及風(fēng)速的變化趨勢較為一致。臺風(fēng)期間海氣界面熱通量發(fā)生強(qiáng)烈交換, 海洋熱量向大氣輸送, 并且潛熱通量遠(yuǎn)大于該海域的感熱通量(“查帕卡”約為4.6倍, “盧碧”約為8.3倍), 表明臺風(fēng)期間海氣界面熱量交換主要貢獻(xiàn)來自于潛熱通量。

2.2.2 近慣性能量 風(fēng)應(yīng)力快速變化是激發(fā)海洋上層近慣性振蕩的主要因素(Sanford, 2011; Pallàs-Sanz, 2016)。如圖7a所示, 在“查帕卡”過境期間, 在慣性頻率附近(1.08)譜能量出現(xiàn)峰值, 量值與M2分潮能量相當(dāng), 此時的近慣性信號為順時針旋轉(zhuǎn); 在臺風(fēng)“盧碧”期間, 也在慣性頻率附近(1.07)能量出現(xiàn)峰值, 但是“盧碧”所激發(fā)的近慣性信號為逆時針旋轉(zhuǎn)(圖7b)。導(dǎo)致這一現(xiàn)象的主要原因是“查帕卡”期間浮標(biāo)位于熱帶氣旋的右側(cè), 風(fēng)場隨時間呈順時針變化, 由于在北半球, 近慣性運(yùn)動是順時針的, 因此, 順時針旋轉(zhuǎn)的風(fēng)場易與海水慣性運(yùn)動共振激發(fā)強(qiáng)烈的近慣性振蕩(Chen, 2015); “盧碧”期間浮標(biāo)位于熱帶氣旋的左側(cè), 風(fēng)場和流場隨時間呈逆時針變化。因此, “盧碧”所激發(fā)的近慣性信號是單純的海水運(yùn)動, 其變化頻率接近慣性頻率, 而非傳統(tǒng)意義上風(fēng)場與慣性運(yùn)動共振激發(fā)的近慣性振蕩。

圖5 “盧碧”登陸過程海表溫度變化分布圖

注: 圖中變化是指臺風(fēng)當(dāng)天與臺風(fēng)生成前3天對比, 實(shí)線框?yàn)榕_風(fēng)中心當(dāng)日位置

圖6 海氣界面潛熱和感熱通量對比圖

注: 灰色陰影表示臺風(fēng)事件

圖7 流速旋轉(zhuǎn)譜

注:為局地科里奧利頻率;2為太陰主要半日分潮的頻率; a為臺風(fēng)“查帕卡”期間, b為臺風(fēng)“盧碧”期間

圖8為近慣性振蕩頻段內(nèi)流速變化的時間序列。在“查帕卡”期間, 緯向流和經(jīng)向流變化不同步, 存在明顯的相位差, 其最大振幅也不同, 分別為8.4和13.1 cm/s, 造成這一現(xiàn)象的原因可能是, 在“查帕卡”期間風(fēng)速橢圓與近慣性流速橢圓長軸正交(圖9a), 故而緯向和經(jīng)向近慣性流在增長的過程中沒有同步發(fā)展, 導(dǎo)致其振幅和相位存在差異?！氨R碧”期間, 緯向流和經(jīng)向流變化趨于一致, 且其極值均為8 cm/s。因?yàn)椤氨R碧”激發(fā)的逆時針近慣性運(yùn)動單純由風(fēng)場驅(qū)動, 其風(fēng)速橢圓與近慣性流速橢圓傾角一致(圖9b), 故而緯向和經(jīng)向流是同步發(fā)展的。

近慣性動能通常隨時間呈指數(shù)狀增長和衰減, 可以利用e折時間尺度來衡量近慣性能量衰減速度。圖8b中的實(shí)心圓點(diǎn)分別代表能量最大值的位置和最大能量衰減至1/e的位置, 計算得出臺風(fēng)期間近慣性動能的e折時間尺度約為38 h和35.5 h, 稍大于局地慣性周期32.14 h, 與陸架及深遠(yuǎn)海相比(e折時間尺度在1周至10 d以上), 淺水深的河口近海區(qū)域的近慣性動能e折時間尺度明顯偏小。迅速衰減的原因可能是, 近海區(qū)域水深較淺, 沒有持續(xù)的近慣性強(qiáng)迫, 受底摩擦的影響, 近慣性能量會在一個慣性周期后快速地耗散。

3 結(jié)論

本文基于浮標(biāo)觀測資料以及再分析數(shù)據(jù), 分析了臺風(fēng)期間珠江口雞啼門海域的海洋動力學(xué)和熱力學(xué)響應(yīng)。臺風(fēng)過境引起緯向流速顯著增強(qiáng), 最大可達(dá)0.8 m/s?！安榕量ā逼陂g, 風(fēng)場隨時間呈順時針變化, 易與海水慣性運(yùn)動發(fā)生共振, 因而激發(fā)了強(qiáng)烈的順時針近慣性振蕩, 但其緯向流速與經(jīng)向流速的響應(yīng)不同步, 存在明顯的相位差。而在“盧碧”期間, 風(fēng)場為逆時針旋轉(zhuǎn), 激發(fā)了接近慣性頻率的逆時針海水運(yùn)動, 并非傳統(tǒng)意義上的近慣性振蕩。受當(dāng)?shù)販\水深的影響, 該區(qū)域的近慣性能量衰減周期小于38 h。在熱力學(xué)方面, 臺風(fēng)引起的海表溫度SST降溫幅度大于2.2 °C同時伴隨著潛熱通量的顯著變化, 這種變化可能是臺風(fēng)過境期間SST晝夜溫差翻轉(zhuǎn)的主要原因之一。

近岸水體對臺風(fēng)的響應(yīng)明顯與開闊大洋不同, 由于近海海洋熱通量及淺水近慣性振蕩還會受到地形、徑流等因素的影響, 其響應(yīng)機(jī)制復(fù)雜, 還需要結(jié)合數(shù)值模式及更多的觀測數(shù)據(jù)進(jìn)行進(jìn)一步分析。

圖8 浮標(biāo)觀測到的近慣性振蕩時間序列

注: 灰色陰影表示臺風(fēng)事件; 滑動時間為16 h; a: 近慣性流速, b: 近慣性動能

圖9 近慣性流速橢圓(藍(lán)線)、擬合橢圓(紅線)和風(fēng)速橢圓(黑線)

注: 圖中紅色虛直線代表近慣性橢圓長軸; 圓點(diǎn)代表旋轉(zhuǎn)橢圓的起點(diǎn); 三角箭頭代表旋轉(zhuǎn)橢圓的終點(diǎn); a為臺風(fēng)“查帕卡”期間, b為臺風(fēng)“盧碧”期間

鄧國通, 劉敏聰, 邢久星, 等, 2021. 深圳近海風(fēng)暴潮影響因素分析[J/OL]. 熱帶海洋學(xué)報: 1-13.

孫啟振, 陳錦年, 閆俊岳, 等, 2010. 2008年南海季風(fēng)爆發(fā)前后西沙海域海氣通量變化特征[J]. 海洋學(xué)報, 32(4): 12-23.

張騫, 廖光洪, 藺飛龍, 等, 2019. 南海西北部上層海洋對臺風(fēng)“杜蘇芮”的響應(yīng)分析[J]. 海洋學(xué)報, 41(7): 22-35.

易香妤, 董文杰, 李劭怡, 等, 2021. 中國東海黑潮海溫變化特征及成因分析[J]. 海洋預(yù)報, 38(3): 38-51.

趙一丁, 尹訓(xùn)強(qiáng), 宋亞娟, 等, 2019. 浪致混合對2016年北太平洋海表溫度季節(jié)性預(yù)測的影響[J]. 海洋學(xué)報, 41(3): 52-61.

趙長進(jìn), 葛建忠, 丁平興, 2015. 長江口及其鄰近海區(qū)無結(jié)構(gòu)網(wǎng)格風(fēng)暴潮預(yù)報系統(tǒng)的研制與分析[J]. 海洋科學(xué)進(jìn)展, 33(2): 182-194.

徐海明, 崔夢雪, 2018. 與冬季北太平洋大范圍海溫異常相聯(lián)系的海氣特征[J]. 大氣科學(xué)學(xué)報, 41(3): 330-343.

梁曉紅, 葛黎麗, 2014. 臺風(fēng)對江蘇近海海表溫度的影響分析與探討[J]. 水產(chǎn)養(yǎng)殖, 35(10): 37-42.

謝玲玲, 何超鳳, 李明明, 等, 2017. 瓊東上升流區(qū)海表溫度對臺風(fēng)過境的響應(yīng)[J]. 海洋科學(xué)進(jìn)展, 35(1): 8-19.

賴巧珍, 吳立廣, SHIE C L, 2013. 0908號臺風(fēng)“莫拉克”登陸過程中海表溫度變化特點(diǎn)及其對“莫拉克”的影響[J]. 熱帶氣象學(xué)報, 29(2): 221-234.

管守德, 2014. 南海北部近慣性振蕩研究[D]. 青島: 中國海洋大學(xué): 27-44.

BINGHAM F M, 2007. Physical response of the coastal ocean to Hurricane Isabel near landfall [J]. Ocean Science, 3: 159-171.

BLACK W J, DICKEY T D, 2008. Observations and analyses of upper ocean responses to tropical storms and hurricanes in the vicinity of Bermuda [J]. Journal of Geophysical Research: Oceans, 113(C8): C08009.

BROOKS D A, 1983. The wake of hurricane Allen in the western Gulf of Mexico [J]. Journal of Physical Oceanography, 13(1): 117-129.

CHEN S L, CHEN D Y, XING J X, 2017. A study on some basic features of inertial oscillations and near-inertial internal waves [J]. Ocean Science, 13(5): 829-836.

CHEN S L, HU J Y, POLTON J A, 2015. Features of near-inertial motions observed on the northern South China Sea shelf during the passage of two typhoons [J]. Acta Oceanologica Sinica, 34(1): 38-43.

DE FREITAS P, DE LUCA LOPES AMORIM F, MILL G N,, 2019. Observations of near-inertial oscillations along the Brazilian continental shelf break [J]. Ocean Dynamics, 69(10): 1203-1215.

GLENN S M, MILES T N, SEROKA G N,, 2016. Stratified coastal ocean interactions with tropical cyclones [J]. Nature Communications, 7: 10887.

HUANG S M, OEY L Y, 2015. Right-side cooling and phytoplankton bloom in the wake of a tropical cyclone [J]. Journal of Geophysical Research: Oceans, 120(8): 5735-5748.

IGETA Y, WATANABE T, YAMADA H,, 2011. Coastal currents caused by superposition of coastal-trapped waves and near-inertial oscillations observed near the Noto Peninsula, Japan [J]. Continental Shelf Research, 31(16): 1739-1749.

KUNZE E, SANFORD T B, 1984. Observations of near-inertial waves in a front [J]. Journal of Physical Oceanography, 14(3): 566-581.

LI R X, CHEN C S, DONG W J,, 2021. Slope-intensified storm-induced near-inertial oscillations in the South China Sea [J]. Journal of Geophysical Research: Oceans, 126(3): e2020JC016713.

LIN I, LIU W T, WU C C,, 2003. New evidence for enhanced ocean primary production triggered by tropical cyclone [J]. Geophysical Research Letters, 30(13): 1718.

MACKINNON J A, ZHAO Z X, WHALEN C B,, 2017. Climate process team on internal wave-driven ocean mixing [J]. Bulletin of the American Meteorological Society, 98(11): 2429-2454.

MILLOT C, CRéPON M, 1981. Inertial oscillations on the continental shelf of the Gulf of Lions-Observations and theory [J]. Journal of Physical Oceanography, 11(5): 639-657.

OEY L Y, INOUE M, LAI R,, 2008. Stalling of near-inertial waves in a cyclone [J]. Geophysical Research Letters, 35(12): L12604.

PALLàS-SANZ E, CANDELA J, SHEINBAUM J,, 2016. Mooring observations of the near-inertial wave wake of Hurricane Ida (2009) [J]. Dynamics of Atmospheres and Oceans, 76: 325-344.

PRICE J F, 1981. Upper ocean response to a hurricane [J]. Journal of Physical Oceanography, 11(2): 153-175.

QI H B, DE SZOEKE R A, PAULSON C A,, 1995. The structure of near-inertial waves during ocean storms [J]. Journal of Physical Oceanography, 25(11): 2853-2871.

SANFORD T B, PRICE J F, GIRTON J B, 2011. Upper-ocean response to Hurricane Frances (2004) observed by profiling EM-APEX floats [J]. Journal of Physical Oceanography, 41(6): 1041-1056.

SEROKA G, MILES T, XU Y,, 2016. Hurricane Irene sensitivity to stratified coastal ocean cooling [J]. Monthly Weather Review, 144(9): 3507-3530.

SHAY L K, MARIANO A J, JACOB S D,, 1998. Mean and near-inertial ocean current response to Hurricane Gilbert [J]. Journal of Physical Oceanography, 28(5): 858-889.

SUN J R, OEY L Y, CHANG R,, 2015. Ocean response to typhoon Nuri (2008) in western Pacific and South China Sea [J]. Ocean Dynamics, 65(5): 735-749.

TEAGUE W J, JAROSZ E, WANG D W,, 2007. Observed oceanic response over the upper continental slope and outer shelf during Hurricane Ivan [J]. Journal of Physical Oceanography, 37(9): 2181-2206.

ZHANG Q, WU L G, LIU Q F, 2009. Tropical cyclone damages in China 1983-2006 [J]. Bulletin of the American Meteorological Society, 90(4): 489-496.

RESPONSE TO TYPHOONS IN COASTAL WATERS AT JITIMEN IN ZHUJIANG RIVER ESTUARY

CHEN Xing-Wen1, 2, QIU Chun-Hua1, 2, 3, ZHANG Heng1, 2, LIU Dong4, SHI Hui1,5, MA Yong-Gui1, 2, WANG Dong-Xiao1, 2

(1. School of Marine Sciences, Sun Yat-sen University, Zhuhai 519082, China; 2. Pearl River Estuary Marine Ecosystem Research Station, Ministry of Education, Zhuhai 519082, China; 3. Key Laboratory of Marine Environmental Detection Technology and Application, Ministry of Natural Resources, Guangzhou510310, China; 4. Guangzhou Ruihai Ocean Technology Co., Ltd, Guangzhou 511400, China; 5. China Water Resources, Pearl River Planning, Surveying, and Designing Co. Ltd, Guangzhou 510610, China)

Understanding the oceanis response to typhoon is of great significance for typhoon forecasting. Due to the complexity of the dynamic structure, the response of coastal water to typhoon is still unclear. Based on observation data of a buoy (July—August 2021) deployed at Jitimen in Zhujiang (Pearl) River estuary, we analyzed the dynamical and thermodynamic responses of the ocean waters to typhoons Cempaka and Lubit passing the study areas during the observation. The two typhoons caused a significant increase in zonal velocity of water. Cempaka intensified remarkably the clockwise near-inertial oscillations, while Lupit motivated the counterclockwise seawater motion that close to inertial frequency. Affected by the local shallow water, the near-inertia kinetic energy dissipated rapidly after an inertial cycle in 32 hours. The cooling amplitude of sea surface temperature (SST) caused by the two typhoons was greater than 2.2 °C, which was apparently affected by latent heat flux. At the same time, the SST during the night was higher than that during the day after typhoon passing, which was also influenced by latent heat flux.

typhoon; near-inertial oscillation; sea surface temperature (SST); heat flux; Jitimen; the Zhujiang (Pearl) River estuary

* 自然資源部海洋環(huán)境探測技術(shù)與應(yīng)用重點(diǎn)實(shí)驗(yàn)室開放基金課題, MESTA-2020-A006號; 國家自然科學(xué)重大項(xiàng)目, 41890851號; 國家自然科學(xué)基金青年項(xiàng)目, 41506102號。陳杏文, 博士研究生, E-mail: chenxw78@mail2.sysu.edu.cn

王東曉, 博士生導(dǎo)師, 教授, E-mail: dxwang@mail.sysu.edu.cn

2022-01-01,

2022-03-01

P731

10.11693/hyhz20220100001