黃艷松,宋金寶

(1.中國(guó)科學(xué)院 海洋研究所 環(huán)流與波動(dòng)重點(diǎn)實(shí)驗(yàn)室,山東 青島 266071; 2.中國(guó)科學(xué)院 研究生院,北京100039)

海-氣熱交換在海氣相互作用過(guò)程中起著非常重要的作用。海表凈熱通量(Fnh)包括向下短波輻射通量(Fds)、向下長(zhǎng)波輻射通量(Fdl)、向上短波輻射通量(Fus)、向上長(zhǎng)波輻射通量(Ful)、感熱通量(Hs)和潛熱通量(HL)。這些通量被用來(lái)驅(qū)動(dòng)海洋模式并檢驗(yàn)氣象預(yù)報(bào)模式結(jié)果的可信度。但由于現(xiàn)場(chǎng)觀測(cè)資料的匱乏,我們對(duì)這些通量參數(shù)了解很少。目前使用最廣泛的用于為海洋、大氣模式提供初始場(chǎng)且驗(yàn)證模式結(jié)果正確與否的通量產(chǎn)品來(lái)自于 NCEP再分析資料。NCEP再分析資料是由美國(guó)國(guó)家環(huán)境預(yù)報(bào)中心(NCEP)和國(guó)家大氣研究中心(NCAR)協(xié)作,對(duì)來(lái)源于地面、船舶、無(wú)線電探空、探空氣球、飛機(jī)、衛(wèi)星等氣象觀測(cè)資料進(jìn)行同化處理后,研制的全球氣象資料數(shù)據(jù)庫(kù),包含一天4次、日平均和月平均3種時(shí)間格式的數(shù)據(jù)。NCEP再分析資料分為兩個(gè)版本,即NCEP1(NCEP/NCAR Reanalysis 1)和NCEP2(NCEP/DOE Reanalysis 2)。NCEP2是NCEP1的升級(jí)版,兩者具有相同的分辨率且用的是相同的原始觀測(cè)數(shù)據(jù),其感熱和潛熱通量算法相同,但兩者在涉及到輻射通量的一些參數(shù)上、同化數(shù)據(jù)源和同化方法是不一樣的[1-3]。關(guān)于NCEP再分析資料在世界大洋各個(gè)區(qū)域的可信性評(píng)估前人已經(jīng)作了很多工作,包括：阿拉伯海[4]、熱帶太平洋[5-7]、東北大西洋[8-9]、拉布拉多海[10]、北大西洋的西邊界流[11]、南半球[12]、印度洋[13-14],中國(guó)南海[15]和中國(guó)東海[16]。NCEP再分析資料在不同區(qū)域表現(xiàn)不同,文獻(xiàn)[17]根據(jù)東南阿拉伯海的浮標(biāo)資料結(jié)果指出 NCEP再分析資料低估了海表氣象學(xué)參數(shù),進(jìn)而低估了海表的感熱通量和潛熱通量; 文獻(xiàn)[4]應(yīng)用阿拉伯海上另一點(diǎn)觀測(cè)資料,文獻(xiàn)[9]應(yīng)用東北大西洋的浮標(biāo)觀測(cè)資料分別指出NCEP再分析資料計(jì)算的湍流熱通量存在高估現(xiàn)象。比較研究表明 NCEP再分析資料與現(xiàn)場(chǎng)觀測(cè)資料計(jì)算結(jié)果之間的差異在不同的區(qū)域甚至同一區(qū)域不同時(shí)間內(nèi)都是表現(xiàn)不同的[18-19]。

在本文研究中,我們對(duì) NCEP再分析資料中的海表氣象學(xué)參數(shù)及其海表熱通量與布放在黃海北部的浮標(biāo)觀測(cè)資料及通量計(jì)算結(jié)果進(jìn)行了比較。對(duì)其可信度進(jìn)行了初步檢驗(yàn)和評(píng)估,以便為研究中國(guó)近海選用合適的再分析資料提供借鑒和參考。

1 觀測(cè)數(shù)據(jù)及數(shù)據(jù)處理方法

