北黃海冬季溫鹽年際變化時空模態與氣候響應

石強

（1.山東省海洋生態環境與防災減災重點實驗室，山東 青島，266033；2.國家海洋局北海環境監測中心，山東 青島 266033；3.海洋溢油鑒別與損害評估技術國家海洋局重點實驗室，山東 青島 266033）

黃海是一個被中國大陸和朝鮮半島圍繞的半封閉陸架淺海，黃海南側以濟州島至長江口連線為界與東海相鄰；北側以旅順口至蓬萊連線為界與渤海相鄰；以山東半島東端成山頭為界分為南、北黃海。黃海表面積約40×104km2，平均水深44 m，最大水深140 m。在黃海中部分布著南北方向伸展的黃海槽，水深約50～80 m。黃海主要分布有5種水團：黃海暖流水團、黃海底層冷水團、朝鮮沿岸水團、中國沿岸水團和長江沖淡水團。黃海環流分為兩種季節性流系：冬季環流和夏季環流。

在山東半島東端的成山頭與大連之間有一個自1928年開始監測的海洋調查斷面，張以懇等（1996）研究了該斷面自1928-1990年期間選取的42年觀測資料，分析了該斷面夏季8月份冷水團年際變化特征。趙保仁等（1996）、江蓓潔等（2007）研究了北黃海冷水團的溫度、鹽度以及海流系統的季節變化。王宗山等（1992）、徐玲玲等（2008）研究了北黃海溫鹽垂直分布。黃江嬋等（2012）根據2006-2007年期間對北黃海19條斷面調查資料，研究了北黃海四季水團分布特征。

北黃海冬季溫鹽結構受到北黃海環流和西北季風的共同影響，冬季渤南沿岸水勢力最強，該水團在西北季風作用下可以越過山東半島東端成山頭向南黃海伸展（李繁華等，1989）。北黃海魯北沿岸流冬季最強，春季減弱，夏秋季消失（鮑獻文等，2009）。在秋冬季節山東半島東及北沿岸，海水交換又呈上下不同的兩層模式：上層黃海沿岸流自西向東南下，而在深底層，黃海暖流余脈的深底層流卻自南而北再向西流動，致使表層溫度低于下層溫度，形成冬秋季逆溫層現象（顏文彬等，1991；管秉賢，2000）。

隨著全球氣候變暖，黃海海域冬季氣候也有所改變。冬季2月份平均氣溫增暖最快（張黎紅等，2005；陳少勇等，2009），東亞季風明顯變弱（黃榮輝等，2008）。冬季氣候要素的變化對北黃海冬季溫鹽結構也產生顯著影響，而對北黃海冬季溫鹽結構的氣候響應仍然缺少研究，特別是對北黃海斷面冬季溫度、鹽度連續長期觀測資料分析有重要的科學意義。

采用旋轉經驗正交函數（REOF）等方法研究近36年來北黃海冬季斷面溫度、鹽度年際變化時空結構與氣候響應，以及冬季北黃海與渤海之間溫鹽年際變化的相關性。這些工作有助了解北黃海冬季溫鹽年際時空變化趨勢和演變過程，并進一步研究北黃海冬季溫鹽年際變化與生化要素的影響關系。

1 材料與方法

1.1 分析資料

北黃海斷面8個監測站分布在北黃海成山頭至大連之間，簡稱CD斷面。該斷面南端水深約31m；中段水深約52 m；北端水深約48 m，見圖1。國家海洋局北海分局海洋調查隊和北海環境監測中心于1977-2012年歷年1-2月（多數是2月份）對斷面表底層海水溫度、鹽度監測，分析方法按照海洋調查規范一篇2部分：海洋水文觀測（GB/T 12763.2-2007）執行，其中1984年、1996年各站缺測，分析中缺測數據由時間相鄰站點資料內插填補。

圖1 北黃海斷面站位分布

