吳瀚逸，熊俊峰，侯渲，林晨，許金朵，馬榮華

1.可持續(xù)發(fā)展大數(shù)據(jù)國際研究中心，北京 100094

2.北京師范大學(xué)遙感科學(xué)國家重點(diǎn)實(shí)驗(yàn)室，北京 100875

3.北京師范大學(xué)環(huán)境演變與自然災(zāi)害教育部重點(diǎn)實(shí)驗(yàn)室，北京 100875

4.中國科學(xué)院南京地理與湖泊研究所，南京 210008

5.國家科技資源共享服務(wù)平臺(tái)，國家地球系統(tǒng)科學(xué)數(shù)據(jù)中心湖泊–流域分中心，南京 210008

引 言

土壤侵蝕是土地退化的重要原因，引發(fā)的一系列生態(tài)環(huán)境問題已嚴(yán)重威脅到區(qū)域的糧食生產(chǎn)與水土保持[1-3]。東北黑土區(qū)是我國重要的糧食產(chǎn)區(qū)和生態(tài)保障區(qū)。近年來，由于氣候變化和過度墾殖，導(dǎo)致了嚴(yán)重的水土流失，威脅到東北地區(qū)糧食安全和水土資源安全[4-6]。因此對(duì)于土壤侵蝕進(jìn)行準(zhǔn)確的監(jiān)測、評(píng)估以及侵蝕量計(jì)算，有助于水土保持措施的配置與優(yōu)化。

由于精準(zhǔn)的實(shí)地測量和高精度土壤水蝕數(shù)據(jù)獲取方式在大區(qū)域尺度內(nèi)難以實(shí)施且耗費(fèi)的人力物力較大，因此利用基于遙感和GIS 開發(fā)的土壤侵蝕模型對(duì)于土壤侵蝕進(jìn)行定量模擬計(jì)算，是當(dāng)前土壤侵蝕計(jì)算以及水土流失動(dòng)態(tài)監(jiān)測的重要途徑之一。在大尺度區(qū)域的土壤水蝕研究中最為常用的模型有通用土壤侵蝕模型（Universal Soil Loss Equation，USLE）和其修訂模型RUSLE（Revised Universal Soil Loss Equation），以及經(jīng)改進(jìn)后適用于中國土壤侵蝕特征的CSLE（Chinese Soil Loss Equation）模型。USLE/RUSLE 系列侵蝕模型將復(fù)雜的土壤侵蝕過程解構(gòu)為土壤特性、地形、植被覆蓋、土地管理以及降雨和徑流等多個(gè)影響因子進(jìn)行計(jì)算，可以較好地對(duì)土壤侵蝕進(jìn)行定量描述，故本研究采用RUSLE 模型進(jìn)行水蝕量計(jì)算。如表1 所示。

表1 部分土壤水蝕數(shù)據(jù)產(chǎn)品Table 1 Partial data products of soil water erosion

以上土壤水蝕數(shù)據(jù)產(chǎn)品在大尺度土壤水蝕估算研究中通常利用單一土壤水蝕影響因子算法，對(duì)于具有不同地理分異規(guī)律的地區(qū)，缺少區(qū)域性優(yōu)化，時(shí)空分辨率不夠精細(xì)；并且，對(duì)于數(shù)據(jù)精度缺少直接驗(yàn)證。同時(shí)考慮到東北地區(qū)與全國其他區(qū)域存在地理環(huán)境差異，基于其他地理?xiàng)l件開發(fā)的模型算法在東北地區(qū)的適用性有待驗(yàn)證，因此本數(shù)據(jù)集將當(dāng)前研究中常用的模型針對(duì)東北區(qū)域做出優(yōu)選與結(jié)果的直接驗(yàn)證。

針對(duì)當(dāng)前土壤水蝕模型應(yīng)用過程中區(qū)域性優(yōu)化和實(shí)測數(shù)據(jù)驗(yàn)證缺失的短板，本研究基于多源遙感數(shù)據(jù)，解構(gòu)RUSLE 模型，對(duì)模型每個(gè)因子選取多種經(jīng)典算法進(jìn)行排列組合，結(jié)合水文站數(shù)據(jù)直接驗(yàn)證篩選最優(yōu)算法組合，模擬東北地區(qū)2001–2020 年土壤水蝕模數(shù)時(shí)空分布，建立東北地區(qū)土壤水蝕時(shí)空數(shù)據(jù)集，對(duì)實(shí)現(xiàn)東北地區(qū)農(nóng)業(yè)持續(xù)發(fā)展和土地退化監(jiān)測等可持續(xù)發(fā)展目標(biāo)有著重要意義。

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

