南黃海海氣熱通量觀測及其與OAflux數(shù)據(jù)集比較研究＊

肖 斌，喬方利，呂連港

（1.國家海洋局 第一海洋研究所，山東 青島266061；2.海洋環(huán)境科學(xué)和數(shù)值模擬國家海洋局重點實驗室，山東 青島266061）

海洋和大氣之間的熱量傳遞是研究海—氣相互作用中的關(guān)鍵問題，其中海氣界面處的海氣熱通量強烈影響著上層海洋以及大氣邊界層結(jié)構(gòu)。研究表明，中國近海感熱通量冬、秋季較強，春、夏季較弱，其地理分布特點是冬季感熱通量的分布隨緯度變化十分明顯，緯度越高感熱通量越大，且等值線分布密集；潛熱通量冬、秋季在臺灣東北部、日本南部和東南部海域形成極大值區(qū)，等值線呈SW—NE走向［1］。在黃海海域水深較淺，受海氣熱交換影響尤其顯著。錢粉蘭等針對黃海、東海海域氣旋發(fā)展過程的熱通量進行過走航觀測［2］，閆俊岳等在南海海域開展了較為全面的海氣通量觀測［3］，然而，在黃海海域連續(xù)的覆蓋海域較廣的熱通量觀測非常有限。Na等［4］通過分析1978—1995年的天氣圖，計算了東亞地區(qū)0.5°×0.5°分辨率的海表熱通量，對黃海海域進行了凈熱通量估計。Hirose等［5］對COARDS（Cooperative Ocean/Atmosphere Research Data Service），NODC（National Oceanographic Data Center），SD2（Oceanographic Station Data），JODC（Japan Oceanographic Data Center）航次報告的通量相關(guān)數(shù)據(jù)進行了分析，計算了東海、黃海海域的熱通量。Chu等［6］利用 MOODS（U.S.Navy’s Master Oceanographic Observation Data Set）和COARDS數(shù)據(jù)集，結(jié)合黃海流系、地形、水團因素將黃海分為若干區(qū)域，分析黃海海域內(nèi)各個區(qū)域內(nèi)熱通量的特征。與東海海域相比，黃海海域可用的數(shù)據(jù)明顯偏少，從而導(dǎo)致在此海域計算得到的熱通量數(shù)據(jù)不準(zhǔn)確。2007年3個航次的熱通量觀測加深了我們對南黃海海域熱通量季節(jié)變化特征的準(zhǔn)確了解。

熱通量的數(shù)據(jù)基本有下面4種：走航觀測、衛(wèi)星數(shù)據(jù)、模式再分析數(shù)據(jù)和浮標(biāo)觀測。各種數(shù)據(jù)都有不同特點。走航觀測數(shù)據(jù)的優(yōu)點是精確度高，一直以來計算海氣通量氣候態(tài)的主要方法［7］，缺點是成本高，數(shù)據(jù)覆蓋海域有限，不能滿足時間和空間連續(xù)性要求，因此利用有限的觀測資料對已有的數(shù)據(jù)集進行比對就成為更有效利用觀測資料的一種方法［8－9］。隨著衛(wèi)星遙感技術(shù)的不斷發(fā)展和提高，國內(nèi)利用衛(wèi)星遙感資料反演中國近海海面熱通量和亞印太邊緣海海氣熱通量產(chǎn)品［10－16］增多。

利用2007年南黃海通量走航觀測中得到的大氣與海洋數(shù)據(jù)（包括近海表空氣溫度、相對濕度、海表溫度和風(fēng)速等），用整體參數(shù)法［17］計算相應(yīng)的海氣熱通量。針對感熱通量、潛熱通量以及通量相關(guān)變量的季節(jié)特點進行分析，并與常用的通量數(shù)據(jù)集進行比對，為黃海海域熱通量研究提供參考。

1 觀測與方法

1.1 觀 測

