全球海氣界面潛熱通量產(chǎn)品對(duì)比研究*

楊新飛， 管　磊

(中國海洋大學(xué)信息科學(xué)與工程學(xué)院, 山東 青島 266100)

?

研究簡報(bào)

全球海氣界面潛熱通量產(chǎn)品對(duì)比研究*

楊新飛， 管磊**

(中國海洋大學(xué)信息科學(xué)與工程學(xué)院, 山東 青島 266100)

摘要：潛熱是大氣能量的重要來源，海氣界面潛熱通量對(duì)于海氣相互作用的研究具有重要意義。本文比較了1996—2005年J-OFURO2、HOAPS3、GSSTF3、OAFLUX、NCEP/NCAR再分析資料5種產(chǎn)品，研究結(jié)果表明5種月平均潛熱通量產(chǎn)品基本一致，其中J-OFURO2、HOAPS3、GSSTF3較為接近，OAFLUX與NCEP/NCAR較為接近。比較不同的月平均潛熱通量產(chǎn)品的統(tǒng)計(jì)特征，均方根誤差最小的是GSSTF3與HOAPS3(14.08 W/m2)，最大的是GSSTF3與NCEP/NCAR(25.98 W/m2)，相關(guān)系數(shù)最小的是HOAPS3與NCEP/NCAR(0.87)、GSSTF3與OAFLUX(0.87)、最大的是GSSTF3與HOAPS3(0.97)，5種潛熱通量產(chǎn)品多年平均的緯向分布雖然基本相似但仍然存在差異，基于遙感數(shù)據(jù)的3種產(chǎn)品與OAFLUX、NCEP/NCAR再分析資料之間較大差值主要分布在太平洋和大西洋赤道區(qū)域、西北太平洋、印度洋北部。

關(guān)鍵詞：潛熱通量；熱通量；海氣界面；衛(wèi)星遙感

潛熱是大氣能量的重要來源，海氣界面熱通量是海氣相互作用的關(guān)鍵因子，包括太陽短波輻射、海面長波輻射，潛熱通量和感熱通量，其中短波輻射將熱量由大氣傳遞到海洋，另外3種主要將熱量由海洋傳遞到大氣[1]。由海洋傳遞到大氣相當(dāng)部分的熱量是以潛熱的形式進(jìn)行的[2]，因此潛熱對(duì)于研究海-氣相互作用十分必要。

潛熱通量一般采用塊體參數(shù)法計(jì)算，其計(jì)算公式如下：

Hl=ρa(bǔ)LeCeU(qs-qa)。

其中：ρa(bǔ)為空氣密度；Le為蒸發(fā)潛熱；Ce為潛熱交換系數(shù)。U為海表面風(fēng)速；qs為海表面飽和比濕；qa為近海面空氣比濕。

目前海氣潛熱通量產(chǎn)品主要有4種類型：基于實(shí)測數(shù)據(jù)(實(shí)驗(yàn)船，浮標(biāo)和觀測站)的產(chǎn)品、基于衛(wèi)星遙感數(shù)據(jù)的產(chǎn)品、基于數(shù)值模型的產(chǎn)品和基于合成數(shù)據(jù)(再分析資料和衛(wèi)星遙感數(shù)據(jù)的融合數(shù)據(jù))的產(chǎn)品。在過去的十幾年中，一些科研機(jī)構(gòu)用不同的算法獲得潛熱通量產(chǎn)品，其中具有代表性的包括：J-OFURO2(the Japanese Ocean Fluxes Data Sets with Use of Remote Sensing Observations version 2)、HOAPS3(the Hamburg Ocean Atmosphere Parameters and Fluxes from Satellite Data version 3)、GSSTF3(the Goddard Satellite-Based Surface Turbulent Fluxes data set version 3)、OAFLUX(the Objectively Analyzed Air Sea Fluxes for the Global Oceans data set)、ECMWF再分析數(shù)據(jù)(The European Centre for Medium-Range Weather Forecasts)、NCEP/NCAR再分析數(shù)據(jù)(the National Centers for Environmental Prediction and the National Center for Atmospheric Research Reanalysis dataset)等。這些產(chǎn)品用于研究海氣相互作用、海洋環(huán)流模型等。然而不同的潛熱通量產(chǎn)品之間存在一定的差異，因此近年來，很多學(xué)者對(duì)此進(jìn)行了研究。

