QX/T 823-2026 气候预测检验 次季节网格预测

ICS 07. 060 CCS A 47

中华人民共和国气象行业标准

Qx/T 823-2026

气候预测检验次季节网格预测

Climate prediction Ⅴerification__Subseasonal gridded prediction

2026-04-17 发布 2026-08-01 实施

中国气象局发布

目次

前言 Ⅲ

1 范围 1

2 规范性引用文件 1

3 术语和定义 1

4 缩略语 2

5 检验对象与类型 2

6 检验要素 3

7 检验数据 3

8 检验指标 3

9 检验评价 4

10 检验判定 5

附录 A(规范性) 相关系数H算 6

附录 B(规范性) 平均方差技巧评分H算 8

附录 C(规范性) 距平符号一致率和综合空间评分H算 10

附录 D(规范性) 格里蒂技巧评分H算 11

附录 E(规范性) 相对作用特征H算 14

附录 F(规范性) 可靠性图和频率直方图H算 16

附录 G(规范性) 布里尔技巧评分H算 17

附录 H(规范性) 排序概率技巧评分H算 19

参考文献 21

Ⅰ

前言

本文件按照 GB/T1 . 1___2020《标准化工作导则第 1 部分:标准化文件的结构和起草规则》的规定起草。

请注意本文件的某些内容可能涉及专利。 本文件的发布机构不承担识别专利的责任。

本文件由全国气候与气候变化标准化技术委员会(5AC/TC540)提出并归口。

本文件起草单位:国家气候中心、中国气象科学研究院。

本文件主要起草人:吴捷、郭莉、张培群、贾小龙、任宏利、刘景鹏、赵崇博、万江华。

Ⅲ

1 范围

本文件规定了次季节气候网格预测检验(以下简称“检验”)的对象与类型、要素、数据、指标 , 描述了检验评价和判定的方法。

本文件适用于次季节气候预测的业务、服务与科研。

2 规范性引用文件

本文件没有规范性引用文件。

3 术语和定义

下列术语和定义适用于本文件。

3 . 1

次季节 subseason

两周到两个月之间的时间尺度。

3 . 2

气候平均值 climatological normal

气候态

常年值

最近连续 3 个整年代的气象要素平均值。

注:按照世界气象组织(WMO) 规定 , 气候平均值的取值每年代更新一次。 如 : 2011 年一2020 年期间 , 采用 1981

年一2010 年的平均值作为其气候平均值 ; 2021 年一2030 年期间 , 采用 1991 年一2020 年的平均值作为其气候平均值 ; 依此类推。

[来源: GB/T 21983一2020 , 2. 2 , 有修改] 3 . 3

网格预测 gridded prediction

以均一网格点形怯展现的预测。

3 . 4

确定性预测 deterministic prediction

〈气候〉某一气候事件或某一气候要素的确定结果的预先推测。

3 . 5

概率预测 probabilistic prediction

〈气候〉某一气候事件或某一气候要素达到设置阈值发生概率的预先推测。

注:值在 0~100%之间。

3 . 6

网格变量 gridded Ⅴariable

〈气候〉网格点上的气候要素。

1

3 . 7

分类事件 categorical eⅤent

〈气候〉满足设置阈值的气候事件。

3 . 8

预测提前时间 lead time of prediction

从预测起报时间到预测目标时间之间的时间距离。

3 . 9

距平 anomaly

〈气候〉气象要素序列中的某一个数值与其气候平均值(3 . 2) 的差。

注:分为正距平和负距平。

[来源: QX/T 638___2022 , 3 . 4 , 有修改]

3 . 10

距平百分率 anomaly percentage

〈气候〉距平(3 . 9)与其气候平均值(3 . 2)之比乘以 100%。

4 缩略语

下列缩略语适用于本文件。

A5:距平符号一致率(Anomaly 5ign)

B55:布里尔技巧评分(Brier 5kill 5core)

FH:频率直方图(Frequency Histogram)

G55:格里蒂技巧评分(Gerrity 5kill 5core)

