不同分辨率再分析资料对浙江省气温刻画能力的对比评估

第36卷第3期2020年6月
热带气象学报
JOURNAL OF TROPICAL METEOROLOGY
Vol.36,No.3Jun.,2020
殷悦，马浩，葛敬文，等.不同分辨率再分析资料对浙江省气温刻画能力的对比评估[J].热带气象学报，2020，36(3)：377-388.
文章编号：1004-4965(2020)03-0377-12
不同分辨率再分析资料对浙江省气温刻画能力的
对比评估
殷悦1，马浩1，葛敬文1，严睿恺2，高大伟1，孟仲1，金希3
(1.浙江省气候中心，浙江杭州310017；2.浙江省气象台，浙江杭州310017；3.杭州市气象局，浙江杭州310051)
摘要：利用浙江省66个基本气象站1979—2010年的日平均气温数据，系统评估了三套再分析资料R1、
R2和CFSR 对浙江省气温的刻画能力。

结果表明：三套再分析资料的气候平均态与观测均存在一定差异，其中R1、R2的空间分布型与观测较为接近，CFSR 与观测差异较大；三套再分析资料均存在系统性冷偏差且这一偏差在32年中稳定存在，其中CFSR 的冷偏差更显著，浙南地区是其冷偏差的重要来源。

三套资料的均方根误差均存在季节变化：冬季（特别是1月）误差较小而夏季（特别是7-8月）误差较大，R1和R2的季节差异强于CFSR 。

CFSR 对浙江省气温变率的把握能力优于R1和R2，其距平场EOF 分解前三模态的空间型态和时间系数与观测更为接近。

系统误差订正后，三套再分析资料的可信度得到显著改善，CFSR 的改善效果最明显，说明系统性误差是三套再分析资料偏差的重要来源。

改善后三套再分析资料的均方根误差和空间相关系数大体相当。

CFSR 网格点气温插值到观测站点时因海拔差异导致的误差以及CFSR 在浙江省的模式地形偏高可能是其有较大冷偏差的重要原因。

关
键
词：资料偏差；平均气温；CFSR ；R1；R2
中图分类号：P423；
文献标志码：A
Doi ：10.16032/j.issn.1004-4965.2020.036
收稿日期：2019-07-24；修订日期：2020-02-18
基金项目：浙江省基础公益研究计划项目（LGF19D050001）；中国气象局预报员专项项目（CMAYBY2019-048）；国家重点研发计划项目
（2018YFC1505600）；浙江省气象科技计划项目（2018QN05）共同资助
通信作者：马浩，男，安徽省人，博士，主要从事气候动力学、气候预测和极端天气气候事件研究。

E-mail ：************************
1引言
大气再分析资料是融合多源观测数据和数值模拟结果、基于资料同化技术建立的一套网格化长序列气象数据集，它有效弥补了站点观测时空分布不均的缺陷，是气候诊断分析的重要基石[1]。

目前，全球大气再分析资料自1990年代中期以来已历经三代的发展[2]。

第一代产品主要有美国国家环境预测中心（National Center of Environment Prediction ，NCEP ）和美国国家大气科学研究中心（National
Center
of
Atmospheric
Research ，
NCAR ）研发的NCEP-NCAR Reanalysis Ⅰ（以下简
称R1）[3-4]
和NCEP-DOE Reanalysis Ⅱ（以下简称R2）[5]；第二代产品主要有欧洲中期天气预报中心
（European Center for Medium -Range Weather Forecasts ，ECMWF ）研发的ERA-40[6]和日本气象厅（Japan Meteorological Agency ，JMA ）研发的JRA -25[7]；第三代产品主要有NCEP 研发的CFSR [8]、ECMWF 研发的ERA-Interim [9]、JMA 研发的JRA -55[10]以及美国航空航天局（National Aeronautics and Space Administration,NASA ）研发的MERRA [11]。

再分析资料的质量受到观测数据来源、数值预报模式和同化方法的影响，因此对
热带气象学报第36卷
再分析资料进行可信度评估尤为重要[12]。

R1自1990年诞生以来，历经近30年的发展已成为目前国际上使用最广泛的再分析资料。

