《土壤侵蚀的估算方法》由会员分享,可在线阅读,更多相关《土壤侵蚀的估算方法(11页珍藏版)》请在金锄头文库上搜索。
1、侵蚀的估算方法数据处理作者:牛健平时间:2011年10月11日北京天合数维科技有限公司目录(CONTENT)一、所需数据与参数31、所需数据32、所需中间参数32.1、水土保持因子P 32.2、地标覆盖因子C 32.3、地形因子LS 42.4、土壤可视性因子K 42.5、降水侵蚀因子R 43、所需参数53.1、潜在土壤侵蚀量Ap 53.2、现实土壤侵蚀量Ar 53.3、土壤保持量Ac 54、指标结果参数54.1、保护土壤肥力的经济效益必64.2、减少土地废弃的经济效益励64.3、减轻泥沙淤积的经济效益助6二、处理流程71、DEM数据的处理81. 1、坡长 L 81.2、百分比坡度a 81.3、
2、地形因子LS 92、气象数据92. 1、月降雨量Pi的计算92.2、土壤侵蚀力指标R 103、土壤类型数据104、遥感影像数据105、土壤理性化数据11三、所需参数的计算11四、指标结果参数计算11一、所需数据与参数在计算的过程中,总共涉及到的数据有地形数据、遥感影像数据、气象数据、 土壤类型数据、土壤理性化数据以及统计数据,涉及到的中间参数有水土保持因 子P,地标覆盖因子C,地形因子LS, 土壤可视性因子K,降水侵蚀因子R,所 需要的参数有潜在土壤侵蚀量Ap,现实土壤侵蚀量Ar, 土壤保持量Ac,指标结 果参数有保护土壤肥力的经济效益Ef,减少土地废弃的经济效益Es,减轻泥沙 淤积的经济效益
3、En。1、所需数据在进行土壤侵蚀的估算过程中,需要以下数据:A、地形数据;B、遥感影像数据;C、气象数据,主要是降雨量数据;D、土壤类型数据;E、土壤理性化数据;F、统计数据。2、所需中间参数在数据处理的过程中,所涉及到的中间参数与计算公式如下。2.1、水土保持因子P按照游松财的方法,水田的P值取0.15,其他土地利用方式基本没有采取水 土保持措施,因此取值为l.OOo2.2、地标覆盖因子C地表覆盖因子是根据地面植被覆盖状况不同而反映植被对土壤侵蚀影响的因素,与土地利用类型、覆盖度密切相关。C值的估算采用如下公式:c = l,fC = Oc = c = 0.6508-0.34361g fC,0
4、 fC 78.3%式中,fC表示植被覆盖度,由遥感影像通过NDVI计算。2.3、地形因子LS通过数字高程模型(DEM),计算获得坡长和坡度,然后根据黄炎和等建 立的方程式,获得LS的空间分布特征。LS = 0.08L0356t05式中,L,坡长(m); a,百分比坡度。2.4、土壤可视性因子K釆用陈明华等人建立的土壤可蚀性K值的计算公式:/r = 10 3 (160. 80 - 2. 31 尤 1 + 0. 38尤 2 + 2. 26尤 3 + 1. 31 尤 4 + 14. 67Z5)式中:可蚀性值(美国习用单位);尤1、尤2、尤3、尤4、尤5分别表示 细砾(1 3mm )%、细砂(0. 0
5、50. 25mm )%、粗粉粒(0. 01 0. 05 mm )%、 细粉粒(0. 0050. 01mm )%、有机质(10g/kg) 02.5、降水侵蚀因子R采用周伏建和黄炎和等人(1997,2000,2002)根据实测数据提出的R值计算 式,该值考虑了月降水量。12/?=工(1.5527 + 0.1792H)1=1式中,R为降雨侵蚀力指标(j. cm/m. h), Pi为月降雨量(mm)。月降水量Pr. 根据气象站点观测的降水点位及数据,进行克里金插值获得。3、所需参数在计算的过程中主要涉及到潜在土壤侵蚀量Ap,现实土壤侵蚀量Ar以及土 壤保持量Ac三个参数,可以通过以上5个中间参数来获得
6、,其计算公式如下。3.1、潜在土壤侵蚀量Ap潜在土壤侵蚀量不考虑地表覆盖类型和土地管理因素,即c=l, p=l,此时, USLE的形式为:Ar = RxKx LS式中,Ap为潜在土壤侵蚀量(吨/平方千米),R为降雨侵蚀力指标,K为土壤可 侵蚀因子,LS为坡长坡度因子。3.2、现实土壤侵蚀量Ar现实土壤侵蚀量主要考虑了地表覆盖类型和土地管理因素,其计算式为:At = RxKx LS x Cx P式中,Ar为现实土壤侵蚀量(吨/平方千米),R为降雨侵蚀力指标,K为土壤可 侵蚀因子,LS为坡长坡度因子,C为地表覆盖因子,P为土壤保持措施因子。3.3、土壤保持量Ac由上两式可得到土壤保持量:AC_AP
7、-A,.式中,A;,为潜在土壤侵蚀量(吨/平方千米),A为现实土壤侵蚀量(吨/平方千米), A为土壤保持量(吨/平方千米)。4、指标结果参数指标结果参数总共涉及到的数据有保护土壤肥力的经济效益侈;减少土地废 弃的经济效益厉和减轻泥沙淤积的经济效益劭三个参数,其计算公式与所需参 数值见下。4.1、保护土壤肥力的经济效益耐土壤侵蚀使大量土壤营养物质流失,主要表现在N、P、K等营养物质的 流失。Ef = (A x De x C, x R x C) + 工 A C /10000(/ = N、P、K)9I式中,E为保护土壤肥力的经济效益(元/ a); A( 土壤保持量(t/hm2a); Dt., 土壤密
8、度(t/m3); Co, 土壤有机质含量(kg/t); R,薪柴转化为土壤有机质的 比例,一般为2; C,薪柴的平均价格(元/kg); Ci, 土壤中N、P、K的纯含量; P., N、P、K的价格。注:N、P、K的价格:参照国家林业局2008年4月28日发布的中华人民共和 国林业行业标准一森林生态系统服务功能评估规范,氯化钾价格为2200元/吨, 磷酸二铁价格为2400元/吨。元素名称价格(元/吨)N509. 04P563. 52K1151.7表1、NPK价格表4.2、减少土地废弃的经济效益根据土壤保持量和土壤表土平均厚度(0.6m)来推算因土壤侵蚀而造成的废 弃土地面积,再根据机会成本法计算
9、因土地废弃而失去的年经济价值。Es=Ac4-P4-0. 6XB4-10000式中,Es,减少土地废弃的经济效益(元/a); Ac, 土壤保持量(t/a); P, 土壤 的容重(t/m3); B,湿地年均收益(元/血)。注:土壤容量取26.5吨/立方米,湿地年均收益根据崔丽娟(2004年)研究推 算,湿地年均收益为245.5元/平方米。4.3、减轻泥沙淤积的经济效益按照我国主要流域的泥沙运动规律,全国土壤侵蚀流失的泥沙有24%淤积于水库、江河、湖泊,这部分泥沙直接造成了水库江河、湖泊蓄水量的下降,在一定程度上增加了干旱、洪涝灾害发生的机会,因此可根据蓄水成本计算损失价值。En=Ac4-pX24%
10、XC式中,En,减轻泥沙淤积的经济效益(元/a); Acf 土壤保持量(t/a); P, 土壤 的容重(t/m;1); C,水库工程费用(元/m)。注:土壤容量取26.5吨/立方米,核算过程中采用的水库清淤工程费用为6.94 元/立方米,此数据来源于硕士论文。二、处理流程总体处理流程如下图所示:图1、数据处理流程图1. DEM数据的处理运用DEM数据得到的最后结果是求出了地形因子LS,在求LS的过程中,需 要坡长L和百分比坡度a,其求解主要在Arcmap中完成。1.1.坡长L坡长的求解主要参考汤国安老师的ArcGIS地理信息系统空间分析实验教程中的第十一章中的水文分析来完成,Arcmap中建立
11、模型如下图所示:图2、无洼地DEM生成模型图3、生成坡长L模型1.2. 百分比坡度a百分比坡度a是用slope来实现的。图4、百分比坡度的生成13、地形因子LS地形因子LS计算的模型为:图5、地形因子LS的计算模型2、气象数据根据气象站点观测的降水点位及数据,进行克里金插值获得月降雨量数据Pi,再根据Pi计算土壤侵蚀力指标R。2.1、月降雨量Pi的计算2.2、土壤侵蚀力指标R土壤侵蚀力指标 P 是在 Arcmap 下 Spatial Analyst | Raster Caculator中计算而来。Spatial Analyst Layer: |lsDista neeBDgn sity.Inte
12、rpolate to RasterSurface AnalysisCell Statistics. Neighborhood Statistics.Zonal Statistics.Zonal Histogram Reclassify.Raster Calculator.图6、栅格计算3、土壤类型数据土壤可蚀性因子K是通过土壤类型数据直接对其进行赋值。4、遥感影像数据通过遥感影像数据获得地表覆盖因子C,首先求NDVI, NDVI的求解在ENVI 中实现,主菜单中transformndvi,算归一化植被指数。接着根据NDVI数据 来求解植被覆盖度FC,其计算公式为:FC=(NDVI-NDVI_M
13、IN)/(NDVI_MAX-NDVI_MIN)式中:NDVI_MIN为NDVI的最小值,DVI_MAX为NDVI的最大值。由此,我们可以根据公式c = l,/C = 0c = c = 0.6508-0.3436 lg /C, 0/C78.3%来计算 C 值,在 Arcmap 中 Spatial Analyst | Raster Caculator实现。Is snk snkdep waters 卜 shink zonalmaxl zonalmin 0 & fcimg *0.781 0.6508 - 0 3436 * LoglO(fcimg)? 0)5、土壤理性化数据按照游松财的方法,水田的P值取
14、0.15,其他土地利用方式基本没有采取水 土保持措施,因此取值为1.00。三. 所需参数的计算所需要的参数有潜在土壤侵蚀量Ap,现实土壤侵蚀量Ac 土壤保持量Ac, 在进行了上述的计算以后,其计算比较简单,根据公式,在Spatial Analyst I Ra航er Caculator中输入对应公示即可实现。最后对Ac进行汇总统计,在Spatial Analyst | Zonal statistics中进行。Spatial AnalysLwyer: isDistanceIDensity.Interpolate to RasterISurface AnalysisICell Statistics.Neighborhood Statistics.Zonzl Statistics.四、指标结果参数计算这一步只是简单的数值计算,在此不再详述,但是在计算的过程中注意单位 的统一。注:在进行统计Ac总和的时候,要乘以图像分辨率。还有单位一定要统一!