基于多源观测资料对张家口一次雨雪天气降水相态特征的分析

姬雪帅，王丽婧，郭 宏，康博思，黄山江，张曦丹，郭旭晖

（1. 河北省张家口市气象局，河北 张家口 075000；2. 宁夏中卫市气象局，宁夏 中卫 755000）

引 言

冬季降水相态比较复杂，主要有雨、雨夹雪、雪、冰粒和冻雨等，其中液态和固态降水的转换一直是预报中的难点，特别是在秋冬和冬春季转换期间容易出现降水相态转变，相同降水量级对应不同降水相态时会造成不同影响。对于气象服务而言，降水量达10 mm，相态为雨时刚达中雨量级，对城市交通、社会生产和人民生活不会产生较大影响，但相态为雪时则达暴雪量级，对城市交通等会带来很大影响。因此降水相态的准确预报对气象服务质量与效果有决定性作用。张家口地区冬季相对较长，季节交替时降水过程通常会伴随相态转换，特别当强冷空气活动频繁时，前期水汽积累充足，后期降温幅度大，降水相态转变天气经常发生，降水相态及其转换时间预报是当地天气预报的难点和重点。

有关雨雪相态转换的研究由来已久，目前已经取得了很多研究成果［1-3］。如通过分析不同地区气象要素分布与锋生变化特征，归纳不同降水相态下各高度层温度阈值，确定本地化的雨雪相态判定指标，建立适宜于本地的降水相态预报方程［4-16］。随着对已有判据的不断使用和对降水相态天气特征的不断认识，位势高度和相对湿度等也逐渐用于降水相态预报方程的建立［16-17］。

受观测资料限制，对降水相态转变的云微物理结构认识还不够深刻。近年来随着风廓线雷达、微波辐射计、云雷达等探测设备增加，多源探测资料在雨雪相态分析中逐渐得到了应用［18-26］。微波辐射计、云雷达等新型探测设备可以探测雨雪相态转换期间的大气层结和云微物理结构，资料时空分辨率高，对降水相态转换特征有一定的指示意义。

2020 年11 月17—19 日，张家口地区出现了一次雨转雪天气过程，降水过程以雨开始、以雪结束。本文利用多源观测资料对此次过程降水相态变化特征进行分析，探索新资料在雨雪相态预报中的应用，以期加深对张家口地区雨雪相态转换特征的认识，为张家口地区雨雪相态预报提供参考依据。

1 仪器介绍与资料说明

1.1 仪器介绍

毫米波云雷达（CR）和微波辐射计（MR）位于河北省张家口观测站（114.53°E，40.47°N），SA双偏振多普勒雷达位于张家口康保（114.54°E，41.47°N）。

毫米波云雷达型号为KPS-HMB，全固态单通道体制，Ka 波段（35.5 GHz），可连续探测垂直方向的回波强度、径向速度及速度谱宽，主要用于分析云宏微观结构，性能参数如表1所示。

表1 云雷达主要参数Tab.1 Key parameters of cloud radar

微波辐射计为美国RADIOMETRICS 公司生产的35 通道MP-3000A 型微波辐射计，提供从地面到10 km 高度内温度、相对湿度、水汽密度和液态水廓线数据，用于分析探测环境的温度、湿度场。

天气雷达为北京敏视达公司生产的S波段双偏振全相参多普勒天气雷达（简称“SA 双偏振多普勒雷达”），可以提供种类更多、质量更高的基数据和产品，除包含基本反射率因子、基本速度和谱宽外，还包含差分反射率因子、零滞后相关系数和差分相位移等，性能参数如表2所示。

表2 天气雷达主要参数Tab.2 Key parameters of weather radar

1.2 资料说明

（1）欧洲中心提供的全球再分析资料（简称“ERA5”资料），时间间隔为1 h，水平分辨率为0.25°×0.25°，该资料用于环流背景和物理量诊断。

（2）探空资料为张家口探空站秒级探空数据，时间分辨率为1 s，每天08：00（北京时，下同）和20：00 进行观测，用于对比分析微波辐射计和ERA5的数据。

（3）云雷达数据主要使用17 日23：20 至18 日07：20和18日23：20至19日07：20两个时段数据，主要用于分析反射率因子、液态水含量和基本速度特征。

（4）微波辐射计数据主要用于分析17、18 日温度、水汽密度、积分水汽含量和积分液态水含量特征，该设备垂直方向为58 层，0~0.5 km 高度分辨率为0.05 km，0.5~2 km 高度分辨率为0.1 km，2~10 km高度分辨率为0.25 km。

