基于耦合气候模式的始新世全球季风数值模拟

郭衍游，王成善，杨 杰，张 健，蔡惠慧

（1.中国地质大学 地球科学与资源学院，北京100083；2.生物地质与环境地质国家重点实验室（中国地质大学），武汉430074；3.成都理工大学 地球科学学院，成都610059）

季风是全球气候系统的一个重要组成部分，是全球大气能量和水汽循环最显著和最重要的区域，全球超过一半人口生活在季风区［1］。季风是国际“气候变异与可预报性研究（CLIVAR）”计划、“全球能量和水循环（GEWEX）”计划和“世界气候研究（WCRP）”计划的重要组成部分，是气候研究领域最重要的科学问题之一。

当前，全球变暖是国际社会共同关注的热点问题。随着CO2排放量的持续增加，地球可能进入温室气候时期。在温室气候条件下，全球季风会有哪些显著变化成了人们极为关心的问题。古气候研究为这一问题提供了很好的途径。目前的研究多集中于现代季风格局的形成及其与青藏高原隆起等因素的关系上，对更早的温室地球时期的季风特征研究很少。地质证据表明，地球历史上最强盛的季风正发生于温室气候时期。在二叠纪至早侏罗世的温室地球时期，全球范围内存在强盛的“超级季风（megamonsoon）”［2，3］。研究结果表明，这种超级季风主要与当时的海陆地形有关。当时，北半球的劳亚大陆和南半球的冈瓦纳大陆在赤道附近连接拼合形成联合古陆，形成大范围强烈的海陆差异，雨量集中在特提斯洋附近，联合古陆内部降雨量几乎为零［2］。联合古陆解体后，地球仍处于温室气候时期，数值模拟结果表明，白垩纪时期东亚地区的盛行风向有明显的季风反转，表现出季风特性［4］。但是，由于温室时期地质记录相对较少、距今时间久远，人们对温室气候下全球季风的范围、强度和变化机制等知之甚少，亟待从地质记录和数值模拟等方面加强研究工作［5］。

耦合气候模式考虑了不同气候子系统之间的相互作用，是当前气候模式研究发展的方向。本文利用目前最先进的耦合气候模式（Community Earth System Model，简称CESM），模拟始新世（40Ma B.P.）地球气候环境，并与现今气候作对比，探讨始新世这一典型温室气候时期全球季风的分布特征。

1 现代全球季风

传统上，季风被定义为盛行风向和降水的季节性反转，是低纬度地区的主要气候系统［1］。现在一般认为季风还表现为整个大气环流格局发生季节性反转［6］。季风由海陆热力差异及行星风带的季节性差异所引起。B.Wang和Q.H.Ding根据夏、冬两季的降水差异，定义了一个季风降水指数［7］（Monsoon Precipitation Index，MPI）

根据MPI值，将现代季风划分为6个季风系统：北美季风、南美季风、北非季风、南非季风、亚洲季风和澳大利亚－印尼季风（图1）。其结果与“全球降水气候学计划”（Global Precipitation Climatology Project，简称GPCP）的全球季风区降水率差异＞3.75×10－8m/s的区域相符（图2）。

全球范围内，降水率的季节差异很好地反映了季风的分布和强度。因此，本文主要采用降水这一指标来判断季风的存在，将中低纬度区域（30°N～30°S）干、湿季降水率差异＞3.75×10－8m/s的地区划为季风区。

2 模式设置

图1 全球季风Fig.1 The global monsoon（据文献［7］修改，转引自文献［5］）

图2 GPCP观测的全球干、湿季降水率差Fig.2 The precipitation difference between the wet season and the dry season by GPCP observations

CESM数值模式是在美国国家大气科学研究中心主导下，集中全球科学家共同开发和维护的气候模式。该模式从CCSM（Community Climate System Model）演变而来，其目标是将地球表层圈层视为一个系统，通过耦合技术，模拟大气圈、水圈、生物圈的相互作用，从而预测现代气候变化或重现地质历史时期古气候/古环境［8，9］。本文的模拟采用的CESM模式版本号为1.0.3，包含大气、海洋、陆地、陆冰和海冰等5个模块，各模块之间通过耦合器（coupler）来实现质量、能量的交互传输。其中，大气模块为NCAR开发的CAM（community atmosphere model），采用 T31谱模式，空间网格约为2°×2°，大气分层为26层，垂直坐标为sigma－p混合坐标；海洋模块为改进的POP（parallel ocean program）模式，空间网格约为1°×1°，海水分层为60层，垂直坐标为z坐标；陆地模块为CLM 2.1（Community Land Model version 2.1），植被参数根据 Dickinson等的统计数据按纬向分布设置［10］。更详细的模式信息请参照CESM网页。