本文資料對(duì)應(yīng)時(shí)間為2009年9月26日到2009年11月20日(世界時(shí))?，F(xiàn)場(chǎng)觀測(cè)資料來(lái)自于布放在黃海北部 122°44.355′E 和 39°02.908′N 的錨定浮標(biāo)觀測(cè)數(shù)據(jù),關(guān)于該浮標(biāo)的具體描述可見文獻(xiàn)[20]。浮標(biāo)上距離海面6.8 m左右的常規(guī)氣象站可以每分鐘輸出一次風(fēng)速(U)、風(fēng)向、大氣溫度(Tair)、大氣壓力(P)和相對(duì)濕度,應(yīng)用Tair、P和相對(duì)濕度,我們計(jì)算出比濕Q; 海面以下0.2 m左右的水溫儀每5 s觀測(cè)一次水溫,我們將其作為海表溫度Tsea; 另外,在距離海面4 m左右的浮標(biāo)桁架上還安裝了4個(gè)輻射儀,可以分別觀測(cè)Fds、Fdl、Fus和Ful,輻射通量也是每分鐘觀測(cè)一次。所有的數(shù)據(jù)在使用之前都經(jīng)過(guò)了質(zhì)量控制。6.8 m處常規(guī)風(fēng)速儀觀測(cè)到的風(fēng)速通過(guò)中性條件下對(duì)數(shù)定律[23]轉(zhuǎn)換成 10 m 風(fēng)速U10; 感熱通量和潛熱通量是通過(guò)COARE3.0塊體通量程序[22]計(jì)算得到的。為了與NCEP再分析資料進(jìn)行對(duì)比,我們對(duì)浮標(biāo)觀測(cè)到的海表氣象參數(shù)、輻射通量和湍流熱通量進(jìn)行日平均。

關(guān)于 NCEP 再分析資料日平均U10、Q、Tair、Tsea、Fds、Fdl、Fus、Ful、Hs和HL來(lái)自于兩個(gè)數(shù)據(jù)庫(kù)NCEP/NCAR Reanalysis 1與NCEP/DOE Reanalysis 2,分別命名為NCEP1和NCEP2。為了與浮標(biāo)數(shù)據(jù)進(jìn)行對(duì)比,我們采用最鄰近插值法將 NCEP再分析資料空間插值到錨定浮標(biāo)的位置。

為了更細(xì)致地比較 NCEP再分析資料與浮標(biāo)資料給出的海表氣象學(xué)參數(shù)和熱通量間的異同,我們考慮平均偏差、相對(duì)誤差、均方根誤差和相關(guān)系數(shù)4個(gè)統(tǒng)計(jì)量,計(jì)算公式如下：

其中Xi表示浮標(biāo)資料時(shí)間序列,Yi代表NCEP再分析資料時(shí)間序列,上劃線代表平均值,“| |”表示取絕對(duì)值,σ代表均方差。

2 結(jié)果

2.1 氣象學(xué)參數(shù)比對(duì)

圖1給出了來(lái)自于浮標(biāo)、NCEP1和 NCEP2的U10、Q、Tair、Tsea數(shù)據(jù)的日平均時(shí)間序列圖。可以看出 NCEP資料與浮標(biāo)觀測(cè)資料之間還是存在一定差異的。表1給出了NCEP資料與浮標(biāo)觀測(cè)資料中海表氣象參數(shù)及各種熱通量的對(duì)比,表 2給出了NCEP再計(jì)算熱通量(應(yīng)用NCEP1、NCEP2海表氣象資料而采用 COARE3.0塊體通量程序計(jì)算出的熱通量)與浮標(biāo)觀測(cè)資料計(jì)算結(jié)果的對(duì)比,包括平均偏差、相對(duì)誤差、均方根誤差和相關(guān)系數(shù),計(jì)算方法見公式(1)～(4),表格中的Fns代表凈短波輻射通量(Fns=Fus-Fds),Fnl代表凈長(zhǎng)波輻射通量(Fnl=Ful-Fdl),Fnr代表凈輻射通量(Fnr=Fns+Fnl)。

圖1 海表氣象參數(shù)日平均值的時(shí)間序列Fig.1 Time series plots of daily averaged values of marine meteorological parameters

圖1(a)為海面U10時(shí)間序列?？梢钥闯?NCEP再分析資料與浮標(biāo)觀測(cè)風(fēng)速變化趨勢(shì)基本一致(NCEP1、NCEP2與浮標(biāo)觀測(cè)值的相關(guān)系數(shù)分別為0.76,0.79); 在風(fēng)速小于11 m/s時(shí),NCEP1和NCEP2風(fēng)速資料一般要低于浮標(biāo)觀測(cè)風(fēng)速,NCEP2與浮標(biāo)觀測(cè)的風(fēng)速更接近; 當(dāng)風(fēng)速大于11 m/s時(shí),NCEP1和NCEP2一般要高于浮標(biāo)觀測(cè)風(fēng)速,NCEP1與浮標(biāo)觀測(cè)的風(fēng)速更接近。該結(jié)果與文獻(xiàn)[15,16]存在相似性。文獻(xiàn)[15]給出了中國(guó)南?，F(xiàn)場(chǎng)觀測(cè)資料和 NCEP2資料的對(duì)比,除科學(xué)3號(hào)比較結(jié)果外其他4組比較結(jié)果均在風(fēng)速小于10 m/s時(shí),NCEP2風(fēng)速資料要低于浮標(biāo)觀測(cè)風(fēng)速。文獻(xiàn)[16]應(yīng)用東海浮標(biāo)資料指出當(dāng)風(fēng)速處于10～20 m/s之間時(shí),NCEP再分析資料中的動(dòng)量通量要高于浮標(biāo)資料計(jì)算的結(jié)果,因?yàn)閯?dòng)量通量高也代表著風(fēng)速高,所以這一點(diǎn)與本文風(fēng)速的比較結(jié)果是一致的。對(duì)于整個(gè)比較期間,NCEP1風(fēng)速與浮標(biāo)風(fēng)速之間的平均誤差為-1.22 m/s,平均低估風(fēng)速達(dá)17.69%,NCEP2風(fēng)速與浮標(biāo)風(fēng)速之間的平均誤差為0.22 m/s,平均高估風(fēng)速為3.19%,NCEP2風(fēng)速的長(zhǎng)期平均值更接近于浮標(biāo)風(fēng)速的長(zhǎng)期平均值。

