張緒東，張菀倫，李云波

（1.中國(guó)海洋大學(xué)，山東青島266003；2.中國(guó)科學(xué)院海洋環(huán)流與波動(dòng)重點(diǎn)實(shí)驗(yàn)室，山東青島266071；3.海軍海洋水文氣象中心，北京100161）

北黃海海溫分布變化的數(shù)值模擬分析

張緒東1,3，張菀倫2，李云波3

（1.中國(guó)海洋大學(xué)，山東青島266003；2.中國(guó)科學(xué)院海洋環(huán)流與波動(dòng)重點(diǎn)實(shí)驗(yàn)室，山東青島266071；3.海軍海洋水文氣象中心，北京100161）

通過(guò)對(duì)各月北黃海海溫統(tǒng)計(jì)產(chǎn)品的分析，研究了北黃海海溫的分布特征及其變化特點(diǎn)，在對(duì)北黃海冷水現(xiàn)象的分析研究后，對(duì)于海溫的季節(jié)變化特點(diǎn)也進(jìn)行了分析闡述。通過(guò)ROMS海洋模型對(duì)海溫進(jìn)行數(shù)值模擬分析研究，我們對(duì)于北黃海海溫的詳細(xì)特征和變化機(jī)制有了更深刻的把握，認(rèn)識(shí)到了熱平流、上升流、混合機(jī)制等因素是形成北黃海海溫區(qū)域分布特點(diǎn)的重要因素。

北黃海；平均水溫；冷水團(tuán)；黃海暖流；ROMS海洋模型；變化機(jī)制

1 引言

北黃海位于渤海海區(qū)的東側(cè)，成山角至朝鮮半島的長(zhǎng)山島以北的海區(qū)，面積約7.1×104km2。黃海地貌形態(tài)較復(fù)雜，最突出的特征有黃海槽、潮流脊和水下階地。注入北黃海的河流，中國(guó)沿岸主要有鴨綠江，朝鮮沿岸有大同江等。

北黃海屬溫帶海洋性氣候，季風(fēng)明顯，潮汐性質(zhì)復(fù)雜，潮流顯著，加之特殊的海底地形、陸地徑流和受外部海域環(huán)境影響等因素造成了北黃海海洋環(huán)境狀況的復(fù)雜，但其海洋環(huán)境要素仍具有顯著的季節(jié)分布特征和年變化規(guī)律。

從20世紀(jì)70年代至今，國(guó)內(nèi)對(duì)黃海海溫的研究，從海洋調(diào)查資料分析開(kāi)始，完成了海溫分布特征和黑潮分支黃海暖流對(duì)該海區(qū)溫度分布的影響研究，重點(diǎn)是對(duì)黃海冷水團(tuán)現(xiàn)象進(jìn)行了研究。80年代以后開(kāi)始了數(shù)值模擬的研究，對(duì)黃海溫鹽變化的研究考慮了較多的因素，受觀測(cè)資料的限制，有的研究側(cè)重于某類動(dòng)力過(guò)程對(duì)溫度的影響，如黃海水溫與環(huán)流變化特征研究[1]、風(fēng)浪及地形對(duì)渤、黃海夏季溫度垂直結(jié)構(gòu)的影響[2]等，有的研究側(cè)重于研究對(duì)溫度水平分布或垂直結(jié)構(gòu)分布的變化等，多限于針對(duì)某個(gè)季節(jié)的分析研究，如渤黃海溫度躍層的分布與變化[3]、黃海西部沿岸冷水在夏季南黃海西部底層冷水形成和季節(jié)演變過(guò)程中的作用[4]等。

近年來(lái)利用衛(wèi)星資料和更優(yōu)化的模式，對(duì)黃海溫度也有新的研究，主要為對(duì)黃海海表溫度的年際與年代際變化、潮生陸架鋒及斷面水溫分布等數(shù)值模擬研究，特別是風(fēng)、浪、流耦合數(shù)值模式的應(yīng)用使得近海三維溫鹽流的預(yù)報(bào)進(jìn)入了業(yè)務(wù)預(yù)報(bào)時(shí)期。

