• CN 62-1175/P
  • ISSN 1006-7639
  • 双月刊
  • 中国科技核心期刊
  • 中国学术期刊综合评价数据库统计源期刊
  • 中文科技期刊数据库收录期刊

干旱气象, 2026, 44(3): 398-411 DOI: 10.11755/j.issn.1006-7639-2026-03-0398

论文

机器学习在中国区域极端气候指数集合预估中的应用

刘明铭,1,2,3, 徐影,2,3

1 中国气象科学研究院北京 100081

2 中国气象局国家气候中心北京 100081

3 中国气象局气候研究开放实验室北京 100081

Application of machine learning in the ensemble projection of regional extreme climate indices over China

LIU Mingming,1,2,3, XU Ying,2,3

1 Chinese Academy of Meteorological SciencesBeijing 100081China

2 National Climate CenterChina Meteorological AdministrationBeijing 100081China

3 Open Laboratory for Climate StudiesChina Meteorological AdministrationBeijing 100081China

通讯作者: 徐影(1967—),女,研究员,主要从事气候变化未来预估研究。E-mail:xuying@cma.gov.cn

责任编辑: 王涓力;校对:邓祖琴

收稿日期: 2026-01-8   修回日期: 2026-04-5  

基金资助: 西藏自治区重大科技专项(XZ202402ZD0006)
国家科技重大专项(2025ZD1208301)
国家气候中心重点创新团队“第三极气候变化监测预估”(NCCCXTD007)

Received: 2026-01-8   Revised: 2026-04-5  

作者简介 About authors

刘明铭(2000—),男,硕士研究生,主要从事气候变化未来预估研究。E-mail: 2545646576@qq.com

摘要

在气候变暖背景下,极端气候事件的变化趋势备受关注。基于3种机器学习(随机森林、极端随机树和岭回归)模型评估其对中国区域极端气候指数[暖昼指数(TX90p)、冷夜指数(TN10p)、日最大降水量(RX1day)及5 d最大降水量(RX5day)]的模拟能力,并与传统全局偏差订正后的多模式集合方法对比,确定最优模型方案;进一步分析所选取的极端气候指数在不同排放情景(SSP1-2.6、SSP2-4.5、SSP5-8.5)下2024—2100年相对于基准期(1961—1990年)的空间分布及其变化趋势。结果表明:基于机器学习的极端气候指数模拟方案能够在不同程度上提升对极端气候事件的模拟能力,有效减小模拟偏差;2024—2100年中国地区极端暖事件在各排放情景下均显著上升,且高排放(SSP5-8.5)下增幅最大(相较基准期上升约52%),极端冷事件显著减少,且随排放增加减少更明显;空间分布上,TX90p在华北、长江中下游、四川盆地及华南部分地区增幅相对明显,而青藏高原及西北部分高海拔地区增幅相对较小;极端降水指数(RX1day、RX5day)在所有排放情景下均呈增加趋势,且高排放情景下增加最明显,其中RX5day的增强幅度整体高于RX1day;东北和华北地区在极端降水(特别RX5day)上的响应最强,对气候变暖的敏感性更高。

关键词: 机器学习; 极端气候指数; 未来预估

Abstract

Under the background of climate warming, changes in extreme climate events have received increasing attention. In this study, three machine learning models, namely Random Forest, Extremely Randomized Trees, and Ridge Regression, were used to evaluate their capability in simulating extreme climate indices over China, including the warm day index (TX90p), cold night index (TN10p), maximum 1-day precipitation (RX1day), and maximum consecutive 5-day precipitation (RX5day). Their performances were compared with that of the traditional globally bias-corrected multi-model ensemble method to determine the optimal model schemes. Furthermore, the spatial distributions and temporal trends of the selected extreme climate indices during 2024-2100 under different emission scenarios, including SSP1-2.6, SSP2-4.5, and SSP5-8.5, were analyzed relative to the baseline period of 1961-1990. The results show that the machine-learning-based schemes can improve the simulation capability for extreme climate events to varying degrees and effectively reduce simulation biases. During 2024-2100, extremely warm events over China increase significantly under all emission scenarios, with the largest increase under the high-emission scenario SSP5-8.5, reaching approximately 52% relative to the baseline period. In contrast, extremely cold events decrease significantly, and the reduction becomes more pronounced with increasing emissions. Spatially, TX90p shows relatively large increases in North China, the middle and lower reaches of the Yangtze River, the Sichuan Basin, and parts of South China, while the increases are relatively smaller over the Qinghai-Xizang Plateau and some high-elevation areas in Northwest China. Extreme precipitation indices, including RX1day and RX5day, show increasing trends under all emission scenarios, with the most pronounced increase under the high-emission scenario. The enhancement of RX5day is generally stronger than that of RX1day. Northeast China and North China show the strongest responses in extreme precipitation, especially for RX5day, indicating higher sensitivity to climate warming.

Keywords: machine learning; extreme climate indices; future projection

PDF (22997KB) 元数据 多维度评价 相关文章 导出 EndNote| Ris| Bibtex  收藏本文

本文引用格式

刘明铭, 徐影. 机器学习在中国区域极端气候指数集合预估中的应用[J]. 干旱气象, 2026, 44(3): 398-411 DOI:10.11755/j.issn.1006-7639-2026-03-0398

LIU Mingming, XU Ying. Application of machine learning in the ensemble projection of regional extreme climate indices over China[J]. Arid Meteorology, 2026, 44(3): 398-411 DOI:10.11755/j.issn.1006-7639-2026-03-0398

0 引言

全球气候变暖背景下,极端高温、暴雨和干旱等极端气候事件发生频率和强度均呈显著变化趋势,其引发的风险不容忽视(Li and Fang,2016;Cook et al.,2020)。气候变化通过影响降水、蒸发和径流等气象要素,直接或间接地改变水循环过程,导致区域极端水文事件演变更加复杂,对公共卫生、粮食安全、防洪、供水及生态环境安全构成严重威胁(Meaurio et al.,2017;Zheng et al.,2018;Schilling et al.,2020)。全球变暖已成为不争的事实,并将诱发多种极端气候事件(徐影等,2017),中国区域极端气候事件已造成巨大经济损失和人员伤亡(Zhang et al.,2022),准确预估未来极端降水和极端温度变化,对于减轻高温干旱、低温冷害风险以及防洪减灾具有重要意义(徐影,2015;Li et al.,2020)。

全球气候模式(Global Climate Models,GCMs)是研究气候变化最重要的工具之一。自1995年全球耦合模式比较计划(Coupled Model Intercomparison Project,CMIP)启动以来,GCMs被广泛应用于气候变化研究(Eyring et al.,2016;Kim et al.,2020a)。尽管国际耦合模式比较计划已进入CMIP7阶段,但CMIP6仍是当前气候变化研究和未来预估中被广泛应用的重要国际模式比较计划之一(Dunne et al.,2025)。与早期阶段相比,CMIP6在空间分辨率、模式复杂度及关键物理过程刻画等方面均有明显改进(Eyring et al.,2016)。然而,由于对气候系统认识的局限、模式结构简化、参数化方案不确定性以及自然变率的影响,CMIP6在区域尺度上的模拟仍存在较大偏差(Reichler and Kim,2008;Chen and Frauenfeld,2014;Eyring et al.,2016;Wang et al.,2022;Zhang et al.,2022;Zhao et al.,2022)。例如,在中国区域,CMIP6对极端降水的模拟仍难以准确再现其空间分布特征及年际变化(高学杰,2007;任宏利等,2015;童尧等,2017;Abdelmoaty et al.,2021;Wang et al.,2021;Zhu and Yang,2026)。这些偏差显著影响了气候变化情景下极端事件风险评估的可靠性。

为提高气候模式结果在区域气候变化研究和影响评估中的适用性,以往研究发展了多种偏差订正方法,用于对模式输出结果中的系统性误差进行统计订正,常用方法包括线性尺度变换(Linear Scaling,LS)、分位数映射(Quantile Mapping,QM)、累积分布函数转换(Cumulative Distribution Function Transform,CDF-t)、等距累积分布函数匹配(Equidistant Cumulative Distribution Function Matching,EDCDFm)以及分位数增量映射(Quantile Delta Mapping,QDM)等(Teutschbein and Seibert,2012;韩振宇等,2018;Guo et al.,2018;吉璐莹,2021;Tong et al.,2021;Patel et al.,2022)。这些方法在订正模式均值态和分布特征方面已得到广泛应用,并在区域气候变化及影响评估研究中发挥了重要作用(Themeßl et al.,2012;Maraun,2016)。然而,传统偏差订正方法大多依赖预设的统计变换关系,在处理不同模式之间复杂的非线性误差结构以及极端事件模拟时仍存在一定局限,尤其在未来气候情景下,其趋势保持能力和适用性仍需谨慎评估(Maraun,2016;Ranasinghe et al.,2021)。因此,如何在多模式信息整合的基础上进一步提高极端气候指数的模拟精度,仍是当前研究需要关注的问题。

