徐記亮 張偉 蘇艷華 顧雷

(中國(guó)民用航空局空中交通管理局航空氣象中心，北京 100018)

引言

低云是危及飛行安全的危險(xiǎn)性天氣之一。在低云覆蓋機(jī)場(chǎng)的情況下著陸，如果飛機(jī)出云后離地面的高度較低，且又未對(duì)準(zhǔn)跑道，往往來(lái)不及修正，容易造成復(fù)飛。有時(shí)，由于指揮或操作不當(dāng)，還可能造成飛機(jī)與地面障礙物相撞或失速的事故[1]。因此，云的觀測(cè)預(yù)報(bào)是機(jī)場(chǎng)天氣報(bào)告中必須包含的氣象要素，云高、云量更是關(guān)注的重點(diǎn)。

目前，云觀測(cè)主要為地面人工觀測(cè)、飛機(jī)探測(cè)和衛(wèi)星遙感等，每種觀測(cè)方式只能反映云的一部分信息[2]。綜合利用各類(lèi)觀測(cè)及模式輸出結(jié)果的云分析方法，可以獲取云量、云狀、云高等信息，是云預(yù)報(bào)和分析的基礎(chǔ)[3]。其中，探空識(shí)別法可以利用相對(duì)濕度確定云的垂直分層。Wang和Rossow[4]利用相對(duì)濕度廓線提出具體的閾值標(biāo)準(zhǔn)，稱為WR95算法。趙仕偉等[5]將WR95算法運(yùn)用到數(shù)值模式中，以模式輸出的數(shù)據(jù)代替探空數(shù)據(jù)判定云底高，結(jié)果表明該方法對(duì)中低云較為適用。蔡淼等[6]在WR95基礎(chǔ)上對(duì)相對(duì)濕度閾值進(jìn)行了優(yōu)化，優(yōu)化后的準(zhǔn)確率和TS評(píng)分略有提高。以上探空識(shí)別方法可以識(shí)別云的垂直分層，但無(wú)法確定云量。衛(wèi)星遙感可以檢測(cè)云并通過(guò)反演獲取云量[7- 9]。但衛(wèi)星產(chǎn)品中的云量是天空中的總云量，與地面觀測(cè)總云量具有不同的物理意義[8]，且與地面觀測(cè)云量存在一定偏差[9]，就地面觀測(cè)點(diǎn)的總云量而言，地面觀測(cè)值的可靠性相對(duì)較高[8]。在開(kāi)展的云量診斷預(yù)報(bào)研究中[10-14]， Xu和Randall[15]基于CRM模擬結(jié)果提出云量診斷方法，此云量診斷方法被應(yīng)用于NCEP GFS模式中計(jì)算云量[16-17]。

目前，中國(guó)氣象局建立了三維云融合分析業(yè)務(wù)系統(tǒng)[18]，生成云底、云頂、云量等三維產(chǎn)品，用于天氣分析。但機(jī)場(chǎng)云高、云量的預(yù)報(bào)還不能從數(shù)值模式中直接獲得，更多地依靠經(jīng)驗(yàn)進(jìn)行預(yù)報(bào)[19]，部分機(jī)場(chǎng)建立了MOS系統(tǒng)預(yù)報(bào)云量[20]，但仍無(wú)法達(dá)到定量計(jì)算、準(zhǔn)確預(yù)報(bào)的水平，且針對(duì)機(jī)場(chǎng)云的預(yù)報(bào)診斷研究仍較少。

根據(jù)2019年中國(guó)大陸地區(qū)機(jī)場(chǎng)航班量統(tǒng)計(jì)，最繁忙的前十大機(jī)場(chǎng)分別為北京首都、上海浦東、廣州白云、成都雙流、深圳寶安、昆明長(zhǎng)水、西安咸陽(yáng)、上海虹橋、重慶江北和杭州蕭山機(jī)場(chǎng)。十大機(jī)場(chǎng)的云診斷預(yù)報(bào)分析可增強(qiáng)民航氣象部門(mén)對(duì)航班運(yùn)行的保障能力，提高航班正常率和航空服務(wù)質(zhì)量，對(duì)于民航部門(mén)來(lái)說(shuō)非常重要。本文基于數(shù)值模式的云高、云量診斷算法，對(duì)比分析云診斷方法在中國(guó)十大機(jī)場(chǎng)的應(yīng)用情況，以探討適用于某一機(jī)場(chǎng)的云診斷算法，為中國(guó)機(jī)場(chǎng)的云預(yù)報(bào)提供參考。

