基于降雨资料的北京山区坡面土壤侵蚀方程研究.docxVIP

基于降雨资料的北京山区坡面土壤侵蚀方程研究 土壤侵蚀是全球的问题，越来越受到重视。北京山区也存在较为严重的土壤侵蚀。北京山区面积为1.04万km2,占全市总面积的62%。根据2004年北京市水土流失监测公报,山区土壤侵蚀面积为4 089 km2。严重的水土流失不仅造成当地生态环境的恶化,而且土壤中的养分、农药和杀虫剂等也随水土流失进入下游水体,造成下游水体的水质污染。据统计,在密云水库周边点源污染得到初步控制的情况下,畜禽养殖、农用化肥、农药污染已成为水源地的主要污染源。根据监测资料分析结果,密云水库上游73%COD,71%BOD,94%NH3-N,75%TN和95%TP都来自非点源污染。为了防治水源污染,治理和控制北京山区水土流失具有重要意义。坡面土壤侵蚀模型是进行水土资源管理、定量评价水土保持措施和进行土壤侵蚀灾害评价的技术工具。 美国通用土壤流失方程USLE基本包括了影响坡面土壤流失的主要因素,公式建立所用资料范围较广,并且统一了侵蚀模型形式,从而使USLE方程在世界各国得到了广泛应用。我国对坡面土壤侵蚀模型也有大量的研究。为了进行区域的土壤侵蚀评价和预报,一些研究者已结合区域实际情况,采用USLE方程的形式,建立了各自的坡面土壤侵蚀预报模型[9,10,11,12,13,14,15]。到目前为止,对北京坡面土壤侵蚀的预报研究还较少。自1987年以来,北京市水土保持工作总站已建立了大量坡面径流试验小区,获取了约1 000个小区的年降雨径流和泥沙资料。这些试验观测资料为北京坡面土壤侵蚀预报方程的建立提供了很好的基本数据资料。 1 年降雨径流泥沙量计算 美国通用土壤流失方程已在许多国家和地区得到了应用,因此,北京土壤流失方程也采用USLE的形式,并考虑北京的实际情况,采用的估算方程形式如下: A=RKLSCPRF(1) 式中:A为土壤侵蚀模数,t·hm-2·a-1;R为降雨侵蚀力,MJ·mm·h-1·hm-2·a-1;K为土壤可蚀性因子,t·hm2·h·MJ-1·hm-2·mm-1;L为坡长因子;S为坡度因子;C为田间覆盖与管理因子;P为水土保持措施因子;RF为砾石覆盖因子。 自1987年以来,北京市水土保持工作总站分别在平谷挂甲峪、密云石匣和大关桥、怀柔汤河口和三渡河、延庆上辛庄、昌平下口、门头沟清水和田寺以及房山蒲洼建立了大量坡面径流小区,获取了近1 000个小区的年降雨径流泥沙资料,以观测不同土壤、土地利用、植被覆盖对土壤侵蚀的影响。 本研究利用北京市水土保持工作总站资料年限较长的密云石匣、怀柔汤河口、延庆上辛庄、房山蒲洼设置的75个径流试验小区降雨和泥沙资料,研究不同土壤、土地利用、植被覆盖对土壤侵蚀的影响。具体情况见表1。 径流试验小区用水泥板围圈,两侧设0.3 m宽隔离区,上部设排水沟,下部连接集流桶。试验区的作物田间管理与当地水平一致,种植作物品种为当地普遍种植的玉米和大豆。观测降雨量采用自记雨量计,根据记录结果分析每次降雨的雨量、降雨时间和雨强。 观测流失泥沙量方法是在每次降雨后即刻观测记载集流桶中水深,并计算出地表径流总量,然后将桶中泥水充分搅匀,取2 000 mL水样带回实验室过滤烘干,计算出含沙量,从而得出每次降雨各小区土壤流失量,将一年中各次土壤流失量累加计算得到各小区年土壤流失量。 在进行天然径流小区观测的同时,也进行了人工降雨试验,用以研究砾石覆盖对土壤侵蚀的影响。试验考虑了0、5%、10%、20%、40%和60%6个不同砾石覆盖度,7°和20° 2种不同坡度和30 mm/h、60 mm/h和90 mm/h 3种不同雨强,用长1 m宽0.5 m深0.5 m的试验土盘共进行了60场人工降雨模拟试验,用以得到砾石覆盖因子值。 为了客观评价土壤流失方程及其参数的合理性,本文选用了模型确定性系数Ef这一指标来评价土壤流失方程预测结果的精度,该指标由Nash和Sutcliffe1970年提出,用以评价计算值和实测值的接近程度。计算式为 Ef=1-∑(Robsi-Rcali)2∑(Robsi-ˉRobs)2∑(Robsi?Rcali)2∑(Robsi?Rˉˉˉobs)2(2) 式中:Robs为实测土壤流失量的平均值;Robsi为实测土壤流失量;Rcali为计算土壤流失量。 如果模型确定性系数为1,表示模型计算结果与实测值非常接近;模型确定性系数为0,表示计算值与实测值之差的平方和刚好等于实测值与其平均值的平方和,这表明可以用实测值的平均值来进行预报,其预报效果与模型预报的效果一样;模型确定性系数为负,则表示实测值的平均值比模型预报值能更好地进行土壤流失量的估计。 2 结果与分析 2.1 日降雨侵蚀力计算 一般将降雨总动能E和最大30 min雨强I30的乘积EI30作为降雨侵蚀力指标。该指标在计算的过程中涉及到