随着社会经济与人口的快速发展及其对生态系统服务的过度需求, 全球生态系统有60%的功能项正在退化^[1], 影响了人类的生存发展和区域的生态安全, 其主要原因之一是缺乏生态系统服务的有效管理^[2]。由于生态系统类型的多样化及其相互关系的复杂性, 过于强调某项生态系统服务时, 可能导致其他服务的降低^[3-4]。因此, 权衡和协调各生态系统服务, 减轻人类对自然环境的干扰, 促进生态系统功能恢复, 是提高生态系统服务管理水平的有效手段之一^[5-7]。在生态系统服务权衡的方法中, 阈值分析法^[8]可分析不同生态系统类型之间生态系统服务价值的权衡关系^[9], 但该方法本身存在不确定性, 转换自然资本为经济资本的过程较为片面静态；极值分析法假设两个或多个存在冲突的对象存在线性关系, 通过设置不同系数值使之求和最大化, 但目前考虑因素有限, 主要应用于水资源管理^[10]；InVEST模型分析法可预测不同情况下各生态系统服务的潜在变化, 以减少其与生物多样性之间的冲突^[11], 适用于大尺度研究；多目标分析法^[12]可将复杂的冲突简化, 从而达到经济、生态、环境及社会多方平衡的目标, 适用于复杂的系统分析^[6]。以上方法都可归结为土地利用类型的改造与调整^[13-14]。

流域作为社会、经济与自然相互作用强烈的区域, 水土流失、生态系统退化等问题突出^[15]。实现流域生态系统能流、物流良性循环，自然资源可持续利用，从而达到社会经济可持续发展的目的是流域生态系统管理的目标^[16]。我国小流域水土流失综合治理模式主要以治理试验观测为主, 随着实地治理经验的积累^[17-19], 形成了一系列试验示范成果与治理模式^[20-21]。但这种模式存在治理周期长、投入较大和措施推广性较差等局限性^[22]。过程模拟与情景分析可权衡并预测小流域未来状态, 计算各类治理措施下的经济、生态与社会效益^[23-24], 能快速经济地为制定不确定性的政策提供重要参考^[25]。此外, 从生态系统服务的角度优化小流域治理方案已成为当前的研究热点^[26]。主要集中在分析人类活动(治理措施、农业措施和土地利用类型变化等)对流域生态系统服务的影响, 权衡生态与经济的关系^{[8, 27-29]}等方面。

黄土高原丘壑区高低起伏, 沟壑纵横, 干旱少雨, 地形破碎, 水土流失严峻^[30], 严重制约着农业经济发展与生态系统服务供给。同时, 土地资源长期不合理利用进一步导致水土流失日渐严重及土地生产力的下降^[31]。甲积峪小流域于1998—2004年实施了黄土高原水土保持世行二期贷款项目, 2007—2011年又实施了英国赠款小流域治理管理项目^[32], 水土流失得到了有效的治理。然而, 在后续治理中各类生态系统服务的权衡问题逐步突显。据此, 本文以甲积峪小流域为研究对象, 采用多目标线性规划方法, 在生态系统服务价值评估的基础上, 设置供给服务、水土保持服务和其他服务3个目标, 根据小流域土地利用优化方向和现有治理措施设计规划变量并建立约束方程, 通过规划求解实现生态系统服务的权衡优化, 为小流域综合治理提供参考依据。

1 研究区概况

甲积峪小流域(106°57′02″—107°01′59″E, 35°14′18″—35°18′00″N)^[32]是泾河的一级支流, 属黄土高原残塬沟壑区, 位于甘肃省平凉市崆峒区东南部(图 1)。海拔1382—1785m, 面积28.2km², 平均气温8.6℃, 年降雨量511mm, 2012年总人口7713人, 以回民为主。在甘肃省水土保持区划中, 该流域属黄土高原沟壑侵蚀区。

