汪晉，廖前鋒，閆申，王毅，項杰

(1.南京工業(yè)大學浦江學院，南京211222；2.陸軍工程大學教務(wù)處，南京210007；3.國防科技大學前沿交叉學科學院，長沙410028；4.國防科技大學氣象海洋學院，長沙410028)

引 言

大氣邊界層(Atmospheric Boundary Layer，ABL)位于對流層底部，是最接近地球表面的一層大氣。由于受到地表摩擦和非絕熱加熱等的影響，ABL具有明顯不同于其上自由大氣的特征，如，有別于自由大氣中的平流特征，ABL 內(nèi)大氣的運動形式主要為湍流，會對動量、熱量和水汽等的輸送過程產(chǎn)生重要影響(Stull，1988；趙鳴和苗曼倩，1992；Garratt，1992；盛裴軒等，2013)。研究ABL過程對了解大尺度天氣過程的演變、研究長期預報和氣候理論等問題，具有重要的意義(高登義，1994；周明煜等，2000；張強和胡隱樵，2001；胡非等，2003；Medeiros et al.，2005)。此外，人類主要生活于ABL 內(nèi)，ABL 對人們的生產(chǎn)、生活等各種活動有重要影響。如ABL與空氣污染密切相關(guān)(任桂萍等，2021)，ABL 中的湍流運動對于空氣污染物具有擴散和混合的作用(因而在空氣污染氣象學中ABL又稱為混合層)；ABL高度決定了污染物擴散稀釋的潛在空間范圍，ABL 高度高時，污染物可以在更大的空間內(nèi)擴散稀釋，從而降低濃度。因此，ABL 高度與近地面污染物濃度呈負相關(guān)(Du et al.，2013)，重污染過程往往伴隨較低的ABL高度(Qu et al.，2017)；同時，作為空氣污染物主要成分的氣溶膠對ABL 也有反饋作用(Yu et al.，2002；Li et al.，2017)；對于都市ABL來說，高的顆粒物濃度會加強ABL的穩(wěn)定性，這會進一步地降低ABL 的高度，導致顆粒物濃度進一步增加(Petaja et al.，2016)。所以，研究ABL 過程對于保護環(huán)境以及保障人民的身心健康具有重要的作用。

ABL 與降水(尤其是強降水，如暴雨)過程有密切的關(guān)系。研究表明，ABL 中的物理過程，如ABL 中的動量、熱量、水汽的輸送、ABL中的摩擦效應(yīng)以及地形的作用等對于暴雨的發(fā)生發(fā)展有重要的影響(Benjamin and Carlson，1986；趙鳴，2008；陳鵬等，2014)。由于ABL中湍流運動的尺度很小(屬于次網(wǎng)格尺度的運動)，大氣數(shù)值天氣(或氣候)模式難以解析這種尺度的運動，通常需要進行參數(shù)化處理。因此，ABL 過程的參數(shù)化一直是數(shù)值天氣預報模式及數(shù)值氣候模式研發(fā)的關(guān)鍵問題之一，也是暴雨模擬研究中的熱點問題(蔡薌寧等，2006；盛杰和林永輝，2010；沈新勇等，2017；程銳等，2019)。ABL高度(或深度)是ABL的一個重要參數(shù)，估計準確的ABL高度數(shù)值對于數(shù)值天氣(或氣候)模式具有重要的意義。一方面，ABL高度數(shù)值是ABL 本身的參數(shù)化方案所必須的；另一方面，它也是數(shù)值天氣(或氣候)模式中其他的物理過程參數(shù)化所需要的。同時，ABL高度信息也是空氣污染物預報模式的輸入?yún)?shù)。準確的ABL 高度數(shù)據(jù)對于提高數(shù)值天氣預報精度、化學成分預報精度等有重要影響。因此，確定ABL高度在氣象學上具有重要意義(潘云仙和蔣維楣，1982；洪鐘祥等，1998；李茂善等，2006；伍大洲等，2006；張璐等，2011；孔揚和張?zhí)K平，2014)。

