黑潮入侵東海陸架途徑及其影響研究進展*

楊德周 尹寶樹, 侯一筠, 孫 松, 俞志明,

宋金明1, 3, 4, 5 許靈靜1, 2, 5 于 非1, 2, 3, 5 宋秀賢1, 3, 4, 5 于仁成1, 3, 4, 5袁華茂1, 3, 4, 5 李新正1, 3, 4, 5 袁涌銓1, 3, 4

(1. 中國科學院海洋研究所 青島 266071; 2. 中國科學院海洋環流與波動重點實驗室 青島 266071; 3. 青島海洋科學與技術國家實驗室 青島 266237; 4.中國科學院海洋生態與環境科學重點實驗室 青島 266071; 5. 中國科學院大學 北京100049)

在全球變化和人類活動共同影響下, 近年來我國近海生態系統動蕩加劇, 生物資源衰退, 赤潮、綠潮、水母旺發等生態災害頻發, 危及社會經濟發展和近海生態安全。國際上對海洋生態系統的研究主要圍繞著環境變化、生態系統演變、生態安全以及它們之間的關系展開(Chavezet al, 2003)。發起并組織實施的相關大型研究計劃包括以研究赤潮為主的“全球有害藻華生態學與海洋學研究(GEOHAB)”和“有害藻華與環境響應(HARRNESS)”; 強調陸海相互作用對生態系統影響的“海岸帶陸海相互作用研究(LOICZ)”; 以及全面研究生態系統的“全球海洋生態系統動力學研究(GLOBEC)”、“海洋生物地球化學與生態系統整合研究(IMBER)”等。以上研究計劃加深了對近海生態環境現狀的了解, 也為理解生態環境變化機制提供了基礎。

我國近海生態環境領域已經開展了大量的調查和研究工作, 對我國近海環境演變和生態系統的基本結構與功能的研究有了一定的科研積累。例如, “渤海生態系統動力學與生物資源持續利用”重大基金項目和“東、黃海生態系統動力學與生物資源可持續利用”國家“973”計劃項目等, 研究主要集中在陸源輸入、人類活動和近海本身的變化, 但對大洋輸入影響近海的過程認識不足, 對環境演變的主要驅動因子和生態系統變化規律缺乏系統的了解和認識。

我國鄰近海域處于北太平洋環流系統中, 外海環流及其變異是影響我國近海生態系統的穩定性和演變趨勢的重要驅動器之一。在北太平洋熱帶和副熱帶海域, 赤道東風帶和西風帶的存在導致了風應力負渦度的輸入, 根據 Sverdrup關系(Pedlosky, 1996; Huang,2010), 該負渦度的輸入必然導致太平洋副熱帶海域海水向赤道方向的水體輸送; 根據質量守恒和渦度守恒,向南輸送的水體必然通過西邊界流——黑潮, 向北輸送。由于東海大陸架和沖繩海槽的存在, 向北輸送的海水在臺灣東北進入沖繩海槽后為保持位渦守恒其流動路徑基本平行于東海陸架的大陸坡, 從而黑潮主軸主要出現在200m等深線附近。然后, 在30°N附近黑潮主軸通過吐噶喇海峽流出沖繩海槽, 離開東海進入太平洋, 從而構成了熱帶、副熱帶海域閉合的海洋環流系統。黑潮源源不斷地把熱帶太平洋海水向北輸送到副熱帶海區, 攜帶了大量的熱能、機械能和營養鹽,對東海陸架的環流和生態系統具有重要的調控作用(Su, 1998; Isobe, 2008; Guoet al, 2012; Yanget al, 2012;Guoet al, 2013; Yanget al, 2013)。黑潮帶來的高溫、高鹽的表層水及高磷酸鹽的次表層水, 可以改變近海的水團結構, 改變近海生源要素, 從而影響近海的生物多樣性和生態系統的結構與功能。