多模式集合(Multi-Model Ensemble, MME)被认为是降低单一模式不确定性的有效手段,可显著提高未来气候变化预估的可靠性(徐影等,2017;Daron et al.,2018;Jose et al.,2022)。传统集合方法主要包括等权平均和加权集合(Weighted Ensemble Mean,WEM),相关研究表明其整体性能优于单一模式(Oh and Suh,2017;Wang et al.,2021;Yue et al.,2021;Wang et al.,2022)。近年来,随着机器学习技术的发展,基于人工智能的多模式集合方法逐渐受到关注。随机森林(Random Forest,RF)、支持向量机(Support Vector Machine,SVM)、人工神经网络(Artificial Neural Network,ANN)等算法能够有效挖掘非线性特征,在极端气候事件模拟中表现出更强的优势(匡志远等,2020;刘新等,2020;Seo and Ahn,2023;Menapace et al.,2025;Spuler et al.,2026)。

研究表明,基于机器学习方法构建的多模式集合在降水和温度模拟方面整体上显著优于传统多模式集合方法和单一全球气候模式。机器学习算法能够更有效地反映气候变量之间的非线性关系,从而提升区域降水和温度模拟的精度(Acharya et al.,2014;Wang et al.,2018;Crawford et al.,2019;Ahmed et al.,2020;Jose et al.,2022)。其中,不同机器学习算法在多模式集合中的表现存在一定差异,但总体上均优于传统等权集合方法,尤其在降水模拟方面优势更为显著(Ahmed et al.,2020;Jose et al.,2022)。

综上所述,尽管机器学习在气候模拟中的应用已取得一定进展,但将机器学习系统应用于中国区域的GCMs偏差订正及极端气候事件评估仍处于起步阶段,尤其是在极端气候指数模拟方面。因此,本文基于CMIP6的19个全球气候模式,采用随机森林(RF)、极端随机树(Extra Trees,ET)和岭回归(Ridge Regression,Ridge)三种机器学习方法,对中国区域未来极端温度和极端降水的时空变化特征进行预估,并与传统的多模式集合平均方法(MME)进行对比,筛选最优模型,进一步提升极端气候事件预估的可靠性,为气候变化影响评估及适应决策提供科学依据。

1 数据来源

1.1 观测数据

观测数据来源于基于全国2 416个国家基准气候站和基本气象站观测数据,采用“距平逼近法”与“角距权重法”插值构建的逐日温度和降水格点化产品(CN05.1)(Xu et al.,2009;吴佳和高学杰,2013),此数据集的空间分辨率为0.25°×0.25°,时间段为1961—2022年,已广泛用于中国区域多个领域的气候变化研究(吴佳等,2015;王苗等,2023;朱欢欢,2023;鞠琴等,2024;赖雨曈和徐影,2024)。基于此逐日温度和降水数据,根据气候变化检测和指数专家组(Expert Team on Climate Change Detection and Indices,ETCCDI)的标准定义(表1),进一步计算极端温度和降水指数:暖昼指数(TX90p)、冷夜指数(TN10p)、日最大降水量(RX1day)以及5 d最大降水量(RX5day),空间分辨率为0.25°×0.25°,时间段为1961—2014年。

表1   4个极端气候指数的定义(ETCCDI)

Tab.1  Definitions of the four extreme climate indices recommended by the ETCCDI

指数缩写指数名称定义单位
TX90p暖昼指数每月日最高气温大于基准期内90%分位值的天数百分率%
TN10p冷夜指数每月日最低气温小于基准期内10%分位值的天数百分率%
RX5day5 d最大降水量每月最大的连续5 d降水量mm
RX1day日最大降水量每月最大日降水量mm

新窗口打开| 下载CSV


1.2 CMIP6全球气候模式数据

全球气候模式数据来源于CMIP6的19个全球气候模式(表2),包括历史模拟和未来预估。历史阶段采用全强迫试验(Historical,HIST)结果,未来预估采用3种共享社会经济路径情景下的模拟结果,即SSP1-2.6、SSP2-4.5和SSP5-8.5,分别对应不同温室气体排放水平下的未来气候变化。

表2   19个CMIP6模式的基本信息

Tab.2  Basic information of the 19 CMIP6 global climate models (GCMs)

序号模式名称模式单位(国家)分辨率 (经度×纬度)
1ACCESS-CM2澳大利亚气候科学研究中心(澳大利亚)1.25°×1.875°
2ACCESS-ESM1-5澳大利亚气候科学研究中心(澳大利亚)1.25°×1.875°
3BCC-CSM2-MR北京气候中心(中国)1.125°×1.125°
4CanESM5加拿大气候模型与分析中心(加拿大)2.8°×2.8°
5CNRM-CM6-1法国国家气象研究中心(法国)1.4°×1.4°
6CNRM-ESM2-1法国国家气象研究中心(法国)1.4°×1.4°
7EC-Earth3EC-Earth联盟(欧洲)0.7°×0.7°
8FGOALS-g3中国科学院大气物理研究所(中国)2.25°×2.25°
9GFDL-ESM4地球物理流体动力学实验室(美国)1.0°×1.25°
10HadGEM3-GC31-LL英国气象局(英国)1.0°×1.0°
11INM-CM4-8俄罗斯数值数学研究所(俄罗斯)1.5°×1.5°
12KIOST-ESM韩国海洋科学技术院(韩国)2.0°×2.0°
13MIROC6日本海洋地球科学技术研究所(日本)1.0°×1.0°
14MPI-ESM1-2-HR马克斯·普朗克气象研究所(德国)1.5°×1.5°
15MPI-ESM1-2-LR马克斯·普朗克气象研究所(德国)1.5°×1.5°
16MRI-ESM2-0气象研究所(日本)1.125°×1.125°
17NESM3国家气候中心(中国)0.75°×0.75°
18NorESM2-LM挪威气候中心(挪威)1.25°×1.25°
19UKESM1-0-LL英国气象局(英国)0.5°×0.5°

新窗口打开| 下载CSV


使用CMIP6的逐日温度和降水数据,依据ETCCDI的标准定义,计算逐月暖昼指数(TX90p)、冷夜指数(TN10p)、日最大降水量(RX1day)以及5 d最大降水量(RX5day),其中TX90p、TN10p分别代表极端高温、极端低温,RX1day、RX5day则代表极端降水事件。为了便于CMIP6模拟结果与观测进行对比,将CMIP6和CN05.1的温度与降水数据均通过双线性插值方法插值至0.5°×0.5°分辨率。由于CN05.1观测资料仅覆盖1961年及之后的时间段,因此选择CMIP6历史模拟数据与CN05.1观测数据重叠的时期进行多模式集合,同时设置基准期为1961—1990年,根据数据的时间范围,将前80%作为机器学习的训练期(1961—2004年),后20%作为验证期(2005—2014年),其中后20%也用来评估19个全球气候模式的优劣。

2 技术路线与研究方法

2.1 技术路线

采用3种机器学习方法(RF、ET和Ridge)进行多模式集合,由于这3种方法都基于CN05.1观测数据进行校正,为了对比的公平性,也采用基于观测数据进行订正的多模式全局偏差订正法(Arithmetic Mean,AM)进行对比。

首先利用1961—2004年的训练期数据对机器学习模型进行训练。在模型训练过程中,将每个像素点在不同时间上的数据分别作为独立样本,以19个CMIP6模式历史阶段的模拟结果作为输入特征,以同期CN05.1观测数据作为目标变量进行训练。这里采用19个模式结果作为特征输入,目的并非简单增加样本量,而是利用不同模式对同一气候信号响应及误差结构的差异,提取其中相对有效的信息,进而构建更优的多模式集成关系,并且相较于传统简单集合有一定优势(Jose et al.,2022;Bilbao-Barrenetxea et al.,2024)。训练完成后,使用CMIP6验证期的数据作为特征输入到训练好的机器学习方案中进行测试,得到验证期各机器学习方案的输出结果。通过对比机器学习输出结果、传统集合方案和CN05.1验证期数据,评估机器学习方案与传统集合方式的优劣并选出性能最佳的机器学习模型。最后,基于选定的最优模型,将未来3种SSP温室气体排放情景下的CMIP6预估数据输入该模型,得到未来的气候预估结果。具体研究框架如图1所示。

图1

图1   研究框架

Fig.1   Research framework


2.2 多模式集合方法

2.2.1 多模式全局偏差订正法

多模式集合平均(MME)是最直观且应用最广泛的多模式集成方法,由于未引入观测数据,容易保留各模式的共同系统偏差,为消除此偏差并与其他方案进行公平比较,在标准算术平均基础上,按照Lenderink等(2007)和Teutschbein与Seibert(2012)的思路,加入了基于观测的全局线性偏差订正(AM)。

偏差订正是通过观测值与原始集合平均值之差对模拟结果进行线性调整。温度和降水都采用绝对偏差校正,订正公式如下:

