1. 引言
近百年来全球温度经历了以气候增暖为主要特征的变化,气候变化成为国际上关注的重大问题。我国专家学者亦积极关注中国的气候变化问题以及相关研究。在建立中国区域气温序列的基础上,对中国的气温变化进行了较为系统的研究分析,得出中国区域气温变暖趋势显著且高于全球平均水平 [1] 。Cao et al. (2013) [2] 基于均一化的气温序列发现中国东部近百年气温增暖速率高达1.52℃/100a。但在百年尺度上,无论是全球变暖还是中国气候变暖都存在不确定性 [3] 。一方面由于台站空间分布不均匀,观测时次和时制不统一,观测仪器变更,站点迁移或城市化引起的观测场环境变化等造成气象资料存在较大的不确定性;另一方面,气候变暖的原因亦存在不确定性,主要是模拟气候变暖的模式存在不确定性、情景设置存在不确定性、气候序列归因分析也存在不确定性。气候时间序列往往是非平稳的并且存在多种时间尺度信息 [4] [5] [6] [7] 。在包含多种时间尺度信号和模态的时间序列中,仅仅关注其线性特征,显然其分析结果会存在不确定性。为了得到更加客观的结果和结论,需对气候序列的相关特征结构以及非线性特征进行分析。
在气候序列的相关中,既有短程相关亦有长程相关。大量研究表明,大气环流系统要素、气候序列以及极端气候序列存在显著的长程相关性 [8] [9] [10] [11] [12] 。本论文在探讨中国东部百年气温线性变化规律的基础上,研究了中国东部18个百年序列站的单站长程相关性,并且探索了其年代际尺度长程相关性规律,揭示了中国百年气温在不同时空尺度上的相关结构特征,为分析中国东部百年尺度气候特征提供新的方法和思路。
2. 资料和方法
本文使用国家气象信息中心(Chinese National Meteorological Information Center, NMIC)整理的“中国中东部百年均一化气温月值序列(V1.0)”中的中国东部16个站点1916.01~2015.12共100年的月平均气温数据。该数据集1951年以前的数据主要来源于国家气象信息中心在2002年建立的“中国长年代温度、降水数据集”,在2008年建立的“中国24城市长年代气温数据集”,2009年建立的“全国60个重点城市长时间序列气温和降水数据集”,以及中国科学院大气物理研究所在1991年建立的“中华人民共和国的两个长期仪器记录的气候数据库”的气温资料。1951年之后的资料来自于经过严格质量控制的中国地面气象台站观测资料。此外,中国香港与澳门的气温资料分别来自各自网站。该数据集收集整编不同来源的百年尺度器测气温资料,经过资料融合拼接,严格质量控制,综合利用标准序列法、偏最小二乘回归、多元线性回归及梯度距离权重反比法等插补方法进行缺测值插补,在建立中国近百年气温月值序列基础上,挑取1951年之前资料序列超过30年且质量较好的18个气象台站资料进行均一化检验与订正。本文选用的中国东部16个台站的经纬度和资料时长如表1 (18站点中贵阳和昆明资料时长相对较短,且处于西南地区);其站点空间分布如图1所示。
Table 1. The location and time of the century-old station in eastern China
表1. 中国东部百年台站的位置及时长
Figure 1. The distribution diagram of 16 site locations
图1. 16站点位置分布示意图
去趋势涨落分析方法(detrended fluctuation analysis,DFA)是研究非线性时间序列长程相关性的工具。它可以有效滤去各阶趋势成份,消除其中的伪相关现象,检测非平稳时间序列中的相关性,已广泛应用到自然科学甚至社会科学各个领域 [13] [14] [15] 。其计算过程如下:
第一步,计算时间序列(F(s), s)的累积离差,
(1)
其中
。
第二步,把Y(i)分成Ns个不重叠的等间隔s的区间v,其中
。
第三步,对于每个区间v,用最小二乘法拟合数据,得到局部趋势。滤去该趋势后的时间序列记为YS(i),表示原序列与拟合值之差,
(2)
其中
为第v区间的拟合多项式。
第四步,计算每个区间滤去趋势后的方差,
(3)
第五步,对所有等长度区间的方差求均值并开方,计算标准DFA波动函数,
(4)
若时间序列是长程相关,则DFA波动函数F(s)与滞后时间s成幂律关系
(5)
在双对数坐标(F(s), s)中,用最小二乘法拟合,其直线部分的斜率即为标度指数
。当
时,表示时间序列是非持久的,只有短期记忆性,当前事件不会对长期的未来事件产生影响;当
时,表示原序列是白噪声,例如布朗运动;当
时,表示序列具有长程相关特征,即当前发生的事件和未来事件之间存在长程相关性,时间序列具有长期记忆性,且标度指数
越大记忆性越好。
3. 中国东部百年温度变化规律分析
为便于分析和比较,本文选取了自1916年1月至2015年12月共100年的月平均温度序列进行研究。首先,对16站点的百年温度均值、年际方差及年均值的趋势进行了统计,如表2所示。
Table 2. Centennial mean, variance and trend of temperature in stations in eastern China
表2. 中国东部台站温度百年均值、方差及趋势
从表2中可以看出16站百年温度序列呈现以增温为主要特征的趋势,但是长沙和南京两站百年温度趋势为负值,主要是因为该两个站点的温度在1980s前主要是降温趋势,之后虽然增温趋势显著,但整体上趋势为负,其中长沙站温度资料演变如图2(c)所示。本文选择香港、沈阳、长沙、徐家汇四个台站为例,展示了站点100年的年平均温度随时间的变化情况,从图中可以看出,香港、沈阳两站在100年时间尺度内年平均温度总体呈现上升趋势;长沙站在1980s前温度呈现下降趋势,1980年后才明显上升(图2(c));徐家汇站在100年时间内发生了两次较为明显的转折,一次在1940s左右,另一次出现了1980s前后(图2(d))。从这四个站点的变化情况可知,不同区域的温度在100年时间尺度内变化规律存在显著差异,近百年的全球增暖需要结合区域特征来分析对气候系统的影响。
Figure 2. Annual average temperature of 100 years of temperature at different sites: (a) Hong Kong; (b) Shenyang; (c) Changsha; (d) Xujiahui
图2. 不同站点百年温度年均值状况:(a)香港;(b)沈阳;(c)长沙;(d)徐家汇
4. 中国东部温度长程相关性分析
本文首先计算了16个台站百年尺度的温度标度指数,用以分析中国东部各个台站温度百年尺度的长程相关性;鉴于台站温度变化存在转折的情况,且绝大多数仅有一次转折,因此计算了前后50年的标度指数以及滑动窗口为10年的各个时期50年的标度指数,用以分析各个时期台站温度长程相关性。
4.1. 中国东部百年和50年温度序列长程相关性
图3所示为香港、沈阳、长沙、徐家汇四个代表站点波动函数幂率分布,计算得到四个站点月平均温度序列标度指数分布为0.73、0.69、0.75、0.72,均表现出一定的长程相关特征。
Figure 3. Long-range correlation of monthly temperature-month series of different sites: (a) Hong Kong; (b) Shenyang; (c) Changsha; (d) Xujiahui
图3. 不同站点百年月温度月值序列长程相关性:(a)香港;(b)沈阳;(c)长沙;(d)徐家汇
4.2. 中国东部百年温度序列长程相关性
进一步计算出全部16站点百年温度月均值序列的标度指数,根据部分站点温度变化的转折,同时将百年尺度的时间分为前50年和后50进行了分别计算。结果如表3所示。
Table 3. 100-year temperature monthly mean sequence long-range correlation scale index of 16 sites
表3. 16站点百年温度月均值序列长程相关标度指数
从表中可以看出,各站标度指数均大于0.5,表明都具备长程相关特征,存在较为明显的可预测性。
4.3. 50年滑动月平均序列DFA
进一步,选取了50年时间窗口以10年为步长进行滑动平均,计算不同时段各站点标度指数分布,如图4所示。从图中可以发现,在100年时间尺度内,中国东部16站点的温度序列长程相关性存在转折,具体而言对于大部分站点,其长程相关性在前70年(1916~1986年)相对较为稳定,但是进入1980s,大部分站点长程相关性明显减弱,直到最近10年标度指数又出现整体增长趋势,这可能与我国1980s初改革开放,大规模工业化及城市建设有关,很可能在一定程度上反映了人为因素对区域气候系统状态的影响。此种转折关系在许多站点的温度序列中往往没有反映(如图5(a)、图5(b))。
Figure 4. 100-year temperature monthly average sequence 50-year sliding scale index of 16 sites (sliding interval is 10 years)
图4. 16站点百年温度月均值序列50年滑动标度指数(滑动间隔为10年)
Figure 5. 50-year sliding scale index (sliding interval is 10 years) for a hundred-year temperature-average sequence of different regional sites: (a) 3 sites in northeast China; (b) 3 sites in the pearl river delta region
图5. 不同区域站点百年温度月均值序列50年滑动标度指数(滑动间隔为10年):(a)东北地区3站点;(b)珠江三角洲地区3站点
5. 结论和讨论
本文对中国16个具有百年以上温度观测序列的站点的长程相关性进行了对比研究。结果发现,从区域分布来看,北方地区站点的长程相关性要弱于南方地区,从时间分段来看,前50年的长程相关性总体要比后50年弱。进一步计算了50年窗口10年步长的标度指数演变情况,发现绝大部分站点标度指数呈现先减少后增加的趋势,这反映了中国东部地区百年温度变化具有年代际特征,并且在增暖背景下存在一定的区域变化差异。
基金项目
国家自然科学基金(批准号:41675107、41775081和41530531)资助。