评估一种局部化粒子滤波器与全天空FY-4 AGRI红外辐射同化的结合效果,用于模拟强降雨事件中的对流过程

《Atmospheric Research》:Evaluation of a localized particle filter with all-sky FY-4 AGRI infrared radiance assimilation for the convection-permitting simulation of an extremely heavy rainfall event

【字体: 时间:2026年02月18日 来源:Atmospheric Research 4.4

编辑推荐:

  本研究基于WRF-DART框架,利用支持全天空观测的CRTM模型,结合局部粒子滤波(LPF)方法同化FY-4A AGRI水汽通道数据。针对“23.7”极端降雨事件,实验表明LPF在非高斯、非线性问题中优于EnKF,降水预报评分提高5.3%-15.2%,且物理场协调性更优。

  
陈浩宇|王远兵|陈耀登|迪迪|孟德明|李欣|朱丽健
国家气候系统预测与风险管理重点实验室/教育部气象灾害重点实验室/南京信息科技大学气象灾害预报与评估协同创新中心,中国南京210044

摘要

风云-4(FY-4)先进地球静止轨道辐射成像仪(AGRI)提供了包含云结构和分布信息的高分辨率观测数据。先前的研究表明,对AGRI数据的全天同化可以提高数值天气预报的准确性。然而,其中涉及的非线性和非高斯特性仍然构成挑战。本研究基于天气研究与预报模型(WRF)和数据同化研究测试平台(DART),采用具有全局卫星观测支持能力的社区辐射传输模型(CRTM),开发了一个同化框架,并应用了局部粒子滤波器(LPF)方法来同化FY-4 A AGRI的水汽通道观测数据。针对中国北部的“23.7”极端强降雨事件,进行了一系列循环同化实验。敏感性实验表明,局部化和膨胀强度对LPF的性能至关重要。在这种情况下,LPF在全天观测条件下的表现优于晴空条件,并且在湿气调节和降雨模拟方面优于集合卡尔曼滤波器(EnKF),这体现了LPF在处理非高斯问题时的优势。具体而言,与EnKF相比,LPF在公平威胁分数和分数技能分数方面的降水技能得分分别提高了5.3%(24.8%)和2.6%(15.2%)(仅限晴空条件)。这些改进伴随着更协调的水分、热力学和动态条件,有利于降水的发展。这些结果展示了LPF在全天红外数据同化中的潜力。

引言