2007年在中國南黃海進行了3個航次的多學(xué)科綜合海洋調(diào)查（圖1）。調(diào)查中獲得與海氣熱通量相關(guān)的氣象和海洋觀測數(shù)據(jù)。圖1中黑點為3個航次觀測航線的 OAflux（Objectively Analyzed air－sea Flux）數(shù)據(jù)格點及編號。這3個航次都由青島沿岸海域的觀測開始，其中冬季、秋季航次走航順序是在調(diào)查海域由南向北，春季航次中走航順序則是由北向南。冬季航次時間是2007－01－08—02－05，春季航次時間是04－03—27，秋季航次時間為10－06—11－05。關(guān)于此次調(diào)查航次的詳細信息可參見文獻［18］。

主要測量要素包括溫度、濕度探頭（HMP45C）測量空氣溫度和濕度；紅外皮溫傳感器（IRR－P）測量海表皮溫；日輻射強度計測量來自太陽和大氣漫射的短波輻射，波段為310～2 800nm；長波輻射表用來測量長波輻射，波段4.5～40μm。風(fēng)速、氣壓等氣象資料來自船舶氣象儀。

圖1 2007年3個航次的觀測數(shù)據(jù)點及OAflux數(shù)據(jù)格點編號Fig.1 The observed data points of the three cruises in 2007and the grid points of OAflux dataset

1.2 數(shù)據(jù)處理方法

資料在后期經(jīng)過分析匯總處理成了30min平均數(shù)據(jù)，原始的觀測數(shù)據(jù)中包括大霧、下雨、大風(fēng)靠港和設(shè)備維護等非工作狀態(tài)，為提高數(shù)據(jù)可信度，在后期數(shù)據(jù)處理中根據(jù)觀測日志把這些數(shù)據(jù)都去除掉。其中風(fēng)速、氣溫和相對濕度的測量高度均為海面以上15m。

本次觀測數(shù)據(jù)應(yīng)用 COARE3.0（COARE Bulk Flux Algorithm version 3.0b）［19］計算海氣熱通量。COARE3.0 起源于 1992年的 COARE（The Coupled Ocean－Atmosphere Response Experiment）實驗，針對全球、同時兼顧高風(fēng)速的海氣通量版本的發(fā)展過程中進行過大量的實際觀測，包括走航、浮標(biāo)、飛機、衛(wèi)星等，到目前該公式在1～10m/s的風(fēng)速區(qū)域計算精度已較為可靠［19］，對高風(fēng)速條件下的熱通量計算精度也有較大進步。因此，COARE3.0已成為估算海氣通量較可靠的公式之一。其潛熱通量（QLH）與感熱通量（QSH）的表達式：

式中，ρ表示空氣密度；Le表示蒸發(fā)潛熱；cp表示定壓比熱容；U表示在海面高度Zr的風(fēng)速，需要另外說明的是，由于OAflux數(shù)據(jù)集的風(fēng)速是衛(wèi)星反演數(shù)據(jù)，是等價中性風(fēng)速（Equivalent Neutral Wind），雖然實測風(fēng)速已訂正的與OAflux相同的高度，但是沒有考慮實測風(fēng)速與等價中性風(fēng)速的差別，這個差別被認(rèn)為是次要的；ce和ch分別表示潛熱和感熱湍流交換系數(shù)，其為風(fēng)速、測量高度、大氣穩(wěn)定度的函數(shù)；qs和qa分別表示對應(yīng)于海表溫度的大氣飽和比濕和近海表大氣比濕；Ts和θ分別表示海表皮溫和近海表氣溫。需要指出的是，通常提到的海表溫度有2種：一種是海表的浮標(biāo)或CTD測得的海表溫度數(shù)據(jù)（the mixed layer or the bulk SST），另一種則是皮溫（the skin SST）包括衛(wèi)星遙測和紅外皮溫探測等，2種海表溫度存在一定的差別［20－21］，OAflux與本次走航觀測中的海表溫度都屬于皮溫Ts。在式（1）和（2）中只有U，Ts，qa和θ是獨立的變量，這意味著其它的變量包括湍流交換系數(shù)均可通過U，Ts，qa和θ計算得到。