流域内以农业种植为主, 主要农作物为小麦、玉米、胡麻、大豆, 少量紫花苜蓿, 下游地区有苹果和梨等果树。主要的畜牧养殖为牛、羊和鸡等。土地自然坡度在5°以下面积为454hm², 占16.1%；5°—15°面积1194hm², 占42.3%；15°—25°面积713hm², 占25.3%；25°—35°面积296hm², 占10.5%；35°以上面积163hm², 占5.8%。地貌变化大, 沟壑纵横, 植被稀少, 地表为黄土覆盖。由于长期冲刷侵蚀, 地貌被分割成塬、梁、峁和台等多级阶状。“Ⅴ”型沟极为发育, 沟壑密度1.34km/km²。多年平均土壤侵蚀模数6600t km^-2 a^-1, 径流模数63000m³ km^-2 a^-1。

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

遥感数据为2012年GeoEye影像, 分辨率为0.5m。DEM与国际标准分幅地形图来源于国家基础地理信息中心, 比例尺为1:50000。遥感影像已经过辐射校正、几何校正等预处理。由于小流域整体为长条形, 面积小, 采用监督分类和野外检验相结合的方法解译获得各土地利用类型空间分布。粮食单产数据来源于《平凉市崆峒区统计年鉴(2012)》、文献资料及2012年甲积峪居民问卷调查等。

2.2 生态系统服务价值评估方法

小流域生态系统服务可分为供给服务(产品供给)、调节服务(固碳释氧、涵养水源、净化环境)、支持服务(保持土壤、维持营养物质循环)和文化服务(休闲旅游和科研教育)^[1]。如表 1所示, 供给服务价值采用市场价值法评估, 调节服务和支持服务采用影子价格法、影子工程法、替代成本法和机会成本法评估。由于甲积峪水域和未利用地研究资料较少, 结合研究区实际情况, 采用谢高地当量表^[33]对其生态系统服务进行评估。研究区域无旅游景区, 且文化服务价值缺少可靠的计算数据, 本文未对其进行评估。

2.2.1 产品供给

产品供给服务价值采用市场价值法评估, 市场价格采用2012年当地收购价。主要产品包括玉米、小麦、牧草、仁用杏、苹果、农作物秸秆和活立木(乔木林)。出于水土保持的需要, 甲积峪小流域内荒草地和封山育林禁伐禁牧, 本文未计算其产品供给服务价值。

2.2.2 固碳释氧

根据光合作用方程, 按照各类植被净初级生产力计算其固碳释氧量^[34]。采用造林成本法计算固碳释氧价值, O₂单价按照352.93元/t计算, CO₂单价按照260.9元/t计算^[35]。乔木林和经济林年净初级生产力分别为2.98^[36]、1.445t hm^-2 a^-1[37]；果园和封山育林分别采用绿当量换算, 折算为乔木林的系数分别为0.47、0.73^[38]；根据“英国赠款中国小流域治理管理项目”监测结果, 荒草地年净初级生产力为1.1t hm^-2 a^-1[39]；根据调查结果, 人工草地年净初级生产力为12.50t hm^-2 a^-1(按每年刈割3次计)。研究区内粮食主要供当地居民食用, 秸秆与人工草用来饲养当地牛羊。农田作物为一年生, 人工草地的牧草为10年生, 从权衡优化的时间尺度考虑, 本文未计算农田固碳释氧价值, 人工草地以其地下部生物量进行评估。

2.2.3 涵养水源

根据各水土保持措施的保水率计算其拦蓄径流量, 采用影子工程法计算其价值。因缺少裸地径流深监测数据, 参照《平凉市崆峒区甲积峪示范小流域规划报告》中甲积峪多年平均径流深63mm进行保守估算, 各措施保水率取当地经验值。蓄水成本为3.75元/m³。计算公式为^[40]：

(1)

式中, W为保水量(m³)；P为保水定额或保水率(%)；H为径流深(mm)；Q为措施面积(hm²)。

2.2.4 净化环境

研究区净化环境服务以滞尘为主, 采用替代成本法计算。农田(坡耕地和梯田)滞尘能力为0.95t hm^-2 a^-1[41], 乔木林滞尘能力为10.11t hm^-2 a^-1[42], 果园、经济林和封山育林折算为乔木林系数分别为0.47、0.73和0.73^[38], 草地(荒草地和人工草地)滞尘能力为0.12t hm^-2 a^-1[43], 滞尘成本为170元/t^[44]。

2.2.5 土壤保持

