黃、渤海斷面及海洋站的鹽度分布特征與變化趨勢(shì)

馬 超, 鞠 霞, 吳德星, 林霄沛

(1. 中國(guó)海洋大學(xué) 海洋環(huán)境學(xué)院, 山東 青島 266100; 2. 國(guó)家海洋局 第一海洋研究所, 山東 青島 266061)

黃、渤海斷面及海洋站的鹽度分布特征與變化趨勢(shì)

馬 超1, 鞠 霞2, 吳德星1, 林霄沛1

(1. 中國(guó)海洋大學(xué) 海洋環(huán)境學(xué)院, 山東 青島 266100; 2. 國(guó)家海洋局 第一海洋研究所, 山東 青島 266061)

采用經(jīng)驗(yàn)正交函數(shù)(EOF)分解、超前滯后相關(guān)分析等方法, 分析了 36°N斷面、大連—成山頭斷面和渤海中部斷面鹽度的空間分布特征和時(shí)間變化特征, 討論了環(huán)黃、渤海海洋站之間鹽度的關(guān)系。結(jié)果表明, 36°N斷面、大連—成山頭斷面和渤海中部斷面鹽度的整體趨勢(shì)都是升高, 黃河流量的減小是造成渤海鹽度升高的重要原因。黃海千里巖站 1960～2000年間的鹽度整體趨勢(shì)也是升高, 年變化率為0.01/a, 造成鹽度升高的主要原因是降水的減少, 并且其短期震蕩可能和El Ni?o有關(guān), 長(zhǎng)期變化則與太平洋年代際震蕩(PDO)有關(guān)。此外, 千里巖站和其他海洋站之間的鹽度有較好的正相關(guān)關(guān)系, 只是在時(shí)間上存在著 1～3個(gè)月的超前。對(duì)黃、渤海鹽度長(zhǎng)期變化的系統(tǒng)分析表明, 近海鹽度可能受到外洋的影響, 有助于加深對(duì)中國(guó)近海環(huán)境演變的認(rèn)識(shí)。

黃海; 渤海; 鹽度; 變化趨勢(shì)

黃海是全部位于大陸架上的一個(gè)半封閉的淺海,地勢(shì)較為平坦, 夏季底層水溫低于8 ℃, 是有名的黃海冷水團(tuán)盤(pán)踞的海域。渤海是中國(guó)的內(nèi)海, 是深入中國(guó)大陸的近封閉型的一個(gè)淺海, 平均水深 18 m。兩者僅通過(guò)渤海海峽進(jìn)行交換, 表層有暖水舌從南黃海經(jīng)北黃海直指渤海海峽[1]。

作為中國(guó)重要的邊緣海, 黃、渤海海洋要素被廣泛研究, 發(fā)現(xiàn)其存在明顯的季節(jié)變化。劉哲等[2]認(rèn)為外強(qiáng)迫、岸形和地形所共同誘導(dǎo)的海水平均環(huán)流, 是形成渤海水平分布和季節(jié)變化的決定因素之一。

近年來(lái), 黃、渤海海洋要素的年際以上尺度的中長(zhǎng)期變化成為研究的熱門(mén)問(wèn)題。Lin等[3]認(rèn)為渤海海表鹽度, 氣溫和海表溫度在 1960～1997年間都是升高的、變化率分別是 0.074 /a, 0.024 ℃/a和0.011 ℃/a。方國(guó)洪等[4]認(rèn)為渤海在1965～1997年間,海表溫度升高 0.48 ℃, 海表鹽度升高 1.34, 海溫與氣溫及鹽度與降水都有較好的相關(guān)關(guān)系。白虹等[5]對(duì)黃海冷水團(tuán)的長(zhǎng)期變化作了研究, 認(rèn)為溫躍層頂?shù)奈恢糜心觌H變化, 且與El Nin?有關(guān)。這些研究大多是針對(duì)溫度的變化, 而鹽度對(duì)海洋內(nèi)部動(dòng)力結(jié)構(gòu)有重要影響。吳德星等[6]認(rèn)為1961～1996年渤海鹽度升高近2, 并認(rèn)為黃河入海流量持續(xù)銳減是導(dǎo)致渤海鹽度升高的主導(dǎo)原因。吳德星等[7]對(duì)渤海1958年和2000年夏季溫鹽場(chǎng)和流場(chǎng)結(jié)構(gòu)作了比較, 認(rèn)為鹽度場(chǎng)發(fā)生了根本性變異, 渤海灣及灣口區(qū)成為高鹽度區(qū), 從而導(dǎo)致渤海流場(chǎng)結(jié)構(gòu)發(fā)生變化。