根據(jù)大氣的層結(jié)穩(wěn)定性ABL分為不穩(wěn)定ABL、穩(wěn)定ABL 和中性ABL。不穩(wěn)定ABL 又稱為對流ABL(Garratt，1992)，通常出現(xiàn)在地表受熱產(chǎn)生對流或云頂輻射冷卻產(chǎn)生對流的情形下，其頂部往往伴隨著逆溫和(或)濕度的急劇下降；穩(wěn)定ABL通常在夜間出現(xiàn)，常常伴隨著以地表為底的逆溫層(Garratt，1992)；中性ABL 較少出現(xiàn)，如陸地夜間ABL 中的剩余層(Residual Layer，RL)就屬于中性ABL(Stull，1988)。不穩(wěn)定ABL的頂(其高度即為ABL高度)常?？梢愿鶕?jù)一些物理量在ABL頂具有顯著的梯度(如逆溫、水汽含量的急劇下降等)而確定，而穩(wěn)定ABL 頂?shù)拇_定卻沒有統(tǒng)一的標準。在實際應(yīng)用中，穩(wěn)定ABL頂可以根據(jù)需要采用不同的方式定義，如定義為以地表為底的逆溫層的頂(Seidel et al.，2010)，或者是湍流強度占近地面湍流強度的一小部分比例的位置(Stull，1988)，或者是利用Richardson 數(shù) 定 義(Vogelezang and Holtslag，1996)。不 穩(wěn) 定ABL 的高度從數(shù)百米到數(shù)千米不等，而由于穩(wěn)定ABL中的湍流擴散弱，穩(wěn)定ABL 的高度最多不超過幾百米(Garratt，1994)。

確定ABL 高度的研究工作有很多。張宏昇等(2020)對ABL 高度的確定作了很好的總結(jié)，介紹了兩種類型的確定ABL高度的研究工作，其一是利用遙感觀測數(shù)據(jù)確定ABL高度，重點是地基遙感觀測數(shù)據(jù)(如激光雷達/云高儀、地基微波輻射計、聲雷達、風廓線雷達、多普勒測風激光雷達等數(shù)據(jù))，對于衛(wèi)星遙感數(shù)據(jù)只是簡單介紹了掩星探測數(shù)據(jù)的應(yīng)用；其二是利用預報方程、診斷公式通過參數(shù)化的方法計算ABL 高度，其優(yōu)點是只需輸入少量數(shù)據(jù)即可計算得出ABL高度，盡管方法簡便實用，但所得結(jié)果精度相對較差，在缺乏大氣廓線觀測值的地區(qū)參數(shù)化方法有較高的應(yīng)用價值。至于ABL高度數(shù)據(jù)的應(yīng)用，該文獻只在引言中籠統(tǒng)地提及了在天氣預報和空氣質(zhì)量預報模式中有應(yīng)用，并沒有詳細的研究工作介紹。本文從不同角度介紹了ABL高度的確定及其應(yīng)用，即針對四種常見類型的觀測資料，包括常規(guī)探空資料、衛(wèi)星遙感探測資料、激光雷達探測資料及全球?qū)Ш叫l(wèi)星系統(tǒng)GNSS(Global Navigation Satellite System)掩星探測資料，介紹確定ABL高度的主要方法及其進展情況；通過對不同方法和結(jié)果進行比較，提出在同時具有多種觀測資料時如何處理ABL 高度的三個基本原則；最后，以目前廣泛應(yīng)用的數(shù)值天氣預報模式WRF (Weather Research and Forecasting)為例，重點介紹ABL 高度數(shù)據(jù)在數(shù)值天氣(或氣候)預報模式中的應(yīng)用，以期為改進、完善數(shù)值天氣(或氣候)預報模式中的ABL參數(shù)化方案提供有益的參考依據(jù)。

1 利用常規(guī)探空資料確定ABL高度

對于不穩(wěn)定ABL來說，ABL頂往往伴隨著一個或多個物理量存在明顯的梯度，如溫度(逆溫)、濕度、氣溶膠粒子濃度、臭氧濃度、湍流強度和折射率等。因此，利用這些物理量的單個探測廓線或模式輸出廓線就可以確定ABL 頂所在的高度，即ABL高度。常規(guī)探空數(shù)據(jù)一直是氣象學家確定ABL 高度的主要數(shù)據(jù)源(Seibert et al.，2000；徐桂榮等，2014；Dai et al.，2014)，基于常規(guī)探空廓線確定ABL 高度的方法主要包括氣塊法、最大梯度法、湍流法及Richardson 數(shù)法等四種，下面分別進行簡要介紹。

1.1 氣塊法

Holzworth(1964)提出了“氣塊法”的基本思想，把從地面(具有最高溫度的)某點出發(fā)的干絕熱廓線與該處(或臨近的)探空溫度廓線的交點處的高度定義為ABL高度，即在這個高度上與在地面具有相同的虛位溫值。利用該方法確定的ABL 高度也稱為“混合高度”，通常用于空氣污染物的預報以評估在ABL 中釋放的污染物的濃度。然而用這種方法確定的結(jié)果強烈依賴于地面溫度，會出現(xiàn)對流ABL頂沒有逆溫的情形(Seibert et al.，2000)。因此，Beljaars和Betts(1992)提出了改進的氣塊法，并且被Wotawa 等(1996)和Seibert等(2000)所應(yīng)用。氣塊法是一種基本的方法，主要用于不穩(wěn)定ABL的情形。

1.2 最大梯度法