土壤保持价值体现在减轻泥沙淤积、减少土地废弃和土壤肥力保持三方面。采用修正后的通用土壤流失方程(Revised Universal Soil Loss Equation, RUSLE)^[45]估算潜在土壤侵蚀量与现实土壤侵蚀量, 两者之差为土壤保持量。计算公式为：

(2)

式中, A_c为单位面积土壤保持量(t/hm²)；A_p为单位面积潜在土壤侵蚀量(t/hm²)；A_r为单位面积实际土壤侵蚀量(t/hm²)；R为降雨侵蚀力因子(MJ mm hm^-2 h^-1)；K为土壤可蚀性因子(t h MJ^-1 mm^-1)；LS为坡长、坡度因子(无量纲)；C为植被覆盖因子(无量纲)；P为水土保持措施因子(无量纲)。

依据相关研究成果, 获取各土地利用类型的C、P因子^[46-47]；计算得出R为51.58MJ mm hm^-2 h^-1；K为0.470t h MJ^-1 mm^-1[46]；通过GIS软件水文分析模块获得LS因子图层。

(3)

式中, E₁为单位面积减轻泥沙淤积价值(元/hm²)；24%为泥沙淤积在河道占总流失量的百分比；A_c为单位面积的土壤保持量(t/hm²)；P为减沙成本(18元/t)；

(4)

式中, E₂为单位面积减少土地废弃价值(元/hm²)；A_c为单位面积土壤保持量(t/hm²)；B为单位面积土地的机会成本(元/hm²)；ρ为土壤容重(t/m³)；D为土层厚度(0.5m)。

(5)

式中, E₃为单位面积保持土壤营养物质价值(元/hm²)；A_c为单位面积土壤保持量(t/hm²)；P为养分价格(元/t)；R为养分含量(g/kg)。氮、磷、钾养分价格采用我国化肥平均价2549元/t^[44], 有机质价格为320元/t^[48]。

2.2.6 营养物质循环

根据各植被年净增长量(gC/m²)和生态系统营养物质分配率, 计算植物体内氮磷钾积累量, 采用影子价格法计算其价值。农田年净增长量由经济产量、经济系数和含水量推算^[49-52]。人工草地种植年限较长, 主要以地上部生长与还田方式维持营养物质循环, 本文以其地上部年净增长量的氮磷钾含量进行计算。

2.3 多目标线性规划模型 2.3.1 参数选取

如图 2所示, 甲积峪较严重的土壤侵蚀主要发生在上游与中游, 这与其坡度和土地利用类型的空间分布特征一致。小流域上游 < 15°、15°—25°和>25°面积分别占上游总面积的23.09%、58.16%和18.75%, 均以封山育林为主, 坡耕地和荒草地面积分别占上游总面积的8.71%和9.24%。中游 < 15°和15°—25°面积分别占中游总面积的32.93%和46.58%, 以梯田和乔木林为主, >25°面积占中游总面积的20.49%, 以梯田和封山育林为主, 坡耕地和荒草地面积分别占中游总面积的7.38%和7.92%。下游 < 15°和15°—25°面积分别占下游总面积的56.44%和39.39%, 以梯田、建设用地和经济林为主, >25°面积占下游总面积的4.17%, 以经济林和人工草为主, 其中坡耕地和荒草地面积分别占下游总面积的6.81%和6.69%。坡度越大, 可能导致潜在的土壤侵蚀量越大, 同时中下游的坡耕地和荒草地已得到了集中治理, 余下的坡耕地和荒草地则主要集中在上游和中游, 因此上游和中游土壤侵蚀程度相对较为严重。

以甲积峪在治理中形成的上、中、下游生态功能区定位为基础, 依据坡度适宜性^[53], 分析了坡度约束条件下甲积峪土地利用的可优化方向(表 2)。以各优化方向对应的土地利用类型面积, 作为约束方程中可优化的各类生态系统面积变量。