1 資料與方法

1.1 資料來(lái)源

民航機(jī)場(chǎng)云的觀測(cè)以人工目測(cè)方式進(jìn)行，云量采用八分量制(中國(guó)民用航空局空管行業(yè)管理辦公室.《民用航空飛行氣象情報(bào)發(fā)布與交換辦法》，2012)：1—2個(gè)八分量為少云(FEW)，3—4個(gè)八分量為疏云(SCT)，5—7個(gè)八分量為多云(BKN)，8個(gè)八分量為滿天云(OVC)。民航機(jī)場(chǎng)氣象臺(tái)按規(guī)定每日24 h于整點(diǎn)或半點(diǎn)連續(xù)發(fā)布機(jī)場(chǎng)例行天氣報(bào)告(METAR報(bào))(中國(guó)民用航空局空管行業(yè)管理辦公室，中國(guó)民用航空局空中交通管理局.《民用航空飛行氣象情報(bào)發(fā)布與交換辦法》，2009)。METAR報(bào)中，云高以30 m為增量等級(jí)編報(bào)。

民航運(yùn)行中，對(duì)飛行有重要影響的云是指云底高為1500 m的云，或最高的最低扇區(qū)高度(兩者取其大)以下的云，或任何高度的積雨云、濃積云。另外，距離機(jī)場(chǎng)更近的最低層的云層或云塊對(duì)飛機(jī)的正常起降影響更大。

選取2016年5月至2019年4月共3 a的NCEP再分析資料FNL(0.25°×0.25°)世界時(shí)每天00時(shí)、06時(shí)、12時(shí)、18時(shí)的6 h數(shù)值預(yù)報(bào)資料及相應(yīng)時(shí)次的METAR報(bào)資料，分析云底、云高在1500 m以下最低層云的情況。

1.2 研究方法

1.2.1 探空識(shí)別方法

1.2.1.1 云高計(jì)算

WR95算法[4-5]計(jì)算云底高方法：從地表向上，當(dāng)相對(duì)濕度(RH)≥87%時(shí)即為入云，以此判斷為第一個(gè)云層的云底。當(dāng)84%≤RH<87%且高度層RH的變化大于3%時(shí)，也判斷為第一個(gè)云層的云底。當(dāng)RH<84%即為出云，以此判斷為第一個(gè)云層的云頂。第N(N≥2)個(gè)云層的相對(duì)濕度閾值為84%，當(dāng)RH≥84%時(shí)即為入云，以此判斷為云底；當(dāng)RH<84%時(shí)即為出云，以此判斷為云頂。

優(yōu)化的WR95方法(簡(jiǎn)稱WR95opt)各高度相對(duì)濕度閾值與高度的關(guān)系見(jiàn)式(1)[6]。

(1)

式(1)中，Hr為相對(duì)濕度閾值；H為高度，單位為km。同時(shí)參考Zhang等[21]的方法，考慮了云層數(shù)量?jī)?yōu)化。

1.2.1.2 云高推算云量

以上探空識(shí)別方法可以識(shí)別云的垂直分層，但無(wú)法確定云的云量。本文設(shè)計(jì)的由云高推算云量的模型，具體算法如下：

在得出云底高h(yuǎn)后，首先計(jì)算在機(jī)場(chǎng)A處，觀測(cè)者觀測(cè)到高度為h的視場(chǎng)范圍，即半徑為r的圓形區(qū)域；再計(jì)算此圓形區(qū)域內(nèi)包含的網(wǎng)格點(diǎn)數(shù)，此網(wǎng)格點(diǎn)數(shù)同一高度h處有云的比例即為機(jī)場(chǎng)A處觀測(cè)到的云底高為h的云的云量。示意如圖1。

