南四湖流域产水量空间格局与驱动因素分析
孙小银, 郭洪伟*, 廉丽姝, 刘飞, 李宝富
曲阜师范大学地理与旅游学院/南四湖湿地生态与环境保护山东省高校重点实验室,山东 日照 276826

第一作者简介:孙小银(1978- ),男,副教授,重庆开县人,硕士生导师,主要研究方向为流域污染生态过程与生态修复、流域生态服务功能评估与区划。E-mail:xiaoyin-sky@163.com

*通信作者简介:郭洪伟(1988- ),男,山东昌邑人,硕士研究生,主要研究方向为流域生态服务功能评估与区划。E-mail:hongwei_guo2015@163.com

摘要

流域生态系统产水服务功能的空间化和定量化评估,对流域水资源管理、优化配置以及提高流域水生态保护效率具有重要意义。以山东省南四湖流域为研究对象,基于1990—2013年土地利用、降水、蒸散以及土壤属性等基础地理数据,以InVEST模型为基础,评估和模拟南四湖流域近25 a的产水量,并采用ArcGIS分析产水量的空间分布格局以及变化趋势,探讨了降水、地形等自然地理要素以及人口、土地利用和国内生产总值(GDP)等社会经济因素与产水量空间格局动态变化之间的关系,并在此基础上划分出南四湖流域生态系统产水功能区。研究结果表明:流域产水量在空间格局上呈现出由东向西递减的趋势,东部、东北部等山区、丘陵地区产水量高,西部平原地区产水相对较低;受自然地理要素影响,流域产水量空间分布与社会经济发展水平即GDP、人口密度的空间分布格局有较大差异。近25 a来,流域产水量呈现减少趋势,且产水量峰值区域由东北部向偏南地区转移,最低值区由西部向中部地区转移。降水、海拔和坡度等地理环境与产水量的空间变化呈显著正相关,其中降水量的相关程度最强;人口、GDP等社会经济数据与产水量变化也呈显著正相关,主要原因在于城市化的发展,城市建设用地等不透水层增加,促进了流域产水量。研究结果可以为流域水资源政策制定以及社会经济发展规划等宏观决策提供科学支撑。

关键词: 生态系统服务; 产水功能; InVEST模型; 南四湖流域
中图分类号:Q149 文献标志码:A 文章编号:1000-3037(2017)04-0669-11
The Spatial Pattern of Water Yield and Its Driving Factors in Nansi Lake Basin
SUN Xiao-yin, GUO Hong-wei, LIAN Li-shu, LIU Fei, LI Bao-fu
Geography and Tourism College of Qufu Normal University/ Shandong Key Laboratory of Nansi Lake Wetland Ecological and Environmental Protection, Rizhao 276826, China
Abstract

The quantitative assessment on spatial pattern of ecosystem service of water yield in a basin has important significance for management and optimal allocation of water resource, and improvement of ecological and environmental conservation capacity. Based on basic geographical data of land use, meteorology and soil attribute in the study area, the paper used InVEST model to evaluate the amount of water yield in Nansi Lake Basin, quantitatively assessing the spatial distribution feature of water supply ability from 1990 to 2013 in the whole basin. The relationships between the dynamic spatial pattern of water yield and natural geographic elements such as precipitation, landform, and socioeconomic factors such as GDP and population were discussed. The basin was divided into several zones according to the water resources function. The result showed that the spatial pattern of water yield presents a decreasing trend from east to west, water yield being high in the east and northeast, and water yield being relatively low in west and other regions in Nansi Lake Basin. Affected by natural geographical elements, the spatial distribution of water yield is inconsistent with the spatial distribution of social economic development level, namely GDP and population. Water yield in the basin was in decreasing trend in past 25 years, and the peak area of water yield shifted from the northeast to the south area, but low yield area shifted from the west to the central region. There are significant positive correlations between the variation of water yield and physical geographical environment factors, such as precipitation, DEM and slope, and the correlation between the variation of water yield and precipitation is the strongest. There are significant positive correlations between variation of water yield and social economic factors of GDP and population, since urbanization caused the increase of urban construction land and other impervious surfaces which promoted water yield. The results of the study provide scientific support for water resources policy, social and economic development planning, and macro decision making.

Key words: ecosystem service; water yield function; InVEST modle; Nansi Lake Basin