参考Costanza等^[54]与谢高地等^[33]对生态系统类型的划分与《土地利用现状分类》(GB/T 21010—2007), 结合研究区土地利用情况与生态系统服务优化目标, 将研究区土地利用类型分为耕地、草地、果园、林地、水域、建设用地和未利用地七大类, 耕地细分为坡耕地与梯田, 草地细分为荒草地与人工草地, 林地细分为乔木林、经济林与封山育林。建设用地相关研究资料较少, 本文未计算其生态系统服务价值^[38]。根据表 2的坡度适宜性约束结果, 结合区域主要治理措施, 进一步分析甲积峪粮油产量需求、人口和规划面积约束。选取以下35个决策变量(表 3), 包括各土地利用类型规划后的面积, 以及上、中、下游不同坡度区间内土地可调整的面积。

2.3.2 规划目标

由于过度强化某类生态系统服务可能导致同一区域其他服务衰减, 在生态系统服务优化时需要对不同生态系统服务进行权衡, 如供给服务与其他服务的权衡^[55]。本文结合研究区水土流失突出的特点, 将生态系统供给服务价值、水土保持服务价值和其他服务价值分别设为目标函数F₁(X)、F₂(X)和F₃(X), 则有：

(6)

式中, V_1i为第i类土地利用类型的单位面积供给服务价值(元/hm²)；V_2i为第i类土地利用类型的单位面积水土保持服务价值(元/hm²)；V_3i为第i类土地利用类型的单位面积其他服务价值(元/hm²)；X_i为第i类土地利用类型的面积(hm²)。

对不同生态系统服务价值进行归一化处理, 得到MaxF′₁(X)、MaxF′₂(X)与MaxF′₃(X), 将3个目标函数的权重λ₁、λ₂和λ₃进行线性加权求和(λ₁+λ₂+λ₃=1), 3个权重系数均取1/3。将多目标转换为单目标求解, 则有：

(7)

2.3.3 约束条件

粮油产量需求约束。经调查得知, 2012年梯田小麦产量2632kg/hm², 小流域内居民以小麦为主食, 人均小麦占有量125.16kg/a, 自留种子用粮平均124.37kg/hm²；玉米产量7806.08kg/hm², 人均玉米占有量为391.90kg/a, 种子由购买获得；人均胡麻油占有量为30kg/a, 产量1683.82kg/hm²。

人口约束。2012年甲积峪人口7713人, 以崆峒区“十二五”规划中人口自然增长率8‰计算, 规划期末人口宜控制在8353人以内。

面积约束。优化后各土地利用类型增加面积与可调整面积(坡耕地212.45hm²、荒草地224.63hm²)相等。

依据上述粮油产量需求、人口和土地面积约束, 结合各土地利用类型调整时对坡度的适宜性约束(表 2), 设置17个约束方程(表 4)。

2.4 敏感性分析

为验证生态系统服务价值系数的准确性, 将各土地利用类型的生态系统服务价值系数分别上下调整50%, 计算其敏感度(Coefficient Sensibility, CS), 以反映生态系统服务价值对各生态系统服务价值系数的依赖程度, 计算公式^[56]为：

(8)

式中, CS为生态系统服务价值敏感性指数；ESV_i和ESV_j分别为调整前后的生态系统服务价值(万元)；VC_ik和VC_jk分别为第k类土地利用类型调整前后的生态系统服务价值系数(元/hm²)。若CS>1, ESV对VC富有弹性；CS < 1, ESV对VC缺乏弹性。CS越大, 表明VC的准确性对ESV评估的影响越大。

2.5 小流域综合治理效益评价方法 2.5.1 保水保土效益评价

根据治理所减少的水土流失量(t), 按照甲积峪实际工程成本计算其价值(元)。依据研究区多年治理经验, 保水效益按照3.75元/t计算, 减沙效益按照18元/t计算。

2.5.2 经济效益评价

经济效益为坡改梯、种植经济林和果园等土地利用调整措施产生的经济收益(元/hm²)。根据调查结果，粮食、秸秆、经济林果实和果园果实价格分别按2.00、0.60、2.50、3.50元/kg计算。

2.5.3 生态效益评价

森林覆盖率。森林植被面积占整个小流域面积的比例(%)^[57]。

固碳综合能力。指单位流域面积固定CO₂的量, 用于衡量流域内植被固定CO₂能力。公式为^[57]：

(9)