R為地球平均半徑；r為觀測(cè)視場(chǎng)半徑；h為云高度圖1 機(jī)場(chǎng)A處觀測(cè)到云的視場(chǎng)范圍示意Fig.1 Schematic diagram of the field of view of the cloud observed at the airport A

1.2.2 云量診斷方法

NCEP GFS模式中診斷云量的計(jì)算方法[15-16]見(jiàn)式(2)～式(4)，為C云量診斷法(簡(jiǎn)稱C方法)。

(2)

(3)

(4)

式(2)～式(4)中，C為云量(成);RH為相對(duì)濕度(%);q*為飽和比濕(kg/kg);qc為云水/云冰混合比(kg/kg)；e為水汽壓(Pa)；p為氣壓(Pa)；T為溫度(K)。

計(jì)算出各高度上的云量，可以推算不同云量的云底高。

1.2.3 診斷效果檢驗(yàn)

參考晴雨預(yù)報(bào)的混淆(誤差)矩陣分別計(jì)算十大機(jī)場(chǎng)診斷云的準(zhǔn)確率、漏報(bào)率及TS評(píng)分，見(jiàn)表1。

表1 云診斷的混淆(誤差)矩陣Table 1 Confusion (error) matrix for cloud diagnostics

準(zhǔn)確率=TP/(TP+FP)，漏報(bào)率=FN/(TP+FN)，TS評(píng)分=TP/(TP+FN+FP)。

準(zhǔn)確率是相對(duì)預(yù)報(bào)而言的，數(shù)值越大表示預(yù)報(bào)可信度越高。漏報(bào)率是相對(duì)實(shí)況而言的，數(shù)值越大表示出現(xiàn)漏報(bào)的可能性越大。TS評(píng)分綜合考慮了準(zhǔn)確率和漏報(bào)率，數(shù)值越大表示預(yù)報(bào)效果越好。云底高的診斷效果根據(jù)診斷的云的平均云底高與相應(yīng)時(shí)次機(jī)場(chǎng)觀測(cè)的平均云底高的相對(duì)誤差來(lái)檢驗(yàn)。

2 結(jié)果分析

2.1 不分云量

不分云量時(shí)， 2016年5月至2019年4月中國(guó)十大機(jī)場(chǎng)利用C、WR95、WR95opt方法計(jì)算的1500 m以下6 h預(yù)報(bào)出現(xiàn)最低層云的準(zhǔn)確率、漏報(bào)率及TS評(píng)分見(jiàn)表2。由表2可知，在廣州、成都、深圳、昆明、重慶、杭州機(jī)場(chǎng)，C方法云量計(jì)算結(jié)果的準(zhǔn)確率均較高，為0.60以上，其中昆明最高，為0.85；北京和西安機(jī)場(chǎng)較低，為0.30以下。成都、重慶、西安、杭州、昆明、廣州機(jī)場(chǎng)WR95和WR95opt方法云量計(jì)算結(jié)果的準(zhǔn)確率較高，為0.80以上，北京最低。所有機(jī)場(chǎng)WR95opt的準(zhǔn)確率均比WR95稍高。除昆明機(jī)場(chǎng)

表2 2016年5月至2019年4月中國(guó)十大機(jī)場(chǎng)1500 m以下不分云量時(shí)6 h預(yù)報(bào)出現(xiàn)最低層云的準(zhǔn)確率、漏報(bào)率、TS評(píng)分Table 2 The accuracy rate,miss rate,and TS scores of the lowest stratified cloud forecasted at 6 h below 1500 m regardless of cloud cover in the top ten airports in China form May of 2016 to April of 2019

的WR95和WR95opt方法準(zhǔn)確率比C方法稍低外，其他機(jī)場(chǎng)WR95和WR95opt方法的準(zhǔn)確率均比C方法高。

