劉一瑋，王 彥，王慶元

（1.天津市氣象臺，天津 300074；2.天津市人工影響天氣辦公室，天津 300074）

隨著氣象科技的發(fā)展， 天氣預(yù)報的準(zhǔn)確性和實效性不斷提高，并通過發(fā)展數(shù)值預(yù)報模式，逐步實現(xiàn)了從人工分析到定量化、客觀化、自動化、精細(xì)化的轉(zhuǎn)變[1]。 現(xiàn)代氣象預(yù)報業(yè)務(wù)的發(fā)展以無縫隙、全覆蓋、精準(zhǔn)化、智慧型為主要特征，高分辨率的格點(diǎn)要素預(yù)報是中國氣象局主推的天氣預(yù)報業(yè)務(wù)。 隨著數(shù)值預(yù)報技術(shù)的發(fā)展， 越來越豐富的數(shù)值模式產(chǎn)品可以為預(yù)報員提供多種信息， 精細(xì)化的氣象預(yù)報需要預(yù)報員從海量的預(yù)報產(chǎn)品中迅速提取有效信息，并依靠對天氣概念模型的構(gòu)建、 模式預(yù)報的理解進(jìn)行可能的主觀集成和訂正， 進(jìn)而完成預(yù)報產(chǎn)品的加工制作[2]。 因此要兼顧預(yù)報質(zhì)量和工作效率，就需要對高分辨率模式產(chǎn)品的釋用技術(shù)及不同要素之間的協(xié)同技術(shù)進(jìn)行研發(fā)和反復(fù)驗證。

2016 年中國氣象局印發(fā)了“全國精細(xì)化氣象格點(diǎn)預(yù)報業(yè)務(wù)建設(shè)實施方案”，天津市氣象局作為試點(diǎn)單位，目前已經(jīng)開展空間分辨率為1 km、時間分辨率為1 h、預(yù)報時效為10 d 的精細(xì)化格點(diǎn)預(yù)報業(yè)務(wù)，包含降水量、氣溫、風(fēng)向、風(fēng)速、相對濕度、云量、能見度等基礎(chǔ)要素和雷暴、冰雹、霧、沙塵等災(zāi)害性天氣共13 類產(chǎn)品。 從業(yè)務(wù)運(yùn)行現(xiàn)狀來看，存在預(yù)報產(chǎn)品之間不協(xié)調(diào)的問題， 主要原因是相關(guān)模式釋用技術(shù)不足，氣象要素之間沒有建立相互關(guān)聯(lián)[3]。

目前，圍繞模式釋用技術(shù)的研究已逐步開展，在模式溫度、降水、風(fēng)等要素的預(yù)報和檢驗等方面取得不少成果[4-11]。 李佰平和智協(xié)飛[13]采用均方根誤差對歐洲中心（ECWMF）確定性預(yù)報模式進(jìn)行評估，并分別利用一元線性回歸、多元線性回歸、單時效消除偏差和多時效消除偏差平均的訂正方法， 對ECMWF模式地面氣溫預(yù)報進(jìn)行訂正， 結(jié)果表明4 種訂正方法都能有效地減小地面氣溫多個時效預(yù)報的誤差，改進(jìn)幅度約為1 ℃。王婧等[14]采用平均法、雙權(quán)重平均法、 滑動平均法和滑動雙權(quán)重平均法分別對GRAPES_RAFS 系統(tǒng)2 m 溫度預(yù)報產(chǎn)品進(jìn)行偏差訂正，訂正后2m 溫度的平均誤差大部地區(qū)減小到-1～1 ℃。吳振玲等[15]、吳乃庚等[16]開展了多模式氣溫集成預(yù)報方法的研究，方法對最高、最低氣溫的主觀預(yù)報有較明顯的訂正能力。 吳啟樹等[17]應(yīng)用準(zhǔn)對稱混合滑動訓(xùn)練期，設(shè)計了基于TS 評分、最優(yōu)化訂正的OTS算法及基于ETS 評分最優(yōu)化訂正的OETS 算法，以確定預(yù)報日各級降水訂正系數(shù)。 這些方法現(xiàn)在及未來應(yīng)用于格點(diǎn)預(yù)報業(yè)務(wù)， 為格點(diǎn)預(yù)報的智能訂正提供科技支撐。

