鄒雙澤 白愛娟 何科 金海東 黃金穎 金波

摘 要 以復(fù)雜地形的天津薊州為例，通過對比距離薊州最近的大興探空站資料與歐洲中期天氣預(yù)報(bào)中心（European Centre for Medium-Range Weather Forecasts，ECMWF）第五代大氣再分析（ECMWF Reanalysis v5，ERA5）資料的差異，對基于ERA5資料生成的強(qiáng)對流指數(shù)在薊州的適用性進(jìn)行檢驗(yàn)和評估。結(jié)果表明：（1）ERA5資料與大興探空站探測的位勢高度、氣溫和風(fēng)速在對流層高度吻合，說明ERA5資料能夠描述薊州高空氣象條件，且對低空的表現(xiàn)能力比高空準(zhǔn)確，各要素中大氣濕度的表現(xiàn)相對較差；（2）基于ERA5生成的對流指數(shù)中，與強(qiáng)對流天氣密切相關(guān)的對流有效位能（convective available potential energy，CAPE）、K指數(shù)、沙瓦特指數(shù)（Showalter index，SI）和大氣可降水量（precipitable water，PW）與大興探空站對應(yīng)參數(shù)的相關(guān)系數(shù)分別達(dá)到0.66、0.90、0.93和0.99，表明利用ERA5構(gòu)建的對流指數(shù)能夠揭示大氣不穩(wěn)定層結(jié)條件；（3）ERA5對流指數(shù)變化與薊州降水過程相對應(yīng)，ERA5能夠反映天氣的變化和發(fā)展，為強(qiáng)對流潛勢分析提供參考。

關(guān)鍵詞 薊州；大興探空站；ERA5資料；高空探測；強(qiáng)對流指數(shù)

中圖分類號(hào)：P459文獻(xiàn)標(biāo)志碼：A文章編號(hào)：2096-3599（2024）01-0118-11

DOI：10.19513/j.cnki.hyqxxb.20221210001

收稿日期：2022-12-10；修回日期：2023-08-28

基金項(xiàng)目：國家自然科學(xué)基金項(xiàng)目（U2242202）；成都信息工程大學(xué)教師科技創(chuàng)新能力提升計(jì)劃重大項(xiàng)目（KYTD202201）；中國氣象局創(chuàng)新發(fā)展專項(xiàng)（CXFZ2022J012）

第一作者簡介：鄒雙澤，女，碩士，工程師，主要從事災(zāi)害性天氣監(jiān)測預(yù)警研究，317973133@qq.com。

通信作者簡介：白愛娟，女，博士，教授，主要從事天氣動(dòng)力學(xué)研究，baiaj@cuit.edu.cn。

Validation and application of ERA5 in complex terrain of Jizhou

ZOU Shuangze¹， BAI Aijuan²， HE Ke³， JIN Haidong⁴， HUANG Jinying⁴， JIN Bo⁵

（1. Banan Meteorological Bureau of Chongqing， Chongqing 401320， China; 2. School of Atmospheric Sciences， Chengdu University of Information Technology， Sichuan 610225， China; 3. Panzhihua Meteorological Bureau， Sichuan 617000， China; 4. Jizhou Meteorological Bureau of Tianjin， Tianjin 301999， China; 5. Baodi Meteorological Bureau of Tianjin， Tianjin 301800， China）

Abstract By comparing the differences between the data of Daxing sounding station， the nearest to Jizhou， Tianjing， and ECMWF （European Centre for Medium-Range Weather Forecasts） Reanalysis v5 （ERA5） data， the applicability of the ERA5-derived strong convection index under complex terrain conditions of Jizhou is validated and evaluated. The results are listed below. （1） The geopotential height， air temperature， and wind speed detected by Daxing sounding station are consistent with the corresponding ERA5 in the troposphere， which indicates that ERA5 can represent the meteorological characteristics above Jizhou， and it performs better at low altitude than at high altitude， but has poor performance on atmospheric humidity. （2） The convection indexes based on ERA5， such as CAPE （convective available potential energy）， K index， SI （Showalter index）， and PW （precipitable water）， are all closely related with strong convective weather， and the correlation coefficients between them and the corresponding parameters of Daxing sounding station are 0.66， 0.90， 0.93， and 0.99， respectively， indicating that these indexes can reveal the stratification instability conditions. （3） The variation of the convection indexes constructed by ERA5 corresponds with the precipitation process in Jizhou， and ERA5 can reflect the change of weather and provide reference for analyzing strong convective potential.