水源供给是一项重要的生态系统服务, 近年来随着生态系统服务研究的迅速发展, 水源供给服务的空间变化及其影响机制越来越受到人们的关注[1, 2, 3, 4]。流域作为水文水资源过程的基本自然单元, 国内外学者对流域生态环境功能尤其是水源供给功能进行了大量研究[5, 6, 7, 8], 如Terrado等[9]利用 InVEST 模型在高度开发的地中海地区, 模拟研究了极端气候条件对水生态服务价值的影响效应, 潘韬等[10]应用 InVEST 模型定量估算了三江源区过去30 a水源供给量的时空变化特征, 吴哲等[11]应用 InVEST 模型模拟了不同土地利用/覆盖类型变化下海南岛产水量的空间分布, 但是, 较少有涉及流域产水功能的空间格局及其动态变化方面的研究。近年来, “ 精细化、定量化、模型化、空间化” 已经逐渐成为流域生态系统服务功能管理和评估的发展新趋势和研究热点。基于生态过程模拟的生态系统服务与权衡交易的综合评估模型(the Integrate Valuation of Ecosystem Services and Tradeoffs Tools, InVEST), 因为能够适应这一需求而在国内外得到了广泛的应用[12, 13, 14], 该模型通过地理数据空间参数的输入, 使模拟结果得以空间化表达, 更易于识别出生态系统服务的重要区域和热点区域[15], 采用 InVEST模型能够实现对生态系统服务功能多尺度综合的、动态的、可视化的评估[16, 17]。目前, 我国已有不少学者基于InVEST模型的产水模块对区域生态系统产水量空间格局开展了研究, 如林世伟等[18]对横断山“ 三江并流” 区水源供给生态系统服务空间特征进行了研究; 王雅等[19]利用4期土地利用数据研究了黑河流域土地利用变化水文效应的时空变化; 孙艳伟等[1]应用InVEST模型评估了于桥水库流域生态系统的水源供给量, 并定量分析了流域水源供给能力的空间分布格局特征; 包玉斌等[15]基于InVEST水源涵养功能评价模块, 定量评价了退耕还林还草工程背景下土地利用/覆被变化对研究区水源涵养的影响。这些研究为区域土地利用与水资源利用决策等提供了定量和空间化的参考。

南四湖是淮河以北面积最大、结构完整、保存较好的内陆大型淡水、草型湖泊湿地, 是山东省的生物宝库和重要水禽栖息地。同时, 南四湖也是南水北调东线输水干线上的重要调蓄水湖泊。就生态和经济地位而言, 南四湖在维护区域经济发展和生态保护平衡方面具有重大战略意义。然而, 随着近几十年气候变化、社会经济发展和人口的急剧增加, 流域水资源短缺问题日益突出, 研究者和相关利益管理部门迫切需要对南四湖流域产水量的空间格局及其空间结构的动态变化取得全面客观的科学认识, 以为合理地制定流域水资源管理、土地利用以及社会经济发展等方面的政策提供一个基础工具。本研究基于1990、1995、2000、2005、2010、2013年6期土地利用数据, 运用InVEST模型产水模块, 模拟南四湖流域不同时期产水量及其空间分布特征, 并结合土地利用、人均国内生产总值(GDP)、人口等社会经济数据和降水、地形等自然地理环境要素对产水量空间格局分布及变化的原因进行深入探讨, 分析1990— 2013年南四湖流域产水量的空间分布格局及其空间结构动态变化规律, 并据此划定流域水源供给重要功能区域。定量评估南四湖流域生态系统水源供给的空间分布特征及其动态变化, 对于摸清流域生态系统水源供给服务现状、指导流域水资源规划和管理具有重要的意义。

1 模型模拟与数据处理
1.1 研究区概况

南四湖流域地处山东省、河南省、江苏省和安徽省四省的交界处, 位于34° 15′ ~34° 53′ N和114° 50′ ~117° 48′ E之间, 属于淮河流域的重要组成, 行政区上完整包括山东省的菏泽、济宁、枣庄等市以及徐州、商丘、开封、宿州下辖部分县区, 总面积约3.17万km2。全流域地势自北向南降低, 南四湖湖区最低, 入湖河流共计约53条 [图1(a)]。流域地处暖温带半湿润地区, 属暖温带大陆性季风气候, 年平均气温约14 ℃, 年均降水量为700 mm左右, 年均潜在蒸散量为1 400 mm左右。依据传统的“ 土壤发生分类” 系统, 流域分布有盐碱土、人为土、淋溶土、初育土、半淋溶土、半水成土、水成土等7种土纲, 有粗骨土、风沙土等12种土类, 有棕壤性土、褐土、酸性粗骨土等20种亚类土壤, 其中半水成土面积最大, 约占66.08%。2013年全流域69.4 %为农业用地 [图1(b)], 为典型的农业流域, 是我国传统的农业耕作区。