為從協同研究的角度探究太平洋西邊界流與中國近海生態系統演變之間的關系, 刻畫鄰近大洋對我國近海的影響, 中國科學院及時啟動了戰略性先導科技專項(A類)“熱帶西太平洋海洋系統物質能量交換及其影響”, 該專項的項目二為“黑潮及其變異對中國近海生態系統的影響研究”?；谠擁椖繉S、東海大量的現場觀測調查、歷史資料和數值模擬研究,目前對黑潮入侵東海陸架的途徑、機理和影響取得了重要研究進展, 得到了黑潮入侵東海陸架底部分支結構、變異及其物質輸送影響的重要科學結果。

1　黑潮入侵東海的主要研究進展

1.1　構建高精度的海洋環流模式

為了研究黑潮入侵東海陸架及其變化, 基于ROMS (Regional Ocean Modelling System)模式(Shchepetkinet al, 2005), 針對東海和西北太平洋相互影響區域, 構建了高分辨率的數值模式。模擬區域的經度范圍為 105°—136°E, 緯度范圍為 15°—41°N,垂向分為 31 層; 模式的水平分辨率為 3′×3′cos(θ),θ表示緯度。該數值模式的水平分辨率在東?？梢赃_到4km, 可以較好地分辨臺灣以東的中尺度渦旋。利用氣候態驅動場, 對該海域進行了7年的氣候態模擬以便模式 spin-up, 然后在氣候態模擬結果的基礎上,利用熱通量、動量通量、水通量以及10個天文分潮對海洋模式進行了 23年(1993—2016年)的后報模擬(Egbertet al, 2002; Deeet al, 2011), 模擬結果可以很好地再現西邊界流、黑潮入侵呂宋海峽、東海陸架觀測斷面的溫、鹽結構、海洋環流結構等(Guanet al,2006; Guoet al, 2006; Huet al, 2008, 2010; Maturiet al, 2014; Nanet al, 2015)。并且該模式可以很好地再現春季黑潮入侵東海陸架的分支結構等(圖 1), 為該模式結合觀測及生態動力學模式研究黑潮物質輸對東海生態的影響等過程奠定了基礎。

1.2　多學科證實黑潮入侵東海的近岸底部分支(NKBC)是影響近岸的關鍵途徑

為了解長江口鄰近海域底層海流情況和溫度、鹽度特征, 在近海 60m等深線處布放了一套坐底潛標(布放位置見圖 2a黃色三角形), 并成功進行了回收,獲得了2014.6—2014.10近4個月的觀測資料。在濾掉潮流以后, 可以發現該處的底層流速主要指向北和東北(圖2b)。對比模擬得到的2014年6月的月平均底層流可以發現, 該處位于黑潮入侵分支的北端,模擬的結果也表明該處的流場主要以北向和東北向流為主(圖2a)。該資料的觀測數據還表明, 長江口鄰近海域底層水的鹽度高于34.3, 溫度在18°C左右(圖2c), 該溫度、鹽度特征表明其源于黑潮次表層水。非常值得注意的是, 該處底層溫度存在很強的約 40天的振蕩, 該振蕩既可能是黑潮分支水和近岸水局地博弈的結果, 又可能反映了黑潮分支水上游的變化信息。黑潮入侵分支可能受到其他水動力學過程的影響, 如黑潮的高頻擺動、陸架波、黑潮和陸架水界面處的鋒面渦旋等都可能對黑潮入侵分支的高頻變化產生影響(Zhanget al, 2001; Huet al, 2008; Changet al, 2011)。此外, 觀測到的溫度變化范圍高達2°C, 這對該區域底層海洋生態和低氧過程具有潛在影響。