式中, Q为流域单位面积植被的固碳量(t hm^-2 a^-1)；S为流域植被总面积(hm²)；n为植被类型总数(个)；Q_i为第i种植被单位面积的固碳量(t hm^-2 a^-1)；S_i为第i种植被的面积(hm²)。

3 结果与分析 3.1 生态系统服务价值评估

如表 5所示, 2012年果园的单位面积供给服务价值最高, 为26325.74元/hm², 经济林次之；乔木林的单位面积水土保持服务价值最高, 为16425.62元/hm², 封山育林次之；水域的单位面积其他服务价值最高, 为13660.04元/hm²。单位面积生态系统服务价值排序为果园>经济林>梯田>乔木林>水域>人工草地>封山育林>坡耕地>荒草地>未利用地。甲积峪单位面积生态系统服务价值较低, 主要是由于坡耕地与荒草地的生态系统服务价值较低, 且所占面积比例较大所致。坡耕地占总面积的7.53%, 对水土保持价值的贡献仅为3.64%；荒草地占总面积的7.96%, 因水保需要禁止放牧, 短期内不利用其产品供给价值。因此, 需通过生态系统服务权衡优化将其调整为其他土地利用类型。未利用地为河流两岸裸地, 为保证河流生态系统完整性, 不对其进行调整。

表 5所示坡耕地单位面积保持土壤服务价值略高于果园, 这是由于当前果园均位于小流域下游区域, LS因子低所致。进一步结合式2计算结果分析表明, 对坡耕地而言, 改为梯田、人工草地、经济林与果园时, 供给服务、水土保持服务和其他服务价值均明显提高；改为乔木林或采取封山育林措施时, 除供给服务价值略有下降外, 水土保持服务和其他服务价值也明显提高。对荒草地而言, 改为其他各土地利用类型均可提高其供给服务、水土保持服务和其他服务价值。

为提高优化结果的合理性, 排除2012年各土地利用类型坡度坡长因子对模型权衡结果的影响, 分别提取坡耕地和荒草地在不同生态功能区、不同坡度分级斑块的LS值, 评估2022年各土地利用类型在不同斑块的水土保持服务价值(表 6), 带入目标函数(式7)进行规划求解。

3.2 最优解与土地利用变化 3.2.1 线性规划最优解

运用Excel规划求解功能, 分别设定对应目标单元和约束单元, 求得目标方程MaxF(X)最优解时的决策变量值(表 7), 据此分别对坡耕地与荒草地不同斑块优化后的同一土地利用类型面积进行加和, 得到规划前后的土地利用转移关系(表 8)。

3.2.2 土地利用转移矩阵

如表 8所示, 甲积峪生态系统服务权衡优化后, 坡耕地与荒草地全部调整为其他土地利用类型。优化前后梯田所占面积比例均为最大, 分别占小流域总面积的31.40%和34.92%；其次为封山育林, 所占面积比例由17.36%增至19.83%。调整后各土地利用类型面积增加量为乔木林＞梯田＞经济林＞封山育林＞果园, 增幅为经济林＞乔木林＞果园＞封山育林＞梯田。因乔木林可同时提供较高的水土保持服务价值和其他服务价值(表 5), 其面积增加量最大。由于荒草地和坡耕地在不同生态功能区坡度分布不同, 且坡耕地优化需满足梯田的最低面积约束, 使得两者土地利用调整方向的优先顺序不同。优化结果总体表现为坡耕地优先调整为乔木林和梯田, 荒草地优先调整为经济林和乔木林。

坡耕地调整为梯田的面积最大, 其次调整为乔木林, 分别占坡耕地面积的46.78%和40.45%。坡耕地主要集中在上中游区域, 在LS因子相同的条件下, 因乔木林的单位面积水土保持服务价值与其他服务价值最大, 使得坡耕地优化为乔木林的面积仅次于梯田。荒草地调整为经济林、乔木林和封山育林的面积比例大致相当, 分别为35.35%、33.65%和31.00%。荒草地分布在中游的面积最大, 为128.12hm², 占荒草地总面积的57.04%, 由于中游为农业发展区, 对土地利用类型调整的约束较少, 生态系统服务权衡的选择性更大, 优化结果更加多元化。

3.2.3 土地利用空间变化