Keywords Jizhou; Daxing sounding station; ECMWF Reanalysis v5 （ERA5）; aerial exploration; strong convective index

引言

強(qiáng)對流天氣嚴(yán)重威脅人民生命財(cái)產(chǎn)安全，是局地氣象災(zāi)害發(fā)生的重要原因。如北京“7·21”局地特大暴雨^［1-3］，引發(fā)的城市內(nèi)澇和山區(qū)泥石流造成78人遇難。又如河南鄭州“7·20”特大暴雨^［4-5］，造成全市380人因?yàn)?zāi)死亡或失蹤。在特定天氣形勢和環(huán)境條件下產(chǎn)生的強(qiáng)對流天氣是短時(shí)臨近預(yù)報(bào)的重難點(diǎn)，高空探測為其判識(shí)提供了基本思路，探索導(dǎo)致強(qiáng)對流發(fā)生的有利環(huán)境條件是災(zāi)害預(yù)警的重要方法^［6-7］。俞小鼎等^［8］研究表明，低層水汽充沛、層結(jié)不穩(wěn)定和抬升觸發(fā)機(jī)制是強(qiáng)對流發(fā)生的有利環(huán)境條件。趙海軍等^［9］分析表明低層冷暖平流強(qiáng)迫是大氣層結(jié)不穩(wěn)定建立或維持的主導(dǎo)者?；谔娇召Y料的強(qiáng)對流指數(shù)為分析強(qiáng)對流天氣潛勢提供了重要信息，如劉玉玲^［10］在強(qiáng)對流發(fā)生物理機(jī)制的分析中，指出對流有效位能（convective available potential energy，CAPE）、粗里查遜數(shù)等參量在天氣預(yù)測中具有重要作用。Tian等^［11］肯定了強(qiáng)對流指數(shù)對我國中東部降水預(yù)報(bào)具有指示意義。王迪等^［12］利用歐洲中期天氣預(yù)報(bào)中心（European Centre for Medium-Range Weather Forecasts，ECMWF）第五代大氣再分析（ECMWF Reanalysis v5，ERA5）資料計(jì)算了多種物理參數(shù)，結(jié)果表明整層大氣可降水量、地面露點(diǎn)溫度和K指數(shù)等在短時(shí)強(qiáng)降水發(fā)生前后都有明顯的變化。

強(qiáng)對流指數(shù)用于分析對流層大氣不穩(wěn)定性，為強(qiáng)對流分析提供了途徑。高空探測資料獲取的對流層溫濕度和風(fēng)場特征，以及生成強(qiáng)對流指數(shù)，可以作為天氣分析的實(shí)況。但是探空站稀少，大多地區(qū)因缺少探空站而不能直接獲取有效的對流指數(shù)。基于多源觀測和數(shù)值模式的再分析資料，同化了各種觀測資料，具有時(shí)間序列長、覆蓋范圍廣和格點(diǎn)均勻以及分辨率高的優(yōu)點(diǎn)，能夠提供對流層各高度層的溫度、濕度、風(fēng)速和氣壓等，可以間接獲取各種對流特征參數(shù)，為沒有探空站的復(fù)雜地形區(qū)強(qiáng)對流天氣潛勢分析提供重要的資料^［13-15］。ECMWF的ERA5資料在全球得到廣泛應(yīng)用^［16］，但在復(fù)雜地形區(qū)資料的可靠性仍需要驗(yàn)證。

一方面，處于復(fù)雜地形的天津薊州短時(shí)強(qiáng)降水頻繁，夜發(fā)性顯著，近年來多次出現(xiàn)暴雨洪澇災(zāi)害^［17-18］。另一方面，薊州甚至天津本地未布設(shè)高空站，對流層大氣的基本信息量少，嚴(yán)重影響強(qiáng)對流天氣災(zāi)害預(yù)警。西來和南來的天氣系統(tǒng)是影響薊州天氣的主要因素，依據(jù)可參考性和就近原則，實(shí)際業(yè)務(wù)中使用上游北京大興探空站（薊州偏西約80 km處）來分析薊州天氣環(huán)流形勢^［19］。本文以大興探空站的資料為真值，檢驗(yàn)ERA5資料的適用性，嘗試用ERA5資料構(gòu)建薊州強(qiáng)對流天氣指數(shù)，為以后薊州的強(qiáng)降水預(yù)警模型和算法研究奠定基礎(chǔ)。

1 資料與方法

1.1 研究區(qū)域