图1 南四湖流域位置与土地利用现状Fig. 1 Location and current land use of Nansi Lake Basin

1.2 InVEST模型

本研究所用模型为InVEST 3.2.0[20], 研究选取“ Water Yield” 子模型模拟流域产水量及其空间分布格局。“ Water Yield” 模型基于水量平衡原理, 以某个栅格单元内的降雨量减去实际蒸散发量即为产水量。总产水量包括地表产流、土壤含水量、枯落物持水量和冠层截留等, 不区分地表径流、壤中流和地下径流。因此, 一定区域内的地表产水量越多, 水资源供给服务就越多。根据Zhang等[21]提出的基于Budyko曲线在模型中自动计算每个栅格的实际蒸散发量。潜在蒸散发参数, 本研究采用InVEST模型提供的“ 改进的Hargreaves公式” [17]进行计算。土壤有效含水量, 根据土壤有效水含量经验公式[22]进行计算。

1.3 数据来源及处理

本研究所采用的数据包括土地利用数据、气象要素数据、土壤基础数据、DEM数据等。土地利用数据包括6期数据, 其中1990、1995、2000、2005、2010年土地利用数据来源于中国科学院资源环境数据中心, 2013年土地利用数据 [图1(b)] 来源于遥感影像解译。以南四湖流域Landsat TM 影像为基础, 利用ENVI 4.7 进行人机交互式目视解译获得, 解译结果的精度检验采取交叉检验的方式, 精度> 90%。数据空间分辨率均为30 m× 30 m。上述土地利用数据均根据刘纪远等[23]提出的中国土地利用/土地覆被遥感分类系统对土地利用类型进行分类, 全流域共划分为6种一级地类22种二级地类。降水和气温数据来自中国气象数据网(http://data.cma.cn/)。辐射数据来自NASA官网(http://www.nasa.gov/)。土壤深度和有效含水量来源于联合国粮农组织(FAO)和维也纳国际应用系统研究所(IIASA) 构建的世界土壤数据库(HWSD)中的中国土壤数据集(1∶ 100万)。DEM数据(空间分辨率30 m)由美国马里兰大学全球土地覆被数据库提供。1995、2000、2005、2010年GDP和人口数据来自中国科学院资源环境数据中心。流域水系数据通过南四湖水利管理局提供的流域水系图(1∶ 20万)矢量化获得。

流域1990— 2013年产水量空间格局分布变化运用ArcMap-Raster calculator功能将两期产水量图层相减得到。近25 a产水量空间变化影响因素采用以下方法分析, 即将各期产水量、降水量、高程、坡度等数据重采样成1 km × 1 km分辨率, 分区统计各县区的数据。在SPSS 21.0统计软件中对每个县区产水量和降水、社会经济数据进行相关性分析。

使用ArcMap空间分析工具里面的multivariate, 选择band collection statistics, 计算当前产水量与驱动因素的相关系数矩阵, 分析当前产水量空间格局影响因素。

2 结果与讨论
2.1 模型率定

以2010年的南四湖流域土地利用、降水、蒸散发数据为基础, 调整z值(表征多年平均降水和其他水文地质特征的常数), 多次进行产水量的模拟, 并利用ArcGIS分区统计功能提取枣庄、济宁、菏泽三市的产水量, 计算其平均产水模数, 建立产水模数(y)与z值(x)之间的相关关系, 即y=-0.472x+35.21, 并与《山东省水资源公报(2010)》[24]中相同地区的水资源量进行比较。当z值等于29.4时, 模拟产水模数与实际产水模数吻合, R2高达0.996, 产水量相对误差-0.004 6%。因此, 可以认为当z=29.4时, InVEST模型产水子模块在南四湖流域的模拟效果最好。

2.2 1990— 2013年南四湖流域生态系统产水功能变化

模拟结果显示, 当前(2013年)南四湖流域平均产水深度138.67 mm, 产水总量40.82× 108 m3。1990— 2013年南四湖流域年平均产水量变化趋势比较明显, 呈现出先降后升再降低、总体趋势下降的波动变化趋势, 1990— 1995年呈现出下降趋势, 减少幅度较大。1995— 2005年产水量逐步上升。2005年以后产水量持续下降, 下降幅度在34.97%~42.46%之间。总体看来, 流域产水功能出现衰退趋势(图2)。当前流域平均降水深度为597.03 mm, 是近25 a流域降雨量最低值。从总体趋势来看, 1990— 2013年流域平均降雨量波动明显, 1990年降雨量最大, 达到1 026.09 mm, 除了2005年降雨量急剧上升, 流域近25 a降雨量呈现下降趋势。流域产水量与降雨量变化趋势一致, 降雨量为流域生态系统产水功能空间分布的重要驱动因素。

图2 1990— 2013年南四湖流域产水量与降雨量Fig. 2 Water yield and precipitation histograms in Nansi Lake Basin during 1990-2013

2.3 当前南四湖流域产水功能空间格局分析

流域生态系统产水量空间分布格局与自然地理和社会经济要素的相关性并不密切(表1图3), 其中, 产水量空间分布与人口、降水和GDP空间部分格局的相关性相对较高, 但是仅在0.208~0.275之间, 表明除上述因素之外, 还有其他影响因素起着重要作用, 如土地利用类型。总体来看, 当前南四湖流域产水空间格局具有东高西低的特点, 湖东山地、丘陵地区的新泰市、泗水县、山亭区、薛城区、邹城市、滕州市等县区产水量较大, 湖区西部平原的东明县、兰考县、定陶县、曹县、菏泽市、牡丹区等县区水量较低, 流域中部嘉祥县、金乡县、丰县、鱼台县等产水量处于中等水平(图4)。考虑到土地利用类型对流域产水格局的影响, 可以看出, 建设用地或者城区分布集中的地区产水量较高, 林地、草地等土地利用类型产水量较低。城市化的发展增加了不透水地面的面积, 改变了水量平衡, 使得降水入渗减少、洪峰流量增加[25], 因此建设用地产水量高, 而林地、草地则对地表径流具有拦蓄作用, 同时土壤入渗量增加, 延迟了流域降水量汇流时间, 削减了洪峰流量[26], 故其产水量较低。

图3 南四湖流域当前产水与各因素分布Fig. 3 The spatial patterns of current water yield and driving factors in Nansi Lake Basin

表1 当前产水量与各影响因素相关性分析 Table 1 Correlation analysis of current water yield and various influencing factors

流域生态系统产水量空间分布格局与经济发展水平、人口密度分布并不一致, 分布格局与二者的相关性分别为0.208、0.275。湖区西部产水量高的地方多是经济发达、人口密集的地区, 如砀山县, 生态系统产水量符合经济发展需求, 在该区域经济发展与产水分布相一致。然而, 新泰市、泗水县、山亭区等县区虽然产水量很高, 但是经济发展水平却不高, 人口也比较稀疏, 原因在于这几个县区在山地丘陵地区, 海拔较高。巨野县、曹县、金乡县地处平原, 产水量较低, 但是人口密度却很高。铜山县、丰县等县区降水量高, 因而产水量较高, 但是经济发展水平较低。总体来看, 南四湖流域生态系统产水功能空间分布不均匀。

图4 2013年南四湖流域各县区产水量Fig. 4 Water yield of every district in Nansi Lake Basin in 2013

2.4 流域产水量空间格局变化分析

图5可知, 1990— 1995年南四湖流域大部分地区产水量呈现下降的趋势, 中部、东部地区产水减少量在200~400 mm 之间, 西部地区产水减少量相对较小, 在0~200 mm之间, 全流域仅有部分地区产水量增加。1995— 2000年流域产水量略有增加, 基本与1995年产水量持平, 南部地区产水量增加, 增加量在0~200 mm之间, 西南部部分地区产水增加量达到200 mm以上; 但是, 北部地区产水量呈降低趋势, 减少量在0~100 mm之间。2000— 2005年流域整体产水量同样呈增加趋势, 仅有西南部部分地区产水量减少。2005— 2010年流域产水量有较大幅度降低, 东部部分地区减少量达到200 mm以上。2010— 2013年流域产水量仍然呈现下降趋势, 大部分地区产水量减少量在100 mm 左右, 仅有东南部部分地区产水量在增加。此外, 产水量峰值区域由东北部的新泰市、泗水县向偏南地区山亭区、薛城区转移, 最低值区由西部定陶县、曹县向中部地区成武县、巨野县转移。总体来看南四湖流域产水量空间格局处于动态变化之中, 并且与降水量空间分布格局一致, 可以看出, 产水量空间分布格局变化主要受降水因素的影响(图5图6)。

图5 1990— 2013南四湖流域产水变化量 (mm)Fig. 5 Variation of the water yield during 1990-2013 in Nansi Lake Basin

图6 1990— 2013年南四湖流域降水量 (mm) 空间分布格局Fig. 6 Spatial distribution pattern of precipation in Nansi Lake Basin during 1990-2013

2.5 产水量变化驱动因素分析

表2所示, 1990— 2013年南四湖流域产水量空间变化趋势与气候、社会经济和土地利用变化因素显著相关。产水量与降雨量变化趋势相同, 显著正相关。因1990年流域发生洪水灾害, 受高强度降雨影响, 产水量为近20多年来最高。而2000年山东省平均降水量比多年平均偏少10.2%, 故流域产水量下降, 属偏枯水年。2005年山东省平均降雨量超出多年平均降水量19.3%, 产水量升高, 属丰水年。2010年山东省比多年平均降水量(679.5 mm)高出2.5%, 产水量变化不大, 属于持平年份。2013年菏泽、济宁、枣庄等地降雨量比多年平均降水量偏少13.6%~23.4%, 流域属枯水年。

表2 1990— 2013年产水量与各影响因素的相关性分析 Table 2 Correlation analysis of water yield and various influencing factors during 1990-2013

社会经济活动与土地利用变化紧密联系, 1990— 2013年土地利用类型变化趋势中城市建设用地、草地、林地的面积变化与产水量变化趋势显著相关。城市建设用地与社会经济发展水平是一致的, 因此, GDP高值区域一般位于城市建设用地密集的地区, 人口密度较大, 不透水层面积较大, 容易形成地表径流, 因此, 流域产水量变化与GDP、人口呈现显著的正相关(表2)。

海拔低的低洼地区对地表径流有一定的蓄积作用, 形成径流潜力大, 因此产水量变化与海拔高度呈显著正相关(r=0.321, P< 0.01)。坡度与产水量变化的相关性与海拔相同(r=0.321, P< 0.01)。受坡度和海拔的影响, 流域东部山地、丘陵地区产水量较大, 西部平原地区产水量较低。

2.6 水源供给重要功能区划分

南四湖流域水资源区属全国水资源三级分区, 属于淮河区(一级分区)、沂沭泗河区(二级分区)。为了筛选出三级分区水源供给的重要功能区, 本研究依据流域单位面积产水量的数据分布情况, 提出了南四湖流域水源供给功能重要功能区的划分标准(表3), 将流域170个子汇流单元划分为4个级别, 分别为极重要(Ⅳ )、中等重要(Ⅲ )和比较重要(Ⅱ )以及不重要(Ⅰ )(表4图7)。流域产水分区各县区的分布情如表4。I区以耕地和建设用地为主, 虽然建设用地面积分布较大, 但由于降雨量偏低, 加上大面积农田水分耗散, 因此为不重要产水区, 面积占36.1%, 但是水资源供给量仅占22.1%。II区虽然耕地面积分布较大, 但是由于受到丰沛的降水影响, 其产水深度要比Ⅰ 区高, 其面积占流域面积的24.4%, 产水量占流域的25.8%, 是水资源的比较重要功能区域。Ⅲ 型产水区降雨量较大, 且分布有大面积的城镇建设用地, 不透水层的增加使其产水量急剧增加, 在暴雨期间易形成城市内涝以及洪峰风险, 虽然仅占流域总面积的21.2%, 但Ⅲ 型产水区产水总量在所有分区中最大, 占流域总产水量的27.9%, 既是流域洪涝灾害预防的关键区域, 也是水资源供给的中等重要功能区域。林地、草地等大部分位于Ⅳ 区, 处于山东丘陵的边缘地带。山地、丘陵产水区受到地形影响, 降雨量高, 占流域总面积的14.4%, 占流域产水总量的23.7%, 因此单位面积产水量最大, 但是由于受到地形的阻拦以及林地、草地对径流的拦截作用, 洪峰风险相对较低, 对于流域水资源的涵养具有重要作用, 是水资源的极重要功能区域。

表3 流域水源供给重要功能区划分标准 Table 3 Division criterion of water resources function zoning in the basin
表4 水资源分区中的县区分布情况 Table 4 The distribution of the county districts in water resources

图7 南四湖流域水资源功能区Fig. 7 Water resources zoning in Nansi Lake Basin

3 结论与建议

本研究基于土地利用数据以及自然地理和社会经济等数据, 运用InVEST模型研究了南四湖流域产水量空间分布及其变化特征, 得出了以下结论:

1)通过InVEST模型模拟发现, 当前南四湖流域生态系统的产水量在空间格局上呈现出由东向西递减的趋势, 东部、东北部等山区丘陵地区产水量高, 西部平原地区产水相对较低, 产水量空间分布与经济发展水平即GDP、人口密度的空间分布格局有较大差别, 南四湖流域生态系统产水功能空间分布不均匀。