它的优势是时间跨度最长（1948年1月至今）、使用群体庞大从而成为许多研究的参照基准。

然而R1存在不可忽视的人为误差，且在南半球的可信度较差，因此对R1进行优化改进十分必要。

R2则修正了R1中的人为误差，进一步改进了模式并同化了更多的观测资料，使得R2在土壤湿度、冬季降水、地面温度和高纬度地面通量等方面较R1得到较大改善[5]。

CFSR作为第三代全球再分析资料，其优势包括全面改进了模式框架和同化技术，对大气、陆地、海洋、海冰采用完全耦合方案，直接同化了卫星观测的原始数据，分辨率在R1、R2基础上显著提升（由水平方向2.5°×2.5°、竖直方向28层提升为水平方向0.5°×0.5°、竖直方向64层）等[8]。

然而，目前对CFSR进行可信度评估的研究相对较少，特别是目前对其气象要素的刻画能力所知有限。

目前国际上广泛开展的对再分析资料的评估表明，第三代再分析资料的性能较前两代有明显改进，特别是在较大的空间尺度上（如全球和中国尺度），第三代再分析资料能够给出更多的细节特征，同化结果与观测更加接近[13-21]。

然而，另一些评估结果却指出，第三代再分析资料在某些方面弱于第一、二代产品。

Wang等[13]总结了CFSR存在的问题，包括高估了西太平洋暖池区域的短波辐射通量、全球大部分地区的潜热通量偏大、1998—2001年期间许多变量发生了跳跃性突变等。

He等[22]利用中国2425个常规气象站数据对R2和CFSR给出的日平均气温进行了评估，结果表明R2的日平均气温在中国中部和东部的大部分地区以及西北地区中东部优于CFSR，CFSR的气温日较差在中国大部分地区可信度较低。

Zhao等[23]基于中国2419个气象站气温资料对比评估了季节尺度上R2和CFSR的表现，发现虽然CFSR在所有季节的日平均气温比R2更加接近观测，但R2在所有季节的日最高/最低气温比CFSR可信度更高。

综上所述，尽管第三代再分析资料在观测数据来源、数值模式分辨率、参数化方案以及同化技术等方面做出了重要改进，其表现是否一定优于
较早产品可能还与研究的区域和对象密切相关。

需要指出的是，相较于全球尺度和中国尺度，在区域尺度上开展再分析资料性能评估的研究较少，而这一工作对于在区域尺度上应用再分析资料开展诊断分析有重要意义。

有研究针对辽宁省海岸带[24]、天山西部[25]、环渤海区域[26]、黄渤海区域[27]再分析资料的适用性进行了评估，但大多基于CFSR、ERA-interim、JRA-55等第三代产品，而非不同年代际的再分析资料进行对比。

较新的高分辨率资料在区域尺度上是否一定优于较早的低分辨率资料并无定论，因此是一个值得探索的科学问题。

浙江省地处我国江南地区东部，下垫面复杂多变，次网格过程参数化难度较大，是检验再分析资料对区域气候刻画能力的试金石。

日平均气温（以下简称气温）作为描述基本气候特征的核心要素之一，其可信度是评估再分析资料性能的重要参考依据。

本研究基于浙江省66个国家基本气象站和R1、R2、CFSR三套同源于NCEP的再分析资料的气温数据，利用多种统计方法对近地面2m 气温场进行可信度评估，以了解高、低分辨率资料对浙江省气温刻画能力的差别，并进一步探讨如何基于站点观测提升再分析资料对区域气温的再现能力。

2资料和方法
2.1观测资料
本文采用的观测资料为浙江省66个国家基本气象站1979—2010年的逐日气温数据。

该66站为浙江省“一类观测站”，观测质量和数据可信度较高。

站点分布和相关信息见图1和表1。

经检查，该66站在研究时段内的气温共有16个缺测值。

考虑到气温在水平方向上大致呈线性变化，利用缺测当天非缺测站点的气温、采用基于径向基函数的Multiquadric插值方法[28]对缺测站点进行水平插补。

2.2再分析资料
本文选取R1、R2、CFSR三套再分析资料1979年1月—2010年12月2m高度月平均气温数据开展分析。

R1、R2的水平分辨率为2.5°×2.5°，CFSR的水平分辨率为0.5°×0.5°。

