长江上游重点水土流失区坡耕地土壤流失速率研究.docxVIP

长江上游重点水土流失区坡耕地土壤流失速率研究 山坡上的土地是长江上游河流的主要水源之一。据统计,实施“退耕还林还草”工程前,长江上游坡耕地面积约666.67万hm2。其中重点流失区(4大片区:金沙江下游及毕节地区,嘉陵江中下游地区,嘉陵江上游陕南陇南地区和三峡库区)坡耕地549.1万hm2,占上游坡耕地总面积82.4%。据有关专家估算,长江上游重点水土流失区坡耕地年土壤侵蚀量约3.8亿 t。坡耕地土壤侵蚀不仅导致了土地生产力下降,同时由于侵蚀产生的物质颗粒细小,归槽率大,大量非点源污染物质随泥沙被搬运到下游,造成下游河道和湖泊淤积和水环境污染。 长江上游坡耕地水土流失严重,已受到广泛关注。80年代以来,有关科研单位、大专院校和水保部门在长江上游地区相继开展了一些坡耕地土壤侵蚀的试验研究,但对广大的长江上游地区来说,坡耕地土壤侵蚀的定量研究仍显得十分贫乏。1989年起实施的“长治”工程,将坡耕地作为治理工作的重点。经过多年努力,“长治”工程取得了明显的成效。但各项治理措施的效益如何,由于缺少土壤侵蚀背景资料,尚难做出定量评价。 在缺少长期水土流失观测资料的情况下,寻求一种简便快捷而又能反映坡面多年土壤侵蚀量的测试方法是非常必要的。国内外大量研究表明,137Cs法是估算坡耕地土壤侵蚀量的一种有效而且可靠的方法。该方法是利用50~70年代大气层核试验产生的137Cs人工同位素在土壤剖面中的残存量,推算自核暴炸以来坡面土壤侵蚀或沉积量。本次利用137Cs技术对长江上游重点流失区坡耕地土壤侵蚀进行定量研究,以期为上游区土壤侵0蚀评价和水土保持效益估算提供依据。 1 研究领域的总结 本次研究以长江上游重点水土流失区坡耕地为研究对象,分别在4大片区采集了6个点6种土壤类型坡耕地土样。研究区概况见表1。 2 学习方法 2.1 样品采集和取样方法 本次样品采集采用地形剖面法,即每一研究地块从坡顶到坡脚顺坡布设取样剖面线,然后沿剖面线按一定间距采集土壤剖面样品。每一地块布设2条以上取样剖面线,每一剖面线至少采集4个土壤剖面样品。 采集的土壤剖面样品,分为分层样和全样两种。分层样采用长20 cm、宽10 cm、高2 cm的矩形框取样,即一个分层样深度为2 cm。取样时将矩形框垂直打入地面,边框顶和地面齐平后,用刮刀刮取出框内土壤,即完成一个分层样。取完一层后,再将取样框垂直向下打入,直至边框顶和上一层取样后的框内地面齐平。全样大部分采用直径7.0cm的取样筒,取样时将取样筒垂直打入地面至一定深度后,拔出取样筒,取出筒内土样。部分全样取样采用直径8.5 cm的取土钻,取样深度一般为35~40 cm,有土壤堆积发生的剖面,取样深度大于40 cm,超过可能含137Cs的土层深度。 2.2 颗粒分析和重测量测定 土样经风干、研磨、过筛(孔径2.0 mm),称重后,粒径2 mm的部分供137Cs含量测试和颗粒分析,样品测重380 g左右。样品137Cs含量采用美国坎培拉公司生产的高纯锗探头多道γ能谱仪测定,测量时间300 00 s。土壤粒径分析采用移液管法。 2.3 计算土壤侵蚀的方法 2.3.1 与耕作土年流失厚度的关系 利用137Cs法计算坡耕地任意点多年土壤平均侵蚀量,一般采用下列计算公式: X=A0×(1?h/H)N?1963(1)X=A0×(1-h/Η)Ν-1963(1) 式中:X——137Cs面积浓度(Bq/m2);A0——137Cs本底值(Bq/m2);H——犁耕层深度(cm);h——土壤年流失厚度(cm);N——取样年份。 此式用于计算湿润地区(如长江上游)坡耕地的侵蚀量,往往偏大。这是由于核爆期间每年犁耕前的耕作土表层137Cs含量高于耕作土的平均值所致(137Cs表面富集)。近年来,137Cs表面富集对侵蚀量计算影响的研究取得了突破,英国研究者提出了复杂模式,但参数值很难确定,不便使用。张信宝对(1)式进行了改进,改进后的公式如下: X=Ae×(1?h/H)N?1963(2)X=Ae×(1-h/Η)Ν-1963(2) 式中:Ae——137Cs有效本底值(Bq/m2)。就坡耕地而言,此值为参与犁耕—侵蚀复合过程的137Cs本底值。 Ae和A0的大致关系为:Ae=A0×(1-R)R为径流系数 (3) 2.3.2 坡耕地土流失量的计算 利用式(2)求得的点土壤侵蚀量是犁耕作用和流水侵蚀的共同结果,不能表征流水侵蚀强度。若计算坡耕地土壤流失量,需消除犁耕作用的影响,以地块为单元进行计算,因为犁耕作用主要发生在地块内部,不将土壤运移出地块,不造成地块137Cs总量的减少。因此,根据某一地块的平均137Cs面积浓度,用(2)式可求得地块的平均侵蚀量。 3 结果与分析 3.11 基本无侵蚀土地非单次降雨期的 137Cs本底值为单位面