西北太平洋臺風對冷渦及葉綠素濃度的影響

王同宇，張書文,2，蔣 晨，劉 瀟，曾偉強，李 韜

西北太平洋臺風對冷渦及葉綠素濃度的影響

王同宇1，張書文1,2，蔣 晨1，劉 瀟1，曾偉強1，李 韜1

（1. 廣東省近海海洋變化與災害預警重點實驗室，廣東海洋大學海洋與氣象學院，廣東 湛江 524088;2. 區域海洋動力學與數值模擬功能實驗室，海洋科學與技術國家實驗室，山東 青島 266237）

【】研究2007—2017年間西北太平洋過境冷渦的11個臺風導致浮游植物生長和冷渦（CCE）變化現象。統計并計算出可能影響葉綠素（Chl-a）濃度變化的因素:臺風性質（強度、移動速度、強迫時間）、臺風前混合層厚度（MLD）、降雨量、海表面溫度（SST）、?？寺槲俾剩‥PV）和兩層約化重力模式下的渦動能（EKE），其中EPV和EKE分別代表上升流和湍流混合強弱。通過線性回歸分析發現，除臺風強度、SST與Chl-a濃度相關性不顯著（＞0.05），移動速度h、強迫時間、降雨和MLD、EPV、EKE與Chl-a均有顯著相關性（＜0.05），并建立了冷渦背景條件下的多元線性回歸模型：= 0.006 - 0.0381+ 0.025 72+ 0.023 83。浮游植物生長主要取決于上升流和湍流混合對營養鹽的輸送作用，慢而尺度大臺風意味著受臺風強迫時間長，足以超過地球自轉調整的時間則會引起強上升流（EPV）以及湍流混合輸送營養鹽，促進Chl-a濃度大幅度增加，強湍流混合同時也需要降雨抑制，避免破壞浮游植物光合作用，臺風前CCE區域MLD(＜25 m)與Chl-a呈現出正相關。

西北太平洋；臺風; 冷渦; 浮游植物生長; 多元線性回歸

西北太平洋是全球臺風最活躍的海域，每年形成約25個臺風，最頻繁出現時期為7~10月。大多數臺風經過此區域（16°―28° N，120°―138° E，圖1），水深大于4 500 m，面積1.21×10-6km2。其海域表現為中尺度渦特征（水平尺度50 ~ 500 km），主導上層海洋動能，比平均海洋動能大兩個量級[1]。

當臺風經過海洋，通過能量的交換發生劇烈的海氣相互作用，海洋表層熱量損失，強風應力向海洋輸入機械能，海表面溫度（SST）降溫并伴隨著浮游植物大量生長。SST降溫主要是由于臺風引起海洋內部垂直混合、夾卷、上升流作用，并將次表層富含營養鹽的冷水抬升至真光層，導致浮游植物大量生長[2-6]，而海氣通量對SST的降溫作用僅5%~15%[7]。中尺度渦可以調節臺風期間的海氣相互作用，尤其是氣旋式冷渦（CCE）。CCE加強，該區域表現出明顯的SST降溫和浮游植物大量生長，其熱力學結構不穩定，臺風過境引起的近慣性震蕩使得混合層內剪切混合強，改變渦旋內部結構，加強上升流，極易將次表層富含營養鹽的冷水抬升[8-12]，CCE引起強海表冷卻通過負反饋機制減弱臺風強度[13]；相反，反氣旋式暖渦（AE）內部較厚的混合層和下沉流則會抑制冷水抬升，SST降溫和浮游植物生長效果并不顯著[14]。對于移動速度較慢，長時間作用同一區域的臺風，海水內部達到?？寺胶?，引起強上升流，會導致CCE加強或新的CCE產生[15]。