表1 NCEP再分析資料與浮標(biāo)觀測(cè)資料之間的比對(duì)Tab.1 Comparison of NCEP reanalysis and buoy data

表2 NCEP再計(jì)算湍流熱通量與浮標(biāo)觀測(cè)資料計(jì)算結(jié)果的對(duì)比Tab.2 Comparison of NCEP recomputed turbulence heat fluxes and that from buoy data

圖1(b)為海面以上2 m處Q日平均值的時(shí)間序列。NCEP1、NCEP2和浮標(biāo)的Q變化趨勢(shì)基本一致(NCEP1、NCEP2與浮標(biāo)觀測(cè)值的相關(guān)系數(shù)分別為0.96,0.96),NCEP1、NCEP2的Q在相對(duì)濕度較低時(shí)存在高估現(xiàn)象,NCEP2的Q更接近于浮標(biāo)測(cè)量值。Q高估會(huì)造成潛熱通量結(jié)果低估。文獻(xiàn)[25]也指出在地球熱帶之外,NCEP1、NCEP2的Q普遍存在高估現(xiàn)象。對(duì)于整個(gè)比較期間,NCEP1的Q與浮標(biāo)Q之間的平均誤差為 1.1 g/kg,平均高估Q達(dá) 16.18%,NCEP2的Q與浮標(biāo)Q之間的平均誤差為0.3 g/kg,平均高估Q為4.64%,NCEP2的Q的長(zhǎng)期平均值更接近于浮標(biāo)Q的長(zhǎng)期平均值。

圖1(c)為海面Tair日平均值時(shí)間序列。與Q比較結(jié)果相類似,NCEP1、NCEP2和浮標(biāo)的Tair結(jié)果變化趨勢(shì)基本一致(NCEP1、NCEP2與浮標(biāo)觀測(cè)值的相關(guān)系數(shù)分別為 0.99,0.99); 在開始風(fēng)速較弱的 6 d里,NCEP1、NCEP2的Tair存在低估現(xiàn)象; 6 d之后,NCEP1、NCEP2的Tair結(jié)果普遍存在高估現(xiàn)象,尤其是相對(duì)濕度比較低時(shí)更明顯; NCEP2的Tair結(jié)果更接近于浮標(biāo)觀測(cè)值。文獻(xiàn)[17]指出在季風(fēng)爆發(fā)前即弱風(fēng)下 NCEP再分析氣溫資料存在低估現(xiàn)象,這與本文在弱風(fēng)下得到的結(jié)論相一致。文獻(xiàn)[23]指出在地球熱帶之外,NCEP1、NCEP2的大氣溫度結(jié)果普遍存在低估現(xiàn)象,這與本文結(jié)果不太一致。對(duì)于整個(gè)比較期間,NCEP1 的Tair與浮標(biāo)Tair之間的平均誤差為 1.22℃,平均高估Tair達(dá) 9.96%,NCEP2 的Tair與浮標(biāo)Tair之間的平均誤差為 0.27℃,平均高估Tair為 2.24%,NCEP2 的Tair長(zhǎng)期平均值更接近于浮標(biāo)Tair的長(zhǎng)期平均值。

圖1(d)為海表水溫日平均值時(shí)間序列?？梢钥闯?NCEP1、NCEP2的Tsea資料與真實(shí)的Tsea總體變化趨勢(shì)相似(NCEP1、NCEP2與浮標(biāo)觀測(cè)值的R分別為 0.97,0.97); 在大部分時(shí)間里 NCEP1、NCEP2的Tsea資料存在高估現(xiàn)象; 9月末NCEP1、NCEP2的Tsea資料存在低估現(xiàn)象; NCEP2的Tsea結(jié)果更接近于浮標(biāo)觀測(cè)值。另外,NCEP再分析Tsea資料指的是皮溫,而我們浮標(biāo)觀測(cè)到的Tsea指的是水下0.2 m左右的海水溫度,關(guān)于這兩種溫度的關(guān)系在文獻(xiàn)[24]中有詳細(xì)介紹。對(duì)于整個(gè)比較期間,NCEP1Tsea與浮標(biāo)Tsea之間的平均誤差為 0.70℃,平均高估Tsea達(dá) 1.08%,NCEP2Tsea與浮標(biāo)Tsea之間的平均誤差為0.27℃,平均高估Tsea為0.34%,NCEP2Tsea的長(zhǎng)期平均值更接近于浮標(biāo)Tsea的長(zhǎng)期平均值。

