高 超，赵 军，李传华，魏 伟，王建邦，康重阳

（西北师范大学地理与环境科学学院，甘肃 兰州 730070）

气候变化和人类开发等多因素影响，可能导致草地生态环境退化，直接影响区域生态、经济和社会的可持续发展[1-2]。分析草地较长时间序列的生长情况，并对同期草地提供的生态系统服务进行评估，探索草地覆盖状况对生态系统服务的影响，可为认识草地生态系统服务规律、提升生态系统科学管理能力提供理论参考[3]。

目前，生态系统服务是学术界的研究热点之一[4-5]，前人从不同角度针对植被展开研究，但已有研究多基于特定时间段来分析植被变化的空间格局，或以较短时序分析植被的变化过程和特征，也有学者从生态恢复与生态效益方面进行探讨[6-9]。植被作为生态系统服务产生和传递的重要基础，气候变化和人为扰动的影响结果直接体现在植被上。因此，植被的生长状况势必对区域内生态系统服务产生重要影响，在分析生态系统服务变化时，有必要明确植被变化对其产生的影响和趋势[10]。

石羊河流域是干旱区内陆河流域，是中国西北生态安全屏障的重要组成部分，水资源、水土保持及植被恢复与保护等问题是区域生态和社会经济发展的重要指标。然而石羊河流域生态问题严重，水土资源作为区域的生态基础至关重要，水土资源不足引起的草地生态环境退化势必影响草地的固碳服务，进而加重人工固碳或减少排放的任务量。因此，本研究选择草地的产水量、土壤保持、固碳服务等主要生态系统服务功能，探讨草地植被覆盖变化对生态系统服务的影响，以期为石羊河流域草地植被生态系统保护和科学管理提供参考。

1 数据来源与评估方法

1.1 研究区概况

石羊河流域地理位置为 101°41′-104°16′ E、36°29′-39°27′ N，总面积 4.16 × 104km2，平均海拔约2 000 m，属于典型的大陆温带干旱气候，其上、中、下游的气候特征具有明显分异：上游祁连山高寒半干旱湿润区，年均气温约0 ℃，常年平均降水量为 300～600 mm，年蒸发量为 700～1 200 mm，是流域主要水源供给地，该区草地主要为高山和亚高山草甸，植被茂密，水土保持及碳固定能力较强；中游走廊平原干旱区，是流域内的主要灌溉农业区，年均温不足7.8 ℃，年降水量150～300 mm，年蒸发量 1 300～2 000 mm，干旱程度明显大于上游区域，且为流域的主要耗水区，该区草地主要为荒漠草原，水土保持能力较差；下游温暖干旱区，主要为民勤绿洲和荒漠，多年平均气温为7.6 ℃，年均降水量不足150 mm，年蒸发量 2 000～2 600 mm，区域降水少且蒸发强，属于重度缺水区，成为生态环境问题最为突出的区域。综合来看，中下游草地均属荒漠草原和干草原植被类型。

1.2 数据来源

植被覆盖度 (fractional vegetation cover, FVC)以MODIS-NDVI产品数据计算获得，归一化植被指数 (normalized difference vegetation index, NDVI)产品空间分辨率为250 m，时间分辨率为16 d，以最大值合成法得到年数据。固碳服务以MODIS-NPP产品数据获得，NPP产品空间分辨率为1 000 m，时间分辨率为1 a。以上MODIS数据均来自https://modis.gsfc.nasa.gov/。产水服务与土壤服务采用InVEST模型计算，所使用的基础数据如表1所列，对获得的原始数据进行处理和计算得到模型的输入参数。

1.3 评估方法

本研究以NDVI计算植被覆盖度，基于InVEST模型的产水量模块、土壤侵蚀模块，基于植物光合作用方程，定量评估石羊河流域草地生态系统的产水服务、土壤保持服务和固碳服务的时空变化特征。

气象数据等实测数据和土地利用等产品数据均满足本研究要求(表1)。植物蒸散系数、土壤保持措施因子及植被覆盖和管理因子等经验参数，以相近研究区、模型参数可借鉴性、文献的新近性为原则，通过分析大量文献获得。Zhang系数、汇水累计阈值和坡度阈值需多次模拟，使其与研究区实际情况相符。多次模拟后分别取9.433、800和25。

1.3.1 草地覆盖变化分析

植被覆盖度基于NDVI数据、以像元二分模型计算获得，公式如下：