C方法云量計(jì)算的漏報(bào)率，昆明最高，為0.99，其他機(jī)場(chǎng)相對(duì)較低，為0.30以下。西安機(jī)場(chǎng)WR95和WR95opt方法的漏報(bào)率最高，為0.70以上，上海虹橋最低。所有機(jī)場(chǎng)WR95opt方法的漏報(bào)率均比WR95方法稍高。除昆明WR95和WR95opt方法計(jì)算的漏報(bào)率比C方法低以外，其他機(jī)場(chǎng)WR95和WR95opt方法計(jì)算的漏報(bào)率均比C方法高。

成都、杭州機(jī)場(chǎng)C方法云量計(jì)算結(jié)果的TS評(píng)分較高，分別為0.65和0.60；西安、北京、昆明機(jī)場(chǎng)較低，其值均小于0.30，其中昆明最低，只有0.01。廣州、上海虹橋、杭州機(jī)場(chǎng)WR95和WR95opt方法云量計(jì)算結(jié)果的TS評(píng)分較高，為0.60左右；西安最低，只有0.26。北京WR95方法云量計(jì)算結(jié)果TS評(píng)分比WR95opt方法稍低，其他機(jī)場(chǎng)的WR95方法比WR95opt方法稍高。北京、上海浦東、昆明、上海虹橋機(jī)場(chǎng)WR95和WR95opt方法云量計(jì)算結(jié)果的TS評(píng)分稍高于C方法；成都、西安、重慶機(jī)場(chǎng)的TS評(píng)分低于C方法；廣州、深圳、杭州機(jī)場(chǎng)C方法云量計(jì)算結(jié)果TS評(píng)分介于WR95方法與WR95opt方法之間。

值得注意的是，昆明機(jī)場(chǎng)C方法云量計(jì)算結(jié)果的漏報(bào)率為0.99、TS評(píng)分為0.01，明顯與其他機(jī)場(chǎng)及WR95和WR95opt方法差異較大，其結(jié)果基本均為無(wú)云的情況。這可能與昆明機(jī)場(chǎng)特殊地形有關(guān)，昆明機(jī)場(chǎng)選址處于山腰，海拔為2103.37 m，機(jī)場(chǎng)東西兩側(cè)為山嶺，南北為低洼地，位于一條東南—西北向的通道中，南部開(kāi)闊，北部稍窄，容易形成山谷風(fēng)。同時(shí)計(jì)算另外兩個(gè)高海拔機(jī)場(chǎng)——西寧曹家堡機(jī)場(chǎng)、拉薩貢嘎機(jī)場(chǎng)的準(zhǔn)確率、漏報(bào)率、TS評(píng)分，其值分別為N、1.00、0.00，即兩個(gè)機(jī)場(chǎng)均沒(méi)有診斷出云，診斷效果均極差。由此可以大致推斷，C方法在高原機(jī)場(chǎng)的診斷效果很差，原因可能是該方法受到高原等復(fù)雜地形或海拔高度的影響較大。

2.2 少云(FEW)

當(dāng)少云量時(shí)， 2016年5月至2019年4月中國(guó)十大機(jī)場(chǎng)利用C、WR95、WR95opt方法計(jì)算的1500 m以下6 h預(yù)報(bào)出現(xiàn)最低層云的準(zhǔn)確率、漏報(bào)率及TS評(píng)分見(jiàn)表3。由表3可知，各機(jī)場(chǎng)不同方法的計(jì)算結(jié)果均明顯變差。FEW少云時(shí)，各機(jī)場(chǎng)準(zhǔn)確率、TS評(píng)分明顯減小，漏報(bào)率增大，尤其WR95及WR95opt方法的漏報(bào)率明顯增大。

除杭州外，其他機(jī)場(chǎng)WR95和WR95opt方法云量計(jì)算結(jié)果的準(zhǔn)確率均比C方法高。成都、重慶三種方法云量計(jì)算結(jié)果的準(zhǔn)確率最高，北京、上海虹橋、杭州機(jī)場(chǎng)最低。

