三江源区生态系统和土壤保持服务对未来气候变化的响应特征

ＡＣＴＡＥＣＯＬＯＧＩＣＡＳＩＮＩＣＡ

生态学报

Ｖｏｌ．４０，Ｎｏ．１８Ｓｅｐ．，２０２０

ＤＯＩ：１０．５８４６／ｓｔｘｂ２０１９０３３００６１７

冯晓玙，黄斌斌，李若男，郑华．三江源区生态系统和土壤保持服务对未来气候变化的响应特征．生态学报，２０２０，４０（１８）：６３５１⁃６３６１．ｒｅｇｉｏｎ，Ｃｈｉｎａ．ＡｃｔａＥｃｏｌｏｇｉｃａＳｉｎｉｃａ，２０２０，４０（１８）：６３５１⁃６３６１．

ＦｅｎｇＸＹ，ＨｕａｎｇＢＢ，ＬｉＲＮ，ＺｈｅｎｇＨ．ＲｅｓｐｏｎｓｅｃｈａｒａｃｔｅｒｉｓｔｉｃｓｏｆｅｃｏｓｙｓｔｅｍｓａｎｄｓｏｉｌｃｏｎｓｅｒｖａｔｉｏｎｓｅｒｖｉｃｅｓｔｏｆｕｔｕｒｅｃｌｉｍａｔｅｃｈａｎｇｅｉｎｔｈｅＴｈｒｅｅ⁃Ｒｉｖｅｒｓ

三江源区生态系统和土壤保持服务对未来气候变化的响应特征

冯晓玙１，２，黄斌斌１，２，李若男１，２，∗，郑　华１，２

１中国科学院生态环境研究中心，城市与区域生态国家重点实验室，北京　１０００８５２中国科学院大学，北京　１０００４９

摘要：气候变化条件下的生态系统响应特征对生态系统服务的提升和生态环境保护具有重要意义。现有气候变化评估多以全球或区域大尺度研究为主，不适合局地小尺度。以２０１５年为基准，根据局地特征改进了综合顺序分类系统（ＣＳＣＳ），模拟了未来不同温室气体排放模式下三江源地区自然植被分布，同时分析了植被覆盖度及土壤保持功能的时空变化。研究结果表明：①在２０８０年的不同排放情境下，三江源地区的降水和气温增长主要发生在生长季（５—９月）和半湿润地区，其中，高排放情境的增幅最大（４２．２１ｍｍ和４．９３℃）；②三江源地区潜在自然植被主要由草地转化为森林，植被覆盖度和土壤保持服务在不同排放情境下均呈现由西北向东南增加的趋势，与水热规律一致，其中高排放情境的增幅最大，低排放最小；③利用三江源地区原生裸地演替特点和多年ＮＰＰ改进的ＣＳＣＳ模型，模拟精度提高了２４％。研究成果为局地小尺度生态系统气候响应特征提供了必要的手段和规律认识。

关键词：修正的综合顺序分类系统（ＣＳＣＳ）；潜在自然植被；植被覆盖度；土壤保持服务；三江源区

ＲｅｓｐｏｎｓｅｃｈａｒａｃｔｅｒｉｓｔｉｃｓｏｆｅｃｏｓｙｓｔｅｍｓａｎｄｓｏｉｌｃｏｎｓｅｒｖａｔｉｏｎｓｅｒｖｉｃｅｓｔｏｆｕｔｕｒｅｃｌｉｍａｔｅｃｈａｎｇｅｉｎｔｈｅＴｈｒｅｅ⁃Ｒｉｖｅｒｓｒｅｇｉｏｎ，Ｃｈｉｎａ

１ＲｅｓｅａｒｃｈＣｅｎｔｅｒｆｏｒＥｃｏ⁃ＥｎｖｉｒｏｎｍｅｎｔａｌＳｃｉｅｎｃｅｓ，ＣｈｉｎｅｓｅＡｃａｄｅｍｙｏｆＳｃｉｅｎｃｅｓ，Ｂｅｉｊｉｎｇ１０００８５，Ｃｈｉｎａ２ＵｎｉｖｅｒｓｉｔｙｏｆＣｈｉｎｅｓｅＡｃａｄｅｍｙｏｆＳｃｉｅｎｃｅｓ，Ｂｅｉｊｉｎｇ１０００４９，Ｃｈｉｎａ

ＦＥＮＧＸｉａｏｙｕ１，２，ＨＵＡＮＧＢｉｎｂｉｎ１，２，ＬＩＲｕｏｎａｎ１，２，∗，ＺＨＥＮＧＨｕａ１，２

Ａｂｓｔｒａｃｔ：Ａｓｓｅｓｓｉｎｇｔｈｅｒｅｓｐｏｎｓｅｃｈａｒａｃｔｅｒｉｓｔｉｃｓｏｆｅｃｏｓｙｓｔｅｍｓｕｎｄｅｒｃｌｉｍａｔｅｃｈａｎｇｅｉｓｃｒｉｔｉｃａｌｔｏｉｍｐｒｏｖｅｅｃｏｓｙｓｔｅｍｓｅｒｖｉｃｅｓａｎｄｐｒｏｔｅｃｔｅｃｏｌｏｇｉｃａｌｅｎｖｉｒｏｎｍｅｎｔｓ．Ｅｘｉｓｔｉｎｇｃｌｉｍａｔｅｃｈａｎｇｅａｓｓｅｓｓｍｅｎｔｍｏｄｅｌｓａｒｅｍｏｓｔｌｙｂａｓｅｄｏｎｌａｒｇｅ⁃ｓｃａｌｅｔｈｉｓｓｔｕｄｙｉｍｐｒｏｖｅｄｔｈｅＣｏｍｐｒｅｈｅｎｓｉｖｅＳｅｑｕｅｎｔｉａｌＣｌａｓｓｉｆｉｃａｔｉｏｎＳｙｓｔｅｍ（ＣＳＣＳ）ｕｓｉｎｇｌｏｃａｌｃｈａｒａｃｔｅｒｉｓｔｉｃｓ，ｓｉｍｕｌａｔｅｄｔｈｅｎａｔｕｒａｌｖｅｇｅｔａｔｉｏｎｄｉｓｔｒｉｂｕｔｉｏｎ，ａｎｄａｎａｌｙｚｅｄｔｈｅｓｐａｔｉａｌａｎｄｔｅｍｐｏｒａｌｃｈａｎｇｅｓｉｎｖｅｇｅｔａｔｉｏｎｃｏｖｅｒａｇｅａｎｄｓｏｉｌｃｏｎｓｅｒｖａｔｉｏｎｓｅｒｖｉｃｅｓｉｎｔｈｅＴｈｒｅｅ⁃Ｒｉｖｅｒｓｒｅｇｉｏｎｕｎｄｅｒｄｉｆｆｅｒｅｎｔｆｕｔｕｒｅｇｒｅｅｎｈｏｕｓｅｇａｓｅｍｉｓｓｉｏｎｐａｔｔｅｒｎｓ．Ｔｈｅｒｅｓｕｌｔｓｉｎｄｉｃａｔｅｄｔｈｅｆｏｌｌｏｗｉｎｇ：①Ｕｎｄｅｒｄｉｆｆｅｒｅｎｔｅｍｉｓｓｉｏｎｓｃｅｎａｒｉｏｓｉｎ２０８０，ｉｎｃｒｅａｓｅｄｐｒｅｃｉｐｉｔａｔｉｏｎａｎｄｔｅｍｐｅｒａｔｕｒｅｍａｉｎｌｙｏｃｃｕｒｒｅｄｉｎｔｈｅｅｍｉｓｓｉｏｎｓｃｅｎａｒｉｏ（４２．２１ｍｍａｎｄ４．９３℃）．②Ｔｈｅｐｏｔｅｎｔｉａｌｎａｔｕｒａｌｖｅｇｅｔａｔｉｏｎｔｙｐｅｗａｓｍａｉｎｌｙｔｒａｎｓｆｏｒｍｅｄｆｒｏｍｇｒａｓｓｌａｎｄｓｏｕｔｈｅａｓｔｕｎｄｅｒｔｈｅｔｈｒｅｅｓｃｅｎａｒｉｏｓ，ｗｈｉｃｈｉｓｃｏｎｓｉｓｔｅｎｔｗｉｔｈｔｈｅｔｒｅｎｄｏｂｓｅｒｖｅｄｆｏｒｗａｔｅｒ⁃ｈｅａｔｃｏｎｄｉｔｉｏｎｓ．Ｔｈｅｈｉｇｈｅｓｔ