通過觀測[16]和數值模擬[17]發現臺風強迫下的CCE區域會出現明顯的降溫和浮游植物大量生長。2003年，在西北太平洋11個臺風中只有2個臺風對SST和浮游植物產生顯著的影響[14]；2000―2008年期間，受到超強臺風影響，僅有20個CCE加強，產生兩個新CCE[9]；Shang等[18]統計1997―2009年南海只有8個臺風過境CCE均導致強烈的SST和浮游植物變化。由此可見，強降溫和浮游植物大量生長事件遠小于預期。雖然臺風性質[強度、移動速度（h）、尺度][4-5, 9, 19-22]和臺風前的海洋層化程度[2-3, 23]對于上層海洋的生態響應十分重要，但是關于不同性質的臺風經過CCE背景條件下對浮游植物生長的影響機制尚未清楚。本研究通過分析西北太平洋2007―2017年期間臺風對CCE及其浮游植物生長的影響，統計Chl-a濃度、SST對不同的臺風性質和CCE響應特征，建立多元回歸模型，以期為進一步揭示復雜背景下臺風與CCE之間生態響應機制提供基礎。

黑線為選取斷面

1 數據與方法

1.1 遙感和臺風數據

Chl-a濃度數據源自哥白尼海洋環境監測服務（CMEMS）提供的遙感產品，該產品融合了SeaWiFS、MODIS-Aqua、MERIS、VIIRSN and OLCI-S3遙感傳感器，空間分辨率為0.036° × 0.036°，時間分辨率為1 d，海表10 m風場數據從遙感系統網站下載（http://www.remss.com/），SST數據空間分辨率為0.083° × 0.083°，時間分辨率為1 d，風場空間分辨率為0.25° × 0.25°，時間分辨率為6 h。海表面高度異常（SLA）數據選自法國研究中心（AVISO）融合TOPEX/Poseidon、Jason-1、ERS-1/2和 Envisat 等多顆衛星數據產品，時間分辨率為1 d，空間分辨率為0.25° × 0.25°?；旌蠈雍穸龋∕LD）是CMEMS的全球觀測海洋溫度、鹽度、高度、MLD、地轉流、海表面鹽度和海表密度產品中下載（http://marine.copernicus.eu/），該數據產品結合了遙感和現場觀測資料[24]，時間分辨率為1周，空間分辨率為0.25° × 0.25°。氣候態硝酸鹽數據源自全球大洋數據集（WOA13），氣候態真光層深度數據從遙感系統網站下載。

臺風路徑、最大風速來自上海臺風所（http://www.typhoon.org.cn/），時間間隔為6 h，h是根據臺風每6 h移動距離計算。

1.2 方法

首先選取在2007—2017年經過該區域的38個臺風，主要集中在7—10月份，并篩選出過境后引起顯著SST異常和高濃度Chl-a的臺風，統計符合條件的臺風引起SST、Chl-a濃度變化情況以及其沿著路徑U、強度變化特征，同時，根據SLA數據來識別CCE，并觀察浮游植物大量生長的區域是否在CCE附近。其中為了突出Chl-a濃度變化，從過境后的兩個星期內選取處高濃度Chl-a所出現的日期，合成平均。臺風強迫時間根據10 m風場計算，將6 h分辨率插值成為30 min，并選取風速＞14 m/s風場范圍作為臺風強迫區域，并對強迫區域時間進行積分得到強迫時間。

臺風引起海洋內部?？寺槲俾剩‥PV）:

上升流引起的等密度面位移可由以下兩個公式計算:

其中，?為海表面高度變化，為重力加速度，′ =（－）是約化重力，為混合層平均密度，為混合層底至水深200 m平均密度，本文取0.015 m/s2。公式（4）、（5）計算出結果雖為同一量級，但數值大小不一致，是（4）基于風應力旋度導致均質上層海洋內部等密度面位移[26-27]，公式（5）則考慮真實海洋表現為中尺度渦特征，海洋內部為第一斜壓模態，兩層約化重力模式可以模擬在實際風應力作用下的自由表面高度變化隨時間的演變[2, 28]，所以，本文將采取兩層約化重力模式。

CCE變化則一方面由SLA來表征，另一方面，通過計算渦動能（EKE）可以反映地轉流場以及中尺度渦強弱：

其中為地轉流異常緯向分量，為地轉流異常緯向分量，A為CCE區域網格點數，h為海洋上層厚度，由全球大洋數據集（WOA13）月平均海洋溫度計算（https://www.nodc.noaa.gov/OC5/woa13/），水平分辨率為0.25o×0.25o，垂直方向從0~1500m共有57層。