2)近25 a来, 流域产水量呈现减少趋势, 且产水量峰值区域由东北部的新泰市、泗水县向偏南山亭区、薛城区转移, 最低值区由西部定陶县、曹县向中部地区成武县、巨野县转移。

3)降水、海拔和坡度等自然地理环境要素与产水量的变化呈显著正相关, 降水量的相关程度最强; 人口、GDP等社会经济数据与产水量变化也呈显著正相关, 原因在于城市化的发展, 城市建设用地等不透水地层增加, 促进了流域产水量。林地、草地的变化趋势则与流域生态系统产水呈显著负相关, 一定程度上抑制了流域的产水, 起到了水源涵养、固持水土的功能。

4)近25 a间, 尽管气候变化对流域的产水量变化有较大的影响, 但经济驱动仍然是导致流域产水量和空间格局变化的主要原因之一, 城市化进程的加快, 人口和GDP的增长对流域产水量变化产生了显著影响。地形因素也在一定程度上影响了南四湖流域产水量的变化, 流域东部山地、丘陵地区产水较高, 西部平原地区产水量较低与流域海拔高度和坡度的分布有一定的关系。

水资源供给是制定水资源利用、生态保护等政策与规划的基础。尽管本研究采用InVEST模型对流域产水量进行了定量评估, 但是由于产水量变化的复杂性, 未来在以下几方面还需要做深入研究:1)产水量与水源涵养息息相关, 通过产水量对流域水资源量进行定量评估进而缓解洪水压力是今后的关键问题; 2)土地利用变化的水文效应会显著影响产水量的变化, 因此揭示土地利用变化水文效应的时空特征, 可为协调水土资源矛盾、改善区域生态安全状况、进行生态系统管理提供科学依据。