薊州位于天津最北端，東西橫跨40 km，南北跨度75 km。地勢北高南低，北部屬于燕山余脈，最高峰九山頂海拔1 078.5 m，南部為沖積平原，屬于典型的一半平原一半山地的復(fù)雜地形區(qū)（圖1）。北京大興探空站距薊州約80 km，選擇ERA5資料在39.75°N、116.5°E的格點(diǎn)值與大興探空站進(jìn)行對比分析。

1.2 研究資料

（1）大興探空站資料。由氣象信息綜合分析處理系統(tǒng)（Meteorological Information Combine Analysis and Process System，MICAPS）獲取2014年4—10月大興探空站（54511）每日08時(shí)和20時(shí)（北京時(shí)，下同）的探測資料，共計(jì)428個(gè)時(shí)次。資料包括925、850、700、500、400、300、250、200、150和100 hPa共10個(gè)氣壓層的位勢高度、氣溫、露點(diǎn)溫度、溫度露點(diǎn)差和風(fēng)向風(fēng)速，以及CAPE值等對流指數(shù)。

（2）薊州國家級地面氣象觀測站降水資料。天津市薊州區(qū)氣象局提供的2014年暖季4—10月國家級地面氣象觀測站逐小時(shí)降水資料。

（3）歐洲中期天氣預(yù)報(bào)中心ERA5再分析資料^［20］。ERA5資料對外開放，通過https：//cds.climate.copernicus.eu/cdsapp#！/home下載，空間分辨率為0.25°×0.25°，時(shí)間分辨率為1 h。本文將用到其中2個(gè)子集：ERA5逐時(shí)氣壓層資料集，包括氣溫、垂直速度和位勢高度等多個(gè)要素，垂直方向有61個(gè)氣壓層可選，選用與探空資料一致的10個(gè)固定氣壓層（925、850、700、500、400、300、250、200、150和100 hPa）；ERA5逐時(shí)單層資料集，包含模式反演的各類物理參數(shù)，分析的強(qiáng)對流指數(shù)來自這個(gè)資料集。

（4）強(qiáng)對流指數(shù)。高空探測可獲取反映大氣溫濕層結(jié)特征的熱力學(xué)參數(shù)，用于分析強(qiáng)對流天氣發(fā)生潛勢的環(huán)境背景。常用指數(shù)有CAPE、K指數(shù)、沙瓦特指數(shù)（Showalter index，SI）和大氣可降水量（precipitable water，PW）。其中，CAPE是反映大氣層結(jié)不穩(wěn)定的重要參數(shù)，定義為從自由對流高度到平衡高度浮力能的垂直積分。K指數(shù)是綜合考慮850～500 hPa平均溫度遞減率、850 hPa露點(diǎn)溫度和700 hPa大氣飽和度的參數(shù)，能夠反映大氣低層的熱力和水汽條件。SI是氣塊從850 hPa開始按干絕熱遞減率，抬升至凝結(jié)高度后，再按濕絕熱遞減率抬升至500 hPa時(shí)，環(huán)境溫度與氣塊的溫度差。PW是單位面積氣柱內(nèi)所含的水汽總量，水汽參與是對流發(fā)生的重要條件。通過ERA5獲取的大氣基本物理量和強(qiáng)對流指數(shù)，與大興探空站進(jìn)行對比，以期對ERA5資料在薊州強(qiáng)對流天氣中的適用性進(jìn)行檢驗(yàn)分析。

1.3 研究方法

采用統(tǒng)計(jì)量平均偏差（Bias）、平均絕對誤差（mean absolute error，MAE）來檢驗(yàn)不同高度層觀測物理量和ERA5再分析資料物理量間的偏差，說明再分析資料對真實(shí)大氣的模擬效果。同時(shí)，采用相關(guān)系數(shù)（R）和均方根誤差（root mean square error，RMSE），評價(jià)再分析資料反演的強(qiáng)對流指數(shù)的適用性，相關(guān)系數(shù)越接近1，再分析資料反演的強(qiáng)對流指數(shù)越準(zhǔn)確。相關(guān)計(jì)算公式如下：

式中：V_Bias為平均偏差值，V_MAE為平均絕對誤差值，V_RMSE為均方根誤差值，R為相關(guān)系數(shù)值，x為大興探空站觀測值，y為ERA5再分析資料值，n為樣本數(shù)。

2 ERA5資料對流層各氣象要素分析

首先選取距離大興探空站最近的一個(gè)格點(diǎn)，其次挑選出與探空資料同時(shí)次的ERA5格點(diǎn)資料。對比分析多時(shí)次大興探空站與ERA5資料的位勢高度、氣溫、露點(diǎn)溫度、溫度露點(diǎn)差、水平風(fēng)速，從地面到對流層頂?shù)拇怪狈植?，確定ERA5資料揭示薊州及周邊對流層環(huán)境條件的能力。