${\mathit{T}}_{\mathrm{c}\mathrm{o}\mathrm{r}}={\mathit{T}}_{\mathrm{o}\mathrm{r}\mathrm{i}}+\mathit{\mu }\left({\mathit{T}}_{\mathrm{o}\mathrm{b}\mathrm{s}}\right)-\mathit{\mu }\left({\mathit{T}}_{\mathrm{o}\mathrm{r}\mathrm{i}}\right)$
${\mathit{P}}_{\mathrm{c}\mathrm{o}\mathrm{r}}={\mathit{P}}_{\mathrm{o}\mathrm{r}\mathrm{i}}+\mathit{\mu }\left({\mathit{P}}_{\mathrm{o}\mathrm{b}\mathrm{s}}\right)-\mathit{\mu }\left({\mathit{P}}_{\mathrm{o}\mathrm{r}\mathrm{i}}\right)$

式中:TP分别为温度和降水,单位分别为℃、mm;下标cor、ori、obs分别表示订正后、订正前、观测值;μ为平均算子,表示对括号内变量在基准期内进行中国区域整体气候态平均,如μ(Tobs)表示观测温度的区域气候态平均值。

2.2.2 随机森林模型

随机森林(RF)是一种集成学习方法,通过构建多棵决策树并结合它们的预估结果来进行回归预估。

RF方案使用Scikit-learn库中的Random Forest Regressor,以CN05.1格点观测为预测目标(y),以19个CMIP6模式在训练期(1961—2004年)的逐月特征为自变量(X)对极端气候指数进行回归建模。为了优化模型性能,使用网格搜索(Grid Search)寻找最佳的超参数组合。具体优化的超参数包括:决策树的数量(n_estimators),选择100、200和500;每棵决策树的最大深度(max_depth),用于控制树的复杂度,防止过拟合,选择None(不限制深度)、10、20和50;每次节点分裂时考虑的最大特征数量(max_features),选择sqrt(总特征数的平方根)、log2(总特征数以2为底的对数)、0.3(总特征数的30%)和None(使用所有特征)。

模型性能通过时间序列交叉验证评估,采用均方误差(Mean Squared Error,MSE)作为评分标准。由于Scikit-learn内部评分机制定义了负均方误差以保持“分数越大越优”的一致性,而MSE是误差指标,数值越小表示模型越优,因此实际计算时取负MSE作为比较依据,结果分析中仍以正的MSE或均方根误差(Root Mean Square Error,RMSE)展示。

2.2.3 极端随机树

极端随机树(ET)与随机森林类似,但在训练过程中,极端随机树在每次节点分裂时采用完全随机的特征选择,而不考虑最佳分裂。这种完全随机化的方式使得极端随机树通常训练速度较快,并且能够减少方差,避免过拟合。

使用Scikit-learn库中的Extra Trees Regressor对极端气候指数进行回归建模。采用网格搜索对以下超参数进行优化:决策树的数量(n_estimators)选择100、200和500;每棵决策树的最大深度(max_depth)选择None、10和20;每次分裂时考虑的特征数量(max_features)选择sqrt(总特征数的平方根)、log2(总特征数以2为底的对数)、0.3(总特征数的30%)和None(使用所有特征)。

超参数优化通过时间序列交叉验证(Time Series Split)进行,评分标准与RF相同。

2.2.4 岭回归

岭回归(Ridge)是一种带有L2正则化(L2 regularization)的线性回归方法,L2正则化通过在损失函数中加入参数平方和惩罚项,约束参数大小,以减弱多重共线性影响并防止过拟合,因此岭回归能够通过正则化限制模型复杂度,提升模型稳定性。岭回归的目标是最小化以下目标函数:

$\widehat{\mathit{\beta }}=\mathit{a}\mathit{r}{\mathit{g}}_{\mathit{\beta }}\mathrm{m}\mathrm{i}\mathrm{n}\left(\sum _{\mathit{i}=1}^{\mathit{n}}{\left({\mathit{y}}_{\mathit{i}}-{\mathit{x}}_{\mathit{i}}^{\mathrm{\top }}\mathit{\beta }\right)}^{2}+\mathit{\alpha }\sum _{\mathit{j}=1}^{\mathit{p}}{\mathit{\beta }}_{\mathit{j}}^{2}\right)$

式中:$\widehat{\mathit{\beta }}$表示通过岭回归估计得到的最优回归系数向量;argβmin表示在所有可能的回归系数β中,寻找使目标函数取最小值的参数;n表示样本数量;p表示输入特征数量;i表示样本序号;j表示特征序号;xi是输入特征向量,yi是目标变量,β是回归系数,α是正则化参数。式中第一项$\sum _{\mathit{i}=1}^{\mathit{n}}{\left({\mathit{y}}_{\mathit{i}}-{\mathit{x}}_{\mathit{i}}^{\mathrm{\top }}\mathit{\beta }\right)}^{2}$为残差平方和,用于衡量模型拟合误差;第二项$\mathit{\alpha }\sum _{\mathit{j}=1}^{\mathit{p}}{\mathit{\beta }}_{\mathit{j}}^{2}$为引入的正则化惩罚项,使得回归系数的大小受到约束,从而防止模型过于复杂,避免过拟合。

采用Scikit-learn库中的Ridge模型对极端气候指数进行回归建模。岭回归通过L2正则化约束回归系数大小,以提高模型稳定性和泛化能力,尤其适用于特征间存在较强相关性的情况。为优化模型性能,对正则化参数α进行网格搜索,取值为0.1、1、10、50、100,并结合时间序列交叉验证,以负均方误差为评分标准筛选最优参数。

2.3 评估方法

为了评估各订正方法的优劣,采用泰勒图和泰勒技巧评分作为评价标准。泰勒图被广泛用于评估气候模式模拟性能(Taylor,2001;Gleckler et al.,2008),是针对空间场模拟能力比较的评估方式之一,包含三方面信息:空间相关系数(R)、空间标准差(Standard Deviation,SD)以及中心化均方根误差(Centered Root Mean Square Error,CRMSE),计算公式如下:

$\mathit{R}=\frac{\frac{1}{\mathit{N}}\sum _{\mathit{n}=1}^{\mathit{N}}\left({\mathit{f}}_{\mathit{n}}-\stackrel{-}{\mathit{f}}\right)\left({\mathit{r}}_{\mathit{n}}-\stackrel{-}{\mathit{r}}\right)}{{\mathit{\sigma }}_{\mathit{f}}{\mathit{\sigma }}_{\mathit{r}}}$
$\mathrm{C}\mathrm{R}\mathrm{M}\mathrm{S}\mathrm{E}={\left\{\frac{1}{\mathit{N}}\sum _{\mathit{n}=1}^{\mathit{N}}{\left[\left({\mathit{f}}_{\mathit{n}}-\stackrel{-}{\mathit{f}}\right)-\left({\mathit{r}}_{\mathit{n}}-\stackrel{-}{\mathit{r}}\right)\right]}^{2}\right\}}^{\frac{1}{2}}$

式中:N为研究区域内参与统计的空间格点总数;fnrn分别为n格点上的模拟值和观测值;$\stackrel{-}{\mathit{f}}$$\stackrel{-}{\mathit{r}}$分别为模拟和观测场在研究区域内的空间平均值;σfσr分别表示模拟和观测场的空间标准差。

泰勒技巧评分是对泰勒图的数值量化总结,计算如下:

$\mathrm{T}\mathrm{S}\mathrm{S}=\frac{4{\left(1+{\mathit{R}}_{\mathit{m}}\right)}^{2}}{{\left(\frac{{\mathit{\sigma }}_{\mathit{m}}}{{\mathit{\sigma }}_{0}}+\frac{{\mathit{\sigma }}_{0}}{{\mathit{\sigma }}_{\mathit{m}}}\right)}^{2}{\left(1+{\mathit{R}}_{0}\right)}^{2}}$

式中:TSS是泰勒技巧评分;Rm是订正结果与观测的空间相关系数,R0是最大空间相关系数,本文取值0.999;σmσ0是预估结果和观测的气候态空间标准差。TSS值越接近1,说明预估结果和观测的一致程度越高。该技巧评分在已有文献中被广泛使用(Wang et al.,2018;Ngoma et al.,2021)。

3 结果分析

3.1 验证期不同机器模型性能对比

为了对比几种机器学习与传统AM方案之间的优劣,首先关注验证期的各方法与观测的气候态偏差,图2给出验证期(2005—2014年)不同方法与观测的气候态偏差空间分布。可以看出,对于极端温度指数TX90p和TN10p,表现最好的分别为Ridge模型和RF模型,相比于AM方案,TSS分别从0.35、0.27提升至0.45、0.76。对于TX90p,机器学习和AM方案的结果分布类似,机器学习在东北和西北区域减少了一定的暖偏差;而TN10p,机器学习主要呈现暖偏差,数值比AM方案更高。在极端降水指数模拟中,机器学习均表现出较好的模拟能力,尤其是中国西南部地区,有效减少了AM方案存在的预估过高问题,东北地区也有所改善。所有的机器学习方案的模拟效果都优于AM方案,并且在极端降水指数的模拟中所有机器学习方案TSS评分都达0.99,其中ET为TSS最高的方案。

图2

图2   中国区域验证期(2005—2014年)不同方案在各极端气候指数模拟中的气候态偏差(模拟值减观测值)空间分布及对应的TSS评分

(红色数字表示此方案TSS评分最高)

