近一个世纪以来，气候变化已对当前生态系统的分布产生了影响。识别及分析影响区域的位置及特征，对于有针对性的制订气候变化下的生态系统保护策略有着重要的作用。

目前，已有的研究表明，使用气候—植被分类模型识别气候变化下生态系统可能发生变化的区域是一种有效的方法^{[1, 2, 3]}。气候—植被模型主要包括了机理模型和非机理模型。机理模型虽然能够较好的将气候变化以及植被响应进行预测，但是此类模型所需数据不易获取，并且气候变化对于生态系统的影响机制尚不完全明确^{[4, 5]}。使得运用非机理性模型评价气候变化对于生态系统的影响仍然是一种有效方法^{[2, 6]}。自德国植物学家柯本建立生物气候分类法以来，生物地理学家已建立了多个非机理性的气候-植被模型，其中使用最为广泛的是柯本气候分类法^[7]。已有大量研究使用此分类法从全球和区域尺度评估气候变化对于植被的影响并且取得了较好的效果^{[2, 8, 9, 10]}。对于中国区域，虽然使用柯本气候分类法评估结果较少^[2]，但仍有一些运用其它气候-植被模型评估气候变化影响的研究,如张新时^[6]，Chen等^[11]，岳天祥等^[12],胥晓等^[13]。然而已有的研究多是基于较大尺度的分析，难以体现气候变化对于山地生态系统的影响。此外，由于山地区域往往海拔变化大，使得不同海拔区间内气候变化的程度存在明显差异，导致对不同海拔区域内生态系统的影响也并不相同。因此，分析不同海拔区间内生物气候类型的变化特征，对于系统的制订保护策略十分重要，但目前的研究并未对此进行详细分析。

随着薄盘样条气候表面插值方法的建立^[14]，能够在较小尺度上获取较高精度的气候表面数据。在此基础上，通过运用柯本气候分类法，可以评估气候变化对于山地生态系统的影响。因此，本研究在区域气候模型(RegCM4.0)对四川省未来气候变化进行预测的基础上，利用薄盘样条插值方法对气候数据进行处理。并使用柯本气候分类法评估未来气候变化对四川省生态系统的影响。同时也对气候变化下各生物气候类型变化区域的海拔特征进行分析。

四川省位于青藏高原东缘，处于中国地势三大阶梯中的第一级和第二级，区域内高差悬殊，地势西高东低，是中国典型山地分布区。四川省独特的地理位置，使其一方面易受气候变化影响，另一方面形成了丰富的生态系统类型。因此，对气候变化下四川省生态系统分布变化的评估，可以较好的反映气候变化对于山地生态系统的影响。

虽然有大量研究对柯本气候分类法进行了修正，但是目前使用比较广泛的分类方法是Kppen-Geiger分类方法(以下使用柯本分类法代表此类方法)。由于已有相关文献对柯本气候分类法进行了详细说明^[1]，因此本文仅作简要描述。按柯本气候分类法分类原则，将全球植被划分为5个主要类型，再根据季节性降水和温度在每个主要类型中划分亚类型。具体划分时，首先根据月最低温，月最高温，月降水划分出5个主要类型，分别为赤道气候带(A)，干燥气候带(B)，暖温带气候带(C)，冷温带气候带(D)，极地气候带(E)，再根据冬季和夏季的降水，月平均温进一步的划分亚类型。A类型中包括的亚类型为常年湿润雨林(Af)，季雨林(Am)，夏干型稀树草原(As)，冬干型稀树草原(Aw)。B类型包括的亚类型为干旱型草原气候(BS)，干旱型沙漠气候(BW)。C类型包括的亚类型为冬湿夏干型暖温带气候(Cs)，夏湿冬干型暖温带气候(Cw)，全年常湿型暖温带气候(Cf)。D类型包括亚类型为冬湿夏干型暖温带气候(Ds)，夏湿冬干型暖温带气候(Dw)，全年常湿型暖温带气候(Df)。E类型包括亚类型为苔原气候(ET)，冰原气候(EF)。其中在C类和D类中又可根据季节性温度，在每个亚类型中再划分出夏季炎热型(a)，夏季温暖型(b)，夏季凉爽型(c)和显著大陆型(d)。