國(guó)外相關(guān)研究較少，常見(jiàn)的研究成果主要是對(duì)黃海冷水團(tuán)的若干研究和海流方面的模擬分析，包括大區(qū)域海洋數(shù)值模擬結(jié)果的改進(jìn)和數(shù)據(jù)同化，對(duì)黃海冷水團(tuán)的變性研究興趣較高。

本文以北黃海海域海溫歷史資料為基礎(chǔ)，利用數(shù)值模式計(jì)算和統(tǒng)計(jì)分析相結(jié)合的技術(shù)方法，對(duì)各類資料進(jìn)行客觀分析，系統(tǒng)、科學(xué)的研究北黃海水文氣象要素及海洋現(xiàn)象的時(shí)空分布規(guī)律和特征。通過(guò)對(duì)資料的處理分析、信息提取，系統(tǒng)的研究了北黃海海水溫度的分布特征及變化規(guī)律，綜合分析北黃海海水溫度區(qū)域分布變化特點(diǎn)的影響機(jī)制，為該海域的各類海洋活動(dòng)提供技術(shù)理論和環(huán)境特征信息，達(dá)到為生產(chǎn)、軍事活動(dòng)服務(wù)保障的目標(biāo)。

2 北黃海典型斷面分布特征

在廣泛收集北黃海海域國(guó)內(nèi)外海洋水文觀測(cè)資料和專項(xiàng)調(diào)查資料的基礎(chǔ)上，采用質(zhì)控處理方法，主要包括：范圍檢驗(yàn)、非法碼檢驗(yàn)、合理性檢驗(yàn)、遞增性檢驗(yàn)、連續(xù)性檢驗(yàn)、穩(wěn)定性檢驗(yàn)、排重檢驗(yàn)等進(jìn)行了資料的標(biāo)準(zhǔn)化處理。垂向數(shù)據(jù)插值計(jì)算，主要采用了阿基瑪（Akima）插值法[5]。

數(shù)據(jù)分析系統(tǒng)的格點(diǎn)化處理方法主要采用了INVERSEDISTANCEA POWER，該方法的基本原理是：

式中：hij是網(wǎng)格點(diǎn)“j”與鄰點(diǎn)“i”之間的有效距離；Z?j是網(wǎng)格點(diǎn)“j”的內(nèi)插值；Zi是相鄰點(diǎn)的要素值；dij是網(wǎng)格點(diǎn)“j”與鄰點(diǎn)“i”之間的距離；δ是權(quán)重參數(shù)；β是平滑參數(shù)。

由于篇幅所限，本節(jié)內(nèi)容選取由遼寧大連至山東成山頭附近海域（經(jīng)北黃海中部）斷面（簡(jiǎn)稱為“斷面A”）水溫分布為代表，分析其季節(jié)變化特征。斷面A的水溫?cái)嗝娣植家?jiàn)圖1。

圖1表明，A斷面在冬季（2月）混合最徹底，由表層可直達(dá)海底，深水區(qū)溫度混合均勻，北黃海中部海槽斷面顯示該處水體平均水溫較高，在5℃左右，向兩端淺水區(qū)，水溫逐漸遞減，成山頭沿岸海域水體溫度在0—2℃左右，西北部淺水區(qū)水體溫度在0℃左右。

圖1 北黃海斷面水溫垂直變化季節(jié)特征

春季（5月）、夏季（8月）A斷面水溫垂直分層結(jié)構(gòu)明顯，但在下層都有冷水團(tuán)盤(pán)踞在海槽洼地處，春季（5月）在30m層以下至槽底有小于6℃的冷水，夏季（8月）在30m層以下至海槽西南側(cè)的底部有小于8℃的冷水盤(pán)踞。

春季（5月）A斷面溫度躍層明顯，但溫度躍層在夏季（8月）強(qiáng)度最強(qiáng)，特別在10—30m層的溫度梯度最大，最大溫差在10℃以上。

秋季（11月）A斷面出現(xiàn)上混合層，30m層以淺溫度分布為垂直均勻型，30m層的中部海區(qū)至海槽偏西南側(cè)底層有溫度低值區(qū)存在，直至冬季該冷水區(qū)消失。

3 北黃海海溫的數(shù)值模擬分析

