应用气象学报  2012, 23 (2): 139-146   PDF    
NCEP再分析资料在强对流环境分析中的应用
王秀明, 俞小鼎, 朱禾     
中国气象局气象干部培训学院,北京 100081
摘要: 为考察NCEP再分析资料在我国强对流天气产生环境分析中的适用性,选取2002—2009年多普勒天气雷达识别的60例超级单体风暴个例,对比分析常规探空资料和NCEP再分析资料提取的温、湿、风垂直廓线,结果表明:NCEP再分析资料计算的对流有效位能因对抬升气团湿度敏感而与观测间差异较大,宜用K指数、温度直减率分析大气层结稳定度;因对流层中高层风与探空差异不大,其中500~700 hPa的风与探空近乎一致,因此NCEP再分析资料计算的深层、中层风垂直切变参量可靠性较高;NCEP再分析资料水汽参数与探空资料差异大,特别是在大气边界层,需用观测资料订正;边界层物理量,特别是风向与探空差异显著,因此不宜用NCEP再分析资料讨论雷暴触发问题;平均而言,NCEP再分析资料湿度廓线低层偏干而中层偏湿,925 hPa以上风速偏小,降低了强对流发生概率。
关键词: NCEP再分析资料    探空资料    超级单体    强对流    
The Applicability of NCEP Reanalysis Data to Severe Convection Environment Analysis
Wang Xiuming, Yu Xiaoding, Zhu He     
China Meteorological Administration Training Center, Beijing 100081
Abstract: Operational sounding observation is carried out twice per day, but four similar vertical profiles can be obtained from NCEP reanalysis data. Therefore, the vertical profiles obtained from NCEP reanalysis data are assessed and its applicability in diagnosing severe convection environment is analyzed. Sixty soundings in close proximity to supercell storms are investigated. These supercell storms are observed in China and identified by Doppler weather radar. The soundings which may be polluted by storm are replaced by soundings from its downstream area or upstream area. The profiles obtained from NCEP data are at the same time and as close to the observation soundings as possible. The vertical profile data that obtained from NCEP reanalysis data are compared with soundings. The ingredients-based methodology is used to discuss the potential for severe convection. Vertical wind shear is the first element to check. NCEP wind in middle and high troposphere is almost consistent with observation, so the deep and middle vertical wind shear (0—6 km, 0—8 km and 0—3 km) and other related dynamical parameters can be calculated from NCEP data except for 0—1 km low vertical wind shear, because wind difference is significant in planet boundary layer (PBL). To study instability of storms, convective available potential energy (CAPE), temperature lapse rate (temperature difference between 850 hPa and 500 hPa), and K index are investigated. The statistic results show that the difference of CAPE between observation and NCEP profiles is significant, because CAPE is sensitive to dewpoint and temperature of lifting air mass. On average, 1℃ temperature increment of lifting airmass brings 200 J·kg-1CAPE augmentation, 1℃ dewpoint increases CAPE by nearly 500 J·kg-1, and the augmentation can vary from 0 to 1000 J·kg-1. Moisture is one of the three ingredients for thunderstorm. The error of NCEP moisture parameter is significant, especially within planet boundary layer, the average dewpoint difference between NCEP and observation is 2℃ in low troposphere, which results in nearly 1000 J·kg-1CAPE difference. To calculate CAPE from NCEP data, moisture should be corrected according to observation. Temperature lapse rate can be used to diagnose atmospheric instability instead of CAPE, as temperature profile can be used to analysis severe convection. The difference of K index is small in most cases. NCEP output variable CAPEsfc (surface CAPE) is unreasonably small, but the tendency can indicate the change of CAPE. Most lifting processes are within PBL, where the difference of the atmospheric parameters especially the wind direction between observation and NCEP data is significant, so it's not suitable to use NCEP data to study lifting mechanism of thunderstorm. On average, the NECP moisture profile is much drier at low level and wetter at middle level than observation, and wind speed above 925 hPa is weaker than observation, which lowers the possibility of severe convection.
Key words: NCEP reanalysis data     sounding     supercell     severe convection    
引言