如圖2所示，西北太平洋區域西北部島嶼眾多，特殊的地理環境使得表層Chl-a濃度四季分布以島鏈為軸，向東南方向遞減。冬季平均葉綠素質量濃度0.1936 mg/m3高于春（0.1594）、夏（0.1530）、秋（0.164 1）季，但在離岸海域Chl-a質量濃度基本上≤0.1 mg/m3。篩選出2007-2017年11個臺風經過CCE，并引起大范圍的浮游植物生長和SST降溫，表1列出11個臺風過境CCE期間的最大風速、h、Chl-a濃度、SST變化特征，表2列出臺風強迫時間、降雨、MLD、EPV、EKE特征。

圖2 2007―2016年西北太平洋典型季節Chl-a濃度空間分布

表1 臺風過境CCE區域臺風特征及Chl-a濃度和SST變化

注：臺風名稱出現兩次表示該臺風過境兩個CCE區域

Note: The name of the typhoon appeared twice, indicating that the typhoon passed through two CCE region

表2 CCE區域臺風強迫時間、降雨、MLD、EPV、EKE特征

注：臺風名稱出現兩次表示，該臺風過境兩個CCE區域

Note: The name of the typhoons appeared twice, indicating that the typhoon passed through two CCE region

2 結果與分析

2.1 特殊路徑臺風盧碧和派比安對Chl-a、SST和SLA的影響

臺風盧碧和派比安的路徑是比較特殊的，盧碧向西運動經過兩個強度相當的CCE（圖3、4），經過第一個CCE時h為5 m/s，強度為60 m/s，降溫幅度達到5.3 ℃（表1），高Chl-a濃度分布在CCE內臺風路徑附近（圖3：b），增長0.22 mg/m3，當經過第二個CCE區域后，路徑轉向東北方向，h為2 m/s，降溫為6.9 ℃，Chl-a質量濃度大幅度增加，增長0.62 mg/m3（圖3：c），另外，在臺風過后1周之后，CCE均加強，并向西傳播，下降幅度-6 cm（圖4：f）。

在派比安臺風過境前，其路徑周圍CCE較弱，SLA=-15 cm。過境時，h分別為1.0、0.8 m/s，強度分別為46、38 m/s，在兩個打轉區域形成一個東北-西南帶狀的強CCE（SLA=-52 cm）（圖4：g-h），同時第一個打轉區域臺風強迫時間大約2 d，SST降溫幅度達到7.3 ℃，Chl-a質量濃度增加0.4 mg/m3，且范圍較小，第二個區域臺風強迫時間＞3 d，SST降溫6.3 ℃，高Chl-a濃度范圍遠大于前者，并分布在路徑右側，增加1.15 mg/m3，浮游植物生長一直維持到第二周（圖3：e-f）。

2.2 連續過境兩個CCE的臺風蘇力和黃蜂對Chl-a、SST和SLA的影響

臺風蘇力和黃蜂路徑均經過2個CCE，均呈現出高濃度Chl-a（圖5），這兩個CCE的SLA變化最為顯著，范圍擴大（圖6）。

蘇力過境第一個CCE時，h為6.1 m/s，強度為50 m/s，導致降溫和Chl-a質量濃度為：6.8 ℃、0.12 mg/m3，高濃度范圍小，CCE強度幾乎未變化；過境第二個CCE時，U為5.8 m/s，強度44 m/s，降溫4.3 ℃，Chl-a質量濃度增長至0.24 mg/m3，CCE內部SLA由-36 cm加強至-41 cm。

黃蜂過境以5.1 m/sh和55 m/s最大風速經過第一個CCE，導致降溫4.5℃和小范圍的Chl-a質量濃度增加0.13 mg/m3，CCE范圍擴大，向西運動，強度變化微弱；而經過第二個CCE時，h為2.5m/s，最大風速達到62 m/s，造成CCE發生顯著加強，SLA減少了16 cm（圖5：c-d），降溫幅度＞7 ℃，Chl-a質量濃度增加0.31 mg/m3（圖5：e-f）。