位于北黃海近岸的大連站自1977-2012年歷年1月、2月月平均氣溫、每日2次（0、12世界時）97 m高度處探空風速資料，由中國氣象局氣象信息中心提供。

1.2 分析方法

1.2.1 旋轉經驗正交函數（REOF）方法

旋轉經驗正交函數（REOF）方法是經驗正交函數（EOF）方法的改進。EOF方法是一種分析矩陣數據中的結構特征和提取數據主要時空特征量的方法，在海洋學科中得到了廣泛的應用（Storch，1999；石強 等，2001，2011）。

EOF方法可以反映分析域整體時空變化，當需要關注分析域局地時空結構時，采用旋轉經驗正交函數（REOF）是一個適宜的方法（Storch，1999）。REOF方法是在保持EOF前K個主成分的總方差不變情況下，進行最大方差正交變換，使變換后的空間型對原變量的方差貢獻差異增大，以反映原變量中的局部結構，第K個特征值方差貢獻用REOFk（%）表示。REOF詳細討論見（吳洪寶等，2010）。

1.2.2 躍變分析方法

采用信噪比公式確定時間序列資料躍變發生的強度與時刻，信噪比公式為：

其中，M1、M2和σ1、σ2分別是基準時刻之前與之后期間（n）內的平均值和標準差，當信噪比值出現局部最大值且大于1.0時，判斷在基準時刻出現躍變現象（Yamamoto et al，1986）。躍變現象可以出現在數年至數十年的多種時間尺度中，本研究主要關注十年尺度的溫鹽模態時間分量躍變，所以取躍變時間尺度n=10年。

1.2.3 正態分布檢驗

當樣本序列符合正態分布時，可以用樣本平均值代表樣本序列的統計特征（陶澍，1994），因此，本文除特別聲明外，計算出的各年斷面溫鹽平均值的觀測序列均通過了正態分布峰度、偏度系數判據檢驗，信度0.001。

2 結果與討論

2.1 氣候要素與表底層溫度、鹽度年際變化

據1961-2006年月平均溫度資料分析，我國東部季風區冬季增溫率為0.39℃/10 a，其中2月增溫最大（張少勇等，2009）。對大連站1977-2012年歷年2月月平均氣溫和歷年1月月緯向、經向偽風應力（Pseudostress）和（（τxτ）y=∑（WxW）y進行5項高斯型低通濾波處理。2月平均氣溫有顯著線性升高趨勢（圖2a）；1月月經向擬風應力和有顯著線性減弱趨勢（圖2b）；月緯向偽風應力和有非顯著線性減弱趨勢（圖2c）。冬季氣溫升高和西北季風強度減弱的氣候趨勢將對北黃海溫鹽結構年際變化產生影響。

將溫度、鹽度數據序列分別沿各站位時間方向做邊界一階、內部二階Shapio濾波，以減少短波噪聲（Shapiro，1970），并繪制溫度、鹽度分布圖 3-圖 6。

圖2 大連站2月平均氣溫（a）和1月經向（b）、緯向（c）偽風應力和

圖3 斷面表層溫度分布（單位：℃）

由圖3、4所示，斷面表底層溫度經歷了10年冷期（1977-1987年）、20年暖期（1988-2008年）和近4年冷期（2009-2012年）的變化，其中2007-2008年是最暖的時期，斷面中部的最高溫度達到5℃以上，次暖期在1995-1999年，暖期海水溫度升高主要在斷面中部海域，并向兩側近岸擴展；冷期最低溫度在1℃以下，出現在1977-1981年斷面的南端海域，次冷期在2010-2012年，冷期海水低溫主要在斷面的南端近岸海域，并向斷面中部海域擴展。

圖4 斷面底層溫度分布（單位：℃）

圖5 斷面表層鹽度分布

圖6 斷面底層鹽度分布