三套再分析资
378
第3期殷悦等：不同分辨率再分析资料对浙江省气温刻画能力的对比评估
料均来自NOAA网站（）。

为了对比再分析资料和站点观测数据的差异，首先将格点尺度的再分析气温资料插值到66个常规观测站上，所用插值方法为三次样条插值[29]。

图1浙江省66个国家基本站分布表1浙江省66个国家基本站信息
站号58443 58446 58448 58449 58450 58451 58452 58453 58454 58455 58456 58457 58458 58459 58464 58467 58468 58472 58477 58537 58542 58543站名
长兴
安吉
临安
富阳
湖州
嘉善
嘉兴
绍兴
德清
海宁
桐乡
杭州
海盐
萧山
平湖
慈溪
余姚
嵊泗
定海
开化
桐庐
淳安
海拔高度/m
25.00
72.20
117.60
47.90
4.10
2.60
4.80
7.90
102.00
5.30
6.00
43.20
4.80
96.50
4.00
5.70
38.00
79.60
35.70
153.60
44.80
172.20
站号
58544
58546
58547
58548
58549
58550
58553
58555
58556
58557
58558
58559
58562
58563
58565
58567
58569
58570
58631
58632
58633
58642
站名
建德
浦江
龙游
兰溪
金华
诸暨
上虞
新昌
嵊州
义乌
东阳
天台
鄞州
北仑
奉化
宁海
石浦
普陀
常山
江山
衢州
武义
海拔高度/m
87.20
85.20
66.20
48.30
64.70
39.10
6.40
115.10
104.30
90.00
91.90
107.60
6.20
5.00
41.50
25.00
129.20
85.20
137.00
126.30
82.40
103.80
站号
58643
58644
58646
58647
58652
58654
58656
58657
58658
58659
58660
58664
58665
58666
58667
58742
58745
58746
58750
58751
58752
58760
站名
永康
遂昌
丽水
龙泉
仙居
缙云
乐清
青田
永嘉
温州
临海
温岭
洪家
大陈岛
玉环
云和
庆元
泰顺
文成
平阳
瑞安
洞头
海拔高度/m
102.90
238.60
59.70
222.10
83.00
179.10
60.80
57.80
34.30
28.30
6.60
35.30
4.60
86.20
95.90
159.60
400.80
538.90
105.40
254.00
39.70
68.60
31°N
30°N
29°N
28°N
27°N
118°E119°E120°E121°E122°E123°E
km
02550100
379
热带气象学报第36卷
2.3主要分析方法2.
3.1时间相关系数
时间相关系数r 能够反映观测资料和再分析资料在时间上的相关性，对于两列长度相同的观测和再分析数据而言，其时间相关系数r 表征如下：
r =
∑i =1
N (X
i
-X
ˉ)(Y i -Y ˉ)∑i =1
N ()X
i
-X
ˉ2
∑i =1
N ()Y -Y
ˉ2
(1)
式中，N 为数据样本个数，
X i 、Y i 分别为气温的观测数据和再分析数据样本，X
ˉ、Y ˉ分别为气温观测数据和再分析数据序列的时间平均值。

2.3.2空间相关系数
将时间相关系数r 应用到空间上，则可以得到空间相关系数（Spatial Correlation Coefficient,SCC ）。

空间相关系数能够反映观测资料和再分析资料空间场之间的相关性。

对于两个站点数目相同的观测资料场和再分析资料场而言，其空间相关系数SCC 表征如下：SCC =
∑i =1
N (X i --X s )(Y i --Y s )
∑i =1N ()
X i --X s 2
∑i =1N
()
Y --Y s 2
(2)
式中，
N 为空间站点个数，X i 、Y i 分别为气温的观测数据和再分析数据样本，
-X s 、-Y s 分别为气温观测数据场和再分析数据场的空间平均值。