（5）SA 双偏振多普勒雷达资料时间分辨率为6 min，主要用于分析0 ℃层亮带变化和粒子相态识别。

2 天气实况与环流背景

2020年11月17日夜间至19日上午河北省张家口市出现全市范围雨雪天气，降水相态由雨转雪，最大降水（33.5 mm）出现在崇礼区。张家口站17日22：00—24：00 为微量降水，19 日07：00 降水过程结束，过程降水量为24.0 mm，小时最大雨量为2.60 mm，其中降雨22.7 mm，降雪1.3 mm；17—18日为降雨过程，19日凌晨为降雪过程（图1）。

图1 2020年11月17日20：00至19日08：00张家口站地面降水逐时变化Fig.1 Hourly variation of surface precipitation at Zhangjiakou station from 20：00 BST 17 to 08：00 BST 19 November 2020

从图2（a）看出，17 日20：00，500 hPa 河北中东部、渤海以及东北地区为高压暖脊控制，高空槽位于内蒙古中西部地区，温度槽落后于高度槽，该高空槽向东移动时不断发展，此时张家口地区处于高空槽前西南气流中，相对湿度逐渐增大；700 hPa（图略），河北中东部有12~14 m·s-1急流存在，张家口处于急流左前侧，有利于水汽的输送和辐合；850 hPa［图2（b）］，山西省内存在明显辐合，山西南部出现风场辐合中心，整个河北地区有风速辐合，张家口温度在0 ℃以上，并处于暖平流中。随着系统逐渐东移，张家口17日夜间开始产生降水，降水相态为雨。

图2 2020年11月17日20：00 500 hPa（a）、850 hPa（b）和18日20：00 500 hPa（c）、850 hPa（d）高度场（黑色等值线，单位：dagpm）、温度场（红色等值线，单位：℃）、相对湿度（填色区，单位：%）及风场（风矢，单位：m·s-1）（五角星位置为张家口）Fig.2 Height field（black isolines，Unit：dagpm），temperature field（red isolines，Unit：℃），relative humidity（color shaded，Unit：%）and wind field（wind vectors，Unit：m·s-1）at 500 hPa（a）and 850 hPa（b）at 20：00 BST 17 November 2020，and 500 hPa（c），850 hPa（d）at 20：00 BST 18 November 2020（The location of the red star is Zhangjiakou）

18 日20：00，500 hPa［图2（c）］高空槽继续东移，700 hPa（图略）和850 hPa［图2（d）］出现低涡，低涡中心位于渤海湾，河北中西部处于涡后，风场与等温线近乎垂直，冷平流强，张家口地区降温明显，850 hPa 温度场上张家口大部分地区温度在0 ℃以下，为18 日后半夜相态转换奠定基础。19 日03：00 850 hPa 温度场（图略）上张家口大部分地区气温低于-4 ℃，降水相态已经发生转换。

从地面图（图略）分析，从17 日20：00 至19 日08：00 地面低压中心位于山东，张家口地区主要受地面倒槽影响，有利于辐合上升运动发展产生降水，19日以后转为冷高压控制，降水结束。

3 热动力结构与微物理特征

利用ERA5 再分析资料，对张家口站上空的温度、湿度、云中冰相粒子、云中液相粒子、散度和垂直速度特征进行分析。从张家口站的小时降水量（图1）可以看出，降水分为2 个时段，17 日22：00 至18日18：00为降雨时段，19 日01：00—07：00为降雪时段。

首先利用张家口探空站秒级探空资料对ERA5再分析资料的温、湿度进行评估［图3（a）］，通过对比18 日08：00 的各层温度，发现探空数据与ERA5再分析资料的温度误差小于1 ℃，变化趋势一致；500 hPa 以下各层相对湿度误差在10%左右。因此可以认为本次过程ERA5 再分析资料数据可以用来分析张家口地区的温度与湿度特征。