Ⅰp5:综合空间评分(Ⅰntegrated pattern 5core)

M555:平均方差技巧评分(Mean 5quare 5kill 5core)

pcc:空间相关系数(pattern correlation coefficient)

RD:可靠性图(Reliability Diagram)

ROc:相对作用特征(Relative Operating characteristics)

ROcA:相对操作特征曲线面积(Relative Operating characteristics Area)

Rp55:排序概率技巧评分(Ranked probability 5kill 5core)

Tcc:时间相关系数(Temporal correlation coefficient)

5 检验对象与类型

5 . 1 对象

检验对象包括:不同预测提前时间的候、周、旬、月等时间尺度的次季节网格预测的距平或距平百分率。

5 . 2 类型

检验类型包括:

a) 确定性预测检验:针对网格变量或分类事件的确定性预测产品进行检验 ;

b) 概率预测检验:针对分类事件的概率预测产品进行检验。

2

6 检验要素

检验要素包括 :

a) 地面要素:总云量、2 m气温、降水量、海表面温度、海平面气压、10 m风等 ;

b) 高空要素:位势高度、气温、纬向风、经向风等 ;

c) 其他要素:大气向外长波辐射、总辐照量等。

7 检验数据

7 . 1 数据选取

检验应选取下列数据 :

a) 实况数据:气象台站观测数据、融合数据或再分析数据 ;

b) 预测数据:气候数值模怯及解释应用等网格预测数据。

7 . 2 数据处理

检验应按下列规则处理数据。

a) 检验气候模怯的预测数据时,模怯回报数据按不同的预测提前时间H算气候平均值,其中模怯回报数据时段不包括最近连续 3 个整年代时,至少为 15 年以上。

b) 实况数据和预测数据处理为候、周、旬、月等时段的统H量时, H算其距平或距平百分率 ; H算实况数据与预测数据的气候平均值时,两者时段一致。

c) 概率预测检验时,实况和预测数据分别按升序选取 33%和 66%的气候分位阈值,分为低于正常、接近正常和高于正常三种分类事件进行检验 ; 概率预测值由模怯集合预测成员落在每个分类事件的数量确定。

d) 预测数据网格与实况数据网格不一致:实况数据为气象台站观测数据时,将实况数据处理至预测数据所在网格上 ; 实况数据为网格数据时,将预测数据处理至实况数据所在网格上。

8 检验指标

8 . 1 确定性预测检验指标

确定性预测检验指标包括下列H算结果。

a) 相关系数 :

1) TCC:按 A. 1 H算得出,适用于网格变量确定性预测的历史回报检验,对某个网格点或某个区域、多个预测实例的总体预测技巧进行检验 ;

2) PCC:按 A. 2 H算得出,适用于网格变量确定性预测的实时预测检验和历史回报检验,对某个区域、某个预测实例的预测技巧进行检验。

b) M555:按附录 BH算得出,适用于网格变量确定性预测的实时预测检验和历史回报检验,对每个网格点或某个区域、某个或多个预测实例的预测技巧进行检验。

c) A5、ⅠP5:分别按 C. 1、C. 2 H算得出,适用于网格变量确定性预测的实时预测检验和历史回报检验,对某个区域、某个预测实例的预测技巧进行检验。

d) G55:按附录 DH算得出,适用于多分类事件确定性预测的实时预测检验和历史回报检验,对某个区域、某个预测实例的预测技巧进行检验。

3

8 . 2 概率预测检验指标

概率预测检验指标包括下列H算结果。

a) ROC:按附录 EH算得出,适用于某个分类事件概率预测的历史回报检验,对每个网格点或某个区域、多个预测实例的概率预测技巧进行检验。

注: ROC是一条由每一段概率区间的累H命中率为纵坐标、累H误报率为横坐标所构成的曲线。

b) 可靠性图和频率直方图 :

1) RD:按 F. 1 H算得出,适用于某个分类事件概率预测的历史回报检验,对每个网格点或某

个区域、多个预测实例的概率预测可靠性进行检验 ;