1.1 數(shù)據(jù)生產(chǎn)流程

本文擬采用技術(shù)路線如下首先基于多源遙感數(shù)據(jù)，利用RUSLE 模型，對(duì)模型每個(gè)因子選取多種經(jīng)典算法進(jìn)行排列組合計(jì)算，結(jié)合水文站數(shù)據(jù)直接驗(yàn)證篩選最優(yōu)算法組合，得到東北地區(qū)不同三級(jí)流域分區(qū)的土壤水蝕模數(shù)估算結(jié)果。本文的技術(shù)路線如圖1 所示，包含以下處理步驟：

圖1 技術(shù)路線圖Figure 1 Technical flow chart

1.2 模型及算法

RULSE 模型是美國農(nóng)業(yè)部基于通用土壤侵蝕模型USLE 修正的一種適用區(qū)域較廣的土壤水蝕模數(shù)估算模型[12-13]。該模型綜合考慮了降雨、土壤理化性質(zhì)、坡度坡長、地表植被覆蓋以及水土保持措施共5 個(gè)因子。其基本公式如下：

其中，A 代表土壤侵蝕模數(shù)，單位為t ha-1y-1；R 是降水侵蝕因子，單位為MJ mm ha-1h-1y-1，代表降水對(duì)于土壤顆粒的沖刷侵蝕作用，作為土壤水蝕的主要?jiǎng)恿σ蜃?；K 是土壤可蝕性因子，代表土壤本身理化性質(zhì)在土壤侵蝕過程中的影響作用；LS 代表地形因子，包括坡度（S）和坡長（L）兩個(gè)因子；C 表示地表植被和作物覆蓋因子，表示在植被覆蓋與種植耕作對(duì)于土壤侵蝕的影響，該因子與植被冠層、地表植被覆蓋以及其他地表植被狀況有關(guān)；P 代表水土保持措施因子，代表著梯田、溝渠、圍欄等水土保持措施對(duì)土壤水蝕的影響，值域?yàn)?–1 之間，越接近1 表示水土保持措施越少。

1.3 數(shù)據(jù)來源

本研究使用的降水?dāng)?shù)據(jù)是來自美國航空航天局（National Aeronautics and Space Administration，NASA）和日本宇宙航空研究開發(fā)機(jī)構(gòu)（Japan Aerospace Exploration Agency，JAXA）合作開展全球衛(wèi)星降水計(jì)劃（Global Precipitation Measurement mission，GPM），其利用多傳感器多衛(wèi)星多算法結(jié)合衛(wèi)星網(wǎng)絡(luò)和雨量計(jì)反演得到更高精度的降水?dāng)?shù)據(jù)，其能夠提供全球范圍基于微波和紅外的半小時(shí)的降水?dāng)?shù)據(jù)產(chǎn)品[14]。GPM 的數(shù)據(jù)產(chǎn)品主要分為：Level-1、Level-2 和Level-3。本研究使用的是Level-3 的數(shù)據(jù)是在Level-2 的基礎(chǔ)上對(duì)固定時(shí)間和空間尺度進(jìn)行插值，數(shù)據(jù)完整性和一致性較高，包括逐半小時(shí)降水以及逐月降水?dāng)?shù)據(jù)。土壤信息數(shù)據(jù)是國際應(yīng)用系統(tǒng)分析研究所（International Institute for Applied Systems Analysis，IIASA）及聯(lián)合國糧農(nóng)組織（Food and Agriculture Organization of the United Nations，F(xiàn)AO）倡導(dǎo)下建立的世界土壤數(shù)據(jù)庫（Harmonized World Soil Database version 1.2，HWSD V1.2）[15]。其中，中國地區(qū)數(shù)據(jù)源為1995 年全國第二次土壤普查由南京土壤所提供的1:1 000 000 土壤數(shù)據(jù)，數(shù)據(jù)格式為grid，投影為WGS84，采用的土壤分類系統(tǒng)主要為FAO-90。DEM 高程數(shù)據(jù)是由美國宇航局NASA 制作由美國地質(zhì)聯(lián)邦調(diào)查局（United States Geological Survey，USGS）免費(fèi)發(fā)布的高程數(shù)據(jù)，包括北緯60°到南緯56°之間陸地表面較高質(zhì)量的海拔數(shù)據(jù)[16]。數(shù)據(jù)來源詳見表2。

