陽海鵬 王 丹

（1.海軍潛艇學院研究生隊 青島 266071）（2.海軍潛艇學院海洋遙感研究所 青島 266071）

1 引言

海洋表面存在一層很薄的熱邊界層，一般它的溫度比相鄰下界水體的溫度低零點幾攝氏度，正是這個特征，人們形象稱之為“冷表皮”［1］。海洋冷表皮是海洋粘性層的一個重要特征，是海－氣界面熱量、水汽等物理量交換的主要場所。因此，了解海洋冷表皮形成的物理機制，獲得海表的準確溫度，對研究海－氣相互作用，尤其是評估表面熱通量具有重要意義［2］。

2 冷表皮方程［3］

研究海面熱邊界層的模型較多，大體分為兩種（如圖1、圖2）：1）有效擴散冷表皮模型；2）表面剪切冷表皮模型［4］。

圖1 有效擴散冷表皮模型

圖2 表面剪切冷表皮模型

其中，TS為海表溫度，TB為冷表皮下水體溫度，δ代表冷表皮厚度。有效擴散冷表皮模型認為溫度在冷表皮內(nèi)成線性關系，表面冷表皮剪切模型考慮了湍流渦對冷表皮的影響。

熱邊界層模型的選取對冷表皮影響不大，本文采用有效擴散冷表皮模型來計算溫差。

海洋冷表皮受海面風速，太陽輻射，界面溫差、濕度等眾多因素影響。為了推導方便，引入一個變量δ（見圖1、圖2），即冷表皮厚度O（1mm），在該深度內(nèi)，認為熱量傳輸主要以分子熱傳導的方式進行。兩個模型主要差別在于冷表皮底邊界的差異，圖1采用［5］

圖2采用

其中，λ是Saunders比例常數(shù)［6］，由實驗觀測確定，ν是海水的運動粘度（m2s－1），u＊是海水的摩擦速度（ms－1）。在海氣界面，從海面輸送到大氣的總熱通量Q可表示如下：

其中Rnl＝R↓l－R↑l，Rnl為凈的長波輻射，Hs為感熱交換，Hl為潛熱交換，Qsol代表不同深度吸收的太陽輻射，在冷表皮內(nèi)約有10%的太陽輻射被吸收。Q為正代表海洋向大氣輻射熱量，從而使表面溫度下降，密度增大，在冷表皮內(nèi)形成自由對流，且由于冷表皮內(nèi)的溫度梯度由分子熱傳導過程確定，所以Q又可以表示如下：k是海水的熱傳導系數(shù)，z代表垂直坐標軸（原點在海表面，向下為正），T代表溫度剖面。由方程（1）和（4）可計算冷表皮內(nèi)的溫差：

3 模型介紹

計算海洋冷表皮溫度波動，關鍵是確定冷表皮厚度。根據(jù)不同Saunders比例常數(shù)λ，確定了四種計算冷表皮厚度的方法［7］。

1）Paulson和Simpsom［8］于1981年提出的模型（簡稱為PS－81），通過實驗觀測簡單認為λ為常數(shù)

2）Wu［9］于1985年提出的模型（簡稱為 W85），考慮了λ與風速的關系

3）Fairall等人與1996年提出的模型（簡稱為F96），考慮了海水摩擦速度和表面冷卻對流過程

g為重力加速度，α為體積熱膨脹系數(shù)，ρw為海水密度，cw為海水的定壓比熱容，TS為凈的太陽熱通量，QE為潛熱通量，S為鹽度，β為鹽擴散系數(shù)，Lv為海水的蒸發(fā)潛熱，Qb為有效表面制冷通量。

4）Artale等人［10］于2002年提出的模型（簡稱為 A02），設定了參考深度h，并考慮了風速

其中，C為86400s，h為參考深度10m，無量綱因子γ根據(jù)風速的不同取值如下

4 仿真分析

因為海洋向大氣輸送的熱通量Q變化沒有風速變化劇烈，所以假定熱通量為常數(shù)，對四種不同模型，仿真所得冷表皮厚度如圖3所示。