基金项目：国家自然科学基金面上项目（４１８７１２１７，４１８７１２１８）；西藏重点研发及转化计划项目；中国科学院战略先导专项（ＸＤＡ２００２０４０２）收稿日期：２０１９⁃０３⁃３０；　　网络出版日期：２０２０⁃０７⁃１１∗通讯作者Ｃｏｒｒｅｓｐｏｎｄｉｎｇａｕｔｈｏｒ．Ｅ⁃ｍａｉｌ：ｒｎｌｉ＠ｒｃｅｅｓ．ａｃ．ｃｎ

（ｇｌｏｂａｌａｎｄｒｅｇｉｏｎａｌ）ｒｅｓｅａｒｃｈ，ａｎｄａｒｅｎｏｔｓｕｉｔａｂｌｅｆｏｒｓｍａｌｌｅｒｇｅｏｇｒａｐｈｉｃａｌｓｃａｌｅｓ（ｌｏｃａｌ）．Ｂａｓｅｄｏｎｔｈｅｄａｔａｆｒｏｍ２０１５，

ｇｒｏｗｉｎｇｓｅａｓｏｎ（ｆｒｏｍＭａｙｔｏＳｅｐｔｅｍｂｅｒ）ａｎｄｉｎｓｅｍｉ⁃ｈｕｍｉｄａｒｅａｓ，ａｎｄｔｈｅｌａｒｇｅｓｔｉｎｃｒｅａｓｅｓｏｃｃｕｒｒｅｄｕｎｄｅｒｔｈｅｈｉｇｈ⁃ｔｏｆｏｒｅｓｔ，ａｎｄｉｎｃｒｅａｓｅｓｉｎｖｅｇｅｔａｔｉｏｎｃｏｖｅｒａｇｅａｎｄｓｏｉｌｃｏｎｓｅｒｖａｔｉｏｎｓｅｒｖｉｃｅｓｆｏｌｌｏｗｅｄａｄｅｓｃｅｎｄｉｎｇｔｒｅｎｄｆｒｏｍｎｏｒｔｈｗｅｓｔｔｏ

ｈｔｔｐ：／／ｗｗｗ．ｅｃｏｌｏｇｉｃａ．ｃｎ

６３５２　生　态　学　报　　　４０卷　

ｉｎｃｒｅａｓｅｓｏｃｃｕｒｒｅｄｕｎｄｅｒｔｈｅｈｉｇｈｅｍｉｓｓｉｏｎｓｃｅｎａｒｉｏａｎｄｔｈｅｌｏｗｅｓｔｉｎｃｒｅａｓｅｓｕｎｄｅｒｔｈｅｌｏｗｅｍｉｓｓｉｏｎｓｃｅｎａｒｉｏ．③Ａｃｃｕｒａｃｙｗａｓｉｍｐｒｏｖｅｄｂｙ２４％ｗｈｅｎｔｈｅＣＳＣＳｗａｓｒｅｖｉｓｅｄｕｓｉｎｇｌｏｃａｌｇｌａｃｉａｌｂａｒｅｌａｎｄｓｕｃｃｅｓｓｉｏｎｃｈａｒａｃｔｅｒｉｓｔｉｃｓａｎｄｍｕｌｔｉ⁃ｙｅａｒＮＰＰ（ＮｅｔＰｒｉｍａｒｙＰｒｏｄｕｃｔｉｏｎ）．Ｔｈｅｓｅｒｅｓｕｌｔｓｐｒｏｖｉｄｅｄｎｅｃｅｓｓａｒｙｍｅｔｈｏｄｓｆｏｒａｓｓｅｓｓｉｎｇｔｈｅｒｅｓｐｏｎｓｅｃｈａｒａｃｔｅｒｉｓｔｉｃｓｏｆｅｃｏｓｙｓｔｅｍｓｕｎｄｅｒｃｌｉｍａｔｅｃｈａｎｇｅａｔｓｍａｌｌｅｒｇｅｏｇｒａｐｈｉｃａｌｓｃａｌｅｓ．

ＫｅｙＷｏｒｄｓ：ｍｏｄｉｆｉｅｄＣｏｍｐｒｅｈｅｎｓｉｖｅＳｅｑｕｅｎｔｉａｌＣｌａｓｓｉｆｉｃａｔｉｏｎＳｙｓｔｅｍ（ＣＳＣＳ）；ｐｏｔｅｎｔｉａｌｎａｔｕｒａｌｖｅｇｅｔａｔｉｏｎ；ｖｅｇｅｔａｔｉｏｎｃｏｖｅｒａｇｅ；ｓｏｉｌｃｏｎｓｅｒｖａｔｉｏｎｓｅｒｖｉｃｅ；Ｔｈｒｅｅ⁃Ｒｉｖｅｒｓｒｅｇｉｏｎ

是气温和降水，是影响陆地生态系统格局和质量的主要因素［２］；同时，陆地生态系统的发育和演替又通过生部分地表植被适应新生境的能力较差，使得其生产力和物候特征都发生了较大的改变［４⁃５］。一些观察结果也表明，自然植被易受到气候变化的影响，并且气候对植被的影响很大程度上取决于气候变化的程度和性质［６⁃７］。虽然未来温室气体的排放引起的全球气候变化存在着不确定性，但分析其对自然植被的影响和由此带来的生态系统服务功能的变化是十分必要的。

气候变化一方面通过降水和气温直接影响陆地生态系统的支撑服务功能，另一方面通过影响土地和其他资源利用强度间接影响陆地生态系统的供给、调节等服务功能［８］。目前，关于气候变化对陆地生态系统影响的研究中，大多基于代表性浓度路径（ＲＣＰｓ）在全球及区域尺度上进行生态系统建模［９⁃１０］，依据局地尺度特征的研究相对较少。

国外对于气候与生态系统及相关服务功能之间的关系的研究相对较早，方法也较成熟，最具代表性的是“霍尔德里奇生命地带（Ｈｏｌｄｒｉｄｇｅｌｉｆｅｚｏｎｅ）”模型，该模型根据气候变量指标（降水、温度、蒸散率）将全球的生统的研究，Ｌｉｕ等［１２］基于“ＨａｄｌｅｙＣｅｎｔｅｒＣｌｉｍａｔｅＭｏｄｅｌ（ＨａｄＣＭ３）”的Ａ２排放情景，发现美国西部干旱地区的Ｓｉｍｕｌａｔｉｏｎ”模型的ＣＯ２高排放情景发现２０５０年全球渔业收入将比预期减少３５％；Ｋｕｎｉｍｉｔｓｕ［１４］利用“Ｒｅｃｕｒｓｉｖｅ⁃ｄｙｎａｍｉｃＲｅｇｉｏｎａｌＣｏｍｐｕｔａｂｌｅＧｅｎｅｒａｌＥｑｕｉｌｉｂｒｉｕｍ（ＣＧＥ）”模型发现气候变化增加了日本西部地区的水稻产量和农民收入，但东部和北部地区却有所减少；Ｓａｍｐｌｅ等［１５］开发了“Ｈｉｇｈ⁃ｌｅｖｅｌＳｃｒｅｅｎｉｎｇ农业灌溉和水力发电量。

气候变化和陆地生态系统之间的相互作用一直被认为是全球变化研究的主要问题之一［１］。气候，特别

物地球化学循环反作用于气候［３］。现有研究成果表明，在过去几十年持续的气温升高和降水量增加趋势下，

物区系划分为３８个生命带类型和１００多个生命地带［１１］。生态系统服务功能对气候变化的响应也有较为系气候变化使森林生物多样性、空气净化和固碳能力提高；Ｌａｍ等［１３］基于“Ｃｌｉｍａｔｅ⁃ｌｉｖｉｎｇＭａｒｉｎｅＲｅｓｏｕｒｃｅｓ

Ｍｅｔｈｏｄｏｌｏｇｙ”来评估与水有关的生态系统服务对未来气候变化的响应，发现未来气候变化降低了苏格兰地区

国内对于气候—植被关系的研究相对起步较晚，但发展较快，其中影响最广的是任继周在１９５６年基于草原分类理论提出的潜在植被分类系统⁃综合顺序分类系统（Ｃｏｍｐｒｅｈｅｎｓｉｖｅｓｅｑｕｅｎｔｉａｌｃｌａｓｓｉｆｉｃａｔｉｏｎｓｙｓｔｅｍ，ＣＳＣＳ），该方法采用定量分类法，体现了草地类型的地带性及内部发生学关系，其理论几经演变、完善，得到了水气压、日照、风速）之间的相关性，发现ＮＤＶＩ值从西向东逐渐增大，与水分分布一致，且与气温、蒸气压呈正模拟了我国的潜在自然植被类型，并分析了植被类型迁移的重心；张明军和周立华［１９］在未来３种不同气候变化情景下预测我国森林生态系统服务价值，结果表明该服务价值呈增加趋势；Ｗｕ等［２０］在ＲＣＰｓ下预测了我国２０１０—２０５０年的草地生态系统生产力的变化趋势，发现温度和降水是影响草地生产力的主要因素，温度的升高更有利于青藏高原和西南喀斯特地区草地生产力的提高。徐雨晴等［２１］基于ＣＥＶＳ模型研究了气候变化背景下中国草地生态系统服务价值的动态变化，结果表明其空间分布状态与降水、气温的分布一致，都呈现出从西南向东北逐渐增加的趋势。

