基于“上升-平漂-下降”探空资料的长江中下游暴雨同化试验*

张旭鹏 郭启云 杨荣康 马旭林 曹晓钟

1 南京信息工程大学气象灾害教育部重点实验室，南京 210044

2 中国气象局气象探测中心，北京 100081

提 要： 为进一步讨论新型“上升-平漂-下降”探空数据在资料同化与数值预报中的应用效果，基于WRF(Weather Research and Forecast)模式及WRFDA(WRF data assimilation)同化系统进行同化对比试验。在对新型探空试验数据进行质量评估和稀疏化的基础上，将下降段资料与常规观测资料组合同化，并讨论其对于长江中下游地区暴雨预报质量的影响及原因。主要试验结果包括：通过与FNL资料、业务同站探空数据交叉对比验证最新试验数据准确性；使用特性层与规定层结合的方案对新型探空上升、下降段进行稀疏化处理可以得到较优效果；同化下降段数据能够在一定程度上提高暴雨预报技巧；风场及湿度场的调整是暴雨预报技巧有所提高的重要原因之一。

引 言

高空探空观测在当代气象观测中具有极其重要的地位，其结果常用作描述大气状态的相对真值。虽然探空站点数量较为有限，但空间分布相对均匀且其数据中包含丰富的大气垂直信息，能够反映大气三维结构，是高空气象观测中最主要的直接观测方式。探空数据不仅可以作为地基微波辐射计(Xu et al,2015)、毫米波云雷达(赵静等，2017)、风廓线雷达(吴蕾等，2014)、COSMIC掩星(杜明斌等，2009；郭启云等，2020)、及卫星反演(Kwon et al，2012；陈源等，2020)等遥感产品或再分析资料(张思齐等，2018)的检验评估标准，在中短期天气预报中也有极其重要的作用(Laroche and Sarrazin，2013)，是减小预报误差的重要观测手段之一(Singh et al，2014)。李佳英等(2006)通过对比探空资料与模式温湿廓线，检验模式对强对流天气的预报能力。当前业务使用的L波段秒级探空数据可深入揭示大气内部运动规律，例如应用于高层大气温度结构(李刚等，2014)或大气重力波(吴泓锟等，2019)研究中等等。

但固定的常规探空难以满足对移动气象目标跟踪观测的需求。下投探空是指由气球、无人机(肖良华等，2019)、火箭(李金辉等，2020)等飞行器携带配有降落伞的探空仪升空，并在一定高度进行投放，从而获得由下投平台至地(洋)面的大气廓线(陈洪滨和朱彦良，2008)。下投探空可以作为常规探空的补充，弥补特定时刻一些区域内探空资料的空白。目前，下投探空观测已经成为观测台风的主要手段之一，对台风的移动路径及强度预报具有极其重要的意义(曲晓波和Heming，2006；张诚忠等，2012；李杨等，2016)。

探空资料在同化领域有着极为广泛的应用。张利红等(2013)指出同化探空和云迹风资料可以提高我国西南地区降水预报准确率。Yamazaki et al(2015)通过OSE(observing system experiments)验证探空资料可以改善北极区域对流层上层环流模拟情况。莫毅等(2008)、宝兴华和杨舒楠(2015)和孟智勇等(2019)分别使用WRF-3DVar、WRF-EnKF和ETKF方法对探空资料进行同化，证明同化探空资料能够改善暴雨落区及降水强度模拟效果。Choi et al(2015)使用WRF-4DVar同化探空资料时加入气球漂移信息，结果表明降水落区预报效果及QPF(quantitative precipitation forecast)技巧都有所提高。同化探空资料对台风路径及其降水预报也有明显改进作用(董海萍等，2017；龚俊强等，2019)。另外，有多个研究发现：使用加密探空数据可以显著提高数值预报技巧(陈朝平等，2012；徐同等，2016；高笃鸣等，2018；王丹等，2019)，且其正效应能够传播至下游地区(Faccani et al，2009)。L波段秒级探空数据通过插值或稀疏化处理后引入数值模式能够很好地改善模拟效果：郝民等(2014)通过插值到最近模式层的方式将秒级探空数据引入GRAPES模式后，模拟高层系统更为准确且降水预报技巧有所提高；姚爽等(2015)将秒级探空数据从时间上稀疏化为分钟数据并将其加入WRF模式，可以改善高空风场模拟效果及大量级降水预报性能。

