未来气候模式下黄河流域极端降水指数时空分布特征

2024-12-31 00:00:00李小雨张国栋尹昌燕张雪丽
人民黄河 2024年11期
关键词:黄河流域

摘 要:以气候变化敏感区黄河流域为研究区,再现和预估流域内极端降水事件的时空分布特征。采用Delta 降尺度以及多指标评价方法建立CMIP6 气候模式数据集,选取12 个极端降水指数研究极端降水指数的时空分布特征。结果表明,黄河流域历史(1980—2014 年)及SSP126、SSP245、SSP585 情景下未来(2022—2100 年)除连续干日数(CDD)外的11 个极端降水指数多年均值的空间格局相似,整体上呈现出由东南向西北递减的趋势,南北向变化率比东西向的大。年际上,11 个极端降水指数均呈现不显著的波动性上升趋势。整体上各极端降水指数呈现波动性上升趋势,但上升趋势并不显著。

关键词:未来气候模式;极端降水指数;气候情景;黄河流域

中图分类号:P338;TV882.1 文献标志码:A doi:10.3969/ j.issn.1000-1379,2024.11.006

引用格式:李小雨,张国栋,尹昌燕,等.未来气候模式下黄河流域极端降水指数时空分布特征[J].人民黄河,2024,46(11):37-42.

0 引言

联合国气候变化专门委员会(IPCC)第六次评估报告(AR6)指出,全球变暖趋势持续加剧,区域气候变化将更明显[1-6] 。其中,极端降水与气候变暖间的关系是复杂的,温度升高,大气在饱和前储存更多的水分,直接导致短时间内的降水量增多。研究表明,气候变化对极端降水的影响主要表现在改变其频率、强度、持续时间、空间分布等特性上,从而导致洪水、水土流失、泥石流等次生灾害明显增加[7] 。因此,科学认识气候驱动的极端降水变化机理,对于流域生态保护以及水、粮食和能源安全十分重要。

目前对于极端降水事件的评估,最常用的是气候变化监测专家组推荐的极端降水指数。国内外常用的极端降水指数主要包括:Frich 等[8] 提出的极端降水指数,包括持续干日数(CDD)、强降水量(R90P)、极强降水量(R95P)、5 d 最大降水量(RX5Day);IPCC 第四次报告确定的27 个核心极端气候指数[9] ,其中与降水相关的指数有12 个(见表1)。

根据指数所描述降雨特性,所有指数可以分为3类:1)持续日数,反映某一状态持续时间,包括CDD、CWD;2)降水量,反映某一时段内所累计的降水量的大小, 包括PRCPTOT、R90P、R95P、RX1Day 和RX5Day;3)降水强度,反映单位时间内的降水量,包括SDII、SDII(90)、SDII(95)。

从20 世纪90 年代中期开始,气候建模工作主要与IPCC 的耦合模型相互比较项目(CMIP)一起进行。相较于CMIP5,第六阶段耦合比较计划(CMIP6)设计了更多的数值试验,提供了更庞大的模拟数据集,且为了解决分辨率较低的问题,CMIP6 推出了高分辨率区域气候模式( WRF), 但目前所参与的模式并不多[10-11] 。不同于CMIP5,CMIP6 基于不同共享社会经济路径可能发生的能源结构所产生的人为排放及土地利用变化,设计了一系列新的情景预估试验,为未来气候变化机理研究以及气候变化减缓和适应研究提供关键的数据支持。舒章康等[12] 基于CMIP6 的11 种全球气候模式,分析并预测了我国未来极端降水和高温事件的演变特征,结果显示,相较于历史时期,未来我国各地区极端降水和高温事件均呈现上升趋势。Faye等[13] 利用19 个CMIP6 模式模拟西非地区9 个极端降水指数,结果表明CMIP6 模式能够完整地再现研究区极端降水指数时空分布。气候模式对于极端降水的模拟性能存在明显的地区差异,且极端降水相较于极端高温事件,其变化趋势表现出更强的区域差异[14] 。因此,针对不同区域开展深入研究,对于模型的改进具有重要意义。