由圖3可知，冷表皮厚度量級為毫米級，四種模型計算的基本趨勢一致，隨風速增大而減小。仿真結果在低風速下（U＜4m／s）差別較大，在風速大于6m／s時，相差很小，可忽略不計。因此，各模型的主要差別在于低風速下仿真海表溫度，相對PS81，W85模型來說，F(xiàn)96和A02模型更優(yōu)，從其表達式可知，F(xiàn)96和A02模型考慮了在低風速下，由于冷卻導致的密度增加引起的對流效應，使不同溫度的上下水體混合，有效抑制了海表冷卻的進程，理論上更符合實際情況，且在風速趨近于0時，其值更加合理，不會趨向于無窮大等失真值。從仿真及理論分析可知，F(xiàn)96和A02的模型較好，但也需要實測海洋數(shù)據(jù)或?qū)嶒炇襾眚炞C模型的可信度。根據(jù)F96和A02模型計算平均冷表皮厚度如圖4所示。

圖3 冷表皮厚度隨風速的變化

圖4 平均冷表皮厚度隨風速的變化

由圖4可知，平均冷表皮厚度隨10m處的風速增加而減少，其平均值如表1所示。根據(jù)不同風速計算所得平均冷表皮厚度，仿真冷表皮內(nèi)溫差如下：

表1 冷表皮厚度與風速關系

由圖5可得，其平均冷表皮溫差隨10m處風速增加而減少。大的風速會對海洋表面產(chǎn)生擾動，促進冷表皮內(nèi)的對流，從而減弱冷表皮效應，有時甚至會破壞冷表皮。其平均值如表2所示。

圖5 平均冷表皮溫差隨10m處風速的變化

表2 平均冷表皮溫差與風速關系

5 結語

本文主要對海洋冷表皮機理進行介紹，并對四種模型仿真冷表皮厚度進行了簡單分析?？傻萌缦陆Y論：

1）冷表皮厚度范圍為0～5mm，冷表皮內(nèi)溫差范圍為0～1.5℃，與實際海洋情況基本一致；

2）冷表皮厚度和冷表皮內(nèi)溫差都隨風速的增大而減少，四種模型所得冷表皮厚度在高風速條件下，結果一致，在低風速下，相差較大；

3）從理論上來說，F(xiàn)96和A02模型考慮了表面冷卻導致的對流過程，更加合理。

但四個模型都假定海表熱通量為常數(shù)，沒有考慮到風速會對熱通量產(chǎn)生影響，這是模型需改進的地方。

［1］Fabrice Veron，W.Kendall Melville，Luc Lenain，Infrared techniques for measuring ocean surface processes，Journal of atmospheric and oceanic technology，VOL（25），2008：307-326.

［2］Donlon，C.J.，and I.S.Robinson，observation of the ocean thermal skin in the Atlantic Ocean.J.Geophys，Res.，VOL（102），1997：18585-18606.

［3］Fairall，C.W.，Bradley，E.F.，Godfrey，J.S.，Cool-skin and warm-layer effects on sea surface temperature，Journal of Geophysical Research，VOL（101），NO.C1，1996：1295-1308.

［4］Wells，A.J.，Claudia Cenedese，F(xiàn)arrar，J.T.，et al.Variations in ocean surface temperature due to near-surface flow：straining the cool skin layer［J］.Journal of physical ocean-ography，2009，vol（39）：2685-2710.

［5］Osborne，M.F.M.，The effect of convergent and divergent flow patterns on Infrared and Optical Radiation from the sea［J］.Dtsch.Hydrogr.Z.，1965，18：1-25.

［6］Saunders，P.M.，The temperature at the ocean-air interface.J.Atmos.Sci.，vol（24），1967：269-273.

［7］Chia-Ying Tu，Ben-Jei Tsuang，Cool-skin simulation by a onecolumn ocean model.Geophysical Research Letters，VOL.32，2005：237-245.

［8］Paulson，C.A.，and Simpson，J.J.，The temperature difference across cool skin of the ocean，J.Geophys.Res.，86（C11），1981：11044-11054.

［9］Wu，J.，On the cool skin of the ocean，Boundary Layer Meteorol.，31，1985：203-207.

［10］Artale，V.，Iudicone，D.，Santoleri，R.，Role of surface fluxes in ocean general circulation models using satellite sea surface temperature：Validation of and sensitivity to the forcing frequency of the Mediterranean thermohaline circulation，J.Geophys.Res.，107（C8），2002：3120-3125.