生态系统服务是人类在自然生态过程中直接或者间接获得的各种生态产品和生态服务总和^[1]。自生态系统服务提出以来，土地变化被认为是生态系统服务价值损益变动的主要驱动因素^[2]。近十几年，人类经济活动正在把大量的适宜开垦的草地和林地转变为可供生产的耕地^[3]。在人类活动不断加剧背景下，生态系统出现了一定程度的退化问题^[4]，如何实现生态系统持续地向外界提供生态系统产品，已经成为世界各国面临的共性挑战^[5]。系统性地评估土地利用变化与生态系统服务价值损益之间的关系，设计和制定适当的流域土地干预政策，对于保证流域生态安全的土地利用方案具有重要意义。20世纪70年代，LULCC变化对生态系统服务的广泛影响得到诸多学者的印证，尤其是在生物多样性^[6]、支持人类需求的能力方面^[7]。Constanza^[8]提出全球ESV量化方法和模式后，土地利用数量与结构对生态系统服务价值产生的影响得到学术界的重视。土地变化的驱动因素开始用于评价土地利用强度与生态系统服务之间的权衡关系^[9]，这一研究逐渐由全国尺度^[10]向较小的空间尺度进行应用^[11]。现有成果集中在两个方面：①基于土地利用变化趋势对生态系统服务价值的影响进行定量评估^[12]。②根据经济社会驱动因子，选择对区域生态系统服务价值变化具有较大的影响因子，模拟未来土地利用变化趋势，据此分析生态系统服务价值损益变化^[13]。从现有研究方法来看，受到广泛关注的方法分别是经验统计模型^[14]、元细胞自动机、系统动力模型及CLUE-S模型等^[15]，其中CLUE-S模型受到学者广泛应用，在一些中小尺度的地类转换关系和空间分配中取得较好的模拟效果。

额尔齐斯河流域中国段位于“丝绸之路经济带”核心区和北通道，自然生态环境优越，拥有大面积的天然森林和草地，属于我国重要的水源涵养区，对于中国树立生态文明大国形象以及巩固北冰洋权益，确保丝绸之路经济带核心区水安全，支撑中国向全球提供更多生态产品具有重要意义。近年来，由于全球气候变化和人类不合理活动的影响，额尔齐斯河流域的土地利用呈现出怎么样变化？与之相联系的生态系统服务价值呈现出何种损益分异规律？上述问题直接关系到额尔齐斯河流域可持续发展。本研究基于2000年以来的5期土地利用数据，运用CLUE-S模型预测未来额尔齐斯河流域各土地利用类型变化趋势，评估流域生态系统服务价值损益的时空分异规律，旨在为流域土地资源优化管理，巩固生态系统安全提供参考。

1 研究区概况

额尔齐斯河流域位于新疆北部阿勒泰地区，发源于富蕴县阿尔泰山南坡，沿阿尔泰山南麓在哈巴河县西部流入哈萨克斯坦，注入布赫塔尔玛水库(斋桑泊湖)，出湖后流经俄罗斯与额毕河汇流，最终注入北冰洋。额尔齐斯河流域全长4 248 km(中国境内633 km)，流域面积1.64×10⁶ km²(中国境内4.53×10⁴ km²)。额尔齐斯河流域属于温带大陆性气候，受大西洋暖湿气流影响，年均降雨量180 mm，径流主要以融雪和融冰、降水补给。由于人类在流域沿岸开荒耕地、过度放牧以及建坝蓄水，造成河谷森林面积从1959年的4.99×10⁴ hm²减少到1980年的1.88×10⁴ hm²年均减少1 400 hm^2[16]，生态景观破碎化加剧，生境质量下降较为严重。额尔齐斯河流域中国境内部分流经富蕴县、福海县、阿勒泰市、布尔津县、哈巴河县以及北屯市。截至2019年，所在县市的GDP为291.88亿元，人口达到58.84万人，第一产业生产总值达到50.18亿元，粮食产量3.34×10⁵ t，牲畜存栏数量达到311.15，是新疆重要的农牧业生产基地。受到农牧业生产活动以及水资源开发利用影响，土地利用格局发生了较大程度的变化，造成生态系统结构和功能出现了不同程度的退化，严重影响了流域生态系统服务价值的供给。研究额尔齐斯河流域土地利用与生态系统服务价值对于该区域实现“两山”目标意义重大。