如图 3所示, 规划前甲积峪梯田与坡耕地面积分别为885.70hm²和212.45hm², 分别占总面积的31.40%和7.53%。对于坡耕地, 优化后在粮食生产得到基本保障的前提下, 优先调整为具有更高水土保持服务价值和其他服务价值的乔木林。对于荒草地, 上游全部实行封山育林, 中游主要调整为乔木林和经济林(调整面积分别占中游荒草地面积的59.00%、41.00%), 下游全部调整为经济林。从上、中、下游空间分布来看, 上游坡度较陡且水土流失较为严重的坡耕地与荒草地主要优化为乔木林和封山育林, 中下游坡度较缓且土壤肥力相对较高的坡耕地与荒草地主要优化为供给价值较高的经济林和果园等。结合生态功能区划分、土地坡度及土壤肥力对生态系统服务的权衡优化结果, 可同时满足上游水土流失防治和中下游农业经济发展的需要。规划后小流域植被覆盖面积增加, 土地利用格局发生明显正向变化。

3.3 生态系统服务权衡优化效果 3.3.1 生态系统服务价值量变化

小流域生态系统服务优化后, 供给服务价值、水土保持服务价值和其他服务价值分别增加了6.24%、9.81%和17.35%, 生态系统服务总价值增加了9.23%(表 9)。三者增加量分别占生态系统服务总价值增加量的25.96%、54.42%和19.62%。乔木林与封山育林增加面积较大, 由于两者的水土保持服务和其他服务价值相对较高, 但供给价值均低于坡耕地, 使得供给服务价值的增幅相对较小。坡耕地和荒草地的其他服务价值均远低于人工草地、果园、经济林、封山育林和乔木林, 因此将两者调整为上述各类土地利用类型都会明显提高其他服务价值, 使其他服务价值增长率高于供给服务和水土保持服务价值。从生态系统服务结构来看, 调整后供给服务价值占比有所下降, 同时其他服务价值的比例增幅明显, 优化了小流域生态系统服务结构。

3.3.2 生态系统服务价值空间变化

空间分析表明(图 4), 规划目标年上、中、下游供给服务价值变化差异明显, 上游减小了18.14%, 而中、下游分别增加了9.87%和16.24%, 其变化量分别占供给服务价值总变化量的-46.16%、96.91%和49.25%。上游有较大面积土地适宜规划为供给价值低但水土保持服务和其他服务价值高的乔木林和封山育林, 因而其供给服务价值下降。中下游适合发展农业与生态经济, 调整后的土地利用类型中, 经济林和梯田所占比例较大, 因而其供给服务价值明显增加。

优化后上、中、下游水土保持服务价值分别增加8.02%、11.03%和6.83%, 其变化量分别占水土保持服务价值总变化量的20.96%、71.21%和7.82%。中游调整后的经济林、乔木林和梯田的水土保持服务价值均较高, 使其水土保持服务增加量最大。上游均调整为封山育林和乔木林, 下游主要调整为经济林和果园, 由于上、下游可调整面积均小于中游, 且果园水土保持服务价值相对较低, 因而上、下游水土保持服务价值增加量较小。

优化后上、中、下游其他服务价值分别增加23.86%、14.25%和14.76%, 其变化量分别占其他服务价值总变化量的43.54%、46.80%和9.66%。上游和中游优化为乔木林与封山育林面积较大, 两者其他服务价值均高于其他土地利用类型, 因而其他服务价值明显增高。下游可优化面积较小, 且有21.62%的面积优化为其他服务价值较低的果园, 因而其他服务价值增加量最小。

优化后上、中、下游生态系统服务价值分别增加3.00%、10.88%和12.38%, 其变化量分别占生态系统服务总价值变化量的7.40%、73.34%和19.26%。上游有大部分规划为封山育林和乔木林, 两者的水土保持服务价值较高, 但供给服务价值低, 使其总服务价值相对较低。下游尽管有部分土地优化为总服务价值最高的果园, 但优化面积仅有11.73hm², 使其总服务价值增加量所占比例较小。中游有较大面积土地优化为总服务价值仅次于果园的经济林, 使其总服务价值增加量最大。