ｈｔｔｐ：／／ｗｗｗ．ｅｃｏｌｏｇｉｃａ．ｃｎ

广泛的应用［１６］。此外，Ｈｕ等［１７］研究了三江源地区１９８２—２０００年生长季节的ＮＤＶＩ与气候因子（均温、降水、相关，与日照时数和风速呈负相关；车彦军等［１８］基于ＲｅｇＣＭ３模型预测的我国２０７１—２１００年的气温和降水，

　１８期　　　冯晓玙　等：三江源区生态系统和土壤保持服务对未来气候变化的响应特征　６３５３

三江源地区是长江、黄河及澜沧江的发源地和重要补给地（分别占长江、黄河、澜沧江总水量的２５％、４９％和１５％）。由于自然条件恶劣，该区域生态系统极其脆弱，由于近百年来气候变化和人类活动干扰的加剧，三江源地区冰川退缩、水土流失加剧、水源涵养能力退减，威胁到了长江、黄河、澜沧江流域的生态安全［２２］。研究三江源地区气候变化对生态系统特征以及重要生态系统服务的影响不仅有利于局地自然生态系统的保护，更有利于三大流域中、下游地区生态系统服务的可持续利用。本文在三江源地区未来不同气候变化情景下，利用修正的综合顺序分类系统（ＣＳＣＳ）模型和多元线性回归建模模拟局地生态系统潜在自然植被空间分布和植被覆盖度的变化，进而探讨气候及植被变化对三江源地区土壤保持服务功能的影响。１　方法与资料１．１　研究区域

三江源地区（３１°３９′—３７°１０′Ｎ、８９°２４′—１０２°２７′Ｅ）位于我国青海省南部，青藏高原腹地，是长江、黄河和澜沧江三大河流的发源地，被称为“中华水塔”。该区域总面积约３６１２００ｋｍ２，平均海拔３５００—４８００ｍ，年平均气温低于－—１℃，生长季平均日照时数为１１８３小时，年平均降水量约为４７０ｍｍ，其中近８０％的降水集中在６月至９月，形成了单一、短暂的植被生长季节，半湿润、半干旱、干旱地区分别约占研究区总面的２９％、５９％、１２％（图１）。

图１　研究区概况Ｆｉｇ．１　Ｓｔｕｄｙａｒｅａ

１．２　研究方法

１．２．１　未来气候变化模式

研究中选择的未来气候变化模式来自于ＩＰＣＣ第五次报告提出的以稳定浓度为特征的新情景。该情景根据大气辐射强度设定，即“代表性浓度路径（ＲｅｐｒｅｓｅｎｔａｔｉｖｅＣｏｎｃｅｎｔｒａｔｉｏｎＰａｔｈｗａｙｓ，ＲＣＰ）”，并首次考虑到各１．２．２　修正的潜在自然植被模型

国应对气候变化政策对未来排放的影响［２３］，本文选取用其中３种排放情景（低、中、高）进行比较分析（表１）。

本文使用综合顺序分类系统（ＣｏｍｐｒｅｈｅｎｓｉｖｅＳｅｑｕｅｎｔｉａｌＣｌａｓｓｉｆｉｃａｔｉｏｎＳｙｓｔｅｍ，ＣＳＣＳ）来预测潜在自然植被类型（图２）。ＣＳＣＳ主要是通过对具有相似降水和温度特性的单元进行分组或聚类而建立的。该系统以活动积温和湿度作为分类指标，包括７个活动积温级和６个湿度级，共组合为４２种潜在自然植被类型。活动积温由低到高对应自然界中北极向赤道递增的经向地带性规律，湿度由干到湿对应自然界中干旱向湿润地区变化的纬向地带性规律。在众多影响未来植被生长因素难以获取的过程中，ＣＳＣＳ模型具有参数易于定义和测量

ｈｔｔｐ：／／ｗｗｗ．ｅｃｏｌｏｇｉｃａ．ｃｎ

６３５４　生　态　学　报　　　４０卷　

的特点，凸显出仅依靠气温、降水就就能够实现高精度模拟的优势。

表１　代表性浓度路径（ＲＣＰｓ）

Ｔａｂｌｅ１　ＲｅｐｒｅｓｅｎｔａｔｉｖｅＣｏｎｃｅｎｔｒａｔｉｏｎＰａｔｈｗａｙｓ（ＲＣＰｓ）

气候情景

Ｃｌｉｍａｔｅｓｃｅｎａｒｉｏ低排放

Ｌｏｗｅｍｉｓｓｉｏｎ（ＲＣＰ２．６）中排放

Ｍｅｄｉｕｍｅｍｉｓｓｉｏｎ（ＲＣＰ４．５）高排放

Ｈｉｇｈｅｍｉｓｓｉｏｎ（ＲＣＰ８．５）

辐射强迫路径

Ｒａｄｉａｔｉｖｅｆｏｒｃｉｎｇｐａｔｈ

从２０１０年到２１００年累计排放温室气体与基础情景相比下降了７０％，到２１００年辐射强迫稳定在２．６Ｗ／ｍ２

辐射强迫不超过４．５Ｗ／ｍ２并在２１００年之后趋于稳定

长时间高温室气体排放，到２１００年辐射强迫稳定在８．５Ｗ／ｍ２

情景描述

Ｓｃｅｎａｒｉｏｄｅｓｃｒｉｐｔｉｏｎ

人类大力实施节能减排措施，通过利用清洁能源和恢复森林来增加对温室气体的吸收，控制其浓度上升

发展与经济相适应的积极减排手段，介于ＲＣＰ２．６和ＲＣＰ８．５之间

人口数量多，技术革新率不高，能源改善缓慢，收入增长慢，长时间高温室气体排放，浓度增加

刘刚等［２４］研究发现我国植被ＮＰＰ空间上基本呈由东南向西北递减的分布趋势，主要是由于植被分布和气候条件决定的，而且青藏高原及南方地区ＮＰＰ与温度呈正相关，ＮＰＰ与温度、降雨量均存在明显相关性。根据三江源地区冰川裸地的地实际演替变化特征以及当地植被的生长能力，并通过查阅大量文献发现三江源地区净初级生产力小于青藏高原地区，而ＣＳＣＳ模型是基于大尺度范围的预测模型，目前在中国的大尺度范围内已经得到了较好的应用［１８，２５⁃２８］。因此，本文采用如下方法修正ＣＳＣＳ模型：

（１）研究区域内冰川、裸地保持不变

有研究表明冰川消退后的原生裸地需经过上百年时间才能适宜生物生长［２９］，原生演替早期阶段生物种

类少，且其生产力及对环境因子的影响也较微弱［３０］。而青藏高原的冰川缩减只有约３０年的时间，区域的土壤、有机质及影响元素还不足以支持生物生存，故在研究中将冰川、裸地保持不变。

（２）根据ＮｅｔＰｒｉｍａｒｙＰｒｏｄｕｃｔｉｏｎ（ＮＰＰ）调整活动积温和湿度指数

顺序分类体系适应于中国整个地区，其参数对于三江源局地尺度来说适用性较低。因此，本文对ＣＳＣＳ模型原始的活动积温和湿度指数分别乘以三江源地区与青藏高原地区２０００—２０１５年多年ＮＰＰ均值的比值（０．７７）对参数进行调整，调整结果如图２。

湿度指数是年均降水量和活动积温的函数，可以反映陆地生态系统植被的自然发生，表达式为：

Ｋ′＝Ｋ×０．７７

（１）（２）

ＮＰＰ是生态系统生长和繁殖的能量值，是生态系统质量的量化表现和植被生长动态的重要参数。综合

(∑θ)′＝∑θ×０．７７

式中，Ｋ为湿度指数，∑θ为活动积温，Ｋ′和(∑θ)′分别为修正后的湿度指数和活动积温。

被类型的分类。

１．２．３　植被覆盖度预估模型

利用２０１５年的自然植被图（图３）对修正后的ＣＳＣＳ模型进行验证，结果表明修正后的模型达到了７７％，