ABL 頂往往伴隨著逆溫和(或)濕度等的急劇下降，通常利用位溫最大垂直梯度來確定ABL 高度(Stull，1988；Garratt，1992)，這種位溫最大垂直梯度層表示從對流不穩(wěn)定區(qū)域(位于下方)到位于上方的穩(wěn)定區(qū)域的過渡區(qū)。Ao等(2012)直接利用溫度或濕度的最大垂直梯度(絕對值)確定ABL 高度，具體實施時需要首先計算位溫(或溫度、濕度)的垂直梯度，由于觀測數(shù)據(jù)往往含有噪聲，再加上算法的原因，計算得到的垂直梯度可能含有較大的誤差，導致結(jié)果不可靠，甚至是錯誤的。針對這種實際情況，有學者在計算垂直梯度之前先對數(shù)據(jù)進行濾波等平滑處理(Amiridis et al.，2007；楊富燕等，2016)；Zeng 等(2004)在最大梯度法中采用了閾值判據(jù)，即對于海洋上的不穩(wěn)定ABL，把虛位溫的垂直梯度首次大于或等于3 K·km-1(稱為閾值)所在的高度作為ABL高度；Dai等(2014)利用位溫垂直梯度廓線來確定對流邊界層(Convective Boundary Layer，CBL)或云蓋頂邊界層(Cloud-Topped Boundary Layer，CTBL)高度：對于CBL 來說，其頂部逆溫層相對較弱，故設(shè)置相對較小的位溫梯度閾值來確定ABL高度(即在逆溫層底部首次超過閾值的高度)；而對于CTBL來說，其頂部逆溫明顯，故可以設(shè)置相對較大的位溫梯度閾值來確定ABL高度，或者直接取位溫梯度最大值所在的高度為ABL 高度。為了減少觀測誤差對計算梯度的影響，還可以采用數(shù)值微分方法(Cheng et al.，2003；Xiang et al.，2021)。

在最大梯度法中采用閾值判據(jù)時，閾值的選擇具有一定的主觀性，這會產(chǎn)生估計的ABL高度的不確定性。最大梯度法既可用于不穩(wěn)定ABL 也可用于穩(wěn)定ABL，但是對于穩(wěn)定ABL，最大梯度法應(yīng)用起來更復雜(Dai et al.，2014)。目前，最大梯度法仍然是確定ABL高度的主要方法，應(yīng)用廣泛。

1.3 湍流法

ABL 與自由大氣之間的主要區(qū)別在于ABL 具有湍流性，湍流能夠使ABL 內(nèi)產(chǎn)生混合過程，因而ABL能夠響應(yīng)地表的強迫。因此，ABL可定義為湍流擾動層，而ABL 高度就為湍流消失的高度，也即存在連續(xù)湍流的最大高度(Dai et al.，2011)。探測ABL 高度的湍流方法就是利用儀器(如放置于飛機上)直接(高頻)觀測垂直風速，得到垂直風速廓線，然后對風速廓線進行高通濾波(如小波濾波器)，消除低頻分量，得到風速脈動值，風速脈動值顯著減少的高度就是ABL 高 度。Wang 等(1999)、Dai 等(2011)、Zhang 等(2014)均利用了湍流方法對外場觀測數(shù)據(jù)確定了ABL高度。

湍流方法被認為是確定ABL 高度的最準確的方法，對任何類型的大氣層結(jié)(不穩(wěn)定、穩(wěn)定、中性)都適用。但是，湍流直接觀測中所使用的設(shè)備代價昂貴，所以除了在科學實驗中偶爾使用以外，日常觀測很少使用這種設(shè)備，這決定了湍流方法在確定ABL高度的工作中應(yīng)用較少。

1.4 Richardson數(shù)法

Richardson 數(shù)是一個無量綱數(shù)，表示浮力抑制的湍流與風切變產(chǎn)生的湍流之比，是判定湍流產(chǎn)生與否的重要參數(shù)，也是判定流體靜力穩(wěn)定性的指標，可直接用來估算ABL 高度，這就是確定ABL 高度的Richardson數(shù)法。

Richardson 數(shù)法既可應(yīng)用于不穩(wěn)定ABL，也可應(yīng)用于穩(wěn)定ABL(Georgoulias et al.，2009)，特別適合于在數(shù)值天氣預報模式中確定ABL 高度(Straume et al.1998；Zilitinkevich and Baklanov，2002；Jericevic and Grisogono，2006)，因為數(shù)值模式中需要進行自動處理。Richardson數(shù)法也可應(yīng)用于從探空資料確定ABL高度(Zeng et al.，2004；Balsley et al.，2006；Hennemuth and Lammert，2006；Sicard et al.，2006)。如Guo 等(2019)利用中國區(qū)域1976—2016 年的無線電探空資料，采用整體Richardson數(shù)法確定了ABL高度，并分析了ABL高度的時空變化特征。