本研究以黄河流域为研究区,基于优选评估后和多模式集合后的CMIP6 数据集,提取12 个极端降水指数,再现历史及预估未来各情景下黄河流域极端降水事件的时间演变和空间分布特征。

1 研究区概况与研究方法

1.1 研究区概况

黄河流域处于干旱半干旱区,属于气候敏感区,气温升高导致流域内汛期常见大范围暴雨、持续性降水和对流性大风等天气。黄河流域多年平均降水量约为476 mm,主要以暴雨形式降落。流域内降水具有明显的季节性,全年约70%降水主要集中于夏季且呈现出强烈的区域性,而冬季降水少且变化率小。空间上,年降水整体上呈从东南向西北递减趋势,降水高值区主要分布于流域东南部湿润半湿润地区,年降水量超800 mm,低值区主要分布于流域北部干旱地区,年降水量小于200 mm。年际差异大,丰水年与枯水年的降水量之比为1.7~7.5。

1.2 数据的收集与处理

1)水文数据。本文所采用的水文数据主要包括:黄河流域1980—2014 年94 个地面气象站逐日降水量数据,来源于国家气候中心(http:// www.nmic.cn/ );黄河流域1980—2018 年水平分辨率为31 km 的ERA5降水再分析数据,其作为欧洲中期天气预报中心全球气候最新产品,对全球日降水数据具有一定的再现能力。

2) 气候模式数据。基于前人研究, 本文选取CMIP6 气候模式中对于黄河流域极端降水指数模拟较优且资料完备的6 个模式探究历史及未来时期降水量[15-16] ,包括EC -Earth3、EC -Earth3 - Veg、GFDL -ESM4、MPI-ESM1 -2 -HR、MRI-ESM2 -0 和IPSLCM6A-LR 的日降水数据集(见表2),数据来源于ht⁃tps:// esgf-index1.ceda.ac.uk/ projects-cmip6-ceda/ 。

下载的CMIP6 数据集提供了1980—2014 年历史模拟及2015—2100 年未来预估的日降水数据,为与其保持时间信息一致性,本文以1980—2014 年为历史基准期,以2022—2100 年为未来预估时期,选取3 种不同共享社会经济路径( Shared Socio - economicPathways, SSPs),即低、中、高排放情景下的SSP1-2.6(SSP126)、SSP2-4.5(SSP245)和SSP5-8.5(SSP585)进行对比分析。

3)极端降水指数。根据IPCC 第四次报告所确定的27 个核心极端气候指数,本文选择与降水相关的12 个作为评估指标(见表1),从降水的强度、频次、持续性和贡献率等方面全面评估黄河流域降水特性。其中,SDII(90),SDII(95)可由SDII 定义拓展计算。

1.3 研究方法

1)Delta 降尺度法。为保持数据分辨率的一致性,本文采用双线性插值法将所有历史及未来数据统一插值到0.25°×0.25°网格。双线性插值是对线性插值在二维向量上的拓展,即在两个方向上分别进行一次线性插值[17] 。

为弥补低分辨率数据的不足,本文以ERA5 为参考数据,采用Delta 降尺度法对历史模式及未来数据进行偏差修正。该方法主要通过比较历史模式数据与实测数据间的差异求得变化因子,并假设未来情景中该变化因子保持不变,由此达到修正未来数据的目的[18] 。

2)气候模式评估。本文以历史时期为验证期,根据泰勒图(S)[19] 、时间变率技巧评分(IVS)[20] 和综合评级指数(MR)[21] 对历史时期观测的与降尺度模拟的极端降水指数进行误差检验,优选出适合黄河流域极端降水指数分析的气候模式。

2 结果与分析

2.1 气候模式评估及优选

