利用BMA加权集合方法预测年代际极端气温.docVIP

利用BMA加权集合方法预测年代际极端气温 姚遥，罗勇，黄建斌，赵宗慈 清华大学地球系统科学研究中心，100084 摘要 基于12个参加了耦合模式比较计划第五阶段（CMIP5）模式的年代际试验结果，评估了模式对全球和区域年代际尺度极端气温的模拟，探索了贝叶斯模型平均（BMA）加权集合方法在年代际极端气温预测中的应用。用广义帕累托分布拟合气温极值，与模式结果推算出的气温极值相比较，模式趋于高估20年一遇最高（低）气温的趋势，模式对20年一遇最高温的模拟优于最低温。比较模式集合平均（MME）和BMA加权集合结果后发现，BMA加权集合在点预报上并无优势，但由于模式间的不确定性较大，BMA加权集合的概率预报效果更好。 关键词：极端气温，年代际预测，CMIP5，BMA 引言 在利用多模式结果预估未来气候变化时，一个值得考虑的问题是给予各模式相应的权重。其中较为常见的一种选择是等权重模式集合，研究也表明多模式集合平均对现代气候的模拟能力优于多数单个模式的结果[1]。另一方面，可以给予表现较好的模式较高的权重，即应用加权集合方法。 在过去的加权集合方法研究中，有多种定权的原则，如依据模式对当前气候平均态、年际变率或者趋势的模拟表现来直观地确定[2-5]，也有通过贝叶斯方法来进行多模式集合[6]。多模式加权集合方法的假设前提是，模式对当前气候的模拟能力与其对未来的预估能力相关。可是这种关联性现在还无法验证，或只能间接地比较模式间的模拟来验证。 上述研究主要是针对百年尺度的气候预估，所以集合方法与天气预报或者季节预测中常用的集合预报方法有所不同。另一方面，年代际气候预测在时间尺度上介于此二者之间，作为一个新的研究领域，针对年代际试验的加权集合方法研究较少，本文旨在尝试和探索其可行的集合预测方法。 年代际尺度的气候预测既需要考虑外强迫的作用，又需要考虑到气候系统内部的自然变率的影响。年代际的可预测性研究表明，中纬度海洋及大陆尺度的陆地气温都有很高的潜在可预报性[7-8]。已有的年代际试验结果也表明，在初始场同化了海洋信息以后，模式对北大西洋的年代际预测有较大的提升[9]，对全球气温的模拟也有一定的技巧，但同未经初始化的CMIP3模式20世纪模拟试验相比，差别并不显著[10]。 参加CMIP5的模式很多都做了年代际试验，已有的评估结果显示，模式对印度洋、北大西洋和西太平洋温度模拟较好，多数模式对大西洋多年代际涛动（AMO）模拟较好[11]，去除掉ENSO的影响后，模式对全球平均温度的回报有一定的技巧[12]。 已有的研究多数只针对全球平均温度[13]，并未对区域尺度的温度极值做深入的评估，本研究结合多模式集合方法，希望做进一步探讨。 1 CMIP5年代际试验与误差订正 本研究利用了12个CMIP5模式的年代际试验结果，每个模式至少有3组试验结果，每组试验包含10次模拟，分别为以1960年，每隔5年至2005年为起报时刻，向后积分10年。这12个模式中，除了MOHC的HadCM3是以观测海温的距平为初始场外，其余模式均为以观测海温的绝对值为初始场。模式的分辨率等信息参见表1。 表1 12个CMIP5模式基本信息 模式名称 机构 集合数 大气分辨率 bcc-csm1-1 BCC 4 64×128 T42L26 CanCM4 CCCma 10 64×128 T63L35 CNRM-CM5 CNRM-CERFACS 10 128×256 T127L31 FGOALS-s2 LASG-IAP 3 108×128 GEOS-5 NASA-GMAO 3 91×144 HadCM3 MOHC 10 73×96 IPSL-CM5A-LR IPSL 6 96×96 MIROC4h MIROC 3 320×640 T213L56 MIROC5 MIROC 6 128×256 T85L40 MPI-ESM-LR MPI-M 10 96×192 T63L47 MPI-ESM-MR MPI-M 3 96×192 T63L47 MRI-CGCM3 MRI 3 160×320 T159L48 本研究所用变量为逐月2m气温，与之对比的观测为CPC的GHCN_CAMS全球陆地格点逐月观测气温以及ERA再分析资料的2m气温。在所有分析前，先将模式结果、观测和再分析资料双线性插值到128×256高斯格点 对于初始化的数值模拟，模式不可避免地会向其自身的气候态漂移，本研究按照DCPP推荐的方法对模式结果进行误差订正，以消除模式气候漂移的影响。在不同起报时间，各自单独进行误差订正，订正项为，其中为起报时间为t时的模式模拟气候态，为相应时段的观测气候态。 图1 模式简单平均集合（MME）和加权集合（BMA）模拟的极端（上）低温和（下）高温，右图为对应散点图，实线为模式