Richardson 數(shù)法建立在這樣一個基本假設(shè)之上，即當Richardson 數(shù)超過某個臨界值時連續(xù)湍流會終止。這種方法的一個弊端是它的臨界值選取的不確定性，臨界值選取可從0.15變化至0.55，而且一旦選定以后，這個臨界值對于所有類型的ABL(不管是不穩(wěn)定的、穩(wěn)定的還是中性的)都是相同的，這不太合理；同時，Richardson 數(shù)法并不太適用于大多數(shù)被動的空基反演的溫度廓線(紅外和微波)，因為被動獲取的廓線的垂直分辨率最多為1～2 km，對于估計ABL 高度，這是一種非常粗略的分辨率。為了彌補Richardson數(shù)法的不足，Zhang等(2014)利用四個外場觀測實驗的數(shù)據(jù)對Richardson 數(shù)法的臨界值進行了修正，把四個外場實驗中觀測到的ABL分為強穩(wěn)定ABL、弱穩(wěn)定ABL和不穩(wěn)定ABL三類，對每一類都利用線性擬合和統(tǒng)計誤差最小化方法得到各自最優(yōu)的Richardson 數(shù)臨界值。這樣，對不同穩(wěn)定性的ABL采用不同的臨界值就可以得到更準確的ABL 高度估計值。Zhang 等(2014)的Richardson數(shù)臨界值選擇方案值得進一步發(fā)展應(yīng)用。

2 利用常規(guī)衛(wèi)星遙感資料確定ABL高度

利用衛(wèi)星遙感資料確定ABL 高度的優(yōu)勢在于可得到其時空變化的全球和區(qū)域統(tǒng)計特征，許多研究在這方面提出了不同的方法和技巧。Minnis 等(1992)、Wood 和Bretherton (2004)、Zuidema 等(2009)先后利用副熱帶海洋上空層云和層積云頂?shù)臏囟群秃１頊囟?sea surface temperature，SST)來估計ABL 高度，其原理如下：ABL內(nèi)的溫度廓線與干絕熱廓線很接近，并且在大多數(shù)海洋上海氣交界面處的海氣溫度差異是很小的，因此可以使用海表溫度SST結(jié)合干絕熱大氣的溫度直減率得到ABL高度。從氣候?qū)W研究的實用性角度考慮，這個方法只適用于完整云塊的區(qū)域，對于陸地或深度淺的積云區(qū)域效果并不好。因此，有人提出了直接利用云頂高度估計的改進方法，Wu 等(2008)使用空基激光雷達估計云頂高度來確定ABL 高度，而Karlsson等(2010)利用多角度成像光譜輻射計估計云頂高度來確定ABL高度。但是，以上這些方法都面臨著同樣的問題，即塊狀云的頂并不總是與ABL頂一致。

另一些研究人員則直接利用衛(wèi)星遙感反演的溫濕廓線數(shù)據(jù)確定ABL 高度。Fetzer 等(2004)利用大氣紅外探測儀AIRS反演的溫度和濕度廓線資料來確定位于ABL頂?shù)哪鏈匚恢?，因而得到ABL高度。該方法的不足之處是，它僅能給出ABL高度的粗略估計，因為廓線的垂直分辨率是很低的(對溫度來說是約1 km，對濕度來說是約2 km)，不足以解析ABL頂?shù)募毼⒔Y(jié)構(gòu)。

因此，利用常規(guī)衛(wèi)星遙感資料確定ABL高度還存在明顯的不足，這方面的研究工作相對較少。

3 利用激光雷達探測資料確定ABL高度

激光雷達的主要探測數(shù)據(jù)為后向散射回波信號強度，它受氣溶膠和水汽的影響大。許多影響光散射的氣溶膠粒子來自于地球表面(如海鹽、灰塵和土壤顆粒)，因為ABL 頂通常具有逆溫特性(屬于穩(wěn)定層結(jié))，這限制了這些氣溶膠粒子向上的擴散，因此氣溶膠粒子主要分布于ABL內(nèi)，在ABL和其上的清潔大氣之間氣溶膠粒子濃度有一個顯著的下降，反映在激光雷達回波信號中，就是后向散射回波信號強度在ABL頂也有 一 個 顯 著 的 下 降(Boers et al.，1984；Hooper and Eloranta，1986)，這 一 特 征 可 用 來 確 定ABL 高 度。Hooper 和Eloranta(1986)利用地基激光雷達探測ABL風速時也探測了ABL 深度(即ABL 高度)，提出了確定ABL深度的后向散射方差法(backscatter variance，也稱為標準差法)。Steyn等(1999)、Hageli等(2000)、Eresmaa等(2012)提出了曲線擬合方法，用理想的激光雷達后向散射強度廓線擬合實際的后向散射觀測數(shù)據(jù)來確定混合層深度和(或)夾卷層厚度。Davis 等(2000)、Brooks(2003)提出利用小波協(xié)方差變換方法確定ABL高度。Sawyer 和Li (2013)、Hicks 等(2015)還將曲線擬合方法和小波協(xié)方差變換方法結(jié)合起來確定ABL 高度。這些方法本質(zhì)上都是最大梯度方法。由于激光雷達具有高時空分辨率和較高精度的優(yōu)勢，利用激光雷達探測ABL高度近年來受到了廣泛的重視(Wang et al.，2012；王琳等，2012；Compton et al.，2013；李紅等，2015；Lewis et al.，2013；楊富燕等，2016)，它既可以用來探測陸地ABL 高度(McGrath-Spangler and Denning，2012)，也可以用來探測海洋ABL 高度(Palm，2005)，還可以探測全球ABL高度，如星載激光雷達(Zhang et al.，2016)，穩(wěn)定與不穩(wěn)定ABL 都可以探測，是一個前沿的研究方向。