Fig.2   Spatial distribution of climate state bias (simulation values minus observed values) of different schemes in various extreme climate index simulations and the corresponding TSS scores during the verification period (2005-2014) in China

(The highest TSS score in each index is highlighted in red)


图3为验证期4个极端气候指数的泰勒图比较。对于极端温度指数,机器学习对相关系数(R)的提升不如极端降水指数显著,但整体表现仍优于AM,且标准差(σ)更接近观测。对TX90p而言,AM的R为0.42,σ/σobs为0.53,而Ridge方案的σ最接近观测(|∆σ|=0.23,σ/σobs=0.92),明显优于AM(|∆σ|=1.5);对TN10p而言,机器学习的相关系数均高于AM(0.47),其中RF的σ贴近度最佳(σ/σobs=0.98)。总体来看,机器学习方法在温度类指数上的空间相关性提升相对有限,但其标准差更接近观测值,对空间变率幅度的再现能力普遍优于AM。对于极端降水指数(RX1day、RX5day),机器学习显著优于AM,且明显优于单一模式。以RX1day为例,AM的R为0.91、σ/σobs为0.90,而3种机器学习中ET与观测最接近(R为0.99、σ/σobs为0.99),RF与Ridge的R也分别达0.98和0.99。RX5day的结果与RX1day相似,AM的R为0.88、σ/σobs为1.03,而ET的为0.90;Ridge的|∆σ|最小(0.08),ET次之(0.15)。说明在降水极端指数中,机器学习在相关性和方差再现方面整体优于AM,其中ET在相关性上表现最佳。

图3

图3   验证期4个极端气候指数的泰勒图比较

(a)TX90p,(b)TN10p,(c)RX1day,(d)RX5day

[每个子图展示了单个CMIP6模式(灰色标记)及不同方案(彩色标记)与观测的空间相关系数、标准差及中心化均方根误差(红色虚线)的分布]

Fig.3   Taylor diagrams comparison of four extreme climate indices during the verification period

(a) TX90p, (b) TN10p, (c) RX1day, (d) RX5day

(Each subplot shows the distribution of spatial correlation coefficient, standard deviation, and centered root mean square error (CRMSE) (red dashed lines) for a single CMIP6 model (gray markers) and different scenarios (colored markers) compared with observations)


综合来看,TN10p在机器学习中改进更明显,这可能与不同模式对冷夜指数模拟的离散程度较大有关。已有研究表明,温度极端指标在多模式模拟中普遍存在较明显的模式间差异和不确定性,其中冷极端相关指标对模式偏差和区域背景条件更敏感(Sillmann et al.,2013;Kim et al.,2020b)。对于这类受区域气候背景和局地过程影响较强的指数,AM方案难以充分识别不同模式偏差结构的差异,而机器学习能够在观测约束下从多个模式中提取更接近目标变量的有效信息,更容易体现出优势。已有研究也表明,基于机器学习的多模式集成相较于简单平均方法,通常能够更好地改善温度和降水的模拟效果(Ahmed et al.,2020;Jose et al.,2022)。

基于上述客观对比并结合TSS评分,以下针对TX90p采用Ridge方案,TN10p采用RF方案,RX1day和RX5day采用ET方案进行未来预估。

3.2 未来中国地区极端指数时间变化趋势

根据图4表3,中国地区年极端温度指数相对于基准期(1961—1990年)的未来(2024—2100年)变化与温室气体排放情景密切相关。TX90p在所有情景下均呈持续上升趋势,且排放越高,增幅越大。在SSP1-2.6情景下,Ridge方案的增加趋势为1.1%·(10 a)-1,21世纪末增幅约16%;在SSP2-4.5情景下,Ridge方案增加趋势为2.6%·(10 a)-1,到21世纪末增幅约35%;在SSP5-8.5情景下,上升最显著,Ridge方案增加趋势为5.5%·(10 a)-1,到21世纪末增幅约52%。总体来看,AM与Ridge对TX90p的未来变化趋势判断一致,但Ridge方案的预估结果整体低于AM,表明机器学习对极端高温事件增加的估计相对更为保守。

图4

图4   中国区域未来不同SSP情景下极端气候指数距平(相对于1961—1990年)随时间变化趋势

(a) TX90p,(b) TN10p,(c) RX1day,(d) RX5day

(极端温度指数为绝对变化,极端降水指数为相对变化)

Fig.4   Temporal trends of anomalies of extreme climate indices in China under different SSP scenarios (relative to 1961-1990)

(a) TX90p, (b) TN10p, (c) RX1day, (d) RX5day

(The extreme temperature indices are an absolute change, while the extreme precipitation indices are a relative change)


表3   2024—2100年AM与机器学习预估的中国区域极端气候指数不同排放情景下的未来变化趋势(相对于1961—1990年)

Tab.3  Future trends (2024-2100) of extreme climate indices over China estimated by the AM and machine learning under different emission scenarios (relative to 1961-1990)

指数排放情景AM机器学习
TX90pSSP1-2.6+1.07+0.95(Ridge)
SSP2-4.5+2.86+2.62(Ridge)
SSP5-8.5+6.00+5.49(Ridge)
TN10pSSP1-2.6-0.20-0.19(RF)
SSP2-4.5-0.42-0.41(RF)
SSP5-8.5-0.56-0.54(RF)
RX1daySSP1-2.6+0.92+1.01(ET)
SSP2-4.5+1.75+1.63(ET)
SSP5-8.5+3.90+3.76(ET)
RX5daySSP1-2.6+0.98+1.22(ET)
SSP2-4.5+1.71+1.81(ET)
SSP5-8.5+3.73+3.92(ET)

注:趋势值采用线性回归法计算,表示2024—2100年逐年序列的最小二乘回归斜率;所有趋势均通过0.01显著性检验。

新窗口打开| 下载CSV


TN10p则在未来呈持续下降趋势,且下降幅度同样随排放增加而加大。在SSP1-2.6情景下,AM和RF方案的下降趋势均为0.2%·(10 a)-1,21世纪末累计减幅约9%;在SSP2-4.5情景下,RF方案下降趋势为0.4%·(10 a)-1,末期减幅约10%;在SSP5-8.5情景下,下降最明显,RF方案的下降趋势为0.5%·(10 a)-1,到21世纪末下降约12%。总体上,TN10p的机器学习下降幅度略低于AM,说明机器学习对极端冷事件减少的估计也相对更为保守。

RX1day和RX5day在所有排放情景下均表现为持续上升趋势,且增幅随排放情景升高而增加。对于RX1day,在SSP1-2.6、SSP2-4.5和SSP5-8.5情景下,ET方案的变化趋势分别为1.0、1.6和3.8%·(10 a)-1,21世纪末增幅分别约8.5%、13.8%和31.9%。总体上,AM与ET方案对RX1day的趋势判断基本一致,仅在SSP1-2.6情景下ET方案略高于AM,其余情景下ET预估值普遍低于AM,表明机器学习对单日极端降水增强的估计相对更保守。

RX5day的上升幅度整体高于RX1day。在SSP1-2.6、SSP2-4.5和SSP5-8.5情景下,ET方案的变化趋势分别为1.2、1.8、3.9%·(10 a)-1,21世纪末增幅分别约10.4%、15.4%、33.3%。与AM相比,ET方案在3种情景下对RX5day的预估值普遍更高,说明其对持续性极端降水事件增加的响应更强。总体来看,ET方案虽然在早期阶段年际波动较大,但与AM一致地反映出未来极端降水显著增强,并表现出对RX1day较保守、对RX5day更高估的特征。

总体来看,不同机器学习方案对4个极端气候指数未来变化方向的判断较为一致,与AM方案给出的总体趋势基本相同,差异主要体现在变化幅度和局地响应上。这表明机器学习方法并未改变未来极端气候变化的总体判断,而主要是变化强度不同。

3.3 未来中国地区极端气候指数空间变化特征

图5为中国地区在3种温室气体排放情景(SSP1-2.6、SSP2-4.5、SSP5-8.5)下,21世纪末期(2080—2099年),TX90p(Ridge方案)与TN10p(RF方案)的空间距平分布(相对于1961—1990年)。可以看出,TX90p表现为显著增加,且增幅随情景强度明显增大:在SSP1-2.6情景下,多数地区增幅为10%~30%,高值区位于华北平原—黄淮、长江中下游与华南沿海;在SSP2-4.5情景下,增幅上升至15%~50%,内陆部分地区增幅也有所增强;而SSP5-8.5情景下,华北、长江中下游、四川盆地及华南沿海增幅为30%~70%,青藏高原及西北高海拔地区增幅相对较小。TN10p在全国范围普遍降低,且排放情景越高降低幅度越大:在SSP1-2.6情景下,多数地区减幅为2%~11%,东北以及西北干旱区局地减幅可达12%;SSP2-4.5情景下,减幅扩大至4%~12%,西北及其邻近地区减幅相对较明显;SSP5-8.5情景下,华北、东北与西北部分区域出现13%及以上减幅,华中、华南沿海及黄河上游地区变化幅度相对较小。

图5