2.3.3均方根误差
均方根误差（Root-Mean-Square Error ，RMSE ）能够反映两列数据偏离的程度，本文中用于定量计算再分析资料相对于观测资料的偏差，具体公式表征如下：
RMSE =
(3)
2.3.4离差平方和
离差平方和（Sum of Squared Deviations,SSD ）能够反映两列数据之间的累积误差，用其刻画再分析资料相对于观测资料的累积偏差，具体公式表征如下：
SSD =∑i =1
N (Y i -X i )2
(4)
2.3.5经验正交函数分解
经验正交函数（Empirical Orthogonal
Function ，EOF ）分解是一种在气候诊断分析中广泛应用的方法，其主要思想是将复杂的气象变量场分解为一系列相互正交的空间分布型和时间序列，通过计算累积方差贡献提取前几个空间模态和时间系数用以反映变量场主要的时空变化特征，从而达到降维的目的[30]。

EOF 分解的计算方
法详见文献[30]。

3结果与分析
3.1气候平均态
从1979—2010年32年的观测及三套再分析资料中气温的时空分布特征可以看出（图2），观测气温的空间分布大致遵循纬度地带性规律变化[31]，即由北向南气温逐渐升高（图2a ）。

高值中心主要位于浙中、浙南地区，最高值18.6℃出现在丽水青田站；低值中心主要位于浙北北部及衢州地区，最低值16.0℃出现在湖州安吉站。

R1、R2的气温空间分布也呈现出阶梯变化的特点，但并未刻画出纬向分布的特征，而是表现为自西北向
东南气温逐渐升高的分布形态，与观测存在一定差异；此外，R1、R2的气温系统性偏低（图2b 、2c ）。

与西北-东南型分布特征相一致，R1、R2的气温高值区位于东南沿海地区，低值区主要分布在浙西北地区。

与站点观测和R1、R2均不同，CFSR 刻画的气温不仅明显偏低，而且除沿海地区外并不具备阶梯分布特征，在浙西北地区、浙西南地区和浙东南地区出现了三个显著的低值中心（图2d ）。

就气候态而言，尽管三套再分析资料与站点观测的空间分布形态均存在一定差异，但R1、R2的气温分布特征与观测更为接近；虽然CFSR 由于
具有较高分辨率而刻画了更多的细节，却与观测资料反映的南暖北冷的特征差异很大。

考虑到浙江地形以小型山地和丘陵为主，观测站大部分建在平原或低海拔处，因此CFSR 拥有更显著冷偏差的原因可能包括其部分网格点海拔较高，而将网格点气温插值到观测站经纬度时，仅考虑了水平方向插值，未考虑竖直方向上因海拔差异造成的
气温差异，这可能导致插值到站点上的气温偏低（当观测站点的海拔低于周围网格点时）。

相较而言，R1、R2反而由于网格点间距更大，不能捕捉到浙江的地形起伏，从而规避了插值时与观测站点
380
第3期殷悦等：不同分辨率再分析资料对浙江省气温刻画能力的对比评估海拔差异问题。

从全省平均气温来看（图2e~2g ），三套再分析资料的逐年变化趋势与观测非常一致，R1、R2、CFSR 与站点观测时间序列的相关系数分别达到0.97、0.95、0.97。

为排除气温线性趋势对相关系数的影响，去除了四套资料中的线性趋势。

R1、R2、CFSR 与观测时间序列的相关系数分别下降至
0.92、0.92、0.93，但仍呈高相关关系（通过了0.01显著性水平检验）。

然而，如前所述，三套再分析资料均表现出气温系统性偏低的特征。

R1、R2、CFSR 较观测资料平均偏低1.2℃、1.1℃、1.7℃。

R1、R2的偏低幅度相对较小，CFSR 偏低幅度较大，与三套资料的空间分布特征相符。

图21979—2010年32年平均气温的空间分布和年平均气温时间序列与站点观测对比
a.站点观测；b 、e.R1；c 、f.R2；d 、g.CFSR 。

3.2资料偏差时空分布特征3.2.1基本特征
RMSE 是反映资料偏差的良好指标。

由图3可见，R1、R2的RMSE 空间分布型态基本一致，均表现出西部大于东部的特征。

浙北西南部、浙西南大部分地区是RMSE 的相对大值区，最大值出现在衢州常山（2.7℃，R1）和金华兰溪（2.5℃，R2）；东部沿海地区是RMSE 的相对小值区；但总体而言全省范围内RMSE 的差别并不大，各地RMSE 多小于2℃。