Kubota等人[1]對(duì)比了J-OFURO、NCEP/NCAR、ECMWF和da Silva全球潛熱通量月平均產(chǎn)品，由于da Silva只有1993年之前的數(shù)據(jù)，因此da Silva參與對(duì)比的時(shí)間區(qū)間為1992—1993年，其他產(chǎn)品為1992—1994年。結(jié)果顯示就平均值而言，5種產(chǎn)品具有相似的模式，其中J-OFURO與GSSTF較為接近，兩者在熱帶海域的潛熱值都高于HOAPS，J-OFURO與da Silva差異較大，由于南半球海域觀測值稀少，造成J-OFURO與NCEP/NCAR、ECMWF、da Silva三個(gè)產(chǎn)品的時(shí)間相關(guān)性較差。

Chou等人[3]將GSSTF2與GSSTF1、HOAPS、NCEP/NCAR、da Silva 1992—1993年間緯向平均潛熱通量進(jìn)行了對(duì)比，結(jié)果顯示GSSTF2潛熱通量較其他4種產(chǎn)品更為可靠。

Bentamy等人[4]將1996年10月—1997年6月9個(gè)月的基于衛(wèi)星遙感數(shù)據(jù)的周平均潛熱通量分別與船測、浮標(biāo)數(shù)據(jù)、再分析數(shù)據(jù)產(chǎn)品進(jìn)行比較，結(jié)果顯示與船測數(shù)據(jù)的均方根誤差達(dá)34 W/m2，與浮標(biāo)數(shù)據(jù)的均方根誤差低于30 W/m2。再分析數(shù)據(jù)產(chǎn)品NCEP/NCAR、ECMWF的潛熱比衛(wèi)星遙感數(shù)據(jù)分別高出10%和20%。

Brunke等人[5]對(duì)比了11種海面通量產(chǎn)品的潛熱通量和感熱通量，包括6種再分析產(chǎn)品(MERRA、CFSR、ERA-Interim、NCEP/NCAR、NCEP/DOE、ERA-40)，4種衛(wèi)星產(chǎn)品(GSSTF2、GSSTF2b、J-OFURO2、HOAPS2)和1種綜合產(chǎn)品(OAFLUX)。作者根據(jù)與實(shí)測數(shù)據(jù)的總偏差和標(biāo)準(zhǔn)偏差由小到大的次序?qū)?1個(gè)產(chǎn)品分為A、B、C 3個(gè)等級(jí)，其中潛熱通量的A等級(jí)產(chǎn)品有：ERA-40、MERRA、ERA-Interim、GSSTF2b；B等級(jí)的產(chǎn)品有：CFSR、JOFURO、NCEP/NCAR、OAFLUX；C等級(jí)產(chǎn)品有：HOAPS2、NCEP/DOE、GSSTF2。

本文對(duì)比分析了5種全球海氣界面潛熱通量產(chǎn)品，分別是J-OFURO2、HOAPS3、GSSTF3、OAFLUX和NCEP/NCAR再分析數(shù)據(jù)，前3種潛熱通量產(chǎn)品基于衛(wèi)星遙感數(shù)據(jù)，OAFLUX基于合成數(shù)據(jù)，NCEP/NCAR再分析資料基于氣象數(shù)值模型。

1數(shù)據(jù)簡介

本文選用1996—2005年的全球海氣界面的月平均潛熱通量作為研究對(duì)象，本節(jié)簡要介紹5種潛熱通量產(chǎn)品。

J-OFURO2是由Tokai University基于衛(wèi)星遙感數(shù)據(jù)計(jì)算得到全球海氣相互作用參數(shù)的產(chǎn)品[6]，其潛熱通量根據(jù)COARE3.0塊體參數(shù)法[7]推算，相較于J-OFURO的單一數(shù)據(jù)源，J-OFURO2采用了多種衛(wèi)星傳感器的融合數(shù)據(jù)，提高了空間覆蓋率的同時(shí)也減小了插值造成的誤差。風(fēng)速主要來源于微波輻射計(jì)(DMSP/SSMI、Aqua/AMSRE、TRMM/TMI)和微波散射計(jì)(ERS/AMI、QuickScat/SeaWinds)的融合數(shù)據(jù)。近海面空氣比濕來源于SSMI[8]，海面比濕根據(jù)日本氣象廳提供的衛(wèi)星融合海表溫度數(shù)據(jù)(Merged satellite and in situ data Global Daily Sea Surface Temperatures，MGDSST)推算得到。J-OFURO2的空間分辨率為1°×1°，覆蓋范圍為：90°S～90°N，180°W～180°E，時(shí)間跨度為1987年1月—2006年12月，數(shù)據(jù)下載地址為：http://dtsv.scc.u-tokai.ac.jp/j-ofuro/data_retrieval_and_download.html。