历史时期(1980—2014 年)6 个气候模式对于12个极端降水指数的空间格局的模拟能力普遍不高,空间相关系数整体分布于0.3~0.7,标准差整体分布于0.5~1.5,中心均方根误差整体分布于0.5~1.0,对不同指数上显现出不同的模拟效果,其中,所有模式对SDII、SDII(90)、SDII(95)和RX5Day 的模拟能力都较强,对PRCPTOT、R90P、R95P、CDD 和CWD 的模拟能力较弱。对于CWD,6 个模式的S 值均小于0.3,其中IPSL-CM6A-LR 模拟效果最差(S 值= 0.07),MPIESM1-2-HR 模拟效果最优(S 值= 0.29)。各模式整体上对R20MM 的模拟效果最优,有6 个模式的S 值均大于0.4。

IVS 得分显示,6 个CMIP6 气候模式对于黄河流域12 个极端降水指数的年际变化率的模拟能力整体一般。从指数上看,各模式对PRCPTOT、R90P、R95P和R20MM 的年际变化率模拟效果都较好,其中对R95P 的模拟最优,6 个模式的IVS 得分均小于0.4,各模式中EC-Earth3-Veg模拟效果最优,MRI-ESM2-0模拟效果最差。各模式对CWD 的年际变化率的模拟效果最差,6 个模式的IVS 得分均大于0.4。

对于空间格局的再现能力,EC-Earth3-Veg、GFDLESM4和IPSL-CM6A-LR 对于12 个指数的整体模拟性能较好;对于年际变化率的再现能力,MPI-ESM1-2-HR和MRI-ESM2-0 对于12 个指数的整体模拟性能较好。本研究中,对于空间格局表现最佳的模型与年际变化率表现最佳的模型不一致,说明本研究所用的CMIP6 模式对于极端降水指数的空间格局和年际变化率的模拟性能不一致,见图1(图中,MME2 为前两名模式平均集合,MME3 为前三名模式平均集合,MME4 为前四名模式平均集合,MME5 为前五名模式平均集合,MME6 为所有模式平均集合)。

为了降低气候模式的不确定性,基于单个模式的诊断结果,比较MME2、MME3、MME4、MME5、MME6的模拟效果,这几种模式的平均集合对于12 个极端降水指数的模拟性能都有明显的提高,相关系数均能达0.6,且标准差、中心均方根误差较观测值均减小。多模式集合后的数据集对于12 个极端降水指数的年际变化率的模拟性能略有提升,使得IVS 得分更接近于0。MR 得分表明多模式集合能显著改善对12 个极端降水指数整体的空间格局和年际变化率的再现能力。

由各模式和多模式集合数据集的MR 平均得分,获得各模式对于流域内极端降水指数的整体模拟性能排名,结果显示,随着集合平均的模式数量的增加,模拟性能先提高后降低,其中多模式集合平均中前四名模式平均集合MME4 对于12 个极端降水指数的模拟最优(见图2)。因此,本文以MME4 作为模式数据集,用于提取黄河流域历史及未来SSP126、SSP245 和SSP585 情景下12 个极端降水指数。

2.2 历史及各情景下未来极端降水指数时空变化趋势

2.2.1 历史及各情景下未来极端降水指数空间格局

历史及未来时期12 个指数的空间分布格局差异都不明显,除CDD 外,整体上呈现从东南向西北递减的趋势。对于描述降雨持续性的指数CWD、CDD(见图3),历史时期流域内多年均值为21.93、43. 44 d, 未来SSP126、SSP245、SSP585 情景下相较于历史时期的变化率分别为-12.1%、- 11.5%、- 16.2%,4.8%、- 1.8%、-3.3%,空间上呈现南多北少的格局。

对于描述降水量的指数RX1Day、RX5Day、PRCP⁃TOT、R90P、R95P, 历史时期流域多年均值分别为18.91、41.52、529.99、157.08、94.94 mm,相较于历史时期,SSP126、SSP245、SSP585 情景下变化率分别为14.3%、13.6%、15. 4%,12. 3%、12. 0%、13. 8%,19. 7%、15.2%、18.2%,23.6%、19.3%、23.0%,23.7%、20.0%、23.7%。