4 利用GNSS掩星探測資料確定ABL高度

作為一種新型的大氣遙感探測手段，GNSS 掩星探測技術(shù)具有高精度、高垂直分辨率、全球準均勻覆蓋、全天候觀測、長期穩(wěn)定等特點(Kursinski et al.，1996；Kuo et al.，2004)，在國際上受到廣泛的重視。一些國家啟動了眾多的掩星探測項目，這些項目獲得了大量的中性大氣和電離層電子密度廓線數(shù)據(jù)，這些數(shù)據(jù)在數(shù)值天氣預報、氣候及空間天氣研究和預報等方面得到了廣泛的應(yīng)用。如我國臺灣地區(qū)與美國聯(lián)合實施的掩星探測項目“氣象、電離層和氣候星座觀測系統(tǒng)(Constellation Observing System for Meteorology，Ionosphere and Climate，COSMIC，2006-2011)”成為最成功的掩星探測業(yè)務(wù)應(yīng)用系統(tǒng)。在COSMIC 項目之后，美國和臺灣地區(qū)又推出了名為FORMOSAT-7/COSMIC-2的后續(xù)掩星探測項目，COSMIC-2項目每天能跟蹤超過4 000多條高質(zhì)量廓線。

由于GNSS 掩星探測技術(shù)具有獨特的優(yōu)點，特別適合于氣候?qū)W研究，因此不少學者利用掩星探測數(shù)據(jù)產(chǎn)品來確定ABL 高度。von Engeln 等(2005)利用全譜反演中的截斷高度來確定ABL 高度；Sokolovskiy 等(2006)利用線性回歸方法估計折射率廓線的垂直導數(shù)，通過確定折射率垂直梯度的最大下降點(斷點，即不連續(xù)點)來確定ABL 高度；Sokolovskiy 等(2007)又提出利用彎角廓線確定ABL高度，其思想是確定固定長度范圍(一般取300 m)內(nèi)的彎角最大下降值。Guo 等(2011)采用斷點方法，利用一年的COSMIC掩星探測折射率數(shù)據(jù)得到了熱帶和亞熱帶海洋地區(qū)ABL 高度的時空分布，并且空間分布特征與歐洲中期天氣預報中心ECMWF(European Center for Medium-Range Weather Forecasts)全球分析資料的結(jié)果一致，但是平均來說前者的量值要大于后者的量值。Ao等(2012)把折射率和水汽壓垂直梯度的最小值點看作為ABL頂，其高度定義為ABL高度，并且引入了相對最小梯度(或顯著參數(shù))來表征所得ABL 高度的質(zhì)量指標，結(jié)果表明由COSMIC 掩星數(shù)據(jù)確定的ABL 高度的平均值與ECMWF 再分析資料(European Center for Medium-Range Weather Forecasts reanalysis interim，ERA-Int)的結(jié)果具有相似的空間分布和季節(jié)變化特征，但是前者比后者要高500 m，且COSMIC 掩星數(shù)據(jù)的ABL 高度的標準差也高于ERA-Int的相應(yīng)結(jié)果。其他一些學者也進行 了 類 似 的 工 作(Chan and Wood，2013；Ho et al.，2015；廖麒翔等，2015；劉艷等，2015)。Basha 等(2018)利用小波協(xié)方差變換方法基于全球多個掩星探測計劃數(shù)據(jù)得到了全球ABL高度的氣候態(tài)分布，這些分布顯示了ABL 高度具有明顯的季節(jié)變化，其中夏季最高，春秋次之，冬季最低。周文等(2018)利用COSMIC掩星資料，通過計算大氣折射率最小梯度來確定青藏高原地區(qū)ABL高度，并用無線電探空資料對結(jié)果進行了檢驗。而Yan 等(2019)、周建印等(2020)、Xiang 等(2021)則利用數(shù)值微分方法結(jié)合反問題中的正則化策略來計算彎角、折射率廓線的垂直梯度進而確定ABL高度，該方法在理論上更嚴謹，適用范圍更廣。