HOAPS3是Max Planck Institute for Meteorology和the University of Hamburg合作的基于衛(wèi)星遙感數(shù)據(jù)計(jì)算得到全球海氣相互作用參數(shù)產(chǎn)品，潛熱通量采用COARE3.0塊體參數(shù)法推算[9-10]，相較于HOAPS2，時(shí)間跨度由1987—2002年擴(kuò)展為1987—2005年。風(fēng)速、近海面空氣比濕來源于SSMI，海表溫度來源于改進(jìn)型甚高分辨率輻射計(jì)探路者(AVHRR Pathfinder)[11]，其中風(fēng)速采用人工神經(jīng)網(wǎng)絡(luò)反演得到[12]，空氣比濕根據(jù)與SSMI亮溫的線性關(guān)系計(jì)算得到[4]，海表比濕根據(jù)馬格斯公式由海表溫度推算得到[13]。HOAPS3的空間分辨率是0.5°×0.5°，覆蓋范圍為：80°S～80°N，180°W～180°E，時(shí)間跨度為1987年7月—2005年12月，數(shù)據(jù)下載地址為：http://www.hoaps.org/index.php?id=data_access。

GSSTF3是National Aeronautics and Space Administration/Goddard Earth Science Data Information and Services Center基于衛(wèi)星遙感數(shù)據(jù)資料計(jì)算得到全球海氣相互作用參數(shù)產(chǎn)品，潛熱通量采用Chou等人的塊體參數(shù)法計(jì)算得到[3,14-15]，相較于之前的版本GSSTF2、GSSTF2c，GSSTF3改進(jìn)了海表比濕方法的反演算法[15-16]。風(fēng)速、近海面空氣比濕都來源于SSMI[17-18]，海表溫度來源于NCEP/DOE再分析資料。GSSTF3的空間分辨率是0.25°×0.25°。覆蓋范圍為：90°S～90°N，180°W～180°E，時(shí)間跨度為：1987年7月—2008年12月，數(shù)據(jù)下載地址為：ftp://measures.gsfc.nasa.gov/data/s4pa/GSSTF/GSSTF.3/。

OAFLUX是Woods Hole Oceanographic Institution基于融合數(shù)據(jù)計(jì)算得到的全球海氣相互作用參數(shù)產(chǎn)品，潛熱通量采用COARE3.0塊體算法計(jì)算得到[19]。風(fēng)速來源于SSMI、AMSRE、QuickScat、NCEP/NCAR、NCEP/DOE和ERA40的融合數(shù)據(jù)；海表空氣比濕來源于SSMI、NCEP/NCAR、NCEP/DOE和ERA40融合數(shù)據(jù)。OAFLUX的空間分辨率是，覆蓋范圍為：90°S～90°N，180°W～180°E。時(shí)間跨度為：1954年1月至今，數(shù)據(jù)下載地址為http://oaflux.whoi.edu/。

NCEP/NCAR再分析資料是National Center for Environmental Prediction和National Center for Atmospheric Research自1991年開始聯(lián)合進(jìn)行的大氣資料再分析研究項(xiàng)目產(chǎn)品。海表溫度來源于大氣模型比較項(xiàng)目(AMIP-Ⅱ)的分析數(shù)據(jù)[20]。NCEP/NCAR再分析數(shù)據(jù)采用T62高斯網(wǎng)格，空間分辨率約為1.9°×1.875°。覆蓋范圍為：90°S～90°N，180°W～180°E，時(shí)間跨度為1948年1月至今，數(shù)據(jù)下載地址為：ftp://ftp.cdc.noaa.gov/pub/Datasets/ncep.reanalysis/surface_gauss。

2數(shù)據(jù)對(duì)比

由于5種潛熱通量產(chǎn)品空間分辨率不同，為了便于潛熱通量產(chǎn)品對(duì)比，本文首先對(duì)5種產(chǎn)品進(jìn)行了網(wǎng)格統(tǒng)一，將HOAPS3、GSSTF3和NCEP/NCAR插值為1°×1°的網(wǎng)格數(shù)據(jù)，然后進(jìn)行交叉比較分析。