以上研究多是圍繞單一氣象要素的客觀預(yù)報方法而開展， 但針對多種氣象要素的關(guān)聯(lián)協(xié)調(diào)技術(shù)相關(guān)研究開展較少。 未來格點(diǎn)氣象要素間的協(xié)同一致性將是智能網(wǎng)格預(yù)報能力的重要評估手段之一，因此氣象要素間的協(xié)同訂正技術(shù)成為智能網(wǎng)格業(yè)務(wù)急需解決的問題。 本文在對歷史實況數(shù)據(jù)綜合統(tǒng)計分析的基礎(chǔ)上， 就降水與相對濕度的對應(yīng)關(guān)系進(jìn)行研究，得到相對濕度和降水之間的協(xié)同關(guān)系，構(gòu)建降水與相對濕度的協(xié)同方法， 以確保智能網(wǎng)格預(yù)報產(chǎn)品的協(xié)調(diào)性，進(jìn)而有效減輕當(dāng)前格點(diǎn)預(yù)報的工作強(qiáng)度，提升智能網(wǎng)格預(yù)報的科學(xué)性。

1 相對濕度觀測實況特征分析

1.1 資料與方法

在傳統(tǒng)的天氣預(yù)報中， 預(yù)報員只需要預(yù)報一個區(qū)域24 h 內(nèi)最大、最小相對濕度即可，而智能網(wǎng)格預(yù)報業(yè)務(wù)則需要對逐時的、 精細(xì)化網(wǎng)格點(diǎn)上的相對濕度做出差異化預(yù)報， 單憑預(yù)報員主觀訂正是難以實現(xiàn)的。一方面是預(yù)報制作時間有限，訂正精細(xì)化的產(chǎn)品費(fèi)時費(fèi)力； 更重要的是前期沒有相關(guān)的研究基礎(chǔ)，預(yù)報員對相對濕度的變化特征缺乏經(jīng)驗。因此研究協(xié)調(diào)訂正方法首先應(yīng)分析相對濕度的實況特征。

為確保計算效率及實況資料的可靠性，本文選取2016 年1 月—2017 年12 月天津市10 個國家氣象站逐小時地面觀測數(shù)據(jù)， 按照12—2 月為冬季、3—5月為春季、6—8 月為夏季、9—11 月為秋季進(jìn)行統(tǒng)計，得到相對濕度平均態(tài)及在降水期間的變化特征。

1.2 不同季節(jié)相對濕度逐時分布特征

4 個季節(jié)的平均相對濕度逐小時變化（圖1），其分布有明顯的季節(jié)特征和日變化特征。 春季的平均相對濕度最小，相對濕度在30%～62%；夏季最大，平均相對濕度在67%～84%；秋季和冬季的相對濕度居中； 相對濕度的這種季節(jié)差異主要受大氣環(huán)流和水汽的影響。相對濕度日變化特征突出，最小相對濕度出現(xiàn)在14—15 時，最大相對濕度出現(xiàn)在05—07 時，季節(jié)不同略有差異。 其中一天中最小相對濕度夏季出現(xiàn)的時間最晚，出現(xiàn)在15 時，其他季節(jié)出現(xiàn)在14時；最大相對濕度夏季則出現(xiàn)得最早，為05 時；冬季的最大相對濕度出現(xiàn)在07 時，在4 個季節(jié)中最晚。

圖1 不同季節(jié)逐時平均相對濕度

1.3 不同強(qiáng)度降水天氣的相對濕度分布特征

常規(guī)的預(yù)報訂正技術(shù)都是針對單一氣象要素進(jìn)行，當(dāng)某個氣象要素被訂正后，其他要素也應(yīng)該隨之改變。如某個格點(diǎn)上原先預(yù)報無降水，當(dāng)預(yù)報員訂正有降水后， 格點(diǎn)上相對濕度等其他氣象要素由于與降水有某種物理約束，其數(shù)值也應(yīng)該發(fā)生改變，從而達(dá)到與降水天氣現(xiàn)象一致。受傳統(tǒng)預(yù)報的影響，預(yù)報員對降水、氣溫等氣象要素的訂正能力較強(qiáng)，而對于相對濕度缺乏經(jīng)驗。

