基于GIS的USLE方程各因子计算方法研究.docVIP

基于GIS的USLE方程各因子计算方法研究 摘要 传统方法估算土壤侵蚀量限制因素过多，usle方程为土壤侵蚀量估算提供了一般算法，但是usle方程中各因子的求算没有通用方法，在各个研究区域的运用难以量化。对基于gis的usle方程各因子计算方法进行探讨，给出了基于gis技术usle估算流域土壤侵蚀量的一般方法，以提高usle方程在我国的可用性。 关键词 土壤侵蚀量；gis；usle；计算方法 中图分类号 s157 文献标识码 a 文章编号 1007-5739（2013）01-0206-02 土壤侵蚀是破坏土地资源的重要因素，导致土地退化，严重影响耕地的数量和质量。我国是土壤侵蚀最为严重的国家之一，约有38%的国土属于土壤侵蚀区，开展土壤侵蚀研究非常重要，预测土壤侵蚀量对预防和治理土壤侵蚀有着重要意义。 1 usle方程概述 用传统的测量学方法估算土壤侵蚀量有诸多弊端，如高程实测法、直接丈量法等方法要求样本数量过大，工作周期过长，限制因素过多。美国学者wischmeier和smith于20世纪50年代提出了通用土壤流失方程（universal soilloss equation，usle）[1]，该模型建立在土壤侵蚀理论及大量实地观测数据统计分析的基础上，方程对各因子的解释都具有物理量化意义，至今仍是预测土壤侵蚀量最为广泛的模型。通用土壤流失方程的表达式为： a=r×k×l×s×p×c（1） 式（1）中，a为年平均土壤流失量，单位t/hm2；r为降雨侵蚀力因子，单位m·t·cm/（hm2·h）；k为土壤侵蚀性因子，单位t·hm2·mj·mm；l为坡长因子，无量纲；s为坡度因子，无量纲；p为水土保持措施因子，无量纲；c为地表植被覆盖及经营管理因子，无量纲。 自通用土壤流失方程提出以来，国内外很多学者运用该模型成功预测了很多流域的土壤侵蚀量，但是usle各因子求算方法没有通用的方法，只能根据试验研究区域自己定义，难以推广到我国普遍区域。本文对基于gis的usle方程各因子计算方法进行探讨，给出了基于gis技术usle估算流域土壤侵蚀量的一般方法。 2 数据来源与处理 研究用数据包括研究区域tm遥感影像、数字高程模型（dem）、土壤类型分布图、土地利用现状图。 运用arc/info的空间数据管理和分析功能，建立研究流域的数字高程模型（dem）、土壤类型分布图、土地利用现状图等矢量图，对其属性数据进行相应的数据编码操作，运用arcview对以上各图件进行grid格式栅格化处理，结合流域尺度和影像精度，进行栅格大小的划分，作出usle方程的因子图，通过arcview将各因子栅格图作连乘运算，获得研究流域的土壤侵蚀量及分布图。 3 usle方程各因子的确定 3.1 降雨侵蚀力因子r值的确定 美国学者wischmeier和smith利用美国35个土壤保持试验站8 250个试验分区的降雨侵蚀资料统计得出降雨侵蚀力的经典计算方法，简称威式经典公式[2]： r=1/100∑（e·i30）（2） 式（2）中，r为1次降雨的侵蚀力值，单位m·t·cm/（hm2·h）；e为1次降雨的总动能，单位m·t/hm2；i30为1次降雨过程中最大30 min的降雨强度，单位cm/h。 由于我国地域辽阔，气候因子变化很大，很难直接应用国外学者的r因子计算式。王万忠等统计分析了全国12个站点的313个站年资料后，提出了我国估算年r值的计算公式[3]： r=0.207（p·i60/100）1.205 （3） 式（3）中，p为年降雨量，单位mm；i60为年最大60 min降雨量，单位mm。该式提供了计算我国年r值的简便方式，具有良好的准确率。 根据式（3）计算出研究流域降雨侵蚀力因子，进而运用arcview生成r因子图。 3.2 土壤侵蚀性因子k值的确定 目前的研究确定k值的计算方法可分为3种：诺谟图法、粒级计算法和查表法。wischmeier于1971年提出了诺谟图法，计算公式为： k=[2.1×10-4（12-a）m1.14+3.25（b-2）+2.5（c-3）]/100（4） 式（4）中，m为（粉砂+极细砂百分数）乘以（100-黏土百分数）；a为有机质分数；b为土壤结构编号；c为剖面渗透性等级。 romkens于1986年分析土壤颗粒粒级和k值相关性后，建立了颗粒粒级计算法。1987年，kirkby建立了查表法，通过查看表中对应的土壤质地以及有机质含量得到k值[4]。由于运用诺谟图法和粒级计算法计算k值均需要较多参数，并且这些参数难以精确，因此查表法是适用性最好的取k值法。 1994年，卜兆宏选取具体的研究对象，对kirkby查表法进行修正并提出k值修正系数表[5]，修正原表中只列出有机质小于等于4%的土壤而缺乏高有机质含量土壤的k值。为使其适用于高有机质含量的土壤