本文采用區(qū)域海洋海流模式ROMS（Regional Ocean M odel System)，ROMS是自由表面的、基于流體力學(xué)原始方程的靜力平衡模式。

3.1 模式配置

3.1.1 資料的選取與處理

在ROMS模式基礎(chǔ)上，我們構(gòu)建了適合中國(guó)陸架邊緣海的高分辨率數(shù)值模型。模式計(jì)算的海區(qū)為東中國(guó)海海域（24°20′—40°55′N,117°5′—130°20′E），水平分辨率為5′×5′，網(wǎng)格個(gè)數(shù)為160×200。模式垂向采用S坐標(biāo)，垂直分16層。模式中外模時(shí)間步長(zhǎng)為30 s，內(nèi)模時(shí)間步長(zhǎng)為600 s，模式積分時(shí)間從1958年1月1日—2005年12月31日。潮汐作用只考慮M 2分潮，通過(guò)潮汐水位和潮汐流速?gòu)拈_(kāi)邊界引入。模式的東南邊界為開(kāi)邊界，西北邊界為閉邊界。模式水深為實(shí)際海底地形，取自方國(guó)洪2002年調(diào)查的2′×2′的水深資料，插值到模式網(wǎng)格中。模式中保留了所有的島嶼，將水深大于1200 m處作為1200m處理。

3.1.2 邊界條件

初始溫度、鹽度場(chǎng)取自全球簡(jiǎn)單海洋資料同化分析系統(tǒng)（Simple Ocean Data Assim ilation,SODA) 1958年資料。上表面風(fēng)應(yīng)力數(shù)據(jù)取NCEP1958—2005年每天的數(shù)據(jù)資料插值到模式網(wǎng)格點(diǎn)上。短波輻射和向下的長(zhǎng)波輻射數(shù)據(jù)（NCEP 1958—2005年每天的數(shù)據(jù)資料插值）計(jì)算上表面熱通量。實(shí)驗(yàn)結(jié)果表明用這種方法計(jì)算出的溫度優(yōu)于直接用凈長(zhǎng)波和短波輻射作為強(qiáng)迫時(shí)的模擬結(jié)果。在模式中考慮了作為東中國(guó)海最主要的淡水注入的黃河和長(zhǎng)江這兩大河流，資料來(lái)源于利津和大通站的河口流量檢測(cè)數(shù)據(jù)。模式輸出為1958年1月1日至2005年12月31日，模擬結(jié)果用于分析中國(guó)陸架邊緣海的環(huán)流及其變異。

3.2 模式結(jié)果驗(yàn)證

3.2.1 衛(wèi)星遙感觀測(cè)數(shù)據(jù)驗(yàn)證

對(duì)于海表SST的大面模擬結(jié)果驗(yàn)證，我們選取NOAA/NASA高分辨率輻射計(jì)（Advanced Very High Resolution Radiometers,AVHRR）衛(wèi)星觀測(cè)進(jìn)行對(duì)比分析（數(shù)據(jù)下載網(wǎng)址:http://iridl.ldeo. columbia.edu/SOURCES/.NOAA/.NCDC/.OISST/. version2/.AVHRR/.sst/）。由于數(shù)值模擬的精度為1/ 12°，我們選擇9 km精度的相等角度網(wǎng)格觀測(cè)數(shù)據(jù)，這樣遙感數(shù)據(jù)精度和模式的數(shù)據(jù)精度較為匹配[6]。盡管如此，NOAA系列衛(wèi)星攜帶的是紅外遙感輻射計(jì)，除了受觀測(cè)海域上空云的影響，還存在缺測(cè)現(xiàn)象，雖然對(duì)觀測(cè)數(shù)據(jù)進(jìn)行了插值處理，但資料仍出現(xiàn)小面積SST偏低的現(xiàn)象。