图3 2020年11月18日08：00 ERA5再分析资料与探空资料的温度、湿度散点图（a），ERA5再分析资料温度、相对湿度的时间-高度剖面（绿线为相对湿度，单位：%；红线为温度，单位：℃；填色区域表示温度-18~-12 ℃范围）（b），云中冰相粒子和液相粒子含量的时间-高度剖面（红线为温度，单位：℃；绿线为云中液态水粒子含量，黑线为云中冰相粒子含量，单位：10-2 g·kg-1）（c），散度和垂直速度的时间-高度剖面（蓝线为散度，单位：10-5 s-1；红线为垂直速度，单位：Pa·s-1）（d）Fig.3 The scattering plots of temperature and relative humidity between ERA5 reanalysis data and the sounding data at 08：00 BST 18 November 2020（a），the time-height profile of temperature and relative humidity from ERA5（the green line for relative humidity，Unit：%；the red line for temperature，Unit：℃；the shaded for temperature ranged from-18 to-12 ℃）（b），the time-height profile of cloud ice particles and cloud liquid particles content（the red line for temperature，Unit：℃；the green line for cloud liquid particles content and the black line for cloud ice particles content，Unit：10-2 g·kg-1）（c），the time-height profile of divergence and vertical velocity（the blue line for divergence，Unit：10-5 s-1；the red line for vertical velocity，Unit：Pa·s-1）（d）

从图3（b）可以看出，第一阶段降雨开始时，0 ℃等温线位于700 hPa 附近，近地面温度大于0 ℃，而-18~-12 ℃温度范围所在高度基本处于600~500 hPa 之内，而-18～-12 ℃温度范围对树枝状雪花的形成至关重要，雪花形成以后，降落过程中融化层厚度比较深厚，降落到地面为液态降水；第二阶段降雪开始时，整个大气层温度基本处于0 ℃以下，且-18~-12 ℃之间厚度增加，有利于雪花形成，雪花形成后降落过程中都处于0 ℃以下环境中，降落到地面后为固态降水。从相对湿度分布来看，第一阶段降雨时，从500 hPa 到近地面相对湿度都在90%以上；第二阶段降雪时，相对湿度在90%以上的厚度明显要薄且近地面相对湿度在70%以下，因此有利于降雪粒子相态的保持。

HOBBS 等［27］研究表明：温度高于-10 ℃，降水粒子以过冷水滴为主；温度低于-40 ℃，则基本以冰相粒子为主。由图3（c）可知，第一阶段降雨时，云体的伸展高度在-40 ℃等温线所在高度以上，云体上部基本以冰相粒子为主，云体下部位于-10 ℃等温线所在高度以下，700 hPa 以下都处于0 ℃等温线所在高度以下，暖层深厚，因此云体内的降水粒子降落到地面后形成液态降水；第二阶段降雪时，云体伸展高度降低，云体上部600~500 hPa 以冰相粒子为主，850~600 hPa为过冷水滴，相比于第一阶段，冰相粒子和液相粒子含量均有明显下降，降水效率降低。

从图3（d）可以看出，第一阶段降雨时，负散度中心位于650 hPa附近，强度为-3×10-5s-1，正散度中心位于300 hPa 附近，无辐散层位于500 hPa 附近，从低层到高层呈现“正、负、正”分布，上升运动主要位于中高层，最大上升速度达-0.6 Pa·s-1，最大上升速度所在高度位于500 hPa，表明降雨阶段的动力强迫主要位于高层；第二阶段降雪时，负散度中心位于850 hPa 附近，强度为-9×10-5s-1，正散度中心位于700 hPa附近，无辐散层位于800 hPa附近，上升运动主要位于低层，最大上升速度为-0.4 Pa·s-1，其所在高度位于850 hPa 附近，表明降雪阶段的动力强迫主要位于低层。

4 多源观测获取的雨雪转换特征

4.1 云雷达回波特征与液态水含量演变

从Ka 波段云雷达观测（图4）发现，17 日23：40开始，反射率因子从低层到高层开始增强，融化层高度在2.1~2.4 km；降雨初期［图4（a）］，18 日00：00—01：00，高反射率因子质心下降，18 日00：24，大于10 dBZ 的质心及地，地面出现降水，夜间出现2 个较强降水时段，分别为01：00—02：10 和05：00—06：10，2 个阶段云雷达反射率因子均出现大于10 dBZ 大值区，与地面自动站降水观测一致。18 日10：00—14：00（图略），也出现强质心及地现象，5 h降水量达10.6 mm。从液态水含量［图4（c）］来看，冷层液态水含量为0.1~0.5 g·m-3，降水较强时，暖层液态水含量达0.5~1.0 g·m-3。液态水含量大于0.5 g·m-3大值区与反射率因子大于10 dBZ 大值区分布相一致。从基本速度［图4（e）］可以看出，降雨时段最大速度可达6.0~8.0 m·s-1。降雪时段［图4（b）］，19 日01：00—06：00，01：00 在3.0~3.5 km 高度出现强质心，最大反射率因子在30 dBZ左右，云顶高度超过6.0 km，降雪强度最大达0.4 mm·h-1，降雪结束以后云顶高度下降，云顶最大高度在3.0 km；云雷达反射率因子大于10 dBZ 质心及地时间与降雪持续时间一致。降雪时段液态水含量在0.5~1.0 g·m-3［图4（d）］，同时液态水含量大于0.5 g·m-3的大值区与反射率因子大于10 dBZ 大值区分布一致，基本速度在2 m·s-1以下［图4（f）］，较降雨时段大大减小。

