本研究使用的Parsivel2激光雨滴谱仪布设在重庆市沙坪坝气象站，由德国OTT公司制造，是一种现代化的以激光技术为基础的新型光学测量系统.主要由激光发射器、接收器、控制运算单元、存储器等部分组成，激光测量尺寸为180 mm×30 mm，可测量区域为54 cm²，液态降水粒径测量范围为0.2~8 mm，固态降水粒径测量范围为0.2~25 mm，速度测量范围为0.2~20 m/s，根据粒径大小及速度可推测出毛毛雨、小雨/雨、雨、雨夹雪、雪、米雪、冻雨、冰雹等8种降水类型，与气象观测员观测的降雨类型吻合度超过97%，因此能很好识别相应的天气现象.仪器采样时间为60 s，粒子直径和速度都分为32档，每1 min采样的粒子谱数据为32×32=1 024个，由于粒子直径的前两档精度不够，因此在使用时剔除，对于液态粒子下落时的形变问题本版本仪器在设计时已考虑订正，因此不再对液态降水粒子进行形变订正，但对固态降水粒子还需要根据类别进行订正.

微雨雷达布设在沙坪坝气象站，与激光雨滴谱仪同位置，该仪器由德国Metek公司生产制造，是一种小型的垂直指向雷达，采用连续调频(FW-CW)技术，通过雨滴的尺寸与雷达散射截面、雨滴的尺寸与下落速度之间的关系来判断降雨的状况及雨滴的大小分布，工作频率为24.230 GHz，可用于观测降水的垂直结构.本研究采用的微雨雷达探测数据垂直分辨率为100 m，共30层，时间分辨率为1 min，粒子直径测量范围为0.109~6 mm，速度测量范围为0~12.192 m/s，速度分辨率为0.190 5 m/s.仪器提供雷达反射率因子Z、雨强R、液态水含量LWC、雨滴下落速度W等物理量，主要物理量计算公式如下：

雷达反射率因子：

雨强：

液态水含量：

雨滴平均直径：

雨滴总数浓度：

本研究主要针对2016年1月22-24日的一次雨转雪天气过程进行研究，首先利用同站布设的自动站雨量、OTT雨滴谱仪观测雨量，验证微雨雷达对此次过程观测的可靠性，然后通过OTT观测的天气现象把降水过程分为不同时段，分析不同时段的MRR雷达反射率、雨强值、雨滴下落速度、粒子平均直径、谱变化等来研究降水的垂直分布特征.其中雨滴谱仪及MRR微雨雷达观测数据时间分辨率都为1 min，自动站雨量数据分辨率为1 h.

2016年1月19日夜间至24日，重庆地区出现了一次明显的降温降雨(雪)天气过程. 19日夜间至22日白天，大部地区以间断小雨为主，中部部分地区达中雨，城口、巫溪、巫山、黔江、酉阳、秀山等地小雪到中雪，酉阳有短暂冻雨；22日夜间至24日白天，重庆市中西部和东南部先后出现降雪，其中中西部普降中雪到大雪，局地暴雪，东南部小雪到中雪.沙坪坝站地面雨滴谱仪在23日01：21观测到雨夹雪，05：17开始转为中小雪，13：22后转为零星小雪、雨夹雪、小雨混合等，过程持续到24日2：30结束.从雷达回波来看，整个降雪时间段观测站点上空回波都在25 dBZ以内，整个观测范围内无较强回波，雨滴谱仪资料判断为稳定性降水.天气形势分析可知，此次重庆市大范围降雪过程产生的主要原因为东北冷涡西侧横槽转竖引导冷空气南下，且过程前期中低层储备了较为充沛的水汽. 1月21日20时，500 hPa蒙古东部有冷涡，冷中心达-48 ℃，冷涡底部有一横槽，延伸至新疆中部，青藏高原上有高原槽，重庆地区受西风波动气流影响；700 hPa重庆长江沿线以北有偏东风与偏西风的切变，水汽较为充沛，比湿达到4g/kg；850 hPa重庆地区为较弱的偏南风，温度露点差较小，0 ℃温度线位于四川南部—贵州中部—湖南中部. 1月23日08时实况500 hPa上我国东北地区的冷涡中心略有南压，横槽已经转竖，引导冷空气南下，南支槽压过重庆，重庆受槽后西北气流控制. 850 hPa受高压底部的偏东冷空气回流影响，温度-4℃线压至贵州中部地区，同时在我市西部地区形成较强辐合，有动力抬升作用，且前期中低层储备了较为充沛的水汽，有利于增大西部的降雪量.截至24日，中西部地区日平均气温下降了6~9 ℃，其余地区下降4~5 ℃，24日夜间全市大部地区最低气温在0 ℃以下，城口最低为-9.5 ℃.