黃蜂過境以5.1 m/s 的h和55 m/s的最大風速經過第一個CCE，導致降溫4.5 ℃和小范圍的Chl-a質量濃度增加0.13 mg/m3，CCE范圍擴大，向西運動，強度變化微弱；而經過第二個CCE時，h為2.5 m/s，最大風速達到62 m/s，造成CCE發生顯著加強，SLA減少了16 cm（圖5：c-d），降溫幅度＞7 ℃，同時Chl-a質量濃度增加0.31 mg/m3（圖5：e-f）。

2.3 過境一個CCE的臺風對Chl-a、SST和SLA的影響

除上述4個臺風外，還有7個臺風：羅莎、森克拉、薔薇、莫拉克、杰拉華、巴蓬、鯰魚，其中森克拉和薔薇連續經過同一個CCE，鯰魚、杰拉華、莫拉克路徑位于CCE中心，其余路徑在CCE邊緣附近，用相同的方法計算SST、Chl-a、SLA變化，結果如圖7所示。

（a-c）：盧碧臺風；（d-f）：派比安；（a）：盧碧臺風過境前、（b）：過境第一個CCE后、（c）：第二個CCE后；（d）：派比安臺風過境前、（e）：過境第一個CCE后、（f）：第二個CCE后.紅線代表臺風路徑，黑線代表CCE內SLA等值線

a、b、e、f：盧碧臺風；c、d、g、h：派比安臺風；a和c：過境第一個CCE后SST降溫；b和d：第二個CCE后SST降溫；e和g：過境第一個CCE后SLA；f和h：第二個CCE后SLA；粗黑線代表臺風路徑，細黑線代表CCE內SLA等值線

（a-c）：蘇力臺風；（d-f）：黃蜂臺風；（a，d）：臺風過境前Chl-a濃度、（b，e）：臺風過境后Chl-a濃度、（c，f）：臺風過境后海表溫度降溫；粗黑線代表臺風路徑，細黑線代表CCE內SLA等值線

（a，b）：蘇力臺風；（c，d）：黃蜂臺風；（a，c）：臺風過境前SLA、（b，d）：臺風過境后SLA.粗黑線代表臺風路徑，細黑線代表冷渦內SLA等值線

7個臺風過境CCE時最大風速＞42 m/s，除森拉克h＜3 m/s，其他臺風h較快，連續臺風森克拉和薔薇過境CCE位于近岸區域，臺風前CCE區域Chl-a質量濃度分別為0.11、0.18 mg/m3，與秋季該區域Chl-a濃度大致相同，但臺風過境后Chl-a質量濃度增加幅度為0.30、0.23 mg/m3，且分布在CCE不同區域，森克拉SST降溫（5.7 ℃）和范圍比薔薇顯著（2 ℃）（圖7：d-i）。鯰魚、杰拉華、莫拉克h為7.2、4.8、6.0 m/s，最大風速為46、57、43 m/s，莫拉克h和強度弱相對較弱，而且臺風前CCE較弱（SLA=-10cm），過境后CCE強度和范圍加強，但浮游植物生長卻出現在CCE區域之外，而CCE內Chl-a濃度增加0.12 mg/m3，SST降溫3.8 ℃; 鯰魚過境后，SST降溫3.9 ℃，與莫拉克相近，但引起Chl-a質量濃度增加了0.3 mg/m3，分布在CCE內部，SLA減少至-32 cm; 杰拉華過境后導致CCE發展成為沿臺風路徑的西北東南帶狀特征，SLA從-28 cm減少至-37 cm，高Chl-a濃度也呈現出帶狀結構，增加了0.46 mg/m3。羅莎和巴蓬臺風最大風速相近（46、48 m/s），羅莎h（4.2 m/s）大于巴蓬（5.2 m/s），羅莎過境前，CCE呈現出東北―西南帶狀特征，而過境后CCE中心則位于臺風路徑周圍，高濃度Chl-a卻出現在CCE邊緣，增加0.56 mg/m3，SST降溫6 ℃；巴蓬過境前，CCE強度弱（SLA=-13 cm），過境后SLA減少至-18 cm（圖8），SST降溫也達到了5.6 ℃，CCE區域內只出現小范圍高濃度Chl-a，增加了0.26 mg/m3。

