鄒海天，李子軒，張彩云

（1.海河水利委員會海河流域水土保持監(jiān)測中心站，天津 300170；2.懷來縣水務(wù)局，河北 懷來 075400）

隨著國家深入推進長江經(jīng)濟帶發(fā)展、黃河流域生態(tài)保護和高質(zhì)量發(fā)展以及京津冀協(xié)同發(fā)展等戰(zhàn)略的實施，水土流失等生態(tài)環(huán)境問題進一步得到重視，水土流失綜合治理力度逐年加大，迫切需要水土流失監(jiān)測數(shù)據(jù)為水土流失綜合治理、規(guī)劃設(shè)計和效果評價提供數(shù)據(jù)支撐。近年來，全國采用遙感解譯、定位觀測與柵格模型計算等方法，開展水土流失動態(tài)監(jiān)測，獲取了區(qū)域的水土流失面積、強度和分布[1]?；跂鸥裣裨嬎?，能夠快速掌握區(qū)域內(nèi)水土流失及分布情況，但在同一地塊內(nèi)存在不同土壤侵蝕強度等級的柵格像元，這就導(dǎo)致無法直接確定地塊土壤侵蝕強度等級，結(jié)果難以用于水土流失綜合治理圖斑設(shè)計。本文選取河北省懷來縣外井溝小流域，采用基于柵格像元和矢量地塊2種方法分別計算獲取小流域內(nèi)水土流失狀況，對比分析2種方法計算結(jié)果差異性和適用性，為水土流失監(jiān)測成果應(yīng)用提供參考。

1 研究區(qū)概況

外井溝小流域位于河北省懷來縣南部，距縣城37 km，面積56.93 km2，處于永定河上游，屬燕山山地丘陵水源涵養(yǎng)生態(tài)維護區(qū)，以水力侵蝕為主。地勢南高北低、東高西低，最高海拔1 653 m，最低海拔810 m，相對高差843 m，地貌形態(tài)主要為山地。屬溫帶大陸性季風(fēng)氣候，四季分明，年均氣溫6℃，7月溫度最高，1月溫度最低，無霜期125 d，凍土深大于1 m。多年平均年降雨量550 mm，多集中在6—9月，多年平均徑流深85 mm。土壤類型以褐土為主，巖石以白云巖為主，兼有少量石灰?guī)r和石英砂巖。植物類型以次生類植物為主，多為森林植被以及灌叢。近年來，小流域內(nèi)實施了京津風(fēng)沙源治理工程等水土流失綜合治理，具有梯田、水平階等水土保持工程措施和人造林等水土保持生物措施。

2 材料與方法

2.1 數(shù)據(jù)及處理

（1）土地利用和水土保持措施數(shù)據(jù)。利用空間分辨率優(yōu)于2 m的遙感影像開展解譯工作，所有解譯最小地塊面積為400 m2。

（2）地形數(shù)據(jù)。采用空間分辨率為12.5 m的數(shù)字高程模型，通過坡度坡長因子計算工具[2]，提取坡度、坡長和CSLE模型中涉及的坡度因子、坡長因子數(shù)據(jù)。

（3）土壤數(shù)據(jù)。源于1∶50萬土壤類型圖。

（4）植被覆蓋度數(shù)據(jù)。采用2016—2018年每年24期MOD13Q1 NDVI數(shù)據(jù)，依據(jù)章文波[3]等提出的計算方法計算植被覆蓋度數(shù)據(jù)，確定CSLE模型中的植被覆蓋與生物措施因子。

（5）CSLE模型中涉及的降雨侵蝕力因子、土壤可蝕性因子和耕作措施因子。采用全國水土流失動態(tài)監(jiān)測項目數(shù)據(jù)。

2.2 土壤侵蝕模數(shù)計算方程

采用CSLE模型計算土壤侵蝕模數(shù)[4]。

式中：A為土壤侵蝕模數(shù)[t/（hm2·a）]；R為降雨侵蝕力因子[MJ·mm/（hm2·h·a）]；K為土壤可蝕性因子[t·hm2·h/（hm2·MJ·mm）]；L為坡長因子；S為坡度因子；B為植被覆蓋與生物措施因子；E為工程措施因子；T為耕作措施因子。

