生态系统服务是指人类从生态系统获得的各种惠益^[1-4]，它是连接自然生态系统与社会经济系统之间的桥梁^[5-6]，主要包括供给、调节、文化和支持服务四大类^[4]。生态系统服务功能已经成为评估生态系统状况和变化的重要指标，并能够为有效保护与科学利用生态系统、权衡和选择不同的生态系统服务提供依据^[6]。评估生态系统服务价值是认识生态价值，评估生态资产，使绿水青山转变为金山银山的基础环节^[7]，对推动生态效益纳入经济社会发展评价体系、制定满足可持续发展和生态文明要求的政策措施都具有重要意义^[8-9]。

灌丛是以灌木为主体的植被类型，灌木植株通常为簇生，无明显主干，群落高度多小于5 m，盖度多大于30%^[10-11]，具有种类多、分布广、生产力高、生命力强等特点^[12]。在我国，灌丛主要分布在较干旱和海拔较高的地区，覆盖了约1/5的陆地面积^[12-13]。灌丛具有重要的生态功能，对我国陆地生态系统碳收支的研究显示，灌丛储存了约30%的碳，是除森林外最重要的碳汇^[14]。此外，在全球气候变暖和大规模的植被恢复行动（例如退耕还林、“三北”工程等^[15]）的影响下，灌丛分布范围的扩张和灌木入侵使得灌丛在群落演替中发挥着越来越重要的作用^{[12-14, 16]}。虽然灌丛具有多种生态价值，对保护区域生态环境也很重要^[12]，但与森林、草地等陆地生态系统相比，针对灌丛的研究较少，这可能与其定义不清晰、群落分布变化快、基础生物量数据获取困难有关^[13]。早期的研究有部分将灌丛归类到森林或是草地生态系统中进行评估^{[11, 13, 15, 17]}。例如王兵等^[15]将灌木林与经济林、竹林并列，归类为我国的重要森林资源；赵海珍等^[17]选择了山地灌丛草地生态系统为研究对象进行服务功能价值评价，将山地灌丛归为一种草地类型。受可获取数据的限制，现有的研究或是仅对某一区域的灌丛进行服务功能评估，如王伟民等^[11]选择我国西部为研究区域对灌丛的服务功能进行评估；或是在对退耕还林、生态林建设等工程的成效进行评估时，对工程范围内的灌丛进行测算^[18-19]。在国家尺度上，针对灌丛的研究很少。王庆慧等^[20]对我国2009年和2015年灌丛生态系统服务的价值量进行了评估，但其研究数据来源于各省（直辖市、自治区）政府的统计公报，不具有统一性，且缺少功能量和空间分布研究。

基于对我国灌丛生态系统服务评估研究现状的整理，本研究以生态环境部和中国科学院资助的“全国生态环境变化（2010—2015年）调查评估”项目（于2018年4月交由生态环境部审议通过并发布）数据为来源，分类结果经过验证，并采用《陆地生态系统生产总值（GEP）核算技术指南》^[21]中的GEP方法（gross ecosystem product，简称 GEP），评估2015年全国灌丛生态系统4项调节服务（水源涵养、土壤保持、防风固沙、碳固定）的功能及价值，以期为正确认识我国灌丛生态系统服务功能及空间分布情况，深入理解灌丛生态系统的重要生态价值提供数据支撑。

一.   研究区域概况

我国大部分位于北温带，小部分位于热带，南北纬度跨越近50°，气候复杂多样，植被类型丰富，是全球灌丛生态系统分布最广泛的国家之一^[12]。我国的灌丛种类繁多、生活型各异、区系组成复杂、生态适应范围也很广^[13]。本研究采用30 m空间分辨率的光学卫星遥感数据和生态系统特征相结合的分类方法，将灌丛生态系统分为常绿阔叶灌木林、落叶阔叶灌木林、常绿针叶灌木林和稀疏灌木林4类。2015年，4类灌丛生态系统的面积分别为16.52万、44.45万、0.87万、5.66 万 km²，具体的空间分布情况如图1所示。

图  1  全国灌丛生态系统空间分布情况

下载: 全尺寸图片 幻灯片

二.   研究数据与方法

本研究以“全国生态环境变化（2010—2015年）调查评估”项目为基础，根据生态系统分类结果，以遥感参量如地上生物量、植被覆盖度等为主要数据源，结合地面调查数据与统计资料等^[22]，评估全国灌丛生态系统水源涵养、土壤保持、防风固沙和碳固定4项调节服务的功能量和价值量，并明确其空间分布状况。港澳台地区由于数据缺乏没有评估。