与R1、R2相比，CFSR 的RMSE 显著增长，浙南地区RMSE 普遍在3℃以上，最大达4.2℃，远大于R1、R2的最大误差；次大值中心位于浙西北地区，与R1、R2的分布型有相似之处。

整体而言CFSR 浙北地区的RMSE 较
小，多分布在2℃以下，但仍明显大于R1、R2。

比较三套再分析资料，可知CFSR 的RMSE 远大于R1和R2，特别是在浙南地区。

RMSE 的年际变化具有相对稳定性。

从全省RMSE 的时间序列来看（图3d ），1979—2010年R1、R2的RMSE 年际变化相近，主要分布在1.2~1.8℃之间；CFSR 的RMSE 明显大于R1和R2，主要分布在1.9~2.2℃之间，平均值达到2.1℃，而R1和R2的平均值仅为1.5℃和1.4℃。

由此可知，与R1和R2相比，CFSR 的RMSE 偏大并非少数年份出现的个例现象，而是在32年中稳定存在，这一系统性偏差的存在极大地影响了CFSR 气温资料在浙江地区的适用性。

31°N 30°N 29°N 28°N 27°N
118.5°E 120.5°E 122.5°E 118.5°E
120.5°E 122.5°E 118.5°E 120.5°E 122.5°E118.5°E
120.5°E 122.5°E
13.013.514.014.515.015.516.016.517.017.518.018.5
气温/℃
OBS
R1
R2
CFSR
1918171615
14
气温/℃198019851990199520002005201019801985199019952000200520101980198519901995200020052010
时间/年
R1VS OBS
R2VS OBS
CFSR VS OBS
(a)(b)
(c)
(d)
(e)r =0.97
(f)r =0.95
(g)r =0.97
381
热带气象学报
第36卷
为了进一步探究CFSR 误差的空间分布特征，选取1979—2010年RMSE 大于2.8℃的9个站点（图4a 中黄点，以下简称浙南9站），考察浙南地区误差对浙江全省误差的贡献。

经计算可知，浙南9站的逐年RMSE 分布在3.3~3.9℃之间，远大于其它地区的RMSE （图4a ）；进一步计算浙南9站SSD 与全省SSD 之比，发现这一地区的站点数目虽然仅占全省的14%，却贡献了全省SSD 的39%~45%（图4b ），再次说明浙南误差是全省误差的重要来源，南
部RMSE 明显偏大是全省RMSE 整体较高的重要原因。

进一步考察浙江地形分布发现，偏差较大的区域恰好位于海拔较高地带。

这在一定程度上证实了前述猜想，即CFSR 冷偏差较大可能由于在将网格点气温值插值到站点时，未对因两地海拔差异造成的气温差异进行修正。

此外，前人研究指出，模式地形偏差可能是造成再分析资料误差的重要原因[33-34]，因此可能CFSR 在这些区域本身的模式地形偏高，从而产生了较大的系统性冷偏差。

图3R1（a ）、R2（b ）、CFSR （c ）1979—2010年32年平均气温RMSE 的空间分布及逐年气温RMSE 时间序列（d ）
图4浙江南部气温RMSE >2.8℃的站点位置（a ，黄点）；1979—2010年浙南9站气温RMSE 及SSD 全省占比时间序列（b ）
3.2.2月季变化
为进一步了解三套再分析资料在月季尺度上的表现，分别计算了三者在春季（3—5月）、夏季（6—8月）、秋季（9—11月）、冬季（12月—次年1月）
以及各月的RMSE 。

可以看出，R1、R2的四季RMSE 演变特征相似：秋冬季节RMSE 较小、春夏季节RMSE 较大，RMSE 在夏季和冬季分别达到最大值和最小值。