本研究使用区域气候模式RegCM4.0，单向嵌套BCC_CSM1.1全球气候系统模式，在中等排放情景rcp4.5下，以50 km×50 km的尺度，模拟1981—2099年逐月的平均气温及降水。在此基础上使用204个省内外站点(包括四川省内119个气象站点以及周边区域85个气象站点)数据与模拟数据进行对比，去除预估数据的系统偏差。

在气象数据校正的基础上，使用薄盘样条插值软件Aunsplin 4.36进行插值。在插值时，选择插值模型以及样条次数的不同可能对插值结果产生一定影响，而本研究又需要对较长时间序列的气象数据进行插值。因此需要选取的模型既要保证表面精度又要考虑模型的稳定性，使得不同月份间可以相互比较。在进行气候表面生成前，分别选取薄盘样条以及局部薄盘样条两个插值函数，以海拔作自变量或协变量，在样条次数为2，3，4，5次条件下进行插值。在Anusplin中，可通过对比不同插值函数以及不同样条次数生成的广义交叉验证(GCV)或是最大似然法(GML)值的大小确定最佳插值函数及其样条次数。通过对比，本研究选取局部薄盘模型，在样条次数为2次，海拔为协变量的条件下对四川省1981—2099年逐月平均气温和降水空间分布进行插值，插值尺度为1 km×1 km。

根据世界气象组织(WMO)对于气候的定义，应选取至少30a为1个气候周期。同时，考虑到气候变化的不确定性，本研究选取1981—2010年，2011—2040年，2031—2060年和2070—2099年4个气候时段分别代表当前(2010)，近期(2040)，中期(2060)以及远期(2099)的气候时段。按照柯本气候分类法的划分原则，利用Arcgis软件与Python语言对这4个气候时段内的柯本气候类型分布进行划分。在划分时，本研究使用1个气候时段内各月平均温度，月降水，各年最暖(冷)月温度的平均值分别代表柯本气候分类法中的月平均温度，月降水，最暖(冷)月温度。

本研究所使用逐月气温与降水的站点数据来源于四川省气象局。

为检验模拟数据与站点数据间的差异，随机选取4个不同海拔范围内的气象站点(图2)，对1981—2010时段逐月气温与降水的实测值与模拟值进行回归分析。结果表明，选取的4个站点中，气温的模拟值与实测值间差异较小，而降水预测值与实测值间差异明显(图1)。由于柯本气候分类法主要是通过气温对生物气候类型进行划分。因此，使用模式输出结果是可以满足分析要求。

图 1 预估数据与气象站点数据对比 Fig.1 The difference between projection data and observations data a,b,c,d代表不同气像站点，1代表月平均气温，2代表月降水量

图选项

图 2 四川省1981—2010年月平均气温空间分布图(1 km×1 km) Fig.2 The mean temperature with 1 km spatial resolution in 1981—2010 period

图选项

如图2所示，使用薄盘样条插值方法将站点数据插值到1 km×1 km尺度是可以较好地反映地形对于区域气候的影响。为了验证插值结果与站点数值差异是否显著，本文以100 m的海拔间隔为标准，按照海拔梯度将四川省划分为55个区域，其中第55个区域代表海拔5500 m及以上的区域。选取1981—2010时段，分别对插值生成的气候表面数据以及预估的气象站点数据在各海拔区间内的月平均温度进行统计，并进行回归分析。由于并不是所有海拔区间均有气象站点分布，因此本研究仅选取有站点分布的海拔区间进行分析。分析结果表明，在选取的时段内气候表面数据与站点的预估数据有较好的一致性(图3)。

图 3 1981—2010时段气候表面平均温度与预估站点平均温度对比 Fig.3 Difference between projection temperature data and climate surface data

图选项

