数据可用性声明：本文关联实体数据集已在《全球变化数据仓储电子杂志（中英文）》出版，可获取：https://doi.org/10.3974/geodb.2020.03.18.V1.
1 前言
生态系统服务是人类从生态系统中直接或间接获得的各种惠益[1]，联合国千年生态系统服务评估（MEA）将其分为供给、调节、支持和文化服务四类[2]。生态系统的产水服务大致隶属于调节和供给服务，具体指生态系统对水的利用、过滤等作用后，水在时空及数量等方面发生变化的现象及过程[3]。作为重要的调节、供给服务，生态系统产水服务是多种生态过程和生态服务功能的基础[4]，其在拦蓄降水、调节径流、净化水质、改善水文状况、调节区域水循环中发挥着重要作用[5]。开展产水服务的量化评估是合理利用水资源，保障区域可持续发展的基础[6]。因此，目前针对产水服务的研究已经成为水文学和流域管理领域的研究热点[7, 8]。
本数据集是在国家重点研发计划的支持下完成的，主要目的是为构建国家屏障区产水服务的时间序列产品，开展屏障区生态系统服务权衡与协同的研究，进而巩固和优化区域生态服务功能，确保国家生态安全及可持续发展。通过推求栅格尺度地表最大潜在可利用水量来衡量水文调节能力，即将降水量与蒸散量作差，获得屏障区的产水服务量。国家生态屏障区1-km分辨率产水数据集的计算使用了基于水量平衡原理的InVEST模型的产水模块，并将结果以栅格地图的形式进行直观表达。产水数据是屏障区生态系统服务科学中的重要产出成果，是该区域生态环境产水量演变监测评估、生态系统服务权衡和协同及生态系统可持续发展的重要数字化资源。本文介绍该数据的详细信息、数据算法基本原理及数据结果。
2 数据集元数据简介
《国家生态屏障区产水服务量1-km栅格逐年数据集（2000‒2015）》[9]的名称、短名名称、作者、地理区域、数据年代、数据分辨率、数据格式、出版与共享服务平台、数据共享政策等信息一并列于表1。
| 条目 | 描述 |
|---|---|
| 数据集名称 | 国家生态屏障区产水服务量1-km栅格逐年数据集（2000‒2015） |
| 数据集短名 | NBZ_water_1km_2000-2015 |
| 作者信息 | 尹礼唱 AAS-4914-2020，中国科学院地理科学与自然资源研究所陆地表层格局与模拟重点实验室，yinlichang3064@163.com 王晓峰 AAS-5271-2020，长安大学土地工程学院；陕西省土地工程重点实验室，wangxf@chd.edu.cn 王怡 AAS-5036-2020，长安大学地球科学与资源学院，wangyichangan134@163.com |
| 地理区域 | 涉及的省份[10]包括：黑龙江、吉林、青海、甘肃、四川、新疆、内蒙古、河北、辽宁、西藏、宁夏、云南、广西、广东、贵州、湖南、江西、山西北方防沙带（3645N‒4506N，7550E‒12418E）青藏高原生态屏障（2940N‒3810N，8250E‒1055E）川滇-黄土高原生态屏障（2410N‒3850N，9905E‒11425E）南方丘陵山地带（2245N‒2710N，10310E‒11915E）东北森林带（4052N‒5334N，11848E‒13422E） |
| 数据年代 | 2000‒2015 |
| 时间分辨率 | 1年 |
| 空间分辨率 | 1 km |
| 数据格式 | .tif |
| 数据量 | 169 MB（压缩后） |
| 数据集组成 | 16个.tif数据, 命名规则是 ‘NBZ_water_’+空间分辨率+年份，如‘NBZ_water_1km_2000.tif’表明是国家屏障区2000年1-km空间分辨率的产水数据。年份依次是从2000年到2015年 |
| 投影坐标 | WGS_1984_Albers |
续表
| 条目 | 描述 |
|---|---|
| 基金项目 | 中华人民共和国科学技术部（2019QZKK0405）；中国科学院（XDA2002040201） |
| 出版与共享服务平台 | 全球变化科学研究数据出版系统 http://www.geodoi.ac.cn |
| 地址 | 北京市朝阳区大屯路甲11号100101，中国科学院地理科学与资源研究所 |
| 数据共享政策 | 全球变化科学研究数据出版系统的“数据”包括元数据（中英文）、通过《全球变化数据仓储电子杂志（中英文）》发表的实体数据和通过《全球变化数据学报（中英文）》发表的数据论文。其共享政策如下：（1）“数据”以最便利的方式通过互联网系统免费向全社会开放，用户免费浏览、免费下载；（2）最终用户使用“数据”需要按照引用格式在参考文献或适当的位置标注数据来源；（3）增值服务用户或以任何形式散发和传播（包括通过计算机服务器）“数据”的用户需要与《全球变化数据学报（中英文）》编辑部签署书面协议，获得许可；（4）摘取“数据”中的部分记录创作新数据的作者需要遵循10%引用原则，即从本数据集中摘取的数据记录少于新数据集总记录量的10%，同时需要对摘取的数据记录标注数据来源[11] |
| 数据和论文检索系统 | DOI，DCI，CSCD，WDS/ISC，GEOSS，China GEOSS，Crossref |
3 数据研发方法
3.1 数据源及预处理
（1）气象数据是中国地面气候资料日值数据集（V3.0），来自中国气象数据网。采用专业气象插值软件ANUSPLIN[12]批量插值全国气象要素，时空分辨率均为1 kmd–1。
（2）中国土地利用数据（2000、2005、2010和2015年）来自中国土地利用现状遥感监测数据库，空间分辨率为1 km。
（3）NDVI数据（2000‒2015年）是MOD13A2 1 km植被指数16天合成产品，下载自地理空间数据云。经过批量格式转换、年最大值合成、批量裁剪和投影转换得到研究区分辨率为1 km的年NDVI栅格数据集。
3.2 InVEST模型
本数据集基于水量平衡原理，采用InVEST模型的产水模块计算国家生态屏障区2000‒2015年的产水量，以表征区域水资源供给量。利用区域水分的输入量降水与输出量蒸散发的插值，同时结合气候、地形、植被、土壤等因素[13]，定量评估基于栅格单元的不同景观组分的产水能力，具体计算公式如下：
式中：WY(x)为某景观类型x在栅格单元的年产水量（mm）；AET(x)为单元x的年实际蒸发量（mm）；P(x)为单元x的年降雨量（mm）。PET(x)为单元x的潜在蒸散量（mm），通过彭曼公式计算[14]；AWC(x)为植物可利用含水量；w(x)为经验参数；Z为Zhang系数[15]。
4 数据结果与验证
4.1 数据集组成
数据NBZ_water_1 km_2000-2015是ArcGIS TIFF格式的20002015年国家生态屏障区的年产水数据集，空间分辨率为1 km，单位是mm，投影坐标系为WGS_1984_Albers。压缩后的总数据量共计169 MB，数据解压后在ArcGIS软件下应用此数据。
4.2 数据结果
国家生态屏障区2000‒2015年1-km分辨率产水的空间分布见图1。在2000‒2015年，研究区的平均产水量为206.63 mm，总体上呈现东南高而西北低的格局。高值区集中于川滇-黄土高原生态屏障和南方丘陵山地带，中间值在东北森林带和青藏高原生态屏障的东南部，而低值区位于青藏高原生态屏障的西北部和北方防沙带。在2000‒2015年，国家屏障区86.3%的区域产水服务增加，增幅较大区域集中在东北森林带的西北部及南方丘陵山地带的中部，而产水量降低的区域集中于川滇-黄土高原的东南部。
（注：本图基于国家测绘地理信息局标准地图服务网站下载的审图号为GS(2020)4619号的
标准地图制作，底图无修改）
4.3 数据结果验证
30个站点上的验证结果表明本文2000‒2015年平均产水量与基于综合蓄水能力法计算的多年平均产水[16]具有高度的相关性，决定系数R2达到了0.81，表明本文的产水量能够较好的表征国家屏障区的产水空间分布。
（注：站点数据为观测资料）
5 结论
为构建国家屏障区产水服务的时间序列产品，使用遥感、气象、土地利用等数据基于InVEST模型的产水模块进行了建模研究。结果表明研究区的平均产水量为206.63 mm，总体上呈现东南高而西北低的格局。与基于站点观测资料的综合蓄水能力法计算的产水结果对比表明本文的产水量能够较好的表征国家屏障区的产水空间分布。
2000‒2015年1-km空间分辨率的国家生态屏障区产水数据集展示了“两屏三带”产水量的时空分布状况，是揭示近年来屏障区生态系统服务及权衡协同关系时空演变过程中的重要输入因子。它可为研究全球变化背景下国家生态屏障区的生态系统服务功能，全面深入了解和把握国家生态屏障区安全态势，为生态系统的可持续发展提供可靠的基础数据和信息。
作者分工：尹礼唱对数据集的开发做了总体设计；尹礼唱、王晓峰采集和处理了产水数据；尹礼唱设计了模型和算法；王晓峰、王怡撰写了数据论文等。