卫星红外辐射(IR)观测由于其高空间和时间分辨率,以及捕捉云系垂直结构和水平分布的能力,在全球和区域数值天气预报(NWP)中发挥着重要作用(Le Marshall等人,2006;Li等人,2019;Li等人,2022;Wang等人,2025)。值得注意的是,多云和降雨区域的卫星观测通常包含与恶劣天气发生和演变密切相关的信息,而预报敏感区域也与云覆盖度密切相关(McNally,2002;Geer等人,2018)。这意味着充分利用多云区域的红外观测可以提高数值模型表示天气过程的能力(Honda等人,2018;Wang等人,2022)。因此,如何有效同化此类观测数据已成为数据同化(DA)中的一个重要问题。
然而,卫星IR观测的全天同化仍然面临许多困难,包括数值模型和辐射传输模型在多云区域的不确定性,以及由云相关物理过程引入的非高斯和非线性问题(Errico等人,2007;Kumar和Shukla,2019)。集合卡尔曼滤波器(EnKF)和变分方法都依赖于高斯先验分布和线性化假设。因此,当处理具有显著非线性和非高斯特性的问题时,它们可能存在局限性,尤其是在对流尺度的高网格分辨率下(Van Leeuwen等人,2019)。
与上述方法依赖于先验概率密度函数的高斯假设和观测算子的线性不同,粒子滤波器(PF)基于蒙特卡洛方法和贝叶斯定理。它使用一组粒子的加权狄拉克δ函数之和来近似集合先验分布,无需对先验分布或观测似然进行高斯假设。该方法的核心思想是根据观测似然为每个粒子分配权重,然后通过重采样过程丢弃低权重粒子并复制高权重粒子来更新集合。通过避免强参数假设,当粒子数量足够大时,该方法在大多数应用中表现良好(Doucet等人,2001;Poterjoy和Anderson,2016),这使得PF成为解决DA中非高斯和非线性问题的一个有吸引力的候选方法。
然而,尽管PF在理论上有优势,但其应用在DA中仍然受到限制。例如,这种方法通常需要一个极大的粒子集合,且集合大小往往随着系统维度的增加而呈指数级增长(Snyder等人,2008),这意味着在高维系统中计算上不可行。当集合大小不足时,可能会发生滤波器发散,即粒子权重逐渐集中在少数几个粒子上,从而降低集合的代表性。这长期以来一直是限制PF应用于高维地球科学问题的关键因素之一(Poterjoy,2016)。
最近,随着计算资源的进步和各种解决PF滤波器发散问题的策略的发展,如局部化(Van Leeuwen,2009;Poterjoy,2016;Potthast等人,2019)、提议密度方法(Doucet等人,2001)以及与EnKFs的混合同化(Frei和Künsch,2013;Potthast等人,2019;Poterjoy,2022),PF在NWP模型中的可行性明显提高。其中,引入局部化以将观测的影响限制在局部区域,从地球科学应用的角度来看是一种更直接的方法(Van Leeuwen等人,2019)。其目的是在同化大量数据时抑制观测之间的虚假相关性。Penny和Miyoshi(2016)提出了一种基于观测空间局部化的局部粒子滤波器,并使用40变量Lorenz-96(L96)模型证明了其在非高斯条件下的优势。在另一种方法中,Poterjoy(2016)引入了局部粒子滤波器(LPF),该滤波器将粒子重采样权重向量化并加入了局部化函数,使用L96模型表明该方法能有效避免滤波器发散,并在非线性情况下表现优于EnKFs。随后,Poterjoy等人(2019)进一步提出了LPF的改进措施,包括优化的权重计算和概率映射策略,以及在新观测稀疏但准确的情况下增强滤波器稳定性的新方法。使用L96模型和天气研究与预报(WRF)模型的数值实验证明了这些改进的积极效果,并突出了LPF在高维模型中的潜力。在此基础上,Poterjoy(2022)引入了正则化和自适应温度调节策略,进一步增强了LPF在采样误差较大的区域中的有效性。
中国新一代地球静止轨道卫星风云-4(FY-4)上的先进地球静止轨道辐射成像仪(AGRI)可以提供高空间和时间分辨率的大气观测数据(Yang等人,2017)。先前的研究表明,同化FY-4 A AGRI IR数据可以有效提高NWP的准确性。具体来说,Zhang等人(2022)使用3D-Var在同晴空条件下同化了FY-4 A AGRI水汽通道数据,从而获得了对台风“In-fa”的改进预报。随后,Xu等人(2023)首次尝试在全天条件下使用3D-Var同化FY-4 A AGRI水汽通道数据,获得了更接近观测结果的湿度分析场,并提高了模型的降水预报能力。最近,Zhong等人(2026)在不同同化频率下进行了FY-4 A AGRI全天DA实验,发现每小时同化频率的预报性能更优,因为它有效减少了由于频率过快或过慢造成的不协调问题。然而,尽管越来越多的研究关注AGRI同化,但大多数研究仍然依赖于基于高斯假设的变分或EnKF方法。
现有的LPF研究主要集中在传统观测或雷达数据上。近年来,一些研究开始探索将PF方法应用于与卫星相关的DA问题。例如,Kumar和Shukla(2019)应用PF处理卫星观测数据,以选择适当的模型物理参数化方案,从而减少与模型物理相关的不确定性。Shen等人(2024)开发了一种基于PF的云检测方案,以便于识别适合同化的晴空辐射。此外,Zhou等人(2023)在OSSE框架内使用LPF同化卫星可见辐射。总体而言,这些研究要么将PF视为辅助工具而不直接同化卫星辐射观测数据,要么专注于在模拟观测条件下进行的DA实验。
与以往的工作不同,本研究是首次尝试直接使用LPF在同流尺度模拟中同化全天真实红外亮度温度(BT)。以中国北部的“23.7”极端降雨事件为例,进行了一系列实验,以检验LPF对局部化半径和膨胀强度的敏感性,并将其性能与EnKF进行比较,从而展示了LPF在全天条件下的分析和预报场中的影响。
本文的其余部分组织如下。第2节介绍了本研究中使用的AGRI辐射数据和全天同化框架。第3节介绍了卫星数据的预处理和LPF方法。第4节介绍了选定的案例、区域模型和DA系统的配置以及实验设计。第5节展示了实验结果和分析,第6节给出了本研究的主要结论和讨论。