分析兩類資料各高空氣象要素的平均廓線，發(fā)現(xiàn)兩者在位勢高度（圖2a）和氣溫（圖2b）的垂直分布上幾乎重疊，各高度層的氣溫差在0.5℃以內(nèi)，位勢高度差值非常小，說明ERA5產(chǎn)品對位勢高度和氣溫的模擬在薊州附近效果好，能夠準(zhǔn)確反映夏季對流層的氣溫和氣壓的垂直分布特征。

分析ERA5資料對大氣層水汽條件的表現(xiàn)，將ERA5相對濕度換算成與探空一致的露點(diǎn)溫度，對應(yīng)的露點(diǎn)溫度（圖2c）和溫度露點(diǎn)差（圖2d）廓線，發(fā)現(xiàn)ERA5與大興探空站露點(diǎn)值在500 hPa以下吻合，即在低層的水汽聚集區(qū)，ERA5的露點(diǎn)溫度能夠揭示大氣濕度環(huán)境條件。在500 hPa高度以上，大氣中的水汽含量相對較少，ERA5反演的露點(diǎn)溫度和探空獲取的露點(diǎn)溫度差異增加，ERA5的露點(diǎn)誤差增大，對濕度的表現(xiàn)能力下降。同時(shí)溫度露點(diǎn)差在700 hPa以下兩種資料數(shù)值接近，偏差在-1℃內(nèi)，說明ERA5對濕度的描述較準(zhǔn)確。從700 hPa向上到200 hPa，ERA5資料的溫度露點(diǎn)差比探空資料偏小，反映的大氣偏濕，結(jié)果將影響到對流層與濕度有關(guān)的對流指數(shù)反演，其中重要的就是CAPE值。

分析ERA5對大氣流場的反演能力，由圖2e—f中水平風(fēng)速的廓線看，ERA5經(jīng)向風(fēng)和緯向風(fēng)與探空觀測曲線重合率高。邊界層內(nèi)夏季平均緯向西風(fēng)風(fēng)速最小為0.7 m·s^-1，隨高度升高，緯向風(fēng)速增加，在200 hPa附近急流帶風(fēng)速超過30.0 m·s^-1，各高度上兩種資料的緯向風(fēng)吻合。經(jīng)向風(fēng)速數(shù)值小，且在850 hPa以下為正值，表明有南風(fēng)發(fā)展，平均值最大為1.5 m·s^-1。850 hPa以上轉(zhuǎn)變?yōu)樨?fù)的偏北風(fēng)，隨高度升高北風(fēng)增強(qiáng)，且在各高度上兩類資料吻合。因此ERA5風(fēng)場資料準(zhǔn)確地反映了薊州周邊大氣流場變化，且可靠性較高，這對分析當(dāng)?shù)貜?qiáng)對流天氣風(fēng)的水平切變和垂直切變非常有利。

進(jìn)一步分析各氣象要素的相關(guān)性垂直分布特征。由多種氣象要素的相關(guān)系數(shù)廓線（圖3a）對比發(fā)現(xiàn)，低層氣象要素的相關(guān)系數(shù)高，500 hPa高度以下均大于0.8。在對流層400 hPa以上，溫濕度物理量相關(guān)系數(shù)降低。對比不同氣象要素，發(fā)現(xiàn)ERA5的經(jīng)向風(fēng)速和緯向風(fēng)速與探空的風(fēng)速相關(guān)性最好，對流層從低空到高空，相關(guān)系數(shù)均達(dá)到0.9以上。氣溫的相關(guān)性較高，200 hPa以下相關(guān)系數(shù)達(dá)到0.8以上，200 hPa以上相關(guān)系數(shù)降低到0.6～0.8。相比而言，露點(diǎn)溫度和溫度露點(diǎn)差的相關(guān)性較差，300 hPa以下相關(guān)系數(shù)大于0.8，300 hPa以上相關(guān)系數(shù)驟降到0左右，甚至出現(xiàn)負(fù)值。由此說明ERA5資料在對流層500 hPa以下的低層效果良好，可以反映大氣溫濕和流場的實(shí)際狀況，但在對流層高層效果降低。