現(xiàn)行業(yè)務(wù)中網(wǎng)格預(yù)報的時間分辨率為1 h，本文將圍繞逐小時降水量和相對濕度展開討論。 首先根據(jù)小時降水量，將降水分為小雨（0.1～1.6 mm）、中雨（1.7～7.0 mm）、大雨（7.1～15.0 mm）、暴雨（＞15.0 mm），并統(tǒng)計不同季節(jié)不同量級降水的相對濕度分布情況（圖2）。

圖2 不同季節(jié)小雨（a）、中雨（b）、大雨（c）、暴雨（d）的相對濕度

當(dāng)出現(xiàn)小雨天氣時， 冬季的相對濕度為39%～99%，春季為23%～99%，夏季為37%～100%，秋季為52%～100%。 4 個季節(jié)相對濕度的離散度較大，原因可能是觀測數(shù)據(jù)中的降水量為1 h 累積量， 而相對濕度是觀測時刻的瞬時量，因此當(dāng)降水量較小時，就容易出現(xiàn)較小的相對濕度， 導(dǎo)致整體離散度較大的情況。冬季的相對濕度最大，秋季其次，春季最小。中雨天氣的相對濕度冬季為90%～99%，春季為69%～99%，夏季為51%～100%，秋季為88%～99%。相對濕度的平均數(shù)值明顯大于小雨天氣， 同時離散度也小于小雨天氣， 分析原因為降水是空氣中水汽壓超過了飽和水汽壓，導(dǎo)致水汽凝結(jié)析出，因此降水量對相對濕度影響明顯，隨著降水量增加，相對濕度相應(yīng)加大。 季節(jié)差異是冬季的相對濕度最大，秋季其次，春季最小。由于冬季的大雨天氣樣本不足，春季和冬季無暴雨天氣樣本，因此不做分析。其他季節(jié)的大雨及暴雨出現(xiàn)時相對濕度的特征類似，其數(shù)值均在73%以上，各季節(jié)之間的差異不大。

當(dāng)對格點(diǎn)上的相對濕度進(jìn)行訂正預(yù)報時， 首先應(yīng)考慮不同量級降水對應(yīng)的相對濕度的差異。

根據(jù)前期的統(tǒng)計結(jié)果， 相對濕度具有日變化的特征，且與降水密切相關(guān)。 如果已經(jīng)產(chǎn)生了降水，相對濕度是否還具有日變化的特征？考慮到樣本數(shù)量，圖3 給出了夏季不同等級降水時相對濕度的逐時變化特征。其中小雨天氣的相對濕度呈現(xiàn)出白天降低、夜晚升高的趨勢，與平均態(tài)類似，但日變化幅度明顯小于平均態(tài)。 11—16 時相對濕度較小，為85%～86%；01—05 時相對濕度較大，為90%～91%；晝夜差異在5%左右。 中雨天氣最小相對濕度出現(xiàn)在14 時，為88%；最大出現(xiàn)在07 時，為94%；日變化特征不明顯。 大雨以上降水的相對濕度日變化特征同樣不明顯。因此做降水與相對濕度協(xié)同訂正時，當(dāng)預(yù)報的降水為小雨或中雨時，應(yīng)考慮日變化，但變化幅度較小，最大相對濕度和最小相對濕度差值為5%～10%；大雨以上量級的降水時，相對濕度可以不考慮日變化。

圖3 夏季不同等級降水時相對濕度的逐時變化

1.4 降水時間、雨強(qiáng)與相對濕度對應(yīng)關(guān)系

圖4a 為降水持續(xù)時間與相對濕度的對應(yīng)關(guān)系。隨著降水時間的持續(xù)增加，相對濕度隨之增大；持續(xù)1 h 降水對應(yīng)的平均相對濕度是85%， 之后降水持續(xù)時間每增加1 h，平均相對濕度增加1.5%～2.0%；當(dāng)降水持續(xù)時間增加5～10 h，平均相對濕度變化不大，穩(wěn)定在91%～92%；當(dāng)持續(xù)時間在11 h 以上時，相對濕度會快速增大到95%； 持續(xù)12 h 以上的降水，相對濕度在97%以上。