超级单体是最强烈的对流风暴,其发生频率虽然不高,但相当比例的极端强对流天气事件 (包括5 cm以上冰雹、33 m/s以上雷暴大风和F2级以上龙卷) 是由超级单体风暴产生的。因此,超级单体风暴发生的环境背景代表了有利于极端强对流事件发生的环境背景。对我国超级单体发生环境背景的分析,有利于在了解超级单体发生环境的基础上做出更准确的潜势预报。现有的探空资料每天两次 (00:00和12:00, 世界时,下同),且空间分辨率一般为200~300 km,而强对流事件常发生在午后,空间尺度常在100 km以下,因此常规探空资料常因资料时空分辨率太低而缺乏代表性。廖晓农等[1]研究表明,增加午后探空能提高北京雷暴潜势预报能力。NCEP再分析资料在模式的基础上加入了大量的观测信息,每天4次,目前其分辨率已达25 km。我国有不少学者利用NCEP再分析资料进行对流暴雨和冰雹、雷暴大风等强对流天气的环境潜势研究[2-7],强对流和对流性暴雨数值模拟研究也常基于NCEP再分析资料[8-9],有不少研究者直接用NCEP再分析资料提取的垂直廓线统计雷暴大风环境探空特征[10-12],并得到了一些有意义的结论,但NCEP再分析资料用于研究强对流、特别是极端强对流事件环境潜势的可行性研究目前尚未涉及。本文通过对比分析探空与NCEP再分析资料描述的60个超级单体强对流事件环境探空的动力、热力参数,利用NCEP再分析资料探讨强对流潜势的适用性。决定对流风暴类型的主要参数是浮力 (一般用对流有效位能表征) 和风垂直切变,观测分析和数值模拟研究表明:风垂直切变对超级单体发展尤为重要,超级单体风暴发生的环境为相对大的热浮力和强的垂直风切变[13-14]。风垂直切变被称为风暴的增强条件,强风垂直切变使得风暴结构更有组织,有利于强风暴的维持。垂直切变是超级单体内中层强烈的旋转 (在雷达图上常表现为中气旋) 的来源,一般用0~6 km风垂直切变表示,国外的观测分析和数值模拟研究表明:6 km以下风垂直切变达到15~20 m·s-1是超级单体发生所必须的,RUC-2分析和预报系统 (Rapid Update Cycle-2 Analysis and Forecast Sys-tem) 同化资料统计结果为超过20 m/s[13]。对流有效位能和0~6 km的风垂直切变是两个常用来诊断超级单体发生潜势的物理量,因此下面主要对NCEP再分析资料计算的这两个物理量进行统计。本文重点在于探讨NCEP再分析资料提取的垂直廓线在强对流环境分析中的适用性、NCEP再分析资料计算的物理量与观测资料差异较大的原因及可能的订正方法。本文考察的对象为超级单体强对流,其结果可推广到利用NCEP再分析资料分析一般雷暴天气潜势。

1 资料与方法 1.1 资料

本文选取2002—2009年我国多普勒天气雷达识别的60个超级单体强对流事件,其产生的天气现象包括冰雹、龙卷、雷雨大风和短时强降水,83%的个例发生在午后到傍晚,早晨的探空资料常常缺乏很好的代表性。88%的个例发生在我国华北和中东部地区,这些地区也是我国的强对流多发地。从时间和空间分布来看,本文选取的超级单体样本能够较好地代表发生在我国的超级单体。探空资料的选取原则为尽量靠近超级单体出现地,且未受对流活动影响。NCEP再分析资料水平分辨率为1°×1°,垂直为26层,每天4次,时间间隔为6 h。与探空资料进行对比的NCEP再分析资料从与观测同一时刻且离测站最近的一个网格点上提取。对比的动力热力参量主要为对流有效位能和水平风垂直切变。为了进一步说明NCEP再分析资料计算的对流有效位能、风垂直切变等参量与观测差异偏大的可能原因,并探讨相关物理量的可用性,本文对60个个例探空资料中基本物理量温、湿、风进行统计,进行对比的垂直层次选取探空的常规观测层:1000,925,850,700,500,400,300,250 hPa和200 hPa。

1.2 方法

对流有效位能的获取方式有多种,本文主要采用GrADS软件附带的tlogps.gs脚本计算,这也是目前国内使用较多的一种计算方式;预报业务上一般采用MICAPS探空分析软件计算;另外,NCEP再分析资料直接输出的变量中也有地面对流有效位能 (CAPEsfc)。下文将通过统计比较说明上述3种常用的对流有效位能获取方式间的差异。