2.1 產(chǎn)品間月平均潛熱通量的偏差、標(biāo)準(zhǔn)偏差、均方根誤差和相關(guān)系數(shù)比較

本文將5種潛熱通量的月平均產(chǎn)品做交叉比較，產(chǎn)品間偏差和標(biāo)準(zhǔn)偏差隨時(shí)間的變化曲線如圖1。圖中a～j分別為J-H、J-G、J-O、J-N、H-G、H-O、H-N、G-O、G-N、N-O(J代表J-OFURO2，H代表HOAPS3，G代表GSSTF3，O代表OAFLUX，N代表NCEP/NCAR，“J-H”代表用J-OFURO2潛熱值減去HOAPS3潛熱值)的偏差和標(biāo)準(zhǔn)偏差，圖中黑線代表偏差，紅色豎線代表標(biāo)準(zhǔn)偏差。表1統(tǒng)計(jì)了5種潛熱通量產(chǎn)品1996—2005年的均方根誤差和相關(guān)系數(shù)。

研究結(jié)果顯示J-OFURO2的月平均潛熱低于HOAPS3和NCEP/NCAR，偏差范圍分別為-9.41～1.57 W/m2、-10.95～4.60 W/m2；而與GSSTF3的偏差呈現(xiàn)周期性，偏差范圍為-3.86～6.26W/m2，夏季J-OFURO2的潛熱低于GSSTF3，春秋季節(jié)則高于GSSTF3。1996—2000年J-OFURO2在夏季時(shí)潛熱高于OAFLUX，2000年之后J-OFURO2與OAFLUX偏差變化幅度減小，J-OFURO2低于OAFLUX，兩者1996—2005年的偏差范圍為-5.13～10.92 W/m2。J-OFURO2與HOAPS3的標(biāo)準(zhǔn)偏差最小，其次是與GSSTF3。

HOAPS3的月平均潛熱值高于J-OFURO2、GSSTF3、OAFLUX，偏差范圍分別為：-1.57～9.41 W/m2、1.58～9.47 W/m2、-1.14～11.73 W/m2。HOAPS3與NCEP/NCAR的偏差范圍為-7.94～6.77 W/m2。HOAPS3與GSSTF3、J-OFURO2的標(biāo)準(zhǔn)偏差較小。

GSSTF3與OAFLUX、NCEP/NCAR的偏差范圍分別為：-6.58～7.17 W/m2、-0.42～14.0 2W/m2，大于GSSTF3與HOAPS3、J-OFURO2的偏差范圍，GSSTF3與HOAPS3的標(biāo)準(zhǔn)偏差最小。

OAFLUX的月平均潛熱值低于HOAPS3和NCEP/NCAR，偏差范圍分別為：-11.73～1.14 W/m2、-10.43～2.45 W/m2；與J-OFURO2、GSSTF3潛熱值的正負(fù)偏差分布均衡，偏差范圍分別為：-10.92～5.13 W/m2，-7.17～6.58 W/m2。OAFLUX與NCEP/NCAR的標(biāo)準(zhǔn)偏差最小，OAFLUX與3種衛(wèi)星產(chǎn)品之間的標(biāo)準(zhǔn)偏差相差不大。

(a. J-H b. J-G c. J-O d. J-N e. H-G f. H-O g. H-N h. G-O i. G-N j. N-O)

J-HJ-GJ-OJ-NH-GH-OH-NG-OG-NN-O均方根誤差①(RMS)/W·m-215.7616.8620.5725.2114.0822.4025.7024.1325.9818.67相關(guān)系數(shù)②(R)0.950.940.910.880.970.900.870.870.880.95

Note：①Root mean square error;②Correlation

NCEP/NCAR的月平均潛熱值高于其他數(shù)據(jù)，與OAFLUX、J-OFURO2、HOAPS3、GSSTF3的月平均偏差范圍為：2.45～10.43 W/m2、-4.60～10.95 W/m2、-6.77～7.94 W/m2、-0.42～14.02 W/m2。