模式中地形和水深数据来自美国PALEOMAP Project的重建资料（图3）［11］。根据浮游有孔虫的硼同位素测量结果，40Ma B.P.大气CO2体积分数为0.8‰～1.4‰［12］，因此，模式中大气CO2体积分数设为工业革命前的4倍，为1.120‰。以该模式模拟了1.2ka，取最后30a的模拟结果进行分析。

3 结果

3.1 始新世温室气候背景

地质记录表明，始新世时期全球地表温度明显高于现今，特别是极地地区比现今温暖得多，温带森林已经扩展到了极地地区［13］。中国主要陆块为热带、亚热带植物覆盖，气候温暖湿润［14］。

本次模拟结果显示，总体上当时全球年平均气温比现今高4.28℃（表1）。其中，低纬度地区年平均温度为27.61℃，比现今温度高2.32℃；但是在南半球高纬度地区，温度高出接近20℃。在局部地区存在温度比现今低的情况，在北纬60°～65°之间，温度比现今低0.1～0.5℃，这可能与局部热量输送以及植被分布有关（图4）。

将模拟的始新世（40Ma B.P.）全球表面气温分布与美国国家环境预报中心（NCEP）观测资料对比，可以进一步看出始新世与现今全球表面气温的变化情况。图5是模拟的始新世（40Ma B.P.）气温值减去NCEP观测的气温值的全球分布图，从图上可以看出，高值主要分布在南极大陆、青藏高原和格陵兰岛，局部升温达30℃。南极大陆和格陵兰岛由于地处高纬度，在地球升温过程中，受地球经向热量输送加强的影响，易于升温，且升温幅度大。青藏高原地区的升温则是由于始新世时期青藏高原还未隆起，其地面气温比现今高原地面气温要高。在赤道及低纬度地区，在温室气体含量高的条件下，热量增加不会直接导致温度升高，而是通过水汽蒸发、对流增强的方式，将热量输送至高纬度地区。在亚欧大陆以及赤道大洋等部分地区，温度比现今低1～4℃，局部低5～6℃。究其原因，亚欧大陆由于其广阔的陆地，在北极冷空气影响下，受全球升温影响较小，表现了与全球性升温不一样的气温变化。赤道大洋局部地区的温度降低则可能与强烈的对流活动有关。与前人的温室气候模拟结果相比，本次模拟出赤道地区总体上有明显的升温（图5），这与地质记录相一致［15］。

图4 始新世模拟年平均气温值与日本气象厅25a（1979—2004年）观测值Fig.4 The simulated annual air temperature of Eocene and the Japanese 25-year reanalysis（JRA-25）observations for the period from 1979to 2004

图5 始新世地表温度与现今观测值的差Fig.5 The difference between the simulated Eocene surface temperature and the modern surface temperature

模拟的降水量纬向年平均值与GPCP观测的现今降水量相比，明显偏高。这是由于温室气候条件下，水汽蒸发强烈，降水也随之增加。因而，始新世的全球年总降水量也高于现今（图6）。

3.2 季风区的分布及特征

图6 始新世模拟降水量与GPCP观测值Fig.6 The simulated annual precipitation of Eocene and GPCP observations

通过对比始新世（40Ma B.P.）与GPCP观测的干、湿季降水分布（6—8月份平均降水率减去12月份至次年2月份平均降水率，图7、表2），可以看出：始新世时期，全球低纬度地区存在明显的季风气候特征，以干、湿季降水率差异3.75×10－8m/s为阈值，也可以划分为6个大季风区。其中由于海陆地形的不同，印尼－澳大利亚季风被古印度季风所取代，其他季风区的位置与现今基本一致。从区域上看，各季风区的范围及强度有所不同。始新世时，北非季风范围和强度比现今大；北美和南美季风则比现今弱；澳大利亚地处高纬度地区，季风特征不明显，可以认为不存在季风；亚洲季风比现今弱，这与当时青藏高原还未全面隆起有关。

从统计的季风区分布面积以及降水季节差异反映的季风强度来看（表2），始新世时期北美季风和北非季风的分布面积与现今有很大差异，全球季风区总面积比现今多20%，季风强度与现今大体相当。

从各季风区的降水来看，始新世时期降水沿纬向分布的特征较明显，北半球的亚洲季风区与北非季风区相连，在低纬度地区形成了广阔的季风带，是当时最重要的季风区；而现今亚洲季风区，季风带来的降水主要沿经向分布于东亚地区。

模拟的盛行风表明，季风区盛行风的季节变化明显（图8）。6—8月份，亚洲东部及南部大部分以南风为主，12月份至次年2月份，则为东北风或西风。非洲北部6—8月份吹南风或西南风，12月份至次年2月份冬季则为东北风。

4 讨论

图7 全球干、湿季降水差异Fig.7 The global precipitation difference between the dry season and the humid season

表2 始新世与现今季风区对比Table 2 Comparison between the Eocene monsoon and the modern monsoon

图8 模拟的始新世全球表面风速Fig.8 The simulated Eocene global surface wind velocity