選取降水持續(xù)時間在5 h 以上的個例， 統(tǒng)計了累積降水量與相對濕度的對應(yīng)關(guān)系， 發(fā)現(xiàn)降水量級與相對濕度也有很好的對應(yīng)（圖4b）。隨著累積降水量的增大， 相對濕度波動增加；1 mm 以下的降水平均相對濕度最小，為85%；降水量每增加1 mm，平均相對濕度增加0.5%～2.0%； 累積量超過5 mm 的降水相對濕度維持在90%以上； 降水達(dá)到25 mm后，相對濕度快速增加，最終穩(wěn)定在96%左右。

圖4 平均相對濕度與降水持續(xù)時間（a）、累積降水量（b）之間的關(guān)系

當(dāng)出現(xiàn)降水時，相對濕度有一個變化的過程，表1 統(tǒng)計了降水前后相對濕度變化的情況， 其中變量為降水時相對濕度與降水前3 h 相對濕度之差，變率表示為相對濕度之差與降水時相對濕度之比。 可以看出，降水時的相對濕度較降水前均有明顯增加；小雨天氣時降水前后相對濕度變化最小， 平均為9%；降水越強(qiáng)，相對濕度增幅越大；當(dāng)達(dá)到暴雨級別時，相對濕度的數(shù)值變化為21%。 變率也有相似的特征，4 種類型的降水相對濕度的變率均在10%以上，且隨著降水的增加而增大，小雨為11%，中雨為15%，大雨為17%，暴雨為22%。

表1 降水前后相對濕度的變化%

2 相對濕度歐洲數(shù)值模式特征分析

2.1 資料與方法

不同的預(yù)報要素間存在一定關(guān)系， 當(dāng)某個要素發(fā)生改變時， 其他相關(guān)要素也需要按照一定規(guī)則改變，以保證數(shù)據(jù)的一致性[3]。 本文使用歐洲中心數(shù)值模式數(shù)據(jù)（降水量和2 m 相對濕度）進(jìn)行特征統(tǒng)計分析。 其中數(shù)據(jù)的水平分辨率為0.125°×0.125°，時間分辨率為3 h， 研究區(qū)域為全國智能網(wǎng)格對天津責(zé)任區(qū)（38.40°～40.35°N，116.60°～118.20°E）。

2.2 不同強(qiáng)度降水天氣的相對濕度分布特征

對模式數(shù)據(jù)進(jìn)行分級統(tǒng)計，結(jié)果如圖5。小雨天氣的相對濕度冬季為29%～95%，春季為29%～95%，夏季為34%～99%，秋季為52%～99%。季節(jié)差異為秋季的相對濕度最大，夏季其次，春季、冬季最?。恢杏晏鞖獾南鄬穸榷緸?5%～99%， 春季為35%～95%，夏季為40%～100%，秋季為57%～98%；大雨天氣冬季無樣本，暴雨天氣冬季、春季無樣本，不做分析；大雨天氣的相對濕度春季為66%～92%，夏季為63%～99%，秋季為82%～98%；暴雨天氣的相對濕度夏季為92%～93%，秋季為74%～99%。 與實況類似，即降水量越小，相對濕度分布的離散度較大，但其離散程度遠(yuǎn)大于實況結(jié)果。 分析原因可能是數(shù)值預(yù)報是使用大氣運(yùn)動方程建立的數(shù)值模式框架， 通過同一積分步長計算出的各種氣象要素之間必然會保持協(xié)同一致性， 由于計算的結(jié)果與最后預(yù)報要素輸出的時間分辨率不一致， 會產(chǎn)生要素間協(xié)同不一致性問題。 涉及本文，模式輸出變量的時間分辨率為3 h，降水要素是個累加量，而降水持續(xù)時間未知，但是相對濕度輸出的則為瞬時量， 這種累加值和瞬時值可能出現(xiàn)協(xié)同不一致的情況， 因此在預(yù)報時不能直接采用數(shù)值模式的結(jié)果。

