梨树县农业非点源污染磷指数风险评价研究.docVIP

梨树县农业非点源污染磷指数风险评价研究 　　【摘 要】随着社会经济的发展，农业非点源污染问题已日益严重。利用磷指数法对梨树县的农业非点源污染磷流失风险进行评价。根据磷素的不同流失程度、风险等级，提出相应的管理措施及建议。 　　【关键词】非点源污染；磷指数；地理信息系统；危险性；风险评价 　　1.研究区域概况 　　研究区域自然条件及社会经济概况。 　　梨树县位于吉林省西南部，辖内21个乡镇。属北温带半湿润季风大陆性气候。地上水资源比较丰富，但开发利用很不充分。农业是梨树县的基础产业，素来就有黄金玉米带之美誉。 　　2.农田磷指数风险评价体系的建立 　　2.1基本数据信息获取 　　本研究用到的数据主要有图形资料、降雨量数据、研究区域调查资料等。 　　2.2空间数据库的建立 　　通过地理信息系统软件对信息数据的提取，数据赋值，空间插值，公式计算，图层叠加等一系列步骤，建立评价指标体系，得到磷指数值，确定研究区磷流失的不同危险等级，从而进行研究区域内的磷指数风险评价。 　　3.磷指数法的具体应用 　　本文进行评价的方法，磷指数法采用宾夕法尼亚的索引版本8/2002。 　　 　　第一部分：筛选工具 　　当研究区域满足土壤有效磷（Mehlich-3）200mg P kg-1，或满足污染源距河流的贡献距离45.7m时，即满足其中任意一条则进入第二部分评价。 　　第二部分：因子 　　A部分：来源因子 　　表3-1 来源因子 　　Table 3-1 source factors 　　土壤有效磷比率=0.20×（Mehlich-3）P（mg Pkg-1） （3-1） 　　化肥磷比率=化肥P含量×方法 （3-2） 　　有机肥磷比率=有机肥P含量×方法×有效性 （3-3） 　　来源因子=土壤有效磷比率+化肥磷比率+有机肥磷比率 （3-4） 　　B部分：迁移因子 　　表3-2 迁移因子 　　Table 3-2 Transport factor 　　*或在小河附近迅速的渗透性土壤 　　运输总和=土壤侵蚀+径流潜力+地表水体连通性+径流贡献距离 （3-5） 　　迁移因子=修改连接因子×（运输总和/22） （3-6） 　　磷指数（PI）值=2×来源因子×迁移因子 （3-7） 　　表3-3 磷指数的危险等级和管理指导 　　Table 3-3 Risk classes and management guidance for the phosphorus index 　　3.1来源因子的获取 　　3.1.1土壤有效磷 　　将研究区域各采样点获得的土壤速效磷含量转化为（Mehlich-3）P[1]，通过公式（3-1），转化为土壤有效磷比率，生成梨树县的土壤有效磷比率的空间分布数据图。经计算，梨树县土壤有效磷比率的均值为26.82mg/kg。 　　3.1.2施肥因子 　　通过实地调查，研究区域的施肥时间比较集中，化肥施用方式多以基肥为主，施肥时间主要集中在4-8月，一般施在土壤以下2cm左右。施肥多依赖化肥，所以本次评价不考虑有机肥的因素。结合当地施肥现状，将施肥因子中各相关方面对于不同地点给予相应赋值。 　　3.2迁移因子的获取 　　3.2.1土壤侵蚀量的计算 　　美国土壤流失通用方程（USLE）是目前估算高地侵蚀造成的土壤流失最普遍的估算式，也是当今世界上应用最为广泛的水蚀模型。表达式为： 　　A=R×K×L×S×C×P （3-8） 　　（1）降雨侵蚀力因子R。 　　本研究应用的是由梨树县气象局提供的近十年的降雨数据。整理得出多年月平均降雨量及年平均降雨量，利用降雨侵蚀力因子R的计算公式，求得该区的年降雨量为468.36mm，R值为26.16。 　　（2）土壤可蚀性因子K。 　　土壤可蚀性因子K反映不同土壤的潜在侵蚀性，K值越大，侵蚀敏感度越高，越易被侵蚀。相反，K值越小，侵蚀敏感度越低，也越不容易被侵蚀。K值大小与土壤粒度结构、土壤物理性质、机械组成、土壤有机质含量、土壤渗透性等息息相关。本文采用1995年，陈明华等人[2]建立的土壤可蚀性K值的计算公式： 　　K=10-3（164.8-2.31X1+0.38X2+2.26X3+1.31X4-14.67X5） （3-9） 　　利用ArcGIS软件，以梨树县土壤类型图作为底图，对研究区域gird图层进行相应的赋值，在计算机上可以获得土壤可蚀性K值空间分布。 　　（3）坡长、坡度因子LS。 　　有关地形因子LS的计算，本研究采用的是由W.H.维希迈尔和D.D.史密斯提出的，在我国应用较多[57-59]的计算公式。利用数字高程图确定坡度因子S以及坡长因子L。 　　LS=（L/22）0.3（θ/5.16）1.3 （