各权衡目标价值在中游的增加量占总增加量的比例均最高, 这是由于中游主要将坡耕地与荒草地调整为梯田、乔木林和经济林, 三者都具有较高水土保持服务价值或供给服务价值。

生态系统服务不变区域为未规划区域, 包括水域、建设用地和未利用地, 以及经过前期治理、土地利用适宜的梯田、人工草地和林地(乔木林、经济林和封山育林)。

3.4 敏感性分析

由表 10可知, 规划前后各土地利用类型生态系统服务价值的敏感度均小于1。2012年敏感度排序为梯田＞封山育林＞乔木林＞人工草地＞经济林＞坡耕地＞荒草地＞果园＞水域＞未利用地, 以梯田的敏感度最大(0.41), 即当梯田生态价值系数增加1%时, 生态系统服务总价值增加0.41%；规划后各土地利用类型敏感度的排序变为梯田＞乔木林＞封山育林＞经济林＞人工草地＞果园＞水域＞未利用地, 以经济林和乔木林的敏感度增加较大, 这与经济林和乔木林面积增加较大有关^[56]。敏感性分析结果还反映了各土地利用类型对生态系统服务总价值的影响程度, 即敏感度越高, 其影响程度越大。表明, 研究区的生态系统服务总价值对各土地利用类型的生态价值系数缺乏弹性, 即研究结果可信。

3.5 综合治理效益评价

通过小流域生态系统服务与土地利用类型的权衡优化, 以2012年为基础, 到2022年研究区年径流减少11.02万t, 保水价值增加41.33万元；减少产沙4406.67t, 减少河道清淤成本7.93万元；森林覆盖率由32.66%提高到44.22%；坡改梯单产由3865.95kg/hm²提高到5219.04kg/hm², 增加产值3453.07元/hm²；新增果品601.63t, 产值达159.23万元；固碳综合能力由1.14t hm^-2 a^-1提高到1.37t hm^-2 a^-1, 提高了20.19%。此外，坡改梯和果园种植可为农田机械化耕作提供便利条件。表明, 生态系统服务和经济、生态、社会三大效益均得到了明显提高。

4 结论与讨论

本文综合了研究区土地利用的坡度适宜性、生态功能分区和黄土高原主要的小流域综合治理措施, 在生态系统服务价值评估的基础上, 采用多目标线性规划方法, 通过设置供给服务价值、水土保持服务价值和其他服务价值3个决策目标, 探讨了甲积峪小流域的生态系统服务权衡优化。

坡耕地和荒草地的生态系统服务价值较低且面积比例较大, 是研究区生态系统服务价值偏低的主要原因。通过生态系统服务权衡优化, 在满足土地利用类型的坡度适宜性和人口粮油需求等的约束条件下, 将坡耕地与荒草地主要调整为乔木林、梯田、经济林和封山育林, 少量调整为果园, 与相关研究结果一致^{[22, 58-60]}, 使得3个目标函数值均明显增加。其中, 其他服务价值增长率最高, 其次为水土保持服务价值, 供给服务价值增长率相对较低。空间分析表明, 其他服务价值在上游的增长率最高, 供给服务和总服务价值在下游的增长率最高, 水土保持服务价值在中游的增长率最高。仅有上游部分区域供给服务价值和总服务价值略有下降。优化后, 小流域的森林覆盖率、减少的水土流失量、固碳综合指数和农业产值等指标均有所提高, 为研究区坡耕地和荒草地的治理方向提供了参考依据。