綜上所述, 該潛標資料證實了 Yang等( 2012,2013)基于夏季數值模擬結果和斷面觀測資料提出的黑潮入侵分支確實存在, 并且黑潮入侵分支可以在從6月到10月持續存在。由于地球行星β效應和地形效應, 黑潮主要沿著大陸坡流動, 在東海陸架的大部分海域, 黑潮的存在像一個屏障把陸架水和大洋水隔開, 使得外海大洋水難以直接進入東海陸架。該黑潮底部入侵分支是聯系外海大洋和近海的一個重要通道, 通過該通道外海大洋水可以進入到長江口鄰近海域的底層, 并通過上升流上升到海洋的表層(Lvet al, 2006; Qiaoet al, 2006), 從而對該生態災害高發區的溶解氧和赤潮產生重要的調控作用。

圖1　春季(5月)表層(a)和底層(b)流場Fig.1 The surface(a) and bottom(b) current fields in spring(May)

同時, 基于化學要素觀測, 構建了Ba-T-S黑潮入侵指示體系, 化學指標也證實了黑潮入侵底部分支的存在(圖 3)。航次斷面資料的夏季分析結果表明黑潮次表層水從臺灣東北陸架坡折處沿底部向北偏西方向入侵東海陸架, 至122°E左右(DH9-5站位)底部,黑潮水可占95%。隨后, 在 約27°N, 122°E附近該水體分為近岸分支和離岸分支, 繼續向東海入侵, 其中近岸分支可以一直入侵到浙江近岸(29.5°N), 而離岸分支則沿約100m等深線轉向東北方向。黑潮次表層水入侵流至錢塘江口附近時, 黑潮水占比仍可達到65%左右。垂直方向上, 陸架外側站位受黑潮次表層水的影響范圍更大, 黑潮水占 50%比例位置可延伸至外側 TW0-1站位表層, 而內側靠近大陸的站位則只限于陸架中部位置底層(Liuet al, 2017)。

另外, 由于溶解無機碘(IO3-和 I-)作為一種半保守化學指標在水團辨析上較溫度和鹽度更為靈敏(Wong, 1995; Truesdaleet al, 2002, 2003; Wonget al,2003), 因此能夠指示黑潮次表層水(KSSW)在臺灣島東北側的入侵(Wonget al, 2004)及在陸架上的遷移路徑?；谌芙鉄o機碘、鹽度、水團密度等參數, 構建了多參數水團辨析方法, 并用于指示黑潮次表層水入侵東海陸架海域的路徑。該化學觀測指標也證實了黑潮入侵分支的存在(Zhouet al, 2017)。如圖4所示,2014年5月, KSSW主要從臺灣島東北部入侵東海陸架, 在DH9斷面外海區域(26°N), KSSW的厚度接近80m。而隨著KSSW沿底部向北偏西方向入侵, 其厚度逐漸減小, 影響范圍逐漸變小。隨后, 在 DH6和DH7之間(約27.5°N), KSSW分為近岸分支和離岸分支, 其厚度均為20m左右, 并繼續北向遷移。其中近岸分支可以一直入侵到浙江近岸 50m等深線以內的次表層位置(約 29.5°N)。

總之, 物理、化學等多學科觀測及模擬研究成果均已證實黑潮入侵東海近岸底部分支存在, 并且是連接外海太平洋和東海陸架近岸水的一個關鍵通道,可以把外海大洋的物質能量輸入到東海近海, 使得浙江外海的赤潮、低氧高發區不僅受到人類活動的影響(例如人類活動引起長江入海物質變化等), 而且受到外海大洋變化的調控。

1.3　黑潮入侵東海近岸底部分支(NKBC)存在不同時空尺度的變化

圖2　模擬的月均底層流場(a)及潛標測得的流速、溫、鹽、溶解氧分布(b, c)Fig.2 The modeled monthly bottom current field (a) and the distribution of velocity, temperature, salinity, and dissolved oxygen recorded by the marine submersible buoy (b,c) (its location is shown by the yellow triangle in Fig.2a)