比未修正的模型模拟的准确率提高了２４％。因此，修正后的ＣＳＣＳ模型能够更准确地反映三江源地区不同植

本文根据不同的物候特征，采用多元线性回归模型建立气候与植被覆盖度的关系，用来预测三江源地区的未来植被覆盖度。选取２０００—２０１５年１７个气象站点（图１）的气候、地形因子（均温、年降水、最低气温、最高气温、水气压、日照时数、风速、ＤＥＭ）与相应的植被覆盖度，结合物候（全年、生长季）和气候（整个地区、半湿润地区、干旱和半干旱地区）特征，分析比较不同气候因子组合模式下模型的拟合效果。最终确定了５种气候、地形因子（均温、年降水、最低气温、最高气温、ＤＥＭ），６种物候⁃气候模式的模型用于三江源区域植被覆盖度模拟（表２）。

ｈｔｔｐ：／／ｗｗｗ．ｅｃｏｌｏｇｉｃａ．ｃｎ

　１８期　　　冯晓玙　等：三江源区生态系统和土壤保持服务对未来气候变化的响应特征　６３５５

图２　修正后的综合顺序分类系统（ＣＳＣＳ）

Ｆｉｇ．２　ＭｏｄｉｆｉｅｄＣｏｍｐｒｅｈｅｎｓｉｖｅＳｅｑｕｅｎｔｉａｌＣｌａｓｓｉｆｉｃａｔｉｏｎＳｙｓｔｅｍ（ＣＳＣＳ）

图３　ＣＳＣＳ模型修正前和修正后的２０１５年植被类型图

Ｆｉｇ．３　Ｍａｐｏｆｖｅｇｅｔａｔｉｏｎｔｙｐｅｓｉｎ２０１５ｂｅｆｏｒｅａｎｄａｆｔｅｒＣＳＣＳｍｏｄｅｌｍｏｄｉｆｉｃａｔｉｏｎ

表２　不同模式下三江源区域植被覆盖度的拟合优度

Ｔａｂｌｅ２　ＧｏｏｄｎｅｓｓｏｆｆｉｔｏｆｖｅｇｅｔａｔｉｏｎｃｏｖｅｒａｇｅｕｎｄｅｒｄｉｆｆｅｒｅｎｔｍｏｄｅｌｓｉｎｔｈｅＴｈｒｅｅ⁃Ｒｉｖｅｒｓｒｅｇｉｏｎ

物候

Ｐｈｅｎｏｌｏｇｙ全年Ｗｈｏｌｅｙｅａｒ生长季Ｇｒｏｗｉｎｇｓｅａｓｏｎ

　　不同模式下均为Ｐ＜０．０１，为极显著

整个地区Ａｌｌａｒｅａ０．５７０．５２

半湿润地区Ｓｕｂ⁃ｈｕｍｉｄａｒｅａ

０．６９０．４９

干旱、半干旱地区Ａｒｉｄ，Ｓｅｍｉ⁃ａｒｉｄａｒｅａｓ

０．６９０．４７

在利用２００５—２０１５年的数据建立多元回归模型后，利用２０００—２００５年的数据对模型进行验证，结果表明模型模拟效果较好（图４），可用于三江源地区植被覆盖度的预测。１．２．４　土壤保持服务功能

ＲＵＳＬＥ模型计算土壤保持服务，周来等通过对国内外ＲＵＳＬＥ应用实践进行总结和对比，并分析其科学合理

ｈｔｔｐ：／／ｗｗｗ．ｅｃｏｌｏｇｉｃａ．ｃｎ

土壤保持功能量是潜在土壤侵蚀量和实际土壤侵蚀量之差［３１］。本研究采用周来等［３２］校正和验证好的

６３５６　生　态　学　报　　　４０卷　

图４　三江源区植被覆盖度模型验证

Ｆｉｇ．４　ＶｅｇｅｔａｔｉｏｎｃｏｖｅｒａｇｅｍｏｄｅｌｖｅｒｉｆｉｃａｔｉｏｎｉｎｔｈｅＴｈｒｅｅ⁃Ｒｉｖｅｒｓｒｅｇｉｏｎ

性，明确不同因子单位类型之间的转化系数，使得ＲＵＳＬＥ在中国具体应用过程中更加科学和便捷，该服务功能受降雨、土壤、地形和植被的强烈影响，该服务模型可表示为：

潜在土壤侵蚀：实际土壤侵蚀：土壤保持：

ＳＥａ＝Ｒ×Ｋ×ＬＳ×ＣＳＣ＝Ｒ×Ｋ×ＬＳ×（１－Ｃ）

ＳＥＰ＝Ｒ×Ｋ×ＬＳ

（３）（４）（５）

式中，ＳＥｐ是潜在的侵蚀量（ｔｈｍ－２ａ－１））；ＳＥａ是实际的侵蚀量（ｔｈｍ－２ａ－１）；ＳＣ是土壤保持能力（ｔｈｍ－２ａ－１）；Ｒ的地形因子；Ｃ是植被覆盖因子。１．３　数据来源

研究中涉及到的数据包括研究区域的生态系统分类、植被覆盖度、高程、土壤、气象等，具体参见表３。

表３　数据来源及精度Ｔａｂｌｅ３　Ｄａｔａｓｏｕｒｃｅａｎｄａｃｃｕｒａｃｙ

数据清单Ｄａｔａｌｉｓｔ

精度Ａｃｃｕｒａｃｙ３０ｍ２５０ｍ９０ｍ１ｋｍ１ｋｍ气象站点０．５°

数据来源Ｄａｔａｓｏｕｒｃｅ

中国生态系统评估与生态安全数据库（ｈｔｔｐ：／／ｗｗｗ．ｅｃｏｓｙｓｔｅｍ．ｃｓｄｂ．ｃｎ）

中国生态系统评估与生态安全数据库（ｈｔｔｐ：／／ｗｗｗ．ｅｃｏｓｙｓｔｅｍ．ｃｓｄｂ．ｃｎ）

地理空间数据云（ｈｔｔｐ：／／ｗｗｗ．ｇｓｃｌｏｕｄ．ｃｎ／）寒区旱区科学数据中心寒区旱区科学数据中心北京师范大学

ＩＰＣＣ的ＳＲＥＳ研究团队利用ＧＣＭ驱动的模拟数据，ＣＣＡＦＳ（ＣｌｉｍａｔｅＣｈａｎｇｅ，ＡｇｒｉｃｕｌｔｕｒｅａｎｄＦｏｏｄＳｅｃｕｒｉｔｙ）组织发布（ｈｔｔｐ：／／ｗｗｗ．ｃｃａｆｓ⁃ｃｌｉｍａｔｅ．ｏｒｇ／）

是雨量的侵蚀因子（ＭＪｍｍｈｍ－２ｈ－１ａ－１）；Ｋ是土壤可蚀性因子（ｔｈｍ２ｈｈｍ－２ＭＪ－１ｍｍ－１）；ＬＳ是表示坡长影响

２０１５年生态系统分类图

Ｅｃｏｓｙｓｔｅｍｃｌａｓｓｉｆｉｃａｔｉｏｎｉｎ２０１５２０１５年植被覆盖度

Ｖｅｇｅｔａｔｉｏｎｃｏｖｅｒａｇｅｉｎ２０１５数字高程模型

Ｄｉｇｉｔａｌｅｌｅｖａｔｉｏｎｍｏｄｅｌ土壤有机质含量Ｓｏｉｌｏｒｇａｎｉｃｍａｔｔｅｒ土壤机械组成Ｓｏｉｌｔｅｘｔｕｒｅ

２０００—２０１５年降雨、气温

Ｐｒｅｃｉｐｉｔａｔｉｏｎａｎｄａｉｒｔｅｍｐｅｒａｔｕｒｅｆｒｏｍ２０００—２０１５２０８０年降雨、气温

Ｐｒｅｃｉｐｉｔａｔｉｏｎａｎｄａｉｒｔｅｍｐｅｒａｔｕｒｅｉｎ２０８０

２　结果与分析

２．１　三江源地区未来气候变化

三江源地区气候模拟结果表明，２０８０年的降水和气温在不同排放模式下均高于２０１５年，增加幅度均表

ｈｔｔｐ：／／ｗｗｗ．ｅｃｏｌｏｇｉｃａ．ｃｎ

　１８期　　　冯晓玙　等：三江源区生态系统和土壤保持服务对未来气候变化的响应特征　６３５７