表1，表2分別給出各航次的熱通量和相關(guān)海洋氣象數(shù)據(jù)的平均值，其中輻射通量數(shù)據(jù)參見文獻［10］。

表1 2007年3個航次的輻射通量和海氣熱通量平均值（單位：W·m－2）Table 1 Averaged values of the radiant heat flux and the air－sea heat flux observed in the three cruises in 2007（unit：W·m－2）

表2 2007年的海氣通量相關(guān)數(shù)據(jù)的航次平均值Table 2 Average values of the flux－related meteorological and oceanic data for the three cruises in 2007

2 觀測結(jié)果

2.1 冬季航次

冬季航次中因為天氣原因（避風(fēng)和大霧等）不可避免地需去掉幾個時間段數(shù)據(jù)。2007－01－23—27和01－29—30因受強風(fēng)影響，調(diào)查船靠港；01－12，01－15和01－16均有降雨；感熱通量的觀測變化范圍是－40.1～127.2W/m2（海洋向大氣傳遞為正），潛熱通量變化范圍是－3.7～292.6W/m2，觀測到風(fēng)速變化范圍是0.3～13.0m/s，海表溫度變化范圍是2.6～16.8℃，氣溫變化范圍－3.2～9.4℃，相對濕度變化范圍21.3%～95.4%。感熱通量、潛熱通量、海表溫度和氣溫的最大值均出現(xiàn)在2007－01－12（圖2）。冬季黃海海域的感熱通量與潛熱通量的變化與海表溫度的變化有很大的相關(guān)性。在冬季，黃海海氣熱通量主要從海洋向大氣輸送能量為主，且以潛熱通量方式為主，潛熱通量和感熱通量的航次平均值為80.7W/m2和32.0 W/m2。相比而言，冬季感熱通量是這3個航次中最顯著的，表2中顯示該航次中的平均海氣溫差達3.3℃，為3個航次最高，這是造成感熱通量最為顯著的主要原因。

圖2 2007年冬季航次的海氣熱通量以及通量相關(guān)氣象和海洋觀測數(shù)據(jù)Fig.2 Air－sea heat flux and the flux－related meteorological and oceanic data obtained in the winter cruise in 2007

2.2 春季航次

在春季航次中，2007－04－07—08大霧，04－15下午有降雨，鑒于不滿足輻射計設(shè)備的使用要求停止了通量觀測。觀測到的風(fēng)速范圍是0.4～16.6m/s，海表溫度變化范圍是5.6～15.1℃，氣溫變化范圍為7.4～18.1℃，平均海氣溫差為－2.0℃。航次平均相對濕度達到82.1%。由于平均氣溫高于平均海表溫度，春季南黃海海域的感熱通量轉(zhuǎn)為主要從大氣向海洋輸送，感熱通量的觀測變化范圍是－62.0～20.0W/m2，平均值為－12.5W/m2（圖3）。近海表相對濕度在這3個航次中最大，較大的海表相對濕度抑制了海水的蒸發(fā)相變導(dǎo)致該航次的潛熱通量為最弱，潛熱通量變化范圍是－32.3～103.6W/m2，平均值為5.6W/m2。春季在南黃海海域海氣熱通量是從大氣向海洋輸送能量（3個航次中唯一的航次），其平均海氣熱通量方向由大氣向海洋。在春季，南黃海區(qū)域近海表空氣濕度較大，且大氣溫度高于海表溫度，是造成大氣向海洋輸送能量的主要原因。

圖3 2007年春季航次的海氣通量以及通量相關(guān)氣象和海洋數(shù)據(jù)Fig.3 Air－sea heat flux and the flux－related meteorological and oceanic data obtained in the spring cruise in 2007

2.3 秋季航次

