李敏，劉婷婷,3*，樊景超,3，孫偉,3，曹珊珊,3，滿芮,3*，高飛

1.中國農(nóng)業(yè)科學(xué)院農(nóng)業(yè)信息研究所，北京 100081

2.農(nóng)業(yè)農(nóng)村部農(nóng)業(yè)大數(shù)據(jù)重點實驗室，北京 100081

3.國家農(nóng)業(yè)科學(xué)數(shù)據(jù)中心，北京 100081

4.中國農(nóng)業(yè)科學(xué)院作物科學(xué)研究所，北京 100081

引 言

生態(tài)系統(tǒng)是人類生存發(fā)展的自然環(huán)境和必然條件，其健康與否在自然資源開發(fā)過度、生境水平退化嚴(yán)重的當(dāng)下顯得至關(guān)重要。國內(nèi)外學(xué)者將人類從生態(tài)系統(tǒng)中獲得的效益定義為生態(tài)系統(tǒng)服務(wù)功能，主要包括供給功能、調(diào)節(jié)功能、文化功能以及支持功能等[1]。對生態(tài)系統(tǒng)服務(wù)功能的客觀監(jiān)測與評價是生態(tài)環(huán)境科學(xué)管理與實現(xiàn)可持續(xù)發(fā)展的先決條件，也是近年來生態(tài)與環(huán)境科學(xué)的研究重點。

延慶作為北京生態(tài)涵養(yǎng)的重要承載主體，在固碳增匯、涵養(yǎng)水源、土壤保育等方面發(fā)揮著舉足輕重的作用。為科學(xué)評估延慶區(qū)生態(tài)系統(tǒng)服務(wù)功能，本研究采用InVEST 模型，以土地利用/覆被數(shù)據(jù)為基礎(chǔ)，對2004 年、2009 年、2014 年延慶區(qū)碳儲存、水源供給、土壤保持這三大關(guān)鍵生態(tài)系統(tǒng)服務(wù)功能的數(shù)量變化、空間分布和時空變化情況進(jìn)行定量分析評估，以期為當(dāng)?shù)赝恋乩靡?guī)劃及生態(tài)管理提供科學(xué)指導(dǎo)及支撐[2]。

本研究采用的InVEST 模型通常利用生產(chǎn)函數(shù)的方法來量化和評估生態(tài)系統(tǒng)服務(wù)，并由一系列特定功能的子模型組成。不同子模型的生產(chǎn)函數(shù)需要對應(yīng)的生態(tài)環(huán)境參數(shù)作為模型輸入，因此本研究的數(shù)據(jù)集由延慶區(qū)遙感數(shù)據(jù)、氣候數(shù)據(jù)、土地數(shù)據(jù)、土壤數(shù)據(jù)、農(nóng)業(yè)專業(yè)數(shù)據(jù)等多源數(shù)據(jù)組成。本數(shù)據(jù)集具有四方面特點及優(yōu)勢：一是可用于生態(tài)系統(tǒng)服務(wù)功能的定量評估，二是適用于大尺度生態(tài)系統(tǒng)服務(wù)功能評估，三是補充了延慶這一特定區(qū)域相關(guān)研究的數(shù)據(jù)空白，四是數(shù)據(jù)種類全面豐富，數(shù)據(jù)源權(quán)威科學(xué)。本數(shù)據(jù)集既可為研究區(qū)生態(tài)系統(tǒng)服務(wù)功能評估提供數(shù)據(jù)基礎(chǔ)，同時對于其他生態(tài)環(huán)境領(lǐng)域的科學(xué)研究也具有一定價值。

1 數(shù)據(jù)采集和處理方法

1.1 數(shù)據(jù)集構(gòu)成及分類

數(shù)據(jù)集由延慶區(qū)遙感數(shù)據(jù)、氣候數(shù)據(jù)、土地數(shù)據(jù)、土壤數(shù)據(jù)、農(nóng)業(yè)專業(yè)數(shù)據(jù)等19 類多源數(shù)據(jù)組成，考慮到本研究涉及數(shù)據(jù)類型較多，可按照數(shù)據(jù)用途及對應(yīng)的子模型對數(shù)據(jù)進(jìn)行分類整理，具體分類及數(shù)據(jù)集構(gòu)成如圖1 所示。

圖1 數(shù)據(jù)集構(gòu)成示意圖

1.2 土地利用/覆被分類圖像采集與處理

1.2.1 數(shù)據(jù)源與數(shù)據(jù)采集