对于描述降水强度的指数SDII、SDII(90)、SDII(95)(见图4),历史时期流域多年均值为3.86、11.63、13.93 mm/ d,相较于历史时期,SSP126、SSP245、SSP585情景下变化率分别为8.7%、7.5%、8.5%,10.8%、10.6%、11.7%,11.4%、11.1%、12.5%。这些指数在空间分布上具有相似的格局,整体上呈现从东南向西北递减的趋势,高值区出现在三门峡至花园口段,低值区出现于河源至头道拐段,最大值与最小值之间差异明显。

2.2.2 历史及各情景下未来极端降水指数年际变化

CDD 变化率最大,其历史及SSP126、SSP245 和SSP585 情景下未来多年均值分别为21.9、23.0、21.5、21.2 d,标准差分别为3.0、3.2、2.8、2.4 d,CDD 上下波动幅度较大,历史时期最大值51.05 d 出现于2005 年,最小值36.02 d 出现于1983 年。对于描述降水量的指数,PRCPTOT 的增长较为明显,其历史及各情景下未来多年均值分别为530.0、634.2、610.7、626.3 mm,标准差分别为33.7、44.2、50.6、49.3 mm,PRCPTOT 历史时期最大值586.66 mm 出现于2013 年,最小值451.09 mm出现于1992 年。对于描述降水强度的指数,其变化率的大小规律为SDII<SDII(90)<SDII(95),即95%分位值下极强降水强度的年际变化率最大,其历史及各情景下未来多年均值分别为13.9、15.5、15.5、15.7 mm/ d,标准差分别为0.6、0.4、0.5、0.6 mm/ d,SDII 历史时期最大值15.65 mm/ d 出现于2007 年,最小值12.71 mm/ d出现于1986 年。对于描述降水日数的指数,R10MM变化率最大,其历史及SSP126、SSP245 和SSP585 情景下未来多年均值分别为7.86、11.40、10.73、11.29 d,标准差分别为1.14、1.57、1.83、1.72 d,R10MM 上下波动幅度较大,历史时期最大值10.27 d 出现于1982 年,最小值5.52 d 出现于2003 年。历史及各情景下未来黄河流域极端降水指数SDII 年际变化如图5 所示。

3 讨论

模式EC - Earth3 - Veg、GFDL - ESM4 和IPSL -CM6A-LR 对于12 个指数的空间分布再现性能相对较好,MPI-ESM1-2-HR 和MRI-ESM2-0 对于12 个指数的年际变化率的整体模拟较好,表明模式对于空间格局和年际变化率的捕捉能力并不一致[22] 。原因可能是黄河流域横跨多个地貌单元,本身的降雨时空分布特征较为复杂,导致模式较难捕捉其降雨时空分布特征。同时,模式对于不同指数的模拟能力不同,整体上看, 模式对于SDII、SDII ( 90)、SDII ( 95) 和RX5Day 的模拟能力较强,而对于PRCPTOT、R90P、R95P、CWD 的模拟能力较弱,表现出对于平均态势指数的模拟能力较强,而对于描述降水量年内分布的指数的模拟能力较差,存在明显的偏高估[15] ,这可能与研究区所处的地理位置有关,CMIP6 气候模式相较于CMIP5 气候模式,其对于干旱、半干旱地区的平均极端降水的模拟能力显著改善,对于湿润区的改善并不明显,其物理过程、动力学框架还需要更深入的研究。

冯安兰等[23] 研究发现黄河流域极端降水指数空间分布格局呈现西北低、东南高的格局,这与本研究的结果是一致的。相较于原始的数据集,MME2、MME3、MME4、MME5、MME6 对于12 个极端降水指数的年际变化率和空间格局的模拟性能均有所提高。综合得分随着模式平均个数的增加先增大后减小,文中前四名模式集合平均生成的数据集(MME4)最能真实地模拟黄河流域的极端降水指数。原因在于多模式集合使得空间格局和年际变化率间的模拟误差相互抵消,可提高模式的模拟能力,但随着模式样本容量的增大,达到一定数量后会趋于稳定[24-27] 。由此,多模式集合中,通过改变集合模式个数获得最优的数据集,而不是利用所有模式的集合,可进一步降低模式的不确定性。