由圖5、6所示，斷面鹽度分布表底層有差異，表層鹽度南、北端有低鹽度水體出現，底層鹽度相對表層分布比較均勻。高鹽度水體主要出現在斷面中部海域，低鹽度水體出現在斷面南、北端近岸海域。表底層高鹽度期在1981-1984年、1989-1990年、1999-2003年和2008-2010年期間；低鹽度期在1977-1980年、1985-1988年、1991-1993年和2011-2012年。

斷面表層平均溫度年際變化呈現顯著線性升高趨勢（圖7a），斷面表層平均鹽度年際變化呈現準平衡狀態（圖7b），底層平均溫鹽值與表層平均溫鹽是相似的年際變化。

圖7 斷面表層平均溫度（a）、平均鹽度（b）

北黃海冬季西北季風將表層水體向南輸送，在魯北沿岸堆積的南移水體將攜帶低溫、高鹽效應的表層水向深層移動；在遼南沿岸表層水體被吹離，將引起底層相對高溫、高鹽水體向遼南沿岸補償流動，由此，在北黃海形成一種西北季風驅動的風動-補償環流。由于氣候變暖，北黃海冬季西北季風強度逐漸減弱（圖2b、2c），北黃海冬季風動-補償環流也逐漸減弱，從斷面表底層溫鹽差年際分布分析，表底層溫鹽差均存在顯著的線性趨勢（圖8a、8b），

表底層溫度顯著升高，表層平均溫度升高幅度比底層大，表底層平均溫差減小；表底層平均鹽度有升高趨勢（未通過顯著性判據），表層平均鹽度升高幅度比底層大，表底層平均鹽度差顯著減小。因此，斷面冬季表底層平均溫鹽差顯著線性正距平趨勢，底層溫度、鹽度出現高于表層的情形將減少。

圖8 斷面平均表底層溫度差（a）、鹽度差 (b)分布

2.2 表底層溫度年際時空變化

采用REOF分析斷面表底層溫度距平序列資料，取EOF前2項做旋轉變換，時間分量采用5項高斯型低通濾波處理，消除短期擾動。

根據斷面表底層溫度模態空間分量高值分布，斷面表底層溫度模態空間分量有兩種型式：開闊海型（圖 9d、圖 10c）和南北近岸型（圖 9c、圖10d）。表層溫度模態中，兩種模態方差貢獻差別不大（相差2%）；底層溫度模態中，開闊海型方差貢獻大于南北近岸型（相差20%）。

圖9 斷面表層溫度時空模態（a，c.REOF1=47%，b，d.REOF2=45%），↓躍變時刻

圖10 斷面底層溫度時空模態

斷面表層溫度模態第二時間分量線性擬合相關R=0.23，與信度0.05相關系數臨界值0.28十分接近，存在近似的線性升高趨勢（圖9b），斷面其他表底層溫度模態時間分量均有顯著線性升高趨勢（圖9a、圖10a、10b）。開闊海型溫度模態中底層模態線性升溫趨勢比表層快，南北近岸型溫度模態中表層模態線性升溫趨勢比底層快。最大熵譜分析結果：表層第二模態和底層第一模態時間分量顯著周期5.5年、9.0年和36.0年；表層第一時間分量顯著周期5.1年、6.5年和14.4年；底層第二時間分量顯著周期5.5年、10.3年。表底層開闊海型溫度模態有相同的顯著年際周期變化，表底層南北近岸型溫度模態顯著年際周期不同。

表底層溫度南北近岸型模態時間分量存在十年尺度躍變現象；開闊海型模態無躍變現象（圖9a，圖 10b）。

2.3 斷面溫度的氣候響應

延遲分析表明：斷面2月溫度模態與大連站歷年1月經向偽風應力和序列顯著相關；與2月經向偽風應力和序列無顯著相關，溫度模態位相落后1個月對風應力響應的情況與其他學者分析結果是一致的（劉哲等，2003）。2月溫度模態與2月平均氣溫序列有顯著同步正相關，與1月平均氣溫模態顯著正相關分布有峰值偏離現象（圖略）。