2.3 基于柵格像元的土壤侵蝕計算和統(tǒng)計

運用基礎(chǔ)數(shù)據(jù)計算各因子，經(jīng)重采樣，生成10 m×10 m分辨率柵格，運用ArcGIS軟件，對7個圖層進行乘積運算，得到每個柵格的土壤侵蝕模數(shù)，按照《土壤侵蝕分類分級標(biāo)準(zhǔn)》（SL 190-2007），判斷每個柵格的土壤侵蝕強度，統(tǒng)計不同土壤侵蝕強度的面積和比例。

2.4 基于矢量地塊的土壤侵蝕計算和統(tǒng)計

根據(jù)CSLE模型計算原理可知，計算結(jié)果受土地利用類型、水土保持措施類型、降水、土壤、地形、植被覆蓋等因素影響。為盡可能保證同一地塊內(nèi)計算因子相對唯一，在遙感解譯的土地利用和水土保持措施基礎(chǔ)上，依據(jù)土壤類型、植被覆蓋度、地形數(shù)據(jù)進一步細(xì)化確定矢量地塊。依據(jù)影響因素確定矢量地塊，獲取每個地塊的R、K、L、S、B、E和T因子值，在矢量文件屬性表中新建屬性列，運用ArcGIS軟件對7個因子進行乘積運算，得到每個地塊的土壤侵蝕模數(shù)，按照《土壤侵蝕分類分級標(biāo)準(zhǔn)》（SL 190-2007），判定每個地塊土壤侵蝕強度，統(tǒng)計不同土壤侵蝕強度的面積和比例。

3 結(jié)果與分析

3.1 基于柵格像元計算的外井溝小流域土壤侵蝕結(jié)果

基于柵格像元計算的外井溝小流域共有水土流失面積31.13 km2，占小流域面積的55.02%，其中輕度侵蝕31.05 km2，占小流域面積的54.88%；中度侵蝕0.03 km2，占小流域面積的0.05%；強烈侵蝕0.04 km2，占小流域面積的0.07%；極強烈侵蝕0.01 km2，占小流域面積的0.02%。

按土地利用類型劃分，小流域內(nèi)林地水土流失面積占總水土流失面積的66.65%，草地28.22%，耕地2.70%，建設(shè)用地、園地和交通用地分別為1.35%、1.04%和0.04%。

3.2 基于矢量地塊計算的外井溝小流域土壤侵蝕

（1）基于土壤侵蝕影響因素確定的外井溝小流域矢量地塊解譯結(jié)果。外井溝小流域解譯土地利用和水土保持措施矢量地塊912個，平均16個/km2，最大地塊為有林地，面積2.05 km2；最小地塊為灌木林地，面積429 m2。疊加土壤類別數(shù)據(jù)后，矢量地塊數(shù)量為1 117個；疊加植被覆蓋度分級數(shù)據(jù)并消除面積小于400 m2地塊后，矢量地塊數(shù)量為2 512個；疊加地形數(shù)據(jù)并消除面積小于400 m2地塊后，矢量地塊數(shù)量為7 833個。

結(jié)果顯示，植被覆蓋和地形對地塊提取數(shù)量影響較大，小流域內(nèi)提取的園地、林地和草地內(nèi)不同等級植被覆蓋地塊有1 545個，提取的不同坡面和坡度地塊有1 808個，在疊加植被覆蓋和地形因素后，矢量地塊較未疊加該因素前分別增加了2.25倍和3.12倍。不同影響因素確定矢量地塊情況，詳見表1。

表1 不同影響因素確定矢量地塊情況

（2）基于矢量地塊計算的外井溝小流域土壤侵蝕結(jié)果?；诟魍寥狼治g影響因素確定后的矢量地塊進行土壤侵蝕強度計算，結(jié)果顯示：外井溝小流域共有水土流失面積33.77 km2，占小流域面積的59.68%，其中輕度侵蝕33.71 km2，占小流域面積的59.58%；中度侵蝕0.02 km2，占小流域面積的0.04%；強烈侵蝕0.04 km2，占小流域面積的0.06%；無極強烈侵蝕。