由表1可知：基于衛(wèi)星遙感數(shù)據(jù)的3種產(chǎn)品J-OFURO2、HOAPS3和GSSTF3之間相關(guān)系數(shù)均高于0.9，均方根誤差均小于20 W/m2；NCEP/NCAR再分析資料產(chǎn)品與3種基于衛(wèi)星遙感數(shù)據(jù)的產(chǎn)品相關(guān)系數(shù)均低于0.9，均方根誤差均大于20 W/m2；基于融合數(shù)據(jù)的OAFLUX與NCEP/NCAR再分析資料之間的相關(guān)系數(shù)大于其與3種衛(wèi)星數(shù)據(jù)產(chǎn)品的相關(guān)系數(shù)，相應(yīng)地OAFLUX與NCEP/NCAR之間的均方根誤差值小于其與3種衛(wèi)星數(shù)據(jù)產(chǎn)品的。相關(guān)系數(shù)最大值存在于HOAPS3和GSSTF3之間(0.97)，最小值存在于GSSTF3和OAFLUX之間(0.87)、HOAPS3和NCEP/NCAR之間(0.87)。均方根誤差最大值存在于NCEP/NCAR和GSSTF3之間(25.98 W/m2)，最小值存在于HOAPS3和GSSTF3之間(14.08 W/m2)。

2.2 產(chǎn)品間潛熱通量值緯向分布比較

將5種潛熱通量的月平均數(shù)據(jù)沿緯度方向進(jìn)行平均獲得全球海氣界面潛熱通量緯向變化趨勢，結(jié)果如圖2。圖2表明：5種產(chǎn)品的全球緯向分布具有相同模式，潛熱極大值出現(xiàn)在南北信風(fēng)帶，赤道的潛熱顯著低于信風(fēng)帶海域；在南半球中高緯度海域和北半球中緯度海域，潛熱通量隨著緯度的升高降低，而在60°N附近潛熱降低幅度變慢，GSSTF3和J-OFURO2有小幅的回升。

5種產(chǎn)品在0°～10°S、30°S、40°S、40°N附近海域潛熱值差異較小，極大值附近、30°N～40°N和高緯度海域差異較大，其中在極大值區(qū)域OAFLUX的潛熱值與其他4種產(chǎn)品差異明顯，潛熱通量最低，NCEP/NCAR潛熱通量最高。在南北半球高緯度海域和40°S～60°S之間海域NCEP/NCAR和OAFLUX均低于其他數(shù)據(jù)。3種基于衛(wèi)星遙感數(shù)據(jù)的產(chǎn)品整體上接近，差異主要分布在南半球中高緯度和北半球高緯度海域：在50°S～60°S海域J-OFURO2和GSSTF3高于HOAPS3，在50°N～60°N海域GSSTF3和HOAPS3高于J-OFURO2。NCEP/NCAR和OAFLUX在潛熱極大值分布區(qū)以外的海域很接近。

圖2　1996—2005年平均潛熱通量緯向分布

2.3 產(chǎn)品間潛熱通量平均差值的全球分布

將5種產(chǎn)品月平均數(shù)據(jù)兩兩比較，獲得產(chǎn)品間1996—2005年潛熱通量平均差值的全球分布，結(jié)果顯示3種基于衛(wèi)星遙感數(shù)據(jù)的產(chǎn)品之間、NCEP/NCAR與OAFLUX之間平均差值在全球范圍內(nèi)較小，3種基于衛(wèi)星遙感數(shù)據(jù)的產(chǎn)品與NCEP/NCAR、OAFLUX之間平均差值較大。圖3給出了J-OFURO2、HOAPS3、GSSTF3分別于與NCEP/NCAR、OAFLUX的全球潛熱通量平均差值分布，圖中a～f分別為J-N、J-O、H-N、H-O、G-N、G-O，符號(hào)意義同圖1。

由圖3顯示：J-OFURO2，HOAPS3，GSSTF3與NCEP/NCAR、OAFLUX相比，在南半球太平洋赤道地區(qū)前者明顯高于后者，而在西北太平洋、印度洋北部、西北大西洋以及大西洋赤道海域則相反。在50°S～60°S海域J-OFURO2、HOAPS3、GSSTF3與NCEP/NCAR、OAFLUX的潛熱差值為正。另外，J-OFURO2、HOAPS3、GSSTF3與OAFLUX的潛熱差值在北太平洋還有一個(gè)較大范圍的高值區(qū)，如圖3中右側(cè)3幅圖所示，這個(gè)區(qū)域內(nèi)J-OFURO2與OAFLUX的差值明顯低于HOAPS3與OAFLUX，GSSTF3與OAFLUX的差值。