0~6 km水平风垂直切变用上下两层风矢量差大小除以两层间的距离计算,单位为1.0×10-3s-1。上层取500 hPa,下层探空资料取底层,NCEP再分析资料取1000 hPa。本文选取的个例大多在东部平原地区,65%的个例探空第1层资料在985~1015 hPa,80%的个例在925 hPa以上,所以尽管下层高度,探空资料与NCEP再分析资料略有差异,但差异不大。因此,上述求取的虽不是严格意义上的0~6 km风垂直切变,但两类资料求取的风垂直切变具有可比性。

2 风垂直切变参数特征

图 1表明,由NCEP再分析资料计算的0~6 km风垂直切变与探空差异小,两者几乎相等,75%的个例差异在0.7×10-3s-1以下。NCEP再分析资料计算的风垂直切变平均值为3.28×10-3s-1,略小于探空值 (3.41×10-3s-1),两者仅相差0.14×10-3s-1。用NCEP再分析资料计算的0~6 km风垂直切变参数是可靠的,可用于强对流环境风垂直切变分析。下面通过逐层风对比分析NCEP再分析资料风垂直切变差异小的原因,进而讨论其他动力参量适用性。

图 1. NCEP再分析资料与探空资料的0~6 km风垂直切变散点图 (a) 及风垂直切变差分布图 (b) Fig 1. The scatter diagram of 0—6 km vertical wind shear (a) and the shear difference (b) between NCEP reanalysis data and observations

逐层风速差统计表明 (图 2a):500~925 hPa的NCEP再分析资料风速与探空资料差值多为2~3 m·s-1,75%的个例差值在3 m·s-1以下,且各层风速偏差值相对集中。美国探空观测风速的精度为1.5 m·s-1[15],从这个意义上看,NCEP再分析资料风速的代表性较好。平均而言,NCEP再分析资料风速比探空资料小 (图 3a)。NCEP再分析资料风向有两个显著特点 (图 2b):700 hPa以上NCEP再分析资料风向与探空资料差异小,边界层的风向差异大。700 hPa以上层次超过75%的个例风向差在20°以内,其中500 hPa风向差最小,75%的个例风向差小于10°,几乎所有的个例风向差均在20°以下。美国探空资料风向的观测精度一般为5°,从这个意义上说,可以用NCEP再分析资料的500 hPa风向代替探空风向。可以推断,500 hPa的NCEP再分析资料环流型与探空资料的环流型基本一致。另一方面,边界层的风向差异显著,特别是1000 hPa,50%个例的风向差近45°,NCEP再分析资料边界层风向几乎不可用,因此,由边界层风向造成的低层辐合强迫抬升可靠性差,不宜用NCEP再分析资料边界层辐散、辐合讨论对流触发问题。

图 2. NCEP再分析资料与探空资料的风速差值 (a) 和风向 (b) 分布 Fig 2. The differences of wind direction (a) and speed (b) between NCEP reanalysis data and observations

图 3. NCEP再分析资料与探空资料各层风速 (a) 和温度、露点 (b) 的平均差 Fig 3. Mean differences of temerature, dew point and wind speed between NCEP reanalysis data and observations

一般低层风速比高层风速小,因此0~6 km的风垂直切变很大程度上依赖于高层风,分析表明:500 hPa附近NCEP再分析资料风场与探空资料差异小,由此计算的风垂直切变与探空资料差异自然也小。NCEP再分析资料700 hPa风场与探空资料差异小,NCEP再分析资料给出的0~3 km风垂直切变也应是可用的。NCEP再分析资料中边界层风场与探空资料差异显著,可知0~1 km风垂直切变量与探空资料差异较大。上文逐层风向、风速差统计表明,中高层风与探空几乎一致,因而可以推断,与其有关的诊断量,如风暴相对螺旋度等均可用NCEP再分析资料计算。另外,因500 hPa风向与探空近乎一致,风速差异也较小,因此可以用NCEP再分析资料的500 hPa环流场来考察产生雷暴的天气系统演变,这已经被广泛应用到短时强降水、雷暴大风等强对流天气分析中[2-7]

3 对流有效位能等热力参数特征