為了研究中國(guó)近海鹽度的變化趨勢(shì), 探討外洋對(duì)近海鹽度的影響, 加深對(duì)中國(guó)近海環(huán)境演變的認(rèn)識(shí), 作者利用最新收集到的環(huán)黃、渤海海洋站和斷面調(diào)查資料, 針對(duì)黃、渤海鹽度的長(zhǎng)期變化進(jìn)行系統(tǒng)的研究, 并對(duì)其成因作簡(jiǎn)要分析。

1 資料來(lái)源和處理

斷面資料取自1976～2003年36°N大連—成山頭斷面, 1978～1999年渤海中部斷面的溫鹽深儀(CTD)調(diào)查資料。海洋站資料為鲅魚(yú)圈、葫蘆島、秦皇島、塘沽、龍口、煙臺(tái)、老虎灘7個(gè)沿岸站。1965～1995年的溫、鹽、降水資料和千里巖站 1960～2000年的鹽度資料, 站位見(jiàn)圖1。此外, 還有青島站1960～2000年的蒸發(fā)降水資料和利津站 1950～2000年的黃河入海徑流資料。

作者繪制了各個(gè)斷面的鹽度等值線(xiàn)圖, 分析了斷面鹽度的年際變化。對(duì)千里巖站鹽度的長(zhǎng)期變化趨勢(shì)進(jìn)行了研究, 并與其他環(huán)黃、渤海海洋站鹽度之間作了超前滯后分析。

圖1 北黃海與渤海地形及觀測(cè)站位Fig. 1 Topography of the northern Yellow Sea, the Bohai Sea, and the observation stations

對(duì)于黃、渤海海溫, 月最低值一般出現(xiàn)在 2月,月最高值一般出現(xiàn)在8月, 因此作者選用2月和8月分別代表冬、夏兩季。

對(duì)黃、渤海鹽度的長(zhǎng)期變化進(jìn)行的系統(tǒng)分析, 可以說(shuō)明近海鹽度變化的可能原因, 有助于加深對(duì)中國(guó)近海環(huán)境演變的認(rèn)識(shí)。

2 結(jié)果和討論

2.1 各斷面的鹽度分布特征

由圖2可知, 36°N斷面夏季上10 m垂直混合比較均勻, 鹽度東高西低, 等鹽線(xiàn)在陸坡處向下傾斜,在斷面東側(cè)20 m以下形成高鹽區(qū)。冬季從表層到底層垂直混合均勻, 等鹽度線(xiàn)豎直分布, 鹽度東高西 低。

圖2 鹽度斷面圖Fig. 2 Distribution of salinity at the sectionsa. 夏季36°N斷面; b. 夏季大連—成山頭斷面; c. 夏季渤海中部斷面; d. 冬季36°N斷面; e. 冬季大連—成山頭斷面; f. 冬季渤海中部斷面a. 36 °N section in summer; b. Dalian to Chengshantou section in summer; c. central Bohai section in summer; d. 36°N section in winter; e. Dalian to Chengshantou section in winter; f. central Bohai section in winter

大連—成山頭斷面夏季上10 m垂直混合比較均勻, 鹽度大致南高北低, 海峽北部10～15 m鹽度垂直梯度較大, 等鹽度線(xiàn)較密集, 并在中部海底突起處向下延伸, 在海峽北部15 m以下形成一個(gè)高鹽中心,位置與冷水中心一致。冬季垂直混合較強(qiáng), 等鹽度線(xiàn)豎直分布, 鹽度大致北高南低, 一定程度上反映出北黃海水從海峽北部進(jìn)入渤海, 萊州灣較低鹽度水沿海峽南側(cè)流出的特性, 最高值出現(xiàn)在中部偏北的狹長(zhǎng)區(qū)域。