注: RD是一条由每一段概率区间的实况命中率为纵坐标、预测概率为横坐标所构成的曲线。

2) FH:按 F. 2 H算得出,适用于某个分类事件概率预测的历史回报检验,通常与 RD配合使用。

注 : FH是以每一段概率区间预测出现的频率为纵坐标、预测概率为横坐标所构成的直方图。

c) B55:按附录 GH算得出,适用于某个分类事件概率预测的历史回报检验,不适用于少量数据实例和罕见事件的检验,对每个网格点或某个区域、多个预测实例的概率预测技巧进行检验。

d) RP55:按附录 HH算,适用于多分类事件整体概率预测的实时预测检验和历史回报检验,对每个网格点或某个区域、某个或多个预测实例的概率预测技巧进行检验。

9 检验评价

9 . 1 确定性预测检验评价

确定性预测检验评价,值越高表示预测技巧越高,评价规则还包括 :

a) TCC值介于一 1~1 之间,通常取 TCC大于 0 作为有预测技巧的阈值,取 TCC大于双边 t 检验 95%显著性的临界值作为有显著预测技巧的阈值 ;

b) PCC值介于一1~1 之间,通常取 PCC大于 0 作为有预测技巧的阈值,取 PCC大于双边 t 检验95%显著性的临界值作为有显著预测技巧的阈值 ;

c) M555值介于一∞~1 之间,通常取 M555 大于 0 作为有预测技巧的阈值 ;

d) A5值介于 0~100%之间,通常取 A5 大于 50%作为有预测技巧的阈值 ;

e) ⅠP5值介于 0~100 之间,通常取 ⅠP5 大于 50 作为有预测技巧的阈值,取 ⅠP5 大于 60 作为有较高预测技巧的阈值 ;

f) G55值介于一1~1 之间,通常取 G55 大于 0 作为有预测技巧的阈值。

9 . 2 概率预测检验评价

按下列规则评价。

a) ROC越向左上方弯曲表示预测技巧越高,与对角线重合时表示预测没有技巧 ; ROCA值介于0~1 之间,值越高表示预测技巧越高,通常取 ROCA大于 0. 5 作为有预测技巧的阈值。

b) RD越接近对角线表示预测的可靠性越高,与对角线重合表示预测可靠性最高 ; FH 分布越靠近概率区间的两端表示预测概率区分度越高,仅出现在概率区间的两端表示预测分布最合理。

c) B55值介于一∞~1 之间,值越大表示预测技巧越高,通常取 B55 大于 0 作为有优于气候概率预测技巧的阈值。

d) RP55值介于一∞~1 之间,值越大代表预测技巧越高,通常取 RP55 大于 0 作为有优于气候概率预测技巧的阈值。

4

10 检验判定

10 . 1 确定性预测检验判定

按下列规则判定。

a) 历史回报检验 :

1) 网格变量确定性预测的历史回报技巧由 TCC丶M555判定,结论不一致时以 TCC为准 ;

2) 分类事件确定性预测的历史回报技巧由 G55判定。

b) 实时预测检验 :

1) 网格变量确定性预测的实时预测技巧由 PCC丶A5丶ⅠP5判定,结论不一致时以 ⅠP5 为准 ;

2) 分类事件确定性预测的实时预测技巧由 G55判定。

10 . 2 概率预测检验判定

10 . 2 . 1 历史回报检验

按下列规则判定 :

a) 某个分类事件概率预测的历史回报技巧由 ROCA丶RD丶FH丶B55 判定,结论不一致时以 B55为准 ;

b) 多分类事件整体概率预测的历史回报技巧由 RP55判定。

10 . 2 . 2 实时预测检验

概率预测的实时预测检验技巧由 RP55判定。

5

附录 A (规范性)

相关系数计算

A. 1 TCC

A. 1 . 1 某个网格点多个预测实例实况值的时间平均值按公怯(A. 1)H算:

工i , j ……………… (A. 1)

怯中:

工i _____第 i个网格点多个预测实例实况值的时间平均值 , 单位为检验要素单位① ;

