2021年夏季CMA-MESO模式降水预报评估_蔡怡.pdf

第 41 卷 第 3 期2023 年 6 月Vol.41 No.3June，2023干旱气象Journal of Arid Meteorology2021年夏季CMA-MESO模式降水预报评估蔡怡1，2，徐枝芳2，3，龚玺4，钟若嵋5，黄观胜6，龙海川7（1.浙江省长兴县气象局，浙江 长兴 313100；2.中国气象局地球系统数值预报重点开放实验室，北京 100081；3.中国气象科学研究院灾害天气国家重点实验室，北京 100081；4.湖南省株洲市气象局，湖南 株洲 412000；5.中国民用航空中南地区空中交通管理局广西分局，广西 南宁 530048；6.湖南省郴州市气象局，湖南 郴州 423000；7.中国船舶集团有限公司，重庆 401120）摘要：利用中国2021年夏季（68月）逐3 h降水预报资料和地面气象站观测资料，从多个角度诊断分析CMA-MESO（China Meteorological Administration Mesoscale Model）3 km系统模式降水预报性能，为预报员提供参考，为模式系统改进提供依据。结果表明：CMA-MESO 3 km模式能较好地预报出不同地区3 h累积降水量和有效降水频率时空分布特征，区域降水预报能力强于单站，持续性降水预报效果好于局地短时强降水。不同预报时效结果显示，CMA-MESO 3 km模式3 h降水预报值最大且大于观测值，6、9、12 h降水预报最接近实况；短时强降水个例结果显示CMA-MESO 3 km模式短临预报效果较好，3、6 h降水预报接近观测，降水量及其时间变化特征预报与观测基本一致，且8个预报时效的区域降水预报平均值接近观测值，预报效果较好。关键词：CMA-MESO 3 km；降水量；降水频率；偏差文章编号：1006-7639（2023）03-0503-13 DOI：10.11755/j.issn.1006-7639（2023）-03-0503中图分类号：P456；P426.6 文献标志码：A引 言暴雨洪涝是中国最常见、最严重的自然灾害之一，约三分之二的国土面积都可能发生不同类型和不同程度的暴雨洪涝灾害，给人民群众生命安全、经济社会发展、国家安全等造成严重影响（李泽椿等，2015）。强对流天气具有突发性、局地性和极端性等特征，对流系统往往发展剧烈，易在短时间内造成极端灾害天气，如何对其准确预报是目前气象领域公认的难题（齐艳军等，2016；翟盘茂等，1999；Stensrud and Fritsch，1994；Stensrud et al.，1999；Stensrud et al.，2000）。近年来，随着数值模式的不断发展和高性能计算机能力的提升，高分辨率的区域快速更新循环同化系统发展迅速（Benjamin et al.，2004；薛纪善等，2008；陈子通等，2010；王叶红等，2011；傅娜等，2013）。2001年中国气象局开始自主研发新一代全球/区域通用数值天气预报系统GRAPES（Global/Regional Assimilation and Prediction Enhanced System），GRAPES区域数值预报业务系统（GRAPES-MESO，2021年更名为CMA-MESO）于2006年正式投入业务化运行（Ma et al.，2009），经过不断升级改进，要素场和降水预报准确性明显提升（黄丽萍等，2017；徐枝芳等，2013），2015年建成了覆盖中国东部区域3 km水平分辨率的试验系统。2020年覆盖全中国区域业务系统的分辨率从10 km升级为3 km（庄照荣等，2020），每日开展8次同化分析和预报，资料同化种类和资料量大幅增加，CMA-MESO模式在短临预报中发挥越来越重要的作用，极大地提高了降水天气过程的预报准确率，但对极端暴雨和强降水的捕捉能力仍有较大的提升空间。模式降水预报检验评估不仅可为预报员预报提供参考，还可为模式系统改进提供依据。有研究采用降水确定性预报评分如预报成功指数评分（Threat Score，TS）、公平成功指数评分（Equitable 蔡 怡，徐枝芳，龚 玺，等.2021年夏季CMA-MESO模式降水预报评估 J.干旱气象，2023，41（3）：503-515，CAI Yi，XU Zhifang，GONG Xi，et al.Evaluation of precipitation forecast of CMA-MESO model in summer of 2021 J.Journal of Arid Meteorology，2023，41（3）：503-515，DOI：10.11755/j.issn.1006-7639（2023）-03-0503收稿日期：2022-08-23；改回日期：2022-10-27基金项目：河南省重大科技专项项目（201400210800）和国家重点研发计划项目（2018YFC1507605）共同资助作者简介：蔡怡（1996），女，浙江杭州人，助理工程师，主要从事数值模式与资料同化研究。E-mail：。通信作者：徐枝芳（1977），女，北京人，研究员，主要从事数值模式与资料同化研究。E-mail：。41 卷干旱气象Threat Sore，ETS）和 预 报 偏 差 Bias，检 验 评 估 了CMA-GFS（China Meteorological Administration Global Forecast System）、ECMWF（European Centre for Medium-Range Weather Forecasts）、SMS-WARMS（Shanghai Meteorological Service-WRF ADAS Rapid Refresh System）等不同模式的降水预报性能，发现模式对小雨、中雨、大雨量级预报评分较高，而对暴雨预报评分总体较低（宫宇等，2018；梅钦等，2018；吴秋霞等，2007；张立凤和罗雨，2010）。另外，还有研究采用概率预报检验评分如Brier评分、相对作用特征（Relative Operative Characteristic，ROC）、可靠性 评 分（Reliability Score，RS）和 邻 域 空 间 检 验（Fraction Skill Score，FSS）对模式降水进行不确定性分析（孙建华等，2018；曹越等，2019；徐同等，2019），发现对于临近降水预报，概率预报效果较集合平均更具有优势，但对极端强降水的指示意义较差（徐致真等，2019；韦青等，2020；吴晶等，2020；王婧卓等，2021）。针对 CMA-MESO模式降水预报的定量统计检验已开展了大量研究：如史文茹等（2021）对“720”特大暴雨天气过程的模式预报性能进行检验评估，肯定了数值模式对此次特大暴雨具有一定预报能力；张小雯等（2020）利用降水确定性预报评分评估了CMA-MESO 3 km模式强对流天气预报性能，发现该模式对风暴等强对流天气的预报明显优于其他模式。另外，许晨璐等（2017）、谢漪云和王建捷（2021）对CMA-MESO模式预报降水量、降水频率、降水强度特征等综合表现也进行了检验评估，指出CMA-MESO 模式能够很好地捕捉日均降水量和降水频率与地域分布特征的关系，弱降水频率偏低，强降水频率与实况吻合，但强度偏弱。多项研究结果均表明CMA-MESO模式对中国区域强对流天气精细化预报具有明显优势（陈静等，2022；黄丽萍等，2017；刘雪晴等，2020）。但上述研究均基于某一定特定区域范围内（如青藏高原东部及四川西部的高原复杂地形、华南沿海地区等）24 h降水预报产品开展评估，对不同区域范围降水预报的对比分析较少，且对不同起报时间不同预报时效的模式预报评估分析尚未深入展开。因此，本文从中国2021年夏季（68月）每日逐3 h降水量和有效降水频率入手，通过划分不同区域及不同降水量级，从多视角定性和定量评估CMA-MESO 3 km模式8个起报时次预报结果，为CMA-MESO 3 km模式预报产品的深度应用和系统改进提供参考。1资料和方法1.1资 料所用资料为：2021年夏季（68月）10 796个国家级地面气象观测站逐3 h累积降水量资料，由于一些观测站降水资料存在缺失，对缺测率大于5%的观测站进行剔除以保证所用资料的完整性；同期CMA-MESO 3 km 模式每日 8 个起报时次（00：00、03：00、06：00、09：00、12：00、15：00、18：00、21：00）逐3 h地面降水预报产品。CMA-MESO 3 km模式的水平分辨率为 0.030.03，预报范围 70.00E145.00E、15.00N65.00N，每隔 3 h 同化一次资料，每日00：00和12：00冷启动一次，00：00和12：00启动的预报时效为 72 h，其他时次的预报时效为36 h，为确保诊断分析数据样本一致性，本文仅对前24 h预报产品做诊断分析。文中所有时间均为世界时。文中附图涉及地图均基于中华人民共和国民政部官网提供的审图号为GS（2022）1873号的标准地图制作，底图无修改。采用双线性插值法将降水预报产品插值到观测站点上，该方法对月降水量处理的结果比双三次插值法、最近邻点插值法和三次样条插值法有更高的精度（税军峰等，2019）。参考降水量等级划分的国家标准GB/T 28592-2012（国家气象中心，2012），将3 h累积降水量（RC）划分为7个等级，即：RC0.1 mm（无降水）、0.1RC3.0 mm（小雨）、3.0RC10.0 mm（中雨）、10.0RC20.0 mm（大雨）、20.0RC50.0 mm（暴雨）、50.0RC70.0 mm（大暴雨）、RC70.0 mm（特大暴雨）。1.2有效降水频率计算方法站点平均有效降水频率定义为：将每日00：0000：00站点观测及模式预报3 h累积降水量分为无降水（RC0.1 mm）、有效降水（RC0.1）和强降水（RC10.0 mm）3种情况进行统计。若RC0.1 mm，则有效降水次数Mn,k计为1（n=1，2，3，N；k=1，2，3，K；N、K 分别为总天数、总站数），若 RC0.1 mm，则Mn,k记为0，强降水次数统计同理。分别计算观测时段内每个站点夏季有效降水频率fk，再对所选区域内所有站点的有效降水频率累加平均，得到区域内站点平均有效降水频率-fsta，具体公式如下：fk=n=1NMn,kN（1）-fsta=k=1KfkK（2）504第 3 期蔡怡等：2021年夏季CMA-MESO模式降水预报评估区域平均有效降水频率定义为：对所选取区域内站点3 h累积降水量Rn,k进行累加平均，即为每日每个观测时段区域平均3 h累积降水量-Rn。然后对-Rn分无降水（-Rn0.1 mm）和有效降水（-Rn0.1 mm）2种情况进行统计，若-Rn0.1 mm，则有效降水次数On记为1，若-Rn0.1 mm，则On记为0。区域平均有效降水频率-fday为统计时段符合条件降水次数与统计时段总天数N之比。具体公式如下：-Rn=k=1KRn,kK（3）-fday=n=1NOnN（4）2夏季降水实况与预报特征分析2.1夏季3 h降水量及降水频率实况时空分布由2021年夏季国家级地面气象观测站逐3 h累积降水量（图1）和有效降水频率（图2）空间分布可知，站点3 h累积降水量和有效降水频率空间分布不均匀，南方多北方少、东部多西部少，下午至傍晚（06：0012：00）降水范围最广，且降水量明显大于其他观测时段；白天（21：00至次日09：00）中部地区（28.00N35.00N）沿长江流域有一条明显的雨带。东南沿海地区午后降水次数较多，同时段 3 h累积降水量远大于其他地区，部分地区降水达中雨以上量级，说明该地区主要以午后局地短时强降水为 主。距 离 临 近 的 川 渝 地 区 东 部（103.00E108.00E，28.00N33.00N）和青藏高原东部及四川西部（97.00E102.00E，28.00N33.00N）呈现出不同的降水特征，川渝地区东部降水主要集中在凌晨至上午（18：00至次日03：00），青藏高原东部和四川西部降水主要集中在深夜至清晨（15：0021：00），有效降水频率大。云南西南地区部分站点有效降水频率高于沿海地区，但降水量远小于东南沿海，说明云南西南地区降水具有区域性、连续性和降水强度小的特征。2.2模式预报夏季3 h降水时空分布图3是2021年夏季逐3 h累积降水量预报偏差图12021年夏季逐3 h累积降水量空间分布（单位：mm）Fig.1The spatial distribution of 3-hour cumulative precipitation in summer of 2021（Unit：mm）50541 卷干旱气象分布。可以看出，模式对不同地区的3 h降水预报表现不一致，以黑河腾冲为分界，降水量偏差东部较大、西部较小。大部分地区预报的3 h累积降水量高于观测值，而云贵高原和京津冀地区的降水量预报小于观测值，其中云南西南部少数站点的负偏差值较大，其余地区负偏差均在1.0 mm以内。对不同观测时段的预报偏差分析发现，降水偏差小值区在各观测时段的分布特征基本一致，而偏差大值区在各观测时段的分布略有不同。长江中下游地区降水预报偏差大值集中在白天（21：00至次日09：00），东南沿海和海南下午至傍晚（06：0012：00）的降水量预报偏差大，川西和青藏高原东部预报偏差大值集中在傍晚至凌晨（09：0015：00），而四川盆地及重庆地区偏差大值出现时间为凌晨至上午（15：00至次日03：00），且为正偏差。结合图1发现，预报偏差大值区与降水观测大值区基本对应。2021年夏季逐3 h累积降水量预报均方根误差空间分布（图略）与3 h累积降水预报偏差相似，均方根误差大值区偏差也较大，如东南沿海、青藏高原东部、海南及川渝地区，而长江以北内陆夏季降水量偏少地区，其偏差和均方根误差值均较小，说明模式不同预报时效预报均对无雨或弱降水有较好的预报效果。图4为2021年夏季逐3 h累积降水量观测和不同预报时效（3、6、9、12、15、18、21、24 h）预报值对比。可以看出，全国不同观测时段所有预报时效平均降水量均大于观测值，且随着预报时效增长，平均降水量均呈现先减少后增加趋势，其中3 h预报降水量最大，均比观测值大1.0 mm以上，即3 h预报降水量偏多，6、9、12 h预报降水量接近观测值。黄丽萍等（2017）指出 CMA-MESO模式预报的中雨以上量级偏差较大，可能是模式预报的水汽含量增加所致。朱立娟等（2017）发现采用云分析的短时降水对背景场有很大的依赖性，而云分析对背景场原因引起的虚假降水很难剔除。张小雯等（2020）通图22021年夏季逐3 h有效降水频率空间分布Fig.2The spatial distribution of 3-hour effective precipitation frequency in summer of 2021506第 3 期蔡怡等：2021年夏季CMA-MESO模式降水预报评估过评估CMA-MESO km模式预报降水发现，因加入云分析，系统spin-up问题不明显，3 h预报TS值最高。由此推断，云分析系统的不完善和模式预报的虚假云水信息可能是造成3 h降水量偏大的原因。2.3夏季3 h有效降水频率预报时空分布图5为2021年夏季逐3 h有效降水频率预报偏差空间分布。可以看出，06：0012：00全国大部分站点夏季3 h有效降水频率预报值高于观测值，说明该时段模式预报的降水日数比观测多。松花江流域东部不同观测时段夏季3 h有效降水频率预报均偏高，云南在凌晨至早上（18：00至次日03：00）有效降水频率偏低，川渝及青藏高原东部边缘地区的有效降水频率偏差在不同时段的分布与降水量预报偏差相似。2021年夏季逐3 h站点有效降水频率图42021年夏季逐3 h累积降水量观测和不同预报时效降水预报值对比Fig.4Comparison of observed 3-hour cumulative precipitation and forecast values with different timeliness in summer of 2021图32021年夏季逐3 h累积降水量预报偏差（预报减去观测）（单位：mm）Fig.3The forecast bias of 3-hour cumulative precipitation in summer in 2021（forecast minus observation）（Unit：mm）50741 卷干旱气象均方根误差大值区空间分布与有效降水频率偏差基本一致（图略）。不同预报时效的3 h有效降水频率接近 图6（a），除中午至傍晚（03：0012：00）3 h有效降水频率预报偏大外，其余时段预报与观测差距不大，说明模式对是否存在有效降水预报较好。强降水频率预报与观测差异明显 图6（b），尤其3 h预报，3 h强降水的频率预报远大于观测，说明模式3 h强降水的预报次数偏多，而6、9、12 h强降水频率预报接近观测。2.4不同区域夏季3 h降水量和降水频率特征选取前文分析的降水特征存在明显差异的5个区域做进一步分析：（1）川西区（97.00E102.00E，28.00N33.00N），即青藏高原东部边缘和四川西部 高 海 拔 复 杂 地 形 区；（2）川 东 区（103.00E108.00E，28.00N33.00N），即四川东部盆地区；（3）云南区（97.00E102.50E，21.50N26.00N）；（4）华南区（110.00E116.50E，21.50N26.00N）；（5）华东区（117.00E122.00E，28.00N33.00N）。综合2021年夏季中国不同区域平均3 h累积降水量及站点平均、区域平均有效降水频率（图7）分析，发现华东和华南地区午后降水量大，同时间站点平均有效降水频率小而区域平均有效降水频率大，表明华东和华南地区主要降水特点为午后局地强降水；川西和川东降水主要发生在深夜至清晨，同时间站点和区域平均有效降水频率都相对较大，说明川西和川东降水为区域性夜雨；云南降水量、站点和区域平均有效降水频率日变化相对较小，为区域连续性降水。不同地区不同预报时效预报的3 h累积降水量及站点平均、区域平均有效降水频率日变化（图8）显示，模式（除川西地区外）基本能预报出不同时效降水量的日变化特征，川西地区预报峰值比观测早约9 h。不同预报时效在不同区域表现不一，总体而言，平均预报效果优于单次预报。华南和华东 3 h降水量预报和站点平均有效降水频率远大于观测，表明3 h预报相对较差。从站点平均和区域平均有图52021年夏季逐3 h有效降水频率预报偏差空间分布Fig.5The spatial distribution of forecast bias of 3-hour effective precipitation frequency in summer of 2021508第 3 期蔡怡等：2021年夏季CMA-MESO模式降水预报评估效降水频率结果看，模式预报的区域平均有效降水频率更接近观测，说明模式能准确预报出区域降水，但精准性略差。云南区域性连续降水模式预报效果相对较好。川西和川东降水过程模式预报效果相对差一些，绝大部分时效预报的降水量比观测值大，其中川西地区预报的降水量和有效降水频次大值均出现在傍晚，和观测有明显差异，而川东地区预报的降水量和有效降水频次比观测高，但日变化趋势和观测较为一致。结合前文可知，西南地区降水局部差异大、变化复杂，受季风环流和复杂地形的共同影响，常有强降水发生，因此选取川西和川东地区2021年7月1724日持续一周的降水过程进行深入分析。从图9可以看出，模式3 h累积降水预报值基本都高于观测值，且川东地区降水预报效果优于川西地区。川东地区降水增加和峰值出现时刻预报相对准确，但3 h降水量预报偏大，6、9、12 h降水预报与观测最接近，这与夏季该地区平均降水日变化表现一致。图62021年夏季逐3 h有效降水（a）及强降水（b）频率观测和不同预报时效预报值对比Fig.6Comparison of observed 3-hour effective precipitation（a）and heavy precipitation（b）frequency and forecast values with different timeliness in summer of 2021图72021年夏季中国不同区域平均3 h累积降水量（a）及区域平均（b）、站点平均（c）有效降水频率的日变化Fig.7The diurnal variation of 3-hour cumulative precipitation（a），regional mean（b）and station mean（c）effective precipitation frequency in different regions of China in summer of 202150941 卷干旱气象图82021年夏季中国不同区域站点观测和不同时效预报的3 h累积降水量（左）及站点平均（中）、区域平均（右）有效降水频率日变化Fig.8Diurnal variation of observed and forecasted 3-hour cumulative precipitation（the left），station average（the center）and regional average（the right）effective precipitation frequency with different forecasting timeliness in different regions of China in summer of 2021510第 3 期蔡怡等：2021年夏季CMA-MESO模式降水预报评估川西地区夜雨较明显，模式对降水峰值出现时间把握不太准确，峰值出现时间较观测明显提前，7月19日提前了12 h。由此可见，模式有较好的降水预报能力，但对川西复杂大地形附近的预报还存在不足，其原因可能是川西地区地形复杂，多高原山地，因此需要进一步优化和完善该区域的资料同化、模式系统的动力框架和物理过程。3典型强降水过程预报性能评估选取 2021 年 3 个典型暴雨过程考察 CMA-MESO 3 km模式对短时区域强降水及单站极端降水的捕捉能力。3.17月1113日京津冀地区强降水2021年夏季京津冀地区（114.00E117.50E，37.00N41.40N）整体降水较少，7月1113日出现一次较强降水过程（图10），7月1112日，京津冀等地出现入汛以来最强降雨，大部地区出现暴雨到大暴雨。此次降雨具有影响范围广、累积降水量大、大风天气明显、局地雨强破极值等特点。模式较好地预报出了此次强降水过程，不同时效预报平均值与观测值基本一致，3 h降水预报略偏大，7月13日6、9 h预报与观测最接近。此次降水过程中，保定站 3 h累积降水量有 2次达到大暴雨量级，前3 h预报较好地抓住了降水极值出现时间，其预报趋势与观测基本一致，其余预报时效对于第一次降水峰值未预报出来，12 h及以上时效预报在11日第一次降水峰值出现了漏报，12日和13日出现了空报，可见对于大雨以上量级的短时强降水，应多关注9 h以内的短临预报结果。3.27 20郑州极端强降水2021年7月1722日，河南地区出现罕见的极端强降水天气。由7月1921日3 h累积降水量演变可知，此次降水主体分布在河南中北部和西部山地的东侧，郑州市及周边地区普遍出现短时超强降水，部分站点3 h累积降水量超过25.0 mm，郑州站出现该区域内降水极大值，3 h累积降水量最大达174.8 mm，打破建站以来历史记录。选取河南地区（113.00E114.50E，34.25N35.75N）及郑州站气象观测和模式预报（图11）进行分析。河南地区降水过程从7月19日06：00起，降水量级逐渐增加，21：00达到暴雨，至20日09：00降水量达到顶峰，强降水一直持续到21日21：00才逐渐减弱为大雨量级。不同预报时效降水预报平均值与区域3 h降水观测值趋势较为接近，均表现为先增后减。在3 d降水过程中，3 h降水量等级预报与实况基本一致，但大部分时段3 h降水量预报值略高于观测值，除极端降水出现时段预报量较观测值偏低20.0 mm外，其余均在6.0 mm以内。同时间段不同时效预报结果差异较大，最大差异可达60.0 mm，且降水量级越高，不同时效间预报结果差异越大。总体而言，3、6 h预报最接近观测。郑州站逐3 h累积降水量时序结果表明，不同图92021年7月1724日川西（a）和川东（b）地区3 h累积降水量观测值和不同预报时效预报值Fig.9The 3-hour cumulative precipitation of observation and forecast values with different timeliness in western（a）and eastern（b）Sichuan Province from 17 to 24 July 202151141 卷干旱气象图102021年7月1113日京津冀地区（a）及保定站（b）3 h累积降水量观测值和不同预报时效预报值Fig.10The 3-hour cumulative precipitation of observation and forecast values with different leading times in Beijing-Tianjin-Hebei（a）and Baoding station（b）from 11 to 13 July 2021图112021年7月1921日河南地区（a）及郑州站（b）3 h累积降水量观测值和不同预报时效预报值Fig.11The observed 3-hour cumulative precipitation and forecast values with different leading times in Henan（a）and Zhengzhou station（b）from 19 to 21 July 2021512第 3 期蔡怡等：2021年夏季CMA-MESO模式降水预报评估预报时效的模式预报均难以准确预报极端降水量值。3 h预报时效可以较为准确预报出郑州站的特大暴雨量级，短临（3、6 h）预报对降水趋势预报也相对比较准确，而其他预报时效的降水预报效果稍微差一些，且随着预报时效增加预报效果相对更差。由此可见，对于短时强降水过程，CMA-MESO 3 km模式的短临（3、6 h）预报效果更好。3.3台风烟花强降水为对比模式在不同天气系统下强降水过程的预报能力，选取2021年7月2426日台风“烟花”登陆前后降水过程进行分析（图12）。台风“烟花”移动速度慢、风雨强度大、影响范围广，长时间陆上滞留更是打破纪录。与郑州极端强降水过程相比，此次受台风影响的降水过程雨量略小。从区域（121.00E122.00E，30.75N31.50N）平均降水可知，这次降水过程模式预报效果较好，不同时效预报降水量平均值与观测值差异较小，且部分时段预报量级与观测量级一致，不同时效预报差异较小，都能很好地预报出降水量变化趋势。浦东站降水时序结果和区域降水平均结果较相似，3、6 h短临预报能更好地捕捉降水量短时变化。由此可见，CMA-MESO 3 km模式对不同降水过程的预报能力均较强，模式的区域降水预报能力强于单站，持续性降水预报效果好于局地短时强降水，局地短时强降水 3、6 h 等短临降水预报效果较好。4结 论本文利用2021年夏季（68月）国家级地面气象观测站每日00：0000：00逐3 h累积降水量观测资料，从多个角度对CMA-MESO 3 km模式降水预报进行评估。得到以下结论：（1）2021年夏季3 h降水量和有效降水频率时空分布不均，南方多北方少、东部多西部少，下午至傍晚（06：0012：00）降水范围广、降水量大。云南地区是连续区域性弱降水；东南沿海地区午后降水次数较多，以短时局地强降水为主；川渝地区多为区域性夜雨。（2）除青藏高原复杂地形附近川西地区外，模式可以较好地预报3 h降水量和降水频率的时空分布特征，区域降水预报能力强于单站，持续性降水预报效果好于局地短时强降水，且不同预报时效均对无雨或弱降水的预报较好。川西地区，模式降水图122021年7月2426日上海地区（a）及浦东站（b）3 h累积降水量观测值和不同预报时效预报值Fig.12The observed 3-hour precipitation and forecast values with different leading times in Shanghai（a）and Pudong station（b）from 24 to 26 July 202151341 卷干旱气象预报时段过长，降水峰值较观测提前约9 h。（3）CMA-MESO 3 km模式对不同降水过程的预报能力较强。不同预报时效、不同预报量级的降水预报效果存在差异。统计结果显示，对于非短时强降水过程，模式3 h降水量预报最大且远高于观测值，6、9、12 h预报结果接近观测值，因此可重点关注6、9、12 h的模式预报以及8个不同时效平均降水量预报。短时强降水个例结果显示，模式3、6 h降水预报更接近观测值，因此对于短时强降水需重点关注短临（3、6 h预报）。参考文献：曹越，赵琳娜，巩远发，等，2019.ECMWF高分辨率模式降水预报能力评估与误差分析 J.暴雨灾害，38（3）：249-258.陈静，庞波，吴政秋，等，2022.华南复杂地形下GRAPES_Meso 3 km 对流尺度模式前汛期精细化降水预报评估J.大气科学学报，45（1）：99-111.陈子通，黄燕燕，万齐林，等，2010.快速更新循环同化预报系统的汛期试验与分析 J.热带气象学报，26（1）：49-54.傅娜，陈葆德，谭燕，等，2013.基于快速更新同化的滞后短时集合预报试验及检验 J.气象，39（10）：1 247-1 256.宫宇，代刊，徐珺，等，2018.GRAPES-GFS模式暴雨预报天气学检验特征 J.气象，44（9）：1 148-1 159.国家气象中心，2012.降水量等级：GB/T 285922012 S.北京：中国标准出版社.黄丽萍，陈德辉，邓莲堂，等，2017.GRAPES_MesoV4.0主要技术改进和预报效果检验 J.应用气象学报，28（1）：25-37.李泽椿，谌芸，张芳华，等，2015.由河南 758 特大暴雨引发的思考 J.气象与环境科学，38（3）：1-12.刘雪晴，陈静，陈法敬，等，2020.降水邻域集合概率方法尺度敏感性试验 J.大气科学，44（2）：282-296梅钦，智协飞，王佳，2018.WRF模式不同云参数化方案的暴雨预报能力检验及集成试验 J.大气科学学报，41（6）：731-742.齐艳军，张人禾，LI T，等，2016.1998年夏季长江流域大气季节内振荡的结构演变及其对降水的影响 J.大气科学，40（3）：451462.史文茹，李昕，曾明剑，等，2021.“720”郑州特大暴雨的多模式对比及高分辨率区域模式预报分析 J.大气科学学报，44（5）：688-702.税军峰，任婧宇，彭守璋，等，2019.19012014年黄土高原1 km分辨率月均气温和月降水量数据集 J.中国科学数据（中英文网络版），4（4）：133-142.孙建华，卫捷，傅慎明，等，2018.江淮流域持续性强降水过程的多尺度物理模型 J.大气科学，42（4）：741-754.王婧卓，陈法敬，陈静，等，2021.GRAPES区域集合预报对2019年中国汛期降水预报评估 J.大气科学，45（3）：664-682.王叶红，彭菊香，公颖，等，2011.AREM-RUC 3 h快速更新同化预报系统的建立与实时预报对比检验 J.暴雨灾害，30（4）：296-304.韦青，代刊，林建，等，2020.20162018年全国智能网格降水及温度预报检验评估 J.气象，46（10）：1 272-1 285.吴晶，李照荣，颜鹏程，等，2020.西北四省（区）GRAPES模式降水预报的定量评估 J.气象，46（3）：346-356.吴秋霞，史历，翁永辉，等，2007.AREMS/973模式系统对2004年中国汛期降水实时预报检验 J.大气科学，31（2）：298310.谢漪云，王建捷，2021.GRAPES千米尺度模式在西南复杂地形区降水预报偏差与成因初探 J.气象学报，79（5）：732-749.徐同，杨玉华，李佳，等，2019.SMS-WARMS V2.0模式对中国西南地区降水预报能力的客观检验 J.气象，45（8）：1 065-1 074.徐枝芳，郝民，朱立娟，等，2013.GRAPES_RAFS系统研发J.气象，39（4）：466-477.徐致真，陈静，王勇，等，2019.中尺度降水集合预报随机参数扰动方法敏感性试验 J.气象学报，77（5）：849-868.许 晨 璐，王 建 捷，黄 丽 萍，2017.千 米 尺 度 分 辨 率 下GRAPES-Meso 4.0模式定量降水预报性能评估 J.气象学报，75（6）：851-876.薛纪善，庄世宇，朱国富，等，2008.GRAPES新一代全球/区域变分同化系统研究 J.科学通报，53（20）：2 408-2 417.翟盘茂，任福民，张强，1999.中国降水极值变化趋势检测J.气象学报，57（2）：208-216.张立凤，罗雨，2010.初始场对暴雨数值预报的影响及集合预报试验 J.气象科学，30（5）：650656.张小雯，唐文苑，郑永光，等，2020.GRAPES_3 km数值模式对流风暴预报能力的多方法综合评估 J.气象，46（3）：367-380.朱立娟，龚建东，黄丽萍，等，2017.GRAPES三维云初始场形成及在短临预报中的应用 J.应用气象学报，28（1）：38-51.庄照荣，王瑞春，李兴良，2020.全球大尺度信息在 3 km GRAPES-RAFS 系统中的应用 J.气象学报，78（1）：33-46.BENJAMIN S G，DVNYI D，WEYGANDT S S，et al，2004.An hourly assimilation-forecast cycle：the RUC J.Monthly Weather Review，132（2）：495-518.STENSRUD D J，BROOKS H E，DU J，et al，1999.Using ensembles for short-range forecastingJ.Monthly Weather Review，127（4）：433-446.514第 3 期蔡怡等：2021年夏季CMA-MESO模式降水预报评估STENSRUD D J，BAO J W，WARNER T T，2000.Using initial condition and model physics perturbations in short-range ensemble simulations of mesoscale convective systems J.Monthly Weather Review，128（7）：2 077-2 107.STENSRUD D J，FRITSCH J M，1994.Mesoscale convective systems in weakly forced large-scale environments.Part II：generation of a mesoscale initial conditionJ.Mon