基于遥感和GIS的密云县土壤侵蚀时空变化分析.pdfVIP

基于遥感和GIS的密云县土壤侵蚀 时空变化分析 张有全∞宫辉力@赵文吉。李小娟。 (1．中国科学院东北地理与农业生态研究所，长春130012：2．中国科学院研究生 院，北京100039；3．资源环境与地理信息系统北京市重点实验室，北京100037) 摘要：将地理信息系统与通用土壤流失力程(USLE)相结合进行密云县土壤侵蚀量的预测。在全 面调查，收集研究区气候、地形、土壤等基础资料的基础上，通过TM遥感影像提取土地利用信息、 植被覆盖度信息，并建立土壤侵蚀空间数据库。运用通用土壤侵蚀模型USLE定量计算出土壤侵 蚀量，对密云水库周边区域10年来的土壤侵蚀强度的时空变化进行分析．为该区域水土保持和北 京市的安仝供水提供科学依据。 关键词：土壤侵蚀；USLE；GIS；遥感 0引言 土壤侵蚀是指地球表面的土壤及其母质受水力、风力、冻融、重力等外力的作用，在各 种自然因素和人为因素的影响下发生的各种破坏、分离、搬运和沉积的现象川。土壤侵蚀不 仅破坏土地资源，而且将大量携带养分、重金属和化肥、农药的泥沙进入江河湖库，造成水 体富营养化”J，同时导致土壤肥力下降，生态环境恶化，是世界E的主要环境问题之一。 密云县境内的密云水库是北京市最大的饮用水源基地，保障着首都1／3的供水。由于近 年人们对其周边土地不合理的开发利用，导致了库区周围尤其是库东北、北部地区的水土流 失严重。据1999年北京市遥感调查，密云县水土流失总面积1000．282km2，占全县总面积的 45％[31随着北京市经济的发展，用水量增加，加之近几年的持续干旱，北京市水资源供 求矛盾日趋紧张。进行密云水库周边10年来的土壤侵蚀强度的时空变化分析，对保护和 改善水库周边生态环境，保证首都人民用水以及用水安全，以及合理规划水库周边土地利用有 着重要意义。 1研究区概况 风区，地处暖温带，属大陆性半湿润半干旱气候。年均降水量660mm。地势自北向西南倾 斜，东西两侧高，中部为断陷燕落盆地，西南为山前冲积平原。山地占全县面积的46．7％； 丘陵占36．6％；平原占8．3％；水库水域占8．4％，其中密云水库占7．9％。耕地主要分布在 西南部，山地以林地利用为主，中部为密云水库。 2研究方法与技术路线 在广泛收集研究气象资料、地形资料、土壤资料和有关水土保持研究 成果的基础上，采用卫星遥感、地理信息系统和野外实地调查相结合的方法，完成研究区土 地利用分布图、植被覆盖图等；依据通用土壤侵蚀方程USLE及国内众多学者提出的修正意 见和参数测定规范，结合研究区的实际情况及野外调查数据，将方程中各因子的算法进行本地 化处理，建立密云县土壤水蚀模型。以30m大小的像元为单位，利用地理信息系统的空问分 析等功能，对收集、提取的信息进行栅格运算，根据侵蚀强度分级标准对运算结果进行分级， ·86· 对密云水库周边区域10年来的土壤侵蚀强度的时空变化进行了分析，在时问序列上确定 研究区域内水土流失强度在空间和等级的变化情况。确定相应的保护措施，技术路线如 图1所示。 巨垂垂堕圃 图1土壤侵蚀强度预测流程图 3土壤水蚀预报模型 计算方法，结合研究地区的实际情况，确定各因子计算方法，并进行本地化赋值。 3．1降雨侵蚀力因子(向 雨滴击溅作用和因降雨产生的径流，是最主要的土壤侵蚀动力。降雨侵蚀力的经典计算 方法： 豆=专摅㈨] ㈤ 式中：R一多年平均年降雨侵蚀力，单位MJ·mm(hm2．h)庐l，2，…是计算平均值使用的年 数。降雨侵蚀力尺的简易计算方法： R=O．207(P*160／100)“…(2) P为年降雨量，‰为年最大60min降雨量(mm)。本文采用年R值的简易算式，首先， 对密云水文观测站测得的1990--2000年观测点数据‰、以及多年降雨等值线进行插值，内 插生成研究区域P、，神栅格图，利用GIS栅格运算功能按照公式2计算该流域降雨侵蚀力 栅格分布图。 3．2土壤可蚀性因子(∞ 土壤可蚀性因子(目为标准小区上单位降雨侵蚀力引起的土壤流失率。密云水库所在的 密云县土壤共分为三大类：棕壤、褐土、潮土。根据土壤各因子指标值查诺谟图，并参照我 国学