2 数据来源与研究方法 2.1 数据来源

本文以2000，2005，2010，2015和2018年连续5期土地利用/覆被数据为主要数据源，该数据来源于中国科学院(http://www.resdc.cn)，空间分辨率为30 m，土地利用数据格式为栅格类型，一级类型数据精度可达到95%以上，将土地划分为耕地、林地、草地、水域、建设用地、未利用土地6大类。基于草地和森林资源调查、固定样地等实测数据，所有影像的总体分类精度均高于95%，能够满足研究需求。

海拔及地形数据来自地理空间数据云的SRTMDEM数据，空间分辨率为30 m。该区城镇居民点、道路数据来源于国家基础地理信息中心。阿勒泰地区主要农作物类型、总产量、播种面积及价格等数据来源于《阿勒泰地区统计年鉴》，并结合2000—2019年《全国农产品成本收益资料汇编》估算生态系统服务价值当量。城镇居民点和道路数据，来源国家基础地理信息中心数据库。

2.2 研究方法 2.2.1 CLUE-S模型介绍

CLUE-S模型是由荷兰瓦赫特宁根大学开发，基于不同土地利用类型与多个驱动因子线性关系经验，对区域所有土地利用进行精确空间分配，适宜较小尺度的土地利用变化与预测研究，主要研究步骤如下：

(1) 土地利用数据整理。基于2000—2018年额尔齐斯河流域土地利用历史变化趋势，应用马尔科夫链模型估算了2025和2030年各类土地利用类型的面积，马尔科夫链模型是基于时间序列的随机过程，广泛用于地理学的数值未来预测，根据额尔齐斯河流域现有土地利用类型集合X(t)预测未来各地类集合状态X(t+n)，基本形式为：

$ X(t+n)=X(t) \times \boldsymbol{P}_{i j}^n $

(1)

式中：n为时间步长；P_ij为土地利用转移概率矩阵。

(2) 设置土地利用转化规则及转化弹性。将额尔齐斯河流域历史时段时内不同土地类型间发生过转化情况赋值为1，将不同地类间未发生过转化情况赋值为0^[14]。弹性系数描述了不同地类转化的稳定性，参考已有研究成果^[17]，并经过多次模拟调试，确定水域、建设用地、未利用地、耕地的弹性系数分别为0.9，0.95，0.95，0.4，0.9。选择地形因子(海拔、坡度、坡向)和区位因子(距离道路距离、距城镇距离)作为驱动力对各地类转换进行考察。表 1的回归分析显示，得到各地类在空间某点的分布概率(P_i，u)。根据分布概率、弹性系数以及迭代系数得到的总概率对各地类进行空间分配，计算公式为：

$ \operatorname{TPTOPA}_{i, u}=P_{i, u}+\operatorname{ELSA}_u+\operatorname{ITER}_u $

(2)

式中：TPTOPA_{i, u}为栅格i是土地利用类型u的总概率；P_{i, u}为Logistics回归求得的空间分布概率；ELSA_u为土地利用类型u的弹性系数；ITER_u为迭代系数。

由表 1可知，土地利用驱动与坡度和海拔关系密切，同时ROC值曲线大于0.5，表明逐步回归方法能够较好地解释研究区各土地利用类型的转换概率。

(3) 马尔科夫链模拟检验及CLUE-S模拟结果检验。应用卡方检验对马尔科夫链预测的未来各种土地利用结果进行验证。使用历史土地利用转移概率矩阵估算2018年额尔齐斯河流域不同土地利用面积，并与其现实结果进行比较(图 1)。卡方检验结果和kappa系数计算公式分别为：

$ \chi^2=\sum\limits_{i=1}^n \frac{O_i-E_i}{\bar{E}_i} $

(3)

式中：O_i为观测值；E_i为模拟值；E_i为模拟期望

$ \text { kappa }=\left(P_a-P_c\right) /\left(P_p-P_c\right) $

(4)

式中：P_a为正确模拟的比例; P_c为随机情况期望模拟比例; P_p为理想情况下正确模拟比例。

根据公式(4)得到模型卡方值0.003 4，小于χ_0.94(4)²=9.57，即预测值与观测值差异并没有得到不显著验证，马尔科夫链预测的未来各种土地利用类型面积结果准确且可靠，计算得到kappa系数为0.814 6，即模拟精度为81.46%，高于80%的检验标准。

CLUE-S能够较准确地预测额尔齐斯河流域未来土地利用变化的空间格局，结果合理可信。

2.2.2 流域生态系统服务价值评估

本研究参考谢高地^[10]、Constanza等人^[5]的研究结果，构建了额尔齐斯河流域生态系统服务价值当量因子，由于额尔齐斯河流域属于新疆维吾尔自治区重要的旅游目的地，一些特色建筑具有较高的美学价值景观，根据相关成果^[18]，将建设用地参照荒漠设定生态系统服务价值当量因子。

借鉴谢高地等人的成果，将额尔齐斯河流域生态系统服务当量因子的经济价值为流域平均粮食市场价值的1/7，相关数据由《阿勒泰地区统计年鉴》《阿勒泰地区统计公报(2000—2019年)》和《全国农产品成本收益资料汇编2000—2019年》得到，2000—2018年全国粮食平均出售价格为1 098.80~2 313元/t，同时得到全国和额尔齐斯河流域单位面积粮食产量范围分别为4.16~7.01 t/hm²，4.58~8.71 t/hm²，获得额尔齐斯河流域生物量修正因子，借鉴Zhang等人^[17]成果，额尔齐斯河流域不同土地利用类型生态系统服务价值范围在802~2 502元(表 2)。

$ V_{i, k}=\frac{1}{7} \times P \cdot Q \cdot D_{i, k} \cdot B $

(5)

式中：V_i，k为额尔齐斯河流域生态系统服务价值当量(元)；P为粮食平均价格(元/kg)；Q为粮食平均产量(kg/hm²)；D_i，k为第i类生态系统服务第k项服务功能的当量因子；B为生态系统服务当量因子调整系数，额尔齐斯河流域为0.724；

$ \mathrm{ESV}=\sum\left(A_i \cdot V_{i, k}\right) $

(6)

式中：ESV为生态系统服务价值；A_i为土地利用类型i面积(hm²)。

2.2.3 土地利用变化与额尔齐斯河流域生态系统服务损益价值变动

已有研究成果中，使用单一动态度和综合动态度分别描述了不同土地类型及区域整体土地利用变化速率趋势和剧烈程度。

(1) 单一土地利用动态度。一般情况下，单一动态度表征一个区域某一时间段内土地利用/覆被类型变化速率和剧烈程度。

$ K=\frac{U_b-U_a}{U_a} \times \frac{1}{T} \times 100 \% $

(7)

式中：k为某种土地利用的单一动态度; U_a, U_b为区域研究期初和期末该生态系统的面积; T为研究时间段长, 研究中T为年时间段; K表示土地利用/覆被类型的年变化率。

(2) 综合动态度。综合土地利用类型动态度侧重研究区域土地利用变化幅度大小，计算公式为：

$ \mathrm{LC}=\left[\frac{\sum\limits_{i=1}^n \Delta \mathrm{LW}_{i-j}}{2 \sum\limits_{i=1}^n \mathrm{LW}_i}\right] \times \frac{1}{T} \times 100 \% $

(8)

式中：LC表示研究区一定时间段内土地利用的年综合变化率；ΔLW_i-j为研究时间段内第i类土地利用类型转为非i类土地利用类型面积的绝对值；T为研究时间段(a)。

2.3 土地利用变化与生态系统服务价值损益判别

本研究利用土地利用转移矩阵^[19]和ESV损益计算公式，对研究区不同土地利用类型的土地相互间转换而形成的ESV价值损益进行定量化评估。在此基础上，运用格网对额尔齐斯河流域生态系统服务价值损益进行空间格局进行分析，刻画未来损益价值的空间规律。现有成果中，常用格网单元有250 m×250 m，500 m×500 m，1 km×1 km，3 km×3 km等。现有成果解释研究尺度精细以及最小单元对研究结果产生较大影响^[20]，研究发现相比其他大尺度，250 m×250 m网格尺度的分析结果能够充分体现额尔齐斯河流域生态系统服务价值损益的空间分异规律。本研究将额尔齐斯河流域按照250 m×250 m大小进行划分，按照自然断点法将生态系统服务价值分为5级，使用空间计量软件Geoda软件进行生态系统服务价值损益的全局自相关分析，计算全局、局部Moran’s指数以及局部G_i指数。

$ \mathrm{LP}_i=\sum A_i \times\left(\mathrm{VC}_{i, k}-\mathrm{VC}_{i, k}^{\prime}\right) $

(9)

式中：VC_{i, k}-VC′_{i, k}分别是第i种土地利用类型的第k种生态系统服务在土地利用转变前和转变后的生态系统服务价值当量价值(元/a)。

$ \mathrm{DY}_t=\frac{L_{i, t_2}-L_{i, t_1}}{L_{i, t_1}} \times \frac{1}{t_2-t_1} \times 100 \% $

(10)

式中：DY_t为土地利用动态度(%)；L_{i, t1}, L_{i, t2}分别为土地利用类型i在t₁和t₂年的面积(hm²)。

$ \mathrm{DY}_t=\left|\frac{\sum\limits_{i=1}^n \Delta L_{i-j}}{2 \sum\limits_{i=1}^n L_i}\right| \times \frac{1}{t_2-t_1} \times 100 \% $

(11)

$ \mathrm{IM}_t=\frac{\sum\left(\mathrm{LP}_{i, j}+\mathrm{LP}_{j, i}\right)}{2 \mathrm{LP}} $

(12)

式中：IM_t为土地利用变化对生态系统服务价值影响程度(元); LP_i，j为第i类土地利用类型转化为第j土地利用类型后生态系统服务价值损益(元/hm²)；LP为生态系统服务价值损益(元)。

$ I=\frac{\sum\limits_{i=1}^n \sum\limits_{i \neq 1}^n w_{i j}\left(\chi_i-\bar{\chi}\right)\left(\chi_j-\bar{\chi}\right)}{S^2 \sum\limits_{i=1}^n \sum\limits_{i \neq 1}^n w_{i j}} $

(13)

$ I_i=\frac{\left(\chi_i-\bar{\chi}\right) \sum\limits_{i \neq 1}^n\left(\chi_j-\bar{\chi}\right)}{S^2} $

(14)

$ S^2=\frac{1}{n} \sum\limits_i\left(\chi_i-\bar{\chi}\right)^2 $

(15)

$ G_i=\frac{\sum\limits_{i \neq 1}^n w_{i j} \chi_j}{\sum\limits_{i \neq 1}^n \chi_j} $

(16)

式中：I，I_i分别为全局和局部Moran's I; G_i为局部G_i统计量；x_i为网格i从2000—2030年生态系统服务价值损益(元/hm²)；χ为所有网格内生态系统服务价值损益平均值(元/hm²)；w_ij为网格i与网格j之间的空间权重；n为网格个数；S为网格i的生态系统服务价值损益方差。

3 结果与分析 3.1 额尔齐斯河流域土地利用变化趋势

2000—2018年额尔齐斯河流域耕地面积由1 719.02 km²增加到3 174.65 km²(图 2)，年均增长3.47%，增加耕地主要由草地、未利用土地转化而来，分别增加了848.48和678.57 km²。由于农业生产的利益驱动，到2030年耕地面积持续增加到4 070 km²，草地和未利用土地向耕地转化的土地类型，分别为1 279.01和951 km²。2000—2018年，林地面积呈现快速下降趋势，由2000年的7 547.26 km²减少到2018年的5 424.36 km²，年均增长率为-2.14%，到2030年林地面积减少到3 945 km²，43%的林地转化为草地，13.57%的林地转化为未利用土地，其次是3.86%的林地转化为耕地。2000—2018年，草地呈波动上升趋势，年均增长速度0.03%，但是草地与耕地、未利用土地之间的转换幅度较为剧烈。经预测，到2030年草地面积呈现波动下降趋势，分别有12%，4.64%，3.78%的草地转换为未利用土地，耕地和林地，同时到2030年11%的林地、10%的未利用土地转化为草地。2000—2018年，水域面积由664.08 km²增长到697.42 km²，年均增长速度0.19%。到2030年，水域面积为697.42 km²。水域面积大部分转向了未利用土地和草地，转出面积分别占到水域总面积的27.14%和19.98%。建设用地由2000年的80.04 km²增长到220.95 km²，年均增长速度5.8%。到2030年，建设用地规模预计达到304.83 km²，未利用土地、草地、耕地的贡献最高，其贡献比例分别达到44%，24%和13%。从未利用土地来看，由2000年的10 988.38 km²增加到2018年的11 346.37 km²，年均增速为0.18%。到2030年，未利用土地规模预计达到11 372.42 km²，林地和草地是转化为未利用土地的主要类型。

2000—2010年，土地利用综合动态度由0.13%上升到0.40%，这一时期耕地和建设用地的快速变化，其中建设用地单一动态度由2.28%上升到8.66%，耕地由0.16%上升到6.03%。在此之后，土地利用综合动态度呈趋缓态势，综合动态度在2025年达到最低(0.19%)，到2030年，综合动态度呈缓慢上升趋势(0.23%)。2000—2030年，建设用地、耕地、水域、未利用土地单一动态度呈现正增长，分别为9.36%，4.56%，0.17%和0.12%，林地和草地的单一动态度始终小于0，说明林地和草地持续地向其他土地利用类型进行转换(图 3)。

3.2 额尔齐斯河流域生态系统服务价值变化

2000—2018年，额尔齐斯河流域生态系统服务价值518.16亿元增长到1 832.44亿元，到2030年，流域生态系统服务价值波动缓慢下降趋势维持在1 777.86亿元。耕地的生态系统服务价值一直处于增加趋势，由2000年的7.74亿元增加到2018年的52.97亿元，经预测耕地生态系统服务价值到2030年达到76.13亿元。林地生态系统服务价值在2018年达到最大值346.38亿元，此后逐渐降低到2030年的256.24亿元。草地生态系统服务价值处于持续增长趋势，到2030年可达1 059.89亿元。水域生态系统服务价值由2000年的25.30亿元增加到2030年的100.57亿元。建设用地生态系统服务价值由2000年的0.02亿元增加到2030年的0.28亿元。未利用土地由2000年的72.93亿元，增加到2030年的284.76亿元(表 4)。

2000—2030年的各类生态系统服务价值呈现出正增长趋势(表 5)，其中水文调节的生态系统服务增长量最高，由153.27亿元增长到553.79亿元，其余依次为：是土壤保持(149.75亿元)＞生物多样性(131.11亿元)＞气候调节(127.67亿元)＞净化环境(112.65亿元)＞气体调节(118.81亿元)＞水资源供给(75.33亿元)＞美学景观(68.32亿元)＞原料生产(37.42亿元)＞食物生产(33.77亿元)＞维持养分循环(11.65亿元)。

然而，自2000年以来，除了食物生产、水资源供给之外，其余各类生态系统服务正增长趋势并没有一直持续，尤其自2015年以后，气体调节(12.24亿元)、原料生产(减少11.24亿元)、生物多样性维持(减少8.61亿元)、土壤保持(减少6.91亿元)、气候调节(减少7.97亿元)出现了较大浮动的减少，造成这一现象的主要原因可能是由于耕地面积的持续的增长，引起的食物生产价值持续上升，而林地面积持续萎缩则带来了的是森林生态系统提供的服务价值下降。

3.3 额尔齐斯河流域生态系统服务价值损益时空变化

2000—2030年，流域生态系统服务损益价值为263.97亿元，耕地转换为林地(13.74亿元)，耕地转换为草地(26.61亿元)，草地转换为林地(85.66亿元)，草地转换为水域(13.87亿元)，未利用土地转换为林地、草地和水域(40.16亿元、41.82亿元、21.87亿元)水域转换为草地(16.89亿元)，水域转换为未利用土地18.31亿元。2000—2030年，流域生态系统服务价值损失为108.80亿元，其中草地转换为未利用土地损失最大(34.39亿元)，所占比例达到31.61%，其次是林地转换为草地(28.49亿元)，所占比例26.18%，再次是水域转换为未利用土地(18.31亿元)，所占比例16.83%。额尔齐斯河流域生态系统服务价值的增值主要是源于其他土地类型向林地和草地转换，而损失主要是源于林地、草地以及水域向其他土地类型进行转换而造成(图 4)。

2000年额尔齐斯河流域生态系统服务价值呈现出较为均匀的空间分布，北部海拔较高的阿尔泰山系的生态系统服务价值较高。到2018年，耕地面积呈现快速上升，林地面积持续萎缩，农业生产使得南部河谷地区林地和草地转换为耕地，使流域生态系统服务价值的分布格局差异逐渐显现(图 5)。探索性空间分析结果显示，流域生态系统服务价值损益全局Moran's I指数及p值分别为0.123 002和0.000 0(表 6)，说明生态系统服务价值损益表现为显著的正相关关系。

由图 5可以看出，流域生态系统服务价值损益的热点主要集中在哈巴河县加依勒玛乡、布尔津县禾木哈纳斯蒙古民族乡、齐巴尔乡、布尔津镇、阿勒泰市阿克吐木斯克牧场、富蕴县可可托海镇、克孜勒希力克乡，青河县阿热勒托别镇等区域，表现为高—高、高—低集聚特征和热点。而生态系统服务价值损益冷点主要分布在布尔津县喀拉希力克乡、哈巴河县加依勒玛乡、萨尔塔木乡、萨尔布拉克乡，阿勒泰市切木尔切克镇、兵团181团、福海解特阿热勒乡、阿尔边乡、阔克阿尕什乡、汗德尕特蒙古族乡，富蕴县克孜勒希力克乡、库额尔齐斯镇呈现低—低、低—高的集聚特点。

4 讨论与结论 4.1 讨论

(1) 2000—2030年，额尔齐斯河流域土地利用变化导致的生态系统服务损失价值为263.97亿元，耕地向林地、草地、草地向林地转换是主要贡献的类型。随着耕地向林草的转换，原因是农业生产扩张，林地和草地向耕地转换以及缺少有效灌溉而造成弃荒行为，价值损益在全局和局部尺度上呈现出较为显著的正向空间自相关特点。额尔齐斯河流域生态系统服务损益变化与其他区域存在较大差别^{[9, 21]}，额尔齐斯河流域土地变化主要驱动力是源于农业生产无序扩张，大量林地和草地转化耕地和用来放牧。一直以来，人类活动在额尔齐斯河流域LUCC之间的联系中扮演了非常重要的角色：一方面，阿勒泰地区人口从1990年的51.28万人(自然增长率15‰，新疆地区平均水平为11.17‰)到2018年的65.95万人(自然增长率4.63‰，新疆地区平均水平为3.69‰)，人口增长直接加大了对土地的生计需求；另一方面，市场对农牧民的土地行为决策影响：面对食物，肉类的巨大市场需求，刺激农牧民将越来越的林地和草地用作种植经济作物和放牧，然而由于农产品周期性因素造成的供过于求，又造成农牧民土地弃耕和牲畜无法及时出售，引发了土地荒漠化问题。

(2) 2000—2018年额尔齐斯河流域耕地面积的持续扩大，林地持续地转化为其他土地利用类型，到2030年耕地占到流域总面积的8.41%。林地面积呈现不断萎缩的趋势，年均增长率为-2.14%(2000—2018年)，这也直接印证刘汉初等人^[22]的结论，即森林种群死亡、人为砍伐以及矿产资源开发，造成森林生态系统受到极大的生存威胁。本研究还进一步弥补孙庆祥等^[23]人并没有对森林生态系统服务价值时空变动趋势进行分析的不足之处，即得到额尔齐斯河流域林地生态系统服务价值呈倒U形变化趋势，由2000年的129.70亿元增加到2018年的346.38亿元，预测2030年处于256.24亿元左右。2018年林地生态系统服务价值大约相当于2018年阿勒泰地区GDP(284.09亿元)的1.22倍左右。

(3) 土地利用变化驱动的预测不仅要基于地表农田、森林、草地以及湿地的自然地理现状，还需要结合生态修复工程；相比其他流域^[24]，2019年以来额尔齐斯河流域山水林田湖草生态保护修复工程试点项目在该地区进行试点实施，该项工程极大地推进流域生态系统整体保护，更加重视山水林田湖草生命共同体建设。那么，生态系统服务价值损益预测还需要建立一个生态保护情景，即更多的土地被用来提供更多的生态系统服务产品。在未来研究中，需要考察生态保护工程对土地利用驱动的影响，将宏观遥感数据与家庭调查、人口数据、土地经济调查等其他微观数据结合起来，以便更好地考察土地利用转型与生态系统服务价值变动，以此来制定长期的土地管理政策。

(4) 从不同土地类型来看，额尔齐斯河流域的水域、森林、草地提供了大量的生态系统服务产品和服务。由于全球气候变暖大背景下，流域北部山区地表覆被发生显著变化，尤其体现在冰川融雪加快，使得额尔齐斯河流域水源涵养服务功能得到一定程度提升。但这一提升有可能会在未来带来较多的自然灾害与极端水文事件，在后期研究中还需要借助更为充实的气温、降雨数据来模拟气候变化情景下冰川积雪变化与生态系统服务增减规律。此外，气候变化与植被覆盖存在紧密联系，尤其是气温与降雨变动带来的森林、草地覆盖变化，在后续跟踪研究中还需要纳入草地、林地类型的NDVI指数等驱动因子对生态系统服务价值进行更加精准模拟和预测。

4.2 结论与建议

(1) 2000—2018年新疆额尔齐斯河流域林地持续地转化为其他土地，耕地面积的大幅度增长，这一趋势一直持续到2030年。额尔齐斯河流域生态系统服务价值在2000—2018年由518.16亿元增长到1 832.44亿元，在此后生态系统服务价值处于波动下降趋势，到2030年，额尔齐斯河流域生态系统服务价值维持在1 777.86亿元。

(2) 2000—2030年，额尔齐斯河流域生态系统服务价值损失为108.80亿元，其中草地转换为未利用土地的损失最大(34.39亿元)，所占比例为31.61%，其次是林地转换为草地(28.49亿元)，所占比例26.18%，再次是水域转换为未利用土地(18.31亿元)，所占比例为16.83%。

(3) 额尔齐斯河流域生态系统服务价值损益与南北地形存在紧密联系，呈现“北高南低”的发展格局，价值损益在全局和局部上均呈现显著的正向空间自相关。

(4) 应严格控制北部山区土地利用强度，严格执行国家设定的生态红线政策，进一步优化土地利用结构，进一步扩大退耕还林还草实施范围，调整农业产业结构，严格控制食葵等高耗水农作物的种植面积和范围。

(5) 在南部河谷地区，应加大耕地生态保护工程，禁止土地开荒行为，扩大草原生态保护奖励补助的范围，纳入流域两岸不适宜放牧的沙化草地进行人工封禁保护，精准实施草地生态修复和保护工程。