表2 數(shù)據(jù)來源表Table 2 Data source table

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

1.4.1 R 因子

首先對(duì)逐半小時(shí)的GPM 降水?dāng)?shù)據(jù)進(jìn)行日合成得到逐日的降水?dāng)?shù)據(jù)，根據(jù)章文波等[17]提出的R因子計(jì)算算法，由逐日降水?dāng)?shù)據(jù)計(jì)算得到逐年的R 因子時(shí)空分布，記為R1；其次，根據(jù)周伏建和黃炎和[18]提出的降水因子計(jì)算模型，由GPM 逐月降水?dāng)?shù)據(jù)進(jìn)行計(jì)算得到逐年降水侵蝕因子，記為R2。具體計(jì)算公式參照表3。

表3 R、K、C 和LS 因子Table 3 R, K, C and LS factors

1.4.2 K 因子

基于HWSD 數(shù)據(jù)選取了土壤不同粒徑顆粒組成含量以及土壤有機(jī)碳含量，通過目前土壤水蝕研究中常用的EPIC[19]、Shrazis[20]以及Torri.D[21]模型計(jì)算表示土壤的可蝕性，得到土壤可蝕性因子空間分布，依次記為K1，K2 和K3。具體計(jì)算公式參照表3。

1.4.3 C 因子

選取MODIS（Moderate-resolution Imaging Spectroradiometer）提供的MOD13A2 植被指數(shù)數(shù)據(jù)的歸一化差值植被指數(shù)（Normalized Difference Vegetation Index，NDVI）數(shù)據(jù)，首先利用蔡崇法和丁樹文[22]提出的分段的非線性轉(zhuǎn)換模型計(jì)算得到C 因子空間分布，記為C1；其次利用Van Leeuwen and Sammons[23]提出的帶參數(shù)的指數(shù)非線性模型計(jì)算得到C 因子空間分布（其中經(jīng)驗(yàn)性參數(shù)取值為=2.5； =1），記為C2。具體計(jì)算公式參照表3。

1.4.4 LS 因子

坡度坡長因子LS 代表著地形條件對(duì)于土壤侵蝕的影響。本研究基于SRTM（Shuttle Radar Topography Mission） 90 米分辨率的數(shù)字高程模型數(shù)據(jù)，利用符素華等[24]的坡度坡長法計(jì)算出坡度坡長因子，記為LS1。具體計(jì)算公式參照表3。

1.4.5 P 因子

水土保持措施因子代表著地表對(duì)于土壤水蝕的保護(hù)措施對(duì)于土壤水蝕的影響。本研究基于文獻(xiàn)調(diào)研統(tǒng)計(jì)以及中國5 年一期的土地利用數(shù)據(jù)，學(xué)者們在東北區(qū)域常用的水土保持措施因子的取值進(jìn)行賦值計(jì)算。詳細(xì)取值如表4，參考文獻(xiàn)見數(shù)據(jù)集中的補(bǔ)充材料。

表4 P 因子賦值表Table 4 P factor value table

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

2.1 數(shù)據(jù)信息及結(jié)構(gòu)

本數(shù)據(jù)集包含2001–2020 年東北區(qū)域土壤水蝕模數(shù)數(shù)據(jù)。數(shù)據(jù)集空間分辨率為250 m，時(shí)間分辨率為1 年，數(shù)值單位為（t ha-1y-1）。數(shù)據(jù)集保存為“2001–2020 年東北區(qū)域土壤水蝕模數(shù).zip”壓縮文件，壓縮文件大小為2.71 GB。其中文件以對(duì)應(yīng)年份進(jìn)行命名，如“NE_2001”代表2001 年東北區(qū)域土壤水蝕模數(shù)，共20 個(gè)文件，背景值為-1。

2.2 數(shù)據(jù)樣本

