重庆微波辐射计资料的评估和探测特征分析

赵美艳，余 君，蒋 镇

（1.重庆市气象信息与技术保障中心，重庆 401147；2.丰都县气象局，重庆 408200；3.云阳县气象局，重庆 404500）

随着气象探测技术的不断进步，探空资料的来源日益增多，由气球携带随风漂移的无线电探空、测风逐渐发展到由飞机、风廓线雷达、地基水汽遥感探测、卫星等构成的综合探空系统。20世纪80年代以来，在探测大气垂直结构的技术领域方面，结构简单、无人值守的被动微波遥感技术发展迅猛[1-3]，尤其是地基微波辐射计，它利用大气在一定频率带中的微波辐射进行测量，能够探测得到分钟量级的气温、水汽等气象要素的廓线信息，已逐渐成为常规探空的重要补充。国内外已有众多专家学者对其探测原理、数据产品的反演、偏差分析等开展了研究和探讨[4-9]。将微波辐射计反演结果与GPS无线电探空和GPS/MET观测对比分析发现，微波辐射计探测的温、湿度廓线与GPS无线电探空结果具有很好的正相关，但其可降水量较GPS/MET偏大[10]。韩珏靖等[11]研究表明，在降水和强对流天气过程中，微波辐射计相对湿度有提前增大的现象，提前时间在强对流天气时要短于降水过程，对天气预报有重要指示意义。

伴随着微波辐射计各种反演技术的发展[12-13]，其观测资料的种类和数量逐渐增多，但其反演产品质量会受到地点、天气等各种因素的影响[14]。以往的研究多采用探空数据对微波辐射计廓线信息的准确性进行分析[15]，但由于受地域的影响和需要花费较高的操作和维护费用，探空观测覆盖范围有限，且观测时次一般每天2次（最多4次），远不能满足气象业务发展的需求。近年来，随着全球导航卫星系统（Global Navigation Satellite Systems，GNSS）的发展，美国GPS（Global Position System）、欧洲Galileo和中国的北斗等，为气象业务和研究提供了高时空分辨率的卫星遥感资料[16-17]。GPS无线电掩星观测（即空基GPS气象学）是基于GPS的探测技术，其无线电信号能够透过厚云层和降水区[18]，探测资料受云雨影响较小，且具有自校准、高精度、低成本等优点，被广泛应用于气象研究和业务应用[19-23]。因此，本研究拟利用Metop-A掩星资料和L波段探空资料分别对重庆沙坪坝MP-3000A型地基微波辐射计探测性能和质量进行统计分析，并针对具体天气个例，对天气过程中的探测特征进行分析，进一步发挥该资料在气象预测预报中的作用。

1 资料和方法

1.1 资料

1.1.1微波辐射计

本文所用微波辐射计资料为沙坪坝站2017年1月—2019年12月MP-3000A型地基微波辐射计观测资料，可反演输出天顶方向、北斜20°、南斜20°及南北平均等多个方向温度、相对湿度、水汽等廓线信息（本文分析为天顶方向），垂直方向从地面至10 km高度共58层，时间分辨率为3 min。

1.1.2Metop-A掩星资料

GPS掩星是一种主动遥感大气探测技术，其利用GPS卫星发射的电磁波信号在穿越地球大气时出现的路径弯曲和信号延迟等信息反演大气温度、压强和水汽等廓线信息[24]。Metop-A掩星为欧洲气象卫星应用组织于2006年发射，每天可提供600条左右大气廓线，其折射率双权重平均值和标准差与COSMIC等几乎无差异，尤其是10 km以下，其折射率偏差较小[25]。本研究使用的Metop-A掩星资料来自于COSMIC数据分析与存储中心（the COSMIC Data Analysis and Archival Center，CDAAC），主要使用掩星反演的Wetprf数据产品，垂直分辨率为200 m，高度范围为地面~40 km，数据时间范围与微波辐射计资料一致。

1.1.3L波段探空资料

探空资料为重庆沙坪坝站2017年1月—2019年12月每日08和20时（北京时，下同）L波段探空秒数据。数据比对前，对气温和相对湿度数据进行异常值检查，剔除缺失和明显不合理的数据。

1.2 方法

为了实现Metop-A掩星资料和探空资料与微波辐射计资料的对比，将掩星及探空气温、湿度廓线用三次样条插值方法插值到与微波辐射计相同的各高度层上。由于微波辐射计站点固定，掩星数据经纬度不固定，为便于分析，选择时间差＜1.5 h、距离差＜150 km范围内的数据与微波辐射计进行匹配[29]（为保证匹配数据量，两种数据在0.6 km以下未纳入统计分析），计算Metop-A、探空与微波辐射计气温与相对湿度数据之间的相关系数、偏差、偏差标准差。公式如下：