全球季风的变化受大陆分合的“威尔逊旋回”、地球运行轨道、太阳活动周期、全球冰量（ice volume）等多种因素调控［5，16］。由于资料缺乏，目前的古季风研究着重探讨现代格局的季风的产生时限、发展变化及其影响因素等，对于新生代早期或“深时（deep time）”季风的研究很少。虽然现代季风的起源被认为出现很晚，但是，季风本质上是由于行星风带的季节性偏移和海陆热力差异所引起的，而这两个因素在地质历史上一直存在；因此，季风是贯穿整个地质历史的永恒现象，只有具体某一区域的季风系统，才可以谈论其形成时间和消亡历史［5］。本次模拟表明，在始新世这一温室时期，就有类似于现今的全球季风存在，其范围及强度与现今气候大体相当。现代气候研究表明，在当前全球变暖背景下，印度夏季季风趋于增强［17］，东亚夏季季风则只发生位置的偏移而不是强度的增强［18］。在始新世时期，出现与现今相似的全球季风是合理的。

5 结论

a.始新世时期，全球平均温度、降水量高于现今，为典型温室气候。

b.全球存在季风，季风分布范围及强度与现代气候相当，但是区域上季风特征有明显差异。

c.亚洲季风沿纬向分布，表现出与现今东亚季风沿经向分布的不同特征。

d.澳大利亚地处高纬度区，季风特征不明显。

e.北非季风比现代强，范围大。

f.北美和南美季风比现今弱。

由于缺乏更精确的古地形、古植被、古CO2含量等数据以及更高分辨率的模拟，目前还难以将模拟结果与地质记录进行准确对比。要弄清地质历史中全球季风的动力学机制和演变规律，还需要在地质记录、数值模拟和理论分析等方面加强研究。

［1］Ramage C.Monsoon Meteorology［M］.San Diego：Academic Press，1971.

［2］Kutzbach J E，Gallimore R G.Pangean climates：megamonsoons of the megacontinent［J］.J Geophys Res，1989，9：3341－3357.

［3］Parrish J T.Climate of the supercontinent Pangea［J］.J Geol，1993，101：215－233.

［4］Chen J M，Zhao P，Wang C S，etal.Modeling east asian climate and impacts of atmospheric CO2concen-tration during the Late Cretaceous（66Ma）［J］.Palaeogeogr Palaeoclimatol Palaeoecol，2013，385：190－201.

［5］汪品先.全球季风的地质演变［J］.科学通报，2009，54（5）：535－556.Wang P X.Global monsoon in a geological perspective［J］.Chinese Sciences Bulletin，2009，54（5）：535－556.（In Chinese）

［6］Trenberth K E，Stepaniak D P，Caron J M.The global monsoon as seen through the divergent atmospheric circulation［J］.J Clim，2000，13：3969－3993.

［7］Wang B，Ding Q H.Global monsoon：Dominant mode of annual variation in the tropics［J］.Dyn Atmos Ocean，2008，44：165－183.

［8］Vertenstein M，Craig T，Middleton A，etal.CESM1.0.3User's Guide［M］.Boulder：NCAR，2011.

［9］CESM home.NCAR CESM（Community Earth System Model）earth· modeling·climate［EB/OL］.［2014-01-01］.http：//www2.cesm.ucar.edu/.

［10］Dickinson R E，Oleson K W，Bonan G，etal.The community land model and its climate statistics as a component of the community climate system model［J］.J Clim，2006，19：2302－2324.

［11］Scotese C R.PALEOMAP project［EB/OL］.［2004-02-20］.http：//www.scotese.com/climate.htm.

［12］Pearson P N，Palmer M R.Atmospheric carbon dioxide concentrations over the past 60million years［J］.Nature，2000，406：695－699.

［13］Huber M，Caballero R.The early Eocene equable climate problem revisited［J］.Clim Past Discuss，2011，7：241－304.

［14］王晓梅，王明镇，张锡麒，等.中国晚始新世－早渐新世地层孢粉组合及其古气候特征［J］.地球科学，2005，30（3）：309－316.Wang X M，Wang M Z，Zhang X Q，etal.Palynology assemblages and paleoclimatic character of the late Eocene to the early Oligocene in China［J］.Earth Science，2005，30（3）：309－316.（In Chinese）

［15］Huber M.Snakes tell a torrid tale［J］.Nature，2009，457：669－671.

［16］An Z S，Steven C.Shen J，etal.Glacial-interglacial Indian summer monsoon dynamics［J］.Science，2011，333（6043）：719－723.

［17］Hu Z Z，Latif M，Roeckner E，etal.Intensified Asian summer monsoon and its variability in a coupled model forced by increasing greenhouse gas concentrations［J］.Geophys Res Lett，2000，27：2681－2684.

［18］Li J P，Wu Z W，Jiang Z H，etal.Can global warming strengthen the East Asian summer monsoon？［J］.Climate，2010，23：6696－6705.