图5   21世纪末期(2080—2099年)中国区域SSP1-2.6(左)、SSP2-4.5(中)、SSP5-8.5(右)情景下TX90p(上,基于Ridge模型)与TN10p(下,基于RF模型)的距平空间分布(相对于1961—1990年基准期)(单位:%)

Fig.5   Spatial distribution of anomalies of the TX90p (the top, based on the Ridge model) and the TN10p (the bottom, based on the RF model) over China under SSP1-2.6 (the left), SSP2-4.5(the middle), and SSP5-8.5 (the right) scenarios for the late 21st century (2080-2099), relative to the 1961-1990 baseline (Unit: %)


通过对比发现,AM方案的空间格局(图略)与机器学习基本一致,均呈“随排放与时间推进而加剧”的趋势;但Ridge与RF在高排放情景下的变化幅度更集中,区域梯度更清晰。

图6为不同温室气体排放情景下,中国区域RX1day和RX5day相对于1961—1990年的空间相对变化分布。在不同温室气体排放情景下,RX1day和RX5day整体均呈增强趋势,且随排放强度升高而逐步加剧,但西北干旱区仍以减弱或变化不明显为主。在低排放情景(SSP1-2.6)下,RX1day在全国大部分地区增加5%~15%,东北、华北及青藏高原增幅较明显,局地超过15%,而青藏高原北部、新疆及河西走廊等西北干旱区则出现轻微下降(-10%~0);RX5day空间分布与RX1day基本一致,但整体增幅略高(5%~20%),在华北平原、长江中下游及东北地区形成连续增强带。在中排放情景(SSP2-4.5)下,增强信号进一步扩大,RX1day在全国大部分地区增幅达10%~25%,东北和华北为高值区,渤海湾沿岸、山东及长江中下游局地超过25%,而新疆北部、青海西部及甘肃西部等西北干旱区仍以下降为主(-10%~0%);RX5day的增强范围和幅度更大,在东北、华北及长江流域形成连续高值带,普遍增加15%~30%,局地可达35%。在高排放情景(SSP5-8.5)下,极端降水增强最明显,RX1day全国多数地区增幅超过20%,高值区进一步扩展,局地可达35%~40%,呈现明显的“东强西弱”空间格局,而青藏高原西部、新疆及河西走廊等西北干旱区仍维持轻微减少或基本不变;RX5day增强趋势更为突出,全国多数地区相对变化超过25%。

图6

图6   21世纪末期(2080—2099 年)中国区域SSP1-2.6(左)、SSP2-4.5(中)、SSP5-8.5(右)情景下,RX1day(上)与RX5day(下)距平空间分布(相对于1961—1990年)(单位:%)

(均基于ET模型)

Fig.6   The spatial distribution of anomalies of the RX1day (the top) and the RX5day (the bottom) over China under SSP1-2.6 (the left), SSP2-4.5 (the middle), and SSP5-8.5 (the right) scenarios for the late 21st century (2080-2099), relative to the 1961-1990 baseline (Unit: %)

(based on the ET model)


总体来看,RX1day与RX5day均随排放情景升高而显著增强,且增强幅度在东部和南部地区最明显,西北干旱区持续表现为偏弱甚至减少。

4 结论和讨论

本文基于19个CMIP6模式输出资料和CN05.1观测数据,围绕中国地区极端气候指数的模拟评估与未来预估问题,构建了基于机器学习的多模式集合方案。通过随机森林、极端随机树和岭回归等方法,对极端温度和极端降水指数进行集合订正,并与传统算术平均集合方案进行对比,分析了不同排放情景下21世纪末期中国地区极端气候指数的空间变化特征,得到如下主要结论。

1)RF、ET、Ridge在订正极端温度和极端降水方面展现出更强的非线性计算能力,集成多个气候模式进行模拟,明显提高了模型的模拟准确性,对于极端降水的TSS评分提高至0.99,显著改善了传统CMIP6模型和AM方案的模拟性能。

2)中国未来极端气候变化对排放情景高度敏感。各排放情景下,TX90p持续上升,TN10p持续降低;RX1day与RX5day均显著增强,且RX5day增幅更明显,ET方案预估普遍高于AM,表明其对持续性强降水响应更强。

3)未来极端气候变化表现出显著的区域差异。TX90p在中国中部地区增幅最大,南方和青藏高原地区变化较小,TN10p在西北和华北地区显著减少;RX1day和RX5day在未来呈递增趋势,东北地区和华北地区的变化尤为显著。

机器学习在极端气候指数偏差订正中具有明显优势,能够更好刻画非线性关系,并提高极端温度和降水模拟精度,尤其在RX5day等指标上表现更优。同时,其可综合多模式信息,增强对气候变化响应的刻画能力,从而提升对未来预估的可靠性。但也需要指出,本文机器学习方案在未来时段的应用,仍建立在历史阶段的统计关系在未来情景下具有一定稳定性的前提之上。换言之,这类方法能够对历史偏差进行订正,但并不能从根本上替代气候模式对未来变化物理过程的演变。已有研究指出,偏差订正通常是在历史参考期建立观测与模式之间的关系,并将其应用于未来气候情景,其合理性取决于原始模式能否较好表征相关气候过程及其变化信号,而统计订正方法本身并不能弥补模式对未来变化机制模拟的根本性缺陷(Maraun,2016;Zhang et al.,2024;Menapace et al.,2025;Spuler et al.,2026)。 IPCC AR6也指出,偏差调整方法虽然已被广泛用于区域气候信息和影响评估研究,但不同方法在趋势保持能力、适用范围及潜在问题方面仍需谨慎评估(Ranasinghe et al.,2021)。

同时,机器学习仍存在一定局限。除对超参数设置和数据质量较为敏感外,其在气候研究中的另一个突出问题是物理机制解释相对不足。虽然机器学习模型能够在统计意义上提高模拟精度,但其结果更多体现为输入与输出之间的拟合关系,对于不同模式权重分配的物理含义、气候系统内部过程之间的相互作用,以及极端事件形成机制的刻画仍较为有限。因此,机器学习结果在应用中仍需结合气候学背景和已有物理认识加以分析,避免仅依据统计性能对结果作出过强解释。此外,在观测稀疏或资料不完整地区,其模拟效果也可能受到限制;部分方法计算成本较高,在局地尺度上还可能出现过拟合问题。未来仍需在优化模型结构、提高数据质量、增强物理约束和提升计算效率等方面进一步改进,以增强其在极端气候研究中的适用性与稳定性。

参考文献

高学杰, 2007.

中国地区极端事件预估研究

[J]. 气候变化研究进展, 3(3): 162-166.

[本文引用: 1]

韩振宇, 童尧, 高学杰, , 2018.

分位数映射法在RegCM4中国气温模拟订正中的应用

[J]. 气候变化研究进展, 14(4): 331-340.

[本文引用: 1]

吉璐莹, 2021. 基于多模式集成和动力降尺度的东亚夏季降水数值预报研究[D]. 南京: 南京信息工程大学.

[本文引用: 1]

鞠琴, 吴金雨, 王兴平, , 2024.

机器学习算法在降水和气温多模式集成中的应用

[J]. 水资源保护, 40(3): 106-115.

[本文引用: 1]

匡志远, 宋振亚, 董昌明, 2020.

基于机器学习订正模型的未来百年全球海表温度预估研究

[J]. 气候变化研究快报, 9(4): 270-284.

[本文引用: 1]

赖雨曈, 徐影, 2024.

CMIP6全球气候模式对中国地区干旱模拟能力评估与预估

[J]. 大气科学, 48(6): 2 157-2 177.

[本文引用: 1]

刘新, 赵宁, 郭金运, , 2020.

基于LSTM神经网络的青藏高原月降水量预测

[J]. 地球信息科学学报, 22(8): 1 617-1 629.

[本文引用: 1]

任宏利, 吴捷, 赵崇博, , 2015.

MJO预报研究进展

[J]. 应用气象学报, 26(6): 658-668.

[本文引用: 1]

童尧, 高学杰, 韩振宇, , 2017.

基于RegCM4模式的中国区域日尺度降水模拟误差订正

[J]. 大气科学, 41(6):1156-1 166.

[本文引用: 1]

王苗, 秦鹏程, 佘晨旭, , 2023.

基于CMIP6多模式集合的华中地区气候变化模拟与预估

[J]. 气象与环境学报, 39(6): 51-60.

[本文引用: 1]

吴佳, 高学杰, 2013.

一套格点化的中国区域逐日观测资料及与其他资料的对比

[J]. 地球物理学报, 56(4): 1 102-1 111.

[本文引用: 1]

吴佳, 周波涛, 徐影, 2015.

中国平均降水和极端降水对气候变暖的响应:CMIP5模式模拟评估和预估

[J]. 地球物理学报, 58(9): 3 048-3 060.

[本文引用: 1]

徐影, 2015. 中国未来极端气候事件变化预估图集[M]. 北京: 气象出版社:59-110.

[本文引用: 1]

徐影, 周波涛, 吴婕, , 2017.

1.5-4 ℃升温阈值下亚洲地区气候变化预估

[J]. 气候变化研究进展, 13(4): 306-315.