表3 2016年5月至2019年4月中國(guó)十大機(jī)場(chǎng)1500 m以下FEW少云時(shí)6 h預(yù)報(bào)出現(xiàn)最低層云的準(zhǔn)確率、漏報(bào)率、TS評(píng)分Table 3 The accuracy rate,miss rate,and TS scores of the lowest stratified cloud forecasted at 6 h below 1500 m when clouds are few in the top ten airports in China from May of 2016 to April of 2019

除昆明機(jī)場(chǎng)外，其他機(jī)場(chǎng)WR95和WR95opt方法云量計(jì)算結(jié)果的漏報(bào)率均高于C方法，為0.84以上。各機(jī)場(chǎng)WR95opt方法云量計(jì)算結(jié)果的漏報(bào)率稍低于WR95方法。昆明機(jī)場(chǎng)三種方法云量計(jì)算結(jié)果的漏報(bào)率最高，為0.96以上，尤其是C方法云量計(jì)算結(jié)果均高于其他機(jī)場(chǎng)。

北京、昆明、上海虹橋機(jī)場(chǎng)C方法云量計(jì)算結(jié)果的TS評(píng)分稍低于或與WR95和WR95opt方法相當(dāng)，其他機(jī)場(chǎng)C方法云量計(jì)算結(jié)果的TS評(píng)分高于WR95和WR95opt方法。WR95和WR95opt方法計(jì)算結(jié)果的TS評(píng)分均較小，為0.10以下。成都、重慶、深圳機(jī)場(chǎng)C方法云量計(jì)算結(jié)果的TS評(píng)分稍高于其他機(jī)場(chǎng)。

2.3 疏云(SCT)

疏云時(shí)， 2016年5月至2019年4月中國(guó)十大機(jī)場(chǎng)利用C、WR95、WR95opt方法計(jì)算的1500 m以下6 h預(yù)報(bào)出現(xiàn)最低層云的準(zhǔn)確率、漏報(bào)率及TS評(píng)分見(jiàn)表4。由表4可知，除杭州外，其他機(jī)場(chǎng)C方法計(jì)算結(jié)果的準(zhǔn)確率均低于WR95和WR95opt方法。北京首都機(jī)場(chǎng)C方法云量計(jì)算結(jié)果準(zhǔn)確率為0.00，或者診斷無(wú)SCT疏云。成都機(jī)場(chǎng)的WR95方法云量計(jì)算結(jié)果準(zhǔn)確率高于WR95opt方法。成都機(jī)場(chǎng)WR95方法云量計(jì)算結(jié)果準(zhǔn)確率最高，而昆明機(jī)場(chǎng)WR95opt方法最高。

C方法計(jì)算結(jié)果的漏報(bào)率均為1.00，WR95和WR95opt方法的漏報(bào)率也較高，均為0.79以上。北京機(jī)場(chǎng)WR95方法云量計(jì)算結(jié)果漏報(bào)率低于WR95opt方法，其他機(jī)場(chǎng)WR95方法云量計(jì)算結(jié)果漏報(bào)率高于WR95opt方法，或兩者相當(dāng)。

C方法云量計(jì)算結(jié)果的TS評(píng)分均為0。WR95方法云量計(jì)算結(jié)果的TS評(píng)分稍低于WR95opt方法。廣州機(jī)場(chǎng)WR95和WR95opt計(jì)算結(jié)果的TS評(píng)分最高，為0.14，其他機(jī)場(chǎng)均為0.10以下。

2.4 多云(BKN)

多云時(shí)， 2016年5月至2019年4月中國(guó)十大機(jī)場(chǎng)利用C、WR95、WR95opt方法計(jì)算的1500 m以下6 h預(yù)報(bào)出現(xiàn)最低層云的準(zhǔn)確率、漏報(bào)率及TS評(píng)分見(jiàn)表5。由表5可知，除上海虹橋、杭州機(jī)場(chǎng)外，其他機(jī)場(chǎng)C方法未診斷出BKN云量。西安WR95方法云量計(jì)算結(jié)果的準(zhǔn)確率高于WR95opt方法，其他機(jī)場(chǎng)兩種方法計(jì)算結(jié)果相當(dāng)，或WR95方法云量計(jì)算結(jié)果低于WR95opt方法。上海虹橋機(jī)場(chǎng)WR95和WR95opt方法云量計(jì)算結(jié)果的準(zhǔn)確率最高，分別為0.45和0.56。