2014—2015年航次的溫鹽觀測數據(圖5)和流場數值模擬結果(圖 6)都表明黑潮入侵分支存在很強的季節變化。通過浙江外海一個斷面(DH-4)的高鹽水(鹽度＞34.3)通量的計算表明, 通過該斷面的高鹽水通量在3—4月份就開始在浙江外海出現并逐漸增強,在7月份水通量達到最強, 可入侵達到近岸區域, 然后逐漸減弱, 在10月份高鹽水逐漸消失, 說明秋、冬季節黑潮入侵分支向外海偏移, 不能到達浙江外海50m等深線以西的近岸區域。化學多參數水團辨析方法同樣證明了入侵存在季節差異: 春季入侵近岸分支可到達近岸區域, 而秋季入侵分支離岸偏向黑潮主軸(圖 7)。

圖3　基于Ba-T-S指示體系的黑潮次表層入侵水占比Fig.3 Proportions of the Kuroshio subsurface intrusion water based on the Ba-T-S index

圖4　2014年5月, 基于多參數水團辨析方法指示的黑潮入侵路徑Fig.4 The constructed Kuroshio intrusion path during May, 2014 determined by the method of multi-parameter water mass analysis注: 顏色代表入侵水團的厚度變化

對比不同年份的斷面觀測結果, 以鹽度為指示可以看出, 黑潮底部分支存在顯著的年際變化: 2009年顯示出到達近岸的入侵, 而 2012年入侵遠離近岸偏向外海(圖 8)。模擬的水通量結果也表明, 黑潮底部入侵分支在浙江近海出現的時間和量值存在年際變化(圖 9)。

圖5　實測的不同季節的溫度(左列, 單位: °C)和鹽度(右)剖面Fig.5 The observed temperature(left) and salinity(right) profiles of four seasons

圖6　模擬的不同季節的底層流場(從上到下: 1、4、7、10月; 流場為月平均的底層流場, 水深大于200m區域的底層流場用200m層的流場代替)Fig.6 The modeled bottom current fields of different seasons (a. January b. April, c. July, d. October)

圖8　不同年份的觀測站位(左)及鹽度剖面(右)Fig.8 The locations of stations(left) and the observed salinity profiles(right) of different years

圖9　斷面位置(a)及模擬的不同年份通過斷面的水通量(b)Fig.9 The location of the transect (left, the black line) and the water volumes that cross the transect in different years

1.4　提出了黑潮入侵東海分支存在和變化的動力機制

黑潮入侵分支存在明顯的季節變化和年際變化,黑潮入侵分支為什么會發生？了解其背后的物理機制, 才可以更好的預測。基于地球流體動力學理論,我們提出了地形β-spiral理論, 該理論是大洋(Stommelet al, 1977) β-spiral在陸架的延伸。該理論提出: 當底層流的方向和地形的梯度方向不垂直的時候(也就是說當入侵流的方向不平行于等深線的時候), 海流就會從表層到底層發生逐層偏轉。如果底層流的方向和地形梯度方向的夾角小于 90°, 入侵陸架的流就會逐漸發生從表層到底層的氣旋式偏轉,這對應于很強的底層入侵; 反之, 如果底層流的方向和地形梯度方向的夾角大于 90°, 入侵陸架的流就會逐漸發生從表層到底層的反氣旋式偏轉, 這時對應于強的表層水入侵?；谠摾碚? 可以預測: 在東海陸架上, 凡是黑潮路徑和等深線不平行的地方都會發生入侵流的垂向偏轉, 并且底層流的方向和地形梯度的角度決定了入侵主要發生在底層還是表層(圖 10)。

數值模擬結果的動量方程診斷結果表明, 黑潮主軸存在季節變化, 當黑潮主軸發生向岸或向海的偏轉時, 就會導致黑潮底層流方向和地形梯度的夾角發生變化, 從而引起黑潮入侵分支的季節變化。此外臺灣以東的渦旋, 會影響黑潮的水平密度結構, 從而影響黑潮底層流的方向(Yinet al, 2017), 根據所提出的地形 β-spiral理論, 其會進一步影響黑潮對東海陸架的入侵。

圖10　解釋黑潮入侵東海分支的地形β-spiral理論Fig.10 The topography β-spiral theory to explain the Kuroshio intrusion branch into East China Sea

