黄土高原丘陵区是我国水资源匮乏和水土流失最严重地区, 水源涵养作为该地区最重要的生态服务功能之一^[1], 对气候、水文、植被、生产力和土壤养分循环等具有直接的影响^[2]。长期以来, 受特殊的自然地理条件和传统农业生产模式的制约, 导致黄土高原丘陵区生态环境不断恶化, 水源涵养能力降低, 生态系统服务功能减弱, 严重影响和制约着黄土高原丘陵区生态和社会经济的可持续发展^[3]。由于黄土高原丘陵区生态环境的脆弱性和水源涵养功能的重要性, 自20世纪80年代以来, 国家在该地区先后开展了小流域治理、退耕还林草工程和生态补偿等一系列生态工程和政策, 该地区的生态条件和水土保持功能得到了显著的改善^[4]。但随着植被结构变化与盖度的增加, 产生了人工植被稳定性低、植物蒸腾量增大、土壤干化、流域产水减少等问题^[5-6], 由此产生的水源涵养功能及其空间变化对区域生态服务的影响引起广泛关注。目前, 国内针对黄土高原生态系统水源涵养功能, 应用多种方法开展了大量的研究。苏常红和傅伯杰^[7]利用InVEST模型对黄土高原在气候和土地利用变化情境下的水源涵养进行评估, 表明陕北黄土高原水源涵养在1990—2008年呈下降趋势。李晶和任志远^[8]基于GIS对陕北黄土高原水源涵养及其价值进行了评估, 表明1990—2000年以年均0.68%递减。包玉斌等^[9]基于InVEST水源涵养功能评价模块, 定量评价退耕还林还草工程背景下土地利用/覆被变化对研究区水源涵养的影响并进行了水源涵养空间分区, 表明2000—2010年, 随着陕北黄土高原草地、灌丛和林地面积的增加, 但水源涵养量整体呈减少趋势。王超和甄霖^[10]利用综合蓄水能力法对宁夏固原市森林水源涵养功能时空动态进行了分析, 表明在2000—2010年灌木林的单位面积水源涵养量最大, 对水源涵养量增长的贡献占93.34%。已有的研究对黄土高原生态建设与保护发挥了积极作用, 但主要集中在黄土高原典型区域退耕前后植被状况发生变化的评价, 时序较早, 没能通过水源涵养功能评估来甄别区域内水源涵养的冷热点范围和确定水源涵养的重要区域。由于黄土高原地貌类型多样^[11], 丘陵区作为主要的地貌单元, 定量地评估黄土高原丘陵区的水源涵养功能其空间分布特征, 确定水源涵养的热点和重要保护区域, 对促进该区域生态保护和水资源的合理配置, 实现可持续发展具有重要科学意义。

对水源涵养水源功能评估主要是基于水量平衡、土壤蓄水力和径流量计算, 评估的模型和方法较多, 但InVEST模型作为整合多种算法的集成式模型^[12], 自2016年发布最新的InVEST3.3.1版本后, 结合土地利用和覆被、降水和土壤深度等数据计算供给水量, 用于分析生态系统的结构和功能变化对生态系统服务及其价值的影响^[13], 在我国主要应用于不同土地利用策略的生态系统服务功能差异及对多种生态系统服务权衡/协同关系的定量辨识, 为区域土地利用规划中潜在的生态系统服务供需提供优化方案^[14]。该模型将GIS技术融合到生态系统服务评估中, 利用函数定量评估现状和未来情景下的多种生态系统服务价值, 是集精细化和定量化为一体的空间化生态系统服务功能评估方法, 在空间格局测度中具有数据易于获取和操作简单的优势, 被国内外广泛的应用^[15-16]。

