文章基于宿迁市4个气象站点1961~2021年逐日降水量、平均气温、最高和最低气温数据,运用线性趋势法、Mann-Kendall检验法、RClimDex模型、小波分析和R/S分析法,分析了18个极端气温和降水指数的变化特征,并进行了趋势预测。结果表明:1) 极端高温和极端低温均呈增加趋势,气温日较差呈减小趋势;极端高温(低温)事件的频率和持续性呈上升(下降)趋势;2) 极端降水事件的强度和频率呈微弱的上升趋势,持续性呈微弱的下降趋势;3) 极端气温冷指数的升温趋势明显高于暖指数,而夜指数的升温幅度也超过昼指数;4) 从周期性来看,极端气温指数存在2~6 a的周期变化;极端降水指数存在2 a左右的周期变化;5) R/S分析表明,未来宿迁地区气候变暖的状态仍继续保持;极端降水指数增减趋势不一致;极端气温指数的平均循环长度多大于15 a,极端降水指数的循环长度多小于10 a。 In this paper, based on daily precipitation, average temperature, and maximum and minimum temperature data of 4 meteorological stations in Suqian City from 1961 to 2021, linear trend method, Mann-Kendall test, RClimDex model, wavelet analysis, and R/S are used to analyze the variation characteristics and trend prediction of 18 extreme temperature and precipitation indexes. The results are as follows: 1) Both extreme high temperature and extreme low temperature showed an increasing trend, while the daily temperature range showed a decreasing trend; the frequency and sustainability of extreme high-temperature (low-temperature) events showed an increasing (decreasing) trend; 2) The intensity and frequency of extreme precipitation events are marginally increasing, while their sustainability has decreased slightly; 3) The increase in temperature for the extreme temperature cold indexes is higher than that of the warm indexes, while the increase in temperature for the night indexes is higher than that of the day indexes; 4) From a periodic perspective, the extreme temperature indexes exhibit a periodic variation of 2~6 a, while the extreme precipitation indexes display a periodic variation of about 2 a; 5) The R/S analysis shows that the state of climate warming in Suqian area will continue to persist, the trends of future increase or decrease of extreme precipitation indexes are inconsistent, and the average cycle length of the extreme temperature indexes is mostly greater than 15 a, while the cycle length of the extreme precipitation indexes are mostly less than 10 a.
文章基于宿迁市4个气象站点1961~2021年逐日降水量、平均气温、最高和最低气温数据,运用线性趋势法、Mann-Kendall检验法、RClimDex模型、小波分析和R/S分析法,分析了18个极端气温和降水指数的变化特征,并进行了趋势预测。结果表明:1) 极端高温和极端低温均呈增加趋势,气温日较差呈减小趋势;极端高温(低温)事件的频率和持续性呈上升(下降)趋势;2) 极端降水事件的强度和频率呈微弱的上升趋势,持续性呈微弱的下降趋势;3) 极端气温冷指数的升温趋势明显高于暖指数,而夜指数的升温幅度也超过昼指数;4) 从周期性来看,极端气温指数存在2~6 a的周期变化;极端降水指数存在2 a左右的周期变化;5) R/S分析表明,未来宿迁地区气候变暖的状态仍继续保持;极端降水指数增减趋势不一致;极端气温指数的平均循环长度多大于15 a,极端降水指数的循环长度多小于10 a。
极端气候,小波分析,变化特征,Hurst指数
En Xu, Hang Zhou, Xin Cheng, Liuxian Song, Xiaomin Ding
Suqian Meteorological Bureau, Suqian Jiangsu
Received: Feb. 20th, 2024; accepted: Mar. 21st, 2024; published: Mar. 28th, 2024
In this paper, based on daily precipitation, average temperature, and maximum and minimum temperature data of 4 meteorological stations in Suqian City from 1961 to 2021, linear trend method, Mann-Kendall test, RClimDex model, wavelet analysis, and R/S are used to analyze the variation characteristics and trend prediction of 18 extreme temperature and precipitation indexes. The results are as follows: 1) Both extreme high temperature and extreme low temperature showed an increasing trend, while the daily temperature range showed a decreasing trend; the frequency and sustainability of extreme high-temperature (low-temperature) events showed an increasing (decreasing) trend; 2) The intensity and frequency of extreme precipitation events are marginally increasing, while their sustainability has decreased slightly; 3) The increase in temperature for the extreme temperature cold indexes is higher than that of the warm indexes, while the increase in temperature for the night indexes is higher than that of the day indexes; 4) From a periodic perspective, the extreme temperature indexes exhibit a periodic variation of 2~6 a, while the extreme precipitation indexes display a periodic variation of about 2 a; 5) The R/S analysis shows that the state of climate warming in Suqian area will continue to persist, the trends of future increase or decrease of extreme precipitation indexes are inconsistent, and the average cycle length of the extreme temperature indexes is mostly greater than 15 a, while the cycle length of the extreme precipitation indexes are mostly less than 10 a.
Keywords:Extreme Climate, Wavelet Analysis, Variation Characteristics, Hurst Index
Copyright © 2024 by author(s) and beplay安卓登录
This work is licensed under the Creative Commons Attribution International License (CC BY 4.0).
http://creativecommons.org/licenses/by/4.0/
2023年3月,根据IPCC发布的第六次评估报告,2011~2020年全球地表平均气温较1850~1900年上升了1.1℃ [
已有研究表明,全球范围内,极端高温事件的频次和强度均呈现出明显的上升趋势,与此同时,极端低温事件的频次则呈现下降趋势,然而,其强度却在不断增强 [
在全球气候持续变暖的背景下,全球水循环受到深刻影响,进而导致全球范围内极端降水事件也不断增多 [
宿迁市地处江苏北部,沂沭泗流域和淮河流域的下游,属暖温带季风气候区,四季分明,季风盛行。宿迁市年均气温在14.7~15.3之间,无霜期长,年平均降水量为898.7 mm,年内降水主要集中于夏季,占全年的60%左右。本文根据宿迁市4个地面气象观测站1961~2021年逐日平均气温、降水量、最高和最低气温数据,选取ETCCDI推荐的18个极端气候指数,其中包括11个极端气温指数和7个极端降水指数,并用线性倾向估计法、M-K突变检验、滑动T检验、RClimDex模型、小波分析和R/S分析等对宿迁市极端气候(气温和降水)事件进行了全方位深入的研究,并对其未来趋势进行了预测,以便更有效地应对极端天气的发生,为宿迁市防灾减灾提供科学的依据。
采用宿迁市4个地面气象观测站1961~2021年逐日平均气温、降水量、最高和最低气温数据,计算由ETCCDI推荐的11个极端气温指数和7个极端降水指数,选取的18个极端气候指数及其定义见表1。
类型 | 名称 | 代码 | 定义 | 单位 | |
---|---|---|---|---|---|
极端气温指数 | 强度指数 | 极端最高气温 | TXx | 年内最高气温的最大值 | ℃ |
极端最低气温 | TNn | 年内最低气温的最小值 | ℃ | ||
气温日较差 | DTR | 日最高气温与日最低气温的差值 | ℃ | ||
频率指数 | 冷昼日数 | TX10p | 日最高气温 < 10%分位值的天数 | d | |
冷夜日数 | TN10p | 日最低气温 < 10%分位值的天数 | d | ||
暖昼日数 | TX90p | 日最高气温 > 90%分位值的天数 | d | ||
暖夜日数 | TN90p | 日最低气温 > 90%分位值的天数 | d | ||
高温日数 | SU35 | 每年日最高气温 > 35℃的日数 | d | ||
霜冻日数 | FD0 | 每年日最低气温 < 0℃的日数 | d | ||
持续时间指数 | 暖持续日数 | WSDI | 连续6日最高气温 > 90%分位值的日数 | d | |
冷持续日数 | CSDI | 连续6日最低气温 < 10%分位值的日数 | d | ||
极端降水指数 | 强度指数 | 单日最大降水量 | RX1day | 每年内最大日降水量 | mm |
5日最大降水量 | RX5day | 每年内连续5日最大降水量 | mm | ||
强降水量 | R95p | 降水量 > 95%分位值的总量 | mm | ||
降水强度 | SDII | 降水量 ≥ 1 mm的总量与天数的比值 | mm/d | ||
频率指数 | 暴雨日数 | R50 | 一年内日降水量 > 50 mm的日数 | d | |
持续时间指数 | 连续干旱日数 | CDD | 降水量 < 1 mm最长连续日数 | d | |
连续湿润日数 | CWD | 降水量 ≥ 1 mm最长连续日数 | d |
表1. 极端气候指数的定义
运用线性拟合对极端气候指数进行趋势分析;利用小波分析方法 [
等级 | Hurst指数值域 | 持续性强度 | 等级 | Hurst指数值域 | 持续性强度 |
---|---|---|---|---|---|
1 | 0.50 < H ≤ 0.55 | 很弱 | 1 | 0.45 ≤ H < 0.50 | 很弱 |
2 | 0.55 < H ≤ 0.65 | 较弱 | 2 | 0.35 ≤ H < 0.45 | 较弱 |
3 | 0.65 < H ≤ 0.75 | 较强 | 3 | 0.25 ≤ H < 0.35 | 较强 |
4 | 0.75 < H ≤ 0.80 | 强 | 4 | 0.20 ≤ H < 0.25 | 强 |
5 | 0.80 < H ≤ 1.00 | 很强 | 5 | 0.00 ≤ H < 0.20 | 很强 |
表2. Hurst指数分级
1961年以来,宿迁地区极端气温强度指数中,日最高气温极大值(TXx)和日最低气温极小值(TNn)均呈上升趋势(图1(a)~(b)),表明宿迁气温升高,呈增暖趋势,但是TNn的变幅明显高于TXx,呈显著不对称性,TXx近61 a来上升趋势不显著(P > 0.05),气候倾向率为0.08℃/10 a;TNn近61 a来呈显著增加趋势(P < 0.01),气候倾向率为0.6℃/10 a;气温日较差(DTR)具有显著下降趋势(P < 0.01),气候倾向率为−0.26℃/10 a (图1(c))。
极端气温的频率指数变化十分显著。其中,冷相关指数的冷昼(TX10p)、冷夜(TN10p)、霜冻日数(FD0)均呈显著减少趋势(P < 0.01) (图1(d)~(e),图1(i)),其气候倾向率分别为−1.76 d/10 a,−5.13 d/10 a,−6.89 d/10 a。暖相关指数中,暖夜日数(TN90p)以显著趋势上升(P < 0.01) (图1(g)),其增加速率为5.15 d∙10a−1。暖昼日数(TX90p)也呈现显著上升趋势(P < 0.05),其增加速率为1.44 d∙10a−1。超过35℃的高温日数(SU35)变化趋势不显著(图1(h)),近61 a来存在两个峰值区间,分别位于1961~1970年和2014~2015年,1970 s~2010 s之间整体呈小幅波动的态势。
图1. 近61年来(1961~2021)宿迁极端气温指数变化趋势
极端气温持续时间指数中,暖持续日数(WSDI)总体呈缓慢上升趋势(图1(j)),冷持续日数(CSDI)呈显著减少趋势(P < 0.01) (图1(k)),从趋势图上可以看出,1996年至今,除了2003年6天,冷持续日数一直为0 d。
综合以上分析,近61年来,宿迁地区极端高温和极端低温均呈增加趋势,极端低温升温幅度大于极端高温,同时气温日较差呈减小趋势。对于极端高温事件,无论是频率指标,还是持续性指标,均呈上升趋势,而极端低温事件的频率指标和持续性指标都呈下降趋势。在各极端气温指数中,暖夜日数上升速率最快,霜冻日数下降速率最快。
由小波功率谱可知(图2(a)),TXx存在4个显著聚集中心,其中心坐标分别为:(6, 1975)、(3, 1978)、(2, 1992)和(5, 2004)。同时,以每个坐标为中心的小波能量在时域尺度上的强集中影响范围分别是1965~1987年、1976~1981年、1983~1999年和1998~2021年。从全局小波谱中可以看出,TXx存在2 a、6 a、14 a和32 a左右的周期,其中2 a和6 a周期通过了95%显著性水平检验,而6 a震荡周期表现得尤为显著,因此TXx的主周期为6 a左右。
TNn也存在4个主要的聚集中心(图2(b)),其中心坐标分别为:(4, 1972)、(2, 1991)、(7, 1991)和(2, 2012)。以(4, 1972)为中心的气温指数小波能量呈现出强烈的集中态势,其影响范围主要集中在1961~1976年,以(2, 1991)和(7, 1991)为中心的强集中影响范围均是1986~2000年,以(2, 2012)为中心的强集中影响范围均是2006~2020年。从小波谱中看出,TNn具有2 a和32 a的振荡周期,其中2 a左右的振荡周期在95%的显著性水平上通过了检验,因此可以确定TNn的主周期约为2 a。
从图2(c)可以看出,DTR存在4个主要的聚集中心,其中心坐标分别为(4, 1966)、(2, 1982)、(5, 1988)和(2, 2000),这些中心的强集中影响范围为1961~1970、1978~1986、1984~1990和1992~2005年。但以(4, 1966)为中心的边界效应影响时段为1961~1968年。从全局小波谱中看出,DTR存在2 a、16 a和32 a的振荡周期,其中2 a周期通过95%显著水平检验,因此DTR的主周期约为2 a。其他极端气温指数的小波能量谱特征详见表3。
图2. 近61年来(1961~2021)宿迁极端气温强度指数小波能量谱分布
指数 | 集中影响年份 | 尺度中心/a | 存在周期/a | 主周期/a |
---|---|---|---|---|
TX10p | 1961~1982, 1990~2006, 1995~2010 | 2, 4, 6 | 2, 6, 20 | 6 |
TN10p | 1967~1977 | 2 | 20, 32 | 无 |
TX90p | 1990~2018, 2000~2008 | 4, 2 | 4, 10, 20, 32 | 4 |
TN90p | 2004~2008 | 2 | 20, 32 | 无 |
SU35 | 1961~1975, 1985~1995, 2006~2016 | 4, 2, 2 | 2, 6, 12, 32 | 2 |
FD0 | 1967~1972, 2003~2011 | 1, 3 | 20 | 无 |
WSDI | 1961~1974, 1990~2021 | 4, 4 | 4, 6 | 4 |
CSDI | 1961~1996 | 4 | 4, 16, 32 | 4 |
表3. 宿迁地区其他极端气温指数的小波能量谱时频特征
近61 a来宿迁地区极端降水强度指数(RX1day、RX5day、R95p、SDII)均呈上升趋势(图3(a)~(d)),上升速率分别为3.14 mm/10 a、5.09 mm/10 a、9.53 mm/10 a和0.3mm∙d−1/10 a,均未通过0.05显著性水平检验,上升趋势不显著;表征极端降水频率的暴雨日数R50,近61 a呈缓慢上升趋势(图3(e)) (P < 0.05),气候倾向率为0.079 d/10 a;表征降水持续时间的两个指数近61 a来呈下降趋势(图3(f)~(g)),趋势变化不明显(P < 0.05),但年际间波动较大,连续无降水日数(CDD)全市平均在39.97 d,连续降水日数(CWD)平均在5.22天。综上,表征极端降水强度和频率的指数呈上升趋势,持续时间指数呈下降趋势,说明宿迁地区降水集中度提高。
图3. 近61年来(1961~2021)宿迁极端降水指数变化趋势
以RX1day为例,分析宿迁市极端降水指数的周期变化特征,结果见图4。由图可知,RX1day在95%的置信区间下有1个聚集中心,其中心坐标为(3, 1985),影响范围为1961~2012年。从小波谱中可以看出,RX1day存在2 a、5 a和10 a左右的周期,2 a和5 a左右的振荡周期通过0.05显著性水平检验,且它们的震荡强度相当,因此RX1day的主周期为2 a和5 a。其他极端降水指数的小波能量谱特征详见表4。
图4. 近61年来(1961~2021)宿迁极端降水强度指数小波能量谱分布
指数 | 集中影响年份 | 尺度中心/a | 存在周期/a | 主周期/a |
---|---|---|---|---|
RX1day | 1961~2012 | 3 | 2, 5, 10 | 2, 5 |
RX5day | 1961~2015 | 2 | 2, 5, 8 | 2 |
R95p | 1961~1990, 1995~2015 | 4, 2 | 2, 5 | 2 |
SDII | 1961~1988, 1994~2015 | 4, 2 | 2, 5, 14 | 2 |
R50 | 1961~1990, 1995~2015 | 4, 2 | 2, 5 | 2 |
CWD | 1961~1978, 1982~1990, 1995~2021 | 6, 3, 7 | 2, 5, 7 | 7 |
CDD | 1961~1998, 2005~2021 | 2, 3 | 2, 20 | 2 |
表4. 宿迁地区其他极端降水指数的小波能量谱时频特征
应用R/S分析法计算了极端气温指数的Hurst指数(见表5),结果表明,宿迁地区各极端气温指数的Hurst指数均大于0.5,表现为各气温序列具有长程相关性,即未来一段时间范围内,宿迁地区各极端气温指数的变化情况与过去61 a的趋势相同。其中,DTR、TN10p、TX90p、TN90p、FD0的H值为>0.8,表现出了很强的持续性,TXx、TNn、TX10p、SU35、CSDI的H值为0.7~0.8,表现出强或较强的持续性,只有WSDI的H值接近0.5,持续性特征很弱。结合前述分析,可知极端气温持续上升,极端高温事件的强度(TXx、TNn)、频率(TN90p、TX90p、SU35)和持续性(WSDI)持续增大,极端低温事件的频率(TN10p、TX10p、FD0)和持续性(CSDI)持续减小。
极端气温指标 | Hurst指数 | |
---|---|---|
强度指标 | TXx | 0.701 (较强) |
TNn | 0.719 (较强) | |
DTR | 0.916 (很强) | |
频率指标 | TX10p | 0.779 (强) |
TN10p | 0.960 (很强) | |
TX90p | 0.894 (很强) | |
TN90P | 0.991 (很强) | |
SU35 | 0.762 (强) | |
FD0 | 1 (很强) | |
持续时间指标 | WSDI | 0.549 (很弱) |
CSDI | 0.732 (较强) |
表5. 宿迁地区极端气温指标的R/S分析
宿迁地区极端降水指数的Hurst指数如表6所示,RX1day、RX5day、CDD、CWD呈现反持续性,R95p、R50和SDII表现为持续性,结合趋势图可知宿迁地区的R95p、SDII、R50、CWD、CDD未来呈增大趋势,RX1day、RX5day呈减小趋势,但是与极端气温指数相比,极端降水指数在持续性和反持续性方面表现出较弱或很弱的特征。
极端降水指标 | Hurst指数 | |
---|---|---|
强度指标 | RX1day | 0.362 (较弱) |
RX5day | 0.432 (较弱) | |
R95p | 0.599 (较弱) | |
SDII | 0.748 (较强) | |
频率指标 | R50 | 0.618 (较弱) |
持续时间指标 | CWD | 0.428 (较弱) |
CDD | 0.482 (很弱) |
表6. 宿迁地区极端降水指标的R/S分析
Hurst指数仅能揭示时间序列的未来趋势,无法确定这种趋势的持续时间,为此,我们使用统计量V来计算平均循环长度,以反映过去趋势对未来的影响时长。本文以气温日较差指数(DTR)为例,图5中,箭头所示的首个明显转折点,即代表系统的平均循环长度。对于气温日较差指数,其平均循环长度为16 a,意味着在16 a后,该指数将不再依赖初始数据。其他极端气候指数的平均循环长度详见表7。从表5~7中我们可以看到,对于同一指数,其持续性特征与平均循环长度之间存在密切关联,当一个指数的持续性越强时,其平均循环长度也越长。整体上,极端气温指数的平均循环长度多大于15 a,而极端降水指数的循环长度多小于10 a,极端气温指数的平均循环长度普遍大于极端降雨指数,这意味着极端气温指数对未来趋势的影响时间更为持久。
图5. 气温日较差(DTR)的V-ln n曲线
极端气候指标 | 非周期循环长度(年) | ||
---|---|---|---|
极端气温指数 | 强度指标 | TXx | 9 |
TNn | 13 | ||
DTR | 16 | ||
频率指标 | TX10p | 15 | |
TN10p | 17 | ||
TX90p | 13 | ||
TN90P | 15 | ||
SU35 | 15 | ||
FD0 | 22 | ||
持续时间指标 | WSDI | 5 | |
CSDI | 15 | ||
极端降水指数 | 强度指标 | RX1day | 6 |
RX5day | 4 | ||
R95p | 6 | ||
SDII | 11 | ||
频率指标 | R50 | 10 | |
持续时间指标 | CWD | 4 | |
CDD | 5 |
表7. 宿迁地区极端气候要素非周期循环长度
在气候变化作用下,宿迁地区极端气候事件加剧。本研究基于实测气象站点数据资料探讨分析了1961~2021年间宿迁地区极端气候历史变化趋势和演变特征,预测了极端气候事件的未来变化趋势,以期为宿迁地区的防灾减灾工作提供帮助。主要结论如下:
1) 近61 a来,宿迁地区极端高温和极端低温均呈增加趋势,极端低温升温幅度大于极端高温,同时气温日较差呈减小趋势。对于极端高温事件,无论是频率指标,还是持续性指,标均呈上升趋势,而极端低温事件的频率指标和持续性指标均呈下降趋势。此外,极端降水事件的强度和频率呈微弱的上升趋势,持续性呈微弱的下降趋势。
2) 冷暖指数和昼夜指数在变化幅度上存在明显差异和不对称性。具体来说,冷指数的升温趋势明显高于暖指数,而夜指数的升温幅度也超过昼指数。
3) 从周期性来看,极端气温指数多存在2~6 a、20 a和32 a的周期变化,但只有2~6 a的周期通过了95%显著性水平检验;极端降水指数多存在2 a和5 a的周期变化,其中2 a左右的周期通过了95%显著性水平检验。
4) 将极端气候指数的历史变化趋势与Hurst指数叠加,可预测未来宿迁地区极端气温持续上升,极端高温(低温)事件的强度、频率和持续性持续增大(减小),宿迁地区气候变暖的状态仍继续保持;极端降水事件中R95p、SDII、R50、CWD、CDD未来呈增大趋势,RX1day、RX5day呈减小趋势。极端气温指数的平均循环长度多大于15 a,而极端降水指数的循环长度多小于10 a。
本文受到宿迁市气象局科研项目(SQ202205)的资助。
徐 恩,周 航,程 昕,宋柳贤,丁晓敏. 宿迁市极端气候事件变化特征及其未来趋势预测Variation Characteristics and Future Trend Prediction of Extreme Climate Events in Suqian[J]. 气候变化研究快报, 2024, 13(02): 380-390. https://doi.org/10.12677/ccrl.2024.132039
https://doi.org/10.1016/j.atmosres.2012.10.022
https://doi.org/10.1016/j.wace.2014.05.001
https://doi.org/10.1016/j.scitotenv.2018.08.262
https://doi.org/10.1007/s00382-009-0735-0
https://doi.org/10.1016/j.scitotenv.2013.07.050
https://doi.org/10.1175/1520-0477(1998)079<0061:APGTWA>2.0.CO;2
https://doi.org/10.13869/j.cnki.rswc.2019.02.039
https://doi.org/10.13869/j.cnki.rswc.2021.01.027
https://doi.org/10.13866/j.azr.2013.02.021