图 4表明:由探空资料和NCEP再分析资料计算得到的对流有效位能 (CAPE) 有一定的代表性,大多数点分布在y=x附近,50%的对流有效位能差值在700 J·kg-1以下,但两者计算的对流有效位能差异很明显。利用探空资料计算的对流有效位能超过2000 J·kg-1的6个个例,利用NCEP再分析资料的计算结果显著偏大,均超过4000 J·kg-1,即NCEP再分析资料计算的对流有效位能在高值区有偏大的倾向;另一方面,NCEP再分析资料计算的对流有效位能在低值区偏小,2000 J·kg-1以下的部分个例NCEP再分析资料计算的对流有效位能显著偏小,使得NCEP再分析资料计算的超级单体的对流有效位能比探空资料计算结果离散度大。虽然两者平均值相差仅为25 J·kg-1,但平均差距达946 J·kg-1。从图 4b可以看出:25%的个例NCEP再分析资料与探空资料计算结果差值超过1500 J·kg-1,除去5%的奇异点,最大差值仍达2500 J·kg-1。NCEP再分析资料计算的对流有效位能与探空资料计算结果差异显著的原因在于其对抬升气团的温、湿状况敏感,从表 1可知,抬升气团 (最低层) 温度增加1℃,对流有效位能平均增加194 J·kg-1,露点温度增加1℃,对流有效位能平均增加459 J·kg-1。下面将通过NCEP再分析资料温、湿参量与探空资料计算结果差值统计,进一步分析对流有效位能差异大的原因,并分析用NCEP再分析资料计算的其他热力参量的可靠性。

图 4. NCEP再分析资料与探空资料计算的对流有效位能散点图 (a) 及对流有效位能差值分布 (b) Fig 4. Same as in Fig. 1, exept for CAPE

表 1 不同资料及不同方式得到的对流有效位能及其与探空资料计算结果差值 (单位: J·kg-1) Table 1 Average CAPE and CAPE deviation calculated by different dataset and different methods (unit: J·kg-1)

不同方式计算的对流有效位能差异较大,特别是NCEP再分析资料直接输出的地面对流有效位能 (CAPEsfc)。由表 1可知:NCEP再分析资料直接输出的地面对流有效位能与探空资料计算结果相比明显偏小,不到其60%。国内研究者利用NCEP再分析资料的地面对流有效位能统计雷暴天气不稳定能量时也同样偏小[10-12]。因此,NCEP再分析资料直接输出的地面对流有效位能不能代表雷暴,特别是强对流环境不稳定能量大小。研究还表明,其变化对雷暴有一定的指示意义[10-12],其他模式计算的对流有效位能也有同样的特点[16-17]。日常业务中多使用MICAPS探空分析软件包计算对流有效位能,统计表明:该结果较本文计算结果平均偏小约100 J·kg-1,可能是自由对流高度和平衡点计算的截断误差造成的,但这种差异与用不同资料计算的对流有效位能相比几乎可以忽略,比抬升气团湿度影响造成的差异小得多,对估计强对流潜势影响不大。

NCEP再分析资料温度差值统计 (图 5a) 表明:NCEP再分析资料的温度与探空资料接近,特别是400~850 hPa的自由大气,NCEP再分析资料温度差多在1℃以下; 300 hPa以下层次,约60%的个例温度差小于1℃。边界层内温度差异略大,但也多在2℃内。平均而言,NCEP再分析资料温度比探空资料略偏低,500 hPa以下层次偏低多在0.5℃以内 (图 3b);中层的温度差小而高、低层的温度差略有增大,这与陈子通等[18]分析广州中尺度数值模式结论一致,值得一提的是,NCEP再分析资料近地层的温度差与其相比偏小。NCEP再分析资料以数值模式结果为基础,其垂直廓线特点与一般模式有很多相似之处,一般来说,模式对温度的模拟结果比其他物理量好[16-18]。本文统计的NCEP再分析资料温度差比40 km分辨率的RUC-2分析和预报系统给出美国超级单体探空误差偏大2~3倍,这与采用的NCEP再分析资料空间分辨率约100 km而RUC-2分析和预报系统仅为40 km有关,受NCEP再分析资料分辨率限制,本文所选取的NCEP再分析资料与探空资料空间距离只能保证小于50 km,另外RUC-2分析和预报系统同化的观测资料更多。总体而言,NCEP再分析资料温度廓线能较好地描述超级单体环境温度层结状况。

图 5. NCEP再分析资料与探空资料的温度差 (a) 和露点温度差 (b) 分布 Fig 5. Same as in Fig. 2, except for temperature (a) and dew temperature (b)