R1秋冬季节的RMSE 在1℃左
31°N 30°N 29°N 28°N 27°N
118.5°E 120.5°E
122.5°E 118.5°E 120.5°E 122.5°E 118.5°E 120.5°E 122.5°E
温度/℃
0.0
0.5
1.0
1.5
2.0
2.5
3.0
3.5
4.0
4.52.62.21.81.41.0
R M S E /℃
19801983198619891992199519982001200420072010
时间/年
(a)(b)
(c)
(d)
R1
R2
CFSR
R1
R2
CFSR
31°N 30°N
29°N 28°N 27°N
(a)(b)
19801985
1990
1995
2000
2005
2010
时间/年
118.5°E
119.5°E 120.5°E 121.5°E
122.5°E
0.00.40.81.21.62.02.42.83.23.64.04.4
温度/℃
R M S E /℃
4.03.93.83.73.63.53.43.3
0.450.440.430.420.410.40
0.39RMSE
SSD 占比
382
第3期殷悦等：不同分辨率再分析资料对浙江省气温刻画能力的对比评估右，R2则介于1~1.5℃之间，R1较R2偏差更小；春季R1的RMSE 介于1.5~2℃之间，R2则在1.5℃附近，R2的偏差更小；夏季二者的差异更加明显，R1的RMSE 分布在2~3℃之间，R2则分布在1.5~2.5℃之间。

总体而言虽然R1和R2的RMSE 均表现出一定的季节性，但R1的季节差异更为明显。

与R1和R2相比，CFSR 的RMSE 季节变化相对较弱，四季RMSE 均基本在1.5~2.5℃之间；冬季最小，春季次之，夏季与秋季相当且最大。

需要指出的是，CFSR 在冬季、春季、秋季的偏差均大于R1和R2，但其夏季的偏差与R1相当、略高于R2，说明与其它两套再分析资料相较，CFSR 的劣势主要在于冬、春、秋三季。

尽管三套再分析资料误差的季节变化存在一定差异，其共同特点是冬季误差较小、夏季误差较大。

施晓晖等[32]在研究R1在中国区域的可信度时指出其在冬季的可信度最好，而在夏季可信度较差，这与本文研究结果相一致。

三套再分析资料RMSE 的分月变化表现出和季节变化相似的特征。

在所有月份中，误差较小的月份为12月、1月和2月，其中1月误差最小。

CFSR 的月际起伏明显小于R1和R2。

R1、R2的高值月份出现在7月和8月，CFSR 则出现在9月，这也是R1和R2夏季误差最大而CFSR 秋季误差最大的原因。

从逐月误差的上下限（即所有站点当月RMSE 的最大值和最小值）来看，三套资料各月误差的下限相近，但CFSR 的上限明显高于R1和R2，且各月变化较小，R1上限的起伏较R2更明显，二者均表现出7—8月上限较高、其它月份上限相对较低的特点，说明夏季7—8月R1和R2中各站RMSE 的差异较大。

综上所述，三套资料均表现出冬季偏差最小而夏季偏差较大（特别是7—8月）的特征；CFSR 在冬、春、秋三个季节中RMSE 均大于R1和R2，而在夏季与R1和R2相当，说明CFSR 的偏差主要来源于于冬、春、秋三季。

图51979—2010年R1（a 、d ）、R2（b 、e ）、CFSR （c 、f ）春季、夏季、秋季、冬季气温RMSE 时间序列（a 、b 、c ）
（d 、e 、f ）和32年平均1—12月气温RMSE （黑色实线）及其上下限分布（灰色误差线）
3.2.3空间相关性
再分析资料与站点观测资料在空间分布上的相似程度也是资料偏差的重要方面。

由三套资料与站点观测的逐年SCC 可知，R1、R2与站点观测之间存在较好的相关，SCC 始终保持在0.25以上
（通过了0.05显著性水平检验），且多数年份达到0.3以上（通过了0.02的显著性水平检验）。

二者的SCC 时间序列较为接近，平均值分别为0.42和0.44；1979年、1995年、2005年为三个低值年份，1991年、1999年、2001年为三个高值年份。

SCC
3.02.52.01.51.00.5
R M S E /℃
3.02.52.01.51.00.5R M S E /℃
19801985
1990
19952000
2005
2010
时间/年R1
R1
时间/月
1
2
3
4
5678
9
101112
R2
R2
3.02.52.01.51.00.5R M S E /℃
3.02.52.01.51.00.5
R M S E /℃
3.02.52.01.51.00.5R M S E /℃
3.02.52.01.51.00.5
R M S E /℃
19801985
1990
19952000
2005
2010
时间/年时间/月1
2
3
4
5
678
9
101112
CFSR
CFSR
1980
1985
199019952000
20052010
时间/年
时间/月
12345
6789101112
(a)
(d)
(b)
(e)
(c)
(f)
春季夏季秋季冬季
383
热带气象学报第36卷
表现出较明显的阶段性变化：1980年代—1990年代前期SCC 相对较高，1990年代中后期相对较低，1990年代末—20世纪初相对较高，之后转为相对较低。