目前国内对生态服务方面的研究很多, 但针对黄土高原丘陵区水源涵养的研究不多, 且主要是以功能和价值评估为主的单一性研究, 缺乏结合功能评估、热点分析和重要性分区的综合研究。因此, 本研究应用InVEST模型, 选择地处黄土高原中部定西市安定区开展水源涵养功能的定量评估, 应用空间统计方法对水源涵养功能的冷热点区的空间分布特征进行分析和热点区识别, 以明确水源涵养功能的强弱分布, 在此基础上对水源涵养功能进行重要性分级和空间分区, 确定水源涵养优先和重点保护的区域, 为科学合理地制定黄土高原丘陵区生态恢复与水源涵养保护决策提供依据。

1 研究区概况

定西市安定区位于黄土高原中部, 东经104°12′48″—105°01′06″之间, 北纬35°17′54″—36°02′40″；地貌以梁峁、沟壑、丘陵为主, 海拔在1662—2570 m之间；年均气温3—9℃之间, 多年平均降水量515.2 mm, 干燥指数在1—2之间, 属于中温带半干旱大陆性气候, 自然植被覆盖率较低, 土地利用类型以草地、耕地、林地、灌草地、建设用地和未利用土地为主, 是典型的干旱半干旱地区^[17]。总面积为363871 hm², 境内河流分布为北部的关川河, 西部的称钩河, 西南部的西河, 东南部的东河, 东部的西巩河。其中, 关川河为境内最大河流, 干流全长180 km。该区多年平均土壤侵蚀模数达到3500—6000 t km^-2 a^-1以上, 水土流失面积达3.65×10⁵ hm², 占土地总面积的80.64%, 是我国水土流失最严重的区域之一^[18]。2017年总人口达43.05万人, 农村人口22.31万人, 生产总值达到80.08亿元, 是以农牧业生产为主导产业的传统农业地区, 脆弱的自然生态环境与落后社会经济水平相互作用, 使该地区的生态经济问题十分严重^[19]。研究区地理位置图(图 1)。

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

数据主要包括土地利用、植被覆盖、土壤、水分、气候、DEM及社会经济等数据。土地利用和植被覆盖数据来源于地面监测和地理空间数据云的Landsat-8影像数据, 空间分辨率为30 m；降水、温度、日照和蒸发等数据来源于中国气象局西北区域气象中心在定西市安定区内14个气象站点2017年数据插值和中国气象科学数据共享网(http:/data.cma.cn)；DEM数据下载自地理空间数据云, 空间分辨率为30 m；土壤含水量、砂粒、粉粒、粘粒、有机质、氮磷、容重和植物根系深度等数据是在2016年10月开始到2017年10月期间, 定西市安定区33个土壤监测点获取的, 分别在阳坡和阴坡沿顺坡方向设置样带(间隔15 m), 采用网格布点法设置5个水分和土壤监测点, 平均每15 d测定一次土壤水分, 使用时域反射仪TDR (TRIME, IMKO Micromodultechnik, Germany)测定0—180 cm土层的土壤水分；土壤的其他理化指标于2016年10月和2017年10月取样2次, 在实验室进行分析。DEM数据下载自地理空间数据云, 空间分辨率为30 m；社会经济数据来源于《定西市统计年鉴》, 行政区域、道路、水系等矢量数据, 来自国家基础地理数据库。

2.2 研究方法

本研究利用InVEST模型进行分析与评价。InVEST水源涵养模型是基于栅格单元的降水、植物蒸腾、地表蒸发、植物根系深度、植物有效可利用水和土壤最大根系埋藏深度等参数计算出产水量, 包括地表产流、土壤含水量、枯落物持水量和冠层截留量^[20]。计算出年产水量后, 根据DEM计算径流路径地形指数, 利用土壤渗透性、地表径流流速系数计算径流在栅格上停留时间, 计算出水源涵养量^[21]。该模型综合考虑了不同土地利用类型的土壤渗透性、地形差异、地表粗糙程度等对地表径流的影响, 能够较好地表达流域水源涵养量的空间分布状况及其影响水源涵养量的主要因素^[22]。

(1) 产水量计算模型

(1)

(2)

