基于Landsat8衛(wèi)星影像的廈門市水土流失動(dòng)態(tài)變化分析

肖惠珍，陳志霞

(廈門精圖信息技術(shù)有限公司，福建 廈門 363000)

1 研究區(qū)概況

廈門市地處福建省東南沿海九龍江入海處，背山面海，地形自西北向東南傾斜，地貌構(gòu)成以臺(tái)地、丘陵、平原、山地為主，平均海拔63.2 m；屬南亞熱帶海洋性季風(fēng)氣候區(qū)，受太平洋暖濕氣流影響，每年一般有4～5次臺(tái)風(fēng)過境，年均氣溫21 ℃，年均降水量約1 200 mm，終年溫暖濕潤(rùn)，盛行東北風(fēng)；植被種類繁多，以次生植被和人工植被為主；土壤類型以紅壤、水稻土和赤紅壤為主，不同類型土壤分布與海拔密切相關(guān)。隨著經(jīng)濟(jì)的迅速發(fā)展，人多地少、人地矛盾突出的問題愈發(fā)嚴(yán)重。在自然和人為因素的共同影響下，廈門市已成為閩南沿海地區(qū)水土流失較為嚴(yán)重的區(qū)域之一。

2 水土流失影響因子

參考許木土等[1]的研究成果，確定了廈門市區(qū)范圍內(nèi)對(duì)水土流失影響較大且便于動(dòng)態(tài)監(jiān)測(cè)的3個(gè)因子，即土地利用類型、植被覆蓋度和坡度。

(1)土地利用類型。不合理的土地利用是導(dǎo)致水土流失的主要原因之一，用土地利用轉(zhuǎn)移矩陣可以很直觀地反映土地利用變化情況。本研究采用兩個(gè)時(shí)相的土地利用類型圖進(jìn)行疊加計(jì)算，得到2014—2015、2015—2016年土地利用轉(zhuǎn)移矩陣，見表1、2。

由表1、2知，廈門市林地主要轉(zhuǎn)化為耕地，其次是建筑用地。2014—2015年有11.97%的林地轉(zhuǎn)化成了耕地，3.62%轉(zhuǎn)化為建筑用地；2015—2016年有6.15%的林地轉(zhuǎn)化為耕地，2.27%轉(zhuǎn)化為建筑用地。耕地主要轉(zhuǎn)化為林地和建筑用地，而建筑用地又主要轉(zhuǎn)化為耕地。未利用地主要轉(zhuǎn)化為建筑用地和耕地，其中城市用地范圍的未利用地主要轉(zhuǎn)化為建筑用地，非城市用地范圍的未利用地主要轉(zhuǎn)化為耕地。灘涂和水體面積變化主要受潮位變化的影響。

表1 2014—2015年廈門市土地利用轉(zhuǎn)移矩陣 km2

表2 2015—2016年廈門市土地利用轉(zhuǎn)移矩陣 km2

(2)植被覆蓋度。植被覆蓋可以削弱水滴的重力勢(shì)能，使得雨水不直接沖擊地面，從而減輕水土流失。當(dāng)達(dá)到一定郁閉度時(shí)，植被能有效保護(hù)土壤不被侵蝕。本研究基于像元二分模型的歸一化植被指數(shù)(NDVI)估算植被覆蓋度(fc)[2]。植被覆蓋度計(jì)算公式為

fc=(NDVI-NDVIsoil)/(NDVIveg-NDVIsoil)

式中：NDVI為歸一化植被指數(shù)；NDVIsoil為沒有植被覆蓋區(qū)域的NDVI值；NDVIveg為完全被植被覆蓋區(qū)域的NDVI值。

NDVI的值介于-1～1：當(dāng)取值小于零時(shí)，表示地物對(duì)可見光高反射，一般為水體、云、雪等；等于零時(shí)，表示地表為裸地，一般為巖石、裸土等；大于零時(shí)，數(shù)值越大，表示植被覆蓋度越高?？紤]到?jīng)]有實(shí)測(cè)數(shù)據(jù)，NDVIveg和NDVIsoil取一定置信度范圍內(nèi)的NDVImax和NDVImin，而這兩個(gè)值可以通過統(tǒng)計(jì)研究區(qū)NDVI值的累積概率分布來獲取，取累積概率為95%的值為NDVImax，累積概率為5%的值為NDVImin。

由表3知，廈門市植被覆蓋度以30%～45%為主，其次是>75%。調(diào)查結(jié)果顯示，植被覆蓋度>75%的地區(qū)主要是有植被覆蓋的山區(qū)，覆蓋度45%～75%的主要是耕地和海拔較低的林地，覆蓋度0～45%的主要為城區(qū)用地。2014—2016年，廈門市植被覆蓋度有所增加，尤其是高植被覆蓋區(qū)域面積增加較明顯。