可见，云雷达高时空分辨率资料能对降雨降雪系统进行精细化探测，反射率因子大于10 dBZ 大值区能反映云系中液态水含量大于0.5 g·m-3的大值区，同时反射率因子大于10 dBZ 质心及地持续时间与降水持续时间呈现明显正相关，降雨阶段垂直速度大于降雪阶段。

4.2 微波辐射计温度与水汽特征演变

从17日20：00微波辐射计与张家口探空站秒级探空资料反演的温度、湿度廓线［图5（a）、图5（b）］可以看出，8 km 以下两者温度变化趋势较一致，微波辐射计反演的温度比探空温度偏低；而7 km 以上两者相对湿度差距较大，但5 km 以下相对湿度变化趋势整体较一致，差值在15%以内，因此可以认为该微波辐射计资料的参考性比较大。降雨阶段［图5（c）］，从温度廓线可以看出，0 ℃层在2.0~3.0 km 之间，稍高于云雷达观测到的亮带位置，因为冰晶越过0 ℃层在暖层的融化需要一定时间，因此可认为云雷达与微波辐射计的观测比较一致，并且从18 日00：00—15：00，0 ℃层高度波动不大，15：00 以后温度迅速下降，0 ℃层高度及地；从水汽密度廓线可以看出，17：00（图略），低层水汽密度出现跃增，23：00 高层开始出现大值区，降雨时大值区分为两层，高层位于5.0~7.0 km 之间，为云中冰相粒子与液态水相粒子的混合物，最大值为6.0~8.0 g·m-3；低层位于2.0 km 以下高度，主要为液态水相粒子，与ERA5 再分析资料产品保持一致，最大值为10.0 g·m-3。降雪阶段［图5（d）］，从温度廓线可以看出，整层大气基本处于0 ℃以下；降雪阶段的水汽主要集中在2.0 km 以下，水汽密度较降雨阶段明显降低，最大值为4.0 g·m-3。从张家口微波辐射计反演的积分水汽含量和积分液态水含量分析可知，17 日17：00—24：00（图略），积分水汽含量从1.5 cm 升高至3.0 cm，降雨开始时积分水汽含量为2.0 cm；积分液态水含量在19：00 出现一次峰值并且在17：00—19：00 出现跃升，从0.2 mm 跃升至0.8 mm，较降雨开始时刻提前3 h 左右。降雨时段［图5（e）］，积分水汽含量和积分液态水含量有3 次跃增，与3次较强降水时刻分别对应，降雨时积分水汽含量为2.0~3.0 cm，积分液态水含量为2.0~4.0 mm，15：00 以后积分水汽含量迅速下降至2.0 cm 以下，积分液态水含量迅速下降至1.0 mm以下，降雨结束；降雪时段［图5（f）］，积分水汽含量在0.6~1.0 cm，积分液态水含量在0~0.2 mm。

图5 2020年11月17日20：00微波辐射计与探空温度（a）及湿度（b）廓线，18日（c、e）和19日（d、f）温度（上，单位：℃）与水汽密度（下，单位：g·m-3）（c、d）及积分水汽含量、积分液态水含量（e、f）Fig.5 The temperature（a）and relative humidity（b）profiles from the microwave radiometer and the sounding data at 20：00 BST 17 November 2020，the temperature（the top，Unit：℃）and water vapor density（the bottom，Unit：g·m-3）（c，d）and integrated water vapor content and integrated liquid water content（e，f）on 18（c，e）and 19（d，f）November 2020

综上所述，微波辐射计反演的温湿廓线、积分水汽含量和积分液态水含量在雨雪过程的2个阶段均有明显变化，根据大气整层温度变化对云体结构判断，特别是0 ℃层高度迅速降低并及地是雨转雪的重要判据，此外雨转雪过程中积分液态水含量迅速降低。同时积分水汽含量的迅速降低对降水结束有重要指示意义。