(3)

(4)

式中, Y_x为年产水量(mm)；P_x表示栅格单元x的年降水量；AET_x(annual evapotranspiration)为不同土地利用类型栅格单元x的年平均蒸散发量(mm)；R_x为Bydyko干燥指数; ω_x为自然气候-土壤性质的非物理参数；Z为经验常数(ZHANG系数), 取值范围为1—10；AWC_x(available water content)为栅格单元 x 的土壤有效含水量(mm)；K_x为栅格单元 x中不同土地覆被类型的植被蒸散系数；ET₀(evapotranspiration)为参考作物蒸散系数。

(2) 水源涵养量计算模型

(5)

(6)

式中, WR(water resource)为水源涵养量(mm)；Velocity为流速系数；TI(terrain index)为地形指数；Ksat为土壤饱和导水率(cm/d)；Drainage area为流域单元集水区的栅格数量；Soil depth为土层厚度(mm)；Percent slope为百分比坡度(%)。

(3) 热点分析法

水源涵养功能的热点或冷点区域是指水源涵养功能的特征值分布聚集度高或者低的区域。热点地区要求水源涵养具有高值, 在空间统计中, 通常采用Gi^*统计方法探查流域水源涵养功能的高值或低值在空间上的聚集情况, 是一种基于距离全矩阵的局部空间自相关指标^[23]。Gi^*的统计意义可用标准化的Z值来检验^[16], 运用ArcGIS中的Hot Spot Analysis with Rendering工具对定西市安定区水源涵养功能的冷热点区域进行分析。

(7)

(8)