2.2 通量對(duì)比

2.2.1 輻射通量

海面輻射通量由Fds、Fdl、Fus和Ful四部分組成。

圖2給出了來(lái)自于NCEP1、NCEP2和浮標(biāo)的輻射通量日平均值時(shí)間序列圖。圖2(a)為Fds日平均值時(shí)間序列,可以看出 NCEP1、NCEP2的Fds與浮標(biāo)觀測(cè)值之間存在一定的差異,普遍大于浮標(biāo)觀測(cè)值,存在高估現(xiàn)象(NCEP1、NCEP2與浮標(biāo)觀測(cè)值的平均誤差分別為 44.56,45.56 W/m2),變化趨勢(shì)大體一致(NCEP1、NCEP2與浮標(biāo)觀測(cè)值的相關(guān)系數(shù)分別為0.77,0.71); NCEP2的Fds相對(duì)于NCEP1改變不大,基本相等(NCEP1、NCEP2與浮標(biāo)觀測(cè)值的相對(duì)誤差分別為34.13%,34.89%)。圖2(b)為Fdl日平均值時(shí)間序列,可以看出NCEP1、NCEP2的Fdl普遍小于浮標(biāo)觀測(cè)值,存在低估現(xiàn)象(NCEP1、NCEP2與浮標(biāo)的平均誤差分別為-13.3,-26.07 W/m2),變化趨勢(shì)大體一致(NCEP1、NCEP2與浮標(biāo)觀測(cè)值的相關(guān)系數(shù)分別為0.90,0.93); NCEP1的Fdl大于NCEP2,與浮標(biāo)觀測(cè)值更接近(NCEP1、NCEP2與浮標(biāo)觀測(cè)值的相對(duì)誤差分別為-4.24%,-8.31%)。圖2(c)為USWR日平均值時(shí)間序列,可以看出 NCEP2的Fus與浮標(biāo)觀測(cè)值比較接近,而 NCEP1要比浮標(biāo)觀測(cè)值高很多 (NCEP1、NCEP2與浮標(biāo)的平均誤差分別為21.47,1.33 W/m2)。圖 2(d)為Ful日平均值時(shí)間序列,可以看出除了在 9月末之外,NCEP1、NCEP2的Ful普遍大于浮標(biāo)觀測(cè)值,存在高估現(xiàn)象(NCEP1、NCEP2與浮標(biāo)觀測(cè)值的平均誤差分別為9.15,10.94 W/m2); NCEP再分析Ful資料不能反映出真實(shí)Ful的巨大變化; NCEP2的Ful相對(duì)于NCEP1改變不大,基本相等(NCEP1、NCEP2與浮標(biāo)的相對(duì)誤差分別為 2.33%,2.78%)。從表 1中可以看出,NCEP再分析資料中的海面吸收的Fns、海面釋放的Fnl與浮標(biāo)觀測(cè)資料相比存在高估現(xiàn)象,這與文獻(xiàn)[9]中表 2b中的結(jié)果相一致; 而針對(duì)于Fnr而言,NCEP1,NCEP2的平均值都要高于浮標(biāo)觀測(cè)值且NCEP1的Fnr更接近浮標(biāo)觀測(cè)值(NCEP1、NCEP2與浮標(biāo)Fnr的平均誤差分別為0.65,7.22 W/m2)。

圖2 輻射通量日平均值的時(shí)間序列Fig.2 Time series plots of daily averaged values of radiative flux

2.2.2 湍流熱通量