利用GNSS掩星探測資料確定ABL高度的不足之處有：掩星探測未必能達到地表附近，因此對于淺薄的ABL(尤其是穩(wěn)定ABL)，掩星探測技術(shù)不太容易探測到它們的高度。

5 確定ABL高度的不同方法和結(jié)果的比較

上文介紹了用于確定ABL 高度的各種數(shù)據(jù)及主要方法。一些學者對不同方法及不同數(shù)據(jù)的結(jié)果進行了比較研究。如Seidel等(2010)利用全球10 a的505個探空站的觀測資料(包括正常分辨率和高分辨率的廓線數(shù)據(jù))，對七種方法(氣塊法、位溫的最大垂直梯度法、比濕的最小垂直梯度法、相對濕度的最小垂直梯度法、折射率的最小垂直梯度法、懸空的逆溫層法、以地面為底的逆溫層法即穩(wěn)定ABL法)確定ABL高度進行了對比分析，探討了ABL高度氣候態(tài)數(shù)據(jù)中的參數(shù)不確定性(指溫度、位溫、虛位溫、相對濕度、比濕、折射率數(shù)據(jù)所得結(jié)果的差異)和結(jié)構(gòu)不確定性(指所用的方法和數(shù)據(jù)分辨率的不同而產(chǎn)生的差異性)，得到了一些有意義的結(jié)果。例如，不同方法所得結(jié)果的差異達到了數(shù)百米，ABL 高度對探空數(shù)據(jù)的分辨率更敏感，標準分辨率數(shù)據(jù)的結(jié)果比高分辨率數(shù)據(jù)的結(jié)果偏大等等。Gu等(2021)用11 a的COSMIC全球無線電掩星數(shù)據(jù)(小波協(xié)方差變換法)、無線電探空數(shù)據(jù)(Integrated Global Radiosonde Archive，IGRA) (整體Richardson 數(shù)法)以及歐洲中心再分析資料ERA-Int(氣塊法)對全球ABL高度進行了估算比較，結(jié)果表明ERA-Int和IGRA數(shù)據(jù)估算得到的全球ABL高度的空間分布較為一致，相關(guān)性較好；由COSMIC 掩星數(shù)據(jù)估算得到的ABL 高度比探空數(shù)據(jù)和再分析數(shù)據(jù)估算的結(jié)果系統(tǒng)性地偏大。

總之，用于確定ABL 高度的各種數(shù)據(jù)、各種方法都有各自的適應(yīng)性和優(yōu)缺點，利用不同方法、不同數(shù)據(jù)所得到的ABL高度具有差異性(不確定性)。如果同時利用多種觀測手段獲得了多種觀測資料，如濕度、溫度、臭氧濃度、湍流強度、示蹤物(如氣溶膠粒子)、折射率等，那么可以利用這些資料得到多個ABL高度的估計值，而且一般來說，這些估計值都是不相同的。如何處理ABL高度同時具有多個估計值的情形，是一個在實際應(yīng)用中必須要考慮的問題?？梢园凑找韵氯齻€原則來處理ABL 高度具有多個估計值的“窘境”：(1)準確性原則。觀測資料的準確性是其被采用的主要原則。如湍流法被認為是確定ABL 高度的最準確的方法，如果有湍流強度的觀測資料，則應(yīng)該首先利用湍流強度觀測資料確定ABL 高度；又如，衛(wèi)星遙感反演的溫度相對濕度來說準確性更高，因此在沒有湍流強度觀測資料的情況下可以利用衛(wèi)星遙感反演的溫度廓線信息估計ABL 高度。如果有無線電探空的溫度濕度廓線，由于探空的溫濕廓線比衛(wèi)星遙感反演的溫濕廓線精度高，則可以利用探空的溫濕廓線確定ABL高度。(2)目的性原則。根據(jù)研究任務(wù)目的的不同而選擇合適的觀測資料來確定ABL 高度。如在空氣質(zhì)量的模擬與預報中，可以采用氣塊法確定ABL 高度；而在ABL 高度的氣候?qū)W研究中，可以利用GNSS 掩星探測折射率廓線數(shù)據(jù)來確定ABL 高度，因為折射率同時含有溫度和水汽的信息，而且GNSS 掩星探測廓線數(shù)據(jù)具有高精度、高垂直分辨率、全球準均勻覆蓋等優(yōu)點，特別適合于氣候?qū)W研究。(3)科學性原則。從科學上探討ABL 高度出現(xiàn)多個估計值的原因。

6 ABL高度信息在數(shù)值天氣預報模式中的應(yīng)用