式中, W_ij 为斑块i与斑块j之间的空间权重矩阵；x_j为斑块j的属性值；x为所有属性的平均值；n为总斑块数。根据研究区的自然与社会经济条件, 参考相关文献, 确定冷热点区域的划分规则为：若Z(G_i^*)>2.58, 则定为极显著热点区域, 在1.65＜Z(G_i^*)≤2.58, 则定为显著热点区, 在1＜Z(G_i^*≤1.65, 则定为热点区域。Z(G_i^*)＜-2.58, 则定为极显著冷点区域, 在-2.58≤Z(G_i^*)＜-1.65, 则定为显著冷点区域, 在-1.65≤Z(G_i^*) < 1, 则定为冷点区域, 在-1≤Z(G_i^*)≤1, 则定为不显著区域。

2.3 模型输入参数确定

土地利用/覆被数据选取的成像时间为2017年9月29日, 影像数据从Landsat-8下载自地理空间数据云, 通过遥感影像目视解译的办法获；降水数据是2017年定西市安定区14个气象站点多年降水数据应用ANUSPLIN软件插值获得；土壤与水分数据由33个试验监测点获得；DEM图均由ArcGIS 10.5软件的水文分析工具(Hydrology)操作获得；地形指数(TI)是通过计算流域单元集水区的栅格数量、土壤深度和百分比坡度数据, 集水区的栅格数量数据由ArcGIS的水文分析模块中操作获得；百分比坡度数据由ArcGIS的坡度计算器工具计算得出；壤饱和导水率(K_sat)由Cosby模型计算得出；Z系数是根据研究区产水量情况以及查阅相关参考文献^{[11, 24-25]}, 确定Z系数为5；土壤有效含水量根据Zhou等^[26]的计算结果；参考作物蒸散由FAO56修正的Penman-Monteith公式计算得到^[20]；流域汇流面积、坡度、坡长通过DEM提取；流速系数以USDA-NRCS提供的国家工程手册上的流速—坡度—景观表格为基准, 乘以1000得到；植物蒸散系数是根据相关文献结合试验数据确定^{[2, 9, 27]}。

3 结果分析 3.1 产水量空间分布特征

定西市安定区位于黄河中上游, 分属祖厉河流域和渭河流域, 其中祖厉河流域占全区总面积的99.5%。在该区的北部为关川河, 西部为称钩河, 西南部为西河, 东南部为东河, 东部为西巩河。其中, 关川河为境内最大河流, 干流全长180 km, 水资源短缺和水土流失严重是该区主要的生态问题。

利用公式(1)—(4)计算出定西市安定区2017年的总产水量为3.41×10⁸ m³, 基于栅格单元的产水量介于0—460.825 mm之间；流域单位面积产水量为292.67 m³/hm²。在各流域中, 西河流域的单位面积产水量最高(162.26 m³/hm²), 西巩河流域的单位面积产水量最低(50.54 m³/hm²)；但关川河流域面积最大, 虽然单位面积的年平均产水量仅为74.03 m³/hm², 年产水量占区域总产水量的28.2%。称钩河流域面积最小, 单位面积的年平均产水量为134.42 m³/hm², 年产水量也最小, 仅占区域总产水量的8.8%。；西河子流域的面积仅占总流域面积的18.1%, 但年总产水量在5个流域中最大, 年产水量占该区域总产水量的31.1%(表 1)。由图 2中可知, 定西市安定区产水量的空间分布差异较大, 西南部地区产水量较高, 中部和东部地区的产水量较低, 单位面积产水量最高的流域, 由于面积较小, 产水量占区域总产水量的比例并不高, 流域单位面积产水量与总产水量不具正相关性。水源涵养量较高的区域植被以草地、灌草地和林地为主, 此外, 地形地貌和人类活动对产水量也有影响, 西南部属低山浅山区, 温度较低, 年降雨量相对高于北部地区, 受人类活动干扰相对较少, 对产水和土壤水源涵养过程的影响较小, 而中部和东部地区是丘陵沟壑区, 多属耕地、荒草地、盐碱地和裸地, 年降雨量相对较低, 人类活动干扰较大, 对产水和土壤水源涵养过程的影响较大。整体上产水量呈现为南高北低的分布格局。表明产水量的变化主要受植被、地貌和人类活动的影响较大。

3.2 水源涵养功能空间分布特征

(1) 水源涵养量评估

利用公式(5)—(6)计算出定西市安定区2017年基于栅格单元的水源涵养量介于0—364.541 mm之间, 流域单位面积平均年水源涵养量36.37 m³/hm², 水源涵养总量为11900×10⁴ m³。

西河流域单位面积平均年水源涵养量最高, 达54.64 m³/hm², 西巩河流域最低, 只有17.58 mm, 仅为西河流域的32.2%。西河流域水源涵养总量最高, 达3500×10⁴ m³, 占区域总水源涵养量的29.41%；称钩河流域虽然单位面积平均年水源涵养量较高, 但由于面积较小, 故水源涵养总量最低, 只有1100×10⁴ m³, 仅占区域总水源涵养量的9.24%(图 3)。由于地理位置、气候条件、植被结构在空间上存在差异, 水源涵养也表现出明显的差别。由图 4可知, 定西市安定区水源涵养功能在空间上呈由西南向东北逐渐减少的特征, 水源涵养量与产水量在空间分布态势上基本一致。从图 2与图 3的对比分析, 水源涵养量减少最大的是关川河流域, 与产水量相比下降了37%, 减少最小的是西河流域, 下降了33%。表明在干旱半干旱地区植被覆盖率低的区域, 土壤水分蒸发大, 水源涵养能力低, 而植被覆盖率高的区域, 水源涵养能力高。但这些区域主要以草地、灌草地和林地为主, 地上地下生物量大, 对水分的消耗也大。因此, 该地区水源涵养的分布格局与地形地貌、气候条件、植被覆盖度、植被结构和人类活动密切相关。

3.3 水源涵养功能冷热点空间分布特征

利用公式(7)—(8)计算定西市安定区水源涵养功能的Z(G_i^*)值, 并根据划分规则, 识别出该区水源涵养功能的冷热点在空间上的聚集情况。统计分析表明, 定西市安定区水源涵养功能冷热点区域占该区总面积的比例分别为：极显著热点10.53%、显著热点5.51%、热点4.08%、极显著冷点8.14%、显著冷点13.7%、冷点5.62%、非显著点52.42%。非显著点占总面积的一半以上, 冷点区面积高于热点区面积。由图 5可知, 水源涵养的热点区主要集中分布于西河、东河和关川河流域, 所占面积占总面积的比例达20.12%。该区域植被类型以草地和灌草地为主, 植被覆盖良好, 有利于抑制径流, 具有较好的水源涵养功能。冷点区主要集中分布于西巩河流域, 所占面积占总面积的比例达27.46%。该区域植被类型以耕地和未利用土地为主。总体上, 定西市安定区水源涵养的热点区主要集中在南部地区, 冷点区主要集中在中部地区, 冷热点分布格局呈现与降水和植被覆盖度一致性的分布特征。

3.4 水源涵养功能的重要性等级划分

为了揭示水源涵养重要性特征, 运用ArcGIS将水源涵养功能进行空间叠置分析, 运用Geometrical interval(几何间隔)的原理, 将定西市安定区的水源涵养功能划分为一般重要、较重要、中等重要、高度重要和极重要五个等级。分级结果为：一般重要区、较重要区、中等重要区、高度重要和极重要区面积分别为116007 hm²、31674 hm²、91054 hm²、119309 hm²和5 825 hm², 分别占总面积的31.88%、8.71%、25.02%、32.79%和1.60%。

由图 6可知, 定西市安定区水源涵养高度重要区的面积最大, 主要分布在关川河流域沿岸的周围地区, 土地利用类型以耕地为主, 属水源涵养的冷点区, 植被结构单一, 植被覆盖度较低, 水土流失严重, 受人类活动影响较大, 是退耕还林重点治理区域, 水源涵养功能较弱, 作为该地区水源涵养的重要生态区, 需要优化现有的植被结构, 提高植被覆盖度, 提升水源涵养功能。一般重要区和较重要区属水源涵养的热点区, 主要分布在海拔较低的西河和称钩河流域, 是天然草地和人工灌草地的主要分布区, 但草地退化较为严重, 需要加强退化草地恢复力度, 进一步提升草地的水源涵养能力。

4 结论与讨论

本研究基于InVEST模型评估了定西市安定区2017年的水源涵养功能, 分析了其空间格局及冷热点和重要性分布特征, 得出以下主要结论：

(1) 定西市安定区单位面积平均年水源涵养量为36.37 m³ hm^-2, 水源涵养总量为11900×10⁴ m³, 西河流域水源涵养量最高, 达3500×10⁴ m³, 占该区总水源涵养量的29.41%。该区水源涵养功能在空间分布上呈由西南向东北逐渐减少的特征, 其分布格局与地形地貌、气候条件、植被覆盖度、植被结构和人类活动密切相关。

(2) 定西市安定区水源涵养功能的热点区域主要集中分布在南部地区, 所占面积占该区总面积的20.12%；冷点区主要集中在中部地区, 所占面积占该区总面积的27.46%；非显著点面积占该区总面积的50%以上。冷热点分布格局呈现与降水和植被覆盖度一致性的分布特征。

(3) 定西市安定区水源涵养高度重要和极重要区面积占该区总面积的34.39%, 主要分布在关川河流域, 属水源涵养的冷点区；一般重要区和较重要区面积占该区总面积的40.59%, 主要分布在海拔较低的西河和称钩河流域, 属水源涵养的热点区, 这些区域是该地区重要的水源涵养区, 需要纳入政策重点和优先保护的区域, 优化现有的植被结构, 提高植被覆盖度, 以提升水源涵养功能。

本研究结果表明2017年定西市安定区单位面积平均年水源涵养量达36.37 m³ hm^-2, 与黄土高原其他地区采用相同方法研究的结果有较大的差异。如包玉斌等对2010年延安市和榆林市水源涵养量的评估值在150—200 m³/hm²之间^[9], 其主要原因, 一方面是评估数据来源的差异造成的, 包玉斌等的研究数据来源于面上数据, 而本研究数据主要来源于研究区的定点采样数据, 且时序上存在差异, 植被状况发生了一定的变化。另一方面是不同区域的气候、植被结构、植被恢复年限、和土壤状况有一定的差异。本研究区植被以草地和灌草地为主, 占总土地面积的64%, 耕地占33%。由于近些年来该区域大面积种植高耗水的苜蓿, 加之以柠条为主的灌木处于生命周期的成熟期, 植物耗水量较大, 导致水源涵养量减少^[28]。一些研究也表明在黄土高原地区随着植被覆盖度提高和林、灌、草成分的增加, 植被的生态需水量增加, 土壤含水量出现下降的趋势, 涵养水量的损耗加大^{[5, 9-10]}。表明植被结构和人为干扰对水源涵养具有较大的影响, 这些研究也佐证了本研究的结果。因此, 优化植被结构, 减少耕地, 选择低耗水的牧草和灌木品种是黄土高原丘陵区恢复生态和提高水源涵养能力的重要研究内容。

本研究通过水源涵养量评估, 利用空间统计方法识别该区域水源涵养功能的冷热点空间分布格局, 明确水源涵养功能的强弱分布, 客观认识不同植被类型水源涵养量的差异具有重要价值。热点区要求水源涵养量具有高值, 在该地区水源涵养的热点区主要分布在以草地、灌草地和林地为主的畜牧业生产区域, 是退耕还林还草工程实施的主要区域, 海拔相对较高、气温低、耕地少, 人为活动影响较小, 其空间分布具有随沿河流分布的特征和地形和植被效应^{[8, 29]}。而冷点区是水源涵养量的低区域, 主要分布在海拔相对较低、地势较为平坦、以农业生产为主和人为干扰较大的地区。通过冷热点分析, 对水源涵养功能进行重要性分级和空间分区, 确定水源涵养优先和重点保护的区域, 为区域生态保护和水源涵养管理提供科学指导。水源涵养高度重要区和极重要区围绕着热点区进行分布, 是该区域水源涵养的重要生态区和主要产水区, 降水较为丰富、地表水系较发达、植被覆盖较好, 具有较好的导水性能、拦蓄降水和调节径流的能力, 涵养水分能力较强^[30]。因此, 维护草地、灌草地和林地的稳定与健康, 加强对水源涵养高度重要区和极重要区的保护, 同时, 加强冷点区植被恢复力度, 提升发挥水源涵养功能, 将有利于整体提高区域生态服务水平^[31]。由于定西市安定区属典型的黄土高原丘陵区沟壑区, 地势起伏较大, 水源涵养功能的空间分布呈现异常性, 需要进一步开展水源涵养功能空间异常点的深入分析与研究, 以揭示其空间变化的特殊规律^[2]。本项研究对客观认识该地区水源涵养功能, 优化土地利用结构, 确定水源涵养的重要保护区域, 综合考虑生态系统水分利用效率, 提升水源涵养功能, 实现可持续发展具有重要意义。

InVEST模型在国内外被应用于各类政府和机构的区域规划和区域多种服务功能的综合性评价与空间格局测度分析中^[15]。本研究基于InVEST模型对黄土高原丘陵区水源涵养功能的空间化评估, 并对其冷热点和重要性分级, 其结果具有一定的科学性与合理性, 研究结果验证了InVEST模型在黄土高原丘陵区的适应性^[1]。但由于流域生态系统结构的复杂性和多样性, 在利用InVEST模型对水源涵养功能进行评估时, 模型中一些生物物理参数是以土地利类型为基本单元确定的, 评估结果的精度不高, 在未来的研究中要考虑流域上下游的水资源关系、地表水与地下水进行区分及人类直接消耗的水分, 对模型的本地化和参数验证方面进一步加强研究。