秋季，黃海海域的海氣熱通量完全是由海洋向大氣輸送，海洋—大氣間的熱通量在3個航次中最為強烈，潛熱通量在3個航次中亦是最強，其中感熱通量的變化范圍是5.2～101.4W/m2，平均值為18.9W/m2，潛熱通量（3.7～428.0W/m2）平均值高達142.1W/m2（圖4）。南黃海海域的平均海氣溫差達到2.1℃，其中海表溫度和氣溫的平均值為18.6℃和16.5℃，為3個航次中最高，近海表空氣濕度為最低，平均值為54.4%，這些相關(guān)特征是導(dǎo)致海氣熱通量在3個航次中最為強烈的直接原因。

圖4 2007年秋季航次的海氣通量以及通量相關(guān)氣象和海洋數(shù)據(jù)Fig 4 Air－sea heat flux and the flux－related meteorological and oceanic data obtained in the autumn cruise in 2007

3 對OAflux數(shù)據(jù)集在南黃海海域的評估

OAflux是WHOI的全球通量產(chǎn)品［22］，分辨率是1°×1°。OAflux應(yīng)用了客觀分析方法①YU L S，JIN X，WELLER R A.Multidecade global flux datasets from the objectively analyzed airsea fluxes（OAflux）project：latent and sensible heat fluxes，ocean evaporation，and related surface meteorological variables.Woods Hole Oceanographic Institution，OAFlux Project Technical Report.2008對衛(wèi)星數(shù)據(jù)、模式、再分析數(shù)據(jù)進行了融合，并以COADS（Comprehensive Ocean－Atmosphere Data Set）數(shù)據(jù)作為主要依據(jù)在綜合分析中選擇權(quán)重［23］，應(yīng)用COARE3.0來計算熱通量，是比較常用的海氣通量數(shù)據(jù)集。一些針對中國近海的熱通量研究應(yīng)用了OAflux［24］。

將根據(jù)航次資料計算的海氣熱通量與OAflux數(shù)據(jù)集進行對比，作為對OAflux數(shù)據(jù)集在南黃海海域的評估。由于同樣是應(yīng)用了COARE3.0來計算海氣熱通量，所以可以通過對比通量相關(guān)的氣象和海洋數(shù)據(jù)進一步分析造成偏差的原因。

OAflux數(shù)據(jù)集主要是對潛熱通量和感熱通量的估算，其中長波輻射與短波輻射采用ISCCP（International Satellite Cloud Climatology Project）數(shù)據(jù)［25］。在走航觀測海域內(nèi)有18個 OAflux熱通量數(shù)據(jù)格點，選擇了逐點比對的方式。先選取距離OAflux數(shù)據(jù)格點最近（距離8km左右）的一個走航點，以該走航點為基準(zhǔn)，選擇之前和之后12h的觀測數(shù)據(jù)，將這24h通量的平均值作為該走航點的日平均數(shù)據(jù)與OAflux相應(yīng)數(shù)據(jù)格點上該天的日平均數(shù)據(jù)進行比對。OAflux數(shù)據(jù)集的海表面風(fēng)速、近海表氣溫和近海表大氣濕度的測量高度分別為10m，2m和2m，進行比對之前已經(jīng)將實測數(shù)據(jù)依據(jù)COARE3.0進行了相應(yīng)的高度訂正，其中走航觀測的比濕由相對濕度換算得到。

OAflux數(shù)據(jù)格點編號見圖1，冬季航次的比對結(jié)果見圖5和圖6。由于圖5和圖6中3，8和17號格點附近缺少走航觀測而未能對比。對比結(jié)果顯示，OAflux數(shù)據(jù)集在冬季航次中不論是感熱通量還是潛熱通量，結(jié)果吻合均很好。感熱和潛熱通量的均方差分別為15.3W/m2和21.4W/m2，表明OAflux在2006/2007年冬季南黃海海域有較可靠的海氣熱通量估算。在通量相關(guān)的海洋氣象數(shù)據(jù)的進一步對比中顯示，OAflux數(shù)據(jù)集的比濕較實測數(shù)據(jù)偏低（圖6），但并沒有對熱通量的估算產(chǎn)生較大影響。說明冬季海氣溫差及風(fēng)速對熱通量起著控制作用。