灌丛生态系统调节服务功能评估方法依据2020年生态环境部发布的《陆地生态系统生产总值（GEP）核算技术指南》^[21]，其中功能量核算采用模型法，分别对水源涵养量、土壤保持量、防风固沙量和固碳量进行核算。核算采用的是差别化参数，来源于当地气象、水利、统计等部门，或通过实测获取；价值量的核算基于功能量及其单价，具体的核算方法包括影子工程法、替代成本法和恢复成本法。工程造价及维护成本等数据来自国家发展和改革委员会、水利部等部门发布的工程预算依据，并根据2015年价格指数折算得到^[21，23]。其他核算参数，包括土壤容重、土壤沙化覆沙厚度、泥沙淤积系数、土壤中污染物纯含量来源于已发表的文献^{[7, 8, 24-26]}。

一   水源涵养评估

1   功能量核算

水源涵养量采用水量平衡法，通过水量平衡方程计算。公式为：$ {Q}_{w}=\displaystyle\sum\nolimits_{i=1}^{n}{A}_{i}\times ({P}_{i}-{R}_{i}-{\mathrm{E}\mathrm{T}}_{i})\times {10}^{-3} $。式中：Q_w为水源涵养量（m³/a）；A_i为i类生态系统的面积（m²）；P_i为产流降雨量（mm/a）；R_i为地表径流量（mm/a）；ET_i为蒸散发量（mm/a）；i为核算区内生态系统类型，i = 1，2，3，…，n。

2   价值量核算

水源涵养价值表现在蓄水保水的经济价值，运用影子工程法。公式为：$ {V}_{w} = {Q}_{w} \times {C}_{w} $。式中：V_w为水源涵养价值（元/a）；Q_w为核算区内总的水源涵养量（m³/a）；C_w为水库单位库容的工程造价及维护成本（元/m³）。

二   土壤保持评估

1   功能量核算

土壤保持量的估算基于修正的通用水土流失方程：$ {Q}_{s}=R\cdot K\cdot L\cdot S\cdot \left(1-C\cdot P\right) $。式中：Q_s为土壤保持量（t/a）；R为降雨侵蚀力因子，用多年平均年降雨侵蚀力指数表示；K为土壤可蚀性因子，通常用标准样方上单位降雨侵蚀力所引起的土壤流失量（t/a）来表示；L为坡长因子（无量纲）；S为坡度因子（无量纲）；C为植被覆盖因子（无量纲）；P为水土保持措施因子（无量纲）。

2   价值量核算

土壤保持价值主要包括减少面源污染和减少泥沙淤积两个方面，运用替代成本法核算。公式为：${{{V}}_{{{s}}}} = {V_{1}} + {V_{2},}{V_{1}} = {\rm{\lambda \times (}}{Q_{s}}/\rho ) \cdot {c}$，$ {V_{2}} = \displaystyle\mathop \sum \nolimits_{{{i = 1}}}^{{n}} {Q_{s}} \times {c_i} \times {p_i} $。其中：V_s为生态系统土壤保持价值（元/a）；V₁为减少泥沙淤积价值（元/a）；V₂为减少面源污染价值（元/a）；λ为泥沙淤积系数；Q_s为土壤保持量（t/a）；ρ为土壤容重（t/m³）；c为单位水库清淤工程费用（元/m³）；c_i为土壤中污染物（如氮、磷）的纯含量（%）；p_i为单位污染物处理成本（元/t）；i为土壤中污染物种类数量，i = 1，2，3，…，n。

三   防风固沙评估

1   功能量核算