运用柯本气候分类法对四川省生物气候类型分布的识别结果表明，在当前2010时段，四川省共包括3个柯本气候类型(C类，D类和E类)和11个气候亚类型(图4)。各类型及亚类型详见图4的图例。在当前时段内，分布面积最大的类型为C类型，其分布面积占四川省总面积的54%，分布面积最小的是D类型，其分布面积占四川省总面积的22%。亚类型中分布面积最大的是ET类型，分布面积占四川省总面积的24%。根据各生物气候类型分布的空间特征，从东到西的类型分别为C类，D类和E类。相比川东平原及丘陵地区，川西高原的生物气候类型空间异质性明显较高。

图 4 气候变化情景下四川省柯本气候区分布变化 Fig.4 The distributions of Kppen climates in Sichuan Province under climate change

图选项

由于柯本气候分类法缺少亚热带生态系统的分类，一般认为C类型反映了亚热带和暖温带两个植被类型。在识别出的四川省柯本气候类型中C类型，D类型和E类型分别对应亚热带针阔混交林生态系统，寒温带针叶林生态系统以及高山草甸生态系统。因此，本研究的识别结果与四川省实际的生态系统类型的分布基本相符。并且，使用1 km×1 km的尺度进行分析，也是能够反映较小区域内地形对于植被分布的影响。

在气候变化情景下，四川省柯本气候类型空间分布特征无明显变化(图4)。相比当前2010时段，在未来2040，2060以及2099时段面积增加的气候类型为C类型和D类型，至2099时段这两个类型增加面积占当前分布面积的13%和20%。面积减少的类型为E类，至2070—2099时段E类型减少面积占当前分布面积的48%。为准确的反映气候变化情景下不同生物气候类型的相互转变，本研究使用转移矩阵方法对不同时段生物气候类型间的转变进行统计，结果见表1—表3。表中各行代表某个类型向其它类型转出的面积，各列代表其它类型向某类型转入的面积。表1—表3的结果表明，在气候变化情景下，除了ET亚类，其它亚类型间均有转出和转入，而ET亚类仅有转出。在各柯本气候类型的转变中，最大的转变是由E类型向D类型转变，到2099时段，有55220 km²的E类型转变为了D类型。对比未来不同时段的转变速率表明，近期时段的气候变化造成的气候类型的转变速率要高于中期和远期时段的气候变化。

为识别气候变化下各生物气候类型变化区域的海拔特征。本研究以当前时段为基准，分别与未来3个时段进行对比，识别出各时段内生物气候类型变化区域的空间分布。同时，以100 m的海拔间隔为标准，使用Arcgis软件对各海拔区域内柯本气候类型的变化面积进行统计。由于降水数据存在较大的不确定性，因此本文仅统计由温度产生的变化。结果表明，从海拔900—4900 m区域内均有柯本气候类型的转变，但是相比低海拔区域，高海拔区域的转变面积明显更多(图5)。根据各类型的空间分布，在高海拔区域(川西地区)的转变主要是由E类型向D类型转变，其中发生转变的主要海拔范围是3800—4500 m，转变面积最多的海拔区域是4200—4300 m。在低海拔区域主要是由D类型向C类型转变，以及C类型中各亚类型相互转变，其中发生转变的主要海拔范围是1200—2000 m，转变面积最多的海拔区域是1300—1400 m。对比未来不同的时段，随时间增加各海拔区域内生物气候类型转变的面积逐渐增大。

图 5 不同时段柯本气候类型变化面积随海拔分布特征 Fig.5 The changed area for different time periods in elevation zones of 100m interval

图选项

此外，随着气候变化下各生物气候类型的相互转换，使得各类型分布海拔发生变化。为反映不同类型的海拔变化，分别对未来不同时段各类型分布的平均海拔进行统计，并与当前时段进行比较。结果表明C类，D类和E类分布平均海拔分别上升258，307 m和157 m，平均上移速率分别为2.9，3.4 m/a和1.8 m/a。