为拓展现有探空能力和范围，解决探空加密的关键性技术问题，中国气象局气象探测中心在现有同球双释放技术(郭启云等，2018)和北斗卫星导航定位测风技术基础上研发了新一代探空观测系统，即“上升(1 h)-平漂(4 h)-下降(1 h)”三段式探空观测。该系统在不增加人员投入和经济成本的条件下，实现了探空观测的时空加密，能够显著提升探空观测效益(钱媛，2019；曹晓钟等，2019)。上升段观测与业务探空效果相近；而下降段则可近似看作下投探空，其中包含和上升部分相同的观测要素，且下投地点实时可控，具有很强的机动性，对业务探空有很好的补充作用。当前我国长江中下游地区已建设由六个新型探空试验站点构成的观测网，建立了观测试验数据集，并对其进行了质量控制与偏差订正(包括温度辐射订正、湿度滞后订正、气压平滑等)。曹晓钟等(2019)、钱媛(2019)和王丹等(2020)对试验数据进行质量评估，其结果初步验证了“上升-平漂-下降”探空的可信可用性。杨晨义等(2021)则基于新型探空资料对下平流层重力波特征进行分析。

“上升-平漂-下降”探空作为一种新型探空观测方式，在原有上升基础上加入平漂和下降段观测能力。其中，下降段资料在同化中的应用效果还需要进一步的研究与验证，目前尚且缺乏这一方面的研究。本文将基于“上升-平漂-下降”探空试验数据进行同化对比试验，并根据试验结果讨论新型探空系统对长江中下游地区暴雨预报的改进潜力。

1 试验数据质量评估

同化新型探空数据之前，需先对其质量进行评估。前人在此方面已有一些研究结果，使人们对“上升-平漂-下降”探空数据质量有了初步了解。但随着试验数据集及质量控制方法不断发展完善，该项工作仍有进一步进行的必要。目前对“上升-平漂-下降”探空数据质量的评估方法大致分为两种：一种是利用现有同站业务探空资料对新型探空数据质量进行评估；另一种则是基于高分辨率模式分析场或再分析资料进行质量评估。考虑到业务探空资料的时空分辨率问题，第一种方法可用于上升、下降段评估，但对平漂段评估能力较弱。而第二种方法则可以相对更合理地评估全部三段探空资料，其缺点是由于不同模式资料在分辨率和质量上存在差异，故参考不同资料得到的评估结果可能有所差异。本文在前人研究基础上，选择时间分辨率为6 h、空间分辨率为1°×1°的FNL分析资料以及同站秒级探空资料，对2019年最新试验数据集质量进行评估。

如图1中所示，上升段和下降段温度在日间和夜间都与FNL资料具有较好的一致性，偏差基本在±2℃。日间平漂段温度经过初步辐射订正后在平漂开始和结束阶段仍存有较明显的正偏差，达到+12℃；夜间没有太阳辐射干扰的情况下，平漂段温度质量相对较好，偏差基本与上升段和下降段相当，能够达到±2℃。目前，气象探测中心已初步设计出测量日间平漂段温度“真值”的方案，相信未来日间平漂温度辐射偏差能够得到有效的订正处理。

图1 基于FNL资料的武汉站新型探空温度评估2019年5月12日(a,b)日间和(c，d)夜间(a，c)温度一致性以及(b，d)温度偏差Fig.1 New sounding temperature evaluation of Wuhan Station based on FNL data： (a, c) temperature consistency and (b, d) temperature deviation during (a, b) daytime and (c, d) nighttime on 12 May 2019