對比兩類資料各層氣溫分布（圖3b），發(fā)現(xiàn)氣溫平均偏差小于0.5℃，在邊界層內(nèi)氣溫平均偏差小于0.1℃，平均絕對誤差為0.5～1.4℃，說明ERA5的氣溫值在邊界層反演效果好。隨高度上升，氣溫下降，平均偏差和平均絕對誤差趨于增大，ERA5對低層氣溫表現(xiàn)能力優(yōu)于高層。露點(diǎn)溫度（圖3c）和溫度露點(diǎn)差（圖3d）的平均絕對誤差也表現(xiàn)為隨高度升高而增大，在低空露點(diǎn)的平均偏差和平均絕對誤差小于5.0℃。但到了對流層高層，露點(diǎn)的誤差快速增長，平均偏差超過5.0℃。溫度露點(diǎn)差的平均偏差在700 hPa以下小于-2.0℃，平均絕對誤差在5.0℃。在300 hPa附近溫度露點(diǎn)差的平均偏差達(dá)到最大值4.5℃。

分析兩類資料的水平風(fēng)速差異的廓線（圖3e—f），發(fā)現(xiàn)從邊界層到對流層中層500 hPa，ERA5反演的水平風(fēng)速平均絕對誤差小于1.5 m·s^-1，能夠滿足薊州地區(qū)天氣分析和研究的需求。

綜上所述，大興探空站資料與ERA5資料的對比分析說明ERA5氣象要素產(chǎn)品在薊州臨近地區(qū)能較好描述大氣基本狀態(tài)，其中位勢高度、氣溫和風(fēng)速與探空資料吻合，誤差較小，尤其是各要素在對流層低層的表現(xiàn)能力比高層優(yōu)越，但濕度的誤差相對較大，尤其是在大氣高層。

3 ERA5強(qiáng)對流指數(shù)的檢驗(yàn)分析

3.1 對流有效位能

ERA5直接對外發(fā)布CAPE產(chǎn)品，探空資料的CAPE來自MICAPS系統(tǒng)，對比2014年4—10月大興探空站資料與ERA5資料的CAPE值散點(diǎn)分布（圖4a），兩類資料的相關(guān)系數(shù)為0.66，通過了0.01顯著性檢驗(yàn)，說明ERA5資料的CAPE值在薊州具有可用性。以探空資料CAPE為x，ERA5資料的CAPE值為y，兩類資料序列的回歸方程為y=0.68x+53.20，表明ERA5反演的CAPE值整體略偏低。

大氣的CAPE值較多出現(xiàn)在500 J·kg^-1以下，分析CAPE小于和大于500 J·kg^-1兩種情況以說明ERA5資料的可靠性，其中高CAPE值對強(qiáng)對流天氣分析至關(guān)重要。在CAPE低值區(qū)（圖4b），ERA5的CAPE值較探空偏大；在CAPE高值區(qū)（圖4c），ERA5的CAPE值較探空資料值略偏小。即在不同閾值區(qū)間，CAPE值的誤差不同，且在CAPE的低值區(qū)誤差相對較小，隨CAPE值增大，誤差變大。

由表1中ERA5資料的CAPE與大興探空站對應(yīng)值各月的相關(guān)系數(shù)分布上看，4月兩類資料的CAPE相關(guān)系數(shù)達(dá)0.874，表明相關(guān)性好。其次是8月和9月相關(guān)系數(shù)較高，分別為0.662和0.653，通過0.01的顯著性檢驗(yàn)。5月和6月的相關(guān)系數(shù)值最低。以上分析說明在夏季強(qiáng)對流天氣活躍的6—8月，ERA5格點(diǎn)資料CAPE產(chǎn)品在薊州區(qū)的表現(xiàn)能力優(yōu)越，能夠準(zhǔn)確地反映薊州復(fù)雜地形下的大氣不穩(wěn)定狀態(tài)。CAPE值是對流層低空溫濕度的綜合反映，正是ERA5溫濕度層結(jié)數(shù)據(jù)在低層的良好表現(xiàn)，才是CAPE正確的原因。

3.2 K指數(shù)

K指數(shù)常用在強(qiáng)對流天氣潛勢分析中，在實(shí)際業(yè)務(wù)工作中是一個(gè)非常有用的預(yù)報(bào)參數(shù)，強(qiáng)對流發(fā)生時(shí)，通常K指數(shù)值越大，對流強(qiáng)度越大。對比分析大興探空站與ERA5資料獲取的K指數(shù)的散點(diǎn)分布（圖5），發(fā)現(xiàn)兩種資料的K指數(shù)都集中在20～40 K。ERA5資料的K指數(shù)和探空站的K指數(shù)顯著相關(guān)，相關(guān)系數(shù)為0.90。同上一節(jié)分析類似，兩類資料的回歸方程可表示為y=0.79x+5.56，平均偏差為2.02 K，均方根誤差為7.24 K。且當(dāng)K指數(shù)小于0時(shí)，ERA5較探空站的K指數(shù)偏大，隨著K指數(shù)增大，偏差變小，資料的離散度減小。這說明ERA5資料的K指數(shù)在薊州臨近地區(qū)與探空資料非常吻合，分析強(qiáng)對流天氣可行性高。