第1-7行分別是羅莎（a-c）、森克拉（d-f）、薔薇（g-i）、莫拉克（j-l）、杰拉華（m-o）、巴蓬（p-r）、鯰魚（s-u）臺風過境前、后Chl-a濃度和SST降溫。粗黑線代表臺風路徑，細黑線代表CCE內SLA等值線。

羅莎（a-b）、森克拉（c-d）、薔薇（e-f）、莫拉克（g-h）、杰拉華（i-j）、巴蓬（k-l）、鯰魚（m-n）臺風過境前、后SLA變化。粗黑線代表臺風路徑，細黑線代表CCE內SLA等值線

2.4 Chl-a與多個變量的線性回歸

11個臺風經過15個CCE時，引起的不同程度的浮游植物生長，臺風強度、h、強迫時間、降雨、MLD、EPV、EKE變量如表1和2所示，通過線性回歸分析（表3），表明Chl-a濃度與7個變量之間表現出不同的相關性：h（= 0.67，= 0.006）、強迫時間（= 0.81，= 0.000 3）、降雨（= 0.69，= 0.004 8）、MLD（= 0.690 9，= 0.043）、EPV（= 0.515 0，= 0.049 6）、EKE（= 0.928 6，< 0.000 1）與Chl-a有著顯著相關性（< 0.05），但臺風強度（= -0.266 2，= 0.337 5）和SST（= 0.40，= 0.141）相關性不顯著（> 0.05），并未通過95%置信區間的檢驗（圖9）。

進一步建立多元線性回歸模型，即Chl-a變化與h、降雨、MLD自變量之間的回歸模型：=0+11+22+33+，為隨機誤差，1、2、3分別是自變量h、MLD、降雨量，求解結果如表3所示，復相關系數為2= 0.56，= 0.02 < 0.05，該回歸模型成立。

表3 回歸模型系數和置信區間

（a-h）表示Chl-a與影響因子（Uh、強迫時間、最大風速、ΔSST、MLD、降雨、EKE和EPV）

3 討論

在上層海洋中，光和營養鹽是促進浮游植物大量生長最重要的兩個因素，光照條件和營養鹽濃度分布可以調節Chl-a濃度變化[26]。西北太平洋篩選出的11個臺風主要發生8―10月份，沿著研究區域選取西北―東南向一個斷面，密度躍層基本維持50 m，氣候態營養鹽濃度隨著深度逐漸增加，在上層0~200 m硝酸鹽濃度均＜4 μmol/L，這說明西北太平洋上層海洋營養鹽匱乏，而氣候態真光層深度在90~120 m之間，真光層內雖可以提供給浮游植物光照條件，但低營養鹽濃度也會抑制浮游植物生長，整個西北太平洋開闊洋面上Chl-a質量濃度≤0.1 mg/m3。2007―2017年臺風中僅有11個臺風過境該海域可以導致Chl-a濃度顯著增加，維持7 d，均發生在CCE區域，說明僅有臺風作用不足以將底層高濃度營養鹽攜帶進入真光層，還需要不穩定的熱力學結構才能引起浮游植物生長。

臺風后的CCE區域Chl-a濃度總是高于臺風前，臺風本身性質（h、強度、尺度）對于Chl-a濃度增加和SST降溫程度起著重要作用，主要通過臺風引起的上升流、混合夾卷過程，一方面，次表層高濃度的Chl-a向上輸送至表層，另一方面，底層富含營養鹽的冷水向上輸送，促進浮游植物生長，上升流和湍流混合的強弱決定了為浮游植物提供營養鹽輸送的能力[27-30]，以EPV、EKE代表上升流和湍流混合強度，臺風引起上升流的條件滿足公式=h/（2max）＞1[7]，EPV與Chl-a濃度顯著相關（= 0.515 0，= 0.049 6）；＜1，意味著慢而尺度較大的臺風過境的上層海洋以湍流混合為主，EKE與Chl-a濃度相關性接近于1（= 0.928 6，＜0.000 1）。強迫時間計算過程考慮到臺風尺度和U兩個因素，對EPV和EKE產生重要影響。