基于2019年最新试验数据，选取同站业务秒级探空及FNL资料与之进行交叉对比，结果如图2。上升段湿度廓线与FNL资料一致性较差，但偏差在整体上小于同站业务探空，证明新型探空仪具有更强的湿度探测能力，下降段湿度偏差整体上大于上升段。两类探空数据相对分析资料整体上较为偏干，一般在10 km附近湿度偏差最为明显。上升和下降段测得的气压廓线与FNL资料、同站探空的一致性较高，气压偏差整体在±4 hPa。10 km以下新型探空资料的气压偏差略小于业务探空数据偏差。

图2 基于FNL资料和业务探空资料对2019年5月12日00 UTC武汉站(a～d)上升段和(e～h)下降段的数据评估：(a,e)湿度一致性与(b,f)湿度偏差;(c,g)气压一致性与(d,h)气压偏差Fig.2 Data evaluation of (a-d) ascending section and (e-h) descending section of the sounding at Wuhan Station at 00 UTC 12 May 2019 based on FNL data and operational sounding data： (a, e) humidity consistency and (b, f) humidity deviation， (c, g) pressure consistency and (d, h) pressure deviation

2 “上升-平漂-下降”探空数据稀疏化方案

新型探空观测系统中秒级数据分辨率远超当前模式分辨率，产生的冗余效应会在同化过程中造成负效果(郭欢，2017；马旭林等，2019)。故在同化试验前需先进行稀疏化处理。本文将规定层数据与特性层数据相结合的稀疏化方案(李庆雷等，2018)应用于新型探空系统上升段和下降段中。为避免试验数据集气压与规定层气压在数值上不完全一致而造成规定层遗漏，试验选取规定层一定小阈值内的探空数据。文中设置的上升段与下降段温度特性层阈值为2 ℃。稀疏化对比效果如图3所示，使用稀疏化方案条件下得到的上升、下降段温度廓线，相比只提取规定层邻近数据的温度廓线更能够反映原始秒级数据中温度廓线的变化特征。

图3 两种稀疏化方法得到的(a)上升段、(b)下降段温度廓线对比Fig.3 Comparison of temperature profiles in (a) ascending section, (b) descending section obtained by two thinning methods

由于平漂段中气象要素变化一般较为平缓，规定层结合特性层的稀疏化方法尚不适用于平漂段数据。当前平漂段数据多采用从时间维度上进行稀疏的方法。

3 暴雨个例

2018年7月4—7日，受低空切变线和低空急流影响，湖北中东部、湖南南部大部出现大暴雨，江西局部出现特大暴雨。天气图显示(图略)，6日00 UTC500 hPa上中高纬度呈现“两槽一脊”的形式，东侧槽南伸至长江中下游地区，利于引导冷空气南下；850 hPa上长江中下游地区受低压带控制，存在低空风切变，有利于暴雨的产生。其中，6日18 UTC至7日18 UTC的24 h降水实况如图4。雨带整体呈东北—西南向走势，主要包括三个强降水中心，分别位于江西南部地区、湖南南部和广西北部地区。强降水中心24 h累计雨量能够达到100 mm以上。

图4 2018年7月6日18 UTC至7日18 UTC 24 h降水实况Fig.4 The observed precipitation in 24 hours from 18 UTC 6 to 18 UTC 7 July 2018

4 资料与试验方案

本文同化的观测资料包括业务常规观测资料(Global Telecommunication System,GTS)和长江中下游六个站点(安庆、南昌、长沙、赣州、宜昌和武汉)的新型探空试验资料。常规观测资料中包括：自动站资料、业务探空资料、云迹风资料和飞机报资料。降水实况资料采用中国气象数据网(https:∥data.cma.cn/)提供的中国自动站与CMORPH降水产品融合的逐时降水网格数据集(1.0版)。分析过程中使用欧洲中期预报中心(ECMWF)提供的分辨率为0.25°×0.25°的再分析数据集ERA5作为参考。