海表面的湍流熱通量指的是Hs和HL。圖3給出了湍流熱通量日平均值的時(shí)間序列圖。該圖中正通量代表熱量從海洋傳到大氣; NCEP1和NCEP2代表NCEP再分析資料中直接下載的湍流熱通量值;NCEP1-COARE和 NCEP2-COARE分別代表應(yīng)用NCEP1和NCEP2各自數(shù)據(jù)庫(kù)中的海表面氣象參數(shù),采用 COARE3.0塊體通量計(jì)算程序獲得的感熱通量和潛熱通量,參考文獻(xiàn)[16]我們將其稱為“再計(jì)算感熱通量”和“再計(jì)算潛熱通量”,合稱為“再計(jì)算湍流熱通量”; 浮標(biāo)-COARE代表應(yīng)用浮標(biāo)觀測(cè)到的海表面氣象參數(shù),同樣采用 COARE3.0塊體通量計(jì)算程序獲得的湍流熱通量結(jié)果?？梢钥闯?NCEP再分析資料中的感熱通量和潛熱通量都存在很大程度上的高估現(xiàn)象,很多前人的研究都證實(shí)了這一點(diǎn)[9,23],這與模式計(jì)算中采用了過(guò)大的交換系數(shù)有關(guān)。當(dāng)我們采用塊體通量計(jì)算程序 COARE3.0來(lái)計(jì)算湍流熱通量時(shí),發(fā)現(xiàn)與 NCEP再分析資料中提供湍流熱通量相比較,再計(jì)算湍流熱通量大大地減小了,NCEP1再計(jì)算湍流熱通量與浮標(biāo)計(jì)算的湍流熱通量相比存在低估現(xiàn)象,NCEP2再計(jì)算湍流熱通量與浮標(biāo)計(jì)算的湍流熱通量相比發(fā)現(xiàn)在較低風(fēng)速下存在低估現(xiàn)象,而高風(fēng)速下由于NCEP2的風(fēng)速要遠(yuǎn)大于浮標(biāo)觀測(cè)風(fēng)速,造成此時(shí)NCEP2再計(jì)算湍流熱通量要大于浮標(biāo)計(jì)算的湍流熱通量。文獻(xiàn)[12,17]都指出再計(jì)算湍流熱通量存在低估現(xiàn)象。從表1可知,對(duì)于整個(gè)比較期間,NCEP1與浮標(biāo)Hs、HL之間的平均誤差分別為28.83和50.95 W/m2,平均高估Hs、HL分別為41.41%和32.19%,NCEP2與浮標(biāo)Hs、HL之間的平均誤差分別為44.69和90.96 W/m2,平均高估Hs、HL分別為64.19%和 57.47%,NCEP1 中Hs、HL的長(zhǎng)期平均值更接近于浮標(biāo)Hs、HL的長(zhǎng)期平均值。從表2可知,對(duì)于整個(gè)比較期間,NCEP1-COARE與浮標(biāo)Hs、HL之間的平均誤差分別為-23.01和-52.00W/m2,平均低估Hs、HL分別為-33.06%和-32.85%,NCEP2-COARE與浮標(biāo)Hs、HL之間的平均誤差分別為 9.21和7.39W/m2,平均高估Hs、HL分別為13.22%和4.67%,NCEP2-COARE 中Hs、HL的長(zhǎng)期平均值更接近于浮標(biāo)Hs、HL的長(zhǎng)期平均值。

圖3 湍流熱通量日平均值的時(shí)間序列Fig.3 Time series plots of daily averaged values of turbulent heat flux

圖4 凈熱通量日平均值的時(shí)間序列Fig.4 Time series plots of daily averaged values of net heat flux

2.2.3 凈熱通量

海表面的凈熱通量Fnh=Fus+Ful-Fds-Fdl+Hs+HL。圖 4給出了 5種海表面Fnh日平均值的時(shí)間序列。NCEP1和 NCEP2分別代表各自數(shù)據(jù)庫(kù)中直接下載的Fns、Fnl和湍流熱通量值計(jì)算的Fnh結(jié)果;NCEP1-COARE和 NCEP2-COARE分別代表應(yīng)用NCEP1和NCEP2各自數(shù)據(jù)庫(kù)中的Fns、Fnl和再計(jì)算湍流熱通量得到的再計(jì)算凈熱通量; 浮標(biāo)-COARE代表應(yīng)用浮標(biāo)觀測(cè)到的凈輻射通量和湍流熱通量計(jì)算的浮標(biāo)凈熱通量。可以看出,除了9月末10月初之外,NCEP1和NCEP2的凈熱通量結(jié)果要高于浮標(biāo)計(jì)算結(jié)果; 但用再計(jì)算湍流熱通量替代 NCEP再分析資料中的湍流熱通量之后,再計(jì)算凈熱通量大大地降低,NCEP1再計(jì)算凈熱通量與浮標(biāo)計(jì)算的凈熱通量相比除了 9月末兩者符號(hào)相反外其他時(shí)間存在低估現(xiàn)象,NCEP2再計(jì)算凈熱通量與浮標(biāo)計(jì)算的凈熱通量相比在10月9日左右、10月末和10月10日左右存在高估的現(xiàn)象,9月末兩者符號(hào)相反,其他時(shí)間點(diǎn)兩者比較接近。結(jié)合表1和表2中的Fnr、Hs、HL和Fnh,可知與浮標(biāo)觀測(cè)數(shù)據(jù)計(jì)算得到的海面凈熱通量相比,NCEP1凈熱通量(NCEP1再計(jì)算凈熱通量)高估(低估)了42%(40%); NCEP2凈熱通量(NCEP2再計(jì)算凈熱通量)高估(高估)了 68%(5%); 對(duì)于 NCEP1數(shù)據(jù)庫(kù)來(lái)講,海面凈熱通量的高估主要取決于Hs和HL的高估,海面再計(jì)算凈熱通量的低估主要取決于Hs和HL的低估; 對(duì)于NCEP2數(shù)據(jù)庫(kù)來(lái)講,海面凈熱通量的高估同樣主要取決于Hs和HL的高估,但海面再計(jì)算凈熱通量的高估主要取決于Hs、HL和Fnr的高估。