統計分析發現，強迫時間與Chl-a濃度相關性顯著（= 0.81，= 0.000 3），臺風強度與Chl-a濃度變化之間的相關性弱（= -0.266 2，= 0.337 5），說明決定Chl-a增加程度的主要因素是臺風強迫時間，慢而尺度大的臺風意味著向海洋輸入更多的能量，這些能量主要以近慣性內波的形式向下傳播，造成上層海洋內部的剪切不穩定機制，湍流混合過程使得混合層加深，極易將高濃度營養鹽抬升至真光層內[11, 17]。SST是Chl-a濃度變化的一個重要指標，SST與臺風性質也有一定相關性[7-9]，但SST與Chl-a濃度（= 0.40，= 0.141）未表現出強相關性，可能是SST除海洋內部動力過程外，還受到臺風期間降雨、太陽輻射蒸發、臺風前CCE的層結結構等影響。臺風前15個CCE區域MLD范圍在10~25 m之間，遠小于氣候態MLD（50 m），其上層海洋易受到上升流、湍流混合以及夾卷過程影響，利于MLD加深和營養鹽輸送，但臺風前的MLD與Chl-a濃度表現為正相關（= 0.690 9，= 0.043），說明CCE的熱力學結構越不穩定，強混合夾卷過程會將部分浮游植物位移至真光層之下，破壞光合作用，所以臺風前CEE存在較厚的MLD易出現更高濃度Chl-a[31]。CCE區域累計降雨量與Chl-a濃度表現出顯著正相關（= 0.69，= 0.004 8），主要因為臺風長強迫時間會導致大量降雨，注入海洋出現強分層現象，抑制強湍流混合將更多浮游植物位移至真光層之下的過程，避免破壞光合作用，從而使得浮游植物大量生長[23]。實際上，營養水平因地點而異，營養水平變化對Chl-a濃度增加的影響將在今后的研究中進行，并考慮詳細的情況。

4 結論

本研究基于遙感數據統計了近10 a西北太平洋臺風經過CCE區域引起浮游植物生長程度和CCE變化，發現此區域內氣候態的真光層內營養鹽匱乏＜4 μmol/L，營養鹽成為控制浮游植物生長重要因素，上升流和湍流混合強弱對于營養鹽輸送至關重要。在臺風過境CCE區域的背景條件下，除臺風強度、SST外，h、強迫時間、降雨和MLD、EPV、EKE與Chl-a濃度表現出顯著相關性，建立Chl-a變化與h、MLD、降雨量自變量之間的多元線性回歸模型：= 0.006 - 0.0381+ 0.025 72+ 0.023 83。慢而尺度大的臺風意味著較長的強迫時間，便于形成強上升流，向輸入更多能量，產生強湍流混合，攜帶營養鹽輸送至真光層內；而在熱力結構不穩定的上層海洋，MLD、降雨量與Chl-a變化表現為正相關。

[1] XU C, SHANG X D, HUANG R X. Horizontal eddy energy flux in the world oceans diagnosed from altimetry data[J]. Sci Rep, 2014, 4: 5316.

[2] NAN D W, LEBEN R R, BALASUBRAMANIAN S. Hurricane-forced upwelling and chlorophyll-a enhancement within cold‐core cyclones in the Gulf of Mexico[J]. Geophysical Research Letters, 2005, 32(18): 109-127.

[3] ZHENG Z W, HO C R, KUO N J. Importance of pre-existing oceanic conditions to upper ocean response induced by Super Typhoon Hai-Tang[J]. Geophysical Research Letters, 2008, 35(20): 288-299.

[4] 林靜柔, 唐丹玲, 婁全勝. 超級臺風 “南瑪都” 對南海北部葉綠素a、溫鹽及溶解氧的影響[J]. 生態科學, 2015, 34(4): 9-14.

[5] 楊元建, 傅云飛, 孫亮, 等. 基于多衛星和 Argo 浮標觀測海洋上層對臺風婷婷的響應[J]. 中國科學技術大學學報, 2010, 40(1): 1-7.