K指數(shù)逐月相關(guān)系數(shù)高，說明ERA5反演的K指數(shù)在該地區(qū)應(yīng)用效果良好（表1）。其中10月的相關(guān)系數(shù)高達(dá)0.891。6月、7月和8月的相關(guān)系數(shù)分別為0.722、0.803和0.736，數(shù)值略有降低。6月、7月和8月中K指數(shù)曲線波動(dòng)振幅?。▓D6），也是K指數(shù)平均值最大的3個(gè)月，反映了盛夏季節(jié)大氣水汽含量充沛，層結(jié)穩(wěn)定性下降，對流發(fā)展較頻繁的基本特征。

3.3 沙瓦特指數(shù)

實(shí)際業(yè)務(wù)工作中SI可用于強(qiáng)對流天氣的判識(shí)和分析，具有實(shí)用性。已有研究^［12］表明，SI為負(fù)數(shù)且絕對值大時(shí)，對流強(qiáng)度越大。ERA5資料未發(fā)布SI產(chǎn)品，本文通過天氣學(xué)算法公式^［21］得到SI。SI相關(guān)系數(shù)達(dá)到0.93，比K指數(shù)的相關(guān)性高，數(shù)據(jù)散點(diǎn)分布高度集中在0～10 K（圖7）。兩類資料的回歸方程為y=0.89x+0.42，可以看到兩者接近于y=x，說明ERA5的SI與實(shí)況吻合，應(yīng)用效果佳。且當(dāng)SI小于0 K時(shí)，ERA5較探空資料偏大，當(dāng)SI大于0 K時(shí)，ERA5較探空資料偏小。隨著SI增加，兩類資料的離散度減小，尤其當(dāng)SI大于10 K時(shí)，二者的偏差最小。4月、9月和10月ERA5資料和探空資料SI的相關(guān)系數(shù)分別達(dá)到0.921、0.955和0.924，ERA5反演的SI效果優(yōu)異。

3.4 大氣可降水量

大氣可降水量（PW）表示整層大氣的水汽含量，業(yè)務(wù)中常用來分析水汽條件和識(shí)別對流類型。PW數(shù)值越大，表示可能的降水量越大。ERA5資料直接提供了PW產(chǎn)品，對比結(jié)果表明ERA5和大興探空站PW顯著相關(guān)（圖8），相關(guān)系數(shù)高達(dá)0.99。線性回歸方程為：y=0.97x+1.26。表1中4—10月兩類資料PW的相關(guān)系數(shù)均大于0.95，并通過0.01的顯著性檢驗(yàn)。分析圖6可知，PW的平均值在6月、7月和8月較高，分別為30.67 mm、39.28 mm和34.23 mm。

4 降水分析與強(qiáng)對流指數(shù)

2014年夏季薊州出現(xiàn)了多次強(qiáng)降水事件，獲取距薊州站最近的ERA5格點(diǎn)資料，時(shí)間選擇08時(shí)和20時(shí)，分析強(qiáng)對流指數(shù)與薊州站隨后12 h累計(jì)降水量的對應(yīng)關(guān)系，探究ERA5在薊州強(qiáng)降水分析中的可用性。分析表2中ERA5的CAPE值與薊州站累計(jì)降水量的相關(guān)系數(shù)，發(fā)現(xiàn)各月的相關(guān)系數(shù)不確定，在4月、5月出現(xiàn)了負(fù)值，其余月均為正值，10月相關(guān)系數(shù)最大，達(dá)到0.351。各月K指數(shù)與降水的相關(guān)系數(shù)均為正，符合對流活動(dòng)與K指數(shù)的對應(yīng)關(guān)系，即K指數(shù)越高，降水量越大，其中5月和7月相關(guān)系數(shù)值較高，分別為0.297和0.283；除7月外，PW與降水的相關(guān)系數(shù)在4個(gè)強(qiáng)對流指數(shù)中最大，符合該物理量是大氣水汽含量直接表達(dá)的事實(shí)，充分說明PW在強(qiáng)降水預(yù)報(bào)中具有很好的實(shí)用性。