现为：高排放情景＞中排放情景＞低排放情景。其中，高、中、低排放情景的全年降水增幅（变化率）分别为４２．２１ｍｍ（１０．５７％）、３１．４１ｍｍ（７．８７％）、２４．２３ｍｍ（６．０７％）；全年气温增幅（变化率）分别为４．９３℃（１１５．７３％）、２．７０℃（６３．３８％）、１．８８℃（４４．０１％）。

从时间上看，降水和气温的增加主要发生在生长季节；从空间上看，半湿润地区的降水增加占据主导地

位，干旱地区由于温度的上升导致降水量有所下降，半湿润地区气温增量较大，高排放气候变化情景下年均温度达到了０℃以上。

２．２　三江源地区未来生态系统及植被覆盖状况

依据２０１５年的生态系统现状，三江源地区主要的生态系统类型为森林、灌丛、草地、湿地、农田、人工地、荒漠和裸地８大类别。其中，草地面积最大，为２４６１０４．１ｋｍ２，占总面积的６９．００％，由于降水、气温等气候因素和海拔等地势因素的空间异质性使得草地生态系统类型多样化；其次为荒漠和裸地，占１７．３０％；森林和灌丛占比为４．４３％，主要分布于适宜植被生长的气温高、降水多、海拔相对较低的东南地区；农田和人工地面积占比最小，分别为５８４．９４ｋｍ２（０．１６％）和５２１．３８ｋｍ２（０．１５％），主要分布在三江源东部。

在不同气候变化情景下，２０８０年生态系统类型的整体分布格局在空间上与２０１５年基本保持一致（图５）；

森林和灌丛的面积沿着东南向西北的方向在逐渐增加；少量草地面积转化为森林和灌丛。其中，森林和灌丛的变化率为高排放情景（２９３．４６％）＞中排放情景（１２９．９６％）＞低排放情景（６１．０８％）；相应地，草地面积随着排放浓度的增加而逐渐减小，变化率为高排放情景（－１８．９９％）＞中排放情景（－８．３３％）＞低排放情景（－３．９１％）。

图５　三江源区域２０１５年生态系统类型分布及２０８０年不同气候模式下生态系统类型预测

Ｆｉｇ．５　Ｄｉｓｔｒｉｂｕｔｉｏｎｏｆｅｃｏｓｙｓｔｅｍｔｙｐｅｓｉｎ２０１５ａｎｄｐｒｅｄｉｃｔｉｏｎｏｆｅｃｏｓｙｓｔｅｍｔｙｐｅｓｕｎｄｅｒｄｉｆｆｅｒｅｎｔｃｌｉｍａｔｅｍｏｄｅｌｓｉｎ２０８０ｉｎｔｈｅＴｈｒｅｅ⁃Ｒｉｖｅｒｓｒｅｇｉｏｎ

对于植被覆盖度来说，三江源地区２０１５年和２０８０年３种气候变化情景下的植被覆盖度空间变化趋势一致，整体上从西北到东南方向植被覆盖度逐渐呈现出增加趋势（图６）。平均植被覆盖度分布为高排放情景（８２．０７％）＞中排放情景（６２．６６％）＞低排放情景（５２．２８％）＞２０１５年（４４．５５％）；相比２０１５年，植被覆盖度增长率为高排放情景（８４．２０％）＞中排放情景（４０．６４％）＞低排放情景（１７．３４％）。

由于三江源地区位于青藏高原高海拔的气候敏感区，气温较低、降水较少，生态系统类型和植被生长对水

ｈｔｔｐ：／／ｗｗｗ．ｅｃｏｌｏｇｉｃａ．ｃｎ

６３５８　生　态　学　报　　　４０卷　

热条件极为敏感。本文中采用的ＩＰＣＣ第五次报告中的３种排放情景，其气温和降水均为增加趋势，所以３种排放模式下模拟的生态系统类型以及植被覆盖度均朝着变好的趋势发展。

图６　三江源区域２０１５年植被覆盖度及２０８０年不同排放情景下植被覆盖度预测

Ｆｉｇ．６　Ｖｅｇｅｔａｔｉｏｎｃｏｖｅｒａｇｅｉｎ２０１５ａｎｄｐｒｅｄｉｃｔｉｏｎｏｆｖｅｇｅｔａｔｉｏｎｃｏｖｅｒａｇｅｕｎｄｅｒｄｉｆｆｅｒｅｎｔｅｍｉｓｓｉｏｎｓｃｅｎａｒｉｏｓｉｎ２０８０ｉｎｔｈｅＴｈｒｅｅ⁃Ｒｉｖｅｒｓｒｅｇｉｏｎ

２．３　三江源土壤保持服务

土壤保持服务功能同时受到植被、气候和土壤等因素的共同影响，与生态系统格局和植被覆盖度紧密相关。土壤保持服务在三江源地区的东南方向功能较高，此处温度较高，降水较多，海拔较低，适宜植被生长，而西北方向冰川裸地较多，土壤保持功能相对较低。土壤流失方程能够综合考虑植被、气候和土壤等影响因素，从模拟结果上来看也能够有效的反应研究区的实际情况。

三江源地区２０１５年和２０８０年３种气候变化情景下的土壤保持服务与生态系统植被状况的变化趋势密切相关。土壤保持能力整体上从西北向东南逐渐增强。其中，高排放情景＞中排放情景＞低排放情景＞２０１５年，对应的土壤保持能力（增长率）为８．１４亿吨（２８．１９％）＞７．６９亿吨（２１．１０％）＞７．０９亿吨（１１．６５％）＞６．３５亿吨。

总体而言，三江源地区２０１５年和２０８０年３种气候变化情景下的潜在自然植被分布、植被覆盖度和土壤保持功能关系十分密切。２０１５年的森林生态系统面积最小，主要集中在三江源的东南部；２０８０年３种排放模式的森林生态系统面积呈现由低排放模式至高排放模式逐渐增加趋势，森林分布由三江源东南部向西北扩展，但止步于４００ｍｍ降水线。相应地，植被覆盖度也随着森林生态系统取代草地生态系统而由东南向西北逐渐增加；高排放模式下植被覆盖度达到最高的８２．０７％，是低排放模式的１．５７倍。３种排放模式下，土壤保持服务南高北低、东高西低的空间格局没有变化，但差异有所增加，南北变幅大于东西。３　讨论与结论３．１　讨论

气候变化对生态系统格局、植被覆盖度以及生态系统服务功能产生重要的影响，研究气候⁃植被之间的相互关系有助于应对气候变化给生态系统带来的负面影响，进而有效的保护生态环境。为了使ＣＳＣＳ模型能够

ｈｔｔｐ：／／ｗｗｗ．ｅｃｏｌｏｇｉｃａ．ｃｎ

　１８期　　　冯晓玙　等：三江源区生态系统和土壤保持服务对未来气候变化的响应特征　６３５９

图７　三江源区域２０１５年土壤保持服务分布及２０８０年不同排放情景下土壤保持服务预测

Ｆｉｇ．７　Ｓｏｉｌｃｏｎｓｅｒｖａｔｉｏｎｓｅｒｖｉｃｅｄｉｓｔｒｉｂｕｔｉｏｎｉｎ２０１５ａｎｄｐｒｅｄｉｃｔｉｏｎｏｆｓｏｉｌｃｏｎｓｅｒｖａｔｉｏｎｓｅｒｖｉｃｅｕｎｄｅｒｄｉｆｆｅｒｅｎｔｅｍｉｓｓｉｏｎｓｃｅｎａｒｉｏｓｉｎ２０８０ｉｎｔｈｅＴｈｒｅｅ⁃Ｒｉｖｅｒｓｒｅｇｉｏｎ

适应于三江源地区，本研究利用青藏高原和三江源地区的多年ＮＰＰ以及冰川裸地的演替特点对模型进行修正，修正后的模型比未修正模型模拟精度提高了２４％。但由于影响活动积温和湿度指数的因素较多，本研究仅采用ＮＰＰ对模型进行修正，使得模型精度提升受限，还有待在后续研究中进一步完善。此外，本研究利用线性回归模型预估未来植被覆盖度，基于三江源地区的物候及气候特点建立植被覆盖度与气候因子的多元线性回归模型，虽然模型模拟效果较好，但由于在未来气候变化模式下，物候及气候均存在不确定性，也会给模拟结果带来一定程度的误差。