3 討論

在資料分析過(guò)程中發(fā)現(xiàn)NCEP2的再計(jì)算湍流熱通量與NCEP1相比要更接近于浮標(biāo)觀測(cè)值,但在11月1號(hào)和10號(hào)兩天也存在較大的誤差,試驗(yàn)觀測(cè)數(shù)據(jù)正好在兩個(gè)版本的中間。針對(duì)這一問(wèn)題我們采取了如下處理方式,即計(jì)算兩個(gè)版本的再計(jì)算湍流熱通量平均值,然后與浮標(biāo)計(jì)算得到的湍流熱通量進(jìn)行比較,兩者之間的平均偏差、相對(duì)誤差、均方根誤差和相關(guān)系數(shù) 4個(gè)統(tǒng)計(jì)量分別為-29.21W/m2,-12.81%,65.68 W/m2和 0.95,而 NCEP1(NCEP2)的再計(jì)算湍流熱通量與實(shí)際觀測(cè)到的湍流熱通量之間的 4個(gè)統(tǒng)計(jì)量分別為：-75.01 W/m2(16.60 W/m2)、-32.91%(7.28%)、102.36 W/m2(84.41 W/m2)和0.96(0.94)。從統(tǒng)計(jì)量上來(lái)講,NCEP1和NCEP2之間的平均再計(jì)算湍流熱通量要比NCEP1再計(jì)算湍流熱通量更接近于觀測(cè)值,但其平均偏差要大于 NCEP2再計(jì)算湍流熱通量與浮標(biāo)之間的平均偏差。圖 5為采用新的組合方式計(jì)算的湍流熱通量日平均值和浮標(biāo)觀測(cè)結(jié)果的時(shí)間序列比較圖,這里的混合風(fēng)場(chǎng)指的是在風(fēng)速U10＜11.5 m/s時(shí),風(fēng)速取為NCEP2的風(fēng)速資料,在風(fēng)速U10＞11.5 m/s時(shí),風(fēng)速取為 NCEP1的風(fēng)速資料。從圖5可以直觀地看出,對(duì)于11月1號(hào)和10號(hào)兩天來(lái)講,NCEP1、NCEP 2之間的再計(jì)算湍流熱通量平均值要明顯的好于 NCEP1或 NCEP2的再計(jì)算湍流熱通量。

另外,11月1號(hào)和10號(hào)兩天NCEP2的再計(jì)算湍流熱通量與浮標(biāo)觀測(cè)值誤差較大的原因在于兩者海表氣象參數(shù)之間的差異,尤其是海面上10 m風(fēng)速的差異。通過(guò)比較可以發(fā)現(xiàn),在高風(fēng)速下(U10＞11.5 m/s),NCEP1的風(fēng)速會(huì)更接近于浮標(biāo)測(cè)量風(fēng)速,則我們又考慮了另外一種處理方式,即在高風(fēng)速下,利用NCEP1的風(fēng)速代替NCEP2的風(fēng)速,重新計(jì)算NCEP2的再計(jì)算湍流熱通量,最后得到其與浮標(biāo)數(shù)據(jù)之間的平均偏差、相對(duì)誤差、均方根誤差和相關(guān)系數(shù) 4個(gè)統(tǒng)計(jì)量分別為8.17 W/m2,3.58%,65.39 W/m2和0.95,這一結(jié)果要優(yōu)于 NCEP1或NCEP2的再計(jì)算湍流熱通量,從圖 5中也可以看出使用混合風(fēng)場(chǎng)計(jì)算的再計(jì)算湍流熱通量非常地接近。所以說(shuō),針對(duì)本文的分析結(jié)果來(lái)講,這也許給我們以后應(yīng)用 NCEP數(shù)據(jù)庫(kù)的數(shù)據(jù)提供了一個(gè)新的思路：在使用塊體通量程序COARE3.0計(jì)算湍流熱通量的過(guò)程中,可以考慮使用混合風(fēng)場(chǎng),即高風(fēng)速下采用 NCEP1的風(fēng)速資料,其余情況下使用NCEP2的風(fēng)速資料,當(dāng)然這一想法還需要通過(guò)更多的實(shí)測(cè)資料來(lái)獲得驗(yàn)證。

圖5 采用新的組合方式計(jì)算的湍流熱通量日平均值的時(shí)間序列Fig.5 Time series plots of daily averaged values of turbulent heat flux by new recomposed way

4 結(jié)論