式中：FVC为植被覆盖度，NDVIsoil和NDVIveg分别为纯裸土和纯植被对应的NDVI值，结合前人研究和具体情况，取生长季内植被NDVI最小值0.076和最大值0.818代替NDVI和NDVI[19]。

对2000-2015年草地覆盖变化进行分级，结合已有研究及研究区实地考察，将草地现状分为5类[20]。其中，植被覆盖度之差＜-20%为明显退化，在-20%～-5%为较明显退化，在-5%～5%为无变化，在5%～15%为轻微改善，植被覆盖度之差＞15%为较明显改善。

表1 InVEST 模型参数来源说明Table 1 InVEST data source

1.3.2 生态系统服务评估

1) 产水服务

产水量模块基于Budyko水热耦合平衡假设(1974年)和年平均降水量数据，以水量平衡法计算产水量。公式如下：

式中：AET(x)表示栅格单元x的年实际蒸散量，P(x)表示栅格单元x的年降水量。

水量平衡公式中，植被蒸散发(实际蒸散量与降水的比值)采用Budyko水热耦合平衡假设公式：

式中：PET(x)表示潜在蒸散量、ω(x)表示自然气候-土壤性质的非物理参数，其中PET(x)定义为：

式中：ET0(x)表示栅格单元x的参考作物蒸散，Kc(lx)为栅格单元植物蒸散系数。

式中：AWC(x)为植物可利用含水量(mm)，由土壤质地和土壤有效深度决定，Z为经验常数，反映区域降水分布及其他水文地质特征。

2) 土壤保持服务

土壤保持服务以土壤侵蚀和土壤保持量为指标，土壤保持量定义为裸地条件下的土壤侵蚀量RKLS与植被覆盖情况下的土壤侵蚀量USLE之差。公式如下：

式中：USLE为实际土壤侵蚀，RKLS为潜在土壤侵蚀，R为降水可蚀量因子，K为土壤可蚀性因子，LS为坡度坡长因子，C为作物/植被和管理因子，P为土壤保持措施因子，SD为土壤保持量。

3) 固碳服务

植被吸收空气中的二氧化碳，通过光合作用产生有机物并释放氧气，化学反应方程式为：

根据上式计算得出：植物干物质生产与CO2消耗比例为1∶1.62，因此固碳量为：

2 结果与分析

2.1 石羊河流域草地覆盖度的时空特征

2000-2015年期间石羊河流域草地覆盖度最大值不断减小，最小值却呈现增加趋势，平均值呈现先增加后减少趋势(表2)。结合2000-2015年土地利用数据的草地面积变化统计，造成这种情况的原因是草地面积在不断减少，尤其是低覆盖度草地的不断消失，使得植被覆盖度的最小值在不断增加。植被覆盖度的均值先增加后减少，2005-2015年草地覆盖度整体开始减小(最大值减小，最小值增加，均值减小)，2010-2015 年减小最为严重。上述结果表明：石羊河流域草地生态系统质量下降。

从草地覆盖变化趋势看，石羊河流域草地覆盖有显著减小，且主要集中在上游(图1)。石羊河上游地区，西南区域显著减小，所属子流域为西大河流域和东大河流域，以及西营河流域西部；中部靠南也有显著减小，北部却有较显著改善；杂木河流域、黄羊河流域及古浪河流域也呈较显著或显著减小；东南部大靖河流域呈轻微改善或较显著改善。石羊河流域中游，西部以减小趋势为主，中部及东部呈改善趋势。石羊河流域下游，草地覆盖无明显变化。

表2 2000-2015年石羊河流域草地面积及植被覆盖度变化Table 2 Changes of the grassland area and FVC in Shiyang River basin in 2000-2015

图1 2000-2015 年草地覆盖变化趋势Figure 1 Change trend of FVC in 2000-2015 grassland

2.2 石羊河流域草地生态系统服务的时空特征

2.2.1 产水量的时空特征

在2000-2015年，石羊河流域草地产水量总体呈上升趋势。从流域尺度年际变化来看，产水量呈先减少后增加趋势，在2010年达到最低766.93 mm·hm-2，在 2015 年达到最高值 926.71 mm·hm-2，多年平均值为853.79 mm·hm-2。子流域产水量的变化趋势大致与流域类似，均呈先减少后增加趋势。

从空间分布上来看，草地覆盖茂密的中上游区域差异较大。其中产水量较高的子流域为古浪河、杂木河及黄羊河流域，均达到110 mm·hm-2以上；金塔河、西大河、大靖河及中下游流域产水较少，产水量在70～110 mm·hm-2，其中金塔河流域产水量最小，低于下游流域的产水量。