前人验证微雨雷达探测效果时，通常选择微雨雷达近地面观测数据进行对比分析^{[14, 18]}，也有学者认为受地面湍流影响，最底层观测数据不可靠，选择更高层数据进行对比^[15].本研究将1月22-24日雨转雪过程MRR探测得到的100，200，300，400 m高度雨量数据与同时期OTT雨滴谱仪、自动站观测数据进行对比. 表 1为3类仪器观测雨量，在此次雨转雪过程中各类仪器累积过程雨量有较大差异.与王洪等^[15]、崔云扬^[19]研究不同，此次过程，微雨雷达观测雨量随着观测高度增加偏差越来越大，OTT雨滴谱仪在未进行订正前偏差亦较大，在进行雨雪订正后两者累积雨量较为接近. 图 1给出了整个过程中3类仪器观测的小时雨量及分钟雨强.可以看出，整个过程中自动站观测小时雨量较小，绝大多数时段小时雨量在1 mm以下，最大小时雨量为2.1 mm，出现在22日23时.以自动站雨量为基准作为参照，OTT雨滴谱仪在开始降雪后雨量偏差较大，对其进行雨雪订正后与自动站观测雨量较为接近，偏差较小；微雨雷达观测雨量与自动站偏差主要出现在地面出现降雪前，且高度越高，偏离越大，出现偏差的时间越早，在前期地面弱降水阶段和地面出现降雪阶段，微雨雷达与自动站两者之间趋势变化等较为一致.

为进一步探讨仪器间观测结果的差异，图 2为OTT与MRR近地面(100 m)高度层观测的雨滴谱、平均直径随时间变化图.从图 2a和图 2b看出，在降雪前，对于1 mm以下的雨滴数浓度，微雨雷达与雨滴谱仪观测较为一致，1 mm以上的雨滴数浓度雨滴谱仪观测较微雨雷达为多；降雪开始后，两者有较为明显的差异，雨滴谱仪观测数浓度变化不大，但粒子直径谱不断拓展，而微雨雷达观测雨滴数浓度出现量级上的跃增，比降雪前高1~3个量级，粒子直径谱逐渐收窄，随后趋于稳定. 图 2c为两种仪器观测的粒子平均直径变化图，可以看出雨滴平均直径在降雪前后有明显分化.产生两者仪器观测差异的主要原因可能在于仪器观测原理的不同，对粒子的观测中，雨滴谱仪通过激光束来实测粒子的直径和速度，微雨雷达则利用雨滴的大小与雷达散射截面、雨滴的大小与下落速度之间的关系来判断降雨的状况及雨滴大小分布，仪器导出速度参量后再根据直径速度经验公式来导出直径D，此经验公式适用于液态降水.对于固态降水中的雪花，其下落速度在1~2 m/s，根据经验公式，对应的粒子直径在0.5 mm以下，因此造成两者谱分布及粒子直径等特征出现差异.