圖11 大連站2月平均氣溫與斷面表底層溫度模態延遲相關系數

圖11 所示，斷面表層溫度第一、第二模態和底層溫度第一模態對2月平均氣溫有顯著同步正相關響應，底層溫度第二模態與月平均氣溫顯著正相關響應中存在其他因素的干擾較大，顯著相關峰值有延遲偏離現象。

由圖12所示，斷面表層溫度第一模態和底層溫度第二模態與1月經向偽風應力有顯著正相關響應，其他溫度模態未有顯著相關響應。溫度模態與2月經向偽風應力和未有顯著正相關峰值響應（圖略）。斷面表底層溫度模態對氣溫的顯著響應比對偽風應力的顯著響應快1個月時間，說明西北季風驅動的風動環流對溫度變化位相影響時間大約是1個月時間。表層溫度第一模態和底層溫度第二模態存在對偽風應力顯著延遲1個月響應，可能導致對2月平均氣溫同步相關系數分布中出現峰值弱化和偏移的現象（圖11）。

斷面2月表底層溫度開闊海型模態是對2月平均氣溫年際變化的響應，表層溫度開闊海型模態時間分量線性升高趨勢斜率大于底層；2月表底層溫度南北近岸型模態是對1月經向風應力和與2月平均氣溫的響應，底層溫度南北近岸型模態時間分量線性升高趨勢斜率大于表層。

圖12 大連站1月經向偽風應力和與斷面表底層溫度模態延遲相關系數

渤海與北黃海之間通過渤海海峽產生水體交換，北黃海海域冬季存在黃海水團、渤南（魯北）沿岸水、遼南沿岸水和渤-黃海混合水團（李繁華等，1989；黃江嬋等，2012）。斷面橫跨黃海水團（斷面中部）、渤南（魯北）沿岸水（斷面南部）和渤-黃海混合水團（斷面北部）。將同步的北黃海斷面溫度模態與渤海斷面溫度模態（石強，2013）進行延遲相關分析，研究冬季北黃海與渤海之間溫度模態的相關關系。

圖13 斷面表底層溫度模態與渤海斷面表底層溫度模態延遲相關系數

圖13 所示，冬季北黃海斷面表底層溫度開闊海型、南北近岸型模態與渤海斷面表底層溫度開闊海型模態有顯著同步正相關，斷面底層溫度開闊海型模態與渤海底層溫度遼河口型模態有顯著同步正相關。由此分析，渤海斷面開闊海域表底層溫度主要與北黃海斷面表底層溫度同步年際變化，北黃海斷面底層溫度開闊海模態與渤海斷面底層溫度遼河口型模態有比較顯著的同步正相關，而渤海斷面表底層溫度的黃河口型模態與表層溫度遼河口型模態是顯著局地效應產生的年際變化，與北黃海斷面溫度年際變化無關。

2.4 表底層鹽度年際時空變化

采用REOF分析斷面表底層鹽度距平序列資料，取EOF前2項做旋轉，時間分量采用5項高斯型低通濾波處理，消除短期擾動。

根據斷面表底層鹽度模態空間分量高值分布，斷面表底層鹽度模態空間分量有二種形態：魯北近海型（圖14c、圖15c）和遼南近海型（圖14d、圖15d）。斷面鹽度年際變化主要出現在斷面南北部海域，這種形態說明渤南（魯北）沿岸水和遼南沿岸水鹽度的年際變化是影響北黃海斷面海域鹽度年際變化的主要因素。

表底層鹽度模態時間分量出現兩種相反的線性趨勢變化，表底層鹽度第一模態（圖14a、圖15a）有顯著線性升高趨勢；表底層鹽度第二模態（圖14b、圖15b）有顯著線性降低趨勢。這種形態說明渤南（魯北）沿岸水影響的鹽度模態有顯著線性升高趨勢，遼南沿岸水影響的鹽度模態有顯著線性降低趨勢，并且升高斜率大于降低斜率。最大熵譜分析：表底層鹽度第一模態顯著周期9.0年；表層鹽度第二模態顯著周期5.1年、6.5年；底層鹽度第二模態顯著周期6.5年。