土地利用/覆被分類數(shù)據(jù)是本研究的基礎(chǔ)數(shù)據(jù)，數(shù)據(jù)集中包括2004、2009、2014 年三期的延慶區(qū)土地利用/覆被分類圖像。原始數(shù)據(jù)為來自中國地理空間數(shù)據(jù)云（http://www.gscloud.cn/）landsat 系列數(shù)據(jù)的遙感影像，空間分辨率30 m，時間分辨率16 天。根據(jù)延慶縣的地理位置，將經(jīng)緯度分別定義為116°E 和40°N，從而確定覆蓋延慶區(qū)的遙感影像數(shù)據(jù)范圍。數(shù)據(jù)篩選的原則主要有3 方面：（1）云量小于10%的遙感影像，減少云層覆蓋對后期解譯處理的影響；（2）成像時間在6–9 月生長季節(jié)之間；（3）三期遙感影像月度時間盡可能保持統(tǒng)一。綜合以上因素，最終選取Landsat 5 TM 數(shù)據(jù)兩幅，成像時間分別為2004-09-08 和2009-09-22，Landsat 8 OLI 數(shù)據(jù)一幅，成像時間為2014-09-04。三期遙感影像數(shù)據(jù)的軌道號均為123/32。

1.2.2 數(shù)據(jù)處理

土地利用/覆被分類圖像是在原始遙感影像解譯的基礎(chǔ)上得到的，具體處理過程分為數(shù)據(jù)預(yù)處理和非監(jiān)督分類兩個階段。預(yù)處理主要包括波段疊加、輻射定標(biāo)、大氣校正、影像裁剪等步驟，均利用ENVI 5.1 軟件操作實現(xiàn)。其中，波段疊加選用741 波段組合，其特點是兼容中紅外、近紅外及可見光波段信息的優(yōu)勢，圖面色彩豐富，層次感好，能提高解譯的精度與速度。之后利用ENVI 軟件的Calibration Utilities 工具對疊加后的假彩色圖進(jìn)行輻射定標(biāo)，再利用Flaash 模塊完成大氣校正過程，從而消除遙感傳感器本身的誤差以及大氣散射引起的輻射誤差。預(yù)處理結(jié)束后，利用ENVI 非監(jiān)督分類法對遙感影像進(jìn)行解譯分類，參照全國土地分類的Ⅰ級分類標(biāo)準(zhǔn)及研究區(qū)的實際情況，將延慶區(qū)劃分為城鎮(zhèn)、農(nóng)田、林地、草地、水體、裸地6 大土地利用類型，并最終生成2004–2014 年間的三期土地利用/覆被圖像。

1.3 碳儲存子模型數(shù)據(jù)采集與處理

InVEST 模型中碳儲存子模型可以估算一段時間內(nèi)某一地塊碳的凈儲存量。由于環(huán)境中已有的碳儲存大部分依靠四種基本的碳庫，即地上生物量、地下生物量、枯落物和土壤有機質(zhì)，因此碳儲存模型的基本原理即用各地類的面積與四大碳庫碳密度值進(jìn)行相乘，從而得到不同地類碳儲存總值，再累計得到研究區(qū)總的碳儲存量。本數(shù)據(jù)集中四大碳庫碳密度數(shù)據(jù)自政府間氣候變化專門委員會（Intergovernmental Panel on Climate Change，IPCC），IPCC 于2006 年制定了《2006 年國家溫室氣體清單指南》，其中第四卷則是專門針對確定農(nóng)田、林地及其他類型土地溫室氣體的目錄方法學(xué)[3]。該指南為估算碳儲存提供了全面的碳密度數(shù)據(jù)，但使用時必須明確研究區(qū)所在的氣候域和氣候分區(qū)。通過查閱整理指南有關(guān)數(shù)據(jù)，最終形成本研究所需的不同土地利用類型（6 類）四大碳庫碳密度值數(shù)據(jù)（表），包括地上碳密度、地下碳密度、土壤碳密度、枯落物碳密度，數(shù)據(jù)單位t/hm2。

1.4 水源供給子模型數(shù)據(jù)采集與處理