与R1和R2相比，CFSR 和站点观测在空间上几乎不相关，SCC 在0值上下波动，最大值仅有0.16，2003年之后甚至SCC 进一步下降，说明CFSR 不仅RMSE 远高于R1和R2，在SCC 上也不及R1和R2。

2003年之后，三套再分析资料的SCC 无一例外出现下降，其中的原因值得进一步深入研究。

3.3EOF 分析
为探究三套再分析资料对浙江主要气候变率模态和其随时间演变的刻画能力，对观测资料和三套再分析资料的年平均气温距平进行EOF 分解（图7）。

为避免全球变暖信号的干扰，在EOF 分解之前去除了气温序列的线性趋势。

观测中的第一模态为“全省一致型”，大值区位于浙北东南部，方差贡献高达86.7%（图7a ），说明全省各站同位相变化是浙江气温变化的主要方式。

对应的时间系数表现出显著的年代际变化特征：1980年代为负位相，1990年代—2000年代中后期为正位相，之后重新转为负位相（图7d ）。

R1、R2和CFSR 的第一模态均表现为与观测一致的全省同位相变化（图7g 、7m 、7s ），且三者均较好地刻画了大值区的位置，其中CFSR 刻画的浙西南地区的小值区也与观测一致。

R1、R2、CFSR 第一模态的方差贡献分别达到96.9%、95.9%和90.0%；与观测相比，R1、R2的方差贡献过大，CFSR 的方差贡献与观测相近。

就时间系数而言，三套再分析资料均较好地再现了观测中的年代际变化特征，三者的PC1与观测PC1的相关系数分别达到0.92、0.92、0.93（表2）。

观测中的第二模态为“南北反相型”，约以杭州南部-绍兴南部为界呈现南北反位相变化的特征，方差贡献为4.3%（图7b ），其时间系数表现出显著的年代际转折特征：1980年代前期—1990年代中后期主要为负位相，之后转为正位相（图7e ）。

三套再分析资料同样抓住了南北反位相的空间分布特征，但总体而言，CFSR 与观测更为接近且方差贡献相当（5.4%）（图7h 、7n 、7t ）。

从时间系数来看，R1、R2均未能再现观测中的年代际转折特征、更多地表现为年际振荡（图7k 、7q ），CFSR 给出的时间系数与观测相近（图7w ）。

从时间系数的相关来看，R1、R2、CFSR 的PC2与观测PC2之间的相关系数分别为0.57、0.78、0.86，比较而言CFSR 优于R1和R2（表2）。

观测中的第三模态主要表现为“沿海-内陆反相型”，即沿海地区和内陆地区呈现反位相变化的特征，方差贡献为2.7%（图7c ）。

R1、R2和CFSR 均刻画出了沿海-内陆反相变化的特征，但CFSR 与观测相比零线过于靠近海岸线，R1和R2的空间分布型态与观测更为接近（图7i 、7o 、7u ）。

从时间系数来看，观测中1999年之前下降、1999年之后上升的年代际转折特征在三套再分析资料中均能得到体现（图7l 、7r 、7x ），R1、R2、CFSR 的PC3与观测PC3之间的相关系数分别为0.43、0.49、0.42（表2）；然而叠加在年代际背景之上较强的年际振荡仅在CFSR 中得到体现（图7f 、7l 、7r 、7x ）。

综合EOF 分析结果可知，CFSR 对浙江省气温的主要空间模态和时间系数的刻画优于R1和R2且与观测接近；CFSR 能够展现更多的空间分布细节特征，且其前三模态的时间系数均与观测吻合较好。

总体来看，尽管CFSR 存在显著的系统性冷偏差，但其刻画的气温变率比R1、R2更接近观测。

这提示我们，不同分辨率的再分析资料可能在不同方面各具优势。