進(jìn)一步研究強(qiáng)對流指數(shù)與降水關(guān)聯(lián)性發(fā)現(xiàn)，薊州站有明顯降水過程發(fā)生之前，存在CAPE和K指數(shù)的峰值（圖6），但是降水峰值與強(qiáng)對流指數(shù)峰值并不能完全對應(yīng)。檢查薊州2014年多次降水過程，ERA5的強(qiáng)對流指數(shù)在過程中都有明顯變化。以2014年6月6—7日的強(qiáng)降水過程為例（圖9），對比降水量與強(qiáng)對流指數(shù)逐小時(shí)變化曲線，發(fā)現(xiàn)薊州站過程累計(jì)雨量達(dá)65.6 mm，最大小時(shí)雨量為45.1 mm，對應(yīng)平均CAPE為188.5 J·kg^-1，CAPE最高值為760.6 J·kg^-1。6日11時(shí)開始大氣不穩(wěn)定能量加速累積。14時(shí)降水發(fā)生，CAPE值高達(dá)673.8 J·kg^-1。16時(shí)之后CAPE數(shù)值驟降，大氣不穩(wěn)定能量釋放，降水以小雨量級為主。在降水趨于結(jié)束時(shí)，K指數(shù)、SI和PW開始有較明顯變化。

基于以上個(gè)例可知ERA5的強(qiáng)對流指數(shù)與對流活動(dòng)的發(fā)生有密切關(guān)系，能夠指示和預(yù)警強(qiáng)降水的發(fā)生。隨著降水臨近，PW、CAPE、K指數(shù)數(shù)值不斷增大，說明大氣中水汽充沛，不穩(wěn)定能量不斷積累，有利于強(qiáng)對流天氣的發(fā)生。但是該地強(qiáng)降水發(fā)生并不能與強(qiáng)對流指數(shù)變化完全對應(yīng)，究其原因，該地區(qū)降水多由積層混合云形成，積云降水和層云降水是不同的兩類降水微物理過程的結(jié)果^［22］，對于穩(wěn)定性的層狀云降水，通常強(qiáng)對流指數(shù)在降水前后數(shù)值變化較小。強(qiáng)對流指數(shù)能夠反映大氣層結(jié)不穩(wěn)定和對流環(huán)境條件，但不是對流和強(qiáng)降水發(fā)展的唯一原因，降水的形成還與降水性質(zhì)、風(fēng)切變強(qiáng)度和局地地形等多因素有關(guān)。

5 結(jié)論

通過對比分析ERA5資料與大興探空站資料在對流層高空基本氣象要素以及CAPE、K指數(shù)、SI和PW等相關(guān)的強(qiáng)對流指數(shù)，說明ERA5資料在復(fù)雜地形背景下薊州的適用性，獲得以下結(jié)論。

（1）以探空資料為真值的高空基本氣象要素的檢驗(yàn)結(jié)果表明：位勢高度、氣溫和風(fēng)速產(chǎn)品的平均廓線幾乎重合，ERA5資料能夠反映基本氣象要素的垂直分布特征；由進(jìn)一步的偏差分析看出，各層氣溫的平均偏差小于0.5℃，中低層風(fēng)速平均偏差在0.5 m·s^-1以內(nèi)。而代表濕度條件的物理量溫度露點(diǎn)差的差異顯著，在對流層高層最大達(dá)4.5℃。ERA5資料對大氣濕度的描述能力相對較差，將會(huì)影響對流指數(shù)產(chǎn)品精度。

（2）基于ERA5資料獲取和生成的CAPE、K指數(shù)、SI和PW與探空資料對應(yīng)吻合，相關(guān)系數(shù)分別達(dá)0.66、0.90、0.93和0.99，可以準(zhǔn)確地反映當(dāng)?shù)貙α鲗拥牟环€(wěn)定特征。

（3）研究薊州降水與強(qiáng)對流指數(shù)的聯(lián)系發(fā)現(xiàn)，復(fù)雜地形下薊州降水的發(fā)生與強(qiáng)對流指數(shù)的變化有較好的對應(yīng)關(guān)系。PW的峰值與降水峰值在時(shí)間軸上接近，其次是SI、K指數(shù)和CAPE。各對流指數(shù)與薊州的降水過程變化并不能完全對應(yīng)，原因在于強(qiáng)對流特別是強(qiáng)降水物理過程復(fù)雜，受大氣不穩(wěn)定和地形等多種因素影響，不能由單一的因子決定。