2 偏差分析

2.1 微波辐射计与掩星

2.1.1总体对比分析

由（图1）可知，两种仪器探测的气温和相对湿度的相关系数均在0.8以上，尤其是相对湿度，整层的相关系数在0.9以上（图1a）。而气温在2 km以下相关性最好，系数均在0.9之上，4～5 km则相对较低（图1b）。

从偏差来看，两种探测设备的气温平均偏差为1.08℃。4 km以下，以负偏差为主（图1c），即微波辐射计探测的气温高于掩星。其中，1 km以下的差异最小，基本在1℃以内，2.5 km处差异最大，为3℃。4 km以上，微波辐射计气温则相对偏小。在微波辐射计的整个探测高度，掩星与微波辐射计的偏差分布均较为稳定，整层波动较小。从图1d可知，微波辐射计和掩星相对湿度的差值平均值大致在±5%以内。5 km以下，以负偏差为主，微波辐射计反演的相对湿度偏大，0.5 km处负偏差最大，为-4.5%，这是因为微波辐射计在高湿的环境下，其水汽廓线的反演精度问题所致。5 km以上，二者的偏差较小，基本在±1%以内，且微波辐射计和掩星相对湿度的这种偏差分布，整层均表现得较稳定，只是低层波动稍高。

微波辐射计反演的湿度廓线的质量稍高于气温。在水汽相对集中的中低层，微波辐射计反演温、湿度较掩星偏大。但其探测的温、湿度偏差却相对较小，尤其是气温，这是由于微波辐射计在低层的廓线信息来自于直接探测，而相对高层的廓线信息是来自于网络模型的反演所得，从而使得低空的测量精度要高于高空。

2.1.2冬、夏半年对比分析

图2为微波辐射计和掩星探测的气温、相对湿度在冬、夏半年的各统计量随高度的变化特征（冬半年指10月—次年3月，夏半年指4—9月）。从图2a可以看出，微波辐射计和掩星反演的气温在低层相关性相对较高，且冬半年的相关性要稍高于夏半年。5 km以下，冬、夏半年的气温相关系数随高度呈递减趋势，5 km以上则相反。在整个探测高度，冬、夏半年相对湿度的相关性随高度的变化与气温稍有不同（图2b），1 km以下和5 km以上，夏半年稍高于冬半年，1～5 km，两个半年的湿度廓线的相关系数随高度变化的趋势基本相同。

由图2c可知，微波辐射计气温偏差分布与图1c相似，即低层以负偏差为主，且2～3 km为负偏差的大值区，微波辐射计气温偏高3～4℃，而高层微波辐射计气温则相对偏低。就整层来看，夏半年偏差小于冬半年，即微波辐射计探测的温度廓线的精度，在夏半年优于冬半年。就两种设备探测的湿度廓线来看，其偏差随高度的变化趋势虽然与气温不同，但也呈现出夏半年优于冬半年的特点，尤其在低层（图2d），冬半年的偏差最大时达到-7%，而夏半年基本在±3%以内，但在低层，两种观测设备测量的相对湿度差值的稳定性，冬半年比夏半年好（图2f）。

图1 掩星资料与微波辐射计相关系数、差值平均值、差值标准差

图2 掩星资料与微波辐射计在冬、夏半年相关系数、差值平均值、差值标准差

就两种设备探测的气温来看，冬半年的相关性高于夏半年，但夏半年的差值小于冬半年，即微波辐射计对气温的反演精度在夏半年优于冬半年。就相对湿度而言，不管是从相关性，还是偏差分布，均呈现出夏半年优于冬半年的特性。经统计，重庆主城区夏半年平均总云量约7成，低云量约5成，而冬半年总云量有8.5成，低云量有7成，冬半年的云量明显多于夏半年，而云量的多少也会影响微波辐射计的探测精度[27]。

2.1.3不同天气条件下对比

图3为无降水和有降水时微波辐射计与掩星气温、相对湿度廓线的相关系数、差值平均值及差值标准差分布。由图3可知，有降水时，微波辐射计和掩星的相关性与无降水时的相关性差别较大。4 km以下（图3a），相关系数相差不大，基本在0.9左右，即有、无降水的天气条件下，两种设备反演的气温值在低层相关性均较好。4 km以上时，降水天气对微波辐射计探测的气温影响较大，相关系数降到了0.7左右。对相对湿度而言，降水天气时，微波辐射计与掩星的相对湿度相关性在中、低层均较差（图3b）。6 km以下，系数在0.6～0.9，且变化波动较大，而高层的相关系数与无降水时基本相同。