渤海中部斷面在南部上 5 m最低, 原因是靠近黃河口, 受河流淡水輸入的影響; 渤海中部淺灘兩側(cè)的凹槽處是鹽度最高的區(qū)域, 溫度則對(duì)應(yīng)著兩個(gè)冷水中心; 渤海中部淺灘鹽度上下相對(duì)均勻, 等鹽度線(xiàn)稀疏。冬季垂直混合強(qiáng), 等鹽度線(xiàn)豎直分布, 斷面南部黃河入?？诟浇}度仍較低, 此處等鹽度線(xiàn)密集, 水平梯度最大; 斷面中部鹽度較均勻, 鹽度也較高, 到斷面北部靠岸處鹽度又逐漸降低, 水平梯度變大。

2.2 各斷面鹽度距平的經(jīng)驗(yàn)正交函數(shù)分解

分別對(duì) 3個(gè)斷面的鹽度距平進(jìn)行了經(jīng)驗(yàn)正交函數(shù)(EOF)分解, 表1中是其前3個(gè)模態(tài)的特征值及方差 貢 獻(xiàn) 。 根 據(jù) 相 鄰 特 征 值 滿(mǎn) 足 ：λj-λj+1≥ej,是樣本個(gè)數(shù))時(shí), 認(rèn)為這兩個(gè)特征值對(duì)應(yīng)的 EOF是有價(jià)值的信號(hào)的顯著性檢驗(yàn)條件, 前 3個(gè)特征向量都是可信的。前 3個(gè)模態(tài)累積方差貢獻(xiàn)率達(dá)到 75%以上, 可以代表原始場(chǎng)的主要特征, 是各斷面鹽度距平場(chǎng)的3個(gè)基本形態(tài)。

表1 各斷面鹽度距平的前3個(gè)特征向量對(duì)應(yīng)的特征值及其方差貢獻(xiàn)Tab.1 Eigenvalues and variance contributions of the first three eigenvectors associated with the deviations of salinity

將其第一模態(tài)時(shí)間變化曲線(xiàn)與臨近的海洋站鹽度作比較(圖3), 可以發(fā)現(xiàn)之間有很好的相關(guān)關(guān)系,說(shuō)明其可以反映鹽度的變化趨勢(shì), 即整體趨勢(shì)都是升高。并且渤海的時(shí)間變化曲線(xiàn)與黃河流量變化有很好的反相關(guān)關(guān)系, 說(shuō)明黃河流量的減小是造成渤海鹽度升高的重要原因。

2.3 千里巖鹽度的長(zhǎng)期變化趨勢(shì)及其成因簡(jiǎn)析

對(duì)千里巖站的鹽度進(jìn)行 2年濾波分析(圖 4a)可以發(fā)現(xiàn), 40 a間千里巖站的鹽度長(zhǎng)期變化趨勢(shì)是上升,鹽度年變率為0.01 /a, 并且存在大約5 a和10 a的周期變化。

鹽度與降水有著直接聯(lián)系, 而由于缺少千里巖站的降水資料, 作者與距離該站最近的青島的降水作了比較, 發(fā)現(xiàn)兩者存在比較明顯的負(fù)相關(guān), 并估計(jì)1960～2000年的 40 a間由于降水的減少, 可以造成大約 0.19的鹽度升高, 大概占到鹽度實(shí)際升高的 50%,因此降水的減少是造成其鹽度升高的主要原因。而由于ENSO對(duì)降水有影響, 一般的, El Ni?o年中國(guó)北方降水將減少[10,11], 因此把千里巖站的月平均鹽度和Nin?3指數(shù)標(biāo)準(zhǔn)化之后進(jìn)行比較(圖4b), 可以發(fā)現(xiàn)兩者之間存在著很好的正相關(guān)關(guān)系, 即El Ni?o年鹽度升高,與降水減少也是一致的, 但是經(jīng)過(guò) 2年濾波后則對(duì)應(yīng)得不好, 說(shuō)明千里巖站鹽度的短期震蕩可能和El Ni?o的發(fā)生有關(guān)。再將千里巖站 2年濾波后的鹽度和年平均太平洋年代際震蕩(PDO)指數(shù)比較(圖4c), 發(fā)現(xiàn)其與PDO的變化趨勢(shì)存在較好的正相關(guān)關(guān)系, 這說(shuō)明其年際以上尺度變化和外洋傳來(lái)的信號(hào)有關(guān)系, 而導(dǎo)致這種長(zhǎng)期變化相關(guān)關(guān)系的原因還有待進(jìn)一步研究。