論文選用1989年2月、5月、8月、11月代表各個(gè)季節(jié)SST進(jìn)行對(duì)比分析。從圖2對(duì)比可以看出，模式結(jié)果和衛(wèi)星遙感觀測(cè)數(shù)據(jù)基本一致。冬季渤海、北黃海以及蘇浙沿岸溫度較低，海表溫度在10℃以下，25℃的高溫黑潮水入侵東中國(guó)海，在濟(jì)州島以南入侵黃海，形成黃海暖流，其尾端可達(dá)渤海海峽。溫度從南至北逐漸降低，南北溫差在一年中最大，可達(dá)25℃，這主要是由于冬季太陽(yáng)總輻射量緯向變化最大。春季隨著太陽(yáng)輻射的加強(qiáng)，海表溫度升高，渤海和北黃海海域水溫在15℃左右，仍然低于其他海域。進(jìn)入夏季后，強(qiáng)烈的太陽(yáng)輻射使整個(gè)海區(qū)海表溫度分布趨于一致，全海域的海表溫度高于24℃，海域南北溫差低于5℃。秋季是海域溫度由夏季型分布向冬季型分布特征轉(zhuǎn)換的過(guò)渡期，也是海表溫度下降最快的季節(jié)，海表溫度分布逐漸形成了由南向北遞減的態(tài)勢(shì)。

3.2.2 斷面CTD觀測(cè)數(shù)據(jù)驗(yàn)證

圖2 36°N斷面1983年4個(gè)月份月平均SST

圖2 （續(xù)）

圖3 大連-成山頭斷面1989年4個(gè)月份溫度模擬（左）和觀測(cè)（右）

除了海洋觀測(cè)站SST的對(duì)比，我們也選取了兩個(gè)斷面的觀測(cè)資料與模擬結(jié)果進(jìn)行對(duì)比分析。觀測(cè)

數(shù)據(jù)取自36°N斷面、大連-成山頭斷面（從37.25°N, 122.75°E—38.79°N，121.76°E）1976—1999年的標(biāo)準(zhǔn)層次CTD調(diào)查資料。36°斷面包含11個(gè)觀測(cè)點(diǎn)，大連-成山頭斷面包含9個(gè)觀測(cè)點(diǎn)。首先對(duì)各斷面的資料按照月份進(jìn)行歸類，然后對(duì)其進(jìn)行垂向插值，之后將各月的數(shù)據(jù)做平均，得到各個(gè)斷面月均溫度垂直結(jié)構(gòu)。

（1）大連成山頭斷面溫度

大連-成山頭斷面基本位于渤海和黃海分界處，地理位置十分重要。以1989年為例，選2月、5月、8月、11月分別代表四季，選取1989年的斷面溫度進(jìn)行對(duì)比分析（如圖3所示）。

圖3（續(xù)）

圖3 表明，實(shí)測(cè)資料分析海水溫度在垂直方向上下混合較均勻，趨于一致，水溫在水平方向差異變化也較小，斷面溫度呈柱狀分布特征，海槽中部溫度略高。模式結(jié)果同樣呈現(xiàn)出這樣的分布特征，說(shuō)明模式的物理機(jī)制是正確的，僅溫度模擬結(jié)果略高于實(shí)測(cè)值，存在計(jì)算誤差。

（2）36°N斷面溫度

圖4 36°N斷面1983年4個(gè)月份溫度模擬（左）和觀測(cè)（右）

選用1982年2月、5月、8月、11月的36°斷面觀測(cè)溫度對(duì)模擬結(jié)果進(jìn)行驗(yàn)證，對(duì)比結(jié)果如圖4所示。從圖中可以看出，冬季斷面溫度從表到底混合均勻，大約在5℃左右，隨著氣溫的升高，5月海表溫度開(kāi)始上升至15℃左右，底層溫度在10℃左右，

表底溫差不大。夏季隨著太陽(yáng)輻射的增強(qiáng)，海水升溫至一年的最高值，海表溫度高達(dá)30℃，海底溫度仍舊在10℃左右，海水的層化現(xiàn)象非常明顯。從海面到10m左右為混合層，這一層海水溫度均勻混合?；旌蠈酉聻闇剀S層，溫躍層內(nèi)溫度的垂向梯度很大。溫躍層下是潮混合層，潮流混合在這一層里起主導(dǎo)作用，使得溫度垂向均勻分布，形成略凸起的臺(tái)狀結(jié)構(gòu)。秋季，斷面溫度下降，海表溫度降低至20℃以下，而底層水溫仍保持10℃左右。

圖4 （續(xù)）

4 小結(jié)