本文對(duì) NCEP再分析資料中的海表氣象學(xué)參數(shù)及海表熱通量日平均值與定點(diǎn)布放在黃海北部的浮標(biāo)觀測(cè)資料及通量計(jì)算結(jié)果進(jìn)行了比較分析。比較時(shí)間從2009年9月26日至2009年11月20日。統(tǒng)計(jì)量平均偏差、相對(duì)誤差、均方根誤差和相關(guān)系數(shù)比對(duì)表明：NCEP1和NCEP2再分析資料中海表氣象學(xué)參數(shù)及海表熱通量日平均值的時(shí)間序列與浮標(biāo)結(jié)果變化趨勢(shì)基本一致,但NCEP1對(duì)風(fēng)速存在低估現(xiàn)象,對(duì)大氣比濕、大氣溫度和海表溫度存在高估現(xiàn)象,NCEP2中海表氣象學(xué)參數(shù)都存在高估現(xiàn)象但更接近浮標(biāo)觀測(cè)資料; 雖然在向上短波輻射通量上 NCEP2比NCEP1更接近浮標(biāo)觀測(cè)值,但NCEP1的凈輻射通量值更接近浮標(biāo)觀測(cè)值。對(duì)湍流熱通量而言,與NCEP1和 NCEP2再分析資料中自帶的感熱通量和潛熱通量及NCEP1再計(jì)算湍流熱通量相比,NCEP2再計(jì)算湍流熱通量最接近于浮標(biāo)資料計(jì)算的湍流熱通量。對(duì)于海表凈熱通量而言,應(yīng)用 NCEP1和NCEP2資料提供的湍流熱通量計(jì)算的結(jié)果要普遍高于浮標(biāo)計(jì)算結(jié)果; 但用再計(jì)算湍流熱通量替代NCEP再分析資料中的湍流熱通量之后,再計(jì)算凈熱通量大大地降低,NCEP1的再計(jì)算凈熱通量低估了42%,而NCEP2的再計(jì)算凈熱通量高估了5%; 對(duì)于NCEP1數(shù)據(jù)庫(kù)來(lái)講,海面凈熱通量的高估主要取決于感熱和潛熱通量的高估,海面再計(jì)算凈熱通量的低估主要取決于感熱和潛熱通量的低估; 對(duì)于NCEP2數(shù)據(jù)庫(kù)來(lái)講,海面凈熱通量的高估同樣主要取決于感熱和潛熱通量的高估,但海面再計(jì)算凈熱通量的高估主要取決于感熱、潛熱通量和凈輻射通量的高估。

海面風(fēng)場(chǎng)的比對(duì)表明,在風(fēng)速小于 11 m/s時(shí),NCEP1和 NCEP2風(fēng)速資料一般要低于浮標(biāo)觀測(cè)風(fēng)速,NCEP2與浮標(biāo)觀測(cè)的風(fēng)速更接近; 當(dāng)風(fēng)速大于11 m/s時(shí),NCEP1和NCEP2一般要高于浮標(biāo)觀測(cè)風(fēng)速,NCEP1與浮標(biāo)觀測(cè)的風(fēng)速更接近。針對(duì)該分析結(jié)果來(lái)講,這也許給我們以后應(yīng)用 NCEP數(shù)據(jù)庫(kù)的數(shù)據(jù)提供了一個(gè)新的思路：在使用塊體通量程序COARE3.0計(jì)算湍流熱通量的過(guò)程中,可以考慮使用混合風(fēng)場(chǎng),即高風(fēng)速下采用 NCEP1的風(fēng)速資料,其余情況下使用NCEP2的風(fēng)速資料,計(jì)算海面上的湍流熱通量。當(dāng)然這一想法還需要通過(guò)更多的實(shí)測(cè)資料來(lái)獲得驗(yàn)證。

研究表明,NCEP再分析資料中的海表氣象參數(shù)是可信的,在統(tǒng)計(jì)意義上,NCEP2給出的海表氣象參數(shù)比NCEP1與浮標(biāo)觀測(cè)值更接近,而凈輻射通量則是NCEP1更接近于浮標(biāo)觀測(cè)值; 使用更合理的塊體公式計(jì)算出的湍流熱通量來(lái)代替 NCEP數(shù)據(jù)庫(kù)中已計(jì)算好的湍流熱通量來(lái)驅(qū)動(dòng)海洋模式是很有必要的; NCEP2的再計(jì)算湍流熱通量與NCEP1相比要更接近于浮標(biāo)觀測(cè)值。使用混合風(fēng)場(chǎng)計(jì)算的再計(jì)算湍流熱通量與浮標(biāo)觀測(cè)結(jié)果更加符合。

[1]Kalnay E,Kanamitsu M,Kistler R,et al.The ncep/ncar 40-year reanalysis project [J].Bull Amer Meteor Soc,1996,77：437-471.

[2]Kanamitsu M,Ebisuzaki W,Woolen J,et al.Overview of ncep/doe reanalysis-2 [EB/OL].[2011-01-02].http：//dss.ucar.edu/datasets/ds091.0/docs/ncep-docs/r2sumshort.pdf.