试验采用WRF V3.9模式与WRFDA V3.9同化系统，模拟区域中心为(28°N、113°E)。模式采用两层嵌套，水平分辨率分别为9 km和3 km，垂直方向38层，网格设置如图5。使用FNL资料提供初始场和边界条件。物理参数化方案分别为：微物理参数化方案WSM5，Grell三维集成积云方案，长波辐射方案RRTM，短波辐射方案Dudhia，行星边界层方案YSU，近地面层方案MM5。

图5 WRF模拟区域设计Fig.5 Design of the WRF simulation area

本文设计控制试验CTL和敏感试验Down，以研究原有常规资料基础上加入下降段资料后对暴雨预报的改善效果，具体试验方案如表1。同化方法采用WRF 3D-Var(three-dimensional variational)。由于新型探空气球释放前后经纬度差异较大，本文使用逐点同化的方式来消除气球位置漂移造成的影响。同化要素包括温度、气压、湿度、位势高度、风速风向。Down试验中同化的观测资料分布情况如图6 所示。

图6 Down试验模拟区域内的资料分布(蓝色圆点：自动站资料，绿色圆点：云迹风资料，橙色圆点：飞机报资料，红色圆点：业务探空资料，☆：下降段资料)Fig.6 Data distribution map of down test in the simulation area(blue dots： automatic station data, green dots： cloud trace wind data, orange dots： aircraft report data, red dots： operational sounding data，☆： descending section data)

表1 试验方案Table 1 Test scheme

5 试验结果

5.1 24 h降水预报

控制试验与敏感试验模拟的24 h累计降水分布如图7。可见CTL与Down两组试验对实况中的东北—西南向雨带都具有一定模拟能力，都模拟出了江西、湖南南部的强降水区域。但在控制试验中，江西南部区域出现了较为明显的200 mm以上虚假强降水中心，且湖南西南部也存在较大范围的降水误报，与实况不符。在敏感试验中，江西南部、湖南西南部的虚假强降水区域在范围和强度上都有一定程度减小，降水落区分布相对更加接近于实况。

图7 试验预报的2018年7月6日18 UTC至7日18 UTC 24 h累计降水(a)CTL试验，(b)Down试验Fig.7 The simulated 24 h accumulated precipitation from 18 UTC 6 to 18 UTC 7 July 2018(a) CTL test， (b) Down test

5.2 ETS降水评分

ETS(equitable threat score)客观检验评分相比TS(threat score)评分，增加了对降水空报、漏报进行惩罚的部分，使得评分相对更加公平，其具体计算公式如下：

(1)

(2)

式中：NA为预报正确的格点(次)数，NB为空报格点(次)数，NC为漏报格点(次)数，ND为预报与实况均未达到指定阈值的正确格点(次)数。根据24 h 降水量，划分为小雨(0.1～9.99 mm)，中雨(10～24.99 mm)，大雨(25～49.99 mm)，暴雨(50～99.99 mm)和大暴雨(>100 mm)，共五个量级。

CTL和Down试验在主要强降水区域(24°～28°N、108°～118°E)内的ETS评分如图8所示。可以看出，除了中雨量级上敏感试验预报效果略弱于控制试验外，其他降水量级上敏感试验的预报技巧上相比控制试验中皆有所提高，尤其对小雨、大雨和大暴雨的预报效果改进较为明显。Down试验中大暴雨ETS评分由CTL中的0.006增长至0.012，大雨和暴雨的ETS评分增幅达到29.0%和9.2%。

图8 24 h降水预报ETS评分Fig.8 The 24 h precipitation forecast ETS

总体来看，加入下降段数据后有利于降水预报技巧的提高，尤其对于大雨量级以上的强降水预报技巧改进较为明显。