4 结论

基于CMIP6 中的6 个模式数据,采用双线性插值结合Delta 降尺度对模式数据进行偏差校正,并利用多种指数评估气候模式对于黄河流域12 个极端降水指数的模拟性能,根据优选集合后的模式数据分析了流域内历史及各情景下未来极端降水指数的时空特征。MR 综合得分显示,GFDL-ESM4、EC-Earth3-Veg和EC-Earth3 是对黄河流域12 个极端降水指数的空间格局和年际变化率模拟效果较优的3 个模式。多模式集合能显著改善模式对黄河流域12 个极端降水指数的空间格局和年际变化率的再现性能。除CDD 外的11 个指数,整体上呈现出由东南向西北递减的趋势,南北向比东西向变化率更大,降水高值区位于三门峡至花园口段,低值区出现于兰州至河口段,最大值与最小值之间差异明显。年际上,各极端降水指数呈现波动性上升趋势,但上升趋势并不显著。

参考文献:

[1] 姜克隽.IPCC 第三工作组第六次评估报告:全球减缓走向何方? [J].气候变化研究进展,2020,16(2):251-252.

[2] 周天军,邹立维,陈晓龙.第六次国际耦合模式比较计划(CMIP6)评述[J].气候变化研究进展, 2019,15(5):445-456.

[3] LI C,ZWIERS F,ZHANG X B,et al.Larger Increases in More Extreme Local Precipitation Events as Climate Warms[J].Geo⁃physical Research Letters,2019,46(12):6885-6891.

[4] CHEN L,FORD T W.Effects of 0.5 °C Less Global Warming on Climate Extremes in the Contiguous United States[J].Cli⁃mate Dynamics,2021,57:303-319.

[5] 吕睿喆,赵晓祥,翁白莎,等.乌梁素海流域气候模式适用性评价及旱涝预估[J].人民黄河,2017,39(4):53-57.

[6] 刘昌明,刘璇,杨亚锋,等.水文地理研究发展若干问题商榷[J].地理学报,2022,77(1):3-15.

[7] GAO T,WANG H L.Trends in Precipitation Extremes over the Yellow River Basin in North China:Changing Properties and Causes[J].Hydrological Processes,2017,31(13):2412-2428.

[8] FRICH P, ALEXANDER L, DELLA⁃MARTA P, et al.Observed Coherent Changes in Climatic Extremes During the Second Half of the Twentieth Century[J].Climate Research,2002,19(3):193-212.

[9] 尹红,孙颖.基于ETCCDI 指数2017 年中国极端温度和降水特征分析[J].气候变化研究进展,2019,15(4):363-373.

[10] 王闯,戴长雷,宋成杰.青藏高原气候变化的时空分布特征分析[J].人民黄河,2022,44(9):76-82.

[11] 沙健,路瑞,孙运海,等.基于模型联用评估气候变化对径流特征的影响[J].人民黄河,2017,39(11):57-62,65.

[12] 舒章康,李文鑫,张建云,等.中国极端降水和高温历史变化及未来趋势[J].中国工程科学,2022,24(5):116-125.

[13] FAYE A,AKINSANOLA A A.Evaluation of Extreme Pre⁃cipitation Indices over West Africa in CMIP6 Models[J].Climate Dynamics,2021,58:925-939.

[14] 张磊,王春燕,潘小多.基于区域气候模式未来气候变化研究综述[J].高原气象,2018,37(5):1440-1448.

[15] 王予,李惠心,王会军,等.CMIP6 全球气候模式对中国极端降水模拟能力的评估及其与CMIP5 的比较[J].气象学报,2021,79(3):369-386.

