杜朝正,杨勤科,王春梅,庞国伟

(1 西北大学 城市与环境学院, 陕西 西安 710127;2 菏泽学院 城市建设学院, 山东 菏泽 274015)

土壤侵蚀是全球性生态环境问题之一[1]，被侵蚀的土壤有机碳是温室气体碳成分的重要来源之一[2]，土壤侵蚀速率大小直接影响水土流失区域的土壤侵蚀经济损失状况[3]。治理土壤侵蚀必须以土壤侵蚀定量评价和制图为基础[4-6]，即通常采用土壤侵蚀模型定量计算土壤侵蚀速率[7]。逐步建立的土壤侵蚀模型有美国的WEPP模型[8]、USLE[9]和RUSLE[10]，欧洲的EROSEM模型[11]，荷兰的LISEM模型[12]以及中国的CSLE[13]等。国外较多研究采用USLE/RUSLE计算区域土壤水蚀速率，如Panagos等[14]编制欧盟土壤流失图，Teng等[15]编制澳大利亚土壤侵蚀图。国内学者多采用CSLE计算区域土壤水蚀速率，如章文波等[16]采用CSLE预测坡面和田间土壤侵蚀，程琳等[17]采用CSLE对陕西省土壤侵蚀进行定量评价。

一般进行区域土壤侵蚀定量评价方法可归纳为2种：第一种是基于侵蚀因子图层整体运算的方法，可称为地图代数法，代表性研究有程琳等[17]、Borrelli等[18]、Yang等[19]、周来等[20]；第二种是基于抽样调查结果和插值的方法，在研究区布设基本抽样调查单元，结合土壤侵蚀模型计算抽样单元土壤侵蚀速率，使用空间插值模型来完成区域土壤侵蚀评价制图，简称为空间插值法，代表性研究有刘宝元等[21]、Yin等[22]。地图代数法和空间插值法均已被广泛应用，但这2种方法有何异同还鲜有报道。

本研究以青藏高原典型样区为例，以现有土壤侵蚀因子数据和抽样调查单元数据为数据源，分别采用基于侵蚀因子栅格图层的地图代数法和基于抽样调查单元的空间插值法，计算研究区土壤水蚀速率，并对计算结果进行对比分析，以期为土壤侵蚀定量评价与制图提供方法参考。……