其他数据论文 II 区论文(已发表) 版本 ZH1
下载
基于InVEST模型的北京市延庆区生态系统服务功能评价研究数据集
A dataset of ecosystem services evaluation of Yanqing
 >>
: 2021 - 07 - 15
: 2021 - 08 - 23
: 2021 - 07 - 23
: 2021 - 09 - 27
42 2 0
摘要&关键词
摘要:延庆作为北京生态涵养的重要承载主体,对其生态系统服务功能进行客观科学的定量评估,可为当地土地利用规划及生态管理提供科学指导及支撑。基于InVEST模型中碳储存、水源供给、土壤保持子模型生产函数原理,采集整理2004–2014年延庆区遥感数据、气候数据、土地数据、土壤数据、农业专业数据等19类多源数据。原始数据均来自中国地理空间数据云等专业数据中心,并经过遥感影像解译、非监督分类、空间统计分析等加工处理。本数据集数据类型丰富,部分数据具有创新性,不仅涵盖了生态环境及资源领域内的大部分关键指标数据,同时可作为基础科学数据为今后延庆区相关研究提供有益补充。
关键词:生态系统服务功能;InVEST模型;延庆
Abstract & Keywords
Abstract: Yanqing plays an important role in ecological conservation in Beijing. The objective and scientific quantitative evaluation of its ecosystem services can provide scientific guidance and support for local land use planning and ecological management. Based on the production function principle of carbon storage, water supply and soil conservation sub-model in InVEST model, 19 types of multi-source data were collected, including remote sensing data, climate data, land use data, soil data and agricultural professional data in Yanqing District from 2004 to 2014.The raw data is from China Geospatial Data Cloud and other professional data centers, and processed through remote sensing image interpretation, unsupervised classification, spatial statistical analysis, etc. The dataset contains a variety of data, some of which is innovative. It not only covers most of the key index data in the field of ecological environment and resources, but also can be used as basic science data to provide the foundation and additional support for the future research in Yanqing District.
Keywords: ecosystem services; InVEST model; Yanqing
数据库(集)基本信息简介
数据库(集)名称基于InVEST模型的北京市延庆区生态系统服务功能评价研究数据集
数据作者李敏
数据通信作者刘婷婷(liutingting@caas.cn)
数据时间范围2004–2014年
地理区域北京市延庆区,地理范围为40°16′–40°47′59N,115°44′–116°34′E
空间分辨率30 m(空间栅格数据)
数据量60.8 MB(压缩文件大小)
数据格式*.shp, *.tif, *.xlsx
数据服务系统网址http://dx.doi.org/10.12205/A0007.20210616.00.ds.0386
基金项目中国农业科学院创新工程(2021CX034)
数据库(集)组成由延庆区遥感数据、气候数据、土地数据、土壤数据、农业专业数据等19类多源数据组成。包括空间栅格数据(tiff)12类:土地利用覆被图、年降雨量、气温、潜在蒸散量、根系约束层土壤深度、植物可利用含水率PAWC、DEM数据、降雨侵蚀力因子R、土壤可蚀性因子K、碳储存量、水源供给量、土壤保持量;空间矢量数据(shp)1类:集水区shp;数据表(excel)6类:土地利用覆被变化表、各地类四大碳库碳密度值、最大根系深度、蒸散系数、植被覆盖与管理因子C、土壤保持措施因子P
Dataset Profile
TitleA dataset of ecosystem services evaluation of Yanqing based on InVEST model
Data corresponding authorLIU Tingting (liutingting@caas.cn)
Data authorsLI Min
Time range2004 – 2014
Geographical scopeYangqing Distinct, Beijing, 40°16′–40°47′59N, 115°44′-116°34′E
Spatial resolution30 m (Spatial grid data)
Data volume60.8 MB (After compressed)
Data format.shp, .tif, .xlsx
Data service system<http://dx.doi.org/10.12205/A0007.20210616.00.ds.0386>
Sources of fundingInnovation Project of Chinese Academy of Agricultural Sciences (2021CX034)
Dataset compositionThe dataset consists of 19 types of multi-source data (such as remote sensing data, climate data, land data, soil data, and agricultural data), including 12 types of spatial grid data (tiff), such as land use and cover map, annual rainfall, air temperature, potential evapotranspiration, soil depth of root restraint layer, plant available water content (PAWC), DEM data, rainfall erosivity factor R, soil erodibility factor K, carbon storage, water supply and soil conservation; 1 type of space vector data (shp): catchment shp; 6 types of data sheets (excel): land use and cover change, carbon density values of four carbon pools in different regions, maximum root depth, evapotranspiration coefficient, vegetation cover and management factor C, soil conservation measures factor P.
引 言
生态系统是人类生存发展的自然环境和必然条件,其健康与否在自然资源开发过度、生境水平退化严重的当下显得至关重要。国内外学者将人类从生态系统中获得的效益定义为生态系统服务功能,主要包括供给功能、调节功能、文化功能以及支持功能等[1]。对生态系统服务功能的客观监测与评价是生态环境科学管理与实现可持续发展的先决条件,也是近年来生态与环境科学的研究重点。
延庆作为北京生态涵养的重要承载主体,在固碳增汇、涵养水源、土壤保育等方面发挥着举足轻重的作用。为科学评估延庆区生态系统服务功能,本研究采用InVEST模型,以土地利用/覆被数据为基础,对2004年、2009年、2014年延庆区碳储存、水源供给、土壤保持这三大关键生态系统服务功能的数量变化、空间分布和时空变化情况进行定量分析评估,以期为当地土地利用规划及生态管理提供科学指导及支撑[2]
本研究采用的InVEST模型通常利用生产函数的方法来量化和评估生态系统服务,并由一系列特定功能的子模型组成。不同子模型的生产函数需要对应的生态环境参数作为模型输入,因此本研究的数据集由延庆区遥感数据、气候数据、土地数据、土壤数据、农业专业数据等多源数据组成。本数据集具有四方面特点及优势:一是可用于生态系统服务功能的定量评估,二是适用于大尺度生态系统服务功能评估,三是补充了延庆这一特定区域相关研究的数据空白,四是数据种类全面丰富,数据源权威科学。本数据集既可为研究区生态系统服务功能评估提供数据基础,同时对于其他生态环境领域的科学研究也具有一定价值。
1   数据采集和处理方法
1.1   数据集构成及分类
数据集由延庆区遥感数据、气候数据、土地数据、土壤数据、农业专业数据等19类多源数据组成,考虑到本研究涉及数据类型较多,可按照数据用途及对应的子模型对数据进行分类整理,具体分类及数据集构成如图1所示。