OTT Parsivel雨滴谱仪能较好分辨出毛毛雨、小雨/雨、雨、雨夹雪、雪、米雪、冻雨、冰雹等8种降水类型，与气象观测员观测的降雨类型吻合度超过97%，能很好地识别相应的天气现象. 图 3d给出了2016年1月22日0：30至1月24日2：30间各类天气现象时间演变分布，可以看出，此次天气过程大致可以分为如下8个不同阶段：毛毛雨(22日00：43-17：40，期间夹杂短暂的小雨)、小雨—中雨(22日17：41-23日1：12)、雨夹雪(23日1：13-5：16)、中雪阶段(23日5：17-10：12)、小雪阶段(23日10：13-13：22)、零星小雪(23日13：23-17：13)、间歇期(23日17：14-20：18)、过程末尾的中雪—小雪—雨夹雪—小雨混合阶段(23日20：18-24日2：30). 图 3a~c分别给出了微雨雷达观测的雷达回波反射率Z、液态水含量LWC、粒子下落速度W的随时间演变剖面图，可以看出，整个降水过程雷达回波反射率较低，最大回波反射率为34 dBZ，在毛毛雨阶段回波垂直结构较为平稳，在20 dBZ以下，这一阶段空中液态水含量较低，在0.1 g/m³以下，但从粒子的下落速度来看，粒子的下落速度分布变化较大，前期存在大于4 m/s的粒子；小—中雨阶段，空中回波反射率增强，同一阶段空中液态水含量也较大，大于1 g/m³，回波反射率及液态水含量两者形态较为一致，其高值区底部向下倾斜，逐渐接地，从粒子下落速度来看，同时期液态水含量大值区底部以下粒子下落速度较大(液态水含量大值区粒子下落速度在1~2 m/s)；雨夹雪阶段垂直方向上液态水含量仍然维持大于1 g/m³，回波反射率在10~20 dBZ之间，粒子下落速度1~2 m/s；中雪阶段与雨夹雪阶段类似，变化不大；随着天气过程减弱，到达小雪阶段，液态水含量逐渐降低到0.5 g/m³以下，回波反射率逐渐下降到10 dBZ以下，低层粒子下落速度仍维持1~2 m/s，高层出现大量大于3 m/s的降水粒子；间歇期地面未观测到降水，回波反射率较低，在12~14.7 dBZ之间，各层液态水含量较少；在过程后期，液态水含量再次增大，回波强度增大到15 dBZ左右不再增长，粒子下落速度仍在1~2 m/s，地面降水形态逐渐由雪、雨夹雪向小雨、毛毛雨转变(固态转为液态)，随后过程结束.

从各层粒子随时间演变图(图 4)来看，1 km以上粒子变化不大，形态较为相似，粒子数浓度高值区主要集中在0.5 mm以下，高值区达108 mm/m³.从近地面(100 m)DSD时间演变图(图 5d)来看，毛毛雨阶段粒子直径主要集中在1 mm以下，随着冷空气的加强，降水向中雨转变，粒子平均直径也逐渐增大，向2 mm延伸.降水转为雨夹雪后，粒子直径逐渐减少，且次数粒子数浓度迅速增大1~2个量级.转为纯雪阶段后，粒子数浓度及直径分布与高层相似，高值区主要集中在0.5 mm以下.与图 3对应来看，高空(1 km以上)粒子的数浓度急剧增加的时段在17时后，此时对应于高空液态水含量大值区和回波反射率大值区，1 km以下的DSD分布来看(图略)，粒子数浓度高值区同样存在倾斜特征，一直到地面产生雨夹雪(1：13)，因此，此大值区的出现表明高空粒子已从液态转为固态.

对整个雨转雪过程不同阶段的雷达反射率因子、粒子下落速度分布特征进行分析，计算对应的频率分布及平均垂直廓线(图 5、图 6)，频率分布图中的蓝色虚线曲线为对应物理量的平均垂直廓线，概率分布为同一高度层按一定间隔分布统计所得，全过程样本数为3 001×30个.

从图 5a~h各阶段下不同高度层的雷达反射率频率分布来看，在毛毛雨阶段(图 5a)，雷达反射率Z在1 km以下主要集中在0~10 dBZ之间，1 km以上主要集中在5 dBZ以下；小到中雨阶段(图 5b)，随高度降低Z逐渐增加，1 km以上主要集中0~20 dBZ之间，平均值在10 dBZ附近，1 km以下有所拓展，在10~20 dBZ之间，平均值从10 dBZ增长到20 dBZ；雨夹雪阶段(图 5c)，Z随高度降低逐渐增加，但增长幅度不大，Z的频率分布逐渐趋于向10~20 dBZ区间段集中；转为纯中雪阶段后(图 5d)，Z随高度变化不大，频率分布较为均匀，主要集中在5~20 dBZ之间，平均值为10 dBZ左右；中雪转为小雪后(图 5e)，Z随高度降低呈现波动变化，2~3 km随高度逐渐减小，平均值小于0 dBZ，2 km以下逐渐增大，各层频率分布较中雪阶段更为集中；随着天气过程继续减弱为零星小雪(图 5f)，回波反射率进一步降低，各层Z值在10 dBZ以下，平均值1 km以上在1 dBZ左右，1 km以下在5 dBZ以下，且在近地面逐渐降低；图 5g为短暂间歇期，地面未观测到明显降水，1 km以下Z平均值在0 dBZ以下，但1KM以上仍有较弱回波，50%以上的样本分布在5 dBZ附近；在降水末期(图 5h)，各层Z频率再次趋于均匀分布，地面以小雨、毛毛雨、小雪等多种天气混杂，各层回波反射率平均值变化不大，在2~3 dBZ之间；整个过程来看(图 5i)，Z的平均值在0~10 dBZ之间，1 km以前变化不大，1 km以下有所增加，主要原因是中雨阶段粒子平均直径增长明显，Z值大幅增大，从频率分布来看，2.8 km以下频率分布较为平均.