图 5b可知,除850 hPa露点温度差较小外 (略小于2℃),露点温度差均在2℃以上,且离散度大,资料代表性较差。大气中90%的水汽集中在500 hPa以下,因此本文仅考察500 hPa以下的露点温度。1000 hPa高度最大差值超过4℃,即使是代表性较好的850 hPa,最大露点温度差值也达到6℃。平均而言,NCEP再分析资料显示低层偏干而中层略偏湿 (图 3b),低层偏干使其对流有效位能偏小,中层偏湿不利于产生雷暴大风、冰雹等强对流天气,总体上降低了灾害性强对流天气产生的概率。模式对水汽的模拟能力普遍较差[12-14],NCEP再分析资料虽然同化了大量探空资料,但仍带有模式本身存在的缺陷。

以上温、湿参量分析表明:在低层,NCEP再分析资料与探空资料的温度差为1~2℃,大部分露点差在2℃以上,因而使得两种资料计算的对流有效位能的差异平均超过900 J·kg-1。温度增加1℃导致的对流有效位能增量集中在200 J·kg-1附近 (表 1),因而可以通过观测或预报的地面温度变化定量估计可能造成的对流有效位能变化,从而能更好地估计强对流发生地、强对流发生时的对流有效位能;露点温度导致的对流有效位能变化显著,露点温度增加1℃使75%的个例对流有效位能增加超过500 J·kg-1,且其变化范围大,最大值达1000 J·kg-1。低层露点温度差是导致NCEP再分析资料计算的对流有效位能差异显著的主要原因。为了得到强对流发生环境比较有代表性的对流有效位能,宜用微波辐射计、地面加密观测或雷达反演的水汽资料对NCEP再分析资料低层水汽进行修正。水汽是雷暴发生三要素之一,低层水汽凝结潜热释放是强雷暴发展的主要能量来源,水汽的垂直分布影响对流不稳定状况,因此,NCEP再分析资料用于强对流天气环境分析时,其湿度特别是低层湿度宜用观测订正。

因为同样对低层温度敏感,NCEP再分析资料计算的对流抑制能量 (CIN) 与探空差异显著,几乎不可用 (图略)。图 6表明:NCEP再分析资料与探空资料的850 hPa与500 hPa温度差 (以下简称T85) 几乎相等,925 hPa与700 hPa温度差也有同样的特点 (图略),可以用NCEP再分析资料温度计算温度直减率以表征大气层结稳定度。两种资料计算的超级单体温度直减率T85分布于23~35℃,大多数在30℃以下 (7℃·km-1以下),为相对弱的条件不稳定。图 6表明:NCEP再分析资料与探空资料的K指数相关性明显好于对流有效位能,因其所采用的资料为边界层以上的对流层中层大气;NCEP再分析资料与探空资料的K值大多集中在32~42℃之间,因两者结果相关性好,可以作为超级单体层结不稳定参数阈值的参考。需要说明的是,K指数仍包含水汽参量,其中700 hPa露点差值较大,因而也会出现少数K指数差异显著的个例。

图 6. NCEP再分析资料与探空资料K指数计算的热力参数散点图 (a)850 hPa与500 hPa温度差,(b)K指数 Fig 6. Scattlerplot of thermodynamical parameters between NCEP reanalysis data and observations (a) temperature difference between 850 hPa and 500 hPa, (b)Kindex

4 小结

1) 对流有效位能对抬升层气团温湿敏感,平均而言,温度增加1℃,对流有效位能增加200 J·kg-1,露点温度增加1℃,对流有效位能增加500 J·kg-1。NCEP再分析资料计算的对流有效位能偏差主要源于NCEP再分析资料温、湿廓线中低层湿度 (露点) 与探空差异较大。NCEP再分析资料中近地面露点平均差大于2℃,温度平均差值约为1.5℃,给NCEP再分析资料计算的对流有效位能带来大的不确定性。

2) NCEP再分析资料直接输出的地面对流有效位能显著偏小,其值本身不能表征强对流环境对流有效位能大小,但变化趋势或可用。因温度参数代表性好,中层较之高、低层误差小,NCEP再分析资料计算的K指数和温度直减率能够较好地表征超级单体环境大气层结稳定度。

3) NCEP再分析资料提取的风廓线在自由大气层与探空资料差异小,用NCEP再分析资料计算得到的0~6 km风垂直切变与探空资料几乎一致,可用于强对流天气诊断;NCEP再分析资料700 hPa以上风向与探空资料基本一致。因此,NCEP再分析资料应能给出与探空资料基本一致的环流型,可用NCEP再分析资料进行雷暴环流背景分析。