以下給出2001–2020 年，5 年一期的土壤水蝕結(jié)果時(shí)空分布圖（圖2）?？傮w上東北地區(qū)的土壤侵蝕熱點(diǎn)為西、北、南的大、小興安嶺與長白山地區(qū)，這三個(gè)地區(qū)的海拔都較周圍高，土壤侵蝕敏感性也較強(qiáng)。其中以長白山南側(cè)土壤水蝕模數(shù)較大，主要是因?yàn)樵摰貐^(qū)處在黃海暖流北上所產(chǎn)生的水汽通道上，南側(cè)為迎風(fēng)坡，因而長白山地區(qū)的降水多。整個(gè)流域的降水量非常豐沛，導(dǎo)致該地區(qū)計(jì)算出的土壤水蝕模數(shù)偏大。山區(qū)的侵蝕狀況主要受到植被變化、凍融作用等多營力的影響，隨時(shí)間變化呈現(xiàn)出波動(dòng)的趨勢；平原地區(qū)的侵蝕狀況隨時(shí)間變化呈現(xiàn)出輕度和微度的侵蝕區(qū)域緩慢增大的趨勢，強(qiáng)烈的侵蝕呈現(xiàn)減小趨勢。

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

3.1 質(zhì)量控制

本研究將水利部發(fā)布的2001–2020 年《河流泥沙公報(bào)》[25]中東北區(qū)域主要水文監(jiān)測站測量得到的年輸沙量以及輸沙模數(shù)與模型估算的土壤侵蝕模數(shù)進(jìn)行直接驗(yàn)證，并基于驗(yàn)證結(jié)果選取最優(yōu)的因子算法組合生成最終侵蝕模數(shù)估算結(jié)果。步驟如下：1）首先根據(jù)水文站點(diǎn)分布位置與流域河段位置的相對(duì)關(guān)系以及站點(diǎn)對(duì)應(yīng)的控制流域面積對(duì)初始流域評(píng)價(jià)單元進(jìn)行合并或者遷移，建立一個(gè)站點(diǎn)對(duì)一個(gè)流域單元或一個(gè)站點(diǎn)對(duì)多個(gè)流域單元的流域評(píng)價(jià)單元。2）其次對(duì)每個(gè)評(píng)價(jià)單元建立以站點(diǎn)測量得到的年輸沙量為橫坐標(biāo)，利用土壤水蝕估算模型計(jì)算的流域總土壤水蝕量為縱坐標(biāo)建立散點(diǎn)圖，并結(jié)合時(shí)間序列變化趨勢以及公開資料記載的重大河流泥沙事件剔除異常值。3）最后將水文觀測站點(diǎn)測量得到的年輸沙量與對(duì)應(yīng)流域評(píng)價(jià)單元內(nèi)計(jì)算出的土壤水蝕模數(shù)的總和對(duì)比并計(jì)算二者在控制流域面積尺度上的評(píng)估指標(biāo)，比較不同組合方式在三級(jí)流域尺度上的驗(yàn)證指標(biāo)，選取精度最優(yōu)的組合模型生成最終結(jié)果。

3.2 質(zhì)量評(píng)估

在本研究中將實(shí)測水文站點(diǎn)數(shù)據(jù)與對(duì)應(yīng)流域估算的土壤水蝕模數(shù)進(jìn)行統(tǒng)計(jì)比較，引入相關(guān)系數(shù)（Correlation coefficient，Cov）、均方根誤差（Root Mean Square Error，RMSE）、平均絕對(duì)誤差（Mean Absolute Error，MAE）等統(tǒng)計(jì)指標(biāo)進(jìn)行精度評(píng)估。相關(guān)計(jì)算公式如下：

相關(guān)系數(shù)（Cov）：

均方根誤差（RMSE）：

平均絕對(duì)誤差（MAE）：

各因子組合排列方案如表5 所示，每種因子組合方案的相關(guān)系數(shù)，RMSE 與MAE 三個(gè)評(píng)估指標(biāo)統(tǒng)計(jì)如圖3 所示，最終得到12 種土壤水蝕估算結(jié)果及驗(yàn)證結(jié)果。由于因子計(jì)算方法的適用條件以及地理要素之間的差異，根據(jù)R、K、C 因子的不同算法計(jì)算出的水蝕結(jié)果主要會(huì)根據(jù)降雨、土壤理化性質(zhì)以及地表植被覆蓋狀況的不同，而產(chǎn)生不同的結(jié)果。由圖3 可知最佳的組合因子組合方案為A1 組合，其中使用的為R1-K1-C1-LS1-P 的因子組合，其流域估算的土壤水蝕模數(shù)值與站點(diǎn)觀測值相關(guān)性分布主要為0.4–0.5 之間，相較于其他組合方案整體稍高，RMSE 分布位于0.007–0.009 之間相較于其他組合方案整體偏低。A6，A7，A10 組合方式相關(guān)性指標(biāo)明顯偏低且RMSE 與MAE 偏高，說明在東北區(qū)域以上三種因子組合的估算模型較不適用。A2–A5、A8、A9、A11、A12 相關(guān)性分布于0.3–0.45，且分布較為離散，RMSE 與MAE 相較于A1 組合偏高，整體分布于0.007–0.01 之間。