春季航次的對比顯示OAflux潛熱通量明顯比觀測值偏高（圖7），但是隨格點變化有相似的變化趨勢，均方差為28.4W/m2。OAflux感熱通量與觀測數(shù)據(jù)基本吻合，均方差為8.4W/m2。為分析造成潛熱通量偏差的原因，進一步對比了相關(guān)氣象和海洋數(shù)據(jù)（圖8）。可以看出，OAflux的濕度對比觀測中的濕度明顯呈偏低趨勢，均方差為1.0g/kg，這也是造成圖7中OAflux潛熱通量大于實測潛熱通量的主要原因。

秋季，不論是潛熱通量還是感熱通量，OAflux數(shù)值均比走航觀測數(shù)值高（圖9），感熱通量與潛熱通量的均方差分別為20.5W/m2和57.5W/m2。海表溫度的差異是造成上述海氣通量偏差的主要原因（圖10）。海表溫度的對比中顯示了相同的趨勢，2條近似平行的曲線表明，OAflux所用的海表溫度與觀測數(shù)據(jù)存在著系統(tǒng)偏差（均方差為1.4℃），可能對應(yīng)著衛(wèi)星數(shù)據(jù)在近海的系統(tǒng)偏差。為了驗證秋季走航數(shù)據(jù)的可靠性，可以利用測量設(shè)備中長波輻射表的輻射數(shù)據(jù)來推算海表皮溫作為參考［26］。海洋向大氣傳遞的長波輻射包括海洋表面的長波輻射以及向下的長波輻射被海面反射的部分［27］：

式中，ULR為向上長波輻射；DLR為向下的長波輻射；ε為表面輻射系數(shù)，取0.98；σ為Stefan－Boltzmann常數(shù)；α1為長波反照率，取0.02［28］。依據(jù)式（3）利用長波輻射表測量到的向上、向下的長波輻射數(shù)據(jù)可以計算皮溫，計算得到皮溫與紅外皮溫儀的實測數(shù)據(jù)基本一致，沒有系統(tǒng)偏差。將要對比格點的皮溫數(shù)據(jù)挑選出來（圖11），相比OAflux的皮溫數(shù)據(jù)，秋季走航直接觀測和推算的皮溫數(shù)據(jù)保持一致。該結(jié)果進一步表明，秋季南黃海海域OAflux數(shù)據(jù)集的通量與實測通量之間由于海表溫度的偏差導(dǎo)致了通量數(shù)據(jù)差異。3個航次的比對中熱通量和相關(guān)海洋氣象數(shù)據(jù)的均方差見表3。

圖11 紅外皮溫儀觀測、長波輻射換算皮溫和OAflux所用海表溫度數(shù)據(jù)的比較Fig.11 Comparison among the skin SST data observed by IR Radiometer，calculated from long wave radiation and in the OAflux dataset

表3 OAflux數(shù)據(jù)與觀測數(shù)據(jù)的均方根誤差Table 3 The RMS errors between the OAflux dataset and the observed data

4 結(jié)語

基于2007年在南黃海進行的冬、春和秋季3個航次的熱通量觀測，分析了南黃海海氣熱通量的季節(jié)變化特征：春季南黃海有凈熱通量收入，而冬季和秋季為凈熱通量支出。冬、秋季海氣熱通量遠強于春季的。春季南黃海海域近海表空氣濕度遠大于秋季和冬季的，從而導(dǎo)致春季的潛熱通量最弱。秋季南黃海海域海表溫度和近海表空氣溫度較高，而近海表空氣濕度又較低，這些特征導(dǎo)致了秋季海氣熱通量最為強烈。冬季、春季和秋季平均潛熱通量分別為80.7，5.6和142.1W/m2，感熱通量分別為32.0，－12.5和18.9W/m2（海洋向大氣傳遞為正）。