1.5　揭示了黑潮物質輸入對東海陸架生態過程的影響

為了研究黑潮底部分支輸入的磷酸鹽的年際變化對長江口鄰近海域的影響, 在臺灣以東選定的一個固定斷面上(Yanget al, 2013), 每年的1月1日根據觀測資料給定磷酸鹽的濃度, 并在每年的隨后計算中, 保持濃度不變。在數值模擬中, 臺灣東北的黑潮輸入的磷酸鹽是東海磷酸鹽的唯一來源。在每年數值模擬的初始條件中, 整個海區的磷酸鹽濃度設為 0,并且在數值試驗中磷酸鹽的濃度被認為是保守的,即只考慮磷酸鹽的平流輸送和擴散過程。模擬結果表明, 2005年和2014年黑潮輸入的磷酸鹽在長江口鄰近海域表層的出現時間相差了約半個月, 并且其出現的規模也存在較大區別(圖11)。黑潮底部分支在此處出現的早晚意味著該區域底層水溫和磷酸鹽濃度的變化時間存在年際變化。對應于水文和生化環境的年際變化, 相應地赤潮出現的早晚和規模也會發生變化。

同時, 針對 2009年夏季個例情況, 建立了環流與生態耦合模式并進行了生態過程模擬, 模擬結果顯示, 包括黑潮底部入侵分支磷酸鹽輸送的模擬結果與實測葉綠素的分布結果更加一致(圖12), 說明底部分支供給的營養鹽有助于維持長江口外鄰近海域葉綠素生長, 進而影響生態過程。

圖11　不同年份黑潮入侵底部分支對長江口鄰近海域磷酸鹽的輸送量Fig.11 The phosphates that being transported from the Kuroshio bottom branch to Zhejiang coastal region in different years

圖12　模擬和觀測葉綠素結果比較Fig.12 Comparisons between the modeled and observed concentrations of chlorophylls

此外, 基于2014年5月份春季調查資料, 分析了調查期間黃東海原綠球藻分布狀況。值得注意的是在臺灣島東北部存在一個高的原綠球藻分布區, 而在長江口鄰近海域也存在原綠球藻的高濃度分布區,并且兩個高濃度區之間也存在著高濃度條帶分布(圖13)。這也從側面表明了黑潮入侵水的路徑和特征。

圖13　2014年5月春季調查期間黃東海原綠球藻分布狀況Fig.13 Distribution of prochlorococcus in the Yellow and East China Seas during May, 2014

該結果給出了另外一個重大的啟示: 黑潮入侵分支可以把某些熱帶藻種, 輸入到中國近海從而引起近海生態過程的變化。已有研究結果也表明, 某些近岸的藻類, 并不是我國的局地種, 而是外來藻種。此外, 2015年8月底棲生物的調查資料表明, 在黑潮入侵分支的流經區域, 其底棲生物也存在與周邊不一樣的類群(圖14)。

2　急需回答的科學問題和未來的研究展望

國內外已意識到近海生態系統的變化與鄰近大洋乃至更大范圍的環境變化密切相關, 需要從更為廣闊的空間和時間尺度上開展協同研究。例如西北太平洋沿岸上升流區初級生產力異常分布和低氧區異常主要由鄰近大洋的亞極地水異常入侵引起, 東北太平洋海洋環境的大時間尺度(年際、年代際)變異,對其海洋生態系統的年際變化起到主要的控制作用(Granthamet al, 2004)。Diaz等(2008)指出: 近海上升流區的生態系統, 除了受人類活動以及入海徑流的影響外, 其本身還受到外海大尺度海洋環境變化的影響。

圖14　2015年8月底棲生物分布特征Fig.14 Distribution of the benthos in August, 2015