本文中，在２０８０年不同气候模式下，三江源地区的降水、气温以及森林覆盖面积均呈现上升趋势，同时草地面积减少，这与在青藏高原地区的其他相关研究结果基本一致。如，Ｌｉａｎｇ等［３３］利用ＣＳＣＳ模型对青藏高原地区潜在自然植被的研究结果表明为在未来气候变化下（２０００—２０８０年）年平均降水量（ＭｅａｎＡｎｎｕａｌＰｒｅｃｉｐｉｔａｔｉｏｎ，ＭＡＰ）和年平均温度（ｍｅａｎａｎｎｕａｌｔｅｍｐｅｒａｔｕｒｅ，ＭＡＴ）均呈现增加趋势，同时，森林面积将显著增结果也表明，三江源地区植被覆盖度与气温和降水均成正相关关系，其中对气温的相关性大于降水。在土壤保持研究方面，康惠惠等［３４］的研究表明土地植被类型及植被覆盖度对土壤保持功能有重要影响，三江源的生态退化与恢复过程与土壤保持功能变化联系紧密。综上，本研究中选取的研究方法及主要环境因子是适用于三江源地区的，对于模拟结果中生态系统、植被状况和土壤保持服务功能所表现出的差异，主要是由于情景模式以及空间分辨率的不同所致。

研究中还对ＣＳＣＳ模型中的冰川、裸地生态系统的评估进行了修正，其主要依据是生态演替理论中的原生演替过程［３５］。基于该演替过程，认为到２０８０年，三江源地区冰川积雪和冻土裸土面积虽然会有减少，但其从景观生态学边际效应角度来看，随着冰川的消退，冰川、裸地和草地生态系统间会产生过渡区，即生态过渡

ｈｔｔｐ：／／ｗｗｗ．ｅｃｏｌｏｇｉｃａ．ｃｎ

土壤理化性质还不足以支持植被的生长。所以在本研究中，冰川、裸地生态系统的空间分布不随时间改变。加，草地面积则呈一定的下降趋势。此外，Ｈｕ等［１７］对三江源地区的植被覆盖度与气候因子之间关系的研究

６３６０　生　态　学　报　　　４０卷　

带。由于生态交错带对于气候变化的反应较为敏感，气候变化下生态过渡带会发生生物群区的过渡［３６⁃３７］，所以本研究中冰川、裸地生态系统不变的假设会给模型模拟结果带来不确定性。

植被与作物管理因子（Ｃ）反映了土壤被不同的土壤覆盖类型保护的程度，这是ＵＳＬＥ模型中唯一的组分因子，它受自然因素（气温和降水）以及人类活动（如开垦，建筑，采矿，伐木，修路等）或其他事件（如林火，洪水）的影响很大。由于三江源地区人类干扰较少，本研究目前仅考虑自然因子进行预测植被因子，尚未考虑人类活动因子，这也是本研究后续需要深入探讨的内容。在影响植被生长的因子中，气象是最主要的影响因素，因此在分类体系中，将气象作为最高级的分类指标［３８］。本研究在自然植被模拟过程中，除高程因素外，仅考虑了降水和温度两个自然因子。由于气候模型的制约，３种排放模式下三江源区域的水热条件均呈现增加趋势，但从全国范围来看，水热条件变化并不相同。在３种排放模式下（高、中、低），潜在自然植被的模拟结果表明，内蒙古的北部以及与宁夏、甘肃交错地带和新疆的北部地区的气候干燥，湿度下降，自然植被均呈现出逐渐退化演替的状况；在新疆的塔克拉玛干沙漠地带，只有在中排放情景下生态系统是进化演替的；而青藏高原地区，由于未来气候变化使其增温增湿，利于植被生长，生态系统整体呈现进化演替［２１］。此外，三江源区域的植被覆盖度还受到其他气象因子［１７］以及土壤［３９］等因素的影响，在后续研究中还应该沿着“未来气候变化⁃生态系统结构、质量、过程⁃生态系统服务”框架进行深入探讨，提高模型模拟精度。３．２　结论

本文利用综合利用顺序（ＣＳＣＳ）模型、多元线性回归模型以及土壤流失方程模型，模拟了未来不同排放模式对应的气候情景下潜在自然植被、植被覆盖度和土壤保持服务的时空动态变化，提出了适应于三江源局地尺度修正模型，结论如下：

情景＞低排放情景，增加主要表现在生长季节（５—９月）和半湿润地区；

（１）与２０１５年相比，２０８０年３种气候模式下的气温和降水均呈现增加趋势，幅度为高排放情景＞中排放（２）与２０１５年相比，２０８０年３种气候模式下的植被覆盖度和土壤保持服务均增加，幅度为高排放情景＞

中排放情景＞低排放情景＞２０１５年，空间上呈现出自西北向东南方向逐渐增强的趋势；从生态系统变化来看，森林生态系统面积增加，草地生态系统面积减小；

高了２４％，为今后研究局地尺度的潜在自然植被类型和生态系统服务评估提供了修正方法。

参考文献（Ｒｅｆｅｒｅｎｃｅｓ）：

［１］　ＷａｌｋｅｒＢ，ＳｔｅｆｆｅｎＷＬ．ＡｎＯｖｅｒｖｉｅｗｏｆｔｈｅｉｍｐｌｉｃａｔｉｏｎｓｏｆｇｌｏｂａｌｃｈａｎｇｅｆｏｒｎａｔｕｒａｌａｎｄｍａｎａｇｅｄｔｅｒｒｅｓｔｒｉａｌｅｃｏｓｙｓｔｅｍｓ．ＣｏｎｓｅｒｖａｔｉｏｎＥｃｏｌｏｇｙ，［２］　黄银晓．全球气候变化对陆地生态系统的影响．环境科学，１９９３，１４（６）：７９⁃８１．

８（１２）：１１６９⁃１１８２．１９９７，１（２）：２．

（３）利用原生裸地演替特点和多年ＮＰＰ率定了的ＣＳＣＳ模型参数，使得其在三江源区域的模拟精确度提

［３］　ＷａｎｇＧＬ，ＥｌｔａｈｉｒＥＡＢ．ＩｍｐａｃｔｏｆＣＯ２ｃｏｎｃｅｎｔｒａｔｉｏｎｃｈａｎｇｅｓｏｎｔｈｅｂｉｏｓｐｈｅｒｅ⁃ａｔｍｏｓｐｈｅｒｅｓｙｓｔｅｍｏｆＷｅｓｔＡｆｒｉｃａ．ＧｌｏｂａｌＣｈａｎｇｅＢｉｏｌｏｇｙ，２００２，［４］　ＧａｎｇＣＣ，ＺｈｏｕＷ，ＬｉＪＬ，ＣｈｅｎＹＺ，ＭｕＳＪ，ＲｅｎＪＺ，ＣｈｅｎＪＭ，ＧｒｏｉｓｍａｎＰＹ．Ａｓｓｅｓｓｉｎｇｔｈｅｓｐａｔｉｏｔｅｍｐｏｒａｌｖａｒｉａｔｉｏｎｉｎｄｉｓｔｒｉｂｕｔｉｏｎ，ｅｘｔｅｎｔ［５］　ＨｏｒｉｏｎＳ，ＣｏｒｎｅｔＹ，ＥｒｐｉｃｕｍＭ，ＴｙｃｈｏｎＢ．Ｓｔｕｄｙｉｎｇｉｎｔｅｒａｃｔｉｏｎｓｂｅｔｗｅｅｎｃｌｉｍａｔｅｖａｒｉａｂｉｌｉｔｙａｎｄｖｅｇｅｔａｔｉｏｎｄｙｎａｍｉｃｕｓｉｎｇａｐｈｅｎｏｌｏｇｙｂａｓｅｄ［６］　ＢｉｅｒｍａｎｎＦ．‘Ｅａｒｔｈｓｙｓｔｅｍｇｏｖｅｒｎａｎｃｅ’ａｓａｃｒｏｓｓｃｕｔｔｉｎｇｔｈｅｍｅｏｆｇｌｏｂａｌｃｈａｎｇｅｒｅｓｅａｒｃｈ．ＧｌｏｂａｌＥｎｖｉｒｏｎｍｅｎｔａｌＣｈａｎｇｅ，２００７，１７（３／４）：［７］　ＦｉｅｌｄＣＢ，ＬｏｂｅｌｌＤＢ，ＰｅｔｅｒｓＨＡ，ＣｈｉａｒｉｅｌｌｏＮＲ．Ｆｅｅｄｂａｃｋｓｏｆｔｅｒｒｅｓｔｒｉａｌｅｃｏｓｙｓｔｅｍｓｔｏｃｌｉｍａｔｅｃｈａｎｇｅ．ＡｎｎｕａｌＲｅｖｉｅｗｏｆＥｎｖｉｒｏｎｍｅｎｔａｎｄ［８］　ＲｕｎｔｉｎｇＲＫ，ＢｒｙａｎＢＡ，ＤｅｅＬＥ，ＭａｓｅｙｋＦＪＦ，ＭａｎｄｌｅＬ，ＨａｍｅｌＰ，ＷｉｌｓｏｎＫＡ，ＹｅｔｋａＫ，ＰｏｓｓｉｎｇｈａｍＨＰ，ＲｈｏｄｅｓＪＲ．Ｉｎｃｏｒｐｏｒａｔｉｎｇ［９］　ＡｎａｖＡ，ＭａｒｉｏｔｔｉＡ．ＳｅｎｓｉｔｉｖｉｔｙｏｆｎａｔｕｒａｌｖｅｇｅｔａｔｉｏｎｔｏｃｌｉｍａｔｅｃｈａｎｇｅｉｎｔｈｅＥｕｒｏ⁃Ｍｅｄｉｔｅｒｒａｎｅａｎａｒｅａ．ＣｌｉｍａｔｅＲｅｓｅａｒｃｈ，２０１１，４６（３）：２７７⁃２９２．