參考文獻(xiàn)：

［1］ 諶蕓，孫軍，徐珺，等.北京721特大暴雨極端性分析及思考（一）觀測分析及思考［J］.氣象，2012，38（10）：1255-1266.

［2］孫建華，趙思雄，傅慎明，等.2012年7月21日北京特大暴雨的多尺度特征［J］.大氣科學(xué)，2013，37（3）：705-718.

［3］孫軍，諶蕓，楊舒楠，等.北京721特大暴雨極端性分析及思考（二）極端性降水成因初探及思考［J］.氣象，2012，38（10）：1267-1277.

［4］齊道日娜，何立富，王秀明，等.“7·20” 河南極端暴雨精細(xì)觀測及熱動(dòng)力成因［J］.應(yīng)用氣象學(xué)報(bào)，2022，33（1）：1-15.

［5］何秉順.河南鄭州山區(qū)4市2021年“7·20” 特大暴雨災(zāi)害調(diào)查的思考與建議［J］.中國防汛抗旱，2022，32（3）：37-40，51.

［6］王秀明，俞小鼎，周小剛.雷暴潛勢預(yù)報(bào)中幾個(gè)基本問題的討論［J］.氣象，2014，40（4）：389-399.

［7］張澤銘，劉強(qiáng)，張熙，等.一次低渦切變線引發(fā)的暴雨天氣過程分析［J］.山東氣象，2011，31（3）：10-14.

［8］俞小鼎，鄭永光.中國當(dāng)代強(qiáng)對流天氣研究與業(yè)務(wù)進(jìn)展［J］.氣象學(xué)報(bào)，2020，78（3）：391-418.

［9］趙海軍，潘玲，毛子卿.山東省持續(xù)性短時(shí)強(qiáng)降水過程物理量特征分析［J］.海洋氣象學(xué)報(bào)，2023，43（1）：63-74.

［10］劉玉玲.對流參數(shù)在強(qiáng)對流天氣潛勢預(yù)測中的作用［J］.氣象科技，2003，31（3）：147-151.

［11］TIAN F Y， ZHENG Y G， ZHANG T， et al. Statistical characteristics of environmental parameters for warm season short-duration heavy rainfall over central and Eastern China［J］. J Meteorol Res，2015，29（3）：370-384.

［12］王迪，牛淑貞，曾明劍，等.河南省分類強(qiáng)對流環(huán)境物理?xiàng)l件特征分析［J］.氣象，2020，46（5）：618-628.

［13］王秀明，俞小鼎，朱禾.NCEP再分析資料在強(qiáng)對流環(huán)境分析中的應(yīng)用［J］.應(yīng)用氣象學(xué)報(bào)，2012，23（2）：139-146.

［14］孫繼松，陶祖鈺.強(qiáng)對流天氣分析與預(yù)報(bào)中的若干基本問題［J］.氣象，2012，38（2）：164-173.

［15］鄭媛媛，姚晨，郝瑩，等.不同類型大尺度環(huán)流背景下強(qiáng)對流天氣的短時(shí)臨近預(yù)報(bào)預(yù)警研究［J］.氣象，2011，37（7）：795-801.

［16］JIANG Q， LI W Y， FAN Z D， et al. Evaluation of the ERA5 reanalysis precipitation dataset over Chinese Mainland［J］. J Hydrol，2021，595：125660.

［17］熊明明，徐姝，李明財(cái)，等.天津地區(qū)小時(shí)降水特征分析［J］.暴雨災(zāi)害，2016，35（1）：84-90.

［18］鄒雙澤，白愛娟，袁淑杰，等.薊州區(qū)半山地形下汛期短歷時(shí)強(qiáng)降水特征研究［J］.成都信息工程大學(xué)學(xué)報(bào)，2021，36（6）：666-673.

［19］尉英華，花家嘉，王瑩，等.近11年天津冰雹統(tǒng)計(jì)特征及對流參數(shù)指標(biāo)分析［J］.氣象，2023，49（2）：213-223.

［20］HERSBACH H， BELL B， BERRISFORD P， et al. The ERA5 global reanalysis［J］. Quart J Roy Meteor Soc， 2020，146（730）：1999-2049.

［21］陸漢城，楊國祥.中尺度天氣原理和預(yù)報(bào)［M］.2版.北京：氣象出版社，2004.

［22］李芳，李南，萬瑜.臺(tái)風(fēng)“溫比亞” 影響山東期間GPM資料的降水分析［J］.海洋氣象學(xué)報(bào)，2020，40（4）：69-76.