图 6为不同阶段的各层粒子下落速度频率分布图，蓝色虚线为平均速度廓线.可以看出，毛毛雨阶段(图 6a)粒子平均下落速度在1~2 m/s之间，1 km以上频率分布较为均匀，均值变化不大，1 km以下速度有拓宽趋势；小到中雨阶段(图 6b)，1 km以上粒子下落速度进一步向1 m/s附近集中，1 km以下速度最大超过5 m/s；随后的阶段(图 6c~h)开始降雪(雨夹雪、小雪、中雪等)，各层频率分布变化不大，基本集中在1 m/s附近，各高度层变化不大，差异之处在于，雨夹雪阶段低层略大于1 m/s，小雪转零星小雪阶段(图 6e)高层速度有较大变化，间歇期低层主要集中在1 m/s以下；整个过程来看(图 4i)，各层粒子速度变化不大，频率分布较为集中，主要在1 m/s附近. Nakaya等^[20]研究发现，雪花下落速度和尺度之间的关系主要取决于雪花的形状和它是否结霜，各种平板状冰晶速度接近于0.3 m/s左右，雪粒聚集物为0.5 m/s，立体枝状冰晶为0.57 m/s.当冻滴附着在冰晶上时，下落速度增至1 m/s，此时速度随尺度增大.可见此阶段粒子多以雪粒聚集物和立体枝状冰晶为主.此次过程中(图 6b)1 km以上各层粒子下落速度平均分布在1 m/s左右，可能表征了高空已经完成雨向雪的转变.

图 7为不同降水阶段的液态水含量LWC、雨强R、雨滴数浓度N₀以及雨滴平均直径D平均垂直廓线.从液态水含量LWC、雨强R的平均垂直廓线来看，两者具有类似特征，1 km以上变化不大，400~1 000 m高度层液态水含量有所增加，400 m以下随高度降低，从降水阶段的平均廓线来看，毛毛雨、中雪、小雪、零星小雪等阶段液态水含量、雨强在垂直高度上变化不大.天气较强的中雨阶段，LWC、R变化较为明显，在垂直高度上与整个过程的平均廓线特征较为一致；雨滴数浓度(图 7c)在各降水阶段变化特征较为明显，1 km以上各层平均数浓度随高度降低逐渐增大，1 km以下平均数浓度逐渐减小，毛毛雨阶段雨滴数浓度最小，lg(N₀)值在2~4之间，小、中雨阶段数浓度变化最大，平均值在3~5之间，小雪阶段数浓度变化梯度最大，高层数浓度与小雨阶段接近，低层与雨夹雪阶段较为接近.总的来看，随着降雪阶段(固态降水)的不断演变，粒子数浓度存在逐渐降低的趋势；图 7d表明，此次降水过程平均雨滴直径整体较小，前期(毛毛雨、小中雨)随着天气加强，平均半径随高度降低有所增加，转为降雪后各高度层平均粒子半径变化不大，但在小雪阶段高层出现了粒子半径拓宽的现象，此次粒子数浓度剧减，结合粒子下落速度来看有明显增长，表明此高度层粒子可能以液态粒子为主，一定程度上揭示了高空粒子相态的转变.