基于本研究结果, 南四湖流域生态系统产水功能与人口、经济分布格局具有较大差异, 因此实现流域城市化、人口和经济与水资源利用的协调发展, 对于维持流域水资源的可持续发展具有重要意义。在国家和地方政府政策推动下, 合理布局流域城镇建设用地, 控制人口规模, 改变粗放的经济发展方式, 实现人口、水资源、经济和谐发展, 是流域今后发展的一个主要任务。在水资源极重要功能区要加强水源地的保护, 退耕还林, 强化水生态系统的保护与修复; 水资源中等重要功能区域控制好环境污染, 切实保护好流域水质, 完善流域水库防洪体系; 水资源比较重要功能区域, 充分利用区内降水资源, 严格地下水管理和保护, 逐步削减地下水超采量; 水资源不重要功能区域要利用区内调水解决水资源短缺问题; 湖区建立防护林缓冲带, 防患洪峰威胁; 同时建立起符合流域实际的最佳水资源利用与经济发展模式, 为流域生态保护与建设提供科学依据, 也是解决区域水资源供需矛盾最基础的工作之一。

The authors have declared that no competing interests exist.

参考文献
[1] 孙艳伟, 李加林, 马仁锋, . 于桥水库流域水源供给服务的空间分布格局[J]. 水资源与水工程学报, 2015, 26(6): 1-6.
[SUN Y W, LI J L, MA R F, et al. Style of spatial distribution of water supply service in Yuqiao watershed. Journal of Water Resources&Water Engineering, 2015, 26(6): 1-6. ] [本文引用:2]
[2] 孙小银, 周启星, 于宏兵, . 中美生态分区及其分级体系比较研究[J]. 生态学报, 2010, 30(11): 3010-3017.
[SUN X Y, ZHOU Q X, YU H B, et al. Comparative study on ecoregion and its classification systems between China and US. Acta Ecologica Sinica, 2010, 30(11): 3010-3017. ] [本文引用:1]
[3] 吕一河, 胡健, 孙飞翔, . 水源涵养与水文调节和而不同的陆地生态系统水文服务[J]. 生态学报, 2015, 35(15): 5191-5196.
[LV Y H, HU J, SUN F X, et al. Water retention and hydrological regulation: Harmony but not the same in terrestrial hydrological ecosystem services. Acta Ecologica Sinica, 2015, 35(15): 5191-5196. ] [本文引用:1]
[4] DAILY G C, MATSON P A. Ecosystem services: From theory to implementation[J]. Proceedings of National Academy Sciences of the United States of America, 2008, 105: 9455-9456. [本文引用:1]
[5] 孙淑云, 张国玉. 南四湖流域水资源特性研究[J]. 水文, 2013, 33(1): 90-93.
[SUN S Y, ZHANG G Y. Research on characteristics of water resources in Nansihu Lake Basin. Journal of China Hydrology, 2013, 33(1): 90-93. ] [本文引用:1]
[6] 王崴, 许新宜, 王红瑞, . 基于PSR与DCE综合模型的水资源短缺程度及变化趋势分析——以北京市为例[J]. 自然资源学报, 2015, 30(10): 1725-1734.
[WANG W, XU X Y, WANG H R, et al. Water shortage and its trend analysis in Beijing based on PSR and DCE Model. Journal of Natural Resources, 2015, 30(10): 1725-1734. ] [本文引用:1]
[7] 刘斌涛, 陶和平, 孔博. 云南省水资源时空分布格局及综合评价[J]. 自然资源学报, 2014, 29(3): 454-465.
[LIU B T, TAO H P, KONG B, et al. Research on the water resource security and its temporal spatial distributions of Yunnan Province, China. Journal of Natural Resources, 2014, 29(3): 454-465. ] [本文引用:1]
[8] 程怀文, 李玉文, 徐中民. 水资源短缺的社会适应能力理论及实证——以黑河流域为例[J]. 生态学报, 2011, 31(5): 1430-1439.
[CHENG H W, LI Y W, XU Z M . Social adaptive capacity for water resource scarcity in human systems and case study on its measuring. Acta Ecologica Sinica, 2011, 31(5): 1430-1439. ] [本文引用:1]
[9] TETTADO M, ACUÑA V, ENNAANAY D, et al. Impact of climate extremes on hydrological ecosystem services in a heavily humanized Mediterranean Basin[J]. Ecological Indicators, 2014, 37: 199-209. [本文引用:1]
[10] 潘韬, 吴绍洪, 戴尔阜, . 基于InVEST模型的三江源区生态系统水源供给服务时空变化[J]. 应用生态学报, 2013, 24(1): 183-189.
[PAN T, WU S H, DAI E F, et al. Spatiotemporal variation of water source supply service in Three Rivers Source Area of China based on InVEST model. Chinese Journal of Applied Ecology, 2013, 24(1): 183-189. ] [本文引用:1]
[11] 吴哲, 陈歆, 刘贝贝, . 不同土地利用/覆盖类型下海南岛产水量空间分布模拟[J]. 水资源保护, 2014, 30(3): 9-13.
[WU Z, CHEN X, LIU B B, et al. Simulation of spatial distribution of water yield of Hainan Island with different types of land use /land cover. Water Resources Protection, 2014, 30(3): 9-13. ] [本文引用:1]
[12] LI C, ZHENG H, LI S Z, et al. Impacts of conservation and human development policy across stakeholders and scales[J]. Proceedings of National Academy Sciences of the United States of America, 2015, 112(24): 7396-7401. [本文引用:1]
[13] ZHNAG C Q, LI W H, ZHANG B, et al. Water yield of Xitiaoxi River Basin based on InVEST modeling[J]. Journal of Resources and Ecology, 2012, 3(1): 50-54. [本文引用:1]
[14] 王大尚, 李屹峰, 郑华, . 密云水库上游流域生态系统服务功能空间特征及其与居民福祉的关系[J]. 生态学报, 2014, 34(1): 70-81.
[WANG D S, LI Y F, ZHENG H, et al. Ecosystem services’ spatial characteristics and their relationships with residents’ well-being in Miyun Reservoir watershed. Acta Ecologica Sinica, 2014, 34(1): 70-81. ] [本文引用:1]
[15] 包玉斌, 李婷, 柳辉, . 基于InVEST模型的陕北黄土高原水源涵养功能时空变化[J]. 地理研究, 2016, 35(4): 664-676.
[BAO Y B, LI T, LIU H, et al. Spatial and temporal changes of water conservation of Loess Plateau in northern Shaanxi Province by InVEST model. Geographical Research, 2016, 35(4): 664-676. ] [本文引用:2]
[16] SU C H, FU B J. Evolution of ecosystem services in the Chinese Loess Plateau under climatic and land use changes[J]. Global and Planetary Change, 2013, 101: 119-128. [本文引用:1]
[17] 黄博强, 黄金良, 李迅, . 基于GIS和InVEST模型的海岸带生态系统服务价值时空动态变化分析——以龙海市为例[J]. 海洋环境科学, 2015, 34(6): 916-924.
[HUANG B Q, HUANG J Q, LI X, et al. Spatiotemporal dynamics of coastal ecosystem services using GIS and InVEST: A case study in the Longhai City. Marine Environmental Science, 2015, 34(6): 916-924. ] [本文引用:2]
[18] 林世伟, 武瑞东. “三江并流”区生态系统供水服务的空间分布特征[J]. 西部林业科学, 2015, 44(3): 8-15.
[LIN S W, WU R D. Applicable prospect and method of InVEST model in water supply of ecogical system of Shiyang River Basin. Journal of West China Forestry Science, 2015, 44(3): 8-15. ] [本文引用:1]
[19] 王雅, 蒙吉军. 基于InVEST模型的黑河中游土地利用变化水文效应时空分析[J]. 北京大学学报(自然科学版), 2015, 51(6): 1157-1165.
[WANG Y, MENG J J. Spatial-Temporal features of hydrological effects caused by land use changes in middle reaches of Heihe River based on InVEST model. Acta Scientiarum Naturalium Universitatis Pekinensis, 2015, 51(6): 1-9. ] [本文引用:1]
[20] SHARP R, RICKETTS T, GUERRY A D, et al. InVEST 3. 2. 0 Beta User’s Guide[M]. Stanford: The Natural Capital Project, 2015: 114-131. [本文引用:1]
[21] ZHANG L, HICKEL K, DAWES W R. A rational function approach for estimating mean annual evapotranspiration[J]. Water Resources Research, 2004, 40(2): W02502. [本文引用:1]
[22] ZHOU W Z, LIU G H, PAN J J, et al. Distribution of available soil water capacity in China[J]. Journal of Geographical Sciences, 2005, 15(1): 5-14. [本文引用:1]
[23] 刘纪远, 张增祥, 庄大方, . 20世纪90年代中国土地利用变化的遥感时空信息研究 [M]. 北京: 科学出版社, 2005, 241-243.
[LIU J Y, ZHANG Z X, ZHUANG D F, et al. Study on Remote Sensing Spatiotemporal Information in Land Use Change. Beijing: Science Press, 2005. ] [本文引用:1]
[24] 山东省水利厅. 2010年山东省水资源公报 [EB/OL]. http: //www. sdwr. gov. cn/zwgk/szygb/201508/t20150818_194849. html, 2015-08-18/2016-04-30.
[Department of Water Resources of Shangdong Province. Shand ong Provincal Water Resources Bulletin in2010. http://www.sdwr.gov.cn/zwgk/szygb/201508/t20150818_194849.html, 2015-08-18/2016-04-30. ] [本文引用:1]
[25] 李丽娟, 姜德娟, 李九一, . 土地利用/覆被变化的水文效应研究进展[J]. 自然资源学报, 2007, 22(2): 211-224.
[LI L J, JIANG D J, L I J Y, et al. Advances in hydrological response to land use /land cover change. Journal of Natural Resources, 2007, 22(2): 211-224. ] [本文引用:1]
[26] 郑绍伟, 慕长龙, 陈祖铭, . 长江上游森林影响流域水文过程模拟分析[J]. 生态学报, 2010, 30(11): 3046-3056.
[ZHENG S W, MU C L, CHEN Z M, et al. Simulations and analysis on the effects of forest on the hydrological processes in the upper reaches of Yangtze River. Acta Ecologica Sinica, 2010, 30(11): 3046-3056. ] [本文引用:1]