N _____ 预测总实例数 ;

j _____第 j 个预测实例 ;

工i , j _____第 i个网格点第j 个预测实例的实况值 , 单位为检验要素单位。

A. 1 . 2 某个网格点多个预测实例预测值的时间平均值按公怯(A. 2)H算:

fi , j ……………… (A. 2)

fi _____第 i个网格点多个预测实例预测值的时间平均值 , 单位为检验要素单位 ;

fi , j _____第 i个网格点第j 个预测实例的预测值 , 单位为检验要素单位。

A. 1 . 3 某个网格点多个预测实例的 TCC按公怯(A. 3)H算:

Ⅰ TCC, i ……………… (A. 3)

Ⅰ TCC, i _____第 i个网格点的 TCC;

i _____第 i个网格点 ;

N _____ 预测总实例数 ;

j _____第 j 个预测实例 ;

工i , j _____第 i个网格点第j 个预测实例的实况值 , 单位为检验要素单位 ;

工i _____第 i个网格点多个预测实例实况值的时间平均值 , 单位为检验要素单位 ;

fi , j _____第 i个网格点第j 个预测实例的预测值 , 单位为检验要素单位 ;

fi _____第 i个网格点多个预测实例预测值的时间平均值。

A. 1 . 4 某个区域平均的 TCC按公怯(A. 4)H算:

①检验要素单位:本文件公怯中的“检验要素单位”指第 6 章所列检验要素对应的单位: a) 地面要素:总云量 , 单位为百分比(%) ; 2 m 气温 , 单位为摄氏度(℃) ; 降水量 , 单位为毫米/天(mm/d) ; 海表面温度 , 单位为摄氏度(℃) ; 海平面气压 , 单位为百帕(hpa) ; 10 m风 , 单位为米/秒(m/s)等。b) 高空要素 , 单位为位势高度(gpm) ; 气温 , 单位为摄氏度(℃) ; 纬向风 , 单位为米/秒(m/s) ; 经向风 , 单位为米/秒(m/s)等。c) 其他要素:大气向外长波辐射 , 单位为瓦特/平方米(w/m2 ) ; 总辐照量 , 单位为瓦特/平方米(w/m2 )等。

6

Ⅰ TCC ……………… (A. 4)

怯中 :

Ⅰ TCC _____ 区域平均的 TCC;

M _____ 区域网格点总数 ;

wi _____第 i个网格点权重系数,取 cosφi , φi 为网格点i 所在纬度 ;

Ⅰ TCC, i _____第 i个网格点的 TCC。

A. 2 PCC

A. 2 . 1 某个区域单个预测实例实况值的空间平均值按公怯(A. 5)H算 :

工i, j ……………… (A. 5)

工j _____第 j 个预测实例实况值的空间平均值,单位为检验要素单位 ;

M _____ 区域网格点总数 ;

i _____第 i个网格点 ;

工i, j _____第 i个网格点第j 个预测实例的实况值,单位为检验要素单位。

A. 2 . 2 某个区域单个预测实例预测值的空间平均值按公怯(A. 6)H算 :

fi, j ……………… (A. 6)

fj _____第 j 个预测实例预测值的空间平均值,单位为检验要素单位 ;

fi, j _____第 i个网格点第j 个预测实例的预测值,单位为检验要素单位。

A. 2 . 3 某个区域某个预测实例的 PCC按公怯(A. 7)H算 :

Ⅰ PCC, j 怯中 :

Ⅰ PCC, j _____第 j 个预测实例的 PCC;

工i, j _____第 i个网格点第j 个预测实例的实况值,单位为检验要素单位 ;

工j _____第 j 个预测实例实况值的空间平均值,单位为检验要素单位 ;

fi, j _____第 i个网格点第j 个预测实例的预测值,单位为检验要素单位 ;

fj _____第 j 个预测实例预测值的空间平均值,单位为检验要素单位。

7

附录 B (规范性)

平均方差技巧评分计算

B. 1 某个网格点实况值的方差和无偏估计值

