2021地理设计组二等奖:基于GIS的东江源区土壤侵蚀及其影响因素空间分析
一、作品背景
水土保持情況普查對我國具有重要意義。我國目前是世界上水土流失最嚴(yán)重的國家之一,水土流失面積極其廣且量大;嚴(yán)重的水土流失問題是我國生態(tài)環(huán)境問題的重要板塊,若是持續(xù)惡化,將會嚴(yán)重影響我國的生態(tài)安全、飲水安全、防洪安全和糧食安全等,嚴(yán)重阻礙國家經(jīng)濟(jì)社會發(fā)展。以東江源區(qū)為研究對象,基于遙感數(shù)字影像和GIS方法提取影響土壤侵蝕的因子,建立CSLE模型,并利用ARCGIS10.2、GeoScene Pro2.1軟件對各評價(jià)因子指標(biāo)進(jìn)行分級,利用模型確定土壤侵蝕的強(qiáng)度,通過定量分析和分級情況分析不同的土壤侵蝕影響因子對于土壤侵蝕強(qiáng)度的貢獻(xiàn)量。
1.研究區(qū)概況
江西三百山是廣東東江的源頭,簡稱”東江源區(qū)“。東江源區(qū)位于江西省贛州市東南端,包括尋烏、安遠(yuǎn)、定南3縣。地理位置:24°20′—25°13′N,114°47′—115°53′E。流域面積3850.56km2,地形地貌以山地、丘陵為主;地勢北高南抵,東西鈴谷相間,被形象地稱之為”八山半水一分田,半分道路與莊園“。該區(qū)屬于典型的亞熱帶濕潤季風(fēng)氣候,熱量豐富,雨量充沛,年平均氣溫18.8℃,年降雨量1526~1700mm。
2.設(shè)計(jì)思想
圖1 設(shè)計(jì)思想構(gòu)圖
二、模型構(gòu)建與分析
1.數(shù)據(jù)來源及預(yù)處理
1.1 數(shù)據(jù)源
定量研究東江源地區(qū)土壤侵蝕量需從相關(guān)網(wǎng)站收集研究區(qū)域的多種類型的原始數(shù)據(jù),其中有東江源地區(qū)30米分辨率DEM數(shù)據(jù)、Landsat 8 OLI_TIRS遙感影像、土壤類型分布數(shù)據(jù)、降雨數(shù)據(jù)以及東江源范圍等。
1.1.1DEM數(shù)據(jù)及遙感影像
地形圖是在地理空間數(shù)據(jù)云(www.gscloud.cn)免費(fèi)下載的東江源地區(qū)DEM 30M數(shù)據(jù),將獲得的 DEM 在 ArcGIS10.5 軟件中進(jìn)行裁剪獲得研究區(qū)域東江源地區(qū) DEM 數(shù)據(jù)。遙感影像是在地理空間數(shù)據(jù)云(www.gscloud.cn)免費(fèi)下載的東江源地區(qū)的Landsat 8 OL_ TIRS 影像數(shù)據(jù)。
1.1.2降雨數(shù)據(jù)及土壤數(shù)據(jù)
本次研究區(qū)域的降雨數(shù)據(jù)來源是在 WheatA 軟件上下載的東江源地區(qū)所包括的尋烏縣和定南縣的 2016到2020年年降雨量統(tǒng)計(jì)數(shù)據(jù)。土壤數(shù)據(jù)是在資源環(huán)境科學(xué)與數(shù)據(jù)中心(www.resdc.cn)免費(fèi)下載的 1:100中國土壤類型空間分布數(shù)據(jù)。
1.2 數(shù)據(jù)預(yù)處理
獲得的數(shù)據(jù)往往不能完全滿足研究的需求為此研究數(shù)據(jù)處理之前往往需要對數(shù)據(jù)進(jìn)行預(yù)處理,這樣可以使研究更加的方便。本文根據(jù)研究目的和CSLE模型對dem數(shù)據(jù)、遙感影像數(shù)據(jù)等在ENVI中進(jìn)行輻射定標(biāo)、大氣校正、圖像鑲嵌和裁剪等。
2.CSLE模型
水土流失量計(jì)算,采用SL190—2007《土壤侵蝕分類分級標(biāo)準(zhǔn)》推薦的由劉寶元等提出的中國通用土壤流失方程(CSLE模型)
A=R×K×LS×B×E×T
式中,A為單位面積上時間和空間平均的土壤流失量, t/(hm2?a);R為降水侵蝕力因子,K為土壤可蝕性因子,LS為坡長坡度因子,B為植被覆蓋度與生物措施因子,E為工程措施因子,T為耕作措施因子。將R、K、L、S、B、E、T因子的柵格數(shù)據(jù)進(jìn)行空間相乘。根據(jù)國家水利部頒布的《土壤侵蝕分類分級標(biāo)準(zhǔn)》(SL190—2007)將土壤侵蝕模數(shù)分為0~200、200~2500、2500~5000、5000~8000、8000~15000、>15000t?km-2?a-1共 6個等級。
3.CSLE模型因子處理
3.1 土地利用分類
利用ENVI軟件,首先根據(jù)東江源范圍矢量文件,利用ROI裁剪圖像,其次進(jìn)行輻射定標(biāo)和大氣校正;然后利用監(jiān)督分類中的最大似然法進(jìn)行ROI感興趣區(qū)分類,最后進(jìn)行分類合并,最終實(shí)現(xiàn)土地利用類型分類。
3.2 降水侵蝕力因子(R)
降雨侵蝕力因子:是指降雨導(dǎo)致土壤侵蝕發(fā)生的潛在能力,用一次降雨總動能E與該次將于最大30min雨強(qiáng)I30的乘積EI30表示。反映了雨滴對土壤顆粒的擊濺分離以及降雨形成徑流對土壤沖刷的綜合作用。
我們綜合了繆馳遠(yuǎn)等、王明曉和史冬梅等學(xué)者研究提出的基于年降雨量數(shù)據(jù)的R因子計(jì)算模型,計(jì)算公式如下:
式中:
Rn—年降雨量侵蝕力因子,
;Pn—年降雨量,mm。
3.3 土壤可蝕性因子(K)
土壤可蝕性(k)指的是土壤容易受到侵蝕破壞的性質(zhì),土壤可蝕性因子主要受到內(nèi)因(土壤質(zhì)地等)和外因(人類活動等)影響。在內(nèi)外兩個方面的相互作用下,土壤可蝕性因子的值會隨之改變。本項(xiàng)目采用EPIC模型來獲取k值,是一種專門為土壤侵蝕與土地生產(chǎn)力設(shè)計(jì)的模型, 基于土壤有機(jī)碳和粒徑分布進(jìn)行土壤可蝕性K值估算。Williams等在EPIC模型中K值的計(jì)算公式如下:
上式中, S1為砂粒含量, 單位為%;S2為粉粒含量,單位為%;n為黏粒含量,單位為 %;C為有機(jī)碳含量, 單位為%;S3=1-S1/100,單位為 %。
3.4 坡度坡長因子(LS)
地形因子直觀表現(xiàn)由坡長(L)和坡度(S)兩個因子決定。坡長通常是指地面上一點(diǎn)沿水流方向到其流向起點(diǎn)間的最大地面距離在水平面上的投影長度。坡度表示了該局部地表坡面的傾斜程度。坡度越大、坡長越長,其地勢越陡峭,越容易發(fā)生土壤侵蝕。
本次研究,坡度坡長因子(LS)通過遙感技術(shù)獲取的數(shù)字高程模型快速獲取并通過DEM數(shù)據(jù)在GeoScene Pro中完成。根據(jù)東江源地區(qū)的地形地勢,此次研究選用ASTER GDEM 30M數(shù)據(jù),采用符素華、劉寶元等[1]提出的公式。
??
式中:S為坡度因子,θ為坡度,L為坡長因子,λ為坡長(m),22.13為標(biāo)準(zhǔn)小區(qū)坡長,m為坡度坡長指數(shù)
3.5 植被與生物因子(B)?
采用MODIS遙感數(shù)據(jù)的NDVI產(chǎn)品MOD13Q1,利用NDVI數(shù)據(jù)計(jì)算植被覆蓋度FVC,根據(jù)FVC計(jì)算不同地類植被覆蓋與生物措施因子B。歸一化植被指數(shù)(NDVI)是反映土地覆蓋植被狀況的一種遙感指標(biāo),定義為近紅外通道與可見光通道反射率之差與之和的商,即NDVI= (NIR-R)/(NIR+R),NIR、R為近紅外波段、紅光波段的反射值。
3.6 水土保持工程措施因子(E)和耕作措施因子(T)
根據(jù)水土保持工程措施因子賦值表和全國輪作措施及代碼表,對應(yīng)東江源區(qū)土地利用類型分布圖,我們得到了相應(yīng)的水土保持工程措施因子以及耕作措施因子并給予賦值(當(dāng)土地利用類型為耕地時,E=0.335,T=0.338;其他均賦值為1)。
4.地理探測器
地理探測器是一種統(tǒng)計(jì)學(xué)方法,包括因子探測器、交互作用探測器、生態(tài)探測器和風(fēng)險(xiǎn)探測器4個部分,通過探測要素的空間分層異質(zhì)性來揭示某原因背后的驅(qū)動力。
q值計(jì)算公式如下:
式中:h=1,…,m 為變量或因子的分層,即分類或分區(qū);𝑛?和n分別為層h和全區(qū)的單元數(shù);𝜎?2和𝜎2分別為層h和全區(qū)變量的方差。SSW 和 SST 分別為層內(nèi)方差之和和全區(qū)總方差。q 值區(qū)間為[0,1],其 q 值越大,該影響因素對土壤侵蝕空間分布的解釋力越強(qiáng)。此處我們選擇多年平均降雨量、坡度、植被覆蓋度為土壤侵蝕影響因子,多年平均土壤侵蝕數(shù)模作為因變量,模型運(yùn)行前對連續(xù)性變量進(jìn)行分類,按照地理探測器公式對各級別數(shù)值和方差進(jìn)行計(jì)算,得到q值。
5.皮爾遜相關(guān)系數(shù)
本研究使用皮爾遜相關(guān)系數(shù)探究降水、坡度和植被對土壤侵蝕的影響。該方法在研究土壤侵蝕的影響因素分析中廣泛研究,計(jì)算公式如下:
其中:n為總的年數(shù),xi和yi分別表示第i年的土壤侵蝕量和影響因素(降水、坡度、植被),x均和y均分別為土壤侵蝕了量和影響因素的均值,rx.y表征x與y的相關(guān)系數(shù)。??
三、結(jié)果分析
1.各類結(jié)果分布圖
圖2 2016-2020年東江源區(qū)土地利用類型分布圖
圖3 土壤侵蝕各類因子分布圖
圖4 2016-2020東江源區(qū)域土壤侵蝕模數(shù)分布圖
2.土壤侵蝕影響因素分析
2.1 坡度因素
表 1東江源區(qū)域各年份坡度、面積、侵蝕模數(shù)與侵蝕量
研究區(qū)域劃分為5個等級來研究2016~2020年土壤侵蝕與坡度的關(guān)系。從得到不同年份的數(shù)據(jù)可知,侵蝕模數(shù)隨著坡度的增加呈現(xiàn)遞增趨勢,并在第5級達(dá)到最大。在1~2級內(nèi),隨著坡度的增加,土壤侵蝕量呈現(xiàn)大量增加的趨勢,并在第2級達(dá)到最大,在3~5級內(nèi),隨著坡度的增加,土壤侵蝕量呈現(xiàn)遞減的趨勢,并且在第5級達(dá)到最小。由此可以推測該研究區(qū)域土壤侵蝕主要在低山、丘陵等地。
2.2 植被因素
由2016—2020 年東江源區(qū)域不同植被覆蓋度等級土壤侵蝕模數(shù)可知隨著植被 ,覆蓋度的升高,土壤侵蝕模數(shù)降低。當(dāng)植被覆蓋度<30%時,土壤侵蝕模數(shù)最大,易誘發(fā)土壤侵蝕;當(dāng)植被覆蓋度>60%時,土壤侵蝕模數(shù)最小。從不同植被覆蓋度土壤侵蝕模數(shù)總數(shù)來看,當(dāng)植被覆蓋度<30%時,土壤侵蝕模數(shù)總數(shù)為占總土壤侵蝕比重的 40.41%??梢娡寥狼治g主要發(fā)生在植被覆蓋度<30%的區(qū)域,植被覆蓋度處于30%-45%之間與處于45%-60%之間時的土壤侵蝕模數(shù)變化較小可能是因?yàn)?016-2020年間不同植被覆蓋率土壤面積相互轉(zhuǎn)化導(dǎo)致面積波動較大。
表 2東江源區(qū)域各年份植被覆蓋度、面積、侵蝕模數(shù)與侵蝕量
2.3 降水因素
由2016—2020 年東江源區(qū)域降雨量對應(yīng)土壤侵蝕模數(shù)圖表,可知降雨量的升高,土壤侵蝕模數(shù)升高,侵蝕量也隨之升高。在2016-2020年間2016年的降雨量平均值為最高值,與之對應(yīng)的侵蝕模數(shù)也達(dá)到了了年間的最大值。
表 3東江源區(qū)域各年份降水量、侵蝕模數(shù)平均值與侵蝕量
3.皮爾遜指數(shù)結(jié)合地理探測器
利用SPSS軟件進(jìn)行皮爾遜相關(guān)性分析,根據(jù)所得結(jié)果可知,土壤侵蝕模數(shù)與降雨量為高相關(guān)(0.980),土壤侵蝕模數(shù)與植被覆蓋度為中相關(guān)(0.300)。由于DEM數(shù)據(jù)為多年不變,無法用皮爾遜相關(guān)性加以研究,故采用地理探測器的方法加以研究。根據(jù)q值表最終數(shù)據(jù)可知,q值由大到小為:降雨量(0.42)>植被覆蓋度(0.22)>坡度(0.07)因此無論是皮爾遜相關(guān)性分析還是地理探測器,均表明了降雨對土壤侵蝕強(qiáng)度的解釋力最強(qiáng),是影響土壤侵蝕空間格局分布的主要因素;其次為植被覆蓋度,植被覆蓋度與土壤侵蝕密切相關(guān);而坡度對土壤侵蝕的影響則較小。
表 1皮爾遜指數(shù)相關(guān)性及地理探測器最終值表
四、創(chuàng)新之處
1.關(guān)鍵技術(shù)創(chuàng)新
中國CSLE模型相比美國RUSLE模型,其LS因子公式中的角度范圍由10°改為14°,更加符合中國地形特征,也更適用于地形復(fù)雜的東江源區(qū),得出的結(jié)果更有利于了解東江源地區(qū)的土壤侵蝕情況。對于大部分關(guān)于水土流失分析的作品中,都沒有包括對于影響土壤侵蝕的各因素的分析,本作品加入了地理探測器與皮爾遜相關(guān)性的結(jié)合分析,對影響土壤侵蝕的各因素(降水、植被、坡度)進(jìn)行了相關(guān)研究。
2.應(yīng)用范圍創(chuàng)新
東江源流域近年來無人對其進(jìn)行土壤侵蝕研究,通過本次研究,可以了解近年來東江源流域的水土流失情況,為近年來該地區(qū)土壤情況提供補(bǔ)充,為以后的研究者提供研究材料。同時通過本次研究,可以以東江源地區(qū)為例,提出對該地區(qū)的保護(hù)策略,從而為東江源區(qū)域以及南方同類型山地丘陵地域的水土保護(hù)做出參考。
總結(jié)
以上是生活随笔為你收集整理的2021地理设计组二等奖:基于GIS的东江源区土壤侵蚀及其影响因素空间分析的全部內(nèi)容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: c语言--找钱
- 下一篇: Mysql数据备份和数据恢复命令