C方法云量計(jì)算結(jié)果的漏報(bào)率均為1.00。WR95和WR95opt方法云量計(jì)算結(jié)果的漏報(bào)率，除深圳機(jī)場(chǎng)稍低外，其他機(jī)場(chǎng)均較高，為0.60以上。上海浦東、上海虹橋、杭州機(jī)場(chǎng)WR95方法云量計(jì)算結(jié)果的漏報(bào)率高于WR95opt方法，其他機(jī)場(chǎng)WR95方法計(jì)算結(jié)果低于WR95opt方法，或者二者相當(dāng)。

表4 2016年5月至2019年4月中國(guó)十大機(jī)場(chǎng)1500 m以下SCT疏云時(shí)6 h預(yù)報(bào)出現(xiàn)最低層云的準(zhǔn)確率、漏報(bào)率、TS評(píng)分Table 4 The accuracy rate,miss rate,and TS scores of the lowest stratified cloud forecasted at 6 h below 1500 m when clouds are sparse in the top ten airports in China from May of 2016 to April of 2019

表5 2016年5月至2019年4月中國(guó)十大機(jī)場(chǎng)1500 m以下BKN多云時(shí)6 h預(yù)報(bào)出現(xiàn)最低層云的準(zhǔn)確率、漏報(bào)率、TS評(píng)分Table 5 The accuracy rate,miss rate,and TS scores of the lowest stratified cloud forecasted at 6 h below 1500 m when clouds are broken in the top ten airports in China from May of 2016 to April of 2019

C方法云量計(jì)算結(jié)果的TS評(píng)分均為0.00。成都、西安機(jī)場(chǎng)WR95方法云量計(jì)算結(jié)果的TS評(píng)分稍高于WR95opt方法，其他機(jī)場(chǎng)兩者相等，或WR95方法計(jì)算結(jié)果稍低于WR95opt方法。上海浦東、深圳、上海虹橋、杭州機(jī)場(chǎng)WR95和WR95opt方法計(jì)算結(jié)果的TS評(píng)分高于其他機(jī)場(chǎng)，均為0.20以上；北京、成都、昆明、西安、重慶的TS評(píng)分均為0.10以下。

2.5 滿天云(OVC)

滿天云時(shí)，2016年5月至2019年4月中國(guó)十大機(jī)場(chǎng)利用C、WR95、WR95opt方法計(jì)算的1500 m以下6 h預(yù)報(bào)出現(xiàn)最低層云的準(zhǔn)確率、漏報(bào)率及TS評(píng)分見(jiàn)表6。由表6可知，各機(jī)場(chǎng)C方法云量計(jì)算結(jié)果的準(zhǔn)確率均為0.00，或者未診斷出OVC云量。西安、重慶機(jī)場(chǎng)WR95方法云量計(jì)算結(jié)果的準(zhǔn)確率高于WR95opt方法，其他機(jī)場(chǎng)兩者相當(dāng)，或WR95方法云量計(jì)算結(jié)果的準(zhǔn)確率低于WR95opt方法。廣州機(jī)場(chǎng)WR95和WR95opt方法計(jì)算結(jié)果的準(zhǔn)確率最高，其值分別為0.13和0.17，其次為北京機(jī)場(chǎng)，其他機(jī)場(chǎng)均低于0.10。

表6 2016年5月至2019年4月中國(guó)十大機(jī)場(chǎng)1500 m以下OVC滿天云時(shí)6 h預(yù)報(bào)出現(xiàn)最低層云的準(zhǔn)確率、漏報(bào)率、TS評(píng)分Table 6 The accuracy rate,miss rate,and TS scores of the lowest stratified cloud forecasted at 6 h below 1500 m when clouds are overcast in the top ten airports in China from May of 2016 to April of 2019