从差值平均值图（图3c）可知，有降水时，气温差值明显比无降水时大，尤其在2.5 km高度上，气温偏差幅度达到7℃；3 km以上，二者的差值有逐渐减少的趋势，在整个探测高度上，降水天气时，微波辐射计反演气温均比掩星偏大。在中低层，掩星与微波辐射计反演的气温差值在有、无降水天气下均有较稳定的表现（图3e），但在4 km以上，降水条件下两种探测设备反演的气温差值稳定性明显要比无降水时的差。从相对湿度差值图（图3d）可知，降水天气对微波辐射计反演的湿度廓线有一定影响，在1 km以下及4～6 km，降水天气时的负偏差值明显比无降水时大。

图3 掩星资料与微波辐射计在晴、雨天气状况下相关系数、差值平均值、差值标准差

从以上分析可知，降水会对微波辐射计的探测精度产生影响，因为降水时，水汽会附着在微波辐射计的天线罩上，而天线罩上的液态水会导致亮温观测值偏高[14]，从而造成雨天的探测误差偏大。

2.2 微波辐射计与探空

2.2.1总体对比分析

由于探空释放的时间一般为08和20时，且一次过程大概持续正点前后的1 h，因此微波辐射计资料选取探空采样时间段内探测数据的平均值，代表08和20时的值，与对应的探空数据进行对比分析。

从微波辐射计与L波段探空资料的温、湿度廓线的统计分析可以发现，两种设备探测的气温、相对湿度相关性和偏差分布特征与图1相似。在整个探测高度上两种要素均为显著正相关，且相对湿度相关系数稍大于气温，整层均在0.9以上；而气温在低层相关性优于中高层，尤其是2 km以下，相关系数均在0.9以上，之后随着高度的上升，相关系数减少，5 km处最小，为0.8，随后又逐渐增大；在中低层，微波辐射计探测气温、相对湿度均高于L波段探空，高层则相反，且低层偏差幅度小于高层，即微波辐射计在低层探测的温、湿度廓线信息精度优于高层。这跟微波辐射计与掩星的偏差分布基本相同。

2.2.2不同时次对比分析

图4为L波段探空和掩星探测的气温、相对湿度在08、20时两个时次各统计量随高度的变化特征。4.5 km以下，两个时次气温相关性相似（图4a），即随高度上升相关系数呈现递减趋势，在低空的相关性较高，尤其是20时在1 km处相关系数为0.97，4 km处减为0.83。4.5 km以上，两个时次相关系数均随高度递增，但20时相对高于08时，相对湿度的相关性稍高于气温（图4b），尤其是2~5 km，呈递增趋势，5 km处最大，两个时次相关系数均达到0.99。

08时微波辐射计反演气温与探空气温的偏差较20时小（图4c），但两个时次的差值随高度的分布形式相似，即4 km以下，均以负偏差为主，而4 km以上08时的偏差普遍小于20时，但20时的偏差在整个高度层的稳定性要优于08时，尤其是在高层。对相对湿度而言，中低层以负偏差为主，即微波辐射计反演相对湿度偏大，而5.5 km以上，二者的偏差较小，大部分在1%以内。在整个探测高度上，L波段探空与微波辐射计相对湿度的这种偏差分布，均表现得较为稳定（图4f）。因此，无论是从相关性还是偏差分布来看，20时微波辐射计相对湿度值要优于08时。

图4 探空资料与微波辐射计在08、20时相关系数、差值平均值、差值标准差

3 个例分析

3.1 大雾过程

2017年12月24日有一次较强的大雾过程。24日02时以后沙坪坝能见度逐渐降低，尤其是07、08时，能见度只有58、72 m（图5a）。24日01时开始近地面相对湿度达到90%以上，且高湿区向150 m高度扩展（图5b），13时左右，低层的高湿区开始逐渐消失减退。从微波辐射计的温度剖面图（图5c）可知，24日02时开始，近地层的气温逐渐降低，0.15～1 km有一“暖盖”，且温度低值区范围有向高层延伸趋势，07时左右，“暖盖”消失，而后随着能见度的逐渐增大，低层气温开始回升。

图5 2017年12月23日21时—24日16时沙坪坝小时能见度（a，单位：m），微波辐射计相对湿度垂直剖面（b，单位：%），微波辐射计温度垂直剖面（c，单位：℃）

此次大雾正是冬季常见的辐射雾，夜间由于地面辐射冷却作用，低空水汽凝结，尤其是50 m高度以下湿度增大，而微波辐射计相对湿度的增大与实际观测能见度下降几乎同时发生。地面辐射导致近地层气温降低，随着太阳辐射的增强，辐射雾逐渐消失，虽然微波辐射计探测低空相对湿度高值区的消失时间比能见度增大时间稍晚，但其探测到近地面气温回升的同时，能见度开始增加。

3.2 降水过程