[6] FOLTZ G R, BALAGURU K, LEUNG L R. A reassessment of the integrated impact of tropical cyclones on surface chlorophyll in the western subtropical North Atlantic[J]. Geophysical Research Letters, 2015, 42(4): 1158-1164.

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

[8] ZHENG Z W, HO C R, ZHENG Q, et al. Effects of preexisting cyclonic eddies on upper ocean responses to category 5 typhoons in the western North Pacific[J]. Journal of Geophysical Research Oceans, 2010, 115(C9): 1-2.

[9] SUN L, LI Y, YANG Y, et al. Effects of super typhoons on cyclonic ocean eddies in the western North Pacific: A satellite data-based evaluation between 2000 and 2008[J]. Journal of Geophysical Research Oceans, 2015, 119(9): 5585-5598.

[10] JAIMES B, SHAY L K. Near-inertial wave wake of Hurricanes Katrina and Rita over mesoscale oceanic eddies[J]. Journal of Physical Oceanography, 2009, 40(6): 1320-1337.

[11] JAIMES B, SHAY L K. Enhanced wind-driven down welling flow in warm oceanic eddy features during the intensification of tropical cyclone isaac(2012): observations and theory[J].Journal of Physical Oceanography, 2015, 45(6): 1667-1689.

[12] MEI W, LIEN C C, LIN I I, et al. Tropical cyclone-induced ocean response: a comparative study of the south china sea and tropical Northwest Pacific[J]. Journal of Climate, 2015, 28(15): 150424113454002.

[13] SHADE L R. The ocean's effect on the intensity of tropical cyclones: results from a simple coupled atmosphere-ocean model[J]. Journal Atmosphere Science, 1999, 56(4): 642-651.

[14] LIN I I. Typhoon-induced phytoplankton blooms and primary productivity increase in the western North Pacific subtropical ocean[J]. Journal of Geophysical Research Oceans, 2012, 117(C3): C03039.

[15] HIROSHI K. Detection of cyclonic eddy generated by looping tropical cyclone in the northern South China Sea: a case study[J]. Acta oceanologica sinica, 2004, 23(2): 213-224.

[16] LIU F, TANG S. Influence of the interaction between typhoons and oceanic mesoscale eddies on phytoplankton blooms[j]. Journal of geophysical research oceans, 2018, 123(2): 241-256

[17] JAIMES B, SHAY L K, HALLIWELL G R. The response of quasigeostrophic oceanic vortices to tropical cyclone forcing[J]. Journal of Physical Oceanography, 2011, 41(10): 1965-1985.

[18] SHANG X D, ZHU H B, CHEN G Y, et al. Research on cold core eddy change and phytoplankton bloom induced by typhoons: case studies in the South China Sea[J]. Advances in Meteorology, 2015, 2015:42-60.

[19] VINENT E M, LENGAIGNE M, MADEC G, et al. Processes setting the characteristics of sea surface cooling induced by tropical cyclones[J]. Journal of Geophysical Research: Oceans, 2012, 117(C2).

[20] WALKER N D, LEBEN R R, PILLEY C T, et al. Slow translation speed causes rapid collapse of Northeast Pacific Hurricane Kenneth over cold core eddy[J]. Geophysical Research Letters, 2014, 41(21): 7595-7601.

[21] LU Z, WANG G, SHANG X. Response of a preexisting cyclonic ocean eddy to a typhoon[J]. Journal of Physical Oceanography, 2016, 46(8): 2403-2410.

[22] PAN J, HUANG L, DEVLIN A T, et al. Quantification of typhoon-induced phytoplankton blooms using satellite multi-sensor data[J]. Remote Sensing, 2018, 10(2): 318.

[23] LIN Y C, OEY L Y. Rainfall-enhanced blooming in typhoon wakes[J]. Scientific reports, 2016, 6: 31310.

[24] GUINEHUT S, DHOMPS A L, LARNICOL G, et al. High resolution 3-D temperature and salinity fields derived from in situ and satellite observations[J]. Ocean Science, 2012, 8(5): 845-857.

[25] LARGE W G, MCWILLIAMS J C, DONEY S C. Oceanic vertical mixing: A review and a model with a nonlocal boundary layer parameterization[J]. Reviews of Geophysics, 1994, 32(4): 363-403.