按土地利用類型劃分，小流域內(nèi)林地水土流失面積占總水土流失面積的67.98%，草地28.94%，建設(shè)用地1.88%，耕地、園地和交通用地分別為0.57%、0.52%和0.11%。

3.3 基于柵格像元和基于矢量地塊計算方法結(jié)果分析

對比2種方法計算結(jié)果發(fā)現(xiàn)，外井溝小流域基于柵格計算的水土流失面積比基于矢量地塊計算的結(jié)果小2.64 km2，二者相對差異率為8.47%。外井溝小流域2種計算方法的土壤侵蝕結(jié)果，詳見表2。從土壤侵蝕強度等級上看，基于柵格像元計算的結(jié)果在中度及以上侵蝕等級面積大于基于矢量地塊計算的結(jié)果。從空間分布上看，基于柵格像元計算的結(jié)果與基于矢量地塊計算的結(jié)果在小流域內(nèi)的水土流失位置整體沒有明顯變化，但有2點區(qū)別較大：一是基于矢量地塊計算的結(jié)果呈空間連續(xù)狀態(tài)，基于柵格像元計算的結(jié)果多呈離散分布狀態(tài)；二是中度以上侵蝕的位置變化較大，如圖1所示。從土地利用類型上看，二者計算結(jié)果中林地、草地水土流失面積占比較大，但基于柵格計算結(jié)果的耕地水土流失面積大于基于矢量計算結(jié)果。

表2 外井溝小流域2種計算方法的土壤侵蝕結(jié)果

圖1 基于柵格像元和基于矢量地塊計算的土壤侵蝕

4 討論

4.1 土壤侵蝕因子數(shù)據(jù)對計算結(jié)果的影響

基于柵格像元和基于矢量地塊2種計算方法均采用CSLE模型。模型中R、K因子均采用全國統(tǒng)一數(shù)據(jù)，比例尺小于1∶50萬；L、S因子采用空間分辨率為12.5 m的數(shù)字高程模型數(shù)據(jù)計算獲得；B因子采用空間分辨率為250 m的MODIS數(shù)據(jù)計算獲得；E、T因子基于遙感解譯的水土保持措施和土地利用數(shù)據(jù)，依據(jù)《區(qū)域水土流失動態(tài)監(jiān)測技術(shù)規(guī)定（試行）》直接賦值。基于柵格像元計算可直接、精確地反映每個柵格像元的水土流失情況，但也受因子數(shù)據(jù)的空間位置精度、準(zhǔn)確性和時效性影響，在同一地塊上因子值有多種情況?；谑噶康貕K計算采用土地利用類型、水土保持措施類型、降水、土壤、植被覆蓋度、地形等因素確定矢量地塊，能有效地保證同一地塊上各因子值的相對唯一。此處理雖然會概化地塊內(nèi)水土流失情況，但能規(guī)避因柵格像元值不一致導(dǎo)致的同一地塊計算出多種土壤侵蝕強度的問題，實現(xiàn)土壤侵蝕計算結(jié)果分布相對集中，更好地體現(xiàn)地塊總體水土流失情況。

4.2 矢量地塊確定方法

確定矢量地塊實際上就是確定計算土壤侵蝕計算的最小單元，且保證每個單元各屬性相對唯一。確定矢量地塊的方法有多種，比如網(wǎng)格法、地塊法、疊置法等[5]。為盡可能避免計算誤差，本文直接依據(jù)土壤侵蝕模型計算因子來分割矢量地塊，直接獲取矢量地塊上模型計算因子值。外井溝小流域范圍較小，降水空間分布差異不大，矢量地塊確定時暫不考慮降水因素影響。首先，基于高分辨率衛(wèi)星遙感影像目視識別地塊土地利用和水土保持措施類型，判定其屬性。其次，利用土壤類型數(shù)據(jù)，確定屬于同一土種的地塊。再次，在園地、林地和草地等地類上，依據(jù)《土壤侵蝕分類分級標(biāo)準(zhǔn)》（SL190-2007），確定植被高、中高、中、中低和低覆蓋等級，進一步細(xì)分地塊，減少地塊內(nèi)因植被覆蓋度數(shù)據(jù)差異對水土流失計算結(jié)果產(chǎn)生影響。最后，由于地面坡度的大小直接影響水土流失情況[6]，因此，本文提取小流域內(nèi)坡面和坡度，并采用坡面和坡度分級數(shù)據(jù)，確定處于同一地形的矢量地塊，減少因坡度差異較大、地形突變等對水土流失計算結(jié)果的影響。