圖14 斷面表層鹽度時空模態

圖15 斷面底層鹽度時空模態

表底層鹽度魯北近岸型模態時間分量存在十年尺度躍變現象；遼南近岸型模態無躍變現象（圖14a、圖 15a）。

表底層鹽度第一模態時間分量之間有顯著同步正相關系數0.94；第二模態之間有顯著同步正相關系數0.98；依據REOF分析時間分量的正交性，表、底層鹽度第一、二模態時間分量之間無顯著同步正相關系數。

2.5 鹽度模態的氣候響應

冬季遼南沿岸受鴨綠江等河流入海淡水影響，溫度和鹽度較低。遼南沿岸水4個季節都是以低鹽為特征（鮑獻文等，2009）。這種遼南沿岸水低鹽特征是因為遼南沿岸存在一個四季基本不變沿遼南岸流向西南的遼南沿岸流，這種流系夏季流速大，流幅窄；冬季流速小，流幅寬（臧家業等，2001）。北黃海冬季斷面海水鹽度年際變化受到遼南沿岸水、遼南沿岸流以及渤海-北黃海海水交換的影響。。

在冬季西北季風驅動下，渤南水經魯北沿岸流向南黃海，黃海口徑流量逐年減少的趨勢造成渤海冬季斷面鹽度黃河口型模態出現顯著線性升高效應（石強，2013），并且與北黃海斷面表底層鹽度魯北近海型模態顯著同步正相關（圖16）。由于渤海斷面表底層鹽度遼河口型模態滯后1～2年與北黃海斷面表底層鹽度與魯北近海型模態顯著正相關（圖16），因此，渤海斷面表底層鹽度黃河口型模態是影響北黃海斷面表底層鹽度魯北型模態年際變化的主要因素。。表底層鹽度遼南近海型模態與渤海冬季斷面鹽度遼東灣型、黃河口型模態無同步和延遲顯著相關，該模態主要是受遼南沿岸水鹽度的年際變化影響。

圖16 渤海鹽度模態與斷面鹽度模態延遲相關系數

延遲相關分析表明：斷面表底層鹽度魯北近岸模態與2月平均氣溫顯著同步正相關；表底層鹽度遼南近海型模態超前2～3年與2月平均氣溫顯著負相關。斷面表底層鹽度模態與1、2月經向、緯向偽風應力和不存在可以解釋的延遲相關關系（圖略）。由此分析，斷面表底層鹽度魯北近岸模態與2月平均氣溫同步年際變化，表底層鹽度遼南近海型模態形成機制與魯北近岸模態不同，表底層鹽度2個模態年際變化與西北季風強度年際變化無關，主要受到渤南（魯北）、遼南沿岸水鹽度年際變化影響。

為分析斷面溫鹽模態之間年際變化的穩定性，我們分析了表層、底層溫鹽兩個模態之間的延遲相關，結果表明：表（底）層溫度第一模態超前第二模態7～8年顯著正相關；表（底）層鹽度第一模態超前第二模態4年顯著負相關。斷面表層、底層溫（鹽）兩種模態之間年際變化存在顯著延遲相關性。

3 結論

（1）冬季北黃海斷面表底層溫度有兩種時空模態：開闊海型和南北近岸型。兩種模態時間分量主要有準周期和線性趨勢年際變化，開闊海型溫度模態中底層線性升溫斜率比表層大；南北近岸型溫度模態中表層線性升溫斜率比底層大。表底層溫度開闊海型模態是對2月平均氣溫年際變化的響應；表底層溫度南北近岸型模態是對1月經向風應力和與2月平均氣溫的響應。