5.3 分析增量分布

分析增量虽不能直接反映同化效果的好坏，但可以体现加入观测资料后对背景场的调整情况。以下给出同化下降段资料前后850、500、200 hPa上高度场、温度场、湿度场分析增量分布特征。

5.3.1 高度场分析增量

控制试验与敏感试验高度场分析增量如图9。850 hPa低层，两组试验正负增量中心都达到2 m和-5 m，差异较小，敏感试验中负增量略大；500 hPa中层，敏感试验在湖北西北部的-5 m负增量中心相比控制试验中有所增强，且湖南广西一带-2 m 负增量区范围有所增大；200 hPa高层，敏感试验中湖北西部负增量中心范围增大，且在广西北部增加了-5 m负增量中心。从垂直方向上来看，同化下降段资料对高层高度场调整较为明显，对低层高度场调整相对较小。

图9 (a，b，c)CTL试验和(d，e，f)Down试验在(a,d)850 hPa，(b,e)500 hPa，(c,f)200 hPa的高度场分析增量Fig.9 The height field analysis increment of (a， b， c) CTL test and (d， e， f) Down test at (a, d) 850 hPa (b, e) 500 hPa, (c, f) 200 hPa

5.3.2 温度场分析增量

温度分析增量如图10所示。850 hPa上，控制试验中，湖北、湖南、广西一带0.1℃正增量区域在敏感试验中明显减小；500 hPa，敏感试验中，湖南南部-0.5℃负增量中心较控制试验中略有增强；200 hPa控制试验和敏感试验都以正温度分析增量为主，敏感试验中0.3℃正增量范围有所增大。在本次个例中，同化下降段资料对低层和高层温度场改变较为明显，中层温度场变化相对较小。

图10 同图9，但为温度场分析增量Fig.10 Same as Fig.9, but for the temperature field analysis increment

5.3.3 相对湿度场分析增量

相对湿度分析增量分布情况如图11。850 hPa上，控制试验在湖北、湖南、广西一带的负增量在敏感试验中有所减小；500 hPa上，敏感试验中安徽南部、江西北部地区的正增量较控制试验增加了1%～2%；200 hPa上，控制试验和敏感试验都以负相对湿度分析增量为主，二者差异较小。从垂直角度来看，同化下降段资料对中低层相对湿度调整相对明显，对200 hPa高层湿度场调整不明显，这可能与传感器在高层(约12 km)测量湿度质量相对较差有关。

图11 同图9，但为湿度场分析增量Fig.11 Same as Fig.9, but for the humidity field analysis increment

5.4 物理量场分析

2018年7月6日暴雨过程的主要影响系统包括低空急流和与低空切变线。如图12所示，ERA5资料、CTL与Down试验都模拟出湖北、湖南地区上空的低空急流。虽然敏感试验中只加入了三段下降段资料，但由于同化时刻研究区域内原有探空资料偏少，且下降段资料的位置位于关键区域附近，故Down试验中预报的低空急流较CTL试验中有较明显的区别。CTL试验预报的低空急流分为南北两支。在Down试验中，低空急流虽南部风速较CTL试验中有所偏强，但北部急流在空间范围和风速上较CTL试验中都有明显减小，使其在空间范围上整体更接近于ERA5再分析结果。

图12 2018年7月6日18 UTC 850 hPa低空急流和风速(填色)、风向(矢量箭头)(a)ERA5资料,(b)CTL试验,(c)Down试验Fig.12 The 850 hPa low-level jet, wind speed (colored) and wind direction (vector arrow) at 18 UTC 6 July 2018(a) ERA5 data， (b) CTL test， (c) Down test

