en
×

分享给微信好友或者朋友圈

使用微信“扫一扫”功能。
目录contents

    摘要

    邻域法在天气预报中有着广泛的应用,其关键应用领域包括2个方面,一是基于邻域法的高分辨率数值模式检验,二是邻域概率或者集合预报的邻域概率。首先,回顾了邻域法“一对多”和“多对多”的两种邻域法检验框架,归纳了邻域法数据处理方法、常用评分指数的物理意义。其次,总结了网格尺度上的邻域概率和大于网格尺度邻域概率的基本思想和统计意义,重点阐述了与集合预报相结合产生的邻域集合概率(neighborhood ensemble probability,NEP)预报、邻域最大集合概率(neighborhood maximum ensemble probability,NMEP)预报的算法流程和内在含义。第三,进一步结合典型应用个例,分析了邻域法检验和邻域集合概率的优缺点和适用性。总体来说,邻域法检验可以在不同的时空尺度上比较预报产品的性能,具有独特的优势。虽然NEP和NMEP两种邻域概率都可以提高降水的预报评分,但NEP更适合于大尺度、系统性降水预报,NMEP对对流性、极端性降水有更好的应用效果。最后,给出了使用邻域法应注意的问题以及未来研究应用的发展方向。

    Abstract

    The traditional dichotomous contingency table test, which evaluates the objective performance of numerical weather prediction (NWP) based on the point-to-point matching between forecasted and observed events, has notable limitation when applied to high-resolution NWP or convection-allowing models (CAM). The neighborhood method addresses these limitations by relaxing the grid scale matching constraints between forecasted and observed events, making it particularly valuable for evaluating high-resolution numerical weather forecasts and the post-processing of objective probability forecasts. This paper systematically reviews the key applications of the neighborhood method in weather forecasting, focusing on two key aspects: one is the verification of high-resolution numerical models using neighborhood method; and other is the neighborhood probability or neighborhood probability of ensemble forecasts. First, the study outlines the verification frameworks of two neighborhood methods,“one-to-many” and “many-to-many”, and discusses the data processing techniques associated with the neighborhood method, alongside the physical interpretation of common scoring matrices such as FBS (fractions brier score) and FSS (fractions skill score). It is concluded that, in addition to traditional dichotomous contingency table-based verification metrics, the neighborhood method facilitates comparisons of forecast performance across multiple spatial and temporal scales. This enables the derivation of diagnostic metrics for NWP forecast performance based on scale changes, providing unique advantages. Second, it summarizes the fundamental concepts and statistical meaning of the grid scale neighborhood probability and the neighborhood probability at scales larger than the grid. Discussion focuses on expounding the algorithm workflow and internal meaning of neighborhood ensemble probability (NEP) forecast and neighborhood maximum ensemble probability (NMEP) forecast derived from ensemble forecasts. Third, by examining typical application cases, it analyzes the advantages, disadvantages and applicability of the neighborhood method and neighborhood ensemble probability. Results show that both NEP and NMEP enhance precipitation forecast scores. NEP performs better for large-scale and systematic precipitation forecasts, whereas NMEP is more effective for convective and extreme precipitation events. However, the selection of an appropriate neighborhood radius remains a critical technical challenge, as it is influenced by variations in underlying surface conditions and the optimal neighborhood scales of different NWP products. Finally, the paper discusses future directions for the application of the neighborhood method in weather forecasting. Promising areas of research and application include integrating neighborhood ensemble probability with the temporal dimension, developing metrics for the rare-event ensemble neighborhood probability, and exploring synergies between the neighborhood method and artificial intelligence. These directions hold significant potential for advancing the utility and impact of the neighborhood method in weather forecasting.

  • 在过去的几十年中,随着计算资源的快速增加和数值模式的不断发展,数值天气预报(numerical weather prediction,NWP)的水平网格间距越来越小,模式中的积云参数化方案逐渐被能够更加精细、准确描述小尺度大气动力框架的对流可分辨模型(convection allowing model,CAM)所替代。尽管采用CAM的NWP能够比更粗分辨率的对流参数化方案产生更加清晰的对流结构,但传统检验方法通常基于预报与观测事件点对点的匹配来评价数值模式的预报表现,这种检验往往导致客观检验结果与预报员对模式的主观评判不尽相同(Mass et al.,2002; Weisman et al.,2008; 苏翔和袁慧玲,2020; 吴瑞姣等,2020; 张小雯等,2020; 刘侃等,2023)。

  • 空间检验方法是伴随着高分辨率数值预报的发展而发展起来的(Mittermaier and Roberts,2010; Duc et al.,2013; 潘留杰等,2014,2023),发展空间检验方法最初的动力在于平衡客观检验结论和预报员主观判定之间的差异。邻域法、尺度分离法、目标对象法和变形场法是最主要的高分辨率数值模式空间检验方法。与许多空间检验方法类似,考虑到模式对小尺度天气事件的预报能力有限,在网格尺度上非常准确地预报小尺度天气事件并不现实,邻域法通过放宽预报和观测事件在网格尺度上的匹配限制来评价模式的预报表现(潘留杰等,2015,2016a,2017a; Kochasic et al.,2017; Johnson et al.,2020; 栗晗等,2022)。这种将单点预报放到一个较大的背景空间中进行检验的方法,不仅可以加深对数值模式预报性能的理解,而且可以为预报员提供更加全面的视角。

  • 考虑到邻近格点要素的预报信息能够相互补偿,从单一确定性预报可以生成概率预报,Theis et al.(2005)给出了从确定预报生成概率预报的技术指南,并被广泛应用(Ben Bouallègue et al.,2013; Schwartz et al.,2015)。由于其方法直观、物理概念清晰、易于理解和实现,邻域法概率预报很快扩展到集合预报领域(Schwartz et al.,2010)。基于集合预报的成员预报导出每个网格点的概率预报,在此基础上进行平滑,平滑过程可以看作是在概率计算中包含所有成员的空间邻域预报上来扩大集合样本量的计算方法。在集合概率预报基础上按照一定的空间尺度进行平滑的方法称之为网格尺度或模糊逻辑邻域概率(fuzzy logic neighborhood probability)。Ben Bouallègue and Theis(2014)的研究表明,平滑对预报技巧具有积极的影响,特别是在局地小尺度天气事件的可靠性方面表现较好,但在分辨率方面有一定的损失。

  • 天气预报的准确率通常与空间尺度、时间窗等有关,区域范围(Ebert and McBride,2000; Ebert,2008)、时间间隔(Gilleland et al.,2009,2010; 张宏芳等,2014)、点与面之间的关系(Golding,2000; Germann and Zawadzki,2004; Hagen-Zanker et al.,2005; Marsigli et al.,2008)等都会对预报结论产生影响。邻域集合预报概率随定义的空间尺度不同有明显的差异。平滑能够增加预报的可靠性,但不可避免地降低空间分辨率,当预报概率为接近零的低概率或者100%的高概率时,这种做法通常是正确的(Murphy et al.,1980; Nachamkin and Schmidt,2015)。然而,在多数情况下,当概率达到一定的阈值时,就需要发出预警或者给出确定性的预报(Johnson and Wang,2012)。因此产生了另一种大于网格尺度的邻域概率——“升尺度”(upscaling)技术(Clark et al.,2010)。升尺度通过改变预报与输出要素的空间尺度来提高预报能力。空间尺度内通过集合预报成员构建邻域概率预报,邻域概率是该区域内任意点邻域尺度上事件的发生概率,同时能够在区域内保留细网格数值模式预报的优点,具有独特的优势(Clark et al.,2011)。

  • 随着NWP的时空分辨率越来越高,适用于高分辨率数值模式检验的邻域法将会得到更多更广泛的应用。本文从邻域法的基本概念出发,系统地回顾和总结了邻域法检验、邻域概率、集合预报的邻域概率的研究进展和一些典型的应用,详细讨论了邻域法检验和邻域概率应用中的注意事项,以期为更好地在天气预报中应用邻域法提供参考。

  • 1 邻域法简介

  • 1.1 邻域检验框架和确定性预报的邻域概率

  • 邻域法检验是模糊检验的一种,它假定预报的空间位置和观测不完全匹配,但空间位置误差在可接受的范围内,预报仍有意义,即定义这种匹配情况下的检验为邻域法检验。允许位移偏离的程度或大小即为邻域半径,邻域半径内部的区域即为邻域窗。图1给出了传统和邻域检验的概念示意图,可以是“一对多”的一个站点观测匹配周围多个预报格点,也可以是一个格点预报匹配邻域窗内的多个站点,其基本思想是以单点观测或预报为中心搜寻邻域内所有格点的预报或观测(Barthold et al.,2015)。邻域窗可以是圆形(Roberts and Lean,2008)。邻域半径的大小取决于网格间距、时间分辨率和气候背景,因此邻域大小不是一个确定的值,而是变化的,在邻域检验中,通常以格点为单位采用“穷举法”计算预报评分,进而获得最优邻域半径,得到的检验结果通常是邻域半径的函数,其物理意义表示预报相对于观测的空间位置的偏移程度。

  • 观测和预报也可同时在给定的邻域窗内进行匹配,称之为“多对多”邻域检验。这种方法的优点在于去除了观测的不确定性。“多对多”邻域检验最常见的是NM(neighborhood maximum)匹配方案,它根据邻域窗内的最大值是否满足事件标准,定义该格点处的事件发生率(Ben Bouallègue and Theis,2014; Barthold et al.,2015)。如果预报和观测到的事件都发生在半径为r的邻域窗内的任何地方,则记为命中; 如果预报事件发生在r邻域窗内,但观测事件未发生,则为空报; 当邻域窗内没有发生预报,但出现了观测事件,则为漏报。

  • 基于邻域检验匹配方案可以从单一确定性预报生成邻域概率预报。邻域概率检验方法与惯常的概率预报检验方法基本一致。将单一确定性预报转换成概率预报,首先是格点“单点事件标记”,即设定If表示单个格点的预报,根据If是否满足阈值条件转换成1、0预报; 其次是基于邻域窗内的“单点事件标记”,计算邻域窗内以某个格点为中心预报事件发生的概率,即为该点的邻域概率。用〈〉s表示邻域半径为s的邻域窗内的所有值,假定观测数据和预报数据是均匀网格且格距匹配,对于确定性预报邻域半径s内以某个格点为中心的邻域概率〈Pfs则为:

  • 图1 邻域法检验匹配示意图:(a)观测场(黑色方块表示观测站点);(b)传统检验“点对点”和邻域“一对多”检验匹配(黑色方块表示与(a)中“点对点”检验匹配的预报格点,灰色填充表示以3个格点为半径的邻域预报与(a)中“一对多”的匹配);(c)和(d)分别表示邻域半径为3个格点的观测场和预报场的“多对多”的邻域检验匹配

  • Fig.1 Schematic representation of the neighborhood verification method: (a) observation field (black square represents observation site) ; (b) comparison of traditional “point-to-point” and neighborhood “one-to-many” verification approaches (black square represents forecast grid point for “point-to-point” matching, and gray shadings indicate neighborhood of forecasts with a radius of three grid points for“one-to-many” matching; (c) “many-to-many” neighborhood verification approach applied to the observation field, with a neighborhood radius of three grid points; (d) “many-to-many” neighborhood verification approach applied to the forecast field, with a neighborhood radius of three grid points

  • Pfs=1nIf
    (1)
  • 与预报类似,假定Io表示单个格点的观测值,则邻域窗内观测事件发生的频率〈Pos为:

  • Pos=1nIo
    (2)
  • 图2给出了一个暴雨预报的邻域概率转换个例。假定邻域半径设为1,包括中心格点的邻域窗内的格点数为9,在降水分布为图2a的情况下,暴雨预报的单点事件标记为图2b,中心点暴雨单点概率为零,但受周边格点降水的影响,邻域概率为0.33(图2c),仍然有可能出现暴雨。图2c没有计算受边界影响格点的邻域概率,当受边界影响时,参与邻域概率计算的格点数可能会减少。值得注意的是,多大的邻域窗能够获得最好的预报表现,是邻域概率计算的关键因素,对于一个给定的研究区域,基于“穷举法”计算不同邻域半径的预报评分,当增加邻域半径,预报评分不再增加或开始减小,此时的邻域半径即为最优邻域半径。

  • 利用邻域概率预报和观测频率,可检验概率预报的Brier评分(Murphy,1973; Kharin and Zwiers,2003; 潘留杰等,2016b,2017b)或其他概率预报评分。设定邻域半径为s的预报〈Fs和观测〈Os的误差表现为〈Es。邻域法检验的主要框架为:选择半径为ss=1,2,···,S)的邻域和强度为kk=1,2,···,K)的阈值,计算邻域检验结果。这包括2个步骤:对于每一个格点,统计邻域窗内的所有预报,如果观测也是格点的“多对多”检验,同样统计每个格点周围的观测; 对于每个强度阈值,计算强度——尺度统计量的确定性或概率的检验评分。检验结果为随尺度S和强度K而变化的K×S矩阵。

  • 1.2 邻域法检验评分

  • 所有的列联表二分法检验评分都可以基于邻域法进行拓展。主要的对比指标有邻域窗内数值平均(Ebert,2009)、最小覆盖(Yates et al.,2006)、模糊逻辑(Ebert et al.,2003)、组合阈值(Atger,2001)和占比技巧(Roberts and Lean,2008)评分等。其中占比技巧应用最为广泛,该方法直接比较邻域窗中观测和预报满足阈值事件的覆盖率。对于确定性模式,如果预报事件的概率与观测事件的频率趋于一致,则预报是有用的。具体实现时,不是直接比较邻域事件发生的频率,而是计算分数覆盖比率评分(fractions brier score,SFB)。SFB定义如下:

  • SFB=1Nn=1N Pfs-Pos2
    (3)
  • 式中:N是整个区域内以S为尺度的邻域窗的个数; 〈Pfs和〈Pos分别表示邻域窗内预报和观测超过给定阈值事件的概率。由于SFB是一种负向评分指数,即SFB越小,模式预报技巧越高,所以对SFB进行转换,即为FSS(fractions skill score,SFS):

  • SFS=1-SFB1Nn=1N Pfs2+n=1N Pos2
    (4)
  • 式中:〈Pfs和〈Pos含义与式(3)相同。分母表示最坏情况下的预报,即预报和观测的区域没有重叠,SFS的值域范围是0~1。

  • 图2 邻域概率计算示意图:(a)格点降水;(b)单点事件标记;(c)邻域概率

  • Fig.2 Schematic diagram illustrating the neighborhood probability method: (a) grid based precipitation distribution; (b) event making single grid point; (c) probability calculated across the neighborhood

  • 2 集合预报的邻域概率

  • 基于确定性预报,可以从某一个格点与其周边相邻格点的空间关系来建立要素场中事件发生概率的大小。但邻域概率更为广泛的应用是在集合或集成预报多成员的基础上,综合考虑成员预报和要素的空间位置关系,构建集合预报的邻域概率,通常被称之为邻域集合概率。

  • 2.1 网格尺度的邻域概率

  • 集合预报的成员预报在数据处理上相当于一个确定性预报。假定q表示降水量阈值,fij是集合预报系统第jj=1,2,···,N)个集合预报成员的第ii=1,2,···,M)个格点的降水预报,那么对集合预报第j个成员的单点事件标记PqB,ij可简单表示为:

  • P(q)B,ij=1, if fijq;0, if fij<q
    (5)
  • 式中:PqB,ijq的函数,通过选择邻域窗将PqB,ij转换成概率预报。定义第i个格点周围总格点数为NbPqB,ij可以转换为邻域概率PqN,ij:

  • P(q)N,ij=1Nbk=1Nb P(q)B,ij,kSi
    (6)
  • 式中:Si表示第i个格点的邻域窗内的点集。从PqB,ijPqN,ij需要在邻域窗内求平均,由于平均有平滑效果,所以邻域的长度可称作平滑尺度r。严格地说,PqN,ij是集合预报第j个成员中第i个格点邻域窗内预报事件超过阈值的比例,PqN,ijr密切相关,r的大小直接影响概率PqN,ij。所有集合预报成员邻域概率的平均,即为集合预报的邻域概率(neighborhood ensemble probability,NEP,PNE),也被称为模糊逻辑概率(Schwartz et al.,2010)。第i个格点,降水量阈值为qPNE的计算式为:

  • P(q)NE,i=1Nj=1N P(q)N,ij
    (7)
  • PqN,ij一样,PNEq的函数,受邻域半径r的控制。第i个格点的PNE是对所有网格点上的集合预报成员概率PqB,ij计算平均获得的集合预报成员概率PqE,i:

  • P(q)E,i=1Nj=1N P(q)B,ij
    (8)
  • 邻域窗为Nb,第i个格点,降水量阈值为q的集合预报的邻域概率PNE也可以通过下式获得:

  • P(q)NE,i=1Nbk=1Nb P(q)E,k,kSi
    (9)
  • 两种方法获得的集合预报的邻域概率在数学上是相同的,但式(8)和式(9)的计算效率更高。

  • 2.2 大于网格尺度的邻域概率

  • 在集合预报框架内,选取邻域窗内预报成员集合平均的最大值(Schwartz and Sobash,2017)可被看作是一种升尺度,通过该方法能够获得集合预报邻域最大概率(neighborhood maximum ensemble probability,NMEP,PNME)。需要注意的是,NMEP对集合预报成员进行了平均,但NMEP在空间上是离散的。通过平滑可以将离散的NMEP转换成连续的概率场,在一定程度上消除细小的相位误差,提高了预报评分(Skinner et al.,2016)。高斯滤波是常用的平滑方式,转换公式(Schwartz and Sobash,2017)为:

  • PS,NME,i=m=1M 12πσ2exp-12xi-mσ2PNME,i
    (10)
  • 式中:PS,NME,i表示平滑后的PNME,i; xi-m是第i个点到第m个点的距离; 标准差σ是一个可调整的控制平滑权重的尺度。

  • 图3给出了集合预报的邻域概率示意图。对于Nm个成员的集合预报来说,每一个成员预报按照阈值转换成单点概率(图3a),Nm个成员的集合预报概率通过对单点概率计算平均值,即为集合预报概率(图3b)。单点概率计算的集合预报概率并不会受到侧边界的影响。邻域集合概率PNE为邻域半径内集合预报概率的平均值(图3c),但研究区域的侧边界会受到影响。以9个格点的邻域窗为例,中间区域邻域窗内计算时基于完整9个格点,方形邻域窗的4个顶点邻域概率计算基于4个格点,其余为基于6个格点获得的邻域概率。与PNE类似,PNME(图3d)也会受到侧边界的影响。

  • 2.3 自适应邻域概率

  • NEP和NMEP通常使用固定邻域半径,从而忽略了集合传播中的地理和时间变化。Blake et al.(2018)提出了集合一致性尺度(ensemble agreement scale,EAS)技术,当局部一致性减小(增大)时,通过增加(减小)邻域半径来局部调整邻域窗的大小。

  • Dij=(A-B)2(A+B)2, if A>0 and B>0;1, if A=0 or B=0
    (11)
  • 式中:Dij为两个成员之间的一致性或差异性; AB为超过阈值大小的平均值。较小的Dij表示具有较大的一致性,通过设定阈值Dcrit,ij来进行判定,有:

  • DijDcrit ,ij,Dcrit ,ij=α
    (12)
  • 式中:α反映了可容忍的偏差程度。α=0表示在网格尺度上不容忍任何偏差; α=1表示容忍任何偏差。

  • 3 邻域法应用

  • 3.1 邻域法检验

  • 邻域法检验不仅可以提供位移误差和预测技能的真实评估,而且可以将不同分辨率的空间预报场进行比较,从而减少高分辨率模式预报中的“双重惩罚”问题(Mittermaier and Roberts,2010; Mittermaier,2014)。按照预报评分属性的不同,邻域法检验可以分为确定性要素检验和概率预报检验,按照要素预报类型的检验可以分为降水预报检验和其他要素检验,按照要素出现的频次或稀有性可以分为常规要素和极端灾害性天气的检验。当邻域窗较大或比较小时,邻域法检验评分对预报偏差敏感(Duc et al.,2013; Mittermaier,2021; Rempel et al.,2022; Pan et al.,2024),研究区域足够大是邻域法检验的重要基础,如果研究区域太小,则可能无法测量位移误差的全部尺度。此外,针对不同形状邻域窗的检验评估差异,也有一些分析和比较(Nachamkin and Schmidt,2015); 研究表明,圆形和方形邻域窗对检验结果的影响不敏感。

  • 图3 邻域集合概率概念图(N1N2、···、Nm表示集合预报成员;(c)和(d)中概率标记为黑色的表示邻域窗内包括9个格点,蓝色和红色分别表示受边界影响,邻域窗内分别为6个和4个格点):(a)集合成员的单点事件标记;(b)集合预报概率;(c)邻域集合概率;(d)邻域最大集合概率

  • Fig.3 Conceptual depiction of neighborhood ensemble probability methods (N1, N2, ···, Nm represent the ensemble forecast members) : (a) event marking for individual grid points across ensemble members; (b) probabilities of ensemble forecasts; (c) probability of the ensemble forecasts aggregated over a neighborhood window, with 9 grid points, contributing to probabilities marked in black, 6 grid points for boundary-affected probabilities marked in blue, and 4 grid points for boundary-affected probabilities marked in red; (d) maximum neighborhood ensemble probability

  • 传统评分和FSS的联合检验(屠妮妮等,2022)、高分辨率模式降水的空间尺度预报表现(唐文苑等,2018; 刘静等,2019; 唐文苑和郑永光,2019; 李子良等,2021; 俞碧玉和朱科锋,2022)、风速和温度的预报效果分析(张博和赵滨,2019; 林晓霞等,2021)是邻域法检验在国内应用的主要方面。图4给出了基于ECMWF模式预报的一个邻域法FSS和ETS(equitable threat score)检验评分的典型应用。个例中的观测资料为与预报时段一致的CMORPH(NOAA Climate Prediction Center morphing method)卫星与自动站逐小时降水融合资料。由图4可见,在低阈值降水时,随着尺度的增大,FSS和ETS评分同时增大,而在高阈值降水时,增大空间尺度尽管有可能提高面积预报准确率,但使得降水强度预报性能急剧下降。赵滨和张博(2018)研究发现,合适的邻域半径是邻域法检验的核心问题,半径过小难以获得最优可用预报尺度,半径过大则会将不同预报系统合并检验引起评估误导。对于高分辨率模式,随着尺度的增加FSS逼近于2b/(b2+1),这里b为整个研究区域内预报和观测超过给定阈值事件概率的比率。

  • 在预报相对于观测出现偏移的情况下,传统的预报评分不能区分预报的好坏程度。FSS作为一种空间检验评分,能够克服这种缺陷。Skok and Hladnik(2018)通过一个理想化的风场预报试验方案进行了证明。在理想化的风场试验(图5)中包括500×500个格点,实况环境风场为一致的1 m/s西风,假定受山脉或地形的影响,椭圆区域的上半部分和下半部分分别有西南阵风和西北阵风(图5a)。预报的阵风强度以及背景风场与实况完全一致,但阵风的空间位置相对实况的位移逐渐从零增加到100个格点(图5b)。通过计算传统的预报评分RMSE和FSS评分可以看到,RMSE和FSS是位移的函数,当位移为0时,RMSE为0,FSS为1。随着位移的增大,RMSE增大,FSS减小,即预报评分变差。然而,一旦位移达到100个网格点时,RMSE停止增加,表明当空间位移较大时,传统检验评分无法区分不正确预报事件的好坏程度。但FSS分值随着位移的增加继续减小,远超过100个网格点的位移,表明FSS评分即使在空间位移较大的情况下也能区分预报的优劣。

  • 为了减少高分辨率模式相对于全球模式的“双重惩罚”效应,在邻域窗内基于列联表的误差补偿策略被提出(Stein and Stoop,2019),即位于同一邻域窗中的漏报事件对空报、命中和正确预报未出现事件可以相互补偿。这相当于说,如果预报和观测到的相同量级的事件之间的位置误差小于邻域大小,则可以忽略或容忍,从而补偿高分辨率模式的预报相位误差,但补偿过程中观测和预报事件的频率保持不变。对于经典的二分法列联表T事件分类有:

  • 图4 邻域法检验应用示例:(a)不同空间尺度和降水量级上的邻域法FSS;(b)ETS评分检验(引自潘留杰等(2016a))

  • Fig.4 Example application of the neighborhood method for verification: (a) FSS across varying spatial scales and precipitation thresholds; (b) verification of the ETS score (cited from Pan et al. (2016a) )

  • 图5 邻域法检验的风场个例:(a)理想化的观测风场;(b)假定的预报风场;(c)检验评分与预报阵风位移的关系(引自Skok and Hladnik(2018))

  • Fig.5 A case study demonstrating wind field verification using the neighborhood method: (a) idealized observed wind field; (b) hypothetical forecasted wind field; (c) relationship between verification scores and forecast gust displacement (cited from Skok and Hladnik (2018) )

  • T=a bc d
    (13)
  • 邻域窗内采用误差补偿的列联表T事件分类有:

  • T=a+min(b,c)b-min(b,c)c-min(b,c)d+min(b,c)
    (14)
  • 式中:a为事件命中的次数; b为空报的次数; c为漏报的次数; d为正确预报未超过阈值事件的次数。式(14)是以定义的邻域窗为单位进行事件划分,min()为取最小值函数。

  • Stein and Stoop(2019)采用误差补偿列联表,分析了高分辨率数值模式在不同邻域半径下的降水预报评分表现; 结果发现,当降水阈值高于5.0 mm·d-1时,命中率下降较快(图6a),而空报率随着阈值的增加而增加(图6b)。由于邻域内使用了误差补偿策略,邻域越大,这些评分的预报表现越好(图6c)。分析表明,相对于单纯的FSS评分,误差补偿列联表能够进一步抑制高分辨率中尺度模式中因“双重惩罚”效应带来的评分偏低的现象,与主观判断结果更为一致。同时,预报和观测之间的位置偏移距离作为邻域大小的函数,可以从误差补偿策略的评分中反映出来,且误差补偿策略的列联表不会改变预报偏差,对整个研究区域来说,预报偏差也不取决于邻域的大小。

  • 3.2 邻域集合概率

  • NEP在CAM的集合降水预报(Johnson and Wang,2012; Yussouf et al.,2013; Romine et al.,2014; Schwartz,2017)和天气雷达反射率模拟预报(Snook et al.,2012; Hitchcock et al.,2016)等方面有着广泛的应用。这些研究主要是针对小时雨量小于15 mm的降水邻域概率预报或者雷达反射率因子≥19 dBZ的降水概率计算,很少有工作针对小时较大雨强的降水计算NEP并进行检验(Yussouf et al.,2016)。可能的原因是,NEP对稀有事件通常不具有良好的可靠性或分辨率。因此,NEP的这种属性导致其更适合应用于气候概率较大的事件预报。与NEP不同,尽管NMEP可以在空间上进行平滑,但NMEP的生成本身不包含空间平滑,因此它对稀有事件具有更好的分辨能力(Barthold et al.,2015; Schwartz et al.,2015)。

  • 国内学者对集成邻域概率也有广泛的研究。刘雪晴等(2020)利用区域集合预报系统降水预报资料进行邻域集合概率方法试验,并针对邻域概率法的等权重和邻域尺度问题,设计了邻域格点权重修正方案,同时评估了试验方案的预报效果。罗聪等(2021)利用邻域最优概率方法对华南区域GRAPES快速更新循环同化预报系统的24 h预报进行逐时降水订正和检验评估; 结果表明,邻域最优概率方法能有效地提升降水客观预报能力。刘莹等(2022)利用区域集合模式的小时降水产品开展了集合预报邻域法订正试验,认为邻域集合最大概率预报能充分显示大范围的降水中心,可以为预报员提供多视角预报参考。针对极端降水天气,朱科锋等(2022)分析了河南“21.7”极端暴雨天气过程的对流可分辨尺度集合预报的降水概率预报表现(图7); 结果表明,与传统的降水概率相比,邻域集合概率显著提升了日降雨和小时降雨的概率预报技巧。这些研究还表明,NEP和NMEP都可以产生邻域集合预报概率,但两种概率预报产品不能相互替代,一种集合概率预报的阈值也并不适用于另一种产品。Schwartz and Sobash(2017)的研究表明,在暴雨以下量级预报中,NEP有很好的应用效果,但就暴雨及以上量级的预报来说,NMEP对极值预报的精度损失更小,可以更有效预报极端天气事件,从而有效避免漏报。

  • 图6 误差空间补偿应用示例(不同邻域半径的误差补偿列联表的预报评分):(a)命中率;(b)空报率;(c)FSS评分(引自Stein and Stoop(2019))

  • Fig.6 Example application of neighborhood error compensation in verification scores, showing results for varying neighborhood radii: (a) probability of detection; (b) false alarm ratio; (c) FSS score (cited from Stein and Stoop (2019) )

  • 4 讨论

  • 检验评估是邻域法在气象预报中最主要的应用之一。基于邻域窗的预报和观测值可以在定义阈值的基础上,构建“一对多”和“多对多”等多种不同的匹配方案,也可以在邻域窗内将强度超过某个阈值的比例定义为其在邻域窗内发生的概率,并检验邻域概率的预报表现。本文给出的检验评分只是邻域法检验中最常用的一部分,在邻域法检验中想要深入获取所有邻域半径下产品的全部预报表现十分困难,因此,对邻域法检验的使用者来说,首先要明确预报需求,什么是一个好的预报,然后利用好的检验模型,通过计算评分指数对预报做出判断。大多数情况下,现有的方法能够为使用者提供在不同空间尺度下的预报判断信息,如果没有合适的判别模型,也能够方便地利用邻域法特点、预报和观测的强度分布等构造一个或多个评分指数,譬如ROC(receiver operating characteristic)曲线图(Theis et al.,2005)以及与尺度相关的方差等,来评估不同尺度下预报对观测的重现能力。

  • 图7 2021年7月19日12时到20日12时河南“21.7”极端暴雨天气过程的邻域集合降水概率预报(单位:%):(a)24 h累积降水量≥50 mm;(b)24 h累积降水量≥150 mm(引自朱科锋等(2022))

  • Fig.7 Neighborhood ensemble precipitation probability forecast for the “21.7” extreme rainstorm weather process in Henan from 1200 UTC on July 19 to 1200 UTC on July 20, 2021 (units: %) : (a) 24 h cumulative precipitation≥50 mm; (b) 24 h cumulative precipitation≥150 mm (cited from Zhu et al. (2022) )

  • 邻域概率在CAM后处理中有着广泛的应用。确定性预报的邻域概率一般通过给定的阈值条件,在邻域窗内对预报事件进行标记,计算预报事件的发生概率,集合预报概率NEP和邻域最大集合预报概率NMEP是最主要的两种邻域集合预报概率方法。这两种方法产生邻域集合预报概率的差别在于:大多数CAM的集合概率预报偏低(Romine et al,2014),由于NMEP在邻域窗内取最大值,所以它能更好地表现预报的极端信息,提高极端事件的命中率。而对流极端天气通常发生在较小的空间和时间尺度上,NEP算法中固有的平滑使得到的概率接近于零,几乎没有什么意义(Schwartz and Sobash,2017)。但NMEP会对小概率事件造成过度的空报,增大极端事件的预测频率,在使用时一方面可以通过增加成员的数量来进一步增加可靠性,另一方面可以增大阈值参数或利用客观算法进行抑制。结合实际业务,也可对NEP和NMEP进行综合诊断分析,原因在于:NEP对NMEP概率有较好的抑制作用,只有当事件发生的概率较低时,NMEP才与NEP接近。此外,NEP和NMEP方法获得的预报结果在统计学检验上存在显著差异,其预报表现不仅与降水性质、概率阈值相关,而且与预报区域降水的气候态以及评价标准等紧密相关。在实际应用中,一方面应明确描述其使用的邻域方法和事件定义,另一方面应对不同预报产品和研究区域等影响降水预报准确率的因素进行深入分析,并采用合适的邻域半径和集成方案来提高降水预报表现。

  • 需要说明的是,本文以降水、风速为主归纳了邻域法的重点应用,但事实上,在气象要素预报的各种情景下,邻域法都有较好的适用性,譬如灾害性天气上升气流螺旋度(Sobash et al.,2011,2016; Clark et al.,2013; Schwartz et al.,2015)、基于物理量集成的高阈值垂直涡度场(Zhang et al.,2016)、中尺度对流系统(Stratman et al.,2013;Snook et al.,2015,2016)、闪电(Lynn et al.,2015)、低层垂直涡度(Yussouf et al.,2013; Wheatley et al.,2015)等在灾害性天气中的要素或物理量的概率预报。

  • 5 结论和展望

  • 本文系统回顾了邻域法检验的技术框架、常用的检验评分以及两种集成邻域概率的计算流程,给出了邻域法检验和概率的典型应用,详细讨论了邻域法检验和邻域概率应用中的注意事项,得到如下主要结论:

  • 1)邻域法可以在时间或空间的多个尺度比较预报产品的性能。对于不同分辨率的网格预报,通过尺度变换,可获得模式在多大尺度上的预报更加准确等重要信息。

  • 2)邻域法检验对预报和观测的数据处理具有通用性,传统的二分法检验或概率预报在邻域法检验中都是适用的,其独特的优势在于,相对传统的二分法检验,它增加了基于尺度变化获得的一些诊断量。FSS评分在邻域法检验中得到了最为广泛的应用。

  • 3)基于邻域法可以计算确定性预报的邻域概率,但更广泛的应用是利用集合预报多成员或多个确定性预报获得NEP或者NMEP两种邻域集合概率。

  • 4)NEP和NMEP都可以提高降水的预报评分,但两者的应用范围显著不同,NEP更适合于大尺度、系统性降水预报,对于对流性、极端性降水来说,NMEP的应用效果更好。

  • 基于目前的研究进展,本文认为以下方面值得进一步研究:

  • 1)邻域半径的合理选择问题。合适的邻域半径是邻域法检验和邻域集合概率的核心问题。在实际应用中,一般通过“穷举法”计算不同邻域半径的预报评分来获得最优邻域尺度。但是在不断增加的空间尺度上通过评分的大小来捕获产品的预报表现,可能会受到小的、低可预测性事件的影响。当下垫面复杂且研究区域较大时,用单一的邻域尺度来量化产品的预报表现或基于此构建邻域概率,可能会得不到理想的结果。关注研究区域和预报产品的特性,构建差异性或自适应最优邻域半径,是一个重要的研究方向。

  • 2)邻域维度拓展问题。已有邻域法检验的研究主要考虑了二维空间维度,集成预报的邻域概率通过增加成员,形成三维空间,但目前时间维邻域法的应用研究还不多见,为了量化NWP的时间偏差,一般采用时域追踪方法。对高时间分辨率的NWP来说,将邻域法拓展到时间维对邻域法检验及其概率具有重要意义。此外,就一些灾害性天气事件而言,拓展阈值条件也是邻域法检验值得探索的一个方面。

  • 3)邻域法的检验和概率预报的评价指标问题。FSS评分是邻域法检验最主要的评价指标,其他传统的二分法评价指标也可以以邻域窗为基础展开。基于NEP和NMEP获得的概率预报,一般通过已有的概率预报检验方法进行评价。事实上,NEP有空间平滑的作用,获得的概率偏小,而NMEP则突出了局地和单个成员的贡献,对极端事件预报能力较好,但概率偏大。如何构建有效的评价指标使其结果更好地表现邻近格点的补偿作用,或者结合目标检验等一些空间检验方法解析邻域概率预报的属性特征,提高预报技巧,对邻域概率和邻域检验都是重要的研究课题。

  • 4)邻域法与人工智能的结合。邻域法通过相邻格点信息的相互补偿,挖掘预报产品中更多的有用信息,但这些信息的获得,目前为止仍然是一种基于历史资料的统计结果,在邻域窗设定、集合邻域半径优选、不同集合成员最优邻域窗的自洽和耦合等方面仍然存在很多未知、可探索的领域。将邻域法与人工智能算法结合起来,形成正向自反馈,在天气预报中的研究应用前景广阔。

  • 致谢:两位匿名审稿人提出了宝贵的修改意见和建议。谨致谢忱!

  • 参考文献

  • 参考文献

  • 地址:江苏南京宁六路219号南京信息工程大学    邮编:210044

    联系电话:025-58731158    E-mail:xbbjb@nuist.edu.cn    QQ交流群号:344646895

    大气科学学报 ® 2025 版权所有  技术支持:北京勤云科技发展有限公司