圖5 歐洲數(shù)值模式統(tǒng)計的不同季節(jié)小雨（a）、中雨（b）、大雨（c）、暴雨（d）的相對濕度

3 協(xié)同訂正算法檢驗評估

在實際格點(diǎn)預(yù)報工作中， 當(dāng)預(yù)報員應(yīng)用歐洲中心數(shù)值模式降水產(chǎn)品時， 還需對降水時刻的相對濕度產(chǎn)品進(jìn)行協(xié)同訂正。 基于偏差訂正提出降水和相對濕度協(xié)同訂正方法，公式如下：

其中，R 為訂正后的相對濕度；RHmin、RHmax分別為觀測得到的不同季節(jié)、不同量級降水的最小、最大相對濕度（表2，考慮離散度問題，剔除掉10%的最小相對濕度）；Rmin、Rmax為不同級別降水量的最小值、最大值；Xrain為預(yù)報降水量；P 為區(qū)間比例。

對2017 年1—12 月歐洲中心數(shù)值模式的相對濕度預(yù)報產(chǎn)品進(jìn)行協(xié)同計算， 并對訂正后的結(jié)果進(jìn)行檢驗。以平均絕對誤差、平均誤差和均方根誤差作為衡量協(xié)同訂正效果的指標(biāo)，計算結(jié)果顯示小雨、中雨、 大雨、 暴雨時分別平均絕對誤差分別為7.46、10.50、6.80、3.35，平均誤差為5.02、-7.63、-2.87、-1.60，均方根誤差為10.21、12.55、7.89、4.07。 說明利用協(xié)同訂正算法將增加小雨過程偏小的相對濕度，而對于中雨以上訂正誤差偏大。

在實際業(yè)務(wù)應(yīng)用中， 對于降水持續(xù)時間較短的小雨和中雨天氣過程可直接用上述方法協(xié)同； 而對于大雨以上量級的降水或持續(xù)時間超過5 h 的降水過程，可直接將相對濕度協(xié)同到95%以上。 從檢驗結(jié)果看， 使用此種降水和相對濕度的協(xié)同方法比數(shù)值模式直接輸出有一定的提高， 但由于觀測數(shù)據(jù)樣本和模式數(shù)據(jù)的時間分辨率的問題， 協(xié)同統(tǒng)計算法仍有優(yōu)化空間。

4 結(jié)論與討論

統(tǒng)計分析了實況降水與相對濕度的一致性特征關(guān)系， 并檢驗評估了歐洲中心數(shù)值模式預(yù)報的降水與相對濕度的一致性偏差，并基于分析和評估結(jié)果，設(shè)計了降水與相對濕度協(xié)同訂正技術(shù)方法， 可實現(xiàn)降水訂正與相對濕度氣象要素預(yù)報產(chǎn)品協(xié)同輸出，對進(jìn)一步科學(xué)地完善智能網(wǎng)格業(yè)務(wù)具有重要的現(xiàn)實意義。 通過分析得到以下結(jié)論：

（1）實況觀測降水量級越小，其相對濕度的離散度越大，且日變化特征明顯。 降水量越大、降水持續(xù)時間越長，相對濕度越穩(wěn)定，小時降水量在5 mm 以上且持續(xù)時間超過5 h，相對濕度維持在90%以上。

（2）歐洲中心數(shù)值模式預(yù)報的降水和相對濕度由于累計量與瞬時量的輸出時間偏差， 使得降水時對應(yīng)的相對濕度離散度較大，與實況存在偏差，需對降水時的相對濕度進(jìn)行協(xié)同訂正。

（3）基于偏差訂正提出降水和相對濕度協(xié)同訂正算法，并進(jìn)行檢驗評估。結(jié)果表明，訂正后，小雨過程偏小的相對濕度明顯增加，效果較好；對于中雨天氣，訂正后相對濕度較實況偏低，但普遍優(yōu)于數(shù)值預(yù)報輸出產(chǎn)品。 未來可在增加實況統(tǒng)計樣本和提高數(shù)值模式產(chǎn)品時間分辨率的基礎(chǔ)上， 進(jìn)一步優(yōu)化協(xié)同算法。