基于遥感及气象数据土壤侵蚀敏感性评价解决方案.pdfVIP

基于遥感及气象数据的土壤侵蚀敏感性评价的解决方案 土壤侵蚀是指土壤或成土母质在外力（水、风）作用下被破坏剥蚀、搬运和沉积的过程。 土壤侵蚀是全球环境的灾难，他不仅影响了人们的生存和社会发展，更限制了全球经济的可 持续发展。三江源地区由于特定的自然环境和脆弱的生态系统，其水土流失问题的研究具有 重要意义。土壤侵蚀敏感性评价是三江源地区生态环境评价的重要部分，它能够为该地区生 态保护和建设提供科学依据。 本方案通过对三江源地区土壤侵蚀敏感性的评价来识别容易形成土壤侵蚀的区域并进 行分级。评价主要基于三江源地区自1990 年以来生态环境数据积累，包括用于土壤侵蚀敏 感性评价的各项自然 指标数据。在对该地区的生态状况进行全面调查的基础上，制定了较 为合理的技术路线。根据三江源地区土壤侵蚀类型和特点，本次评价将分水蚀区域和风蚀区 域分别进行评价，然后进行综合制图。 1、数据准备及预处理 土壤水蚀的影响因素有降雨侵蚀力（R）、土壤质地、土壤可蚀性、地形起伏度和植被 状况。其中降雨侵蚀力和植被状况指标随时间发生变化，分别有1990、2004 和 2009 年 三期数据。 土壤风蚀敏感性的指标主要有土壤质地、土壤可蚀性、地形起伏度、植被状况、风场强 度和土壤表层湿度，对各因子分别评价，然后进行综合。本次评价中植被状况、风场强度和 土壤表层湿度三项指标随时间而发生变化，分别为1990、2004 和2009 年三期的数据。 2、研究方法 2.1 地形起伏度 利用三江源地区100m 分辨率dem 数据，采用ArcMap 中的空间分析(SpatialAnalyst) 模块中的邻域分析(Neighborhood Statistic)工具，以n×n 像元的圆形为模板算子，对整个研 究区进行移动计算, 先计算出n×n 像元内的格网最大值maximum，然后计算出其领域最小 值minimum，再利用模块中的栅格计算工具(Raster Calculator) 计算最大值与最小值高程 差,就得到了该n×n 窗口的地势起伏度结果值。可表示为如下公式: 式中: RF 是分析窗口内的地势起伏度; Hmax 是分析窗口内的最大高程值; Hmin 是分析窗 口内的最小高程值。 2.2 植被指标 利用1990、2004 和2009 三期TM 遥感卫星图像，结合野外调查和人工目视解译， 得到了三江源地区1990 年以来三期土地覆被/利用和草地退化状况数据。然后融合三江源 地区1:100 万植被覆盖，我们得到了该地区土壤侵蚀敏感性评价的植被状况指标定级。 2.3 风场强度 风场强度的计算采用美国农业部土壤风力侵蚀方程 RWEQ 中的相应公式： 3 3 式中，W 为风能强度因子，单位为 m /s ，U 为离地面 2m 高处的风速（10min 最大 风速），U 是 2m 高处的临界风速，一般设置为5m/s 。计算结果用 ANUSPLINE 方法进 c 行插值得到三江源地区风场强度空间分布。然后参考全国风场强度大小，制定了土壤侵蚀敏 感性评价指的风场强度指标定级值域，对三江源地区的风场强度进行分级。 2.4 土壤表层湿度 采用美国AMSR-E 数据源，经处理得到 25km 分辨率的土壤表层湿度数据，然后插值 成1km 分辨率。根据2009 年全国土壤湿度在 3.69~23.23 范围内，制定了土壤湿度指标 分级值域。 2.5 土壤可蚀性（k）和土壤质地 土壤可蚀性的计算公式为： 其中S 为土壤结构性指数，P 为土壤可渗透性指数，其值在计算过程中根据下表判定， S 和P 的取值，对应不同的土壤类型，取值可能略有差异。 其中 式中 P 为粘粒含量，P 为粉粒含量，P 为细砂粒含量。土壤质地（砂砾、粉粒、粘 1 2 3 粒）和有机质数据来自于邵景安博士2007 年在三江源课题组的工作，依据环保部2003 年 制定的《生态功能区划 暂行规程》，我们确定了土壤质地指标的敏感性分级。 有学者对23 种主要土壤的K 值作了直接的测定，其范围为0.03~0.69，基于此我们对 土壤可蚀性指标进行了敏感性分级。 2.6 降雨侵蚀力（R） 利用日降雨资料，以半月为时段，建立利用日雨量计算半月侵蚀力的