[本文引用: 2]

朱欢欢, 2023. 中国气候未来变化的集合预估比较研究[D]. 南京: 南京信息工程大学.

[本文引用: 1]

ABDELMOATY H M, PAPALEXIOU S M, RAJULAPATI C R, et al, 2021.

Biases beyond the mean in CMIP6 extreme precipitation: A global investigation

[J]. Earth’s Future, 9(10): e2021EF002196. DOI: 10.1029/2021EF002196.

[本文引用: 1]

ACHARYA N, SHRIVASTAVA N A, PANIGRAHI B K, et al, 2014.

Development of an artificial neural network based multi-model ensemble to estimate the northeast monsoon rainfall over south peninsular India: An application of extreme learning machine

[J]. Climate Dynamics, 43(5): 1 303-1 310.

[本文引用: 1]

AHMED K, SACHINDRA D A, SHAHID S, et al, 2020.

Multi-model ensemble predictions of precipitation and temperature using machine learning algorithms

[J]. Atmospheric Research, 236: 104806. DOI:10.1016/j.atmosres.2019.104806.

[本文引用: 3]

BILBAO-BARRENETXEA N, MARTÍNEZ-ESPAÑA R, JIMENO-SÁEZ P, et al, 2024.

Multi-model ensemble machine learning approaches to project climatic scenarios in a river basin in the Pyrenees

[J]. Earth Systems and Environment, 8(4): 1 159-1 177.

[本文引用: 1]

This study employs machine learning algorithms to construct Multi Model Ensembles (MMEs) based on Regional Climate Models (RCMs) within the Esca River basin in the Pyrenees. RCMs are ranked comprehensively based on their performance in simulating precipitation (pr), minimum temperature (tmin), and maximum temperature (tmax), revealing variability across seasons and influenced by the General Circulation Model (GCM) driving each RCM. The top-ranked approach is used to determine the optimal number of RCMs for MME construction, resulting in the selection of seven RCMs. Analysis of MME results demonstrates significant improvements in precipitation on both annual and seasonal scales, while temperature-related enhancements are more subtle at the seasonal level. The effectiveness of the ML–MME technique is highlighted by its impact on hydrological representation using a Temez model, yielding outcomes comparable to climate observations and surpassing results from Simple Ensemble Means (SEMs). The methodology is extended to climate projections under the RCP8.5 scenario, generating more realistic information for precipitation, temperature, and streamflow compared to SEM, thus reducing uncertainty and aiding informed decision-making in hydrological modeling at the basin scale. This study underscores the potential of ML–MME techniques in advancing climate projection accuracy and enhancing the reliability of data for basin-scale impact analyses.

CHEN L, FRAUENFELD O W, 2014.

A comprehensive evaluation of precipitation simulations over China based on CMIP5 multimodel ensemble projections

[J]. Journal of Geophysical Research: Atmospheres, 119(10): 5 767-5 786.

[本文引用: 1]

COOK B I, MANKIN J S, MARVEL K, et al, 2020.

Twenty-first century drought projections in the CMIP6 forcing scenarios

[J]. Earth’s Future, 8(6): e2019EF001461. DOI: 10.1029/2019EF001461.

[本文引用: 1]

CRAWFORD J, VENKATARAMAN K, BOOTH J, 2019.

Developing climate model ensembles:A comparative case study

[J]. Journal of Hydrology, 568: 160-173.

[本文引用: 1]

DARON J, MACADAM I, KANAMARU H, et al, 2018.

Providing future climate projections using multiple models and methods:Insights from the Philippines

[J]. Climatic Change, 148(1): 187-203.

[本文引用: 1]

DUNNE J P, HEWITT H T, ARBLASTER J M, 2025.

An evolving Coupled Model Intercomparison Project phase 7 (CMIP7) and Fast Track in support of future climate assessment

[J]. Geoscientific Model Development, 18(19): 6 671-6 700.

[本文引用: 1]

. The Coupled Model Intercomparison Project (CMIP) coordinates community-based efforts to answer key and timely climate science questions, facilitate delivery of relevant multi-model simulations through shared infrastructure, and support national and international climate assessments. Generations of CMIP have evolved through extensive community engagement from punctuated phasing into more continuous support for the design of experimental protocols, infrastructure for data publication and access, and public delivery of climate information. We identify four fundamental research questions motivating a seventh phase of coupled model intercomparison relating to patterns of sea surface temperature change, changing weather, the water–carbon–climate nexus, and tipping points. Key CMIP7 advances include an expansion of baseline experiments, a focus on CO2-emissions-driven experiments, sustained support for community MIPs, periodic updating of historical forcings and diagnostics requests, and a collection of prioritized experiments, or the “Assessment Fast Track”, drawn from community MIPs to support climate research, assessment, and service goals across prediction and projection, characterization, attribution, and process understanding.

EYRING V, BONY S, MEEHL G A, et al, 2016.

Overview of the Coupled Model Intercomparison Project Phase 6 (CMIP6) experimental design and organization

[J]. Geoscientific Model Development, 9(5): 1 937-1 958.

[本文引用: 3]

. By coordinating the design and distribution of global climate model simulations of the past, current, and future climate, the Coupled Model Intercomparison Project (CMIP) has become one of the foundational elements of climate science. However, the need to address an ever-expanding range of scientific questions arising from more and more research communities has made it necessary to revise the organization of CMIP. After a long and wide community consultation, a new and more federated structure has been put in place. It consists of three major elements: (1) a handful of common experiments, the DECK (Diagnostic, Evaluation and Characterization of Klima) and CMIP historical simulations (1850–near present) that will maintain continuity and help document basic characteristics of models across different phases of CMIP; (2) common standards, coordination, infrastructure, and documentation that will facilitate the distribution of model outputs and the characterization of the model ensemble; and (3) an ensemble of CMIP-Endorsed Model Intercomparison Projects (MIPs) that will be specific to a particular phase of CMIP (now CMIP6) and that will build on the DECK and CMIP historical simulations to address a large range of specific questions and fill the scientific gaps of the previous CMIP phases. The DECK and CMIP historical simulations, together with the use of CMIP data standards, will be the entry cards for models participating in CMIP. Participation in CMIP6-Endorsed MIPs by individual modelling groups will be at their own discretion and will depend on their scientific interests and priorities. With the Grand Science Challenges of the World Climate Research Programme (WCRP) as its scientific backdrop, CMIP6 will address three broad questions: – How does the Earth system respond to forcing? – What are the origins and consequences of systematic model biases? – How can we assess future climate changes given internal climate variability, predictability, and uncertainties in scenarios? This CMIP6 overview paper presents the background and rationale for the new structure of CMIP, provides a detailed description of the DECK and CMIP6 historical simulations, and includes a brief introduction to the 21 CMIP6-Endorsed MIPs.

GLECKLER P J, TAYLOR K E, DOUTRIAUX C, 2008.

Performance metrics for climate models

[J]. Journal of Geophysical Research: Atmospheres, 113(D6): 2007JD008972. DOI: 10.1029/2007JD008972.

[本文引用: 1]

GUO L Y, GAO Q, JIANG Z H, et al, 2018.

Bias correction and projection of surface air temperature in LMDZ multiple simulation over central and Eastern China

[J]. Advances in Climate Change Research, 9(1): 81-92.

[本文引用: 1]

JOSE D M, VINCENT A M, DWARAKISH G S, 2022.

Improving multiple model ensemble predictions of daily precipitation and temperature through machine learning techniques

[J]. Scientific Reports, 12(1):4678. DOI:10.1038/s41598-022-08786-w.

[本文引用: 5]

Multi-Model Ensembles (MMEs) are used for improving the performance of GCM simulations. This study evaluates the performance of MMEs of precipitation, maximum temperature and minimum temperature over a tropical river basin in India developed by various techniques like arithmetic mean, Multiple Linear Regression (MLR), Support Vector Machine (SVM), Extra Tree Regressor (ETR), Random Forest (RF) and long short-term memory (LSTM). The 21 General Circulation Models (GCMs) from National Aeronautics Space Administration (NASA) Earth Exchange Global Daily Downscaled Projections (NEX-GDDP) dataset and 13 GCMs of Coupled Model Inter-comparison Project, Phase 6 (CMIP6) are used for this purpose. The results of the study reveal that the application of a LSTM model for ensembling performs significantly better than models in the case of precipitation with a coefficient of determination (R) value of 0.9. In case of temperature, all the machine learning (ML) methods showed equally good performance, with RF and LSTM performing consistently well in all the cases of temperature with R value ranging from 0.82 to 0.93. Hence, based on this study RF and LSTM methods are recommended for creation of MMEs in the basin. In general, all ML approaches performed better than mean ensemble approach.© 2022. The Author(s).

KIM Y H, MIN S K, ZHANG X B, et al, 2020a.

Evaluation of the CMIP6 multi-model ensemble for climate extreme indices

[J]. Weather and Climate Extremes, 29:100269. DOI: 10.1016/j.wace.2020.100269.

[本文引用: 1]

KIM Y, ROCHETA E, EVANS J P, et al, 2020b.