ｃｌｉｍａｔｅｃｈａｎｇｅｉｎｔｏｅｃｏｓｙｓｔｅｍｓｅｒｖｉｃｅａｓｓｅｓｓｍｅｎｔｓａｎｄｄｅｃｉｓｉｏｎｓ：ａｒｅｖｉｅｗ．ＧｌｏｂａｌＣｈａｎｇｅＢｉｏｌｏｇｙ，２０１７，２３（１）：２８⁃４１．Ｒｅｓｏｕｒｃｅｓ，２００７，３２：１⁃２９．３２６⁃３３７．

ａｐｐｒｏａｃｈ．ＩｎｔｅｒｎａｔｉｏｎａｌＪｏｕｒｎａｌｏｆＡｐｐｌｉｅｄＥａｒｔｈＯｂｓｅｒｖａｔｉｏｎａｎｄＧｅｏｉｎｆｏｒｍａｔｉｏｎ，２０１３，２０：２０⁃３２．

ａｎｄＮＰＰｏｆｔｅｒｒｅｓｔｒｉａｌｅｃｏｓｙｓｔｅｍｓｉｎｒｅｓｐｏｎｓｅｔｏｃｌｉｍａｔｅｃｈａｎｇｅｆｒｏｍ１９１１ｔｏ２０００．ＰＬｏＳＯｎｅ，２０１３，８（１１）：ｅ８０３９４．

［１０］　ＨｉｃｋｌｅｒＴ，ＶｏｈｌａｎｄＫ，ＦｅｅｈａｎＪ，ＭｉｌｌｅｒＰＡ，ＳｍｉｔｈＢ，ＣｏｓｔａＬ，ＧｉｅｓｅｃｋｅＴ，ＦｒｏｎｚｅｋＳ，ＣａｒｔｅｒＴＲ，ＣｒａｍｅｒＷ，ＫüｈｎＩ，ＳｙｋｅｓＭＴ．Ｐｒｏｊｅｃｔｉｎｇ

ｈｔｔｐ：／／ｗｗｗ．ｅｃｏｌｏｇｉｃａ．ｃｎ

　１８期　　　冯晓玙　等：三江源区生态系统和土壤保持服务对未来气候变化的响应特征　６３６１

ｔｈｅｆｕｔｕｒｅｄｉｓｔｒｉｂｕｔｉｏｎｏｆＥｕｒｏｐｅａｎｐｏｔｅｎｔｉａｌｎａｔｕｒａｌｖｅｇｅｔａｔｉｏｎｚｏｎｅｓｗｉｔｈａｇｅｎｅｒａｌｉｚｅｄ，ｔｒｅｅｓｐｅｃｉｅｓ⁃ｂａｓｅｄｄｙｎａｍｉｃｖｅｇｅｔａｔｉｏｎｍｏｄｅｌ．Ｇｌｏｂａｌ

［１１］　ＨｏｌｄｒｉｄｇｅＬＲ．Ｄｅｔｅｒｍｉｎａｔｉｏｎｏｆｗｏｒｌｄｐｌａｎｔｆｏｒｍａｔｉｏｎｓｆｒｏｍｓｉｍｐｌｅｃｌｉｍａｔｉｃｄａｔａ．Ｓｃｉｅｎｃｅ，１９４７，１０５（２７２７）：３６７⁃３６８．

ＥｃｏｌｏｇｙａｎｄＭａｎａｇｅｍｅｎｔ，２０１３，２９４：１２０⁃１３５．Ｒｅｐｏｒｔｓ，２０１６，６：３２６０７．

ＥｃｏｌｏｇｙａｎｄＢｉｏｇｅｏｇｒａｐｈｙ，２０１２，２１（１）：５０⁃６３．

［１２］　ＬｉｕＹＱ，ＧｏｏｄｒｉｃｋＳＬ，ＳｔａｎｔｕｒｆＪＡ．ＦｕｔｕｒｅＵ．Ｓ．ｗｉｌｄｆｉｒｅｐｏｔｅｎｔｉａｌｔｒｅｎｄｓｐｒｏｊｅｃｔｅｄｕｓｉｎｇａｄｙｎａｍｉｃａｌｌｙｄｏｗｎｓｃａｌｅｄｃｌｉｍａｔｅｃｈａｎｇｅｓｃｅｎａｒｉｏ．Ｆｏｒｅｓｔ［１３］　ＬａｍＶＷＹ，ＣｈｅｕｎｇＷＷＬ，ＲｅｙｇｏｎｄｅａｕＧ，ＳｕｍａｉｌａＵＲ．Ｐｒｏｊｅｃｔｅｄｃｈａｎｇｅｉｎｇｌｏｂａｌｆｉｓｈｅｒｉｅｓｒｅｖｅｎｕｅｓｕｎｄｅｒｃｌｉｍａｔｅｃｈａｎｇｅ．Ｓｃｉｅｎｔｉｆｉｃ［１４］　ＫｕｎｉｍｉｔｓｕＹ．Ｒｅｇｉｏｎａｌｉｍｐａｃｔｓｏｆｌｏｎｇ⁃ｔｅｒｍｃｌｉｍａｔｅｃｈａｎｇｅｏｎｒｉｃｅｐｒｏｄｕｃｔｉｏｎａｎｄａｇｒｉｃｕｌｔｕｒａｌｉｎｃｏｍｅ：ｅｖｉｄｅｎｃｅｆｒｏｍｃｏｍｐｕｔａｂｌｅｇｅｎｅｒａｌ［１５］　ＳａｍｐｌｅＪＥ，ＢａｂｅｒＩ，ＢａｄｇｅｒＲ．Ａｓｐａｔｉａｌｌｙｄｉｓｔｒｉｂｕｔｅｄｒｉｓｋｓｃｒｅｅｎｉｎｇｔｏｏｌｔｏａｓｓｅｓｓｃｌｉｍａｔｅａｎｄｌａｎｄｕｓｅｃｈａｎｇｅｉｍｐａｃｔｓｏｎｗａｔｅｒ⁃ｒｅｌａｔｅｄｅｃｏｓｙｓｔｅｍ［１６］　任继周，胡自治，牟新待．我国草原类型第一级分类的生物气候指标．甘肃农业大学学报，１９６５，（２）：４８⁃６４．

ｒｉｖｅｒ⁃ｓｏｕｒｃｅｒｅｇｉｏｎ．ＩｎｔｅｒｎａｔｉｏｎａｌＪｏｕｒｎａｌｏｆＡｐｐｌｉｅｄＥａｒｔｈＯｂｓｅｒｖａｔｉｏｎａｎｄＧｅｏｉｎｆｏｒｍａｔｉｏｎ，２０１１，１３（１）：２４⁃３３．ｓｅｒｖｉｃｅｓ．ＥｎｖｉｒｏｎｍｅｎｔａｌＭｏｄｅｌｌｉｎｇ＆Ｓｏｆｔｗａｒｅ，２０１６，８３：１２⁃２６．

ｅｑｕｉｌｉｂｒｉｕｍａｎａｌｙｓｉｓ．ＪｏｕｒｎａｌｏｆＪａｐａｎＳｏｃｉｅｔｙｏｆＣｉｖｉｌＥｎｇｉｎｅｅｒｓ，２０１４，７０（５）：Ｉ＿１３⁃Ｉ＿１９．