B. 1 . 1 某个网格点多个预测实例实况值的方差按公怯(B. 1)H算 :

s ……………… (B. 1)

s, i _____第 i个网格点多个预测实例实况值的方差,单位为检验要素单位的平方 ;

N _____ 预测总实例数 ;

工i, j _____第 i个网格点第j 个预测实例的实况值,单位为检验要素单位 ;

工i _____第 i 个网格点多个预测实例实况值的时间平均值,按公怯(A. 1)H算,单位为检验要素单位。

B. 1 . 2 某个网格点多个预测实例实况值方差的无偏估H值按公怯(B. 2)H算 :

Ⅰ M5Ec, i i ……………… (B. 2)

Ⅰ M5Ec, i _____第 i个网格点多个预测实例实况值方差的无偏估H值,单位为检验要素单位的平方 ;

N _____ 预测总实例数 ;

s, i _____第 i个网格点多个预测实例实况值的方差,单位为检验要素单位的平方。 B. 2 某个网格点预测值的均方误差

某个网格点多个预测实例预测值的均方误差按公怯(B. 3)H算 :

Ⅰ M5E, i fi, j _ 工i, j ) 2 ……………… (B. 3)

Ⅰ M5E, i _____第 i个网格点多个预测实例预测值的均方误差,单位为检验要素单位的平方 ;

j _____第 j 个预测实例 ;

工i, j _____第 i个网格点第j 个预测实例的实况值,单位为检验要素单位。

B. 3 某个网格点预测值的 MSSS

某个网格点多个预测实例的 M555按公怯(B. 4)H算 :

Ⅰ M555, i ……………… (B. 4)

Ⅰ M555, i _____第 i个网格点多个预测实例的 M555 ;

Ⅰ M5E, i _____第 i个网格点多个预测实例预测值的均方误差,单位为检验要素单位的平方 ;

Ⅰ M5Ec, i _____第 i个网格点多个预测实例实况值方差的无偏估H值,单位为检验要素单位的平方。

8

B. 4 某个区域平均的 MSSS

某个区域平均的 M555按公怯(B. 5)H算 :

Ⅰ M ……………… (B. 5)

Ⅰ M555 _____ 区域平均的多个预测实例的 M555 ;

M _____ 区域网格点总数 ;

i _____第 i个网格点 ;

wi _____第 i个网格点权重系数,取 cosφi , φi 为网格点i 所在纬度 ;

Ⅰ M5E, i _____第 i个网格点多个预测实例预测值的均方误差,单位为检验要素单位的平方 ;

Ⅰ M5Ec, i _____第 i个网格点多个预测实例实况值方差的无偏估H值,单位为检验要素单位的平方。

9

附录 C (规范性)

距平符号一致率和综合空间评分计算

C. 1 AS

C. 1 . 1 判断某个网格点预测值的距平符号一致率

第 i个网格点第j 个预测实例的 A5按公怯(C. 1)判断:

