降尺度方法及其代码实现

降尺度方法

常用的降尺度方法主要包括三种类别：动力降尺度、统计降尺度以及动力和统计相结合的降尺度方法。

1 动力降尺度

动力降尺度，即区域气候模型(RegionalClimateModels,RCMs)，是建立在系统的大气热力动力学方程基础之上，基于GCM提供的初始和侧边界条件，通过动力学方程之间的关系，推导得到的具有更高分辨率的气候变化情景。
动力降尺度方法具有较强的物理机制，且受观测资料的影响较小，可应用于不同的空间分辨率(一般为15~50km)，从而反映影响区域气候的地面属性和未来气候的波动特征。然而，现阶段对气候系统内在的复杂性及其不确定性仍然难以有效刻画和描述，尤其对气候的动力过程缺乏深入认识，动力降尺度模型的输出变量仍然具有较大偏差，例如通常容易低估降水极值和高估降水的发生频率。尽管RCM的空间分辨率已经有了很大提高，仍然不能满足中小流域的实际应用，尤其是站点尺度气候变化的影响评估。同时，RCM的运行是基于GCM提供的侧边界条件，它的输出变量会在一定程度上遗传GCM的系统偏差。另外，由于RCM具有较大的计算量，只能在部分流域产生有限的气候变化情景，难以在实际应用中探究气候变化水文效应评估中的不确定性。

2 统计降尺度

统计降尺度方法通过建立历史阶段气候模式的网格输出变量与站点或区域观测变量之间的统计模型，并将该统计关系应用到未来的气候模拟输出中，从而得到站点或区域尺度的未来气候要素系列。与动力降尺度模型相比，统计降尺度方法具有更小的计算量，同时其模拟性能可以满足要求，由于其使用较为简便在国内外得到广泛应用。
使用统计降尺度方法开展气候变化影响评估时，首先需要对其气候生成情景和响应评估的表现进行评价，可以采用统计检验、相对误差、泰勒图、均方根误差等标准来衡量降尺度方法的可靠性。

2.1 理想预报法(Perfect Prognosis, PP)

理想预报方法是目前应用较为广泛的一种统计降尺度方法，其原理是建立历史时期多个大气环流因子网格输出(预报因子，如大气湿度、太阳辐射和气压等)和站点或区域气象观测要素(预报变量，如气温和降水)之间的统计模型，然后将该函数关系应用于未来时期的大气环流因子网格输出，从而获得站点或区域的未来气候变化情景。
理想预报方法一般假设全球气候模式能够较好地模拟网格尺度的大气环流因子，同时假定大气环流因子与区域或站点尺度的气象要素之间具有较强的相关性，其最主要的假设在于历史时期预报因子与预报变量的统计关系在未来气候变化情景下仍然适用。优选的预报因子通常对理想预报方法能否取得较好的统计降尺度效果具有重要作用，如何选取成为PP 方法的关键。

2.2 模型输出统计法(Model Output Statistics, MOS)

通过建立历史时期降水网格输出与局地或站点实测降水之间的关系，推演未来时期区域或站点降水序列的方法可以称为模型输出统计法（MOS） 方法。类似地，通过建立GCM的气温网格输出与站点或局地的实测气温之间的统计模型，MOS方法也可以用于气温或其他气象要素的统计降尺度。
与PP方法建立网格尺度实测预报因子与站点或局地尺度气象要素关系不同的是，MOS方法直接使用全球气候模式输出的气候变量或与之相对应的观测气候变量之间的关系，用以校正输出气候变量的系统偏差，即常用的偏差校正方法。该方法将历史期实测数据与气候模式输出在统计参数方面的差异定义为模型偏差，而且需要在模拟未来气候系列时应该去除掉该项系统偏差。

2.3 随机天气发生器（Stochastic Weather Generator, SWG）

随机天气发生器的技术要点是首先估计历史时期气象观测数据的统计参数，然后利用随机模拟模型生成与历史阶段气象要素相同统计参数的未来气候序列。历史期的观测数据可能存在较多的缺失值，但是随机天气发生器可以产生任意时长的完整气候系列。对于资料缺失较多的站点或无资料地区，可以通过空间参数插值的方法获得其历史时期的统计特征值，然后通过随机天气发生器产生缺测站点的气象序列。
随机天气发生器的重要功能是生成站点或区域尺度的未来平均气候变化情景，一般可分为参数与非参数模型两种类别。

3 动力和统计相结合的降尺度方法

动力和统计降尺度相结合的方法，其工作流程为：首先基于全球气候模式生成侧边界条件，通过一套热力动力学方程生成区域尺度的气候变量；再使用统计降尺度技术产生站点或区域的气候变化情景，该方法在一定程度上兼顾了动力降尺度和统计降尺度的优点，但实质上只是两种方法的拼接，在气候变化评估的实际应用中存在局限。

参考

1.博士论文-D2019-全球水文气象极值事件对气候变化的响应研究-尹家波