如前所述，數(shù)值天氣預報模式都包括ABL 過程。ABL 過程的湍流特性決定了ABL 運動方程組的閉合主要采用參數(shù)化(即次網(wǎng)格尺度的湍流通量利用預報的平均變量通過垂直擴散方程來表示)來實現(xiàn)。參數(shù)化可以分為不同的形式，如1 階閉合(利用K-理論)、1.5階閉合(閉合后的方程中不僅包括1階閉合情形的標準變量，還包括湍流的位溫方差和湍動能)、2 階閉合(閉合后的方程中只含有湍流的2階矩量)等，也可以分為局地閉合和非局地閉合(根據(jù)參數(shù)化關(guān)系式是否僅使用局地的信息確定)。目前廣泛應(yīng)用的數(shù)值天氣預報模式WRF 包含有多種ABL 參數(shù)化方案(Shin and Hong，2011)。常用的有：局地閉合方案，包括MYJ(Mellor-Yamada-Janjic)，QNSE (quasi-normal scale elimination)，Boulac (Bougeault-LaCarrere)，等 等(為1.5 階 閉合)；非局地閉合方案，包括YSU (Yonsei University)，ACM2 (new version of the asymmetric convective model(ACM))，等等(為1階閉合)。ABL高度信息在數(shù)值天氣預報模式中的作用主要體現(xiàn)在：在ABL過程的參數(shù)化方案中，它是大多數(shù)參數(shù)化方案所必須的參量。

例如，在YSU 參數(shù)化方案中(Hong et al.，2006)，湍流運動的1階非局地閉合產(chǎn)生如下的擴散方程

其中，C是預報變量，C=u，v，θ，q，qθ，qi，Kc是渦旋擴散系數(shù)，γc是對局地梯度的修正項(它代表大尺度渦旋對擴散的貢獻)，是ABL頂逆溫層的通量，h是ABL 高度，z是垂直坐標(從地面開始的高度)。動量渦旋擴散系數(shù)Km(C=u，v時的Kc)取如下形式(Hong et al.，2006)

其中，p是廓線形狀指數(shù)(通常取p=2)，k是馮卡門常數(shù)(k=0.4)，ws是混合層速度尺度，h、z的意義同上。

非局地擴散方程(1)的求解需要事先給定Kc的值，因此需要確定ABL高度h的數(shù)值。在YSU參數(shù)化方案中，h的值根據(jù)最低模式層和其上面的層之間的整體理Richardson 確定為第一個中性層的高度，具體來說是求解下列方程確定的

其中，Ribcr是臨界整體Richardson 數(shù)(取Ribcr=0.5)，U(h)是高度h處的水平風速，θva和θv(h)分別是最低模式層和高度h處的虛位溫，θs是地表附近的合適的溫度值。

在ACM2 參數(shù)化方案中，同樣需要確定ABL 高度的數(shù)值，其方法與YSU的類似。在MYJ、QNSE、Boulac參數(shù)化方案中，ABL高度h的值確定為湍動能TKE達到一個充分小的值(在WRF模式中，這個值取為0.005 m2·s-2)的高度。

對于不同的ABL過程參數(shù)化方案，其ABL內(nèi)的混合方式、夾卷方式以及ABL高度的確定方式等都不一定相同，因此一些學者對ABL參數(shù)化方案的性能進行了評估研究。Hong 和Pan(1996)研究表明，中尺度數(shù)值天氣預報模式在預報降水時，預報能力與非局地的垂直混合方案和臨界Richardson數(shù)密切相關(guān)；Wisse和Vila-Guerau(2004)指出，在降水模擬過程中不同ABL參數(shù)化方案產(chǎn)生的模擬結(jié)果各有特點，并且對流性降水模擬結(jié)果對ABL 參數(shù)化方案十分敏感。Shin 和Hong (2011)通過利用WRF 模式對CASES-99 (the Cooperative Atmosphere-Surface Exchange Study)外場觀測實驗中的某一天的天氣狀態(tài)進行數(shù)值模擬，檢驗了上述五個ABL 參數(shù)化方案的性能。他們除了比較不同的ABL參數(shù)化方案對ABL結(jié)構(gòu)的影響之外，還比較了參數(shù)化方案對ABL 高度的影響，結(jié)果表明，無論白天還是晚上，不同參數(shù)化方案所得到的ABL高度的差別都很大，QNSE參數(shù)化方案給出的ABL最深，而Bou-Lac給出的ABL最淺。Weisman等(2008)、Xie等(2012)的工作表明，局地閉合的ABL參數(shù)化方案診斷的ABL高度比非局地閉合的方案診斷的ABL高度要低。