[3]Kanamitsu M,Ebisuzaki W,Woolen J,et al.The ncep-doe amip-Ⅱ reanalysis (r-2) [J].Bull Amer Meteor Soc,2002,83：1631-1643.

[4]Weller R A,Baumgartner M F,Josey S A,et al.Atmospheric forcing in the Arabian Sea during 1994-1995：observations and comparisons with climatology and models [J].Deep-Sea research Part II-Topical Studies in Oceanography,1998,45B：1961-1999.

[5]Weller R A,Anderson S P.Surface meteorology and air-sea fluxes in the western equatorial Pacific warm pool during the toga coupled ocean-atmosphere response experiment [J].Journal of Climate,1996,9：1959-1999.

[6]Zeng Xubin,Zhao Ming,Dickinson R E.Intercomparison of bulk aerodynamic algorithms for the computation of sea surface fluxes using toga coare and tao data [J].Journal of Climate,1998,11：2628-2644.

[7]Wang W,McPhaden M J.What is the mean seasonal cycle of surface heat flux in the equatorial Pacific? [J].Journal of Geophysical Research-Oceans,2001,106：837-857.

[8]Moyer K A,Weller R A.Observations of surface forcing from the subduction experiment：a comparison with global model products and climatological datasets [J].Journal of Climate,1997,10：2725-2742.

[9]Josey S A.A comparison of ecmwf,ncep-ncar,and soc surface heat fluxes with moored buoy measurements in the subduction region of the Northeast Atlantic [J].Journal of Climate,2001,14：1780-1789.

[10]Renfrew I A,Moore G W K,Guest P S,et al.A comparison of surface layer and surface turbulent flux observations over the Labrador Sea with ecmwf analyses and ncep reanalyses [J].Journal of Physical Oceanography,2002,32：383-400.

[11]Moore G W K,Renfrew I A.An assessment of the surface turbulent heat fluxes from the ncep-ncar reanalysis over the western boundary currents [J].Journal of Climate 2002,15：2020-2037.

[12]Rouault M,Reason C J C,Lutjeharms J R E,et al.Underestimation of latent and sensible heat fluxes above the Agulhas Current in ncep and ecmwf analyses[J].Journal of Climate,2003,16：776-782.

[13]Goswami B N,Sengupta D.A note on the deficiency of ncep/ncar reanalysis surface winds over the equatorial Indian Ocean [J].Journal of Geophysical Research-Oceans,2003,108(C4)：3124.doi：10.1029/2002JC001497.

[14]Goswami B N,Rajagopal E N.Indian Ocean surface winds from ncmrwf analysis as compared to quikscat and moored buoy winds [J].Proceedings of the Indian Academy of Sciences-Earthand Planetary Sciences,2003,112：61-77.

[15]張?jiān)龊?姜華,王輝,等.南海西南季風(fēng)期NCEP2湍流熱通量的質(zhì)量分析[J].海洋學(xué)報(bào),2006,2：29-38.

[16]史劍,趙棟梁,孫建.利用浮標(biāo)資料初步評(píng)估 NCEP再分析湍通量[J].海洋預(yù)報(bào),2007,2：74-82.

[17]Swain D.Comparison of ncep turbulent heat fluxes with in situ observations over the south-eastern Arabian Sea [J].Meteorology and Atmospheric Physics,2009,104：163-175.

[18]Curry J A,Clayson C A,Rossow W B,et al.High-resolution satellite-derived dataset of the surface fluxes of heat,freshwater,and momentum for the toga coare iop [J].Bulletin of The American Meteorological Society,1999,80：2059-2080.

[19]Halpern D,Freilich M H,Weller R A.Ecmwf and ers-1 surface winds over the Arabian Sea during July 1995 [J].Journal of Physical Oceanography,1999,29：1619-1623.

[20]黃艷松,宋金寶,范聰慧.海氣通量渦相關(guān)法計(jì)算中的時(shí)間尺度分析[J].海洋科學(xué),2011,35(11)：114-119.

[21]Panofsky H A,Dutton J A.Atmospheric turbulence [M].New York：Wiley,1984：397.

[22]Fairall C W,Bradley E F,Rogers D P,et al.Bulk parameterization of air-sea fluxes for tropical ocean global atmosphere coupled ocean atmosphere response experiment [J].Journal of Geophysical Research-Oceans,1996,101：3747-3764.

[23]Sun B,Yu L,Weller R A.Comparisons of surface meteorology and turbulent heat fluxes over the Atlantic：nwp model analyses versus moored buoy observations[J].J Clim,2003,16：679-695.

[24]Donlon C W,Minnet P J,Gentemann C,et al.Towards improved validation of satellite sea surface skin temperature measurements for climate research [J].J Clim,2002,15：353-369.