从产水量的变化速率空间分布来看，2000-2015年研究区各时间段内有明显变化，大部分区域呈小幅度减少，但增加区域有较大幅度增长(图2)。2000-2005年，石羊河西部区域产水量有明显增加，流域中心明显减少。2005-2010年，部分区域呈小幅度上涨，涨幅最大的区域与前一期减幅最大的区域相同，减少最严重的区域位于流域西部边界(图2b)。2010-2015年，变化幅度最大。区域年均最大涨幅为40%；年均最大减幅达到10%。从整个时间序列来看，仅大靖河东部及中下游西北部涨幅明显；东大河流域下游及其西部接壤处减幅明显，其余区域均呈小范围浮动变化。

图2 产水量变化速率空间分布Figure 2 Spatial distribution of change rate of water wield

2.2.2 土壤保持时空特征

在2000-2015年，石羊河流域草地土壤保持总体呈波动减少趋势。从流域尺度年际变化来看，土壤保持 2015 年达到最低 115.77 t·hm-2，2000年达到最高值 200.81 t·hm-2，多年平均值为 150.64 t·hm-2。从子流域尺度来看，黄羊河流域、古浪河流域、下游流域及大靖河流域呈不同幅度持续减少，其余子流域呈不同程度的波动减少，其中下游流域因其土壤保持量小，变化不明显。2000-2005年，各子流域土壤保持均大幅减少，其中黄羊河流域、金塔河流域、杂木河流域及西营河流域减幅最大，古浪河流域与大靖河流域减幅较小。2005-2010年，杂木河流域、西营河流域及东大河流域均有明显涨幅，黄羊河流域、金塔河流域及大靖河流域有小幅度减少。2010-2015年，除古浪河流域及中下游有小幅度减少、大靖河流域有小幅度增加外，其余子流域均呈大幅度减少趋势。

从空间分布来看，除中下游流域草地覆被较少，土壤保持量极低外，石羊河上游各子流域之间土壤保持差异明显。从土壤保持变化速率空间分布来看，2000-2015年各时段均有明显变化但总体呈减少趋势(图3)。2000-2005年，石羊河流域北部土壤保持有小幅度增加，南部及东南部有不同程度减少，西大河流域、东大河流域、西营河流域及金塔河流域的中部区域减幅明显，东大河流域、西大河流域的下游及其与下游流域接壤处也有较大减幅。2005-2010年，整个流域土壤保持变化量分区明显，流域西南区域的5个子流域上游均明显增长，流域东南区域明显减少，下游流域小幅度减少。2010-2015年，流域西南部区域与前一期变化相反，流域东南区域无明显分异但总体减少，下游流域与前一期变化相仿。

2.2.3 固碳服务时空特征

2000-2015 年，石羊河流域草地固碳量总体呈减少趋势。从流域尺度年际变化来看，固碳量在2000 年达到最高值 194.25 g·m-2，2010 年达到最低值 139.01 g·m-2，多年平均值为 155.92 g·m-2。从子流域角度看，各子流域固碳量均呈减少趋势且在2000-2005年显著减少。

从空间分布来看，石羊河流域草地固碳量呈整体退化趋势，在空间上退化程度并无太大的分异特征(图4)。这说明影响草地固碳能力的扰动因子持续影响全域草地，并无局部扰动。从固碳能力上看，上游高固碳量区域不断减小，且单位面积固碳能力也在不断减小；中游与下游交界处也有部分区域固碳能力从中等向低等退化。对比同期植被覆盖度分布情况，上游草地固碳能力的退化趋势与草地植被覆盖度变化趋势大致相同：2005年以前东南部地区草地固碳功能退化最强烈，2005年以后西南部地区草地固碳功能退化超过东南部。

图3 土壤保持变化速率空间分布Figure 3 Spatial distribution of change rate of soil conservation

2.3 石羊河流域草地覆盖变化对生态系统服务的影响

为了更好地分析草地覆盖变化对各项生态系统服务的影响，以Pearson相关系数分别分析草地覆盖度与各项生态系统服务的相关性，并分上、中、下游分析草地覆盖度与各项生态系统服务的相关性，探索石羊河流域不同流域段草地覆盖对各项生态系统服务的影响差异。最后综合分析草地植被对草地生态系统服务的综合影响。

2.3.1 草地植被与产水量