圖3 各斷面鹽度距平EOF分解的第一模態(tài)時(shí)間變化曲線(xiàn)(粗實(shí)線(xiàn))Fig. 3 The first temporal EOFs of the deviation of salinity (thick solid line) in each section

2.4 千里巖站與其他海洋站鹽度之間的變化關(guān)系和超前/滯后分析

再將千里巖站和煙臺(tái)、塘沽、秦皇島站的鹽度作比較(圖5), 可以發(fā)現(xiàn)4站的變化趨勢(shì)基本相同,僅存在幅度上的偏差, 4站都在1970年、1984年、1995年附近出現(xiàn)最大值, 而總的趨勢(shì)也都是升高。

將千里巖的鹽度和其他海洋站的鹽度做超前滯后分析, 將千里巖的鹽度記為mx, 在葫蘆島、老虎灘、鲅魚(yú)圈、塘沽、秦皇島、煙臺(tái)海洋站中任取一個(gè)站的鹽度, 記為my, 然后對(duì)mlx+和my作相關(guān)分析, 其中l(wèi)= －12, －11, …, 0, …, 11, 12代表時(shí)間滯后(負(fù)值代表超前)年數(shù), 計(jì)算相關(guān)系數(shù) ：

其中, ∑表示對(duì)xm+l和ym均有數(shù)據(jù)的資料求和,xm+l和ym是這些數(shù)據(jù)的平均值, 所得結(jié)果如圖6。

圖4 千里巖鹽度與青島降水、Ni?o3指數(shù)、PDO指數(shù)時(shí)間變化圖Fig. 4 The correlations of the salinity at Qianliyan between the Qingdao precipitation, Ni?o3 index, and PDO indexa. 青島降水(2年濾波); b. Ni?o3指數(shù); c. PDO指數(shù)(2年濾波)a. Qingdao precipitation (2 years filter); b. Ni?o3 index; c. PDO index (2 years filter)

由圖6可以看到, 千里巖站與環(huán)黃、渤海海洋站的鹽度有較好的相關(guān)關(guān)系。其中, 千里巖與老虎灘、秦皇島和塘沽的相關(guān)系數(shù)以 1l=-時(shí)為最高, 分別是0.352, 0.389和0.374, 與葫蘆島和鲅魚(yú)圈的相關(guān)系數(shù)以 3l=- 時(shí)為最高, 分別是0.319和0.373, 與煙臺(tái)的相關(guān)系數(shù)以 4l=- 時(shí)最高, 為0.299。由此看來(lái), 千里巖站的鹽度比其他環(huán)黃、渤海海洋站的鹽度都有一定的超前, 時(shí)間一般為1～3個(gè)月。這也表明了渤海鹽度的升高可能與黃海的入侵有一定的關(guān)系。

3 結(jié)論

圖5 環(huán)黃、渤海海洋站鹽度時(shí)間變化圖Fig. 5 The variations of the salinity at stations around the Yellow Sea and the Bohai Sea

圖6 千里巖與其他海洋站鹽度超前滯后相關(guān)Fig. 6 The lagging correlation analysis between the salinity values at Qianliyan and other stations

(1)36°N斷面、大連—成山頭斷面、渤海中部斷面的鹽度的整體趨勢(shì)都是升高。黃河流量的減小是造成渤海鹽度升高的重要原因。(2)黃海千里巖站的鹽度40年間總體呈上升趨勢(shì), 年變率為0.01 a-1, 降水減少是造成鹽度上升的主要因素。并且鹽度存在著5a和10a左右的周期變化, 其季節(jié)短期震蕩和El Ni?o有關(guān)系, 而長(zhǎng)期變化與PDO有關(guān), 說(shuō)明年際以上尺度變化和外洋傳來(lái)的信號(hào)有關(guān)系。(3)千里巖站與其他環(huán)渤、黃海海洋站的鹽度成正相關(guān)關(guān)系,它們都有 5～6a的周期, 總體趨勢(shì)都是升高。而且千里巖站的鹽度比鲅魚(yú)圈、葫蘆島、秦皇島、塘沽、煙臺(tái)、老虎灘6站的鹽度都有一定的超前, 表明渤海鹽度的升高可能與黃海的入侵有一定的關(guān)系。