Brunke等人根據(jù)潛熱通量產(chǎn)品評(píng)估結(jié)果將GSSTF2b劃分為A類，J-OFURO、NCEP/NCAR、OAFLUX劃分為B類，HOAPS劃分為C類，本文針對(duì)改進(jìn)后發(fā)布的潛熱通量產(chǎn)品J-OFURO2、HOAPS3、GSSTF3等進(jìn)行了比較，比較結(jié)果可見J-OFURO2、HOAPS3、GSSTF3較為接近，NCEP/NCAR、OAFLUX較為接近。5種產(chǎn)品的差異部分原因是由于計(jì)算潛熱通量過程中輸入數(shù)據(jù)的差異造成的，對(duì)于空氣比濕，3種衛(wèi)星數(shù)據(jù)產(chǎn)品采用SSMI的反演數(shù)據(jù)[4,8,18]，OAFLUX采用SSMI反演數(shù)據(jù)[16]與NCEP/NCAR、NCEP/DOE、ERA-40再分析數(shù)據(jù)的融合數(shù)據(jù)。4種產(chǎn)品的空氣比濕來源中均有SSMI衛(wèi)星數(shù)據(jù)，但是采用的反演算法不盡相同，其中GSSTF3和HOAPS3采用Bentamy等人的算法[4,18]，J-OFURO2采用Schlüssel等人的算法[8]，OAFLUX采用Chou等人的算法[16]。Bentamy等人和Schlüssel等人用SSMI的亮溫值一步計(jì)算獲得空氣比濕，而Chou等人首先用SSMI亮溫值計(jì)算得到大氣邊界層底部500 m的水汽含量WB，然后對(duì)WB進(jìn)行經(jīng)驗(yàn)正交分解(EOF)獲得空氣比濕，比濕數(shù)據(jù)的比較顯示J-OFURO2、HOAPS3、GSSTF3三種產(chǎn)品的空氣比濕較為接近，NCEP/NCAR與OAFLUX的空氣比濕較為接近。

3結(jié)論

本文對(duì)J-OFURO2、HOAPS3、GSSTF3、OAFLUX和NCEP/NCAR再分析資料的全球海氣界面月平均潛熱通量進(jìn)行了交叉比較分析，由于5種產(chǎn)品輸入?yún)?shù)來源、獲取潛熱通量方法不同，其產(chǎn)品存在差異，主要結(jié)果如下：

圖3　1996—2005年潛熱通量平均差值分布

(1)J-OFURO2、GSSTF3和HOAPS3之間偏差和標(biāo)準(zhǔn)偏差均較小，說明3種基于衛(wèi)星遙感數(shù)據(jù)的潛熱產(chǎn)品較為相似，NCEP/NCAR和OAFLUX之間的偏差和標(biāo)準(zhǔn)偏差較NCEP/NCAR與其他3種數(shù)據(jù)的小，表明NCEP/NCAR與OAFLUX較為接近。

(2)5種潛熱通量產(chǎn)品在緯向分布上具有相同的模式，潛熱最大值出現(xiàn)在南北信風(fēng)帶海域，赤道海域潛熱顯著低于此海域的潛熱；在南半球中高緯度海域和北半球中緯度海域，潛熱通量隨著緯度的升高降低。潛熱較小差異值分布在0°～10°S、30°S、40°S、40°N附近海域，差異較大值分布在南北半球信風(fēng)帶和高緯度海域。

(3)J-OFURO2、HOAPS3、GSSTF3的潛熱通量比較接近，這3種數(shù)據(jù)與OAFLUX的平均潛熱差值在南半球中高緯度海域、信風(fēng)區(qū)海域表現(xiàn)為正偏差，負(fù)偏差主要分布在赤道、南北半球30°附近海域。在太平洋中部有兩個(gè)較大的正偏差分布區(qū)域。J-OFURO2、HOAPS3、GSSTF3與NCEP/NCAR的全球平均潛熱差值分布和三者與OAFLUX的分布相似，但是負(fù)偏差范圍更廣，且在太平洋中心只有一個(gè)明顯的正偏差區(qū)域。