Impact of bias correction of regional climate model boundary conditions on the simulation of precipitation extremes

[J]. Climate Dynamics, 55(11): 3 507-3 526.

[本文引用: 1]

LENDERINK G, BUISHAND A, VAN DEURSEN W, 2007.

Estimates of future discharges of the river Rhine using two scenario methodologies: Direct versus delta approach

[J]. Hydrology and Earth System Sciences, 11(3):1145-1 159.

[本文引用: 1]

LI D W, YUAN J C, KOPP R E, 2020.

Escalating global exposure to compound heat-humidity extremes with warming

[J]. Environmental Research Letters, 15(6):064003. DOI:10.1088/1748-9326/ab7d04.

[本文引用: 1]

Heat stress harms human health, agriculture, the economy, and the environment more broadly. Exposure to heat stress is increasing with rising global temperatures. While most studies assessing future heat stress have focused on surface air temperature, compound extremes of heat and humidity are key drivers of heat stress. Here, we use atmospheric reanalysis data and a large initial-condition ensemble of global climate model simulations to evaluate future changes in daily compound heat-humidity extremes as a function of increasing global-mean surface air temperature (GSAT). The changing frequency of heat-humidity extremes, measured using wet bulb globe temperature (WBGT), is strongly related to GSAT and, conditional upon GSAT, nearly independent of forcing pathway. The historical ∼1°C of GSAT increase above preindustrial levels has already increased the population annually exposed to at least one day with WBGT exceeding 33°C (the reference safety value for humans at rest per the ISO-7243 standard) from 97 million to 275 million. Maintaining the current population distribution, this exposure is projected to increase to 508 million with 1.5°C of warming, 789 million with 2.0°C of warming, and 1.22 billion with 3.0°C of warming (similar to late-century warming projected based on current mitigation policies).

LI Z Y, FANG H Y, 2016.

Impacts of climate change on water erosion: A review

[J]. Earth-Science Reviews, 163: 94-117.

[本文引用: 1]

MARAUN D, 2016.

Bias correcting climate change simulations: A critical review

[J]. Current Climate Change Reports, 2(4): 211-220.

[本文引用: 3]

MEAURIO M, ZABALETA A, BOITHIAS L, et al, 2017.

Assessing the hydrological response from an ensemble of CMIP5 climate projections in the transition zone of the Atlantic region (Bay of Biscay)

[J]. Journal of Hydrology, 548: 46-62.

[本文引用: 1]

MENAPACE A, DHAWAN P, DALLA TORRE D, et al, 2025.

Review of bias correction methods for climate model outputs in hydrology

[J]. Journal of Hydrology,660:133213. DOI: 10.1016/j.jhydrol.2025.133213.

[本文引用: 2]

NGOMA H, WEN W, AYUGI B, et al, 2021.

Evaluation of precipitation simulations in CMIP6 models over Uganda

[J]. International Journal of Climatology, 41(9): 4 743-4 768.

[本文引用: 1]

OH S G, SUH M S, 2017.

Comparison of projection skills of deterministic ensemble methods using pseudo-simulation data generated from multivariate Gaussian distribution

[J]. Theoretical and Applied Climatology, 129(1): 243-262.

[本文引用: 1]

PATEL J, GNANASEELAN C, CHOWDARY J S, et al, 2022.

A quantile mapping approach-based bias correction in Coupled Model Intercomparison Project Phase 5 models for decadal temperature predictions over India

[J]. International Journal of Climatology, 42(4): 2 455-2 469.

[本文引用: 1]

RANASINGHE R, RUANE A C, VAUTARD R, et al, 2021. Climate change information for regional impact and for risk assessment[M]//MASSON-DELMOTTE V,ZHAI P,PIRANI A, et al. Climate Change 2021: The Physical Science Basis. Cambridge: Cambridge University Press: 1 767-1 926.

[本文引用: 2]

REICHLER T, KIM J, 2008.

How well do coupled models simulate today’s climate?

[J]. Bulletin of the American Meteorological Society, 89(3): 303-312.

[本文引用: 1]

SCHILLING J, HERTIG E, TRAMBLAY Y, et al, 2020.

Climate change vulnerability, water resources and social implications in North Africa

[J]. Regional Environmental Change, 20(1):15. DOI: 10.1007/s10113-020-01597-7.

[本文引用: 1]

North Africa is considered a climate change hot spot. Existing studies either focus on the physical aspects of climate change or discuss the social ones. The present article aims to address this divide by assessing and comparing the climate change vulnerability of Algeria, Egypt, Libya, Morocco, and Tunisia and linking it to its social implications. The vulnerability assessment focuses on climate change exposure, water resources, sensitivity, and adaptive capacity. The results suggest that all countries are exposed to strong temperature increases and a high drought risk under climate change. Algeria is most vulnerable to climate change, mainly due to the country’s high sensitivity. Across North Africa, the combination of climate change and strong population growth is very likely to further aggravate the already scarce water situation. The so-called Arab Spring has shown that social unrest is partly caused by unmet basic needs of the population for food and water. Thus, climate change may become an indirect driver of social instability in North Africa. To mitigate the impact of climate change, it is important to reduce economic and livelihood dependence on rain-fed agriculture, strengthen sustainable land use practices, and increase the adaptive capacity. Further, increased regional cooperation and sub-national vulnerability assessments are needed.

SEO G Y, AHN J B, 2023.

Comparison of bias correction methods for summertime daily rainfall in South Korea using quantile mapping and machine learning model

[J]. Atmosphere, 14(7): 1057. DOI: 10.3390/atmos14071057.

[本文引用: 1]

This study compares the bias correction techniques of empirical quantile mapping (QM) and the Long Short-Term Memory (LSTM) machine learning model for summertime daily rainfall simulation focusing on precipitation-dependent bias and temporal variation. Numerical experiments using Weather Research and Forecasting (WRF) were conducted over South Korea with lateral boundary conditions of ERA5 reanalysis data. For the spatial distribution of mean summertime rainfall, the bias-uncorrected WRF simulation (WRF_RAW) showed dry bias for most of the region of South Korea. The WRF results corrected by QM and LSTM (WRF_QM and WRF_LSTM, respectively) were improved for the mean summer rainfall simulation with the root mean square error values of 0.17 and 0.69, respectively, which were smaller than those of the WRF_RAW (1.10). Although the WRF_QM performed better than the WRF_LSTM in terms of the summertime mean and monthly precipitation, the WRF_LSTM presented a closer interannual rainfall variation to the observation than the WRF_QM. The coefficient of determination for calendar-day mean rainfall was the highest in the following order: the WRF_LSTM (0.451), WRF_QM (0.230), and WRF_RAW (0.201). However, the WRF_LSTM had a limitation in reproducing extreme rainfall exceeding 50 mm/day due to the few cases of extreme precipitation in training data. Nevertheless, the WRF_LSTM better simulated the observed light-to-moderate precipitation (10–50 mm/day) than the others.

SILLMANN J, KHARIN V V, ZHANG X, et al, 2013.

Climate extremes indices in the CMIP5 multimodel ensemble: Part 1. Model evaluation in the present climate

[J]. Journal of Geophysical Research: Atmospheres, 118(4): 1 716-1 733.

[本文引用: 1]

SPULER F R, WESSEL J B, JEBEILE J, et al, 2026.

Bias adjustment and the question of usable climate information: Methodological assumptions and value judgements

[J]. Bulletin of the American Meteorological Society, 107(1): E79-E102.

[本文引用: 2]

TAYLOR K E, 2001.

Summarizing multiple aspects of model performance in a single diagram

[J]. Journal of Geophysical Research: Atmospheres, 106(D7): 7 183-7 192.

[本文引用: 1]

A diagram has been devised that can provide a concise statistical summary of how well patterns match each other in terms of their correlation, their root‐mean‐square difference, and the ratio of their variances. Although the form of this diagram is general, it is especially useful in evaluating complex models, such as those used to study geophysical phenomena. Examples are given showing that the diagram can be used to summarize the relative merits of a collection of different models or to track changes in performance of a model as it is modified. Methods are suggested for indicating on these diagrams the statistical significance of apparent differences and the degree to which observational uncertainty and unforced internal variability limit the expected agreement between model‐simulated and observed behaviors. The geometric relationship between the statistics plotted on the diagram also provides some guidance for devising skill scores that appropriately weight among the various measures of pattern correspondence.

THEMEßL M J, GOBIET A, HEINRICH G, 2012.

Empirical-statistical downscaling and error correction of regional climate models and its impact on the climate change signal

[J]. Climatic Change, 112(2): 449-468.

[本文引用: 1]

TEUTSCHBEIN C, SEIBERT J, 2012.

Bias correction of regional climate model simulations for hydrological climate-change impact studies: Review and evaluation of different methods

[J]. Journal of Hydrology, 456: 12-29.

[本文引用: 2]

TONG Y, GAO X J, HAN Z Y, et al, 2021.

Bias correction of temperature and precipitation over China for RCM simulations using the QM and QDM methods

[J]. Climate Dynamics, 57(5): 1 425-1 443.

[本文引用: 1]