走航觀測計算得到的熱通量與OAflux數(shù)據(jù)集進行對比可知，OAflux在2006—2007年冬季數(shù)據(jù)最為可靠，感熱通量與潛熱通量的均方差分別是15.3W/m2和21.4W/m2。春季OAflux的近海表空氣濕度數(shù)據(jù)與實測數(shù)據(jù)存在較明顯偏差，均方差為1.0g/kg，造成了OAflux潛熱通量與實測潛熱通量存在偏差，感熱通量與潛熱通量的均方差分別為8.4W/m2和28.4W/m2。秋季OAflux的海表溫度與實測數(shù)據(jù)的系統(tǒng)偏差（均方差為1.4℃）造成了其海氣熱通量與實測存在較大偏差，秋季感熱通量與潛熱通量的均方差分別為20.5W/m2和57.5W/m2。對比中所發(fā)現(xiàn)的OAflux數(shù)據(jù)集中春季的近海表空氣濕度偏差、秋季的海表溫度偏差暗示了衛(wèi)星反演數(shù)據(jù)在南黃海海域所存在的問題，這是我們在選用這些數(shù)據(jù)時需要格外注意的。

（References）：

［1］ZHOU M Y，QIAN F L.Model calculation of air－sea heat fluxes over ocean area near China［J］.Acta Oceanologica Sinica，1998，20（6）：21－30.周明煜，錢粉蘭.中國近海及其鄰近海域海氣通量的模式計算［J］.海洋學(xué)報，1998，20（6）：21－30.

［2］QIAN F L，ZHOU M Y，LI S M，et al.Observational study on the meteorological characteristics of the cyclone developing processes over the Huanghai Sea and the East China Sea area［J］.Acta Oceanologica Sinica，2002，24（S1）：77－83.錢粉蘭，周明煜，李詩明，等.黃海、東海海域出海氣旋發(fā)展過程的氣象場特征和熱通量的觀測研究［J］.海洋學(xué)報，2002，24（S1）：77－83.

［3］YAN J Y，LIU J M，JIANG G R，et al.Advances in the study of air－sea flux exchange over the South China Sea［J］.Advances in Earth Sciences，2007，22（7）：685－697.閆俊岳，劉久萌，蔣國榮，等.南海海－氣通量交換研究進展［J］.地球科學(xué)進展，2007，22（7）：685－697.

［4］NA J，SEO J，LIE H J.Annual and seasonal variations of the sea surface heat fluxes in the East Asian marginal seas［J］.Journal of Oceanography，1999，55（2）：257－270.

［5］HIROSE N，LEE H C，YOON J H.Surface heat flux in the East China Sea and the Yellow Sea［J］.Journal of Physical Oceanography，1999，29（3）：401－417.

［6］CHU P，CHEN Y C，KUNINAKE A.Seasonal variability of the Yellow Sea/East China Sea surface fluxes and thermohaline structure［J］.Advances in Atmospheric Sciences，2005，22（1）：1－20.

［7］JOSEY S A，KENT E C，TAYLOR P K.New insights into the ocean heat budget closure problem from analysis of the SOC air－sea flux climatology［J］.Journal of Climate，1999，12：2856－2880.

［8］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.

［9］SUN B，YU L S，WELLER R A.Comparisons of surface meteorology and turbulent heat fluxes over the Atlantic：NWP model analyses versus moored buoy observations［J］.Journal of Climate，2003，16：679－695.

［10］CHEN J N，WANG H N，LüX Y.Variation feature of the sensible heat flux and latent heat flux in South China sea［J］.Advances in Water Science，2007，18（3）：390－397.陳錦年，王宏娜，呂新艷.南海區(qū)域海氣熱通量的變化特征分析［J］.水科學(xué)進展，2007，18（3）：390－397.

［11］CHU J T，CHEN J N，XU L Y.Improved calculation of turbulent heat fluxes at air－sea interface in maritime China［J］.Oceanologia et Limnologia Sinica，2006，37（6）：481－487.褚健婷，陳錦年，許蘭英，等.海－氣界面熱通量算法的研究及在中國近海的應(yīng)用［J］.海洋與湖沼，2006，37（6）：481－487.