水源供給子模型用于模擬一定區(qū)域內(nèi)的地表水源供給量，水源供給量越多則說明該生態(tài)系統(tǒng)水源供給服務(wù)功能越強。該模型是一個忽略地下水影響的簡化水文循環(huán)模型，利用水量平衡原理，基于Budyko 曲線和區(qū)域年降雨量來計算水源供給量（具體公式見《InVEST 2.6.0 模型用戶指南》）[4]，即各柵格的降雨量減去實際蒸散量即等于該柵格水源供給量。因此，此模型中的水源供給量不僅包括地表產(chǎn)流，還包含了土壤含水量、枯落物持水量和冠層截留量。與模型對應(yīng)的數(shù)據(jù)包括：（1）年均降雨量；（2）年均最高氣溫、年均最低氣溫、年均氣溫；（3）潛在蒸散量；（4）根系約束層土壤深度；（5）植物可利用含水率PAWC，以上均為空間分布tiff 數(shù)據(jù)；（6）不同土地利用類型的最大根系深度和（7）蒸散系數(shù)，以上數(shù)據(jù)為Excel 格式；以及（8）研究區(qū)數(shù)字高程DEM 和集水區(qū)shp數(shù)據(jù)。各類數(shù)據(jù)的來源、采集及處理方法具體如下。

（1）年降雨量的原始數(shù)據(jù)來自戈達(dá)德地球科學(xué)數(shù)據(jù)和信息服務(wù)中心（GES DISC）TRMM_3B43_V6 月降雨數(shù)據(jù)集，空間精度0.25°，數(shù)據(jù)格式為NetCDF，分別下載2004、2009、2014 年1–12 月的月降雨量NC 文件。由原始數(shù)據(jù)生成最終所需年降雨量數(shù)據(jù)，需要經(jīng)過如下步驟：利用ArcGIS10.2 軟件的Multidimension Tools 將原始NetCDF 導(dǎo)出為TIFF 格式并定義投影；對12個月的月降雨數(shù)據(jù)進(jìn)行疊加并計算求得年平均降雨量；此時得到的年降雨量數(shù)據(jù)由于空間分辨率較差，所以需要利用ArcGIS 10.2 的Raster to Point 功能將柵格數(shù)據(jù)轉(zhuǎn)換為點要素shp 數(shù)據(jù)，之后再利用ArcGIS 10.2 Geostatistical Analyst 功能中的Kriging 插值法進(jìn)行空間插值，并完成重采樣、裁剪、空間單位換算等操作，最終生成研究所需的空間分辨率為30 m、單位為mm 的年均降雨量空間數(shù)據(jù)。

（2）年均最高氣溫和最低氣溫原始數(shù)據(jù)來自美國國家氣候數(shù)據(jù)中心（National Climatic Data Center，NCDC），在空間插值的基礎(chǔ)上計算得到年平均氣溫和氣溫溫差數(shù)據(jù)。通過空間計算便可得到延慶地區(qū)三期年潛在蒸散量的空間數(shù)據(jù)。

（3）潛在蒸散量主要指土壤蒸發(fā)和植物蒸騰的總耗水量。在本研究中，該數(shù)據(jù)是通過InVEST模型選用的Modified-Hargreaves 方法，利用年均氣溫、年度最高氣溫與年均最低氣溫的差值等數(shù)據(jù)計算生成的。

（4）根系約束層土壤深度，又可理解為土壤根系最高埋藏深度。數(shù)據(jù)來源于世界土壤數(shù)據(jù)庫（Harmonized World Soil Database version 1.2，HWSD）。該數(shù)據(jù)庫由聯(lián)合國糧農(nóng)組織（FAO）和維也納國際應(yīng)用系統(tǒng)研究所（IIASA）共同構(gòu)建完成，可提供全球土壤數(shù)據(jù)。HWSD 由柵格圖像文件及屬性數(shù)據(jù)庫組成，每個網(wǎng)格單元都能與常用的土壤參數(shù)相關(guān)聯(lián)（如有機碳、PH 值、土壤深度等），因此本研究僅需利用研究區(qū)邊界圖對原始數(shù)據(jù)進(jìn)行裁剪便可得到對應(yīng)的根系約束層土壤深度空間分布數(shù)據(jù)。

（5）植物可利用含水率PAWC用以表征土壤中可以被植物吸收利用的水量占比，一般無法直接獲取。本研究采用周文佐等定義的非線性擬合土壤PAWC估算模型進(jìn)行計算，相關(guān)參數(shù)包括土壤中沙粒（SAN）、粉利（SIL）、粘粒（CLA）、有機碳（C）的含量（%）。以上數(shù)據(jù)同樣來自HWSD，并通過ArcGIS 空間計算最終得到植物可利用含水率PAWC空間分布數(shù)據(jù)。

（6）最大根系深度指植被根系生物量達(dá)到95%時的根系深度。對于水體、城鎮(zhèn)、裸地等沒有植被覆蓋的土地利用/覆被類型，最大根系深度則定義為0。本研究中最大根系深度數(shù)據(jù)同樣來自世界土壤數(shù)據(jù)庫，并整理形成相關(guān)數(shù)據(jù)表（Excel）。