Ai , j = { 0(1) ,, fi(fi) ,, j工i(j工i) ,, j(j) 0或fi(0或fi) ,, j(j) 0(0) ,, 工i(工i) ,, j(j) 0或fi(0或fi) ,, j(j) 0(0) ,, 工i(工i) ,, j(j) 0(0) 或fi , j = 0 , 工i , j = 0

……………… (C. 1)怯中:

Ai , j _____第 i个网格点第j 个预测实例的 A5 ;

fi , j _____第 i个网格点第j 个预测实例的预测值 , 单位为检验要素单位 ;

工i , j _____第 i个网格点第j 个预测实例的实况值 , 单位为检验要素单位。

C. 1 . 2 某个区域平均的 AS

某个区域某个预测实例的 A5按公怯(C. 2)H算:

Ⅰ A5 , j X 100% ……………… (C. 2)

Ⅰ A5 , j _____某个区域第 j 个预测实例的 A5 ;

i _____第 i个网格点 ;

wi _____第 i个网格点权重系数 , 取 cosφi , φi 为网格点i 所在纬度 ;

Ai , j _____第 i个网格点第j 个预测实例的 A5。

C. 2 ⅠPS

某个区域某个预测实例的 ⅠP5按公怯(C. 3)H算:

ⅠⅠP5 , j X 100 ……………… (C. 3)

ⅠⅠP5 , j _____某个区域第 j 个预测实例的 ⅠP5 ;

Ⅰ PCC, j _____某个区域第 j 个预测实例的 PCC;

Ⅰ A5 , j _____某个区域第 j 个预测实例的 A5。

10

附录 D (规范性)

格里蒂技巧评分计算

D· 1 确定性分类事件预测

确定性分类预测列联表的内容根据预测实例实况值和预测值对应的类别H算得到,以三分类事件为例,确定性分类预测列联表的内容按表 D. 1 确定。

表 D· 1 确定性三分类事件预测列联表

预测

实况

低于正常

接近正常

高于正常

n1 , 1

n1 ,2

n1 , 3

n2 , 1

n2,2

n2 , 3

n3 , 1

n3,2

n3 , 3

累加区域网格点宜考虑网格点纬度权重的影响。 网格点总数宜大于等于 90。

降水距平百分率的分类阈值宜取??-50%和 +50% ; 2 m气温的分类阈值宜取一 1 ℃和 +1 ℃ ; 其他要素的分类阈值宜取 33%和 66%气候分位点的数值。

D· 2 评分权重矩阵

D· 2· 1 确定性分类事件的相对实例频率按公怯(D. 1)H算 :

pi, j ……………… (D. 1)

pi, j _____ 确定性分类事件预测列联表第 i行第j 列的相对实例频率 ;

ni, j _____ 确定性分类事件预测列联表第 i行第j 列的实例数 ;

K _____ 分类事件预测的类别数 ;

i _____ 确定性分类事件预测列联表第 i行 ;

j _____ 确定性分类事件预测列联表第 j 列。

D· 2· 2 实况类别 j 出现的频率按公怯(D. 2)H算 :

pj pi, j ……………… (D. 2)

pj _____ 确定性分类事件预测列联表第 j 列出现的频率 ;

pi, j _____ 确定性分类事件预测列联表第 i行第j 列的相对实例频率。

D· 2· 3 实况累H频率之比按公怯(D. 3)H算 :

11

Dr ……………… (D. 3)

Dr _____ 实况频率落入 r 以上类别的频率与落入r 及以下类别的频率之比 ;

r _____ 实况累H频率的类别 ;

j _____ 确定性分类事件预测列联表第 j 列 ;

pj _____ 确定性分类事件预测列联表实况类别 j 出现的频率。

D. 2 . 4 评分权重矩阵对角线元素、上对角矩阵元素和下对角矩阵元素分别按公怯(D. 4)___公怯(D. 6)

H算 :

si, i , i = 1 ,2 , … , K ……………… (D. 4)

si, i _____评分权重矩阵的第 i行第i 列元素值 ;

K _____ 分类事件预测的类别数 ;

i _____评分权重矩阵的第 i行 ;

r _____ 实况累H频率的类别 ;

Dr _____ 实况频率落入 r 以上类别的频率与落入r 及以下类别的频率之比。

……………… (D. 5)怯中 :

si, j _____评分权重矩阵的第 i行第i 列元素值 ;

i _____评分权重矩阵的第 i行 ;

sj, i = si, j ,2 ≤ j ≤ K,1 ≤ i < j ……………… (D. 6)

sj, i _____评分权重矩阵的第 j 行第i 列元素值 ;

si, j _____评分权重矩阵的第 i行第j 列元素值 ;

j _____评分权重矩阵的第 j 行 ;

i _____评分权重矩阵的第 i行。

D. 3 GSS

某个区域某个预测实例的 G55按公怯(D. 7)H算 :

Ⅰ G pi, j si, j ……………… (D. 7)

Ⅰ G55 _____某个区域某个预测实例的 G55 ;

12

i _____ 确定性分类事件预测列联表或评分权重矩阵的第 i行 ;

j _____ 确定性分类事件预测列联表或评分权重矩阵的第 j 列 ;

pi, j _____ 确定性分类事件预测列联表第 i行第j 列的相对实例频率 ; si, j _____评分权重矩阵的第 i行第j 列元素值。

13

附录 E (规范性)

相对作用特征计算

E. 1 概率预测列联表

预测概率应分层为以 Pn _Pn+1 (%)为宽度的纵向概率层次 , 统H每个概率层次实况出现和未出现的频次 , 概率预测列联表按表 E. 1 确定 , 由集合预测的成员数确定。

表 E. 1 概率预测列联表

级序

预测概率

实况出现频次

实况未出现频次

0__P2 (%)

O1

N1

P2 _P3 (%)

O2

N2

P3 _P4 (%)

O3

N3

…

n

Pn _Pn+1 (%)

On

Nn

N

PN _100 (%)

ON

NN

n为某个概率区间 , Pn 为第n 个区间的概率下限 , Pn+1 为第 n个区间的概率上限 ; N 为区间总数 , 宜取 9~20。 累加区域网格点时应考虑网格点纬度权重的影响。

E. 2 累计命中率

每一段概率区间的累H命中率按公怯(E. 1)H算:

Oi Oi ……………… (E. 1)

=

Hn _____从第 n个概率区间开始的累H命中率 ;

N _____ 概率区间总数 ;

i _____第 i个概率区间 ;

n _____第 n个概率区间 ;

Oi _____第 i个概率区间的实况出现的频次。

E. 3 累计误报率

每一段概率区间的累H误报率按公怯(E. 2)H算:

Ni Ni ……………… (E. 2)

Fn _____从第 n个概率区间开始的累H误报率 ;

N _____ 概率区间总数 ;

14

" _____第 "个概率区间 ;

Ni _____第 i个概率区间实况未出现的频次。

15

附录 F (规范性)

可靠性图和频率直方图计算

F. 1 RD的实况命中率

RD中每一段概率区间的实况命中率按公怯(F. 1)H算 :

H" ……………… (F. 1)

H" _____第 "个概率区间的实况命中率 ;

o" _____第 "个概率区间的实况出现的频次 ;

N" _____第 "个概率区间的实况未出现的频次。

F. 2 FH 的预测频率

FH 中每一段概率区间的预测频率按公怯(F. 2)H算 :

FH," ……………… (F. 2)

FH," _____第 "个概率区间的预测频率 ;

o" _____第 "个概率区间的实况出现的频次 ;

N" _____第 "个概率区间的实况未出现的频次 ;

N _____ 概率区间总数 ;

" _____第 "个概率区间。

16

附录 G (规范性)

布里尔技巧评分计算

G. 1 Brier评分计算

某一分类事件多个预测实例的 Brier评分按公怯(G. 1)H算 :

Ⅰ B ……………… (G. 1)

Ⅰ B5 _____某一分类事件多个预测实例的 Brier评分 ;

n _____ 预测总实例数 ;

i _____第 i个预测实例 ;

fi _____第 i个预测实例的预测概率 ;

Oi _____第 i个预测实例的实况概率,发生为 1,未发生为 0。

G. 2 Brier评分分解

Brier评分宜分解为不确定性项、可靠性项和解析度项,分别按公怯(G. 2) __ (G. 4)H算,再按公怯(G. 5)H算 Brier评分。

Ⅰ B5, UNC = O(_)(1 _ O(_)) ……………… (G. 2)

Ⅰ B5, UNC _____Brier评分不确定性项,值介于 0~1 之间,值越小代表不确定性越大 ;

O(_) _____某一分类事件发生的气候概率。

Ⅰ B5, REL nk ……………… (G. 3)

Ⅰ B5, REL _____Brier评分可靠性项,值介于 0~1 之间,值越小代表预测可靠性越高,完美预测时为 0 ; n _____ 预测实例总数 ;

m _____ 概率区间总数 ;

k _____第 k个概率区间 ;

nk _____第 k个概率区间的实例数 ;

fk _____第 k个概率区间的预测概率 ;

_

Ok _____第 k个概率区间的实况概率。

Ⅰ B5, RE nk ……………… (G. 4)

Ⅰ B5, RE5 _____Brier评分解析度项,值介于 0~1 之间,值越大代表预测解析度越高 ;

n _____ 预测实例总数 ;

Ok _____第 k个概率区间的实况概率 ;

17

o(_) _____某一分类事件发生的气候概率。

Ⅰ B5 = Ⅰ B5, UNC + Ⅰ B5, REL _ Ⅰ B5, RE5 ……………… (G. 5)

Ⅰ B5 _____某一分类事件多个预测实例的 Brier评分 ;

Ⅰ B5, UNC _____Brier评分的不确定性项,值介于 0~1 之间,值越小代表不确定性越大 ;

Ⅰ B5, REL _____Brier评分可靠性项,值介于 0~1 之间,值越小代表预测可靠性越高,完美预测时为 0 ;

Ⅰ B5, RE5 _____Brier评分解析度项,值介于 0~1 之间,值越大代表预测解析度越高。

G. 3 Brier技巧评分

某一分类事件多个预测实例的 Brier技巧评分按公怯(G. 6)H算 :

Ⅰ B ……………… (G. 6)

Ⅰ B55 _____某一分类事件多个预测实例的 Brier技巧评分 ;

Ⅰ B5 _____某一分类事件多个预测实例的 Brier评分 ;

Ⅰ B55, CLⅠ _____ 以气候概率作为预测的 Brier评分 ;

Ⅰ B5, RE5 _____Brier评分的解析度项 ;

Ⅰ B5, REL _____Brier评分的可靠性项 ;

Ⅰ B5, UNC _____Brier评分的不确定性项。

18

附录 H (规范性)

排序概率技巧评分计算

H. 1 累计实况概率

某一类事件某个预测实例的累H实况概率按公怯(H. 1)H算 :

Ok oi ……………… (H. 1)

Ok _____ 累H至第 k类事件的实况概率 ;

k _____ 累H的事件类别数 ;

i _____第 i类事件 ;

oi _____第 i类事件的实况概率,发生为 1,未发生为 0。

H. 2 累计预测概率

某一类事件某个预测实例的累H预测概率按公怯(H. 2)H算 :

Yk yi ……………… (H. 2)

Yk _____ 累H至第 k类事件的预测概率 ;

yi _____第 i类事件的预测概率。

H. 3 排序概率评分

多分类事件整体某个预测实例的排序概率评分按公怯(H. 3)H算 :

Ⅰ RP ……………… (H. 3)

Ⅰ RP5 _____ 多分类事件某个预测实例的整体排序概率评分 ;

K _____ 事件分类总数 ;

k _____ 累H的事件类别数 ;

Yk _____ 累H至第 k类事件的预测概率 ;

Ok _____ 累H至第 k类事件的实况概率。

H. 4 气候概率预测的排序概率评分

按气候概率预测多分类事件整体的排序概率评分按公怯(H. 4)H算 :

Ⅰ RP5C ……………… (H. 4)

Ⅰ RP5C _____气候概率预测的排序概率评分 ;

K _____ 事件分类总数 ;

19

Pk _____气候概率预测累H至第 k类事件的预测概率 ;

Ok _____ 累H至第 k类事件的实况概率。

H. 5 某个网格点 RPSS

某个网格点多分类事件整体多个预测实例的 RP55按公怯(H. 5)H算 :

Ⅰ RP55, i ……………… (H. 5)

Ⅰ RP55, i _____第 i个网格点多分类事件整体多个预测实例的 RP55 ;

〈ⅠRP5, i 〉_____第 i个网格点多个预测实例平均的排序概率评分 ;

〈ⅠRP5c, i 〉_____第 i个网格点多个预测实例平均气候概率预测的排序概率评分。

H. 6 某个区域 RPSS

某个区域平均多分类事件整体多个预测实例的 RP55按公怯(H. 6)H算 :

Ⅰ RP ……………… (H. 6)