图1   数据集构成示意图
1.2   土地利用/覆被分类图像采集与处理
1.2.1   数据源与数据采集
土地利用/覆被分类数据是本研究的基础数据,数据集中包括2004、2009、2014年三期的延庆区土地利用/覆被分类图像。原始数据为来自中国地理空间数据云(http://www.gscloud.cn/)landsat系列数据的遥感影像,空间分辨率30 m,时间分辨率16天。根据延庆县的地理位置,将经纬度分别定义为116°E和40°N,从而确定覆盖延庆区的遥感影像数据范围。数据筛选的原则主要有3方面:(1)云量小于10%的遥感影像,减少云层覆盖对后期解译处理的影响;(2)成像时间在6–9月生长季节之间;(3)三期遥感影像月度时间尽可能保持统一。综合以上因素,最终选取Landsat 5 TM数据两幅,成像时间分别为2004-09-08和2009-09-22,Landsat 8 OLI数据一幅,成像时间为2014-09-04。三期遥感影像数据的轨道号均为123/32。
1.2.2   数据处理
土地利用/覆被分类图像是在原始遥感影像解译的基础上得到的,具体处理过程分为数据预处理和非监督分类两个阶段。预处理主要包括波段叠加、辐射定标、大气校正、影像裁剪等步骤,均利用ENVI 5.1软件操作实现。其中,波段叠加选用741波段组合,其特点是兼容中红外、近红外及可见光波段信息的优势,图面色彩丰富,层次感好,能提高解译的精度与速度。之后利用ENVI软件的Calibration Utilities工具对叠加后的假彩色图进行辐射定标,再利用Flaash模块完成大气校正过程,从而消除遥感传感器本身的误差以及大气散射引起的辐射误差。预处理结束后,利用ENVI非监督分类法对遥感影像进行解译分类,参照全国土地分类的Ⅰ级分类标准及研究区的实际情况,将延庆区划分为城镇、农田、林地、草地、水体、裸地6大土地利用类型,并最终生成2004–2014年间的三期土地利用/覆被图像。
1.3   碳储存子模型数据采集与处理
InVEST模型中碳储存子模型可以估算一段时间内某一地块碳的净储存量。由于环境中已有的碳储存大部分依靠四种基本的碳库,即地上生物量、地下生物量、枯落物和土壤有机质,因此碳储存模型的基本原理即用各地类的面积与四大碳库碳密度值进行相乘,从而得到不同地类碳储存总值,再累计得到研究区总的碳储存量。本数据集中四大碳库碳密度数据自政府间气候变化专门委员会(Intergovernmental Panel on Climate Change,IPCC),IPCC于2006年制定了《2006年国家温室气体清单指南》,其中第四卷则是专门针对确定农田、林地及其他类型土地温室气体的目录方法学[3]。该指南为估算碳储存提供了全面的碳密度数据,但使用时必须明确研究区所在的气候域和气候分区。通过查阅整理指南有关数据,最终形成本研究所需的不同土地利用类型(6类)四大碳库碳密度值数据(表),包括地上碳密度、地下碳密度、土壤碳密度、枯落物碳密度,数据单位t/hm2
1.4   水源供给子模型数据采集与处理
水源供给子模型用于模拟一定区域内的地表水源供给量,水源供给量越多则说明该生态系统水源供给服务功能越强。该模型是一个忽略地下水影响的简化水文循环模型,利用水量平衡原理,基于Budyko曲线和区域年降雨量来计算水源供给量(具体公式见《InVEST 2.6.0模型用户指南》)[4],即各栅格的降雨量减去实际蒸散量即等于该栅格水源供给量。因此,此模型中的水源供给量不仅包括地表产流,还包含了土壤含水量、枯落物持水量和冠层截留量。与模型对应的数据包括:(1)年均降雨量;(2)年均最高气温、年均最低气温、年均气温;(3)潜在蒸散量;(4)根系约束层土壤深度;(5)植物可利用含水率PAWC,以上均为空间分布tiff数据;(6)不同土地利用类型的最大根系深度和(7)蒸散系数,以上数据为Excel格式;以及(8)研究区数字高程DEM和集水区shp数据。各类数据的来源、采集及处理方法具体如下。
(1)年降雨量的原始数据来自戈达德地球科学数据和信息服务中心(GES DISC)TRMM_3B43_V6月降雨数据集,空间精度0.25°,数据格式为NetCDF,分别下载2004、2009、2014年1–12月的月降雨量NC文件。由原始数据生成最终所需年降雨量数据,需要经过如下步骤:利用ArcGIS10.2软件的Multidimension Tools将原始NetCDF导出为TIFF格式并定义投影;对12个月的月降雨数据进行叠加并计算求得年平均降雨量;此时得到的年降雨量数据由于空间分辨率较差,所以需要利用ArcGIS 10.2的Raster to Point功能将栅格数据转换为点要素shp数据,之后再利用ArcGIS 10.2 Geostatistical Analyst功能中的Kriging插值法进行空间插值,并完成重采样、裁剪、空间单位换算等操作,最终生成研究所需的空间分辨率为30 m、单位为mm的年均降雨量空间数据。
(2)年均最高气温和最低气温原始数据来自美国国家气候数据中心(National Climatic Data Center,NCDC),在空间插值的基础上计算得到年平均气温和气温温差数据。通过空间计算便可得到延庆地区三期年潜在蒸散量的空间数据。
(3)潜在蒸散量主要指土壤蒸发和植物蒸腾的总耗水量。在本研究中,该数据是通过InVEST模型选用的Modified-Hargreaves方法,利用年均气温、年度最高气温与年均最低气温的差值等数据计算生成的。
(4)根系约束层土壤深度,又可理解为土壤根系最高埋藏深度。数据来源于世界土壤数据库(Harmonized World Soil Database version 1.2,HWSD)。该数据库由联合国粮农组织(FAO)和维也纳国际应用系统研究所(IIASA)共同构建完成,可提供全球土壤数据。HWSD由栅格图像文件及属性数据库组成,每个网格单元都能与常用的土壤参数相关联(如有机碳、PH值、土壤深度等),因此本研究仅需利用研究区边界图对原始数据进行裁剪便可得到对应的根系约束层土壤深度空间分布数据。
(5)植物可利用含水率PAWC用以表征土壤中可以被植物吸收利用的水量占比,一般无法直接获取。本研究采用周文佐等定义的非线性拟合土壤PAWC估算模型进行计算,相关参数包括土壤中沙粒(SAN)、粉利(SIL)、粘粒(CLA)、有机碳(C)的含量(%)。以上数据同样来自HWSD,并通过ArcGIS空间计算最终得到植物可利用含水率PAWC空间分布数据。
(6)最大根系深度指植被根系生物量达到95%时的根系深度。对于水体、城镇、裸地等没有植被覆盖的土地利用/覆被类型,最大根系深度则定义为0。本研究中最大根系深度数据同样来自世界土壤数据库,并整理形成相关数据表(Excel)。
(7)蒸散系数在本研究中是表征不同地面水分蒸发程度的参数,反映自然蒸发和植物蒸腾的综合情况,因此蒸散系数不仅要取决于大气干燥程度、辐射条件、风力强弱,还受到土壤导水能力、植被状况(植被水分输导组织、叶片气孔数量等)等因素的影响。本研究中蒸腾系数由InVEST模型本身指定,无需单独采集处理。
(8)集水区是一个雨水汇合和集中排出的地形单位,也用于区分不同流域及其供应水源地区的界限。本研究的集水区shp数据主要是基于DEM数据进行水文分析生成的。延庆区DEM数据来自美国地质勘探局(United States Geological Survey,USGS)的ASTGTM V002数据集,下载涵盖研究区地面的DEM数据影像共两幅,序号为ASTGTM2_N40E115和ASTGTM2_N40E116。用DEM生成集水区数据则需要利用ArcGIS 10.2软件的hydrology模型,经过提取水流方向,计算汇流累积量、水流长度,提取河流网络[5],最后进行流域分割分析等若干步骤完成。
1.5   土壤保持子模型数据采集与处理
土壤保持子模型主要用于对生态系统的土壤保持量进行量化评估,模型的基本原理是利用土壤潜在侵蚀量减去土壤实际侵蚀量从而得到土壤保持量。模型数据需求包括:(1)降雨侵蚀力因子R;(2)土壤可蚀性因子K;(3)植被覆盖与管理因子C;(4)土壤保持措施因子P。其中(1)(2)为空间分布tiff数据,(3)(4)为excel文件。各类数据的来源、采集及处理方法具体如下。
(1)降雨侵蚀力因子R是表征降雨导致的土壤流失或损失发生的能力大小的指标,本研究中采用Wischmeier公式[6]、利用月平均降雨量和年平均降雨量计算因子数值。降雨数据的数据来源及采集处理方式与第1.4节水源供给子模型中相同,此处省略。
(2)土壤可蚀性因子K被定义为标准小区内单位降雨侵蚀力引起的土壤侵蚀量[7],是衡量土壤抗蚀性的指标。本研究采用EPIC模型,利用沙粒(SAN)、粉利(SIL)、粘粒(CLA)、有机碳(C)含量(%)等土壤属性数据计算土壤可蚀性因子。土壤属性数据来源及采集处理方式与第1.4节水源供给子模型中相同,此处省略。
(3)植被覆盖与管理因子C反映植被覆盖及相应的管理措施对土壤流失的影响,因此不同土地利用/覆被类型、不同的植被覆盖度会形成不同的植被覆盖与管理因子值。规定因子取值介于0–1之间。值越小,表明植被覆盖与管理措施对土壤保持作用越强;相反取值为1,则表明未采取植被覆盖与管理措施。本研究查阅整理路炳军等学者[8]完成的相关研究资料与文献,最终形成不同地类植被覆盖与管理因子C数据表(Excel)。
(4)土壤保持措施因子P是指采取水土保持措施后,土壤流失量与顺坡种植时的土壤流失量的比值。取值范围为0–1[9]。值越小,表明保持措施的作用越强。取值为1时,表示未采取任何水土保持措施;相反,取值为0,则代表无侵蚀发生。本研究中土壤保持措施因子P的取值查阅联合国粮农组织《土地管理》[10]公开目录后整理而成。
2   数据样本描述
本研究数据集包含数据类型较多,数据集构成及各类数据基本信息如表1所示。本数据集是在InVEST模型原理下对研究区进行生态系统服务功能评价研究的数据基础。按照InVEST模型中各子模型生产函数的规定,1–16号数据(除2号以外)是生产函数所需的输入数据,17–19号数据是生产函数计算的结果数据。另外,本研究主要利用InVEST模型对延庆区碳储存、水源供给、土壤保持这三大关键生态系统服务功能进行量化评估,因此也可按不同评估子模型对数据进行分类。其中:1号数据为三个子模型共同的基础输入数据之一;碳储存模型输入数据为3号,输出数据为17号;水源供给模型输入数据为4–12号,输出数据为18号;土壤保持模型输入数据为13–16号,输出数据为19号。数据集共计60.8 MB,部分数据样本示例如图2所示。
表1   数据集构成及各类数据基本信息
序号数据名称格式具体参数
1北京市延庆区土地利用覆被图tif空间分辨率30 m
数量:2004、2009、2014年三期
2北京市延庆区土地利用覆被变化表xlsx单位:km2
3不同土地利用类型四大碳库碳密度值xlsx单位:t/hm2
4北京市延庆区年均降雨量空间分布tif空间分辨率30 m,单位:mm
数量:2004、2009、2014年三期
5北京市延庆区气温空间分布
(包括年均最高气温、年均最低气温、年平均气温)
tif空间分辨率30 m,单位:℃
数量:2004、2009、2014年三期
6北京市延庆区潜在蒸散量空间分布tif空间分辨率30 m,单位:mm
数量:2004、2009、2014年三期
7北京市延庆区根系约束层土壤深度空间分布tif空间分辨率30 m,单位:mm
8不同土地利用类型最大根系深度xlsx单位:mm
9北京市延庆区植物可利用含水率PAWC空间分布tif空间分辨率30 m,单位:%
10不同土地利用类型蒸散系数xlsx无量纲
11北京市延庆区数字高程DEM数据tif空间分辨率:30 m
12北京市延庆区集水区空间数据shp无量纲
13北京市延庆区降雨侵蚀力因子R空间分布tif空间分辨率:30 m
14北京市延庆区土壤可蚀性因子K空间分布tif空间分辨率:30 m
15不同土地利用类型植被覆盖与管理因子Cxlsx取值范围:0–1
16北京市延庆区不同土地利用类型土壤保持措施因子Pxlsx取值范围:0–1
17北京市延庆区碳储存量空间分布tif空间分辨率:30 m,单位:t/hm2
18北京市延庆区水源供给量空间分布tif空间分辨率:30 m,单位:mm
19北京市延庆区土壤保持量空间分布tif空间分辨率:30 m,单位:t








图2   部分数据样本展示
3   数据质量控制和评估
为保证数据质量及InVEST模型数据格式要求,本研究所有原始数据均来自专业数据中心,如中国地理空间数据云、政府间气候变化专门委员会(Intergovernmental Panel on Climate Change,IPCC)、戈达德地球科学数据和信息服务中心(GES DISC)、美国国家气候数据中心(NCDC)、世界土壤数据库(HWSD)、美国地质勘探局(USGS)、联合国粮农组织(FAO)等,部分经查阅整理的数据也都源自领域内专业的公开手册、指南和相关科学文献等。
数据的加工及处理全部在ENVI 5.1、ArcGIS10.2软件平台上进行,对于遥感影像解译、非监督分类、空间统计分析等数据加工过程均按照试验步骤操作完成。而模型计算是在《InVEST 2.6.0模型用户指南》(InVEST 2.6.0 User’s Guide)[3]指导说明下完成,数据格式及质量均严格按照模型需求进行预处理。
4   数据价值
遥感技术、地理信息系统、大数据的迅猛发展,为基于大尺度空间的生态系统服务功能定量直观评估提供了便捷、精确的技术手段。本研究数据集整合了遥感数据、气候数据、土地数据、土壤数据、农业专业数据等多源数据,数据类型丰富,基本涵盖了生态环境及资源领域内的大部分关键指标数据。同时,本数据集也包含了部分专业研究中较难直接获得的创新性数据,如四大碳库碳密度值、根系约束层土壤深度、植物可利用含水率PAWC、最大根系深度、蒸散系数、降雨侵蚀力因子R、土壤可蚀性因子K、植被覆盖与管理因子C、土壤保持措施因子P等。另外,本数据集以北京市延庆区为特定研究区,可作为补充数据为今后相关的区域研究提供科学数据基础。
[1]
张永民, 赵士洞. 全球生态系统服务的状况与趋势[J]. 地球科学进展, 2007, 22(5): 515–520.
[2]
李敏. 基于InVEST模型的生态系统服务功能评价研究——以北京延庆为例[D]. 北京:北京林业大学,2016.
[3]
The Intergovernmental Panel on Climate Change (IPCC). 2006 IPCC Guidelines for National Greenhouse Gas Inventories, Volume 4: Agriculture, Forestry and Other Land Use [EB/OL]. [2021–06–21]. https://www.ipcc-nggip.iges.or.jp/public/2006gl/vol4.html.
[4]
Sharp R, Tallis H, Ricketts T, et al. InVEST 2.6.0 User’ s Guide[M]. The Natural Capital Project, 2012.
[5]
陶艳成, 华璀, 卢远, 等. 基于DEM的钦江流域水文特征提取研究[J]. 广西师范学院学报(自然科学版), 2012, 29(4): 60–64, 82.
[6]
Wischmeier W H, Smith D D. Predicting rainfall erosion losses, USDA Agricultural Handbook No.537[M]. U.S. Department of Agriculture. 1978.
[7]
张科利, 彭文英, 杨红丽. 中国土壤可蚀性值及其估算[J]. 土壤学报, 2007, 44(1): 7–13.路炳军, [8] 刘洪鹄, 张波, 等. 基于GIS的北京市土壤侵蚀量计算[J]. 中国水土保持科学, 2011, 9(2): 24–27.
[9]
刘青泉, 李家春, 陈力, 等. 坡面流及土壤侵蚀动力学(Ⅱ):土壤倾蚀[J]. 力学进展, 2004, 34: 493–506.
[10]
Food and Agriculture Organization of the United Nations (FAO). Land husbandry - Components and Strategy [EB/OL]. [2021–06–21]. http://www.fao.org/3/T1765E/t1765e00.htm#Contents.
数据引用格式
李敏. 基于InVEST模型的北京市延庆区生态系统服务功能评价研究数据集[DB/OL]. 国家农业科学数据中心, 2021. [2021-07-01]. DOI: 10.12205/A0007.20210616.00.ds.0386.
稿件与作者信息
论文引用格式
李敏, 刘婷婷, 樊景超, 孙伟, 曹姗姗, 满芮, 高飞. 基于InVEST模型的北京市延庆区生态系统服务功能评价研究数据集[J/OL]. 中国科学数据, 2021, 6(3). (2021-09-27). DOI: 10.11922/11-6035.nasdc.2021.0003.zh.
李敏
LI Min
主要承担工作:数据采集和数据论文撰写。
女,山西省大同市人,硕士,助理研究员,研究方向为农业信息技术及应用。
刘婷婷
LIU Tingting
主要承担工作:数据整理。
liutingting@caas.cn
女,北京人,硕士,助理研究员,研究方向为科学数据管理。
樊景超
FAN Jingchao
主要承担工作:数据整理。
男,辽宁省沈阳市人,博士,副研究员,研究方向为农业科学数据管理。
孙伟
SUN Wei
主要承担工作:数据分析。
男,山东海阳人,博士,副研究员,研究方向为农林3S技术、空间信息智能分析技术。
曹姗姗
CAO Shanshan
主要承担工作:数据分析。
女,黑龙江哈尔滨人,博士,副研究员,研究方向为农业时空信息智能分析。
满芮
MAN Rui
主要承担工作:数据整理。
女,北京人,博士,助理研究员,研究方向为科学数据管理。
高飞
GAO Fei
主要承担工作:数据整理。
河北石家庄人,博士后,研究方向为农业科学数据管理。
中国农业科学院创新工程(2021CX034)
Innovation Project of Chinese Academy of Agricultural Sciences (2021CX034)
出版历史
II区出版时间:2021年9月27日 ( 版本ZH1
参考文献列表中查看
中国科学数据·农业科学
csdata