圖3 相關(guān)系數(shù)，RMSE 與MAE 統(tǒng)計(jì)箱線圖（其中A1–A12 為12 種算法組合方式，如表5 所示）Figure 3 Statistical boxplot analysis of correlation coefficient, RMSE and MAE ( A1-A12 is a combination of 12 algorithms, as shown in table 5)

表5 模型各因子組合表Table 5 The table of combination for various model factors

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

水土流失對(duì)于我國東北地區(qū)的生態(tài)環(huán)境以及糧食種植影響較為嚴(yán)重。同時(shí)受到氣候、土壤、植被、土地利用等綜合地理環(huán)境因素的影響，東北區(qū)域黑土流失嚴(yán)重，耕層變淺，肥力下降，土壤鹽堿化加劇。水土過程模擬是治理土地退化的關(guān)鍵基礎(chǔ)環(huán)節(jié)。本研究提供的2001–2020 年長時(shí)間序列的高精度、高時(shí)空分辨率的土壤水蝕模數(shù)估算結(jié)果，對(duì)東北地區(qū)水土保持以及生態(tài)環(huán)境監(jiān)測具有重要意義。

本數(shù)據(jù)集是公開發(fā)布的東北區(qū)域2001–2020 年土壤水蝕模數(shù)估算數(shù)據(jù)，具有較高的時(shí)空分辨率和實(shí)際精度，可為東北區(qū)域的土壤侵蝕評(píng)估以及水土保持決策提供基礎(chǔ)數(shù)據(jù)，幫助決策者和生產(chǎn)者制定合理的政策和風(fēng)險(xiǎn)管理策略。同時(shí)，本數(shù)據(jù)集也具有較高的科學(xué)價(jià)值與研究意義，基于模型優(yōu)選思想為空間異質(zhì)性較高區(qū)域的土壤侵蝕估算提供了新的思路與角度。在后續(xù)的研究中，我們將嘗試使用精度更高的源數(shù)據(jù)來計(jì)算各個(gè)影響因子并使用月尺度監(jiān)測數(shù)據(jù)進(jìn)行土壤侵蝕的驗(yàn)證。

數(shù)據(jù)作者分工職責(zé)

吳瀚逸（2000—），男，福建南安人，碩士研究生，研究方向?yàn)檗r(nóng)業(yè)遙感、環(huán)境遙感。主要承擔(dān)工作：研究方案設(shè)計(jì)、算法構(gòu)建、產(chǎn)品生產(chǎn)、精度驗(yàn)證和論文撰寫。

熊俊峰（1993—），男，江西九江人，助理研究員，研究方向?yàn)榄h(huán)境遙感。主要承擔(dān)工作：研究思路與方案設(shè)計(jì)、論文撰寫。

侯渲（1998—），男，甘肅天水人，工程師，研究方向數(shù)據(jù)庫建設(shè)、數(shù)據(jù)共享。主要承擔(dān)工作：數(shù)據(jù)收集整理、數(shù)據(jù)產(chǎn)品評(píng)估。

林晨（1984—），男，江蘇南京人，研究員，研究方向自然資源遙感監(jiān)測。主要承擔(dān)工作：研究思路設(shè)計(jì)。

許金朵（1982—），女，江蘇徐州人，高級(jí)工程師，研究方向數(shù)據(jù)庫建設(shè)、數(shù)據(jù)共享。主要承擔(dān)工作：研究框架設(shè)計(jì)，數(shù)據(jù)提供。

馬榮華（1972—），男，山東臨沂人，研究員，研究方向湖泊環(huán)境遙感。主要承擔(dān)工作：研究方案設(shè)計(jì)。