LeMone等(2013，2014)專門利用ABL高度信息來評估ABL 參數(shù)化方案的性能，其利用CASES-97 (the April-May 1997 Kansas-based Cooperative Atmosphere-Surface Exchange Study) 外場觀測實驗期間獲得的白天和夜晚的觀測數(shù)據(jù)首先確定ABL高度(作為觀測的ABL 高度)，然后利用這些ABL 高度值檢驗WRF 模式中的四種ABL 參數(shù)化方案(YSU、MYJ、Bou-Lac、QNSE)對觀測的ABL高度的模擬能力。對于白天的對流ABL，LeMone等(2013，2014)直接把模式模擬的ABL 高度(即參數(shù)化方案給出的ABL 高度)與觀測的ABL 高度進行比較，結(jié)果表明BouLac、MYJ、QNSE 參數(shù)化方案所得到的ABL 高度通常低于觀測值，而YSU參數(shù)化方案得到的ABL 高度略高于觀測值；對于夜晚的穩(wěn)定ABL，未直接比較模式模擬的ABL 高度和觀測的ABL 高度的偏差，而是比較它們在整個夜間的變化趨勢，結(jié)果表明YSU和MYJ對于穩(wěn)定ABL的高度模擬具有較高的準確度，而BouLac和QNSE的模擬能力則較差，需要進行改進。

可見，ABL高度信息不僅在ABL參數(shù)化中是必不可少的參量，同時也可用來檢驗ABL參數(shù)化方案的性能。由于不同的參數(shù)化方案各有優(yōu)缺點，因此如何進一步改進完善現(xiàn)有ABL參數(shù)化方案仍然是一個挑戰(zhàn)。

7 結(jié)論與討論

本文從不同角度介紹ABL高度的確定及其應(yīng)用，提出在同時具有多種觀測資料時如何處理ABL 高度的三個基本原則，并WRF 模式為例重點介紹ABL 高度數(shù)據(jù)在數(shù)值天氣(或氣候)預報模式中的應(yīng)用，主要結(jié)論如下：

(1)ABL是自由大氣與地球表面(陸地和海洋)的連接通道，承擔著兩者之間動量、熱量、水汽和氣溶膠(污染物)等化學成分的交換任務(wù)，在地球氣候系統(tǒng)中具有關(guān)鍵的作用。ABL高度是ABL的一個重要參數(shù)，數(shù)值天氣(或氣候)模式中的ABL過程的參數(shù)化均需要用到它，低層云的范圍、種類及其演化與ABL 高度也有關(guān)系，準確獲取ABL高度數(shù)據(jù)具有重要的意義。

(2)對估計ABL 高度的主要方法及進展情況進行了簡要介紹。確定ABL高度的方法包括氣塊法、最大梯度法、湍流法、Richardson 數(shù)法等，這些方法本質(zhì)上都是最大梯度方法，涉及到的數(shù)據(jù)主要有常規(guī)探空溫濕廓線數(shù)據(jù)、衛(wèi)星遙感反演溫濕廓線數(shù)據(jù)、氣溶膠濃度廓線數(shù)據(jù)、GNSS 掩星探測廓線數(shù)據(jù)等。關(guān)于ABL高度的研究，研究工作主要集中在兩方面：一是從觀測數(shù)據(jù)出發(fā)估計ABL高度并進行精度驗證；另一方面是結(jié)合ABL參數(shù)化方案應(yīng)用于數(shù)值天氣(或氣候)模式中，檢驗ABL 參數(shù)化方案對暴雨、臺風等高影響天氣過程模擬的影響以及數(shù)值模式對ABL 高度的診斷能力，以評估完善數(shù)值天氣(或氣候)模式和ABL 參數(shù)化方案。

影響ABL高度時空分布的主要因素有兩個：湍流運動和垂直運動(Medeiros et al.，2005；韓美等，2012)；湍流運動與動力作用、熱力作用及云頂輻射冷卻有關(guān)，而垂直運動則與天氣系統(tǒng)或大尺度背景環(huán)流有關(guān)。影響ABL 高度量值大小及時空分布的具體因素包括：大氣穩(wěn)定度、下墊面狀況、太陽輻射、氣溶膠粒子濃度、積云和層積云、大尺度天氣系統(tǒng)等等(廖國蓮，2005)。由于ABL 高度涉及的因素多、問題復雜，關(guān)于ABL 高度確定及相關(guān)的問題仍然是一個前沿性的問題。近年來隨著人工智能技術(shù)的發(fā)展，機器學習技術(shù)廣泛應(yīng)用于各個領(lǐng)域，機器學習技術(shù)為人們研究ABL問題提供了新的手段和途徑(Rieutord et al.，2021；Krishnamurthy et al.，2021)，有望解決ABL相關(guān)科學問題。

致謝：陳飛教授對ABL 高度數(shù)據(jù)在WRF 模式中的應(yīng)用提供了指導，謹致謝忱！