各機(jī)場(chǎng)C方法云量計(jì)算結(jié)果的漏報(bào)率均為1.00。成都、杭州機(jī)場(chǎng)WR95和WR95opt的漏報(bào)率均為1.00，其他機(jī)場(chǎng)WR95方法計(jì)算結(jié)果的漏報(bào)率低于WR95opt方法。

各機(jī)場(chǎng)C方法云量計(jì)算結(jié)果的TS評(píng)分均為0.00。各機(jī)場(chǎng)WR95和WR95opt方法云量計(jì)算結(jié)果的TS評(píng)分均為0.10以下，廣州機(jī)場(chǎng)最高，分別為0.08、0.07，成都、昆明、重慶、杭州機(jī)場(chǎng)均為0.00。北京、廣州、西安機(jī)場(chǎng)WR95方法計(jì)算結(jié)果的TS評(píng)分高于WR95opt方法，深圳機(jī)場(chǎng)低于WR95opt方法，其他機(jī)場(chǎng)兩者相當(dāng)。

2.6 診斷云高相對(duì)誤差

2016年5月至2019年4月中國(guó)十大機(jī)場(chǎng)各方法計(jì)算的1500 m以下最低層云平均云底高與METAR的相對(duì)誤差見(jiàn)圖2。由圖2可知，各方法診斷云高相對(duì)誤差差別較大，相對(duì)誤差有正有負(fù)。昆明機(jī)場(chǎng)C方法云高計(jì)算結(jié)果的相對(duì)誤差為正值，且誤差較大，為163%，表明昆明機(jī)場(chǎng)C方法云高的計(jì)算結(jié)果遠(yuǎn)高于METAR；其余機(jī)場(chǎng)的相對(duì)誤差均為負(fù)值，表明其云高計(jì)算結(jié)果低于METAR。其中重慶機(jī)場(chǎng)相對(duì)誤差最小，為-19%；杭州、深圳、上海虹橋、北京、上海浦東、成都機(jī)場(chǎng)較大，其相對(duì)誤差分別為-73%、-67%、-64%、-61%、-59%、-51%。

成都、西安、重慶機(jī)場(chǎng)WR95方法云高計(jì)算結(jié)果的相對(duì)誤差為正值，分別為13%、63%、51%，其他機(jī)場(chǎng)的相對(duì)誤差為負(fù)，其中昆明機(jī)場(chǎng)最大，為-93%，北京、廣州、深圳機(jī)場(chǎng)較小，其值分別為-17%、-17%、-12%。

正值為高于METAR，負(fù)值為低于METAR圖2 2016年5月至2019年4月中國(guó)十大機(jī)場(chǎng)1500 m以下最低層云平均云底高與METAR的相對(duì)誤差Fig.2 Relative error between METAR and average cloud base height of the lowest layer clouds below 1500 m in top ten airports in China from May of 2016 to April of 2019

上海浦東、昆明、上海虹橋機(jī)場(chǎng)WR95opt方法云高的計(jì)算結(jié)果相對(duì)誤差為負(fù)，其中上海浦東、上海虹橋機(jī)場(chǎng)較小，其值均為-10%，昆明機(jī)場(chǎng)較大，為-43%。其他機(jī)場(chǎng)的相對(duì)誤差為正，其中北京、廣州、深圳、杭州機(jī)場(chǎng)較小，分別為5%、8%、7%、0%；重慶、西安、成都機(jī)場(chǎng)較大，分別為111%、106%、43%。

從相對(duì)誤差絕對(duì)值來(lái)看，北京、上海浦東、廣州、成都、深圳、昆明、上海虹橋、杭州機(jī)場(chǎng)WR95和WR95opt方法云高的計(jì)算結(jié)果好于C方法，西安、重慶機(jī)場(chǎng)的C方法計(jì)算結(jié)果較好。WR95與WR95opt方法比較，北京、上海浦東、廣州、深圳、昆明、上海虹橋、杭州機(jī)場(chǎng)WR95opt方法的云高計(jì)算結(jié)果好于WR95方法，其相對(duì)誤差絕對(duì)值均在10%或以下，其中杭州為0%；成都、西安、重慶機(jī)場(chǎng)WR95方法較好，成都機(jī)場(chǎng)的相對(duì)誤差絕對(duì)值為13%。西安、重慶機(jī)場(chǎng)的C方法云高計(jì)算結(jié)果較好，其相對(duì)誤差的絕對(duì)值相對(duì)較大，分別為31%、19%。