4.3 SA 双偏振多普勒雷达0 ℃层亮带特征与粒子分类识别

张家口SA 双偏振多普勒雷达3.3°仰角反射率［图6（a）］显示：2020 年11 月17 日22：00 降水开始时，整个张家口地区为均匀的层状云降水回波，出现降水的站点均为小雨，并且出现了比较完整的0 ℃层亮带；从反射率因子垂直剖面［图6（b）］可以看出，0 ℃层高度在2.3 km，与云雷达和微波辐射计的观测基本一致。图6（c）黑色方框区域是影响张家口站的云团回波，其强度在20~40 dBZ，与云雷达的观测结果基本一致。通过双偏振雷达的水凝物分类对粒子进行分类［图6（d）］，发现此云系主要以冰晶与聚合物为主，结合图5（c）微波辐射计温度分析可知，02：00 从地面到2.0 km 高度温度都在0 ℃以下，冰晶和聚合物在0 ℃以下的温度层结中仍然会保持固态降落到地面，降水相态转为雪。

图6 2020年11月17日22：00 3.3°仰角反射率因子（a）与沿41.375°N的反射率因子剖面（b）及19日02：00 0.5°仰角反射率因子（c)(单位：dBZ）与水凝物分类图（d）Fig.6 The reflectivity on 3.3°elevation(a)and the profile of reflectivity along 41.375°N(b)at 22:00 BST 17 November 2020 and the reflectivity on 0.5°elevation(c)(Unit:dBZ),the hydrometeor classification(d)at 02:00 BST 19 November 2020

综上所述，云雷达、微波辐射计和SA 双偏振雷达对0 ℃层的观测结果基本一致，而0 ℃层亮带的变化对雨雪天气转换有较好的指示意义；云雷达能对降雨降雪系统垂直结构进行精细化探测，能清晰揭示其发展演变过程，指示降水持续时间；微波辐射计温湿廓线演变对降水的起始时间、结束时间和降水相态均有指示意义；SA 双偏振多普勒雷达探测范围广，可以探测大范围云系的演变，但精细化程度较差，与云雷达结合探测可以对降水起始时间、降水强度精准判断，与微波辐射计的垂直温度结合可以进行粒子相态分类，对降水相态演变精准判断。因此多源观测资料结合可以更好地了解降水云系结构和大气的温湿变化，为云模式发展和人工影响天气提供重要的参考依据。

5 结 论

以张家口2020 年11 月17—19 日雨雪过程为例，应用常规气象观测资料、云雷达资料、微波辐射计资料和SA双偏振多普勒雷达等多源观测资料，分析了降水相态特征。主要结论如下：

（1）此次过程是在高空低槽、中低层低涡与地面倒槽配合下，高空槽后西北气流引导冷空气南下造成气温迅速下降导致降水相态从雨转为雪。

（2）热动力结构与微物理特征诊断表明：降雨时段的动力强迫主要位于高层，-18~-12 ℃温度层高度高，暖层深厚，低层相对湿度大，云体伸展高度高，云中冰相粒子和液相粒子含量都比较高；降雪时段的动力强迫主要位于低层，-18~-12 ℃温度层厚度大，云中冰相和液相粒子含量明显下降，雨雪转换时0 ℃层高度快速下降。

（3）云雷达高分辨率资料可以反映云体中0 ℃层变化，反射率因子大于10 dBZ 大值区能反映云系中液态水含量大于0.5 g·m-3大值区，同时反射率因子大于10 dBZ 质心及地持续时间与降水持续时间呈现明显正相关。

（4）微波辐射计的温湿廓线可以判断整层大气层结结构，大气温度的迅速降低对雨雪转换、积分水汽含量和积分液态水含量的变化对降水强度变化均有重要指示意义。

（5）利用SA双偏振多普勒天气雷达实现对水凝物分类，结合微波辐射计的温湿廓线可以对雨雪进行直观判断，从而实现对雨雪相态转换的短临订正。

本文重点分析不同降水相态时多种新型观测资料的演变特征，对0 ℃层的判断可以看出云雷达、微波辐射计及SA 双偏振多普勒雷达资料的反演结果比较一致，说明精细化的探测资料在冬季短临预报中可以应用。值得强调的是，这里有关的多源探测资料只针对张家口本地设备，所得结论也仅由本个例研究确定，普适性的规律需要收集更多个例进一步研究验证。