致謝:本文J-OFURO2數(shù)據(jù)由School of Marine Science and Technology, Tokai University提供；HOAPS3 數(shù)據(jù)由Max Planck Institute for Meteorology和 The University of Hamburg提供；GSSTF3數(shù)據(jù)由National Aeronautics and Space Administration/Goddard Earth Science Data Information and Services Center提供；OAFLUX數(shù)據(jù)由；Woods Hole Oceanographic Institution提供，NCEP/NCAR再分析資料由National Center for Environmental Prediction和National Center for Atmospheric Research提供。

參考文獻(xiàn)：

[1]Kubota M, Kano A, Muramatsu H, et al. Intercomparison of various surface latent heat flux Fields[J]. Journal of Climate, 2003, 16: 670-678.

[2]Grist J P, Josey S A. Inverse analysis adjustment of the SOC air-Sea flux climatology using ocean heat transport constraints[J]. Journal of Climate, 2003, 16(20): 3274-3295.

[3]Chou S H, Nelkin E, Ardizzone J, et al. Surface turbulent heat and momentum fluxes over global oceans based on the goddard satellite retrievals, Version 2 (GSSTF2)[J]. Journal of Climate, 2003, 16: 3256-3273.

[4]Bentamy A, Katsaros K B, Mestas A M, et al. Satellite estimates of wind speed and latent heat flux over the global oceans[J]. Journal of climate, 2003, 16(4):637-656.

[5]Brunke M A, Wang Z, Zeng X, et al. An assessment of the uncertainties in ocean surface turbulent fluxes in 11 reanalysis, satellite-derived, and combined global datasets[J]. Journal of Climate, 2011, 24(21): 5469-5492.

[6]Kubota M, Tomita H. The present state of the J-OFURO air-sea interaction data product[J]. Flux News, 2007, 4: 13-15.

[7]Fairall C W, Bradley E F, Hare J E, et al. Bulk param-eterization of air-sea fluxes: updates and verification for the COARE algorithm.[J]. Journal of Climate, 2003, 16: 571-591.

[8]Schlüssel P, Schanz L, Englisch G. Retrieval of latent heat flux and longwave irradiance at the sea surface from SSM/I and AVHRR measurements[J]. Advances in Space Research, 1995, 16(10): 107-116.

[9]Andersson A, Fennig K, Klepp C, et al. The Hamburg ocean atmosphere parameters and fluxes from satellite Data-HOAPS-3[J]. Earth System Science Data, 2010, 2(2): 215-234.

[10]Andersson A, Bakan S, Fennig K, et al. Hamburg ocean atmosphere parameters and fluxes from satellite data HOAPS-3 monthly mean, electronic publication[J]. World Data Center for Climate, 2007: 3097,doi:10.1594/WDCC/HOAPS3 _MONTHLY.

[11]Casey K S, Brandon T B, Cornillon P, et al. The past, present, and future of the AVHRR Pathfinder SST program[M]. Oceanogr-aphy from Space，Netherlands： Springer 2010: 273-287.

[12]Krasnopolsky V M, Breaker L C, Gemmill W H. A neural network as a nonlinear transfer function model for retrieving surface wind speeds from the special sensor microwave imager[J]. Journal of Geophysical Research Oceans (1978-2012), 1995, 100(C6): 11033-11045.

[13]Murray F W. On the computation of saturation vapor pressure[J]. Journal of Applied Meteorology, 1967, 6(1): 203-204.

[14]Shie C L. Science background for the reprocessing and Goddard Satellite-based Surface Turbulent Fluxes (GSSTF2c) Data Set for Global Water and Energy Cycle Research[J]. Science Document for the Distributed GSSTF2c via Goddard Earth Sciences (GES) Data and Information Services Center (DISC), 2011: 19.

[15]Shie C L, Hilburn K, Chiu L S, et al. Goddard Satellite Based Surface Turbulent Fluxes[M]// Daily Grid, version 3, Edited by Andrey Savtchenko. Greenbelt, MD, USA: Goddard Earth Science Data and Information Services Center (GES DISC), 2000.

[16]Chou S H, Atlas R M, Shie C L, et al. Estimates of surface humidity and latent heat fluxes over oceans from SSM/I data[J]. Monthly Weather Review, 1995, 123(8): 2405-2425.

[17]Wentz F J, Ricciardulli L, Hilburn K, et al. How much more rain will global warming bring? [J]. Science, 2007, 317(5835): 233-235.

[18]Chou S H, Shie C, Atlas R M, et al. Air-sea fluxes retrieved from special sensor microwave imager data[J]. Journal of geophysical research， 1997, 102(6): 12705-12726.