2019年4月19日—20日，重庆出现了局地短时强降水，雨量主要集中在重庆中西部及中偏北地区。从微波辐射计观测来看，19日14时开始5 km以下相对湿度达到85%以上（图6b），15、16时左右沙坪坝站上空1～5 km探测高度上相对湿度几乎达到100%，而此时对应沙坪坝站小时最强降水（图6a），即16时小时降水量达到29.9 mm，之后高湿区抬升，而中低层相对湿度降低，降水逐渐减弱，20时第一阶段降水结束。20日02时开始，1～8 km出现85%以上的高湿区，此时台站第二阶段降水开始，且在03时达到最大（22.7 mm）。

从温度剖面图发现（图6c），19日13时，温度高值区抬升，15时左右，2.5～4.5 km的高度上出现了较为明显的暖区，且暖区持续到18时左右减弱消失，随后1～2 km为温度低值区，上下形成逆温层。20日01时，3 km左右的暖区再次出现，直至降水结束。3 km高度层暖区出现的时段正好对应两次降水时段。

图6 2019年4月19日11时—20日08时小时降水量（a，单位：mm）和微波辐射计相对湿度（b，单位：%），垂直剖面温度（c，单位：℃）

为进一步检验此次降水过程中的温度特征，选用与微波辐射计同址的沙坪坝站L波段探空数据对此次过程进行分析（（沙坪坝探空站海拔高度是541 m，因此500 m以下不予分析）。图7为4月19日20时和20日08时L波段探空温度垂直剖面图，由图可知，19日20时，在1 km左右的高度上，有一明显的逆温层存在，而20日08时这种逆温现象基本消失，这与微波辐射计探测的逆温现象出现时间相同，只是两种设备探测的逆温层高度稍有不同。

图7 2019年4月19日20时和20日08时沙坪坝L波段探空温度垂直廓线（单位：℃）

综上可知，此次短时强降水过程中，5 km以下出现高湿区时，降水开始，之后随着暖湿气流的上升，低层水汽减少，降水减弱，而水汽到达高层后容易达到饱和而释放凝结潜热[28]，使得2～3 km高度层温度升高，随着低层气团的上升，水汽绝热冷却导致低层降温，在两次降水间歇，低层出现较为明显的逆温层。虽然降水会对微波辐射计测量结果产生影响，但此次降水过程中的温、湿度特征已较好地体现出来。

4 结论与讨论

通过对2017年1月—2019年12月沙坪坝微波辐射计与Metop-A掩星、L波段探空气温、相对湿度两个要素进行对比并分析微波辐射计资料在不同天气条件下的探测特征，得出如下结论：

（1）总体来看，微波辐射计与掩星反演所得气温和相对湿度在整个探测高度上均为显著正相关，低层相关性高于高层，且相对湿度稍高于气温。在中低层，微波辐射计反演温、湿度较掩星偏大。相比于高层，微波辐射计在近地面层探测的温、湿度偏差相对较小。两种设备探测的气温，冬半年的相关性稍高于夏半年，但夏半年偏差小于冬半年。就相对湿度而言，不管是从相关性，还是偏差分布，均呈现出夏半年优于冬半年的特性。

（2）降雨天气对微波辐射计探测精度影响较大。4 km以上，降水天气时气温相关性较低，相关系数仅有0.7。在水汽相对集中的低空区域，微波辐射计反演气温高于掩星，尤其是降雨时，这种偏差更大。1 km以下及4～6 km，降水天气时微波辐射计相对湿度的负偏差值明显比无降水时大。

（3）微波辐射计与L波段探空气温、相对湿度在整个探测高度上也为显著正相关，表现出低层相关性高于高层的特性，且4 km以下微波辐射计探测气温、相对湿度值高于L波段探空值；08时气温平均偏差小于20时，而20时相对湿度探测精度优于08时。

（4）不同的天气过程，微波辐射计相对湿度和气温的变化与天气现象有较好的对应关系。大雾天气时，微波辐射计探测低层相对湿度的增大与实际观测能见度下降几乎同时发生，且地面辐射导致近地层气温降低时段微波辐射计也能较准确地探测出来。通过对2019年4月19—20日短时强降水过程分析发现，降水时段，微波辐射计探测5 km以下为高湿区，随着暖湿气流的上升，降水减弱，暖湿气团上升过程中高层的凝结潜热和低层的绝热冷却作用，致使中低层出现了逆温层。

上述研究表明，微波辐射计反演的温湿度廓线可信度较高，不同的高度、季节和天气条件下，微波辐射计气温和相对湿度廓线均有不同的表现。晴朗少云的天气，其探测精度相对较高，而降水时，水汽对微波辐射计的探测精度有一定的影响。不同的天气条件下，微波辐射计的探测特征不同，后续将考虑加入更多的对比数据源对其探测的各类天气进行系统分析，以便为气象预测预报提供参考依据。