4.3 矢量地塊確定影響因素

矢量地塊確定方法直接影響計算效率和結(jié)果準(zhǔn)確性。在數(shù)據(jù)處理過程中發(fā)現(xiàn)，土地利用、水土保持措施、土壤類型和降水較容易獲取矢量地塊計算所需的數(shù)據(jù)。但植被覆蓋度和地形數(shù)據(jù)因存在像元值空間上的連續(xù)性，在分級后植被覆蓋度、坡度和坡面數(shù)據(jù)進一步細(xì)化矢量地塊時，疊加后的矢量地塊特別破碎，較難實現(xiàn)地塊單元化。本文應(yīng)用了2種方法優(yōu)化植被覆蓋度和地形數(shù)據(jù)：①將分級后植被覆蓋度、坡度和坡面數(shù)據(jù)由柵格數(shù)據(jù)轉(zhuǎn)換為矢量數(shù)據(jù)，再與之前確定的矢量地塊疊加計算，軟件自動優(yōu)化疊加后圖斑邊界和消除碎屑多邊形；②將前序獲取的矢量地塊邊界疊加分級后植被覆蓋度、坡度和坡面柵格數(shù)據(jù)，人工提取植被覆蓋度、坡度和坡面對矢量地塊的分割線，逐個確定新的矢量地塊邊界。這2種方法前者用時較短，但需處理大量碎屑多邊形等問題；后者精準(zhǔn)度較高，但用時較長，在進行大尺度范圍的土壤侵蝕計算時不易實現(xiàn)。

5 結(jié)論

（1）基于柵格像元和基于矢量地塊計算均可實現(xiàn)小流域土壤侵蝕評價且精度相對一致，整體空間分布大致相同，結(jié)果差異為8.47%。產(chǎn)生差異原因與二者的最小計算單元、土壤侵蝕各因子計算方法等因素有關(guān)?；跂鸥裣裨嬎銜r每個計算單元都是10 m×10 m的柵格，但呈現(xiàn)在地塊上的計算結(jié)果存在多樣性；基于矢量地塊計算時計算單元存在不確定性，依據(jù)計算單元所處位置的土地利用、水土保持措施、降水、植被覆蓋度和地形綜合情況而定，但呈現(xiàn)在地塊上的計算結(jié)果是唯一的。

（2）在外井溝小流域，基于柵格像元和基于矢量地塊計算結(jié)果在高等級侵蝕強度、耕地上水土流失面積差異大，主要原因是在確定矢量地塊時對模型計算因子進行均值處理，從而對地塊上土壤侵蝕情況刻畫不如基于柵格像元計算結(jié)果細(xì)致，類別面積越小差別幅度越大。

（3）從計算過程和結(jié)果來看，基于柵格像元計算更快捷、簡單，可快速掌握區(qū)域水土流失及分布情況，但因模型計算因子數(shù)據(jù)精度不一，計算結(jié)果空間上呈現(xiàn)插花狀態(tài)，連續(xù)性不強；基于矢量地塊計算結(jié)果空間分布集中連片，水土流失地塊邊界清晰，有助于實現(xiàn)土壤侵蝕地塊單元化管理，但受計算數(shù)據(jù)精度、質(zhì)量和計算單元劃分影響，數(shù)據(jù)處理分析工作量較大。

（4）在今后的水土流失監(jiān)測工作中，需進一步優(yōu)化模型計算因子，提高數(shù)據(jù)的精度和質(zhì)量，確保不同數(shù)據(jù)層在空間上能有效匹配，使計算結(jié)果既有數(shù)值計算上的準(zhǔn)確性，又有空間分布上的合理性，從而保證水土流失監(jiān)測結(jié)果在小流域水土流失綜合治理規(guī)劃設(shè)計和效益評價中得到有效利用。