在较小的区域尺度上使用再分析资料时，应当针对不同的问题选取不同的资料，
图61979—2010年逐年R1、R2、CFSR 的年平均气温与站
点观测值之间的SCC 时间序列
表2
R1、R2、CFSR 与站点观测EOF 分解前三模态时间系
数之间的相关系数
再分析资料
R1R2CFSR
PC10.920.920.93
PC20.570.780.86
PC30.430.490.42
19801983198619891992199519982001200420072010
时间/年
0.60.40.20.0-0.2
S C C R1VS OBS
R2VS OBS
CFSR VS OBS
384
第3期
殷悦等：不同分辨率再分析资料对浙江省气温刻画能力的对比评估
31°N 30°N 29°N 28°N
27°N
118.5°E 120.5°E 122.5°E (a)119.5°E 121.5°E 118.5°E 120.5°E 122.5°E 119.5°E 121.5°E 118.5°E 120.5°E 122.5°E 119.5°E 121.5°E (b)
(c)
(d)(e)(f)(j)(k)(l)(g)(h)(i)
(m)(n)
(o)
(p)(q)(r)(s)(t)
(u)
(v)(w)
(x)3
1-1
-3
31°N 30°N
29°N 28°N
27°N
198019841988199219962000200420081980198419881992199620002004200819801984198819921996200020042008118.5°E 120.5°E 122.5°E 119.5°E 121.5°E 118.5°E 120.5°E 122.5°E 119.5°E 121.5°E 118.5°E 120.5°E 122.5°E 119.5°E 121.5°E 118.5°E 120.5°E 122.5°E 119.5°E 121.5°E 118.5°E 120.5°E 122.5°E 119.5°E 121.5°E 118.5°E 120.5°E 122.5°E 119.5°E 121.5°E 118.5°E 120.5°E 122.5°E 119.5°E 121.5°E 118.5°E 120.5°E 122.5°E 119.5°E 121.5°E 118.5°E 120.5°E 122.5°E 119.5°E 121.5°E 198019841988199219962000200420081980198419881992199620002004200819801984198819921996200020042008198019841988199219962000200420081980198419881992199620002004200819801984198819921996200020042008198019841988199219962000200420081980198419881992199620002004200819801984198819921996200020042008
3
1-1
-3
31°N 30°N 29°N 28°N
27°N
31-1
-3
31°N 30°N 29°N 28°N
27°N
3
1-1
-3
温度/℃
0.45
0.400.350.300.250.200.150.100.050.00-0.05-0.10-0.15-0.20
图7去除线性趋势之后站点观测与三套再分析资料气温距平EOF 分解前三模态的空间型及时间系数
a~f.站点观测；g~l.R1；m~r.R2；s~x.CFSR 。

3.4系统误差订正效果评估
根据前述分析，三套再分析资料和站点观测之间存在系统性误差，故对系统误差进行订正并分析改善效果。

订正方法如下：计算每个站点上再分析资料和观测资料1979—2010年32年的气候态偏差，再将这一偏差叠加到再分析资料上。

订正系统误差后，三套再分析资料的多年平均RMSE 均得到明显改善（图8），全省大部分区域的RMSE 下降至0.1~0.3℃、仅R2中东部沿海地区存在0.5℃左右的偏差，其中CFSR 的改善较R1和R2更为明显，其RMSE 最大值、平均值分别由4.2℃、1.8℃减小至0.36℃、0.21℃，南部原先偏差较大的区域也得到显著改善。

更重要的是，订
正后CFSR 的RMSE 已与R1和R2相当甚至更小。

从时间演变来看，订正后CFSR 的年平均RMSE 在1996年之前略高于R1和R2，但1996年之后低于R1和R2，特别是2004年之后明显低于R1和R2。

此外，再分析资料与观测之间的相关性也得到了显著改善。

从SCC 逐年演变时间序列可以看出，订正之后R1、R2和CFSR 的逐年SCC 均达到0.9以上。

1980年代—1990年代中期三者SCC 相当，1994年CFSR 明显偏低，此后直至2004年三者再次相当，之后CFSR 高于R1和R2。

整体而言订正后CFSR 的误差水平已与R1和R2相当甚至更低，特别是2004年之后CFSR 已明显优于R1和R2。

385。