节选

AGRI辐射数据

FY-4卫星是中国第二代地球静止轨道卫星系列,目前包括FY-4 A(2016年12月11日发射)和FY-4B(2021年6月3日发射)。AGRI是FY-4上的主要仪器之一,主要负责云成像,同时提供水汽、气溶胶和积雪覆盖的信息。FY-4 A AGRI共有14个通道,其中包括2个对上层和中层敏感的红外水汽通道(通道9和10)

观测误差和辐射数据预处理

AGRI全天IR观测包含大量受云影响的信息。与晴空观测相比,全天观测具有更大的不确定性,更容易将随机和系统误差引入模型,从而增加同化的复杂性(Harnisch等人,2016)。全天红外DA中的一种常见策略是根据云的影响分配观测误差(Stengel等人,2013)。在这里,我们采用了云影响量化方法

案例概述

本研究选择了2023年7月29日至8月1日发生在中国的“23.7”极端强降雨事件。该事件期间最强烈的降水集中在7月29日至30日,7月30日的降雨尤为强烈。7月29日1800 UTC至7月30日0600 UTC的环流演变由欧洲中期天气预报中心(ECMWF)再分析v5(ERA5)显示在图4中。台风Doksuri减弱后,其残余

对局部化半径的敏感性

LPF通过将欧几里得距离函数纳入粒子权重来限制观测的影响范围。然而,应该注意的是,这种方法对采样误差非常敏感,如果局部化区域包含太多独立观测,则可能会发生滤波器退化(Poterjoy等人,2019;Van Leeuwen等人,2019)。这个局部化区域的范围由局部化半径定义。因此,确定一个

结论与讨论

鉴于全天红外卫星观测具有显著的非线性和非高斯特性,有效同化此类数据仍然是一个重要且具有挑战性的问题。局部粒子滤波器(LPF)作为一种不依赖于高斯假设的方法,在解决这一问题方面具有理论优势。在本研究中,我们选择了中国北部的一次强降水事件作为案例,并构建了一个能够

CRediT作者贡献声明

陈浩宇:撰写——原始草稿、可视化、软件、方法论。王远兵:撰写——审阅与编辑、方法论、概念化。陈耀登:监督、方法论。迪迪:软件。孟德明:软件、方法论。李欣:资源、数据管理。朱丽健:软件、方法论。

未引用参考文献

Han等人,2006
Robert等人,2018

利益冲突声明

作者声明他们没有已知的竞争性财务利益或个人关系可能影响本文报告的工作。

致谢

本工作得到了国家自然科学基金(U2142208,42475158)和国家重点研发计划(2022YFC3080500)的共同支持。我们还要感谢南京信息科技大学的高性能计算中心对本工作的支持。
相关新闻
生物通微信公众号
微信
新浪微博
  • 搜索
  • 国际
  • 国内
  • 人物
  • 产业
  • 热点
  • 科普

知名企业招聘

热点排行

    今日动态 | 人才市场 | 新技术专栏 | 中国科学人 | 云展台 | BioHot | 云讲堂直播 | 会展中心 | 特价专栏 | 技术快讯 | 免费试用

    版权所有 生物通

    Copyright© eBiotrade.com, All Rights Reserved

    联系信箱:

    粤ICP备09063491号