表3 2014、2015、2016年廈門市植被覆蓋狀況

(3)地形坡度。坡度能表示地表單元陡緩的程度，在一定范圍內(nèi)坡度越大，徑流速度越快，水流沖刷能力越強(qiáng)，水土流失就越嚴(yán)重。從2016年廈門市的DEM中提取坡度信息可知：廈門市整體坡度較為平緩，67.80%的地區(qū)坡度在5°以下，主要是海域、建筑用地、耕地和灘涂等；坡度>5°的地區(qū)多是林地，主要分布在翔安東部、同安北部和集美西部的山區(qū)。

3 水土流失動(dòng)態(tài)監(jiān)測(cè)結(jié)果

參照《土壤侵蝕分類分級(jí)標(biāo)準(zhǔn)》(SL 190—2007)，結(jié)合研究區(qū)實(shí)際情況，整理得到本研究廈門市水土流失強(qiáng)度分級(jí)標(biāo)準(zhǔn)。將土地利用分類圖、植被覆蓋度分級(jí)圖、地形坡度分級(jí)圖進(jìn)行疊加計(jì)算[3-5]，分別得到2014、2015、2016年3個(gè)時(shí)相的水土流失強(qiáng)度分級(jí)結(jié)果，見表4—6。

從單一時(shí)相的強(qiáng)度分級(jí)結(jié)果來看，2014—2016年廈門市絕大部分地區(qū)是微度侵蝕區(qū)，極強(qiáng)烈及其以上侵蝕面積很小，引起較強(qiáng)侵蝕的地類主要是耕地和未利用地。從3個(gè)時(shí)相侵蝕強(qiáng)度分級(jí)變化情況來看，廈門市輕度及其以上侵蝕面積逐年減少，說明整體侵蝕強(qiáng)度有所降低，尤其是強(qiáng)烈及其以上侵蝕面積減少明顯，由2014年的31.70 km2減少到2016年的21.35 km2，減少了32.65%。

表4 2014年廈門市各地類水土流失強(qiáng)度分級(jí)結(jié)果

表5 2015年廈門市各地類水土流失強(qiáng)度分級(jí)結(jié)果

表6 2016年廈門市各地類水土流失強(qiáng)度分級(jí)結(jié)果

將2014—2016年廈門市的水土流失面積按行政區(qū)進(jìn)行統(tǒng)計(jì)，結(jié)果見表7。其中，同安區(qū)的水土流失情況最為嚴(yán)重，其水土流失面積占全市流失總面積的50%以上，這是由于該區(qū)地形破碎復(fù)雜，特別是境內(nèi)的多個(gè)開山采石點(diǎn)嚴(yán)重破壞地表，導(dǎo)致了嚴(yán)重的水土流失；其次是耕地面積較大的翔安區(qū)；思明區(qū)和湖里區(qū)水土流失面積最少，這是由于兩區(qū)經(jīng)濟(jì)發(fā)達(dá)，城區(qū)綠化面積大，水土保持設(shè)施較完善。2014—2016年廈門市各區(qū)水土流失狀況都得到了改善，水土流失面積有明顯減少。

表7 2014—2016年廈門市各行政區(qū)水土流失面積

4 結(jié) 語(yǔ)

基于廈門市2014—2016年三期Landsat8遙感影像數(shù)據(jù)，利用GIS空間分析功能，依據(jù)《土壤侵蝕分類分級(jí)標(biāo)準(zhǔn)》(SL 190—2007)，對(duì)廈門市的水土流失動(dòng)態(tài)變化進(jìn)行分析，得到以下結(jié)論：

(1)廈門市占地面積最大的地類是林地，占土地總面積的30%左右。從土地利用變化來看，2014—2016年主要表現(xiàn)為建筑用地持續(xù)增加，林地、耕地持續(xù)減少。

(2)廈門市區(qū)植被覆蓋度以30%～45%為主。監(jiān)測(cè)數(shù)據(jù)顯示，近年來廈門市區(qū)低植被覆蓋區(qū)面積逐年減小，而高植被覆蓋區(qū)面積明顯增加，植被覆蓋度變化呈現(xiàn)穩(wěn)中有小幅增長(zhǎng)的趨勢(shì)。

(3)對(duì)廈門市三個(gè)時(shí)相不同強(qiáng)度水土流失面積統(tǒng)計(jì)分析發(fā)現(xiàn)，廈門市區(qū)輕度及其以上侵蝕面積逐年減少，尤其是強(qiáng)烈及其以上侵蝕面積減少明顯。從各個(gè)行政區(qū)來看，水土流失最嚴(yán)重的是開山采石點(diǎn)較多的同安區(qū)，其次是耕地面積較大的翔安區(qū)，思明區(qū)和湖里區(qū)水土流失相對(duì)較輕。