水汽通量可以反映降水过程中水汽的来源和输送情况。CTL与Down试验对水汽通量的模拟效果如图13。可以看出本次降水过程的水汽输送情况很大程度上受到江淮气旋的影响，气旋携带南部洋面的水汽，再将其由北向南输送。Down试验中预报的江西地区水汽通量与CTL试验中相差较小，Down试验中略强。但Down试验在安徽、河南、湖北一带向南输送的水汽通量以及广西地区向北输送的水汽通量都较CTL试验有所减弱，削弱了强降水所需的水汽持续供应，从而有利于减弱湖南西南部的虚假强降水。

图13 2018年7月6日18 UTC 850 hPa水汽通量(填色，单位： kg·m-1·s-1·hPa-1)和风向、风速(矢量箭头)(a)CTL试验，(b)Down试验，(c)Down试验与CTL试验的水汽通量之差Fig.13 The 850 hPa water vapor flux (colored， unit: kg·m-1·s-1·hPa-1), wind direction and wind speed (vector arrow) at 18 UTC on 6 July 2018(a) CTL test, (b) Down test, (c) difference of water vapor flux between Down test and CTL test

通过比较试验预报的整层大气可降水量，可以进一步验证以上结论。如图14所示，控制试验与敏感试验预报的长江中下游地区大气可降水量分布情况整体较为相似，但Down试验预报的湖南西南部大气可降水量较CTL试验预报结果有明显减小，其中负值中心达到-3.5 mm，且湖北南部及湖南北部区域的大气可降水量也有一定程度的减小，与水汽通量中的表现相一致，有利于削弱该地区的虚假强降水。

图14 2018年7月6日20 UTC大气可降水量(a)CTL试验，(b)Down试验，(c)CTL-DownFig.14 Atmospheric precipitable water at 20 UTC 6 July 2018(a) CTL test, (b) Down test, (c) CTL-Down

控制试验与敏感试验预报的对流有效位能(CAPE)也具有一定差异。如图15a、15b所示，两组试验在同化时刻的CAPE值总体分布较为相近，但由图15c可以看出，Down试验预报的湖南西南部以及江西南部地区CAPE值小于CTL试验中的预报结果，一定程度抑制了CTL试验中虚假强降水的发生发展，有利于强降水预报技巧的提高。

图15 同图14,但为7月6日18 UTC对流有效位能Fig.15 Same as Fig.14, but for CAPE at 18 UTC 6 July 2018

6 结论与讨论

本文基于长江中下游新型“上升-平漂-下降”探空观测数据集设计了同化对比试验。在对新型探空数据进行质量评估和稀疏化处理的基础上，将下降段资料与常规观测资料结合并进行同化对比试验，根据试验结果讨论新型探空系统在长江中下游暴雨预报中的应用效果及前景。主要结论有：

(1)通过将“上升-平漂-下降”探空试验数据与FNL资料以及同站业务秒级探空数据进行交叉对比，检验了最新试验探空数据质量。

(2)使用规定层结合特性层的稀疏化方案能够对“上升-平漂-下降”探空系统中上升段、下降段数据进行有效的稀疏化处理。

(3)本文试验中，同化新型探空下降段数据后，能够改善控制试验中出现的降水高估、误报等问题，提高强降水预报技巧。其中大暴雨ETS评分由控制试验中的0.006增长至0.012，大雨和暴雨的ETS评分增幅分别达到29.0%和9.2%。

(4)同化下降段观测资料能够对多个高度上的高度场、温度场和相对湿度场产生一定的调整作用。

(5)在原有常规观测资料基础上加入下降资料并同化后，风场出现了较为明显的调整，进而改善了水汽输送情况，这可能是导致降水预报技巧提高的重要原因之一。

本文试验结果证明“上升-平漂-下降”探空系统在同化方面有很好的应用前景。但由于目前新型探空系统尚未业务化运行，可用于同化试验的数据有限，所得结论尚不具广泛代表性，随着日后数据集的完善，可做更多试验验证当前结论。未来工作中可尝试由高低层风场相互作用及能量输送的角度入手，将更多数据引入到天气、气候模式中，从而更好地体现新型探空系统在资料同化与数值预报领域的价值。