防风固沙量通过潜在风蚀量与实际风蚀量的差值进行估算。公式为：$ {Q}_{f}=0.169 \; 9\times (\mathrm{W}\mathrm{F}\cdot \mathrm{E}\mathrm{F}\cdot \mathrm{S}\mathrm{C}\mathrm{F}\cdot $$ {K}')^{1.371 \; 1}\cdot (1-{C}^{1.371 \; 1}) $。式中：Q_f为防风固沙量（t/a）；WF为气候侵蚀因子（kg/m）；EF为土壤侵蚀因子（无量纲）；SCF为土壤结皮因子（无量纲）；K'为地表糙度因子（无量纲）；C为植被覆盖因子（无量纲）。

2   价值量核算

防风固沙价值由减少的沙化土地面积乘以治沙工程单价得到，运用恢复成本法。减少的沙化土地面积根据防风固沙量和土壤沙化覆沙厚度计算。

式中：V_f为防风固沙价值（元/t）；Q_f为防风固沙量（t/a）；ρ为土壤容重（t/m³）；h为土壤沙化覆沙厚度（m）；c为单位治沙工程的成本或单位植被恢复成本（元/m²）。

四   碳固定评估

1   功能量核算

选用固定CO₂的量作为灌丛生态系统固碳功能的评价指标，固碳速率法评估公式为：$ {Q}_{c}= $$ {M}/ {N}\cdot \mathrm{S}\mathrm{F}\cdot \left(\mathrm{F}\mathrm{C}\mathrm{S}\mathrm{R}\cdot (1+\beta )\right) $。式中：$ {Q}_{c} $为灌丛生态系统CO₂总固定量（t/a）；M代表CO₂的分子质量；N代表碳的分子质量；$ {M}/{N} $ = 44/12是碳转化为CO₂的系数；FCSR为灌丛的固碳速率（t·hm⁻²·a⁻¹）；SF为灌丛面积（hm²）；β为灌丛土壤固碳系数。

2   价值量核算

生态系统固碳价值采用替代成本法（造林成本法）核算。公式为：$ {{V}}_{c} ={Q_c}\cdot {C_c} $。式中：V_c为生态系统固碳价值（元/a）；$ {Q_c} $为生态系统固碳总量（t/a）；$ {C_c} $为CO₂价值（元/t）。

三.   研究结果与分析

一   灌丛生态系统服务功能量

1   全国总体情况

2015年，全国灌丛生态系统水源涵养量为1 652.89 亿m³，单位面积水源涵养量为24.49 万m³/(km²·a)；土壤保持量为267.33 亿t，单位面积土壤保持量为396.05 t/(hm²·a)；固碳量为1 289.96 Tg，固碳能力为1 911.06 g /(m²·a)。这3项生态系统服务的空间分布具有一致性，大体上呈现西南、华中地区高，东北、华北地区低的格局。发挥这3项生态服务的灌丛主要位于太行山脉、黄土高原、秦岭山脉、横断山脉、藏东南和云贵高原等地（见图2a，b，d）。

图  2  全国灌丛生态系统调节服务功能空间格局

下载: 全尺寸图片 幻灯片

全国灌丛生态系统防风固沙量为22.10 亿t，单位面积固沙量为4 437.82 t/(km²·a)。发挥固沙生态服务的灌丛集中分布在浑善达克沙地、吕梁山和太行山所处的山西高原、鄂尔多斯高原、河西走廊、准噶尔盆地和塔里木盆地四周的天山−昆仑山山脉等区域（见图2c）。

2   各省（自治区、直辖市）的情况

1）水源涵养。2015年，在全国31个省（自治区、直辖市）中，灌丛生态系统水源涵养量最高的是四川，为274.14 亿m³，约占全国总量的16.59%；其次是湖南和云南，分别为208.50 亿m³和167.98 亿m³，各占全国总量的12.62%和10.17%；水源涵养量高于100 亿m³的还有广西、贵州和西藏，上述6省（自治区）合计在全国水源涵养量中占比超过60%。水源涵养量较低的省（自治区、直辖市）包括天津、山东、江苏和黑龙江，均小于1 亿m³（见图3a）。

图  3  各省（自治区、直辖市）灌丛生态系统服务功能状况

下载: 全尺寸图片 幻灯片

从单位面积水源涵养量来看，水源涵养能力最强的是江西省，为71.56 万m³/(km²·a)；其次是浙江和福建，分别为70.57 万m³/(km²·a)和69.38 万m³/(km²·a)（见图3a）。

2）土壤保持。2015年，在全国31个省（自治区、直辖市）中，灌丛生态系统土壤保持量最高的是四川，为42.44 亿t，约占全国总量的15.88%；其次是陕西和云南，分别为37.64 亿t和24.74 亿t，各占全国总量的14.08%和9.26%；广西、西藏、湖南等省份的土壤保持量也较高，上述6省（自治区）的土壤保持总量在全国占比超过60%。土壤保持量较低的省（自治区、直辖市）包括天津、江苏、黑龙江和新疆，均小于0.1 亿t（见图3b）。

从单位面积土壤保持量来看，土壤保持能力最强的是海南省，为1 671.30t/(hm²·a)；其次是浙江和福建，分别为1 170.30 t/(hm²·a)和1 091.78 t/(hm²·a)（见图3b）。

3）防风固沙。2015年，灌丛生态系统防风固沙服务分布在15个省（自治区、直辖市）中，其中新疆和内蒙古的防风固沙量最高，分别为14.92 亿t和5.02 亿t，在全国防风固沙量中占比超过90%；其次是陕西、甘肃、河北、山西和宁夏，分别为0.54亿、0.39亿、0.34亿、0.24亿、0.22 亿t；其余省份防风固沙量较低，均低于全国总量的1%（见图3c）。在被核算的省（自治区、直辖市）中，从单位面积固沙量来看，防风固沙能力最强的是新疆，为20 317.54 t/(km²·a)；其次是内蒙古和宁夏，分别为17 388.27 t/(km²·a)和6 870.74 t/(km²·a)（见图3c）。

4）碳固定。2015年，在全国31个省（自治区、直辖市）中，灌丛生态系统固碳量最高的是四川，为229.29 Tg，占全国总量的17.80%；其次是西藏和云南，分别为151.52 Tg和129.97 Tg，各占全国总量的11.76%和10.09%。年固碳量超过5%的还有陕西、广西、湖南、贵州和湖北，上述8省（自治区、直辖市）的固碳总量超过全国的70%。年固碳量较低的省（自治区、直辖市）包括天津、山东、江苏和宁夏，均小于1 Tg（见图3d）。从固碳能力来看，碳固定能力最强的是浙江省，为3 465.34 g /(m²·a)；固碳能力超过3 000 g /(m²·a)的还有广西、江西、福建和湖北（见图3d）。

二   灌丛生态系统服务价值量

1   全国总体情况

2015年，全国灌丛生态系统服务功能价值总量为22 388.15亿元。其中，水源涵养功能的价值量最高，为13 388.37 亿元，占比高达59.80%；其次是碳固定功能，价值量为4 979.26 亿元，占比为22.24%；再次为土壤保持功能，为3 386.00 亿元，占比为15.12%；防风固沙功能的价值量最低，为634.52 亿元，占比仅为2.84%（见表1）。

表  1  全国灌丛生态系统调节服务价值量

服务功能	核算指标	功能量	单位	价值量/亿元	总价值/亿元	比例/%
水源涵养	水源涵养量	1 652.89	亿m3	13 388.37	13 388.37	59.80
土壤保持	减少泥沙淤积	267.33	亿t	846.62	3 386.00
	减少面源污染（氮）	267.33	亿t	1 730.97		15.12
	减少面源污染（磷）	267.33	亿t	808.41
防风固沙	防风固沙量	22.10	亿t	634.52	634.52	2.84
碳固定	固碳量	1 289.96	Tg	4 979.26	4 979.26	22.24
合计					22 388.15	100.00

下载: 导出CSV 

| 显示表格

2   各省（自治区、直辖市）的情况

2015年，在全国31个省（自治区、直辖市）中，灌丛生态系统服务功能总价值量最高的是四川，为3 643.27 亿元，占比达16.28%；总价值量高于2 000亿元的还有湖南和云南，分别为2 282.78 亿元和2 175.73 亿元，占比分别为10.20%和9.72%；除上海外，总价值量较低的省（自治区、直辖市）还有天津、山东和江苏，均低于10 亿元，占比不足0.05%（见图4）。

图  4  各省（自治区、直辖市）灌丛生态系统调节服务价值量及占比

下载: 全尺寸图片 幻灯片

从单项服务功能的价值量分析，水源涵养价值量较高的省（自治区）是四川、湖南、云南、广西和贵州，分别为2 220.55亿、1 688.83亿、1 360.63亿、1 341.80亿、1 057.81亿元；较低的是天津、山东、江苏和黑龙江，分别为1.08亿、4.86亿、6.66亿、7.15亿元。土壤保持价值量较高的省是四川、陕西和云南，分别为537.57亿、476.80亿、313.41 亿元；较低的是天津、江苏和黑龙江，分别为0.34亿、0.65亿、0.90 亿元。防风固沙价值量较高的省（自治区）是新疆和内蒙古，分别为428.46 亿元和144.16 亿元；较低的是云南和四川，分别为0.01亿、0.08 亿元。碳固定价值量较高的省（自治区）是四川、西藏和云南，分别为885.08亿、584.87亿、501.69 亿元；较低的是天津、山东和江苏，分别为1.81亿、8.34亿、8.71亿元（见图4）。

四.   结论与讨论

本研究基于“全国生态系统变化（2010—2015年）调查评估”项目研究结果，对2015年我国灌丛生态系统4项调节服务的功能量和价值量进行评估：2015年全国灌丛生态系统的水源涵养量为1 652.89 亿m³，土壤保持量267.33 亿t，防风固沙量22.10 亿t，固碳量1 289.96 Tg。灌丛生态系统的4项主要调节服务价值总量为22 388.15亿元，其中水源涵养价值13 388.37 亿元，土壤保持价值3 386.00 亿元，防风固沙价值634.52 亿元，碳固定价值4 979.26 亿元。

灌丛在全国生态系统类型中面积占比仅为7.1%，但在提供调节服务功能中占有重要地位。灌丛是我国仅次于森林的陆地生态系统土壤保持和碳固定的主体。2015年灌丛提供的土壤保持量占全国总量的21%左右，提供的固碳量占14%左右，单位面积土壤保持量和固碳能力也是除森林外最高的。灌丛提供的水源涵养量仅次于森林和草地，占全国总量的15%左右，但单位面积水源涵养量较低，高于草地，低于森林和湿地。灌丛在防风固沙功能中占次要地位，防风固沙量占全国总量的7%左右，但单位面积固沙量相对较高，仅低于草地。

2015年，全国灌丛生态系统价值量表现为水源涵养 > 碳固定 > 土壤保持 > 防风固沙，价值量的分布在全国尺度上呈现西部高于东部、南部高于北部的格局，价值量最高的省份集中在西南地区，这种格局是多因素综合作用的结果。生态系统服务功能价值量与生态系统类型、质量、生产力等密切相关^[27]，我国的灌丛生态系统中面积最大的是落叶阔叶灌木林，占比高达65.85%，还有常绿阔叶灌木林（占24.47%）、常绿针叶灌木林（占1.30%）和稀疏灌木林（占8.38%）。灌丛面积和类型在各省（直辖市、自治区）间的不均匀分布是影响价值量分布格局的主要因素。灌丛群落的生物量作为估算灌丛碳储量的重要指标，同时也是影响生态系统服务功能的重要因素，生物量和碳库的空间分布规律与生态系统服务功能的分布应具有一定的相似性^{[20, 27-28]}。杨弦等^[29]对中国北方温带灌丛生物量的分布进行了研究，发现东西部的生物量分布差异较大，东部温带落叶灌丛的生物量高于西部的荒漠灌丛；华北地区的灌丛生物量又略高于东北地区。胡会峰等^[12]估算了我国6种主要灌丛植被的碳储量，发现云、贵、川3省既是我国灌丛的主要分布区域，又是灌丛的重要碳库。海拔、温度、降水、坡度、坡向、土壤结构等自然条件在各省（直辖市、自治区）间差异很大，它们通过影响灌丛的植被质量、群落生产力、健康状况等来间接影响区域生态系统服务功能的强弱和价值量的大小^[27-28]。

需要指出的是，考虑到不同区域的灌丛类型、地理位置、土壤条件等对价值量评估可能造成的影响，本研究在功能量评估的过程中对数据进行了本地化处理，采用了差别化的核算参数。但由于数据可获得性的限制，在全国层面采用统一的而非体现区域差异性的价格参数，可能对评估结果带来一定误差。另外，在选择评估对象时并未涵盖灌丛生态系统提供的所有服务，仅针对主要的4项调节服务，是由于以下两方面的原因：其一，调节服务的价值最高，占各类生态系统服务功能价值的70%^[30]，其评估结果对自然生态系统与社会经济系统协调发展的政策制定具有重要意义，对调节服务评估的需求最为迫切^[6]；其二，供给、文化和支持服务更加依赖于完整的自然生态系统，较难将不同生态系统类型分别评估，其中产品提供、美学、艺术、教育等服务功能受市场价格和公众支付意愿的影响较大，评估易受主观性影响^[6]。因此，确定在国家尺度上具有代表性的核算参数和完善生态系统服务评估指标体系将是后续生态系统服务功能评估工作的重点。同时，本研究也将在“全国生态系统变化（2015—2020年）调查评估”项目完成后，增加对灌丛生态系统服务功能量和价值量的动态变化分析，以更全面地评估灌丛生态系统蕴含的潜力。