4) 平均而言,NCEP再分析资料湿度廓线低层偏干而中、高层偏湿,风垂直切变略小,降低了灾害性强对流天气产生的概率。边界层风向与探空资料差异显著,湿度等参量在边界层误差较大,不宜用NCEP再分析资料边界层风分析对流启动问题。

参考文献
[1] 廖晓农, 俞小鼎, 谭一洲. 14时探空在改进北京地区对流天气潜势预报中的作用. 气象, 2007, 33, (3): 28–32. DOI:10.7519/j.issn.1000-0526.2007.03.004
[2] 何立富, 周庆亮, 陈涛, 等. 北京"7.10"暴雨中尺度对流系统分析. 应用气象学报, 2007, 18, (5): 655–664. DOI:10.11898/1001-7313.20070501
[3] 徐文慧, 倪允琪. 登陆台风环流内的一次中尺度强度对流过程. 应用气象学报, 2009, 20, (3): 267–275. DOI:10.11898/1001-7313.20090302
[4] 王瑾, 蒋建莹, 江吉喜. "7.18"济南突发性大暴雨特征. 应用气象学报, 2009, 20, (3): 295–302. DOI:10.11898/1001-7313.20090305
[5] 何立富, 周庆亮, 陈涛. "05. "05.6"华南暴雨中低纬度系统活动及其相互作用, 2010, 21, (4): 385–394.
[6] 谢建标, 林良勋, 颜文胜, 等. 广州2005年"3.22"强飑线天气过程分析. 应用气象学报, 2007, 18, (3): 321–329.
[7] 纪晓玲, 王式功, 穆建华, 等. 宁夏雷暴天气过程划分及环流分型和环境场特征. 应用气象学报, 2010, 21, (3): 329–334. DOI:10.11898/1001-7313.20100308
[8] 王秀明, 钟青. 环境与强对流 (暴) 云相互作用的个例模拟. 高原气象, 2009, 28, (2): 366–373.
[9] 慕熙昱, 徐琪, 夏文梅, 等. 准线形对流系统中强度-速度负相关性的特征研究. 气象学报, 2009, 67, (4): 631–639. DOI:10.11676/qxxb2009.063
[10] 廖晓农, 王华, 石增云. 北京地区雷暴大风日θe平均廓线特征. 气象, 2004, 30, (11): 35–37. DOI:10.3969/j.issn.1000-0526.2004.11.008
[11] 梁爱民, 张庆红, 申红喜, 等. 北京地区雷暴大风预报研究. 气象, 2006, 32, (11): 73–81. DOI:10.3969/j.issn.1000-0526.2006.11.012
[12] 郑永光, 张喜春, 陈炯, 等. 用NCEP资料分析华北暖季对流性天气的气候背景. 北京大学学报: 自然科学版, 2007, 43(5): 600-608.
[13] Thompson R L, Edwards R, Hart J A, et al. Close proximity soundings within supercell environments obtained from the Rapid Update Cycle. Wea Forecasting, 2003, 18: 1243–1261. DOI:10.1175/1520-0434(2003)018<1243:CPSWSE>2.0.CO;2
[14] Moller A R, Doswell III C A, Foster M P, et al. The operational recognition of supercell thunderstorm environments and storm structures. Wea Forecasting, 1994, 9: 327–347. DOI:10.1175/1520-0434(1994)009<0327:TOROST>2.0.CO;2
[15] [2011-02-01]. http://www.ofcm.gov/fmh3/text/chapter2.htm.
[16] 魏东, 尤凤春, 范水勇, 等. 北京快速更新循环预报系统_BJ_RUC_模式探空质量评估分析. 气象, 2010, 36, (8): 72–80. DOI:10.7519/j.issn.1000-0526.2010.08.010
[17] 李佳英, 俞小鼎, 王迎春. 用探空资料检验中尺度数值模式对强对流天气的诊断分析能力. 气象, 2006, 32, (7): 13–17. DOI:10.7519/j.issn.1000-0526.2006.07.002
[18] 陈子通, 闫敬华, 苏耀墀. 模式探空的评估分析及其在强对流天气预报中的应用研究. 大气科学, 2006, 32, (2): 235–247.