雨滴谱分布特征是了解降水的关键，为进一步分析雨滴在垂直方向上的分布特征，对不同降水阶段的雨滴谱进行平均分析. 图 8给出了此次雨转雪天气过程中不同降水阶段的平均雨滴谱随高度变化情况.可以看出，在毛毛雨阶段(图 8a)粒子直径主要集中在1 mm以下，数浓度最大量级为10⁴，粒子数浓度在垂直方向上变化不大；小、中雨阶段(图 8b)，1 km以上小粒子数浓度变化不大，但量级较毛毛雨阶段有明显增长，达到10⁶.中等尺度、大粒子数浓度随高度降低有所减少. 1 km以下小粒子数浓度变化相对不大，中等尺度、大粒子数浓度随高度降低增长较大，结合雨滴下落速度来看，可能表明此阶段有部分中等尺度、大尺度粒子转化为雪晶，粒子聚合效应明显.而在1 km粒子仍以液态水存在为主，且雨滴间的碰并增长效应较强；降雪各阶段(图 8c~h)雨滴平均谱变化不大，各阶段谱型较为类似，小粒子数浓度在垂直方向上变化不大，中、大尺度粒子随着高度降低略微减少，在2 km以下逐渐趋于平稳，粒子最大数浓度在降雪前期维持在10⁶的量级上，但在降雪后期(图 8f~h)最大值逐渐降低到10⁴量级，整个过程来看(图 8i)，小粒子在垂直高度上变化不大，中、大尺度粒子数浓度随高度降低有先减小后增大的趋势.

综上，降水不同过程粒子谱有着不同的垂直分布特征，前期为稳定弱降水阶段，各层粒子浓度变化不大，转为较强的中雨阶段后，低层雨滴碰并增长明显，随着天气由雨转为雨夹雪，各层雨滴分布逐渐趋于平稳，完全转雪后，各层粒子浓度变化不大，一直到降雪后期.后期过程结束主要以雨滴数浓度减少为主，粒子尺度分布较为稳定.

利用布设在沙坪坝站的雨滴谱仪和微雨雷达获取的连续观测资料，对2016年1月22-24日的一次雨转雪天气过程降水微物理特征进行分析，得出如下结论：

1) 此次天气过程持续时间长达49 h 47 min，其中纯液态降水阶段24 h 29 min，累积降雨量13 mm，累积降雪(含雨夹雪)时长25 h 18 min，降水量13.8 mm.粒子的相态变化与冷空气活动有密切关系.相比雨滴谱仪观测雨强，两者趋势较为一致.

2) 微雨雷达近地面观测物理量与地面雨滴谱仪观测物理量对比表明，两者在对雨强的观测上较为一致，但在雨夹雪阶段存在较大差异，这是由于雨夹雪阶段雨滴谱仪修正公式不同造成的.对粒子谱的观测来看，在纯降雪阶段，受观测原理影响，微雨雷达和雨滴谱仪观测谱有明显不同.

3) 液态降水阶段的雨滴平均特征值和滴谱的垂直特征与固态降水(雪)阶段差异较大，液态降水表现为垂直下落方向上先减小后增大的特征，粒子碰并增长特征明显，固态降水(雪)垂直方向上则较为稳定.

4) 不同降水阶段的雨滴特征量变化差异较大.在毛毛雨阶段，各特征量在垂直方向上变化不大，粒子数浓度较小.在中雨阶段，粒子下落速度、粒子平均直径在高空已经表现为雪花的特征，低层随高度降低而增加，表现为液态粒子碰并增长的特征，其特征量表现为500 m以上(以下)随高度降低增大(减小)的特征.在降雪各阶段，各特征量变化较为一致.

5) 借助微雨雷达以速度来反演粒子大小这一特征，可以较清晰地推断垂直方向上粒子相态的转变.此次过程地面观测到中等强度降水时，高空粒子下落速度均匀分布在1 m/s附近、回波反射率Z分布在10 dBZ、粒子数浓度急剧增加等物理量特征表明高空以冰雪晶为主，随着时间演变，这一特征向地面倾斜，逐渐接地，地面天气现象转变为雨夹雪；而小雪阶段高空粒子相态的突然转变，可能预示天气过程的转变.这一特征可较好地应用于对于降雪的临近及结束预报等.重庆地区微雨雷达建成时间较短，且观测高度仅为地面到3 km高度，观测数据完整且高质量的个例分析较少，针对本地区的观测研究结果有待进一步分析，为相关资料在人工影响天气、天气预报中的应用累积经验.