[26] FURUYA K. Subsurface chlorophyll maximum in the tropical and subtropical western Pacific Ocean: Vertical profiles of phytoplankton biomass and its relationship with chlorophylla and particulate organic carbon[J]. Marine Biology, 1990, 107(3): 529-539.

[27] BABIN S M, CARTON J A, DICKEY T D, et al. Satellite evidence of hurricane-induced phytoplankton blooms in an oceanic desert[J]. Journal of Geophysical Research: Oceans, 2004, 109(C3): 1938-1942.

[28] PRICE J F, SANFORD T B, Forristall G Z. Forced stage response to a moving hurricane[J]. Journal of Physical Oceanography, 1994, 24(2): 233-260.

[29] SHAY L K, GONI G J, BLACK P G. Effects of a warm oceanic feature on Hurricane Opal[J]. Monthly Weather Review, 2000, 128(5): 1366-1383.

[30] ZHAO H, HAN G, ZHANG S, et al. Two phytoplankton blooms near Luzon Strait generated by lingering Typhoon Parma[J]. Journal of Geophysical Research: Biogeosciences, 2013, 118(2): 412-421.

[31] GRAN H H, BRAARUD T. A quantitative study of the phytoplankton in the Bay of Fundy and the Gulf of Maine(including observations on hydrography, chemistry and turbidity)[J]. Journal of the Biological Board of Canada, 1935, 1(5): 279-467.

Cold Core Eddy Change and Phytoplankton Bloom Induced by Typhoons: Case Studies in Northwest Pacific Ocean

WANG Tong-yu1, ZHANG Shu-wen1,2, JIANG Chen1, LIU Xiao1, ZENG Wei-qiang1, LI Tao1

(1.，，524088，; 2.，，266237，)

【】To investigate the cold core eddy（CCE）change and phytoplankton bloom induced by 11 typhoons in northwest Pacific Ocean during the period 2007—2017. 【】The factors that might affect the increase of chlorophyll-a（Chl-a）concentration were calculated: typhoon properties（intensity，translation speedand forcing time)，pre-typhoon mixed layer depth（MLD），cumulative rainfall，sea surface temperature（SST），Ekman pumping velocity（EPV）and eddy kinetic energy（EKE）was estimated from two-layer reduced gravity model. EPV and EKE represent the strength of the upwelling and turbulent mixing，respectively. 【】Linear regression analysis shows that，except for typhoon intensity and SST（> 0.05），h，forcing time，rainfall，MLD，EPV，EKE were significantly correlated with Chl-a under the pre-existing CCE（< 0.05），and establish multiple linear regression model:= 0.006 -0.0381+ 0.02572+ 0.02383was established. Phytoplankton bloom mainly depended on nutrients transport induced by the upwelling and turbulent mixing. Slow typhoons meant that the forcing time was long enough to exceed the geostrophic adjustment time and strong upwelling（EPV）and turbulent mixing greatly increased Chl-a. In the meantime，due to the fact that MLD in CCE region was generally less than 25 m，strong turbulent mixing also needed to cooperate with rainfall inhibition to avoid damaging phytoplankton photosynthesis.

Northwest Pacific Ocean; typhoons; cold core eddy; phytoplankton bloom;multiple linear regression

P733

A

1673-9159（2019）05-0085-11

10.3969/j.issn.1673-9159.2019.05.013

2019-04-20

國家重點研發計劃重點專項（2016YFC1401403）；國家自然科學基金面上項目（41676008和41876005）；國際合作項目（GASI-IPOVI-04）；廣東省自然科學基金（2016A030312004）

王同宇（1995－），男，碩士研究生，研究方向為物理海洋。E-mail：gdou_wty@163.com

張書文（1962－），男，博士，教授，主要從事物理海洋學研究。E-mail：gdouzhangsw@163.com

王同宇，張書文，蔣晨，等. 西北太平洋臺風對冷渦及葉綠素濃度的影響[J]. 廣東海洋大學學報，2019，39（5）：85-95.

（責任編輯：劉朏）