Two different bias correction methods, the quantile mapping (QM) and quantile delta mapping (QDM), are applied to simulated daily temperature and precipitation over China from a set of 21st century regional climate model (the ICTP RegCM4) projections. The RegCM4 is driven by five different general circulation models (GCMs) under the representative concentration pathway RCP4.5 at a grid spacing of 25 km using the CORDEX East Asia domain. The focus is on mean temperature and precipitation in December–January–February (DJF) and June–July–August (JJA). The impacts of the two methods on the present day biases and future change signals are investigated. Results show that both the QM and QDM methods are effective in removing the systematic model biases during the validation period. For the future changes, the QDM preserves the temperature change signals well, in both magnitude and spatial distribution, while the QM artificially modifies the change signal by decreasing the warming and modifying the patterns of change. For precipitation, both methods preserve the change signals well but they produce greater magnitude of the projected increase, especially the QDM. We also show that the effects of bias correction are variable- and season-dependent. Our results show that different bias correction methods can affect in different way the simulated change signals, and therefore care has to be taken in carrying out the bias correction process.

WANG B, ZHENG L H, LIU D L, et al, 2018.

Using multi-model ensembles of CMIP5 global climate models to reproduce observed monthly rainfall and temperature with machine learning methods in Australia

[J]. International Journal of Climatology, 38(13): 4 891-4 902.

[本文引用: 2]

WANG D, LIU J H, SHAO W W, et al, 2021.

Comparison of CMIP5 and CMIP6 multi-model ensemble for precipitation downscaling results and observational data: The case of Hanjiang River Basin

[J]. Atmosphere, 12(7): 867. DOI: 10.3390/atmos12070867.

[本文引用: 2]

Evaluating global climate model (GCM) outputs is essential for accurately simulating future hydrological cycles using hydrological models. The GCM multi-model ensemble (MME) precipitation simulations of the Climate Model Intercomparison Project Phases 5 and 6 (CMIP5 and CMIP6, respectively) were spatially and temporally downscaled according to a multi-site statistical downscaling method for the Hanjiang River Basin (HRB), China. Downscaled precipitation accuracy was assessed using data collected from 14 meteorological stations in the HRB. The spatial performances, temporal performances, and seasonal variations of the downscaled CMIP5-MME and CMIP6-MME were evaluated and compared with observed data from 1970–2005. We found that the multi-site downscaling method accurately downscaled the CMIP5-MME and CMIP6-MME precipitation simulations. The downscaled precipitation of CMIP5-MME and CMIP6-MME captured the spatial pattern, temporal pattern, and seasonal variations; however, precipitation was slightly overestimated in the western and central HRB and precipitation was underestimated in the eastern HRB. The precipitation simulation ability of the downscaled CMIP6-MME relative to the downscaled CMIP5-MME improved because of reduced biases. The downscaled CMIP6-MME better simulated precipitation for most stations compared to the downscaled CMIP5-MME in all seasons except for summer. Both the downscaled CMIP5-MME and CMIP6-MME exhibit poor performance in simulating rainy days in the HRB.

WANG D, LIU J H, WANG H, et al, 2022.

Performance evaluations of CMIP6 and CMIP5 models for precipitation simulation over the Hanjiang River Basin, China

[J]. Journal of Water and Climate Change, 13(5): 2 089-2 106.

[本文引用: 2]

Projecting the climate change impacts on hydrology and water resources relies on the climate scenarios simulated by general circulation models (GCMs), which requires a systematic and comprehensive assessment of the GCMs’ simulation performances at a regional scale. This study evaluates the performances of precipitation simulation over the Hanjiang River Basin (HRB) by six climate models from Phase 6 of the Coupled Model Intercomparison Project (CMIP6), the corresponding six previous models from the CMIP5, and their multi-model ensemble (MME) based on the observational data in the CN05.1. To our knowledge, this is the first preliminary study in the HRB. The Taylor diagram (including standard deviation, root-mean-square difference, and correlation coefficient) and Taylor skill score are used for the evaluation of GCMs’ precipitation simulation performances. The spatial pattern and temporal pattern over the HRB simulated by CMIP6 and CMIP5 models are compared by relative biases. The results of the Taylor diagram and skill score show that CMIP6 models don't necessarily perform better than the corresponding previous CMIP5 models in simulating precipitation over the HRB. The MME exhibits superior performance compared to that of any individual model, and the CMIP6-MME is more skillful than the CMIP5-MME. As to the spatial and temporal variation characteristics, the precipitation biases are both present in CMIP6 and CMIP5 models, and the bias of the CMIP6-MME is lower than that of the CMIP5-MME. The CMIP6 and CMIP5 models overestimate the precipitation from January to June, and simulate larger precipitation biases in the areas and seasons with less precipitation, while they are lower with more precipitation over the HRB. The findings obtained in this study could provide a scientific reference for the research of future hydrological cycle predictions over the HRB.

XU Y, GAO X J, SHEN Y, et al, 2009.

A daily temperature dataset over China and its application in validating a RCM simulation

[J]. Advances in Atmospheric Sciences, 26(4): 763-772.

[本文引用: 1]

YUE Y L, YAN D, YUE Q, et al, 2021.

Future changes in precipitation and temperature over the Yangtze River Basin in China based on CMIP6 GCMs

[J]. Atmospheric Research, 264: 105828. DOI: 10.1016/j.atmosres.2021.105828.

[本文引用: 1]

ZHANG S X, HARROP B, LEUNG L R, et al, 2024.

A machine learning bias correction on large-scale environment of high-impact weather systems in E3SM atmosphere model

[J]. Journal of Advances in Modeling Earth Systems, 16(8): e2023MS004138. DOI: 10.1029/2023MS004138.

[本文引用: 1]

ZHANG X D, REN G Y, YANG Y D, et al, 2022.

Extreme historical droughts and floods in the Hanjiang River Basin, China, since 1426

[J]. Climate of the Past, 18(8): 1 775-1 796.

[本文引用: 2]

. The major droughts and floods in the Hanjiang River\nBasin, central China, have a significant impact on flood prevention and\ncontrol in the middle reaches of the Yangtze River and water resources\nmanagement in the areas of the South–North Water Diversion Middle Line\nProject. However, there is a lack of understanding of the multi-decadal to\ncentennial-scale patterns of extreme droughts and floods in the area.\nApplying the yearly drought and flood records from historical documents and\nprecipitation data in the period of instrumental measurements, this study\nconstructs a time series of extreme droughts and floods in the Hanjiang\nRiver Basin from 1426–2017 and analyzes the temporal and spatial\ncharacteristics of the extreme drought and flood event variations. The results\nshow that there were a total of 45 extreme droughts and 52 extreme floods in the\nbasin over the past 592 years. Extreme droughts and floods were highly\nvariable on a multi-decadal to centennial scale, and the frequencies were\nhigher in the first and last 100 years or so of the study period and\nlower in between. Spatially, the frequencies of extreme droughts and floods\nwere generally higher in the middle and lower reaches than in the upper\nreaches. It was also found that there is a good correlation of drought and\nflood frequencies between the upper Hanjiang River Basin and North China.\nThese results are informative for the study of mechanisms and predictability\nof multi-decadal to centennial-scale variability of extreme hydroclimatic\nevents in the river basin.

ZHAO X, MA X W, CHEN B Y, et al, 2022.

Challenges toward carbon neutrality in China: Strategies and countermeasures

[J]. Resources,Conservation and Recycling,176: 105959. DOI: 10.1016/j.resconrec.2021.105959.

[本文引用: 1]

ZHENG H X, CHIEW F H S, CHARLES S, et al, 2018.

Future climate and runoff projections across South Asia from CMIP5 global climate models and hydrological modelling

[J]. Journal of Hydrology: Regional Studies, 18: 92-109.

[本文引用: 1]

ZHU H H, YANG J N, 2026.

Evaluation of extreme precipitation over East China in CMIP6 models

[J]. Atmosphere, 17(2): 136. DOI: 10.3390/atmos17020136.

[本文引用: 1]

Based on precipitation extremes calculated from high-resolution daily observational data in East China during 1961–2014, the performance of 34 climate models from phase 6 of the Coupled Model Intercomparison Project (CMIP6) are assessed in terms of climatology and interannual variability. Four extreme precipitation indices, including the total precipitation (Prcptot), the total precipitation for events exceeding the 95th percentile (R95p), and the maximum of 1-day (Rx1day) and 5-day (Rx5day) precipitation, are analyzed. Results show that the CMIP6 models demonstrate good performances in reproducing the climatological spatial distribution and interannual variability of precipitation extremes, with the best from Prcptot. Based on an integrated assessment of the above two factors, the models that perform relatively well for all four extreme precipitation indices are GFDL-CM4, MIROC6, EC-Earth3-Veg, EC-Earth3, and EC-Earth3-CC. Furthermore, the optimal multi-model ensemble (A-MME) constructed from a selection of the most skillful models shows improved behavior compared to the all-model ensemble. The wet (dry) biases over the northern (southern) region of East China are all decreased. This may benefit from the improvement that A-MME can reproduce well the characteristics of moisture and vertical velocity.

/