基于LMR的中国石笋氧同位...去千年气候同化中的应用初探_彭友兵.pdf

第 卷 第 期 年 月第 四 纪 研 究，彭友兵，吴怡，程海，等 基于 的中国石笋氧同位素记录在过去千年气候同化中的应用初探 第四纪研究，（）：.，（）：.：文章编号：（）基于 的中国石笋氧同位素记录在过去千年气候同化中的应用初探彭友兵，吴怡，程海，郝赢玺，崔瑜（西安交通大学人居环境与建筑工程学院地球环境科学系，陕西 西安；西安交通大学全球环境变化研究院，陕西 西安）摘要：本研究利用中国石笋氧同位素记录填补了（）的代用资料数据库在中国东部的空白，并分析了代用资料库优化后 对中国降水的同化能力。结果表明，中国石笋氧同位素记录的添加提升了 对中国东部 以南和西北地区降水的同化能力，能较好再现过去 年两个地区发生的旱涝事件。同化数据表明过去 年中国东部发生了干旱和雨涝事件各 次，西北地区发生了 次干旱事件和 次雨涝事件。其中中国东部 年、年干旱事件和 年、年雨涝事件可能主要受北极涛动（）、北大西洋涛动（）和大西洋多年代际振荡（）共同作用，当、和 处于显著正位相则中国东部易发生干旱事件，负位相则易发生雨涝事件；年干旱事件和 年、年和 年雨涝事件则主要受太平洋年代际振荡（）影响，当 处于显著负位相则中国东部易发生干旱事件，正位相则易发生雨涝事件；年干旱事件主要受 负位相、和 正位相影响；而 年干旱事件与大气涛动间的联系并不密切。在中国西北地区，世纪末期之后接连发生的降水异常事件可能与、和 的位相转变存在联系。关键词：；同化；旱涝事件；大气涛动中图分类号：文献标识码：开放科学标识码（）收稿，收修改稿 国家自然科学基金项目（批准号：）资助 第一作者（通讯作者）简介：彭友兵，女，岁，副教授，气候变化研究，：引言过去千年气候变化研究对于我们辨识现代气候变化的自然背景，评估全球变暖自然和人为因素等问题具有重要意义。利用代用资料重建和数值模拟两种研究方法，国内外对过去千年气候变化的研究取得了显著的进展，一方面通过历史文献资料和代用资料获得了大量的反映过去千年气候变化的高分辨率重建数据，另一方面利用数值模式进行模拟，对过去千年气候变化形成机制进行了分析。但上述两种研究过去千年气候变化的方法都有自己的局限性。重建数据的局限性在于：）重建资料数据点稀少及空间分布的不均匀，使得重建资料受地域性限制，故而存在代表性不足的问题；）不同的代用指标对气候变化敏感性不同，导致对过去千年气候变化的认识存在不确定性；）代用资料作为观测数据，其不可避免的存在代表性误差和器测误差；）重建数据缺乏大气环流场的相关数据资料，从而无法对过去气候变化的形成机制进行探讨，记录与记录的简单对比远不足以揭示过去千年气候演变的驱动机制。而数值模拟的局限性在于：模式模拟结果存在高的不确定性，主要是由驱动模式的重建强迫因子、模式巨大的内部变率噪音和各模式间采用不同的物理参数等造成，因此模拟结果不一定能完全再现重建数据所揭示的气候变化事实。重建和数值模拟方法各自的缺陷，使得过去千年气候变化的研究面临了诸多挑战和争议，阻碍了对过去千年气候变化发生机理的理解。因此，随着过去千年气候变化研究工作的快速发展，迫切 期彭友兵，等：基于 的中国石笋氧同位素记录在过去千年气候同化中的应用初探需要既能符合气候变化事实又能符合动力学机制的方法，以便更好地揭示过去千年气候变化的演变规律和机制，提高我们对年际年代际百年尺度气候变化动力学的认识。古气候同化方法正是对国际和国内过去千年气候变化研究在这方面需求的适应，是加强对区域气候动力学研究，理解更加完整和真实物理机制的重要手段。数据同化在 年被引入古气候研究中，它可以将代用资料和气候模式的优势协调统一，通过重建数据约束模型模拟的轨迹，从而产生既符合气候模式物理机制，又与重建数据所反映的过去气候事实一致的结果，还能对重建数据和模拟数据误差进行定量估计，成为目前古气候研究领域的重要研究方向之一。但目前将古气候数据同化技术引入中国气候变化研究尚存不足。如目前应用广泛的 种全球古气候同化数据（）、（）和（）在其采用的代用资料数据库中对中国尤其是中国东部考虑不足，代用资料覆盖不充分，早期资料覆盖严重匮乏。在整个中国只考虑了 个树轮代用记录，和 增加了在中国西部的树轮代用记录，但对中国东部 仅仅只考虑了 个树轮代用记录，而 也只是另外再增加了 个石笋氧同位素 记录。代用资料覆盖不足严重影响了同化数据对中国东部过去千年气候的再现能力。基于此，本研究以 为主要研究工具，通过添加中国石笋 数据，丰富古气候同化代用资料库中中国东部代用记录，来优化 代用资料数据集。从而帮助我们提升同化数据对中国过去千年气候变化的再现能力，进而更深入地理解过去千年中国年际年代际百年尺度气候变化的演变规律和物理机制。同化模型与数据.同化模型与过去千年同化试验 由美国华盛顿大学 等搭建，等通过更新和扩充其代用资料数据库以及采用季节回归的正向模型，将 发展至第二个版本；其采用的是集合卡尔曼滤波同化方法。在传统的集合卡尔曼滤波同化过程中，其对状态变量进行更新的核心公式为：（）（）首先由气候模式提供先验值（）或者叫做背景场，然后根据观测数据（或代用资料如树轮宽度、冰芯 比值、石笋氧同位素等）与先验值 的误差对状态变量进行更新，最终生成一系列最优状态估计同化结果（后验值或者叫做分析值）。（）是一个非线性的向量值观测算子，通常是正向模型，用于将模拟变量 从模拟状态空间映射到观测数据（或代用资料）空间；在古气候领域，它将气候模式模拟的气候状态变量（如温度、降水）转换为代用资料（如树轮宽度、冰芯 等）。以古气候数据同化研究中同化石笋氧同位素资料为例，（）既可以是基于线性拟合的回归统计正向模型，也可以是考虑了物理过程的正向模型。因此（）代表气候模拟数据中没有涵括的来自代用资料的新息向量。在公式（）中，新息向量通过 即卡尔曼增益矩阵来与先验值 进行加权更新状态变量。卡尔曼增益矩阵计算公式为：（）其中，公式（）中，指背景场协方差矩阵，指代用资料误差协方差矩阵，表示公式（）的线性观测算子。在 中，采用的集合平方根滤波同化方法对传统的集合卡尔曼滤波更新方程即公式（）进行了调整，其被分为集合平均和相对集合平均的扰动（相对集合平均的距平）两部分来分别计算。可以把气候模拟数据中没有涵括的来自代用资料的新息向量传递到多个气候变量如 位势高度场和海平面气压场，从而实现多变量同化。新信息向量主要通过 赋予权重来向位势高度场和海平面气压场传递。详细介绍请参考 等。由公式（）可知 同化理论框架由代用资料数据库（即公式（）中的）、气候模式模拟资料数据库（即公式（）中的）和正向模型（即公式（）中）三部分组成。的代用资料数据包括两个数据库，其中之一为最新的国际过去全球变化研究计划（，简称）过去 年气候变化集成代用资料数据库（），另一个数据库来自于 等整 理 的 大 量 树 轮 宽 度 记 录。但 在 版 本.中仅考虑了 数据库。本文在此基础上利用中国石笋氧同位素 记录来扩展 的代用资料数据库。中只考虑具有年和半年精度的代用资料数据，其中半第 四 纪 研 究 年年精度的数据均通过平均处理为年精度数据，其涵括的代用数据类型包括树轮宽度、树轮密度、树轮同位素、冰芯 等。的气候模式模拟数据来自于国际耦合模式比较计划第五阶段（）通用气候系统模式（）过去千年试验（）的输出结果，在 中其分辨率由原来的.降为.，取相对于 年平均的距平值，并采用蒙特卡罗随机扰动法产生模拟气候状态的 个集合成员，从而为进行离线同化提供集合背景场。采用回归统计正向模型将气候模式模拟的气候状态变量（如温度、降水）转换为代用资料形式（如树轮宽度、树轮密度、石笋氧同位素等）。利用中国石笋氧同位素 记录来扩展 的代用资料数据库时，首先对中国石笋氧同位素 记录进行了筛选，包括剔除已经包含在.代用资料库中的数据和选取可以代表降水且具有季节或年精度的数据，最终添加数据如图 和表 所示，共添加 个洞穴的 个石笋 记录。图 添加的中国石笋数据洞穴位置示意图.将 代用资料库扩展优化后，进行同化试验，获取降水量、全球平均气温、空气温度、海表面温度、位势高度场和海平面气压同化数据，得到优化了代用资料库后的新 同化数据集，将其命名为 数据集。.数据及分析方法 为了评估 对中国过去千年气候变化的同化能力，我们使用了已有的 同化数据表 添加的中国石笋数据洞穴位置 序号引用文献位置精度时间跨度 等大鱼洞（，）年 年 等董哥洞（，）年 年 等峨眉洞（，）季节 年 等和尚洞（，）年 年李红春等石花洞（，）年 年 等莲花洞（，）年 年 等万象洞（，）年 年 等香龙洞（，）季节 年 等小白龙洞（，）年 年.（下 载 地 址：）和过去千年中国 重 建 降 水 数 据（）来对比验证，并通过计算皮尔逊相关系数（，用 表示）来评估 的同化能力。由于 数据只能反应过去 年降水变化，因此只对过去 年中国降水变化进行研究。其中 计算方法与 等相同，自由度计算方法与 等相同。为了分析 对过去千年中国旱涝事件的再现能力，我们采用 和 对持续性旱涝事件的定义，认为区域平均降水经过 年滑动平均并标准化后，出现持续 年或以上为负值（正值）且其中至少有连续 年的负值（正值）超过一个标准差则为一次持续性干旱（雨涝）事件。为了分析旱涝事件发生机制，分别计算了太平洋年代际振荡（，简称）指数、.指数、北极涛动（，简称）指数、北大西洋涛动（，简称）指数和大西洋多年代际振荡（，简称）指数。指数被定义为北太平洋从 至北极区域（）异常的第一主成分；.指数定义为 、区域的平均 异常；和 指数计算分别为北大 西 洋 在 至 北 极 区 域 和 、区 域 上 空 海 平 面 气 压（，简称）异常的第一主成分（）；期彭友兵，等：基于 的中国石笋氧同位素记录在过去千年气候同化中的应用初探指数基于北大西洋 的平均异常（，）减去全球（）平均 异常。结果.同化能力评估 将 和 .以 为参照来分别图 （）.与 反映的过去 年中国降水皮尔逊相关系数分布图和（）与 反映的过去 年中国降水皮尔逊相关系数分布图所有数据均进行了 年的滑动平均和标准化；图中黑点表示该区域 通过 置信区间检验.（）.，（）分析其对中国降水的同化能力，结果如图 所示。可以看出，在中国东部 以南和中国西北这两个地区，与 的皮尔逊相关系数 数值较 .与 的更高，说明 代用资料库经优化后的 对这两个地区的同化能力较之 .要强。但在中国东部 以北则优化后的 值比优化前变小，甚至为负数，说明数据集相比 .对这个地区的同化能力反而有所下降。这或许说明该区所选取的石笋氧同位素记录不利于降水同化。.同化数据对过去千年中国年代际旱涝事件的再现能力 因为 同化数据提高了在中国东部和西北地区的同化能力，于是重点分析 对中国东部（，）和西北地区（，）年代际旱涝事件的再现能力。分别将这两个地区的降水进行区域平均，再进行 年滑动平均并标准化处理，然后再与 进行比对。数据揭示过去 年中国东部发生了干旱事件 次，分别发生在 年、年、年、年和 年，其 中 次 干 旱 事 件 与 一 致，即 年、年、年和 年干旱事件（图）。与 对比发现，同化数据严重高估了 年、年干旱事件和严重低估了 反映的 年的干旱事件（图）。过去 年中国东部发生雨涝事件 次，分别为 年、年、年、年和 年，其中 次雨涝事件与 反映一致，分别为 年、年、年和 年雨涝事件（图）。总体而言，相 较 于 模 拟 结 果和 其 它 同 化 数 据，在反映中国东部干旱、雨涝事件的能力有所提升。与 对比发现，同化数据高估了 年雨涝事件的严重程度；另外，反映的 年雨涝事件在 中呈现出完全相反的干旱气候。对比中国东部干旱和雨涝事件可以发现在 世纪和 世纪中国东部降水异常事件发生频繁，如 年干旱事件后紧接着就发生了 年雨涝事件，年雨涝事件后仅 年又发生 年雨涝事件，经过 年又发生 年干旱事件，紧接着是 年雨涝事件。数据揭示过去 年中国西北干旱区发生了干旱事件 次，分别发生在 年、年、年、年和 年，且均能与 反应的干旱事件对第 四 纪 研 究 年图 （）和 反映的过去 年中国东部干旱事件发生图（黄色柱状），以及（）和 反映的 过去 年中国东部雨涝事件发生图（蓝色柱状）所有数据均进行了 年的滑动平均和标准化.（）（），（）（）应，但 揭示的 年这一次的干旱事件在 中 呈 现 出 完 全 相 反 的 雨 涝 气 候（图）。显示的 年的干旱事件对应 的 年，图 中可以看出，高估了此次干旱事件的持续时间。反映过去 年中国西北地区发生了雨涝事件 次，分别发生在 年、年、年和 年（图），且均与 反应的雨涝事件对应，但在 中反映的 年和 年雨涝事件的严重程度在 中被低估了。相较于中国东部，西北地区旱涝事件更易接连发生。年干旱事件后仅 年发生 年雨涝事件。年干旱事件后仅 年又发生 年干图 （）和 反映的过去 年中国西北干旱事件发生图（黄色柱状），以及（）和 反映的 过去 年中国西北雨涝事件发生图（蓝色柱状）所有数据均进行了 年的滑动平均和标准化.（）（），（）（）旱事件。年雨涝事件后马上发生 年干旱事件。年干旱事件后紧接着 年雨涝事件发生，经 年后再次发生 年雨涝事件。表 表明，在 与 对应较好的中国东部 次干旱事件中，除 年外其它 个时期中国西北地区也在同期或近期发生了严重的干旱事件，中国东部的 年、年和 年干旱事件分别与中国西北地区 年、年和 年干旱事件对应。揭示中国东部的另外两次干旱事件时期中国西北地区发生了雨涝事件，即中国东部 年和 年干旱事件与中国西北地区 年和 年雨涝事件 期彭友兵，等：基于 的中国石笋氧同位素记录在过去千年气候同化中的应用初探表 揭示的过去 年中国东部和西北地区旱涝事件（单位：年）（：）世纪东部干旱东部雨涝西北干旱西北雨涝 粗体表示 与 对应一致的事件对应。结果基本与 保持一致。但由图 可知 显示中国西北地区在 年发生了干旱事件，因此 表明中国东部 年干旱时期，中国西北地区也同样发生了严重的干旱事件 这与 结果不同；尽管 显示中国东部 年气候偏干（图），却未达到干旱事件程度，这与 一致，都揭示出该时期中国东部偏干西北地区雨涝的降水空间分布特征。然而与干旱事件不同的是，一个地区的雨涝事件却不一定对应着另一个地区的降水异常事件。显示中国东部的雨涝事件中 年雨涝事件与中国西北地区 年干旱事件同期，年和 年雨涝事件与中国西北地区 年雨涝事件同期，但发生在 世纪的其它两次雨涝事件时期（年和 年）中国西北地区未有降水异常事件发生。结果与 基本保持一致。中国西北地区雨涝事件除前面提到的 年、年和 年雨涝事件分别与中国东部的 年干旱事件、年干旱事件和 年及 年雨涝事件对应外，中国西北地区 年雨涝事件时期中国东部未有降水异常事件发生。以上分析结果与各个降水异常事件降水空间分布结果一致（图略）。.揭示的内部变率与过去 年中国年代际旱涝事件间的联系 通过.、和 对降水的回归分析（图），可以发现这些大气涛动均可以图 （）降水与.指数回归分析、（）降水与 指数回归分析、（）降水与 指数回归分析、（）降水 与 指数回归分析以及（）降水与 指数回归分析所有数据均进行了 年的滑动平均和标准化；图中黑点表示该区域通过 置信区间检验.（）.，（），（），（），（）影响中国东部和西北地区的降水，其中 和 影响最为显著，其次为 和，最弱为.，.仅对中国东部 以南地区降水有显著影响（图）。对于中国东部，和 导致其降水南北一致变化，.、和 则导致其南北降水出现反相变化，当 负位相或 正位相时，中国东部降水减少，当.、和 正位相时，华北降水减少而长江中下游地区降水增加，反之亦然。对于中国西北地区，.、和 导致该地区降水南北一致变化，和 则导致其南北降水出现反相变化，当.、和 正位相时，中国西北地区降水增加，当 和 正位相时，西北北部降水减少而南部第 四 纪 研 究 年降水增加，反之亦然。图 为 反映的过去 年、和 指数变化，由于.指数对两个地区的影响不显著，因此没有分析。图 显示中国东部 年干旱事件对应、和 正位相。可能是、和 正位相的共同作用导致了中国东部此次干旱事件的发生。中国东部 年干旱时期各大气涛动处于正常状态，此次干旱事件有可能受其它因素影响。中国东部 年干旱事件与 负位相、和 正位相对应。高估了这次干旱事件的持续时间（图）或许正缘于 高估了该时期 负位相、和 正位相的持续时间。中国东部在 年干旱事件后紧接着就发生了 年雨涝事件，图 显示两次事件的发生正对应 的显著转变，从 年干旱事件时期的负位相转变成了 年雨涝事件时期的正位相，说明 或许是这两次降水异常事件发生的主导因素。世纪前后是中国东部旱涝事件频繁发生的时期，其中 年、年雨涝事件和 年干旱事件和、和 的位相转变相对应。年、年雨涝事件均对应、和 负位相，之后伴随着、和 显著转变为正位相中国东部发生了 年干旱事件，说明、和 的位相转变在这 次降水异常事件中起了主导作用。而紧接着的 年雨涝事件时期、和 却并未发生显著转变，甚至 和 依然与 年干旱事件时期一样呈显著正位相，但 发生了显著位相变化，在该时期呈正位相。年雨涝事件时期，各大气涛动状态也与 年雨涝事件时期相似，表明 或许是这两次雨涝事件发生的主导因素。综上所述，除中国东部 年干旱事件外，其它中国东部降水异常事件可能都受到大气涛动的影响，中国东部 年、年干旱事件和 年、年雨涝事件可能主要受、和 共同作用，另外 和 对 年干旱事件中也起到了重要作用，当、和 处于显著正位相则中国东部易发生干旱事件，负位相则易发生雨涝事件，或许主导了 年干旱事件和 年、年和 年雨涝事件，其在图 （）反映的过去 年中国东部干旱事件（黄色柱状）和雨涝事件（蓝色柱状）时期、和 指数变化，（）反映的过去 年中国西北地区干旱事件（黄色柱状）和雨涝事件（蓝色柱状）时期、和 指数变化所有数据均进行了 年的滑动平均和标准化.（），（）（），（），（）（）年干旱事件中也起到了重要作用，当 处于显著负位相则中国东部易发生干旱事件，正位相则易发生雨涝事件。这与图 回归分析结果相一致。图 显示前面 次西北降水异常事件与大气涛动间没有明显一致的对应关系，但 年雨涝事件和紧接着发生的 年干旱事件对应着 和 从 年正位相转变为 年负位相，以及 年干旱事件后气候转湿，紧接着发生的 年和 年雨涝事件对应着、和 从 年负位相转变为之后的正位相，说明、和 或许对这些降水异常事件的影响较为重要。结合之前中国东部降水异常事件与大气涛动间的可能联系，发现尽管中国西北 年干旱事件与中国东部 年干旱事件以及 年雨涝事件与中国东部 期彭友兵，等：基于 的中国石笋氧同位素记录在过去千年气候同化中的应用初探年和 年雨涝事件发生在同期，但引起该时期两地区降水异常的主导因子可能并不相同。讨论与总结通过在中国增加石笋氧同位素记录来优化 的代用资料库可以提高对中国东部 以南地区和西北地区的同化能力，但对中国东部 以北地区的同化能力反而有所下降，可能是因为该区所选取的石笋氧同位素记录不利于降水同化，以后还需要通过试验来验证这一可能性。代用资料库优化后的 数据可以较好地反映过去 年中国东部和西北地区的旱涝事件。过去 年中国东部发生 年、年、年、年、年共 次干旱事件和 年、年、年、年和 年共 次雨涝事件，除 年干旱事件和 年雨涝事件外均与 结果对应一致；过去 年中国西北干旱区发生了 年、年、年、年和 年共 次干旱事件和 年、年、年和 年共 次雨涝事件，也都能与 结果对应一致。但 对于过去 年最为严重的崇祯大旱没有反映出来，需要未来通过开展更多的试验来进行提升。过去 年中国东部和西北降水异常事件均存在接连发生的现象。在 世纪早期中国东部接连发生了 年干旱事件和 年雨涝事件，在 世纪末到 世界接连发生了 年雨涝事件、年雨涝事件、年干旱事件和 年雨涝事件。西北地区旱涝事件更易接连发生，几乎没有一次降水异常事件是单独发生的。另外，也出现在同一时期中国东部和西北均发生降水异常事件的现象。在中国东部的 年、年和 年干旱事件时期中国西北地区分别发生了 年、年和 年干旱事件。在中国东部 年和 年干旱事件时期中国西北地区分别发生了 年和 年雨涝事件。在中国东部 年雨涝事件和 年、年雨涝事件时期，中国西北地区分别发生了 年干旱事件和 年雨涝事件。中国东部旱涝事件与大气涛动联系较西北地区更为密切。中国东部除 年干旱事件外，其它降水异常事件均受大气涛动影响，且接连发生的降水异常事件往往对应着大气涛动的显著转变。中国东部 年、年干旱事件和 年、年雨涝事件可能主要受、和 共同作用，当、和 处于显著正位相则中国东部易发生干旱事件，负位相则易发生雨涝事件。年干旱事件和 年、年和 年雨涝事件则主要受 影响。当 处于显著负位相则中国东部易发生干旱事件，正位相则易发生雨涝事件。年干旱事件主要受 负位相、和 正位相影响。而 年干旱则可能受其它因素影响。在中国西北地区，仅在 世纪末期之后接连发生的降水异常事件可能与、和 的位相转变存在联系。另外，在同一时期两地区可能均发生降水异常事件，但各自的主导因子可能并不相同。由于大气涛动无论是对中国东部还是对西北地区降水的影响都存在区域不一致性，所以未来还需要从旱涝事件空间分布特征与大气涛动间的联系来进一步探讨其发生机制。降水异常事件可能受多个大气涛动因子的影响，表明未来也需要开展大气涛动对旱涝事件的协同作用研究。致谢：感谢为本文提出了诸多宝贵修改意见的审稿人和杨美芳编辑！参考文献（）：葛全胜，郑景云，郝志新 过去 年中国气候变化研究的新进展 地理学报，（）：.，（）：.程海，张海伟，赵景耀，等 中国石笋古气候研究的回顾与展望 中国科学：地球科学，（）：.，：，（）：.郑景云，刘洋，郝志新，等 过去 年气候变化的全球集成研究进展与展望 第四纪研究，（）：，（）：，：第 四 纪 研 究 年 ，（）：.周天军，邹立维，陈晓龙 第六次国际耦合模式比较计划（）评述 气候变化研究进展，（）：.，（），（）：.王鑫，杨保，王江林 过去千年北半球中高纬夏季温度多年代际变化的归因分析 第四纪研究，（）：，（）：，：，：周雪城，姜大膀 过去千年年代和百年尺度东亚夏季风变化的模拟分析 第四纪研究，（）：.，（）：.朱信国，严蜜，宁亮，等 过去千年 个特征时期东亚冬夏季风关系的模拟研究 第四纪研究，（）：.，（）：.邢辰，孔瑜，刘飞，等 重建与模拟中过去 年厄尔尼诺和东亚夏季风关系的年代际变化 第四纪研究，（）：.，（）：.，（），.易慧郁，刘健，孙炜毅，等 过去 年典型暖期北半球温度变化特征及成因分析 第四纪研究，（）：，（）：.，：，,（）：.，：.，：.：，：.，（）：.：，：.，：，：，,（）：.，,：.，（）：，,（）：.李红春，顾德隆，等 高分辨率洞穴石笋稳定同位素应用之一 京津地区 来的气候变化 记录 中国科学（辑），（）：.，（），（）：.，：，期彭友兵，等：基于 的中国石笋氧同位素记录在过去千年气候同化中的应用初探，,（）：.，,：，,：，：，：.，,（）：.，（）：.，.彭友兵，徐影 过去 年中国东部旱涝事件模拟研究初探 第四纪研究，（）：.，（）：.黎永杰，张仲石，徐天澳 年至 年中国东部的严重干旱 第四纪研究，（）：.，（）：.，（，；，）（），（）（），（）（）,，（）,，：；