全新世太阳活动对东亚夏季风百年尺度变化影响的模拟研究_陈萌萌.pdf

第 卷 第 期 年 月第 四 纪 研 究，陈萌萌，孙炜毅，刘健，等 全新世太阳活动对东亚夏季风百年尺度变化影响的模拟研究 第四纪研究，（）：，（）：文章编号：（）全新世太阳活动对东亚夏季风百年尺度变化影响的模拟研究陈萌萌，孙炜毅，刘健，宁亮，严蜜，马雅楠，陈少洋（南京师范大学虚拟地理环境教育部重点实验室 江苏省地理环境演化国家重点实验室培育建设点 江苏省地理信息资源开发与利用协同创新中心，南京师范大学地理科学学院，江苏 南京；可持续发展大数据国际研究中心，北京；江苏省大规模复杂系统数值模拟重点实验室，南京师范大学数学科学学院，江苏 南京；青岛海洋科学与技术试点国家实验室海洋气候同位素模拟开放 工作室，山东 青岛；中国科学院地球环境研究所，黄土与第四纪地质国家重点实验室，陕西 西安）摘要：太阳活动对东亚夏季风（）变化有着重要的影响。本研究利用通用地球系统模式（）开展了全新世气候瞬变模拟试验，探究了全新世 降水百年尺度的变化特征及其对太阳活动外强迫的响应机理。结果表明：在变化趋势上，模拟的全新世以来 降水呈减弱趋势，主要受地球轨道参数外强迫的影响。在百年尺度上，中全新世以来太阳活动外强迫引起 降水出现约 年的周期变化特征，这与高分辨率重建资料记录的周期相似。这一百年尺度变化的空间主模态呈南北偶极型分布，北方降水与太阳活动存在显著的同期正相关，当太阳活动增强时东亚海陆热力差异加强，陆地海平面气压降低，引起异常的气旋性环流，有利于北方降水增多。南方降水与太阳活动呈负相关并滞后于太阳活动，太阳活动通过调制赤道太平洋类拉尼娜态的海温梯度，引起西北太平洋风蒸发海温的正反馈机制，触发异常的西北太平洋气旋性环流，从而导致南方降水减少。关键词：东亚夏季风；全新世；太阳活动；古气候模拟；百年尺度中图分类号：文献标识码：开放科学标识码（）收稿，收修改稿 国家自然科学基金项目（批准号：、和）、国家自然科学基金中瑞合作交流项目（批准号：）、可持续发展大数据国际研究中心（批准号：）和江苏省高校优势学科建设项目（批准号：）共同资助 第一作者简介：陈萌萌，女，岁，硕士研究生，地图学与地理信息系统专业，：通讯作者：孙炜毅，：引言季风以近地面盛行风风向随季节有规律变化为主要特征，是全球大气环流的重要组成部分。东亚季风是全球最为活跃的季风环流系统之一，其变化规律对东亚各国社会经济可持续发展等方面具有重要的影响。深入认识东亚季风的变化特征并揭示其成因机制，有助于提高东亚季风的预报能力。全新世（距今 年）是距今最近的一次间冰期，它对人类社会和文明演化有着潜在的影响，它也是深入了解过去气候和环境变化、更好地理解现在和预估未来气候的关键时期，因此，了解全新世气候变化特征及驱动机制至关重要。全新世东亚季风的演变过程与机制是当前学界研究的重点，气候重建工作极大地促进了对全新世东亚季风演变过程的理解，然而，全新世气候代用指标重建结果在反映东亚夏季风（，简称）的长期演变特征上存在争议。中国南部董哥洞、中部三宝洞和北部九仙洞的石笋 记录表明在早全新世东亚夏季风最为强盛，随后东亚夏季风强度呈逐渐下降趋势；同样地，来自中国南方长江流域的神农架大九湖泥炭记录也显示早全新世 最强，气候最为湿润；而公海孢粉、黄土、湖泊等重建记第 四 纪 研 究 年录则显示 最强盛时期出现在中全新世。在百年尺度上，不少学者基于石笋 代用资料发现，全新世以来东亚季风具有明显的约 年的百年周期信号；此外，等对亚洲夏季风区域内的多条高分辨率重建结果进行功率谱分析，发现在中晚全新世期间重建资料显示出 年的百年周期信号。对于东亚季风可能存在的 年的百年尺度周期，已有学者探讨了其成因，比如，等基于和尚洞石笋中的土壤磁性矿物记录重建东亚季风变化，发现其存在着与太阳活动相关的 年周期，认为太阳活动直接影响东亚季风系统；等认为太阳活动通过影响北大西洋和热带太平洋海温变化，二者共同作用引起海气环流变化进而影响东亚季风；通过将济州岛重建记录与太阳总辐照度和热带太平洋沿岸的代用资料进行对比，发现三者具有显著的相关关系，认为太阳活动主要通过改变热带太平洋海温变化影响东亚季风。然而，代用指标时间分辨率存在差异性、空间不均，指示意义也存在不确定性，难以揭示东亚季风百年尺度变化的机理过程。利用地球系统模式进行全新世气候模拟试验可以为解决上述问题提供参考，地球系统模式可以模拟各气候要素场的时空变化特征，定量区分不同外强迫因子的相对贡献。目前针对全新世开展的气候模拟工作主要基于两种方法，其一是固定某一时期的外强迫和边界条件开展平衡态 切片气候模拟试验，如古气候模拟比较计划第二、三、四阶段（，简称 ）开展的中全新世平衡态气候模拟试验；其二是涵盖完整全新世时段的瞬变积分气候模拟试验，目前已发布了多套全新世瞬变模拟数据集，如过去 年的气候瞬变模拟试验（，简称）。然而，基于上述模拟资料，学者们大多侧重于分析全新世长期气候变化趋势或千年尺度气候突变的问题，比如基于 的模拟结果发现，降水自早全新世以来呈下降趋势。另一方面，以往研究发现了太阳活动与 降水之间的关系，比如，靳春寒基于过去千年集合模拟资料（，简称）发现，太阳活动强 年周期会引起类太平洋年代际振荡（，简称）的海温模态，进而在年代际尺度上影响我国北方降水变化；等基于现代器测资料发现，东亚夏季风降水的年代际变化在一定程度上受太阳黑子变化的影响；周天军等利用 模式的过去千年模拟结果探究了东亚夏季风降水的百年尺度变化，并认为太阳活动的百年尺度变化通过改变海陆热力差异来影响东亚夏季风降水。然而，目前的气候模拟工作很少关注全新世以来东亚季风的百年尺度变化过程，已开展的瞬变模拟试验所考虑的外强迫因子主要是地球轨道参数、温室气体、冰盖和淡水注入等，缺乏考虑太阳活动外强迫的作用。尽管有学者基于重建资料结果推测太阳活动影响了全新世东亚季风百年尺度变化，但仍缺乏模拟资料去验证这一猜想，更不清楚太阳活动会如何影响东亚季风的百年尺度变化。为此，本研究将利用通用地球系统模式开展的全新世气候瞬变模拟试验（，简称）结果，探究太阳活动外强迫对全新世 降水百年尺度变化的影响。资料与方法.气候模拟资料 是由美国国家大气研究中心（）推出的通用地球系统模式，我们采用 的 分辨率版本（水平分辨率约为.，垂直方向分为 层；海洋模块在经向 个格点，纬向 个格点，垂直方向分为 层）开展了多个全新世瞬变气候模拟试验。本研究使用其中的 个试验结果进行分析，分别是地球轨道参数试验（，简称）和 太 阳 活 动 试 验（，简称），其中 试验包括了两个效应，分别是太阳常数的变化和地球轨道参数的变化，二者共同作用改变大气顶端接收到的太阳辐射。具体的试验设计见相关文献。需要说明的是，由于受到计算资源和存储空间的限制，模式运行中采用了 倍加速模拟技术，此方法的模拟性能在前人的工作中已得到验证，；此外，太阳活动外强迫数据（）使用的是 等的重建（图 黑线），其变化幅度大约在.至 之间，全时段平均值约为。为了展示全新世太阳活动变化的百年尺度波动变化，我们对这条外强迫序列进行了 年的带通滤波处理（图 红线）。由于该重建的太阳强 期陈萌萌，等：全新世太阳活动对东亚夏季风百年尺度变化影响的模拟研究图 全新世太阳活动变化（单位：）黑色实线为全新世太阳活动重建序列，蓝色虚线为太阳活动全 时段的平均值，红色实线为 年带通滤波处理结果（无量纲）.（：）（）（），（）迫只 包 含 了.的 资 料，因 此试验的时段为.。在这里我们假设太阳活动的作用是基于线性变化的，将与 试验相减结果看作净太阳活动的效应。此外，本文还使用了 开展的过去 年 的 气 候 瞬 变 模 拟 试 验（），。包含地球轨道参数试验、温室气体试验、冰盖试验和淡水注入试验，以及由这 个外强迫因子共同驱动下的全强迫试验。本研究使用 的全强迫试验（，简称）和地球轨道参数试验（）结果进行分析，选择的试验时段为.。.观测 再分析资料 为了验证模式的模拟性能，本文采用全球降水气候计划（，简称）全球降水数据和来自欧洲中期天气预报 中 心 发 布 的 第 五 代 再 分 析 资 料（，简称）的 风场数据，分析时段选取 年。为便于与 和 的模拟结果进行评估验证，本文将观测 再分析数据利用双线性插值到与模式水平分辨率一致，即.。.中国东部南北方夏季降水指数 前人研究发现，全新世以来的不同时间尺度（年代百年千年尺度）中国东部季风区的降水变化均存在南北差异。如代用资料显示中国东部存在全新世降水南北反相变化、或至少是不同步变化的可能。在过去千年特征气候时期（中世纪气候异常期和小冰期），有学者同样发现中国东部降图 年观测 再分析资料和 试验模拟 的夏季气候态（，）降水（填色，单位：）和 风场（箭头，单位：）以及南北降水差异指数的定义范围（蓝色框线）.（，）（，：）（，：）（）水变化呈南北“偶极型”分布特征。在近百年器测时段，等利用来自 个站点 年的降水数据发现，近 年来中国东部夏季降水的年代际变化特征主要由“正负正”三级型转变为南北“偶极型”。可见，中国东部季风区降水的南北差异，即“南涝北旱”“北涝南旱”问题，受到历史学家的广泛关注。因此，本文参考陈建徽等和靳春寒的定义方法，将北方降水区域定义为 和 范围，南方降水定义为 和 范围，并将北方与南方夏季（月，简称）平均降水之差定义为南北降水差异指数，本文采用该指数来衡量 强度。.模拟结果与观测 再分析资料的对比本文使用 的降水资料和 的 风场资料来验证模式模拟的可靠性。图 为观测与模拟的 年东亚夏季降水场和 风场。从图 和 中可以看出，东亚夏季降水总体随纬度增加而减少，模拟的中国东部夏季降水空间分布和 的降水结果相似，其空间相关系数为.，达到了 的置信度；此外，模拟的 风场也与 的第 四 纪 研 究 年环流特征相似，东亚夏季受控于来自西太平洋副热带高压西侧的东南气流，中国东部盛行偏南风。因此，对东亚季风区的模拟是较为合理的。但是，模拟和观测也存在一定的差异，例如，模拟的中国东部偏南风偏强、青藏高原东侧存在虚假的降水大值中心，这些不足普遍存在于现有气候模式中，与水平分辨率等因素有关。表 本文使用的全新世东亚季风区重建资料 序号点位区域资料类型纬度经度时段（）参考文献董哥洞石笋 等，和尚洞石笋 等，乌鸦洞石笋 等，黑风洞石笋 等，狮子洞石笋 等，天门洞石笋 等，祥龙洞石笋 等，公海湖孢粉 等，玛珥湖花粉 等，莲花洞石笋 等，黄土高原黄土磁化率 等，结果分析.降水百年尺度变化特征 重建资料具体来自于表，图 为指代 降水的多条重建资料集合平均结果。由图 可知，代用指标重建的全新世 降水表现为在早全新世先上升，然后在 逐渐下降的趋势。图 为 试验和 试验模拟的全新世东亚夏季降水变化。本文将模式模拟的 降水变化（图 ）与重建资料（图）进行对比，可以发现其趋势演变相似，且我们模拟的全新世 降水变化与 等基于环流指数定义得到的结果较为一致，说明我们模拟结果的合理性和可靠性。在全新世期间，北方 月夏季降水呈下降趋势（图 ），南方降水呈增加趋势（图 ），北方降水和南方降水趋势呈现相反位相变化关系，该特征与学者们基于器测资料发现的中国东部南北方降水位相差异较为相似；而 试验（图）和 试验（图）模拟的中国东部南北方夏季降水指数均在 呈下 降 的 趋 势，且 二 者 的 线 性 趋 势 均 为.，这说明太阳活动外强迫对全新世 降水的长期趋势影响不显著。为了探究太阳活动对 降水百年尺度变化的影响，我们首先比较了在整个全新世时期，和 模拟的中国东部 降水变化的功率谱分析结果。基于 的全强迫试验（图）和地球轨道参数试验（图）的模拟结果显示，全新世 降水变化均无显著的百年尺度周期，这表明地球轨道参数、冰盖、温室气体等外强迫无法引起 降水的百年尺度变化。基于 模拟的结果显示，在 试验中，全新世 降水仅存在约 年的显著周期信号（图）；但在 试验中，在地球轨道参数和太阳活动外强迫的共同影响下，降水 年的百年周期有所加强（图），但仍不显著。我们进一步查看了 期间中国东部 降水变化的功率谱分析结果（图），发现在 试验中，中国东部 降水变化无显著的百年尺度周期，这与全新世全时段的功率谱分析结果一致（图），但 试验的模拟结果（图）表明，中国东部 降水变化存在约 年的显著周期。为了进一步分析太阳活动对 降水百年尺度变化的影响，我们将 与 试验相减结果看作净太阳活动的效应，发现功率谱结果中约 年的周期信号进一步加强（图），这说明中晚全新世以来，太阳活动外强迫对 降水的约 年的周期变化具有显著影响，这与前人基于重建资料发现的 年周期变化特征相似。.太阳活动外强迫对 降水的影响机理 为了分析 降水百年尺度变化与太阳活动外强迫的关系，在净太阳活动影响下，我们对 期陈萌萌，等：全新世太阳活动对东亚夏季风百年尺度变化影响的模拟研究图 （）全新世东亚季风区重建资料集合平均结果（具体资料见表）和（）模拟中国东部夏季北方降水 （，）、南方降水（，）以及南北降水差异（北方降水减南方降水）时间序列均为 年滑动平均结果；（）中蓝色阴影表示重建资料所有样本的 个标准差范围，（，）为 试验，（，）为 试验；单位为 .（）（）（）（，），（，），（）（），（，），（，）；：图 全新世.期间 模拟和 模拟中国东部南北降水差异指数的功率谱（，）分别为 的 试验和 试验，（，）分别为 的 试验和 试验红色虚线代表 的置信度.（）（），（）（），第 四 纪 研 究 年图 中晚全新世 期间 模拟中国东部南北降水差异指数的功率谱（），（）和（）分别为 试验、试验以及 试验减 试验；红色虚线代表 的置信度.（），（），（）降水变化与太阳活动标准化序列进行了滞后图 中晚全新世 期间净太阳活动影响下 降水（，阴影）与太阳活动标准化序列的滞后性回归分析结果来自 减 试验蓝色框线代表南北降水差异指数的定义范围，“”、“”等标签分别表示降水滞后太阳活动 年、年等，打点区域表示达到 的置信度.（，）（）“”“”，性回归分析（图）。同时为了提取百年尺度信号，本文对太阳活动序列和模式模拟结果均进行了 年带通滤波处理。结果表明，中晚全新世以来 降水百年尺度变化与太阳活动外强迫存在明显的相关关系，中国东部夏季降水的空间主模态呈偶极型分布。在太阳峰值年（年，即图），北方地区存在降水最大值中心，南方降水与北方降水存在显著差异，表现为“南旱北涝”的空间分布格局。随着太阳总辐照度逐渐减弱（图 ），北方降水由正异常转为负异常，而南方降水表现为显著减少。进一步，我们计算了中国东部南、北方夏季降水与太阳活动外强迫的相关系数（图）。结果表明，在太阳总辐照度处于极大值（年）时，北方降水与太阳活动的相关系数最高，达到.（通过了显著性检验）；此后，随着太阳总辐照度减弱（滞后 年至 年），北方降水与太阳活动相关性降低；约滞后 年期间北方降水与太阳活动的相关性不显著，这与图 所示的空间场表现 期陈萌萌，等：全新世太阳活动对东亚夏季风百年尺度变化影响的模拟研究图 中国东部南、北方降水与太阳活动外强迫的滞后相 关系数北方（蓝色虚线），南方（蓝色实线），结果来自 试验减 试验，红色圆点为相关系数的峰值，灰色虚线表示 的置信度.（）（）（）一致。然而，在太阳总辐照度处于极大值（年）时，南方降水与太阳活动呈负相关；此后，南方降水与太阳活动的负相关性逐渐增强，并在 年左右达到最大（.，通过了 显著性检验），因此，南方降水明显滞后于太阳活动约 年，这与南方地区重建资料显示的结果类似。为了探究太阳活动对 降水变化影响的机理过程，图 展 示 了 净 太 阳 活 动 影 响 下，期间东亚夏季气候要素与太阳活动标准化序列的滞后性回归分析结果，我们主要选取了 年和 年这两个滞后年份，即各气候要素在太阳活动极大值时期和滞后太阳活动 年，来分析中国东部北方降水和南方降水变化的原因。结果表明，在太阳总辐照度极大值年（年），北方地区存在显 著 的 低 压 中 心（图），上 升 运 动 显 著（图），水汽含量增加（图），北方降水明显增多。前人研究发现，海陆热力差异对季风环流和降水的变化产生影响。为此，我们计算了净太阳活动影响下，东亚季风区的陆地地表气温（，）与西北太平洋海洋区域气温（，）的温度差为.，通过了 的显著性检验，这说明太阳总辐照度的增强导致陆地升温速度超过海洋，东亚海陆热力梯度差加强，东亚陆地气压降低，引起异常的气旋性环流，这有利于增强北方降水。这一过程与前人研究发现在中世纪气候异常期东亚季风区海陆热力差异偏强，其降水呈现华北降水偏多、南方降水偏少的图 中晚全新世 期间净太阳活动影响下中国东部夏季气候要素场与太阳活动标准化序列的滞后性回归分析（，）为表面气压场（，阴影）和 风场（，矢量），仅画出通过 显著性检验的区域；（，）为 垂直速度场 （，阴影）；（，）为 比湿场（，阴影）打点区域表示达到 的置信度，蓝色框线代表南北降水差异指数的定义范围，“”、“”标签分别表示各个气候要素滞后太阳活动 年、年.（，）（，）（，）（，）（，）（，）（，）“”“”特征相似。相比于北方地区，南方地区陆地气压场和风场在太阳总辐照度极大值年均无显著变化，水汽含量第 四 纪 研 究 年变化也不明显，但受到下沉运动影响，南方降水较少。随后太阳总辐照度减弱（滞后 年），我国华北日本南部出现异常的反气旋性环流（图），西北太平洋特别是菲律宾海地区出现异常的气旋性环流（图），这使得我国南方地区出现明显的辐散，下沉运动加强（图），水汽含量变化不明显（图），导致南方降水进一步减少（图）。但是为何南方地区的降水环流变化滞后于太阳活动约 年？此时太阳总辐照度已经减弱并近似回到正常状态，显然不是太阳活动的直接效应造成的，因此我们猜想可能与太阳活动的间接效应有关，比如太阳活动通过调制海温变化进而影响 降水和环流。以往研究发现，太阳活动 年周期通过调制类 的海温模态来影响东亚季风降水的年代际变化，然而太阳活动在百年尺度上会如何影响太平洋海温变化仍不清楚。为了探究这一问题，我们选取了 年、年和 年来分析太平洋海表温度滞后于太阳活动的变化情况（图）。结果表明，在太阳总辐照度极大值时期（年），赤道东太平洋出现降温，西太平洋有一定的增暖（图），该特征与拉尼娜事件发生时的海温异常现象相似，可以被认为是“类拉尼娜”态。前人研究也提到，在太阳活动 年周期变化的峰值时期，赤道太平洋呈现类“拉尼娜”海温模态，同时也和百年尺度上中世纪气候异常期太阳活动增强影响下的“类拉尼娜”态和小冰期太阳活动减弱影响下的“类厄尔尼诺”态特征相似。在海温滞后太阳活动的 年至 年期间，赤道中西太平洋海温进一步变暖，形成异常的低压中心，位于西北太平洋气旋式环流东南侧的异常西南风与背景东北信风相叠加，削弱了气候态风场，从而减少海表面蒸发，加强了菲律宾海以东的增暖，而暖海温异常可以通过激发赤道罗斯贝波响应，进一步增强菲律宾海上空异常的气旋性环流（图 ），这与“风蒸发海温”的正反馈机制相似。而西太平洋增暖达到最大时大约处于 年，此时西北太平洋异常的气旋性环流强度也达到最强，因此对我国南方降水的抑制作用也达到最强。然而需要说明的是，我们模式中采用的是 倍加速模拟技术，可能会延长海温对大气的影响或气候系统内部变率的周期，因此，模拟的太阳活动对赤道太平洋海温影响的过程周期、南方降水滞后于太阳活动的具体年份可能存在不确定性，以后需要基于不加速的瞬变积分模拟来进一步分析。图 中晚全新世 期间净太阳活动影响下海表温度（，阴影）和 风场（，矢量）与太阳活动标 准化序列的滞后性回归分析“”、“”等标签分别表示海表温度滞后太阳活动 年、年，仅画出通过 显著性检验的区域.（，）（，）“”“”，结论与讨论本 文 基 于 模 拟 的 试 验 和试验结果，以及 模拟试验结果，探究了全新世 降水百年尺度的变化特征及其对太阳活动外强迫的响应机理，得出的主要结论如下：（）在变化趋势上，模拟的全新世以来 降水呈减弱趋势，主要受地球轨道参数的影响，太阳活动外强迫对 降水长期趋势变化没有影响。在百年尺度上，中晚全新世以来太阳活动外强迫引起 降水出现约 年的周期变化特征，这与高分辨率重建资料记录的周期相似。（）中晚全新世以来 降水百年尺度变化与太阳活动外强迫存在明显的相关性。在太阳活动百年周期的峰值年，这一百年尺度变化的空间主模态呈南北偶极型分布，表现为“南旱北涝”格局，北 期陈萌萌，等：全新世太阳活动对东亚夏季风百年尺度变化影响的模拟研究方降水与太阳活动存在显著的同期正相关，在太阳总辐照度极大值年其空间相关系数达到最高，当太阳活动增强时，东亚海陆热力差异加强，东亚陆地气压降低，上升运动显著，水汽含量增加，引起异常的气旋性环流，有利于北方降水增多。（）南方降水与太阳活动呈负相关并滞后于太阳活动，太阳活动通过调制赤道太平洋类拉尼娜态的海温梯度，赤道中西太平洋海温增暖，形成异常的低压中心，位于西北太平洋气旋式环流东南侧的异常西南风，通过削弱气候态东北信风令蒸发减弱，加强了菲律宾海以东的增暖，而暖海温异常可以通过激发赤道罗斯贝波响应，引起西北太平洋风蒸发海温的正反馈机制，触发异常的西北太平洋气旋性环流，这导致南方地区出现明显的辐散，下沉运动加强，进一步减少南方降水。由于西太平洋增暖达到最大时大约处于 年，此时西北太平洋异常的气旋性环流强度也达到最强，因此对我国南方降水的抑制作用也达到最强。然而当前研究仍存在以下的不确定性：）本文仅是 气候模式的模拟结果，可能具有一定的模式依赖性，同时模式中采用的是 倍加速模拟技术，也 会 存 在 一 定 的 不 确 定 性；）本 文 将试验与 试验相减的结果作为净太阳活动效应，探究其对 降水百年尺度变化的影响，考虑到本研究主要是基于长时间尺度的变化，我们认为将 试验与 试验相减能够在一定程度上消除内部变率的影响进而分离出太阳活动的净效应，且该结果与 试验模拟结果相似（图略），这表明 降水百年尺度变化主要受太阳活动线性主导，但也会存在其他非线性因子的影响；）代用资料的地理位置、时间覆盖范围等各不相同，其指示意义存在不确定性，将代用资料与模拟结果进行对比验证时可能具有一定的不确定性。因此在未来的工作中需要加强多模式集成研究，有待于基于不加速的瞬变积分模拟进一步分析，并且在今后的研究中考虑其他非线性因子的影响，以及采用更多高分辨率的重建资料结果进行对比验证。此外，针对太阳活动如何影响太平洋海表温度变化，以及太阳活动千年尺度变化周期对东亚季风的影响方面等问题，在未来仍须进一步深入分析。致谢：感谢两位审稿专家和编辑部杨美芳老师提出的宝贵意见和建议！参考文献（）：安芷生，刘晓东 东亚季风气候的历史与变率 科学通报，（）：，（）：，（）：张志平，黄伟，陈建徽，等 全新世东亚夏季风演化的多尺度周期变化及其可能机制探讨 第四纪研究，（）：，（）：，,，（）：，：，,（）：，,，（）：，（）：王梦泽 全新世时期中国北方地区气候变化文献综述 中国资源综合利用，（）：，（）：，孙炜毅，刘健，严蜜，等 全新世亚洲季风百年千年尺度变化的模拟研究进展 地学前缘，（）：，：，（）：万凌峰，刘健，高超超，等 全新世火山喷发对温度变化趋势影响的模拟研究 第四纪研究，（）：，（）：，：，：，第 四 纪 研 究 年 ，,（）：，,（）：，,（）：，（，），（）：余英浩，金映豫，徐德克，等 青海可鲁克湖孢粉记录的.以来植被和气候演化历史 第四纪研究，（）：，（）：冷程程，赵成，崔巧玉，等 大兴安岭阿尔山天池沉积物中正构烷烃记录揭示的全新世古气候变化 第四纪研究，（）：，（）：，（）：，,（）：谭明 环流效应：中国季风区石笋氧同位素短尺度变化的气候意义 古气候记录与现代气候研究的一次对话 第四纪研究，（）：，（）：，（）：马春梅，朱诚，郑朝贵，等 晚冰期以来神农架大九湖泥炭高分辨率气候变化的地球化学记录研究 科学通报，（）：，（）：，（）：，,（）：，：，,：，：，,（）：，（）：，（）：，,（）：，（）：，,：，,（）：，,（）：戴张奇，刘健，孙炜毅，等 过去 年 个典型暖期北半球季风降水的模拟研究 第四纪研究，（）：，（）：李新周，刘晓东 早全新世与未来 后东亚区域气候变化对比：自然强迫和人类活动的影响 第四纪研究，（）：，：，（）：靳立亚，陈发虎 千百年尺度气候快速变化及其数值模拟研究进展 地球科学进展，（）：，期陈萌萌，等：全新世太阳活动对东亚夏季风百年尺度变化影响的模拟研究 ，（）：，：，（）：，,：，,，（）：，,，（）：，,（）：，：，（）：，：，（）：，（）：，：刘珊珊，姜大膀 过去 中国干湿变化的瞬变模拟分析 第四纪研究，（）：，（）：靳春寒 太阳活动对亚洲季风年代际气候变化的影响研究 南京：南京师范大学博士学位论文，：，：，（）：周天军，李博，满文敏，等 过去千年 个特征期气候的 耦合模式模拟 科学通报，（）：，（）：李倩倩，曹剑，赵海坤 过去两万年南北半球季风降水反位相变化特征研究 第四纪研究，（）：，（）：谭亮成，蔡演军，安芷生，等 石笋氧同位素和微量元素记录的陕南地区 高分辨率季风降雨变化 第四纪研究，（）：，（）：，,（）：孙炜毅，刘健，万凌峰，等 全新世北半球中纬度降水变化对不同外强迫响应的模拟研究 第四纪研究，（）：，（）：，：，（）：，：，：，,（）：，（）（），（）：，等 全球再分析产品 气象科技进展，（）：，（）：，,，（）：，（）：，：，,：第 四 纪 研 究 年 ，（）：周秀骥，赵平，刘舸，等 中世纪暖期、小冰期与现代东亚夏季风环流和降水年代百年尺度变化特征分析 科学通报，（）：，（）：，：，（）：黄荣辉，徐予红，周连童 我国夏季降水的年代际变化及华北干旱化趋势 高原气象，（）：，（）：，：，（）：，（）：钱维宏，朱亚芬，汤帅奇 重建千年东亚夏季风干湿分布型指数 科学通报，（）：，（）：邢辰，孔瑜，刘飞，等 重建与模拟中过去 年厄尔尼诺和东亚夏季风关系的年代际变化 第四纪研究，（）：，（）：陈建徽，饶志国，刘建宝，等 全新世东亚夏季风最强盛期出现在何时？兼论中国南方石笋氧同位素的古气候意义 中国科学：地球科学，（）：，？：，（）：，,：，,：，（）：，,：，,，（）：程军，张近童，刘一陶，等 轨道尺度下东亚夏季风降水的北方主导演变模态 地球环境学报，（）：，（）：周雪城，姜大膀 过去千年年代和百年尺度东亚夏季风变化的模拟分析 第四纪研究，（）：，（）：黄平，周士杰 全球变暖下热带降水变化研究回顾与挑战地球科学进展，（）：，（）：，,（）：，（）：，：，,（）：，,（）：，,（）：，：？，（）：期陈萌萌，等：全新世太阳活动对东亚夏季风百年尺度变化影响的模拟研究 ，（，；，；，；，；，）（）（），：；