相比中国其它区域的评估结果^{[11, 12]}，本研究中并未出现新增加的生物气候类型，也没有将消失的生物气候类型。但是，在本研究的气候情景下，未来气候变化的程度足以使得四川省的生物气候类型发生转变，并且转变的主要趋势为由温度低的类型向温度高的类型转变。对于生物气候类型转变的区域，川西高原发生转变的面积明显多于川东丘陵和平原地区。通常高海拔区域的升温程度要高于低海拔区域。因此，分布于高海拔的E类型相比其它类型对气候变化更敏感，发生转换的面积也更多。本研究中，生物气候类型的转变主要集中于E类型向D类型的转变，这一转变趋势与大多数选取中国区域研究气候变化对生物气候类型影响的结果相一致^{[11, 15]}。但是，相比其它区域，未来四川省柯本气候E类型向D类型转变速率相对较低^[16]。根据RUBEL等提供的使用RegCM 3.0模型对不同排放情景下柯本气候类型的分布数据进行分析^[3]，至2076—2100年四川省E类气候区分布面积相比1976—2000年，在最高的排放情景下可能减少91%，在最低的排放情景下可能减少66%，高于本研究中E类型减少48%的结果。形成差异的原因可能是由于(1)本研究使用了更小尺度的气候插值方法。(2)使用气候预测模型RegCM 4.0与其它研究使用RegCM 3.0模型对气候预测存在差异。由于本研究使用气候情景为中等排放情景，E类气候区的变化应在66%—91%范围内。因此，更小尺度的分析是决定气候类型转变的程度的最主要因素。而根据气候变化对川西高原地区雪线以及四川省周边地区树线影响的研究结果表明，过去的气候变化已使得这些区域的雪线，树线向高海拔迁移，但并未对此区域内当前生态系统格局产出明显影响^{[17, 18, 19]}。雪线和树线的上移在一定程度上可以反映E类型向D类型的转变。因此，选择较高精度的气候插值方法分析气候类型的转变结果可能与实际更相符合。由于降水模拟的结果较差，本文不对降水产生的生物气候类型的转变进行讨论。

气候变化下，四川省所有生物气候类型都向高海拔区域迁移。迁移速率由大到小的顺序为D类>C类>E类。对比Lenoir等通过选取171个物种在海拔0—2600m范围内研究气候变化下物种的迁移结果表明，本研究中分布于相似海拔的C类型的迁移速率与Lenoir的研究结果基本一致^[20]。此外，对于气候变化影响比较严重的海拔区域，主要是集中于1200—2000m和3800—4500m区域内。这两个海拔区域内的变化，分别占未来三个时段全部变化的62%，63%和65%。但是不同于低海拔区域以亚类型转变为主，高海拔区域的转变主要为类型的转变。因此，对于生态系统的保护，首先是要排除这些可能受气候变化影响比较严重的区域。然而，若是这些区域内存在重要的物种，则需要对这些区域制订更为严格的保护措施，减少其它干扰因素。因为生态系统的变化是由多个因素共同决定的，虽然气候变化是不可阻止的，但减少其它干扰因素也能在一定程度上维护生态系统的稳定。根据本文的分析结果，在气候变化下对四川省的生态系统进行保护，需针对不同的海拔区域制订相应策略。同时，对近期气候变化造成的影响需要更多的关注。

通过对未来气候变化情景下四川省柯本气候类型分布变化的分析，可以得出以下结论：

(1)在未来气候变化下，四川省柯本气候类型分布格局虽无明显变化，但是未来气候变化的升温程度足以使得四川省柯本气候类型发生变化。在气候变化下，C类型和D类型分布面积增加，E类型分布面积减少。并且在近期气候变化下各类型转变的速率要高于远期的气候变化。通过与大尺度研究结果对比，在本研究所选取的尺度下生物气候类型变化速率与实际更相符合。

(2)在未来气候变化下，高海拔区域受气候变化影响明显高于低海拔区域，其中变化主要集中于海拔3800—450m的范围内。同时，从近期到远期的气候变化，所有生物气候类型分布海拔均呈现上升趋势。相比1981—2010，到2070—2099年C类，D类和E类型平均分布海拔向高海拔区域上移的速率分别为2.9，3.4 m/a和1.8 m/a。因此，使用柯本分类法，是可以在一定程度上反映气候变化下的生态系统类型的垂直迁移。