（7）蒸散系數(shù)在本研究中是表征不同地面水分蒸發(fā)程度的參數(shù)，反映自然蒸發(fā)和植物蒸騰的綜合情況，因此蒸散系數(shù)不僅要取決于大氣干燥程度、輻射條件、風(fēng)力強弱，還受到土壤導(dǎo)水能力、植被狀況（植被水分輸導(dǎo)組織、葉片氣孔數(shù)量等）等因素的影響。本研究中蒸騰系數(shù)由InVEST 模型本身指定，無需單獨采集處理。

（8）集水區(qū)是一個雨水匯合和集中排出的地形單位，也用于區(qū)分不同流域及其供應(yīng)水源地區(qū)的界限。本研究的集水區(qū)shp 數(shù)據(jù)主要是基于DEM 數(shù)據(jù)進(jìn)行水文分析生成的。延慶區(qū)DEM 數(shù)據(jù)來自美國地質(zhì)勘探局（United States Geological Survey，USGS）的ASTGTM V002 數(shù)據(jù)集，下載涵蓋研究區(qū)地面的DEM 數(shù)據(jù)影像共兩幅，序號為ASTGTM2_N40E115 和ASTGTM2_N40E116。用DEM 生成集水區(qū)數(shù)據(jù)則需要利用ArcGIS 10.2 軟件的hydrology 模型，經(jīng)過提取水流方向，計算匯流累積量、水流長度，提取河流網(wǎng)絡(luò)[5]，最后進(jìn)行流域分割分析等若干步驟完成。

1.5 土壤保持子模型數(shù)據(jù)采集與處理

土壤保持子模型主要用于對生態(tài)系統(tǒng)的土壤保持量進(jìn)行量化評估，模型的基本原理是利用土壤潛在侵蝕量減去土壤實際侵蝕量從而得到土壤保持量。模型數(shù)據(jù)需求包括：（1）降雨侵蝕力因子R；（2）土壤可蝕性因子K；（3）植被覆蓋與管理因子C；（4）土壤保持措施因子P。其中（1）（2）為空間分布tiff 數(shù)據(jù)，（3）（4）為excel 文件。各類數(shù)據(jù)的來源、采集及處理方法具體如下。

（1）降雨侵蝕力因子R是表征降雨導(dǎo)致的土壤流失或損失發(fā)生的能力大小的指標(biāo)，本研究中采用Wischmeier 公式[6]、利用月平均降雨量和年平均降雨量計算因子數(shù)值。降雨數(shù)據(jù)的數(shù)據(jù)來源及采集處理方式與第1.4 節(jié)水源供給子模型中相同，此處省略。

（2）土壤可蝕性因子K被定義為標(biāo)準(zhǔn)小區(qū)內(nèi)單位降雨侵蝕力引起的土壤侵蝕量[7]，是衡量土壤抗蝕性的指標(biāo)。本研究采用EPIC 模型，利用沙粒（SAN）、粉利（SIL）、粘粒（CLA）、有機碳（C）含量（%）等土壤屬性數(shù)據(jù)計算土壤可蝕性因子。土壤屬性數(shù)據(jù)來源及采集處理方式與第1.4 節(jié)水源供給子模型中相同，此處省略。

（3）植被覆蓋與管理因子C反映植被覆蓋及相應(yīng)的管理措施對土壤流失的影響，因此不同土地利用/覆被類型、不同的植被覆蓋度會形成不同的植被覆蓋與管理因子值。規(guī)定因子取值介于0–1 之間。值越小，表明植被覆蓋與管理措施對土壤保持作用越強；相反取值為1，則表明未采取植被覆蓋與管理措施。本研究查閱整理路炳軍等學(xué)者[8]完成的相關(guān)研究資料與文獻(xiàn)，最終形成不同地類植被覆蓋與管理因子C數(shù)據(jù)表（Excel）。

（4）土壤保持措施因子P是指采取水土保持措施后，土壤流失量與順坡種植時的土壤流失量的比值。取值范圍為0–1[9]。值越小，表明保持措施的作用越強。取值為1 時，表示未采取任何水土保持措施；相反，取值為0，則代表無侵蝕發(fā)生。本研究中土壤保持措施因子P的取值查閱聯(lián)合國糧農(nóng)組織《土地管理》[10]公開目錄后整理而成。

2 數(shù)據(jù)樣本描述