論文運(yùn)用ROMS海洋模型，加入了潮汐、風(fēng)場(chǎng)、熱通量、邊界流速和水位、徑流以及斜壓作用，對(duì)東中國(guó)海海域進(jìn)行了數(shù)值模擬，模擬時(shí)間為1958—2005年。通過(guò)模式結(jié)果與海洋觀測(cè)站、CTD斷面觀測(cè)資料和衛(wèi)星資料的對(duì)比，結(jié)果顯示模擬結(jié)果和觀測(cè)數(shù)據(jù)在長(zhǎng)期變化趨勢(shì)上具有一致性，并且季節(jié)分布結(jié)構(gòu)相同?？梢钥闯?，在區(qū)域海洋模式ROMS的模擬過(guò)程中，由于模式中考慮了區(qū)域水平對(duì)流、垂相混合，并使得溫躍層和底邊階層的解析度提高，這些在模式結(jié)果與實(shí)測(cè)資料以及統(tǒng)計(jì)分析結(jié)果的比較中都得到了很好的驗(yàn)證。

海水溫度變化盡管量值很小，也常常會(huì)隨著氣溫變化而發(fā)生，但這種變化主要集中在北黃海的春季和夏季。秋季北黃海海水溫度的降低多半與氣溫變化無(wú)關(guān)，因?yàn)楹Ｑ蟮膬魺嵬繛樨?fù)值，其主要原因是由海表面的水和更深冷水層之間發(fā)生了混合。在冬季海溫明顯降低后，黃海暖流的熱平流作用主導(dǎo)了該區(qū)域的溫度分布場(chǎng)。

[1]袁承儀.黃海水溫與環(huán)流變化特征及機(jī)制研究[D].青島:中國(guó)海洋大學(xué),2011.

[2]徐玲玲,林霄沛,吳德星.風(fēng)浪及地形對(duì)渤、黃海夏季溫度垂直結(jié)構(gòu)的影響[J].中國(guó)海洋大學(xué)學(xué)報(bào),2008,38(2):183-188.

[3]屠金釗.渤、黃海溫度躍層的分布與變化[J].海洋通報(bào),1992,11 (4):27-32.

[4]韋久勝,于志剛,葛人峰,等.黃海西部沿岸冷水在夏季南黃海西部底層冷水形成和季節(jié)演變過(guò)程中作用的化學(xué)水文分析[J].海洋與湖沼,2013,44(4):890-905.

[5]陳上及,馬繼瑞.海洋數(shù)據(jù)處理分析方法及其應(yīng)用[M].北京:海洋出版社,1991.

[6]Oey L-Y.Simulation of mesoscale variability in the gulf of Mexico:Sensitivity studies,comparison w ith observations,and trapped wave propagation[J].Journal of Physical Oceanography, 1996,26(2):145-175.

Characteristicsof the sea tem perature in the North Yellow Sea

ZHANG Xu-dong1,3，ZHANGWan-lun2，LIYun-bo3
（1.The Ocean University ofChina,Qingdao 266003 China;2.Key Laboratory ofOcean Circulation and Waves,Institute of Oceano logy, Chinese Academy ofScience,Qingdao 266071 China;3.Navy Marine HydrometeorologicalCenter,Beijing 100161 China）

The distribution characteristics and its variations of sea temperature in the Northern Yellow Sea was studied,the evolution of sea temperature was exam ined through the statistical analysismonthly sea temperature data in many years.Because of the identification of cold watermass in the Yellow Sea,its seasonal variations about the sea temperaturewere led,and the influence of the Yellow SeaWarm Current to the sea temperaturewas recognized.The detailed characteristics and rule of sea temperature was well known.It’s recognized that the warm current,upwelling and m ixing effect etc.are themost important factors resulting in the area distribution charactersof the sea temperature in the North Yellow Sea.

Northern Yellow Sea;average sea temperature;cold watermass;Huanghaiwarm current;ROMS oceanmodel;variations

P731.11

A

1003-0239（2015）05-0089-09

2014-12-17

國(guó)家自然科學(xué)基金（41276026）

張緒東（1975—），男，工程師，主要從事海洋環(huán)境預(yù)報(bào)和研究工作。E-mail：jerryzhangwww@126.com

10.11737/j.issn.1003-0239.2015.05.011