以滇东岩溶地区3个典型气象站1970~2020年逐日降水数据为基础,选取4个极端降水指标,并采用Sen’s斜率法分析、CV值分析、Mann-Kendall突变检验等方法,对比分析了滇东岩溶地区的极端降水事件特征。结果表明:1) 会泽站RX1day呈显著上升趋势且幅度最大,泸西站RX5day和CWD显著下降,3个站点CWD均下降,会泽站和泸西站的CDD均呈上升趋势,砚山站趋势不明显;2) 3个站点四个指标的CV值变化范围分别在0.30~0.35、0.29~0.40、0.30~0.45和0.25~0.35之间,泸西站波动最大,砚山站波动最小;3) Mann-Kendall突变检验中砚山站突变点较多且集中,这是砚山站点与另外两个站点相比最显著的区别,表明砚山站极端降水事件发生的概率较高。综上,总体上降水减少、气候变干旱,所有分析结果均符合气候变化大背景影响下的极端降水时空分布变化趋势。 Based on the daily precipitation data of three typical meteorological stations in the Karst area of eastern Yunnan from 1970 to 2020, four extreme precipitation indexes are selected, and the characteristics of extreme precipitation events are compared and analyzed by using the Sen’s slope and CV value analysis methods as well as Mann Kendall mutation test. The results show that: 1) RX1day of Huize station shows a significant upward trend and the amplitude is the largest, RX5day and CWD of Luxi station decrease significantly, CWD of three stations decreases, CDD of Huize station and Luxi station shows an upward trend, and Yanshan station has no upward trend; 2) The variation range of CV values of the four indicators of the three stations is 0.30~0.35, 0.29~0.40, 0.30~0.45 and 0.25~0.35 respectively. The fluctuation of Luxi station is the largest and that of Yanshan station is the smallest; 3) In Mann Kendall mutation test, there are many concentrated mutation points at Yanshan station, which is the most significant difference between Yanshan station and the other two stations, indicating that the probability of extreme precipitation events at Yanshan station is high. In conclusion, in general, the precipitation decreases and the climate become arid. All the analysis results are in line with the temporal and spatial distribution trend of extreme precipitation under the influence of climate change.
以滇东岩溶地区3个典型气象站1970~2020年逐日降水数据为基础,选取4个极端降水指标,并采用Sen’s斜率法分析、CV值分析、Mann-Kendall突变检验等方法,对比分析了滇东岩溶地区的极端降水事件特征。结果表明:1) 会泽站RX1day呈显著上升趋势且幅度最大,泸西站RX5day和CWD显著下降,3个站点CWD均下降,会泽站和泸西站的CDD均呈上升趋势,砚山站趋势不明显;2) 3个站点四个指标的CV值变化范围分别在0.30~0.35、0.29~0.40、0.30~0.45和0.25~0.35之间,泸西站波动最大,砚山站波动最小;3) Mann-Kendall突变检验中砚山站突变点较多且集中,这是砚山站点与另外两个站点相比最显著的区别,表明砚山站极端降水事件发生的概率较高。综上,总体上降水减少、气候变干旱,所有分析结果均符合气候变化大背景影响下的极端降水时空分布变化趋势。
岩溶石漠化,极端降水,气候变化,滇东南
Jingran Shuai, Wenrong Ding*
Faculty of Geograghy, Yunnan Normal University, Kunming Yunnan
Received: May 12th, 2022; accepted: Jun. 15th, 2022; published: Aug. 31st, 2022
Based on the daily precipitation data of three typical meteorological stations in the Karst area of eastern Yunnan from 1970 to 2020, four extreme precipitation indexes are selected, and the characteristics of extreme precipitation events are compared and analyzed by using the Sen’s slope and CV value analysis methods as well as Mann Kendall mutation test. The results show that: 1) RX1day of Huize station shows a significant upward trend and the amplitude is the largest, RX5day and CWD of Luxi station decrease significantly, CWD of three stations decreases, CDD of Huize station and Luxi station shows an upward trend, and Yanshan station has no upward trend; 2) The variation range of CV values of the four indicators of the three stations is 0.30~0.35, 0.29~0.40, 0.30~0.45 and 0.25~0.35 respectively. The fluctuation of Luxi station is the largest and that of Yanshan station is the smallest; 3) In Mann Kendall mutation test, there are many concentrated mutation points at Yanshan station, which is the most significant difference between Yanshan station and the other two stations, indicating that the probability of extreme precipitation events at Yanshan station is high. In conclusion, in general, the precipitation decreases and the climate become arid. All the analysis results are in line with the temporal and spatial distribution trend of extreme precipitation under the influence of climate change.
Keywords:Karst Rocky Desertification, Extreme Precipitation, Climate Change, Southeast Yunnan
Copyright © 2022 by author(s) and Wuhan University.
This work is licensed under the Creative Commons Attribution International License (CC BY 4.0).
http://creativecommons.org/licenses/by/4.0/
随着社会经济的不断发展,人类活动时刻影响着全球气候的变化,CO2排放量持续增加,气温不断上升,气候变化正在加剧水循环,这会带来更强的降雨和洪水。在许多地区这意味着更严重的干旱 [
岩溶区极端降水事件时空分布变化的研究是当前地学界的一个难点。极端降水事件会造成岩溶区域水土流失严重从而进一步加剧区域气候变化。自20世纪80年代以来,国内外学者 [
因此,为解释滇东岩溶区近50年来极端降水的历史变化特征并服务现今的水资源管理,本文利用云南省滇东岩溶区3个典型站点1970~2020年逐日降水观测数据,选取4个极端降水指标来分析该区域降水事件变化规律,为云南岩溶区水资源规划、防灾减灾,推动区域协调发展等提供科技支撑。
云南省129个县(市、区)总面积为39.4万km2,其中岩溶区119个,总面积为11.1万km2,占云南省总面积的28.2%。滇东地区有51个县(市、区)的岩溶面积超过云南省省总面积的30%,岩溶区面积67,779 km2,占国土总面积的53.9%。岩溶山区地表干旱缺水严重 [
会泽、泸西、砚山站点坐落于滇东岩溶片区,该区位于23˚36'21''N~26˚25'4''N和103˚17'50''E~104˚20'13''E之间,北与滇东北高原相连,南与文山壮族苗族自治州接壤,东临广西,西为滇中高原。整体地势较为平缓,北高南低,自北向南依次为温带高原季风气候、北亚热带季风气候、中亚热带季风气候,高温高湿,雨热同期。年均温12.7℃~16.1℃,降雨量为817~1008 mm。地带性土壤为棕壤或红壤。
图1. 研究区及气象站点
本研究选择的典型站点会泽、泸西、砚山位于云南省东部(图1),属岩溶石漠化区,石漠化比例分别为4.66%、22.45%、37.65%,三个站点的海拔分别为2110.5 m、1704.3 m和1561.1 m (表1)。
站点 | 经度 | 纬度 | 海拔/m | 植被NDVI* | 石漠化比例% |
---|---|---|---|---|---|
会泽 | 103˚17' | 26˚25' | 2110.5 | 0.405 | 4.66 |
泸西 | 103˚46' | 24˚32' | 1704.3 | 0.528 | 22.45 |
砚山 | 104˚20' | 23˚37' | 1561.1 | 0.351 | 37.65 |
表1. 气象站地理位置及周边环境
注:*以气象站为中心的半径3 km内植被NDVI平均值。
本文数据来源于云南省气象局资料中心和中国气象局国家气象信息中心(http://data.cma.cn),选取了会泽、泸西、砚山3个滇东片区典型气象站点的逐日降水量数据资料进行分析,这些资料均经过严格的质量控制,能满足研究的要求。
采用气候变化和指数专家团法(以下简称ETCCI)、Mann-Kendall突变检验法和Sen’s斜率分析法,结合变差系数CV来分析该区域极端降水变化规律。表2为从ETCCI法中选取的最长持续干旱时间(以下简称CDD)、最长持续降水日数(以下简称CWD)、日最大降水量(以下简称RX1day)和连续5日最大降水(以下简称RX5day) 4个极端降水指数。降水数据的计算基于Matlab完成。参考Abdila W. P.,Gabriela-Victoria Harpa,Shawn M. Milrad等 [
指数名称 | 缩写名称与单位 | 定义 |
---|---|---|
最长持续干旱时间 | CDD/d | 日降水量小于1 mm的最大持续天数 |
最长持续降水日数 | CWD/d | 日降水量大于或等于1 mm的最大持续天数续日数 |
日最大降水量 | RX1day/mm | 每年最大1日的降水量 |
连续5日最大降水量 | RX5day/mm | 年内连续5日的最大降水量 |
表2. 选取的各极端降水指数及其定义
1) M-K检验法。曼(H.B. Mann)和肯德尔(M.G. Kendall)提出了原理并发展了这一方法。它可以客观反映时间变化趋势,不受个别异常值干扰。具体计算方法及参数如下 [
构造时间序列X,具n个样本量:
S k = ∑ i = 0 k r i , r i = { 1 , x i > x j j = 1 , 2 , ⋯ , i 0 , e l s e } (1)
秩序列Sk是第i时刻数值大于j时刻数值个数的累计数。时间序列随机独立的情况下,统计量UFk为:
U F k = S k − E ( S k ) V a r ( s k ) , k = 1 , 2 , ⋯ , n (2)
式中:E(Sk)为均值、Var(Sk)为方差。
X 1 , X 2 , ⋯ , X n 相互独立,相同连续分布,计算公式如下:
E ( S k ) = n ( n − 1 ) 4 , V a r ( S k ) = n ( n − 1 ) ( 2 n + 5 ) 72 (3)
UF和UB是按照标准正态分布的算法,将时间序列X分布按照顺序和逆序排列之后得到的统计量。若UF值 > 0,呈上升趋势,反之呈下降趋势。当UF超过的临界值±1.96时,则通过了95%的置信度检验,呈显著上升或下降趋势,超过临界线的范围即为突变时间区域。若UF和UB在临界线之间出现交点,交点所在时刻即为突变时刻。该点即为突变点 [
2) Sen’s斜率法。该方法可以减少噪声对序列变化趋势的干扰。对于时间序列 X ( X 1 , X 2 , ⋯ , X n ) ,计算公式如下 [
β = median ( x j − x i j − i ) , ∀ j > I (4)
式中:Median为中值函数。当β > 0时,呈上升趋势;当β = 0时,序列趋势不明显;当β < 0时,呈下降趋势。
3) 变差系数CV。它可以表示序列总体的离散程度,计算公式为 [
C V = S M × 100 % (5)
式中:S为标准差,M为平均数。
表3为滇东岩溶地区会泽站、泸西站和砚山站三个典型站点的Sen’s斜率结果,可以发现:3个站点的RX5day、RX1day、CWD和CDD四个指数中,会泽站的RX1day、泸西站的RX5day和CWD通过了置信水平为0.05的显著性检验,分别呈显著上升和下降趋势。会泽站的RX1day增加幅度最大,达4.82 mm/10a,泸西站的RX5day减少幅度最大,达−9.31 mm/10a。
Sen’s | RX5day | RX1day | CWD | CDD |
---|---|---|---|---|
会泽 | 0.02 | 0.48* | −0.05 | 0.11 |
泸西 | −0.93* | −0.41 | −0.07* | 0.11 |
砚山 | −0.25 | 0.15 | −0.06 | 0.00 |
表3. 滇东岩溶地区3个典型站点Sen’s斜率
注:*为达到0.05置信水平。
从CWD指数来看,3个站点的最长降水时间均呈下降趋势,其多年平均值分别为−0.48 d/10 a、−0.69 d/10 a和−0.56 d/10 a。从CDD指数来看,会泽站和泸西站的最长干旱时间均呈上升趋势,其多年平均值分别为1.143 d/10 a和1.071 d/10 a,而砚山站却无上升或下降趋势。总体上降水减少,气候变干旱,符合气候变化大环境下的趋势。
表4为滇东岩溶地区会泽站、泸西站和砚山站三个典型站点的不均匀系数结果,可以发现:3个站点的RX5day、RX1day、CWD和CDD四个指数的不均匀系数CV变化范围分别在0.30~0.35、0.29~0.40、0.30~0.45和0.25~0.35之间,其中泸西站CWD指数的不均匀系数最大,达0.432,砚山站的CDD指数的不均匀系数最小,达0.2524。
CV | RX5day | RX1day | CWD | CDD |
---|---|---|---|---|
会泽 | 0.34 | 0.29 | 0.31 | 0.34 |
泸西 | 0.35 | 0.39 | 0.43 | 0.29 |
砚山 | 0.30 | 0.33 | 0.41 | 0.25 |
表4. 滇东岩溶地区3个典型站点CV值
本研究选取滇东岩溶地区3个典型站点会泽站、泸西站和砚山站,利用其1970~2020年逐日降水数据的RX5day、RX1day、CWD和CDD四个指标进行Mann-Kendall突变检验,检验结果显示,发现3个站点的RX5day指标整体上会泽站呈先上升后下降趋势、泸西站呈下降趋势、砚山站呈先下降后上升趋势;RX1day指标会泽站整体呈上升趋势、泸西站呈先上升后下降趋势、砚山站呈先上升后平稳趋势;CWD指标会泽站呈先上升后下降趋势,泸西站和砚山站呈下降趋势;CDD指标会泽站和泸西站呈先下降后上升趋势,砚山站呈先上升后下降趋势。
对比会泽站和泸西站,砚山站的RX5day指数于1983年之后有波动但整体上呈平稳趋势,其突变点较多且比较集中,这是砚山站点与另外两个站点相比最显著的区别,表明砚山站极端气候发生的概率较高。会泽站只有RX1day在1994有一个突变点、泸西站只有RX5day和CWD仅出现一个突变点,分别在1985年和1972年发生突变,这种情况可能是由于大背景下的气候条件的影响。
图2为识别检验会泽站1970~2020年RX5day、RX1day、CWD和CDD四个指标突变检验结果,发现RX5day的UF和UB两条曲线于1976年、1988年、1991年、1994年和2015年附近出现交点,且交点在临界线之间,表明会泽站RX5day序列在上述年份可能发生突变,其中突变年份多发生在1988~1994年间,在此期间会泽站RX5day呈下降趋势,在1976年和2015年发生突变后呈上升趋势。如图2(b)所示,RX1day的UF和UB曲线在1994年附近相交,且交点在临界线之间,即在1994年左右发生突变,突变后呈上升趋势。如图2(c)所示,CWD的UF和UB两条曲线在置信区间内出现多个交点,交点分别发生在1971年、1978年、1983年、1988年和1992年附近,说明会泽站CWD序列在上述年份可能发生突变,其中1988~1994年发生突变比较频繁,1994年突变后会泽站CWD呈下降趋势。如图2(d)所示,CDD的UF和UB曲线在2005年和2019年相交,且交点在临界线之间,说明会泽站CDD序列在2005年和2019年发生突变,突变后会泽站CDD呈上升趋势。
图3为泸西站1970~2020年RX5day、RX1day、CWD和CDD四个指标突变检验结果,发现:RX5day、RX1day、CWD的UF和UB两条曲线分别在临界线之间在1985年、1975年和1985年、1972年附近出现交点,表明泸西站RX5day、RX1day、CWD序列在以上各自所述年份前后发生了突变,均在突变后呈下降趋势。此外,CWD的UF曲线在1975年以后超出置信上限,呈显著下降趋势,1980年后又回复到临界线之间;最后UF曲线在2003年以后又超出置信上限,呈显著下降趋势。如图3(d)所示,CDD的两条曲线在置信区间内出现多个交点,交点分别发生在1975年、1977年、1992年和1998年,表明泸西站CDD序列在1975年、1977年、1992年和1998年可能发生突变,其中1998年突变后泸西站CDD呈上升趋势,其他突变点突变后CDD呈下降趋势。
图4为砚山站1970~2020年RX5day、RX1day、CWD和CDD四个指标突变检验结果,发现:RX5dayUF和UB两条曲线于1973年、1980年、1983年、1988年、1991年、1993年、2005年和2007年附近出现交点,且交点在临界线之间,表明砚山站RX5day序列在上述年份均发生了突变,其中1980~1993年频繁交替出现交点,是由于1980~1993年降水年际间变化较大。如图4(b)所示,RX1day的UF和UB曲线在置信区间内相交,交点发生在1982年、1984年、1992年、1994年、1996年、2003年、2009年、2011年、2015年和2018年附近,说明砚山站RX1day序列在以上年份可能发生突变,其中1992年、1994年和2015年突变后呈下降趋势,其他年份突变后呈上升趋势。如图4(c)所示,CWD的UF和UB曲线在置信区间内相交,交点发生在1972年、1974年、1979~1984年,说明砚山站CWD序列在上述年份可能发生突变,突变后呈下降趋势。此外,UF曲线
图2. 会泽站极端降水指数的Mann-Kendall突变检验
在1993年以后超出置信上限,呈显著下降趋势,1980年后又回复到临界线之间;最后UF曲线在2003年以后又超出置信上限,表明该地区的CWD序列呈显著下降趋势。如图4(d)所示,CDD的两条曲线在置信区间内出现多个交点,交点分别发生在1972年、2000年、2002年、2008年、2011年、2013年、2015年、2017年和2019年,交点在临界线之间,表明砚山站CDD序列在上述年份可能发生突变,其中1972年、2008年突变后砚山站CDD呈上升趋势,其他突变点突变后CDD呈下降趋势。
本文基于滇东岩溶区3个典型站点1970~2017年的逐日降水数据,计算4个极端降水指标,分析其降水变化趋势,得到以下结论:
图3. 泸西站极端降水指数的Mann-Kendall突变检验
1) 对Sen’s斜率的分析,3个站点的RX5day、RX1day、CWD和CDD四个指数中,会泽站的RX1day、泸西站的RX5day和CWD通过了置信水平为0.05的显著性检验,分别呈显著上升和下降趋势。会泽站的RX1day增加幅度最大。从CWD指数来看,3个站点的最长降水时间均呈下降趋势。从CDD指数来看,会泽站和泸西站的最长干旱时间均呈上升趋势,而砚山站却无上升或下降趋势。总体上降水减少,气候变干旱,符合气候变化大环境下的趋势。
2) 对CV值的分析,3个站点的RX5day、RX1day、CWD和CDD四个指数的不均匀系数CV变化范围分别在0.30~0.35、0.29~0.40、0.30~0.45和0.25~0.35之间,其中泸西站CWD指数的不均匀系数最大,达0.43,砚山站的CDD指数的不均匀系数最小,达0.25。
3) 对Mann-Kendall突变检验的结果分析,3个站点的RX5day指标整体上会泽站呈先上升后下降趋势、泸西站呈下降趋势、砚山站呈先下降后上升趋势;RX1day指标会泽站整体呈上升趋势、泸西站呈先上升后下降趋势、砚山站呈先上升后平稳趋势;CWD指标会泽站呈先上升后下降趋势,泸西站和砚山站呈下降趋势。对比会泽站和泸西站,砚山站的RX5day指数于1983年之后有波动但整体上呈平稳趋势,其突变点较多且比较集中,这是砚山站点与另外两个站点相比最显著的区别,表明砚山站极端降水事件发生的概率较高。会泽站只有RX1day在1994有一个突变点、泸西站只有RX5day和CWD仅出现一个突变点,分别在1985年和1972年发生突变,均符合气候变化大背景影响下的极端降水时空分布变化趋势。
图4. 砚山站极端降水指数的Mann-Kendall突变检验
气候变化在宏观上主要表现为全球大气变暖和降水频率变化,当前分析极端降水时空变化对区域的影响已成为重点研究内容。全球增温和人类活动的剧烈作用导致全球极端降水气候事件日益频发,已危害到人类社会经济健康发展。极端降水气候事件自1950年以来愈发频繁 [
全球气候正经历着前所未有的变化,根据现有的评估可以发现,目前在复合型极端事件发生发展机理认识方面还存在不足。同时,未来仍需进一步完善跨学科跨部门跨区域研究,加强对复合型极端事件形成机理、预估及其对生态系统,经济社会影响风险的评估,提高对区域气候变化的适应能力。IPCC最新发布的《气候变化与土地特别报告》(SRCCL),对岩溶区荒漠化展开了探讨。发现随着全球气温不断变暖,岩溶区荒漠化和土地退化威胁着粮食安全。气候不断变化,极端事件发生的概率和严重程度会随之提高。面对气候变化尤其是极端事件给土地带来的巨大压力,我们必须采取相关措施来减轻她所造成的负面影响。如:坚持土地管理可持续性,减少所有行业的CO2等温室气体排放等,才有可能到21世纪末将全球平均升温控制在相对工业化前水平2℃以内。由此减轻气候变化对土地和粮食系统的负面影响 [
国家自然科学基金项目(41761105)资助。
帅静然,丁文荣. 近50年来滇东岩溶地区极端降水变化分析Analysis on Temporal Variation of Extreme Precipitation in Karst Area of Eastern Yunnan in Recent 50 Years[J]. 水资源研究, 2022, 11(04): 346-356. https://doi.org/10.12677/JWRR.2022.114038