草地植被与全流域和上游的产水量之间存在显著(P＜0.001)正相关，且全流域和上游的相关性强于中下游流域(表3)。从InVEST模型产水量输入参数来看，降水是影响产水量的主要因素，这与吴健等[21]研究结果相符。降水同样影响植被覆盖，石羊河流域降水量与草地分布均呈由西南向东北递减。在石羊河流域上游，水源充沛，水资源的消耗在区域环境的承受范围内；石羊河下游重度缺水且草地稀少，多为耐旱性草本植物，因此石羊河流域下游草地植被与产水量的相关性比中游更好。从植被恢复的角度来看，中游水资源的合理利用对全流域更为重要。

图4 固碳量空间分布Figure 4 Spatial distribution of carbon sequestration

2.3.2 草地植被与土壤保持

草地植被与土壤保持在全流域、上游及中游存在显著 (P＜0.001)正相关，在下游存在显著 (P ＜0.001)负相关(表4)。在InVEST土壤保持模块输入参数中，降水侵蚀力、土壤可蚀性、坡度坡长因子在一定时间内基本保持稳定，对于草地也无土壤保持措施，因此植被覆盖因子和经营管理因子易受人类活动干扰，从而影响流域的土壤保持。石羊河流域上游西部，草地覆盖减小显著，上游中部及东部部分区域有较显著改善，对比流域土壤保持变化速率发现，流域上游西部可能为植被覆盖因子和经营管理因子共同影响；流域上游中部及东部区域为经营管理因子影响土壤保持量。石羊河流域中游，西部草地覆盖与土壤保持变化均呈减小趋势；中部及东部变化趋势相反，说明经营管理因子主要影响土壤保持。石羊河流域下游，草地覆盖与土壤保持变化趋势相反，草地分布细碎，经营管理效果微弱，因此影响因子主要为草地覆盖。

2.3.3 草地植被与固碳服务

草地覆盖与固碳服务之间除2000年的中游外其余年份的流域均存在显著(P＜0.001)正相关关系，且相关系数远高于草地覆盖与产水量和土壤保持(表5)。草地覆盖度的计算与固碳服务的计算中均用到NDVI指数，草地覆盖与固碳服务之间存在共线性。因此，石羊河流域草地植被与固碳服务的相关系数均较大。对比2000—2015年石羊河流域草地覆盖与固碳服务的变化趋势，西营河流域和金塔河流域部分区域呈相反变化趋势，草地覆盖有所改善，但固碳服务量仍呈减少趋势。从CASA模型的参数分析，温度和水分可能会影响NPP量，由于产水量的减少幅度不大，温度可能是造成NPP减少的主要原因。

表4 草地植被覆盖度与土壤保持相关系数Table 4 FVC and SDR correlation coefficient fractional

表5 草地植被覆盖度与固碳服务相关系数Table 5 FVC and carbon sequestration correlation coefficient fractional

3 讨论与结论

对石羊河流域2000、2005、2010、2015年的草地覆盖、产水量、土壤保持、固碳能力进行评估，并分析草地覆盖对生态系统服务的影响，研究表明：石羊河流域草地生态系统质量下降，草地覆盖减小明显。西南部上游产水量、土壤保持及固碳能力较强，东北部下游生态系统服务能力较弱。流域内产水量有先减少后增加趋势，土壤保持和固碳能力均呈减少趋势。除草地覆盖度与流域上游产水量相关性呈显著(P＜0.001)正相关关系。石羊河流域上游草地植被覆盖减小显著，引起草地生态系统服务能力下降。闫俊杰等[22]在伊犁河谷的研究也指出，草地覆盖度退化引起了生态服务价值的损失，其中高覆盖度草地退化造成的影响最为严重，与本研究的结论一致。改善上游的草地生态系统对于提高全流域草地生态系统服务至关重要。

本研究以MOD-NDVI产品数据分析石羊河流域草地植被覆盖变化，以InVEST模型的产水量模块和土壤侵蚀模块分析研究区草地产水服务和土壤保持服务变化，以MOD-NPP产品数据分析流域草地固碳服务。在生态系统服务评估上，InVEST模型虽较其他方法有一定优势，但已有的数据很难完全满足其要求。石羊河流域及周边的气象站点较少，且土壤数据较难获得，为保证精度，土地利用数据仅使用了较权威的中国科学院数据(每5年一期)而没有自行解译更小时间间隔的数据。国内外均有研究认为，生态系统服务之间的相互关系随时间变化且存在一定规律，并表现出阶段性和差异性[23-24]。因此，若能缩小研究时间间隔，或许能更好反映生态系统服务的时间变化。