[16] XU H W,CHEN H P,WANG H J.Future Changes in Precipi⁃tation Extremes Across China Based on CMIP6 Models[J].In⁃ternational Journal of Climatology,2022,42(1):635-651.

[17] 徐振亚,任福民,杨修群,等.日最高温度统计降尺度方法的比较研究[J].气象科学,2012,32(4):395-402.

[18] 刘倩,高路,马苗苗,等.辽宁大凌河流域气温和降水降尺度研究[J].水利水电技术(中英文),2021,52(9):16-31.

[19] 向竣文,张利平,邓瑶,等.基于CMIP6 的中国主要地区极端气温/ 降水模拟能力评估及未来情景预估[J].武汉大学学报(工学版),2021,54(1):46-57,81.

[20] CHEN W L,JIANG Z H,LI L.Probabilistic Projections of Climate Change over China Under the SRES A1B Scenario Using 28 AOGCMs[J].Journal of Climate,2011,24(17):4741-4756.

[21] JIANG Z H,LI W,XU J J,et al.Extreme Precipitation Indi⁃ces over China in CMIP5 Models.Part I:Model Evaluation [J].Journal of Climate,2015,28(21):8603-8619.

[22] TANG B,HU W T,DUAN A M.Assessment of Extreme Pre⁃cipitation Indices over Indochina and South China in CMIP6 Models[J].Journal of Climate,2021,34(18):7507-7524.

[23] 冯安兰,张强,宋金帛,等.基于CMIP6 的黄河流域极端降水时空特征分析[J].北京师范大学学报(自然科学版),2024,60(2):270-284.

[24] JIAN S Q,ZHU T S,WANG J Y,et al.The Current and Future Potential Geographical Distribution and Evolution Process of Catalpa bungei in China[J]. Forests,2022,13(1):96.

[25] HAMLET A F,BYUN K,ROBESON S M,et al.Impacts of Climate Change on the State of Indiana:Ensemble Future Projections Based on Statistical Downscaling [J]. Climatic Change,2019,163(4):1881-1895.

[26] INNOCENTI S,MAILHOT A,LEDUC M,et al. Projected Changes in the Probability Distributions, Seasonality, and Spatiotemporal Scaling of Daily and Subdaily Extreme Pre⁃cipitation Simulated by a 50⁃Member Ensemble over North⁃eastern North America[J].Journal of Geophysical Research:Atmospheres,2019,124(19):10427-10449.

[27] CHEN J,BRISSETTE F.Stochastic Generation of Daily Pre⁃cipitation Amounts: Review and Evaluation of Different Models[J].Climate Research,2014,59(3):189-206.

【责任编辑 张 帅】

基金项目:黄委优秀青年人才科技项目(HQK-202305);贵州省科技支撑一般项目(206);国家重点研发计划项目(2022YFC300340502)

猜你喜欢
黄河流域
实施积极财政政策 推进黄河流域生态保护和高质量发展
山西财税(2023年8期)2023-12-28 08:05:28
生态环境部启动新一年度黄河流域“清废行动”
黄河流域文学生态意识的地域透视
文学与文化(2022年3期)2022-11-19 02:32:20
“黄河流域生态保护和高度质量发展”
人民黄河(2022年9期)2022-09-08 13:04:38
“黄河流域生态保护和高质量发展”专题征稿
人民黄河(2022年7期)2022-07-07 06:52:24
“黄河流域生态保护和高质量发展”专题征稿
人民黄河(2022年4期)2022-04-07 09:04:04
黄河流域甘肃段生态治理的法治保障分析
在黄河流域生态保护和高质量发展中展现陕西担当
当代陕西(2019年23期)2020-01-06 12:17:40
黄河流域灿烂的齐家文化——陶器
收藏界(2019年2期)2019-10-12 08:26:10
增强大局意识 提升黄河流域生态保护发展水平
人大建设(2019年12期)2019-05-21 02:55:42