［12］CHU J T，CHEN J N，XIA J，et al.An improved algorithm of simulation on air－sea turbulent heat fluxes in China seas［J］.Chinese Journal of Oceanology and Limnology，2007，25（3）：292－299.

［13］CHE J N，WU Y M，HE Y J.The remote of the sensible and latent heat fluxes in China coast［J］.Acta Oceanologica Sinica，2006，28（4）：26－35.陳錦年，伍玉梅，何宜軍.中國近海海氣界面熱通量的反演［J］.海洋學(xué)報，2006，28（4）：26－35.

［14］WANG H N，CHE J N，HE Y J.Variablity of surface heat budget in China coastal seas［J］.Geoscience and Remote Sensing Symposium，2008，4：946－949.

［15］SUN Q Z，CHEN J N，YAN J Y，et al.The variation characteristics of air－sea fluxes over the Xisha area before and after the onset of the South China Sea monsoon in 2008［J］.Acta Oceanologica Sinica，2010，32（4）：12－23.孫啟振，陳錦年，閆俊岳，等.2008年南海季風(fēng)爆發(fā)前后西沙海域海氣通量變化特征［J］.海洋學(xué)報，2010，32（4）：12－23.

［16］CHEN J N，ZUO T，WANG H N.Remote sensing retrieved of air－sea heat fluxes in the Asian and Indian－Pacific ocean area and its analysis［J］.Geoscience and Remote Sensing Symposium，2011，2093－2096.

［17］BLANC T V.Accuracy of bulk－method－determined flux，stability，and sea surface roughness［J］.JGR，1987，92（c4）：3867－3876.

［18］LüL G，YU F，DIAO X Y，et al.Direct observation of radiative flux in the Southern Yellow Sea［J］.Ocean Science Journal，2008，43（2）：115－126.

［19］FAIRALL C W，BRADLEY E F，HARE J E，et al.Bulk parameterization of air－sea fluxes：Updates and verification for the COARE algorithm［J］.Journal of Climate，2003，16：571－591.

［20］FAIRALL C W，BRADLEY E F，GODFREY J S，et al.Cool－skin and warm－layer effects on sea surface temperature［J］.JGR，1996，101（c1）：1295－1308.

［21］SCHLUSSEL P，SHIN H Y，EMERY W J，et al.Comparison of satellite－derived sea surface temperatures with in situ skin measurements［J］.JGR，1987，92（C3）：2859－2874.

［22］YU L S，WELLER R A.Objectively analyzed air－sea heat fluxes（OAFlux）for the global ice－free oceans（1981—2005）［J］.Bull.A－mer.Meteor.Soc.，2007，88：527－539.

［23］YU L S，WELLER R A，SUN B.Improving latent and sensible heat flux estimates for the Atlantic Ocean（1988—1999）by a Synthesis Approach［J］.Journal of Climate，2004，17：373－393.

［24］LIU N，WU D X，LIN X P.Variations of latent heat flux in the East China Sea during the last 50years［J］.Periodical of Ocean University of China，2010，40（4）：19－25.劉娜，吳德星，林霄沛.近50年東中國海潛熱通量的變化［J］.中國海洋大學(xué)學(xué)報，2010，40（4）：19－25.

［25］ROSSOW W B，SCHIFFER R A.Advances in understanding clouds from ISCCP［J］.Bull.Amer.Meteor.Soc.，1999，80：2261－2287.

［26］WEI C J，YU F，LüL G，et al.Comparison among computational methods for atmospheric longwave radiation above sea surface in the Southern Yellow Sea［J］.Advances in Marine Science，2009，27（3）：302－311.魏傳杰，于非，呂連港，等.南黃海海面大氣長波輻射計算方法的比較［J］.海洋科學(xué)進展，2009，27（3）：302－311.

［27］KANTHA L H，CLAYSON C A.Small－scale processes ingeophysical fluid flows［M］.New York：Academic Press，2000.

［28］PAYNE R E.Albedo of the sea surface［J］.Journal of the Atmospheric Sciences，1972，29：959－970.