［１７］　ＨｕＭＱ，ＭａｏＦ，ＳｕｎＨ，ＨｏｕＹＹ．Ｓｔｕｄｙｏｆｎｏｒｍａｌｉｚｅｄｄｉｆｆｅｒｅｎｃｅｖｅｇｅｔａｔｉｏｎｉｎｄｅｘｖａｒｉａｔｉｏｎａｎｄｉｔｓｃｏｒｒｅｌａｔｉｏｎｗｉｔｈｃｌｉｍａｔｅｆａｃｔｏｒｓｉｎｔｈｅｔｈｒｅｅ⁃［１８］　车彦军，赵军，师银芳，王大为．基于ＣＳＣＳ和ＲｅｇＣＭ３模型的２１世纪末中国潜在植被．生态学杂志，２０１４，３３（２）：４４７⁃４５４．［１９］　张明军，周立华．气候变化对中国森林生态系统服务价值的影响．干旱区资源与环境，２００４，１８（２）：４０⁃４３．

ｉｎＣｈｉｎａ．ＡｄｖａｎｃｅｓｉｎＭｅｔｅｏｒｏｌｏｇｙ，２０１３，２０１３：８１２７２３．

［２０］　ＷｕＦ，ＤｅｎｇＸＺ，ＹｉｎＦ，ＹｕａｎＹＷ．Ｐｒｏｊｅｃｔｅｄｃｈａｎｇｅｓｏｆｇｒａｓｓｌａｎｄｐｒｏｄｕｃｔｉｖｉｔｙａｌｏｎｇｔｈｅｒｅｐｒｅｓｅｎｔａｔｉｖｅｃｏｎｃｅｎｔｒａｔｉｏｎｐａｔｈｗａｙｓｄｕｒｉｎｇ２０１０⁃２０５０［２１］　徐雨晴，於琍，周波涛，石英，徐影．气候变化背景下未来中国草地生态系统服务价值时空动态格局．生态环境学报，２０１７，２６（１０）：［２２］　蒋冲，高艳妮，李芬，王德旺，张林波，李岱青．１９５６—２０１０年三江源区水土流失状况演变．环境科学研究，２０１７，３０（１）：２０⁃２９．

１６４９⁃１６５８．

［２３］　ＭｏｓｓＲＨ，ＥｄｍｏｎｄｓＪＡ，ＨｉｂｂａｒｄＫＡ，ＭａｎｎｉｎｇＭＲ，ＲｏｓｅＳＫ，ｖａｎＶｕｕｒｅｎＤＰ，ＣａｒｔｅｒＴＲ，ＥｍｏｒｉＳ，ＫａｉｎｕｍａＭ，ＫｒａｍＴ，ＭｅｅｈｌＧＡ，

ＭｉｔｃｈｅｌｌＪＦＢ，ＮａｋｉｃｅｎｏｖｉｃＮ，ＲｉａｈｉＫ，ＳｍｉｔｈＳＪ，ＳｔｏｕｆｆｅｒＲＪ，ＴｈｏｍｓｏｎＡＭ，ＷｅｙａｎｔＪＰ，ＷｉｌｂａｎｋｓＴＪ．Ｔｈｅｎｅｘｔｇｅｎｅｒａｔｉｏｎｏｆｓｃｅｎａｒｉｏｓｆｏｒｃｌｉｍａｔｅｃｈａｎｇｅｒｅｓｅａｒｃｈａｎｄａｓｓｅｓｓｍｅｎｔ．Ｎａｔｕｒｅ，２０１０，４６３（７２８２）：７４７⁃７５６．

［２４］　刘刚，孙睿，肖志强，崔天翔．２００１—２０１４年中国植被净初级生产力时空变化及其与气象因素的关系．生态学报，２０１７，３７（１５）：［２５］　张颖，陈怀艳，李建龙．三江源生态系统近１０年净初级生产力估测．天津农业科学，２０１４，２０（１０）：２５⁃２８．

４９３６⁃４９４５．

［２６］　孙庆龄，李宝林，李飞，张志军，丁玲玲，张涛，许丽丽．三江源植被净初级生产力估算研究进展．地理学报，２０１６，７１（９）：１５９６⁃１６１２．［２７］　孙云晓，王思远，常清，薛铸鑫，尹航．青藏高原近３０年植被净初级生产力时空演变研究．广东农业科学，２０１４，４１（１３）：１６０⁃１６６．

２０１３，６８（９）：１１９７⁃１２１１．４２７⁃４４８．

［２８］　张镱锂，祁威，周才平，丁明军，刘林山，高俊刚，摆万奇，王兆锋，郑度．青藏高原高寒草地净初级生产力（ＮＰＰ）时空分异．地理学报，［２９］　ＣｒｏｃｋｅｒＲＬ，ＭａｊｏｒＪ．Ｓｏｉｌｄｅｖｅｌｏｐｍｅｎｔｉｎｒｅｌａｔｉｏｎｔｏｖｅｇｅｔａｔｉｏｎａｎｄｓｕｒｆａｃｅａｇｅａｔｇｌａｃｉｅｒｂａｙ，Ａｌａｓｋａ．ＪｏｕｒｎａｌｏｆＥｃｏｌｏｇｙ，１９５５，４３（２）：［３０］　ＷｈｉｔｔａｋｅｒＲＪ，ＢｕｓｈＭＢ，ＲｉｃｈａｒｄｓＫ．ＰｌａｎｔｒｅｃｏｌｏｎｉｚａｔｉｏｎａｎｄｖｅｇｅｔａｔｉｏｎｓｕｃｃｅｓｓｉｏｎｏｎｔｈｅＫｒａｋａｔａｕＩｓｌａｎｄｓ，Ｉｎｄｏｎｅｓｉａ．ＥｃｏｌｏｇｉｃａｌＭｏｎｏｇｒａｐｈｓ，［３１］　ＲａｏＥＭ，ＯｕｙａｎｇＺＹ，ＹｕＸＸ，ＸｉａｏＹ．ＳｐａｔｉａｌｐａｔｔｅｒｎｓａｎｄｉｍｐａｃｔｓｏｆｓｏｉｌｃｏｎｓｅｒｖａｔｉｏｎｓｅｒｖｉｃｅｉｎＣｈｉｎａ．Ｇｅｏｍｏｒｐｈｏｌｏｇｙ，２０１４，２０７：６４⁃７０．［３２］　周来，李艳洁，孙玉军．修正的通用土壤流失方程中各因子单位的确定．水土保持通报，２０１８，３８（１）：１６９⁃１７４．

［３３］　ＬｉａｎｇＴＧ，ＦｅｎｇＱＳ，ＹｕＨ，ＨｕａｎｇＸＤ，ＬｉｎＨＬ，ＡｎＳＺ，ＲｅｎＪＺ．ＤｙｎａｍｉｃｓｏｆｎａｔｕｒａｌｖｅｇｅｔａｔｉｏｎｏｎｔｈｅＴｉｂｅｔａｎＰｌａｔｅａｕｆｒｏｍｐａｓｔｔｏｆｕｔｕｒｅ［３４］　康惠惠，潘韬，盖艾鸿，刘玉洁．生态退化与恢复对三江源区土壤保持功能的影响．水土保持通报，２０１７，３７（３）：７⁃１４．［３５］　许中旗，李文华，鲍维楷，许晴．植被原生演替研究进展．生态学报，２００５，２５（１２）：３３８３⁃３３８９．［３６］　ＯｄｕｍＥＰ．ＦｕｎｄａｍｅｎｔａｌｓｏｆＥｃｏｌｏｇｙ．３ｒｄｅｄ．Ｐｈｉｌａｄｅｌｐｈｉａ：Ｗ．Ｂ．ＳａｕｎｄｅｒｓＣｏｍｐａｎｙ，１９７１．［３７］　周婷，彭少麟．边缘效应的空间尺度与测度．生态学报，２００８，２８（７）：３３２２⁃３３３３．

ｕｓｉｎｇａｃｏｍｐｒｅｈｅｎｓｉｖｅａｎｄｓｅｑｕｅｎｔｉａｌｃｌａｓｓｉｆｉｃａｔｉｏｎｓｙｓｔｅｍａｎｄｒｅｍｏｔｅｓｅｎｓｉｎｇｄａｔａ．ＧｒａｓｓｌａｎｄＳｃｉｅｎｃｅ，２０１２，５８（４）：２０８⁃２２０．１９８９，５９（２）：５９⁃１２３．

［３８］　杜怀玉，赵军，师银芳，车彦军．气候变化下中国潜在植被演替及其敏感性．生态学杂志，２０１８，３７（５）：１４５９⁃１４６６．

１２４２⁃１２５０．

［３９］　李婷，刘康，胡胜，包玉斌．基于ＩｎＶＥＳＴ模型的秦岭山地土壤流失及土壤保持生态效益评价．长江流域资源与环境，２０１４，２３（９）：

ｈｔｔｐ：／／ｗｗｗ．ｅｃｏｌｏｇｉｃａ．ｃｎ