（2）冬季北黃海斷面表底層鹽度有兩種時空模態：魯北近海型和遼南近海型。兩種模態時間分量主要有準周期和線性趨勢年際變化。表底層鹽度魯北近海型模態時間分量有顯著線性升高趨勢；遼南近海型模態時間分量有顯著線性降低趨勢，其中魯北近海型鹽度模態線性升高斜率大于遼南近海型模態降低斜率。斷面表底層鹽度年際變化與渤南（魯北）沿岸水和遼南沿岸水的鹽度年際變化顯著相關。

（3）表底層溫度南北近岸型模態時間分量存在10年尺度躍變現象；開闊海型溫度模態無躍變現象。表底層鹽度魯北近岸型模態時間分量存在10年尺度躍變現象；遼南近岸型模態無躍變現象。

Shapiro R，1970.Smoothing，filtering，and boundary effects.Reviews of Geophysics and Space Physics，8（2）：359-387.

Storch H V，Zwiers F W，1999.Statistical Analysis in Climate Research.Cambridge University Press.

Yamamoto R，Iwashima T，Sanga N K，1986.Ananalysis of climatic jump J Met Soc Japan，64：273-281.

鮑獻文，李娜，姚志剛，等，2009.北黃海溫鹽分布季節變化特征分析.中國海洋大學學報，39（4）：553-562.

陳少勇，郭忠祥，高蓉，等，2009.我國東部季風區冬季氣溫的氣候變暖特征.應用氣象學報，20（4）：476-485.

管秉賢，2000.山東半島北及東沿岸海域冬季的逆溫躍層及其與深底層逆流的關系.黃渤海海洋，18（3）：1-6.

黃江嬋，臧璐，石曉勇，等，2012.北黃海四季水團劃分及化學特征比較.漁業科學進展，33（1）：17-27.

黃榮輝，顧雷，陳際龍，等.2008.東亞季風系統的時空變化及其對我國氣候異常影響的最近研究進展.大氣科學，4：691-719.

江蓓潔，鮑獻文，吳德星，等，2007.北黃海冷水團溫、鹽多年變化特征及影響因素.海洋學報，29（4）：1-10

李繁華，劉愛菊，趙松鶴，等，1989.山東近海水文狀況.山東省地圖出版社.

劉哲，魏皓，蔣松年，2003.渤海多年月平均溫鹽場的季節變化特征及形成機制的初步分析.青島海洋大學學報，33（1）：7-14.

石強，2013.渤海冬季溫鹽年際變化時空模態與氣候響應.海洋通報，32（5）：505-513.

石強，陳江麟，李崇德，2001.渤海硝酸鹽氮和亞硝酸鹽氮季節循環分析.海洋通報，21（2）：22-29.

石強，楊東方，2011.渤海夏季海水pH值年際時空變化.中國環境科學，31增刊，58-68.

陶澍，1994.應用數理統計方法，中國環境科學出版社，北京.

王宗山，徐伯昌，鄒娥梅，等，1992.黃渤海底層溫、鹽度與其垂直平均值的統計關系.黃渤海海洋，10（1）：1-7.

吳洪寶，吳蕾，2010.氣候變率診斷和預測方法（第2版）.氣象出版社.

徐玲玲，林霄沛，吳德星.2008.風浪及地形對渤、黃海夏季溫度垂直結構的影響.中國海洋大學學報，38（2）：183-188.

顏文彬，陳玉英，1991.渤海和北黃海的逆溫躍層.海洋通報，10（4）：94-98.

臧家業，湯毓祥，鄒娥梅，等，2001.黃海環流分析，科學通報，46增刊：7-15.

張黎紅，王謙謙，2005.大連地區近44a冬季氣溫的變化特征，南京氣象學院學報，28（6）：827-832.

張以懇，楊玉玲.1996.夏季北黃海冷水團多年變化特征分析.海洋預報，13（4）：15-21.

趙保仁，1996.北黃海冷水團環流結構探討-潮混合鋒對環流結構的影響.海洋與湖沼，27（4）：429-435.