[19]Yu L, Weller R A. Objectively analyzed air-sea heat fluxes for the global ice-Free oceans (1981—2005)[J]. Bulletin of the American Meteorological Society, 2007, 88(4): 527-539.

[20]Kalnay E, Kanamitsu M, Kistler R, et al. The NCEP/NCAR 40-year reanalysis project[J]. Bulletin of the American meteorological Society, 1996, 77(3): 437-471.

責(zé)任編輯陳呈超

Comparison of Global Air-Sea Latent Heat Flux Products

YANG Xin-Fei, GUAN Lei

(College of Information Science and Engineering, Ocean University of China, Qingdao 266100, China)

Abstract:The interactions between the ocean and atmosphere play an important role in climate system. The ocean influences the atmosphere by transferring heat energy, while the atmosphere changes the ocean currents and redistribute ocean heat through wind stress. The latent heat (LHF) at the air-sea interface is an important source of the atmospheric energy. It is essential for the study of air-sea interaction. The five LHF monthly mean products including the Japanese Ocean Flux Datasets with Used Remote Sensing Observations version 2 (J-OFURO2), the Hamburg Ocean Atmosphere Parameters and Fluxes from Satellite Dataset version 3 (HOAPS3), the Goddard Satellite-Based Surface Turbulent Fluxes version 3 (GSSTF3), the Objectively Analyzed air-sea Fluxes (OAFLUX) and the National Centre for Environmental Prediction and the National Center for Atmospheric Research Reanalysis dataset (NCEP reanalysis) are compared in this paper. The statistics of the intercomparisons among these monthly LHF products over the global ocean during 1996—2005 are analyzed. The results indicate that five LHF products have a similar patterns of global and zonal distribution, larger differences mainly locate in the trade winds zone and high latitude area. Three satellites-derived latent heat flux datasets, i.e. J-OFURO2, HOAPS3, GSSTF3 are close, with the root mean square error (RMSE) between each group of the two products less than 20 W/m2 and the correlations above 0.90. While OAFLUX is closer to NCEP with the RMSE less than 20 W/m2 and the correlations above 0.90. The lowest RMSE is between GSSTF3 and HOAPS3 with the value of 14.08 W/m2, the highest RMSE is between NCEP and GSSTF3 (25.98 W/m2). The lowest correlation is between NCEP and HOAPS3 (0.87) as well as GSSTF3 and OAFLUX (0.87), the highest correlation is between GSSTF3 and HOAPS3 (0.97). Zonal average values of five products reveal similar pattern with differences in specific zones. By analyzing the global differences among the five datasets, the differences among three satellite-derived latent heat fluxes datasets are smaller. The larger difference between satellite-derived datasets and NCEP as well as OAFLUX mainly locates in the equatorial region of the Pacific and the Atlantic, the Northwest Pacific and the Northern Indian Ocean.

Key words:latent heat flux； heat flux； air-sea interface； remote sensing

中圖法分類號(hào)：P229.2

文獻(xiàn)標(biāo)志碼：A

文章編號(hào)：1672-5174(2016)01-123-08

作者簡介：楊新飛(1987-)，女，碩士生。E-mail:yangfei1105@sina.com

收稿日期：2014-03-17；

修訂日期：2015-10-20

*基金項(xiàng)目：全球變化研究國家重大科學(xué)研究計(jì)劃項(xiàng)目(2015CB953900)；國家自然科學(xué)基金項(xiàng)目(41376105)；高等學(xué)校博士學(xué)科點(diǎn)專項(xiàng)科研基金項(xiàng)目(20130132110013)資助。

DOI：10.6441/j.cnki.hdxb.20140090

引用格式：楊新飛， 管磊.全球海氣界面潛熱通量產(chǎn)品對(duì)比研究[J].中國海洋大學(xué)學(xué)報(bào)(自然科學(xué)版)，2016,46(1)：123-130.

YANG Xin-Fei, GUAN Lei.Comparison of global air-sea latent heat flux products[J]. Periodical of Ocean University of China, 2016, 46(1): 123-130.

Supported by Global Change Research Program of China (2015CB953900); National Natural Science Foundation of China (41376105)； Specialized Research Fund for the Doctoral Program of Higher Education (20130132110013)

**通訊作者：E-mail:leiguan@ouc.edu.cn