在小流域的生态系统供给服务中, 粮食生产是当地居民为生存必须从事的活动, 经济林和果园等食物供给与原材料供给是居民维持生计与区域直接经济效益的体现。片面追求农业生产利益, 会使区域的水土流失问题更加严峻, 生态系统脆弱性也随之增加。各生态系统服务之间存在着复杂的相互作用, 尽可能通过权衡以减少其中的冲突并实现协作, 才能实现生态系统服务的优化。本文将供给服务价值、水土保持服务价值与其他服务价值设置为权衡优化的子目标函数, 同时参考相关研究成果^{[59, 61-64]}, 在权衡优化中侧重三者的平衡关系, 将3个决策目标的权重系数均设定为1/3, 即将供给服务、水土保持服务与其他服务视为同等重要。在约束条件不变的前提下, 将以产品供给、水土保持和其他服务价值3个子目标构建目标函数与仅以单个子目标或总服务价值作为目标函数的优化结果对比发现, 3个子目标的结果和仅其他服务价值为目标的结果极为相似, 即供给服务价值和总服务价值较低, 水土保持服务和其他服务价值较高, 同时林地面积最高, 两种方案都具有明显的优势, 但后者的供给服务价值和总服务价值低于前者。若仅以供给服务或总服务价值为目标, 优化后的供给服务和总服务价值均最高, 但调整为林地(尤其是乔木林)的面积明显减少, 使得水土保持服务和其他服务价值均最低。而仅以水土保持服务价值为目标, 则水土保持服务价值最高但优势不明显, 供给服务和总服务价值中等, 其他服务价值低；供给服务价值居中是由于梯田同时具有很高的保土、保水价值和较高的供给价值, 使得调整为梯田的面积增加而调整为林地的面积减少。因此, 选择3个子目标的结果能在满足产品供给服务价值的条件下, 通过将坡耕地和荒草地更多地调整为林地而获得较高的水土保持服务和其他服务价值, 与预期效果一致, 更符合甲积峪的实际情况。此外, 由于本文从长远效益考虑, 未计入农田的固碳释氧价值, 且人工草地仅计算了其地下部的固碳释氧价值。从黄土高原水土流失的突出问题出发, 净化环境价值仅考虑了滞尘方面, 而未考虑吸收SO₂及NO_x等污染物的价值。因此, 其他服务价值远低于产品供给和水土保持服务价值, 在一定程度上避免了其他服务价值虚高对生态系统服务权衡的影响。

土地利用不合理是导致该地区水土保持服务退化的主要原因之一, 与其他学者的研究结果一致^{[20, 65-66]}。我国在1998—2008年间对土地利用进行了一系列复杂的调整, 使得我国生态系统服务价值下降幅度低于全球平均值^[67], 效果较明显。坡耕地、荒山荒坡地和沟壑地是我国水土流失最严重的土地利用类型^[68], 研究区位于温带季风气候区, 降雨多集中于夏季, 且坡耕地和荒草地面积较大, 加剧了水土流失。如何将其调整为适宜的土地利用方式, 以减少地表径流量, 改善小流域生态环境^{[65, 69]}, 是小流域综合治理的重点和难点^[20]。本文从生态系统服务权衡优化的角度对土地利用类型在数量和空间上进行了调整, 在考虑各土地利用类型坡度适宜性的基础上, 参考《封山(沙)育林技术规程》(GB/T 15163—2004)^[70]中小流域内不适于人工造林的高山、陡坡和水土流失严重地段等地块, 对其进行封育以增加植被盖度。研究区封山育林区域大多未成林, 在计算其生态系统服务价值时生物量取乔木林的0.73^[38], 使其与乔木林有所区别。

除土地利用类型调整外, 小流域治理还可以从沟头防护、工程措施及法律法规方面着手, 以缓解人类活动与自然之间的矛盾, 针对生态系统服务权衡进行多方面综合研究是未来的重要方向^[71]。其中沟头防护措施主要从防治沟头延伸、沟底下切和沟岸扩张三方面遏制水土流失的加剧。当地居民“乱砍滥挖”的情况比较严重, 为防止过度索取导致生态环境遭到更大的破坏, 可成立村民自管小组, 通过政府发布封禁通告实行管护措施, 以此加强村民的自我约束和管理。加强道路建设, 修建小型堤坝和小型拦蓄工程也是小流域治理的有效措施之一。为了保证治理效果, 使水土保持防治机制不断完善, 依法监督和保护是必不可少的, 国家与地方相继制定了《中华人民共和国水土保持法》、《甘肃省水土保持条例》、《水土保持补偿费征收使用管理办法》、《水利工程建设监理规定》等一系列促进水土保持工作法制化、规范化和科学化建设的法律法规和管理办法。当前甲积峪在综合以上监督、防预和治理的基础上, 进行了农田种植结构调整及梯田后续产业综合开发, 力求在提高区域水土保持服务价值和其他服务价值的前提下, 实现供给服务价值最大化。