目前的研究已經回答了大洋西邊流黑潮影響東海陸架的主要途徑及其變化、機制和初步影響效應,但尚有許多科學問題有待回答: 根據繞島理論, 日本以東太平洋大洋環流必然會導致黑潮入侵東海, 其調整機制和過程目前仍然不清楚。黑潮入侵東海陸架底部分支與太平洋氣候變化信號存在什么關系？大洋中尺度渦如何影響黑潮入侵東海陸架分支？黑潮入侵分支和臺灣暖流的相互作用是什么樣的？外海大洋輸入的物質, 有多少進入了我國近岸的生態系統中, 到底對生態過程和赤潮災害的發生、低氧區的形成等起怎樣的作用？哪些區域人類活動比外海入侵更重要, 哪些區域外海入侵起控制作用？近海的藻類有多少是外來物種, 該物種是如何進入到近海的? 這些物種對我國的生態系統的影響是什么樣的?東海原甲藻是否存在孢囊？該孢囊是否可以通過黑潮入侵分支輸入到浙江外海的赤潮高發區？等等。

這些都是未來需要重點回答的科學問題, 這些問題的回答, 將有助于深化外海大洋物質能量輸入對近海生態系統的認識, 提高對近海生態系統變動的預測水平。

致謝本文的數值模擬工作是在中國科學院海洋研究所的高性能計算平臺上完成的,特此致謝。

Chang Y L, Oey L Y, 2011. Interannual and seasonal variations of Kuroshio transport east of Taiwan inferred from 29 years of tide-gauge data. Geophys Res Lett, 38(8): L08603

Chavez F P, Ryan J, Lluch-Cota S Eet al, 2003. From anchovies to sardines and back: multidecadal change in the Pacific Ocean. Science, 299(5604): 217—221

Dee D P, Uppala S M, Simmons A Jet al, 2011. The ERA-interim reanalysis: configuration and performance of the data assimilation system. Q J Roy Meteor Soc, 137(656):553—597

Diaz R J, Rosenberg R, 2008. Spreading dead zones and consequences for marine ecosystems. Science, 321(5891):926—929

Egbert G D, Erofeeva S Y, 2002. Efficient inverse Modeling of barotropic ocean tides. J Atmos Ocean Technol, 19(2):183—204

Grantham B A, Chan F, Nielsen K Jet al, 2004. Upwellingdriven nearshore hypoxia signals ecosystem and oceanographic changes in the northeast Pacific. Nature,429(6993): 749—754

Guan B X, Fang G H, 2006. Winter counter-wind currents off the southeastern China coast: a review. J Oceanogr, 62(1):1—24

Guo X Y, Miyazawa Y, Yamagata T, 2006. The Kuroshio onshore intrusion along the shelf break of the East China Sea: the origin of the Tsushima Warm Current. J Phys Oceanogr,36(12): 2205—2231

Guo X Y, Zhu X H, Wu Q Set al, 2012. The Kuroshio nutrient stream and its temporal variation in the East China Sea. J Geophys Res: Oceans, 117(C1): C01026

Guo X Y, Zhu X H, Long Yet al, 2013. Spatial variations in the Kuroshio nutrient transport from the East China Sea to south of Japan. Biogeosciences, 10(10): 6403—6417

Hu J Y, Kawamura H, Li C Yet al, 2010. Review on current and seawater volume transport through the Taiwan strait. J Oceanogr, 66(5): 591—610

Hu X M, Xiong X J, Qiao F Let al, 2008. Surface current field and seasonal variability in the Kuroshio and adjacent regions derived from satellite-tracked drifter data. Acta Oceanol Sin,27(3): 11—29

Huang R X, 2010. Ocean Circulation: Wind-Driven and Thermohaline Processes. Cambridge: Cambridge University Press, 3302—3305

Isobe A, 2008. Recent advances in ocean-circulation research on the Yellow Sea and East China Sea shelves. J Oceanogr,64(4): 569—584

Liu W, Song J M, Yuan H Met al, 2017. Dissolved barium as a tracer of Kuroshio incursion in the Kuroshio region east of Taiwan Island and the adjacent East China Sea. Sci China Earth Sci, 60(7): 1356—1367, http://dx.doi.org/10.1007/s 11430-016-9039-7