本研究數(shù)據(jù)集包含數(shù)據(jù)類型較多，數(shù)據(jù)集構(gòu)成及各類數(shù)據(jù)基本信息如表1 所示。本數(shù)據(jù)集是在InVEST 模型原理下對研究區(qū)進(jìn)行生態(tài)系統(tǒng)服務(wù)功能評價研究的數(shù)據(jù)基礎(chǔ)。按照InVEST 模型中各子模型生產(chǎn)函數(shù)的規(guī)定，1–16 號數(shù)據(jù)（除2 號以外）是生產(chǎn)函數(shù)所需的輸入數(shù)據(jù)，17–19 號數(shù)據(jù)是生產(chǎn)函數(shù)計算的結(jié)果數(shù)據(jù)。另外，本研究主要利用InVEST 模型對延慶區(qū)碳儲存、水源供給、土壤保持這三大關(guān)鍵生態(tài)系統(tǒng)服務(wù)功能進(jìn)行量化評估，因此也可按不同評估子模型對數(shù)據(jù)進(jìn)行分類。其中：1號數(shù)據(jù)為三個子模型共同的基礎(chǔ)輸入數(shù)據(jù)之一；碳儲存模型輸入數(shù)據(jù)為3 號，輸出數(shù)據(jù)為17 號；水源供給模型輸入數(shù)據(jù)為4–12 號，輸出數(shù)據(jù)為18 號；土壤保持模型輸入數(shù)據(jù)為13–16 號，輸出數(shù)據(jù)為19 號。數(shù)據(jù)集共計60.8 MB，部分?jǐn)?shù)據(jù)樣本示例如圖2 所示。

圖2 部分?jǐn)?shù)據(jù)樣本展示

表1 數(shù)據(jù)集構(gòu)成及各類數(shù)據(jù)基本信息

3 數(shù)據(jù)質(zhì)量控制和評估

為保證數(shù)據(jù)質(zhì)量及InVEST 模型數(shù)據(jù)格式要求，本研究所有原始數(shù)據(jù)均來自專業(yè)數(shù)據(jù)中心，如中國地理空間數(shù)據(jù)云、政府間氣候變化專門委員會（Intergovernmental Panel on Climate Change，IPCC）、戈達(dá)德地球科學(xué)數(shù)據(jù)和信息服務(wù)中心（GES DISC）、美國國家氣候數(shù)據(jù)中心（NCDC）、世界土壤數(shù)據(jù)庫（HWSD）、美國地質(zhì)勘探局（USGS）、聯(lián)合國糧農(nóng)組織（FAO）等，部分經(jīng)查閱整理的數(shù)據(jù)也都源自領(lǐng)域內(nèi)專業(yè)的公開手冊、指南和相關(guān)科學(xué)文獻(xiàn)等。

數(shù)據(jù)的加工及處理全部在ENVI 5.1、ArcGIS10.2 軟件平臺上進(jìn)行，對于遙感影像解譯、非監(jiān)督分類、空間統(tǒng)計分析等數(shù)據(jù)加工過程均按照試驗步驟操作完成。而模型計算是在《InVEST 2.6.0 模型用戶指南》（InVEST 2.6.0 User’s Guide）[3]指導(dǎo)說明下完成，數(shù)據(jù)格式及質(zhì)量均嚴(yán)格按照模型需求進(jìn)行預(yù)處理。

4 數(shù)據(jù)價值

遙感技術(shù)、地理信息系統(tǒng)、大數(shù)據(jù)的迅猛發(fā)展，為基于大尺度空間的生態(tài)系統(tǒng)服務(wù)功能定量直觀評估提供了便捷、精確的技術(shù)手段。本研究數(shù)據(jù)集整合了遙感數(shù)據(jù)、氣候數(shù)據(jù)、土地數(shù)據(jù)、土壤數(shù)據(jù)、農(nóng)業(yè)專業(yè)數(shù)據(jù)等多源數(shù)據(jù)，數(shù)據(jù)類型豐富，基本涵蓋了生態(tài)環(huán)境及資源領(lǐng)域內(nèi)的大部分關(guān)鍵指標(biāo)數(shù)據(jù)。同時，本數(shù)據(jù)集也包含了部分專業(yè)研究中較難直接獲得的創(chuàng)新性數(shù)據(jù)，如四大碳庫碳密度值、根系約束層土壤深度、植物可利用含水率PAWC、最大根系深度、蒸散系數(shù)、降雨侵蝕力因子R、土壤可蝕性因子K、植被覆蓋與管理因子C、土壤保持措施因子P 等。另外，本數(shù)據(jù)集以北京市延慶區(qū)為特定研究區(qū)，可作為補充數(shù)據(jù)為今后相關(guān)的區(qū)域研究提供科學(xué)數(shù)據(jù)基礎(chǔ)。