3 結(jié)論與討論

(1)不分云量時(shí)，除昆明機(jī)場(chǎng)外，其他機(jī)場(chǎng)C方法云量計(jì)算結(jié)果的準(zhǔn)確率低于WR95和WR95opt方法。除昆明、廣州機(jī)場(chǎng)外，其他機(jī)場(chǎng)C方法云量計(jì)算結(jié)果的漏報(bào)率均低于WR95和WR95opt方法。北京、上海浦東、昆明、上海虹橋機(jī)場(chǎng)WR95和WR95opt方法云量計(jì)算結(jié)果的TS評(píng)分稍高于C方法，成都、西安、重慶機(jī)場(chǎng)低于C方法。所有機(jī)場(chǎng)WR95opt方法云量計(jì)算結(jié)果的準(zhǔn)確率、漏報(bào)率均稍高于WR95方法；除北京機(jī)場(chǎng)外，其他機(jī)場(chǎng)WR95方法云量計(jì)算結(jié)果的TS評(píng)分均稍高于WR95opt方法。

(2)分云量時(shí)，各方法的診斷效果明顯變差，準(zhǔn)確率、TS評(píng)分明顯減小，漏報(bào)率明顯增大。各機(jī)場(chǎng)C方法在SCT、BKN、OVC云量時(shí)的云量計(jì)算結(jié)果準(zhǔn)確率基本為0.00，或者診斷無(wú)此類(lèi)云量、漏報(bào)率為1.00、TS評(píng)分均為0.00，其對(duì)于SCT、BKN、OVC云量的診斷效果比FEW云量時(shí)差。大部分機(jī)場(chǎng)C方法云量計(jì)算結(jié)果的準(zhǔn)確率低于WR95和WR95opt方法。FEW時(shí)，C方法云量計(jì)算結(jié)果的TS評(píng)分高于WR95和WR95opt方法，SCT、BKN、OVC時(shí)C方法云量計(jì)算結(jié)果的TS評(píng)分低于WR95和WR95opt方法。

(3)云高相對(duì)誤差有正有負(fù)，昆明機(jī)場(chǎng)C方法云高計(jì)算結(jié)果的相對(duì)誤差為正，其余機(jī)場(chǎng)為負(fù)。成都、西安、重慶機(jī)場(chǎng)WR95方法云高計(jì)算結(jié)果的相對(duì)誤差為正，其余機(jī)場(chǎng)為負(fù)；上海浦東、昆明、上海虹橋機(jī)場(chǎng)WR95opt方法云高計(jì)算結(jié)果的相對(duì)誤差為負(fù)，其余機(jī)場(chǎng)為正。從相對(duì)誤差的絕對(duì)值來(lái)看，北京、上海浦東、廣州、深圳、上海虹橋、杭州、昆明的WR95opt方法診斷效果較好；成都的WR95方法診斷效果較好；西安、重慶的C方法診斷效果較好。

(4)C方法對(duì)云的診斷效果不如探空法，當(dāng)SCT及以上云量時(shí)，大部分沒(méi)有診斷出云量，原因可能是C方法更傾向于診斷少云，可以嘗試通過(guò)提高相對(duì)濕度的比重或加入某一系統(tǒng)偏差項(xiàng)，以提高多云量的診斷效果。大部分機(jī)場(chǎng)的云高相對(duì)誤差為負(fù)值且其絕對(duì)值大于探空法，表明C方法診斷的云高遠(yuǎn)低于實(shí)際，原因可能是C方法能夠診斷出云的“閾值”較小，實(shí)際云層以下容易達(dá)到C方法的“閾值”而被誤診斷為云，可以通過(guò)嘗試提高云的“閾值”以提高診斷的云底高，減小云高相對(duì)誤差，需進(jìn)一步的對(duì)比分析。