Lv X G, Qiao F L, Xia C Set al, 2006. Upwelling off Yangtze River estuary in summer. J Geophys Res: Oceans, 111(C11):C11S08

Maturi E, Sapper J, Harris Aet al, 2014. GHRSST Level 4 OSPO Global Foundation Sea Surface Temperature Analysis (GDS version 2). National Oceanographic Data Center, NOAA,Dataset

Nan F, Xue H J, Yu F, 2015. Kuroshio intrusion into the South China Sea: a review. Prog Oceanogr, 137: 314—333

Pedlosky J, 1996. Ocean Circulation Theory. Berlin, Heidelberg:Springer-Verlag

Qiao F L, Yang Y Z, Lv X Get al, 2006. Coastal upwelling in the East China Sea in winter. J Geophys Res: Oceans, 111(C11):C11S06

Shchepetkin A F, McWilliams J C, 2005. The regional oceanic modeling system (ROMS): a split-explicit, free-surface,topography-following-coordinate oceanic model. Ocean Model, 9(4): 347—404

Stommel H, Schott F, 1977. The beta spiral and the determination of the absolute velocity field from hydrographic station data.Deep Sea Res, 24(3): 325—329

Su J L, 1998. Circulation dynamics of the China seas north of 18°N. In: Robinson A R, Brink K H eds. The Sea, The Global Coastal Ocean, Regional Studies and Syntheses. New York: John Wiley, 483—506

Truesdale V W, Bailey G W, 2002. Iodine distribution in the Southern Benguela system during an upwelling episode.Cont Shelf Res, 22(1): 39—49

Truesdale V W, Upstill-Goddard R, 2003. Dissolved iodate and total iodine along the British east coast. Estuar Coast Shelf Sci, 56(2): 261—270

Wong G T F, 1995. Dissolved iodine across the Gulf Stream front and in the South Atlantic Bight. Deep Sea Res Part I Oceanogr Res Pap, 42(11-12): 2005—2023

Wong G T F, Hung C C, Gong G C, 2004. Dissolved iodine species in the East China Sea—a complementary tracer for upwelling water on the shelf. Cont Shelf Res, 24(13—14):1465—1484

Wong G T F, Zhang L S, 2003. Geochemical dynamics of iodine in marginal seas: the southern East China Sea. Deep Sea Res Part II Top Stud Oceanogr, 50(6—7): 1147—1162

Yang D Z, Yin B S, Liu Z Let al, 2011. Numerical study of the ocean circulation on the East China Sea shelf and a Kuroshio bottom branch northeast of Taiwan in summer. J Geophys Res: Oceans, 116(C5): C05015

Yang D Z, Yin B S, Liu Z Let al, 2012. Numerical study on the pattern and origins of Kuroshio branches in the bottom water of southern East China Sea in summer. J Geophys Res:Oceans, 117(C2): C02014

Yang D Z, Yin B S, Sun J Cet al, 2013. Numerical study on the origins and the forcing mechanism of the phosphate in upwelling areas off the coast of Zhejiang province, China in summer. J Mar Syst, 123-124: 1—18

Yin Y Q, Lin X P, He R Yet al, 2017. Impact of mesoscale eddies on Kuroshio intrusion variability northeast of Taiwan. J Geophys Res: Oceans, 122(4): 3021—3040

Zhang D X, Lee T N, Johns W Eet al, 2001. The Kuroshio east of Taiwan: modes of variability and relationship to interior ocean mesoscale eddies. J Phys Oceanogr, 31(4): 1054—1074

Zhou P, Song X X, Yuan Y Qet al, 2017. Intrusion pattern of the Kuroshio Subsurface Water onto the East China Sea continental shelf traced by dissolved inorganic iodine species during the spring and autumn of 2014. Mar Chem,http://dx.doi.org/10.1016/j.marchem.2017.07.006