[1] 馮士筰, 李鳳岐, 李少菁. 海洋科學(xué)導(dǎo)論 [M]. 北京：高等教育出版社, 1999.

[2] 劉哲, 魏皓, 蔣松年. 渤海多年月平均溫鹽場(chǎng)的季節(jié)變化特征及形成機(jī)制的初步分析 [J].青島海洋大學(xué)學(xué)報(bào), 2003, 33(1)： 7-14.

[3] Lin Chuanlan, Su Jilan, Xu Bingrong,et al. Long-term variations of temperature and salinity of the Bohai Sea and their influence on its ecosystem [J]. Progress in Oceanography, 2001, 49： 7-19.

[4] 方國(guó)洪, 王凱, 郭豐義, 等. 近30年渤海水文和氣象狀況的長(zhǎng)期變化及其相互關(guān)系 [J].海洋與湖沼, 2002,33(5)： 515-525.

[5] Bai Hong, Hu Dunxin, Chen Yongli,et al. Statistic characteristics of thermal structure in the southern Yellow Sea in summer [J].Chinese Journal of Oceanology and Limnology, 2004, 22(3)： 237-243.

[6] 吳德星, 牟林, 李強(qiáng), 等. 渤海鹽度長(zhǎng)期變化特征及可能的主導(dǎo)因素 [J].自然科學(xué)進(jìn)展, 2004, 14(2)：191-195.

[7] 吳德星, 萬(wàn)修全, 鮑獻(xiàn)文, 等. 渤海 1958年和 2000年夏季溫鹽場(chǎng)及環(huán)流結(jié)構(gòu)的比較 [J].科學(xué)通報(bào), 2004,49(3)： 287-292.

[8] 劉用強(qiáng), 丁一匯. ENSO事件對(duì)我國(guó)季節(jié)降水和溫度的影響 [J].大氣科學(xué), 1995, 19(2)： 200-208.

[9] 龔道溢, 王紹武. 近百年ENSO對(duì)全球陸地及中國(guó)降水的影響[J]. 科學(xué)通報(bào), 1999, 44(3)： 315-320.

The spatial distribution and temporal evolution of salinity at the sections and observation stations in the Yellow Sea and Bohai Sea

MA Chao1, JU Xia2, WU De-xing1, LIN Xiao-pei1
(1. College of Physical and Environmental Oceanography, Ocean University of China, Qingdao 266100, China; 2.First Institute of Oceanography, State Oceanic Administration, Qingdao 266061, China)

Dec., 16, 2008

the Yellow Sea, the Bohai Sea, salinity, variation

Empirical Orthogonal Function (EOF) analysis and lagging correlation analysis were adopted to analyze the spatial distribution and temporal evolution of salinity at the 36°N section, Dalian-Chengshantou section and central Bohai section, and to describe the relationship among the observation stations around the Yellow Sea and Bohai Sea. The conclusion is that the salinity at Qianliyan station in the Yellow Sea increased during the period of 1960～2000 with an annual rate of 0.01/a, and it was probably caused by the decrease of precipitation. The short-term oscillation of the salinity might related to the El Ni?o, while the long-term variation related to the Pacific Decal Oscillation (PDO). Besides, there were positive correlations between the salinity at Qianliyan and other stations, and the only difference was a temporal lead of 1～3 months. This study systematically analyzed the long term variation of the salinity in the Yellow and Bohai Seas, and showed that the salinity in coastal sea might be affected by the open ocean, which is helpful to understand the evolution of marine environment in the China coastal sea.

P722.4; P722.5

A

1000-3096(2010)09-0070-06

2008-12-16;

2009-12-30

國(guó)家重點(diǎn)基礎(chǔ)研究發(fā)展計(jì)劃項(xiàng)目(2005CB422302); 國(guó)家自然科學(xué)基金項(xiàng)目(40930844); 國(guó)家科技支撐計(jì)劃項(xiàng)目(2008BAC49B02)

馬超(1983-), 男, 山東濟(jì)寧人, 博士, 主要從事中國(guó)東海環(huán)境演變研究, 電話(huà)： 0532-66782853, E-mail： machao@ouc.edu.cn

(本文編輯： 劉珊珊)