2025年水土流失动态监测技术指南

2025年水土流失动态监测技术指南以支撑国家生态保护修复、国土空间规划及"双碳"战略实施为目标，聚焦全国及重点区域水土流失状况的精准感知与趋势研判，覆盖自然地理单元、行政单元、小流域等多尺度监测场景，整合地面观测、遥感反演、模型模拟等多元技术手段，规范监测指标、数据采集、处理分析、成果应用全流程技术要求，适用于各级水土保持管理部门、科研机构及社会化服务单位开展水土流失动态监测工作。一、监测指标体系构建监测指标需兼顾科学性与可操作性，围绕水土流失成因、过程及影响构建三级指标体系。一级指标包括自然驱动因子、人为扰动因子及流失强度响应因子，二级指标细化至具体量化参数，三级指标为数据采集与计算的最小单元。自然驱动因子中，地形因子需采集坡度（0°-90°分10°间隔）、坡长（5-200m分段统计）、坡型（凸坡、凹坡、直坡、阶状坡）及海拔（按100m梯度分级）；降雨因子需涵盖年降雨总量、侵蚀性降雨频次（日降雨量≥12mm）、降雨侵蚀力R值（采用修正通用土壤流失方程R因子计算方法，逐次降雨动能与最大30分钟雨强乘积求和）；土壤因子包括质地（砂粒、粉粒、黏粒含量）、容重（1.0-1.8g/cm³）、有机质含量（0-10%分档）及土壤可蚀性K值（基于USLE-K公式计算，结合中国土壤数据库本地化修正）；植被因子需获取植被覆盖度（0-100%分5%间隔）、植被类型（乔木、灌木、草本、经济林、农作物等）、生物量（干重0-15kg/m²）及枯落物厚度（0-10cm）。人为扰动因子包含土地利用类型（耕地、园地、林地、草地、建设用地、未利用地）、耕作措施（顺坡耕作、横坡耕作、等高种植、梯田化率）、工程扰动强度（施工面积、挖填土方量、临时防护措施覆盖率）及放牧强度（载畜量0.5-3.0羊单位/hm²）。流失强度响应因子以水力侵蚀为主，兼顾风力、重力侵蚀，分级标准执行《土壤侵蚀分类分级标准》（SL190-2007），具体指标包括土壤流失模数（t/(km²·a)）、侵蚀沟密度（km/km²）、面蚀面积占比（%）、沟蚀面积占比（%）及崩岗、滑坡等重力侵蚀点密度（个/km²）。二、监测技术方法与实施（一）地面监测技术1.径流小区观测：按"标准小区+辅助小区"模式布设，标准小区尺寸20m×5m（长×宽），坡度5°-25°覆盖主要侵蚀坡度，配套自动气象站（监测降雨、温度、湿度）、径流泥沙自动采样仪（0.5小时分辨率）及人工观测桩（每月测量土壤含水量、植被覆盖度）。辅助小区针对特殊下垫面（如裸地、荒坡、梯田）设置，尺寸可调整为10m×5m或5m×5m，重点监测不同措施减沙效益。2.小流域卡口站监测：选择集水面积5-20km²的典型小流域，在出口处布设卡口站，配备超声波流量计（精度±0.5%）、激光泥沙粒度分析仪（测量范围0.1-2000μm）及视频监控系统（24小时记录洪水过程）。每场次洪水后需人工校核泥沙量，同步调查流域内耕作、工程活动变化。3.样地调查：采用分层随机抽样法，按地貌类型（山地、丘陵、平原）、土地利用类型划分抽样层，样地面积0.25-1hm²（林地、草地取大值，耕地取小值）。调查内容包括：①植被：每木检尺（胸径≥5cm乔木）、草本盖度（50cm×50cm样方5个）；②土壤：剖面开挖（1m×1m×1.5m）测量容重（环刀法）、质地（激光粒度仪）；③侵蚀特征：量测浅沟（深0.1-0.5m）、切沟（深0.5-3m）的长度、宽度、深度，计算侵蚀量（体积法）。（二）遥感监测技术1.数据源选择：采用多源卫星数据融合，包括高分一号（2m全色/8m多光谱，重访周期4天）、高分六号（2m全色/8m多光谱，红边波段）、哨兵2号（10m多光谱，覆盖全球）、terrasar-x（1mSAR，穿透云雨）及无人机（0.1-0.5m分辨率，重点区域详查）。2.信息提取流程：①辐射校正：采用FLAASH模型对光学影像进行大气校正，SAR影像采用Gamma软件完成噪声去除与地形校正；②特征提取：植被覆盖度通过NDVI（归一化植被指数）反演（公式：NDVI=(NIR-R)/(NIR+R)），坡度坡向基于DEM（30m分辨率，采用ANUDEM算法插值）计算，土地利用类型通过随机森林分类器（训练样本≥1000个/类）提取；③侵蚀强度分级：结合坡度、植被覆盖度、土地利用及降雨因子，构建决策树模型（规则：耕地+坡度＞25°+覆盖度＜30%→强烈侵蚀），人工修正误判图斑（野外验证点≥总图斑数5%）。（三）模型模拟技术优先选用修正通用土壤流失方程（RUSLE2.0）进行区域尺度模拟，参数本地化校正如下：①R因子：基于全国2100个气象站逐次降雨数据，建立R值空间插值模型（精度R²＞0.85）；②K因子：利用中国土壤数据库（25万比例尺），结合野外采样（每100km²1个样点）建立K值估算公式（K=0.2+0.3exp[-0.0256SAN(1-SIL/100)]×(SIL/CLA+SIL)^0.3×[1-0.25C/(C+exp(3.72-2.95C))]×[1-0.7(1-SAN/100)/(1-SAN/100+exp(-5.51+22.9(1-SAN/100)))]，其中SAN为砂粒含量，SIL为粉粒含量，CLA为黏粒含量，C为有机质含量）；③LS因子：采用Wischmeier-LSmith公式（LS=(λ/22.13)^m×(sinθ/0.0896)^n，m=0.5（θ＞5°）或0.4（3°-5°）或0.3（1°-3°），n=1.4）；④C因子：基于遥感反演的植被覆盖度，采用C=exp[-0.78-0.34×NDVI+0.04×NDVI²]（R²=0.89）；⑤P因子：根据水土保持措施类型赋值（梯田0.1-0.3，等高耕作0.5-0.7，无措施1.0）。模型输出土壤流失模数栅格图（30m×30m），通过地面监测点数据验证（相对误差≤15%），误差超限区域重新校准参数。三、数据采集与质量控制数据采集需建立"天-空-地"协同观测网络，地面观测数据按月上传至国家水土保持监测数据库，遥感数据按季度更新，模型模拟结果年度发布。质量控制贯穿全流程：1.仪器设备质控：径流小区自动采样仪每季度进行流量标定（误差≤2%），泥沙测量仪每月用标准泥沙样校准（误差≤3%）；无人机航摄前需检查POS系统（定位精度≤0.1m）、相机参数（像元分辨率符合设计要求）；气象站传感器每年由计量部门检定（温度±0.5℃，降雨±1%）。2.数据采集质控：外业调查需双人复核，记录观测时间、天气、仪器型号及操作人；遥感影像云量需＜10%（重点区域＜5%），无效数据（云覆盖、条带噪声）需补采；模型输入数据需进行一致性检查（如DEM与遥感影像空间配准误差≤1个像元）。3.数据存储质控：采用"本地+云端"双备份，地面观测数据存储为CSV格式（含时间、地点、指标、单位、备注字段），遥感数据存储为GeoTIFF格式（含投影信息、元数据标签），模型结果存储为GRID格式（附参数说明文档）。元数据需包含数据来源、采集时间、处理方法、质量等级（优、良、中、差），质量等级为"差"的数据需标注问题原因（如仪器故障、天气影响）并限制使用。四、数据处理与分析（一）数据预处理1.格式转换：将不同来源数据统一为WGS84/UTM投影（中国区域采用高斯-克吕格投影，6°分带），分辨率重采样至30m（重点区域10m）。2.噪声剔除：地面观测数据采用3σ准则（超出均值±3倍标准差的异常值）剔除，遥感数据通过邻域均值滤波（3×3窗口）去除椒盐噪声。3.空间配准：以DEM为基准，对遥感影像进行几何校正（地面控制点≥20个，均方根误差≤0.5个像元）。（二）时空分析1.时序变化分析：采用Mann-Kendall趋势检验（α=0.05）判断近10年水土流失强度变化趋势，计算年际变化率（(当前值-基准值)/基准值×100%）；通过Hurst指数（0.5-1.0表示持续性）预测未来演变趋势。2.空间分异分析：利用热点分析（Getis-OrdGi统计量）识别高强度侵蚀聚集区（Z得分＞1.96为显著热点），通过地形位指数（TPI=log[(海拔/区域平均海拔)×(坡度/区域平均坡度)]）分析不同地形部位的侵蚀特征。3.驱动因子识别：采用地理探测器模型（q统计量）量化各因子对侵蚀的解释力（q值越大，影响越强），重点分析植被覆盖度（q常＞0.6）、坡度（q常＞0.5）、土地利用（q常＞0.4）的交互作用（双因子增强或非线性增强）。（三）模型模拟与验证基于RUSLE模型输出年土壤流失量，结合小流域卡口站实测输沙量进行验证（纳什效率系数NSE≥0.65为可接受，NSE≥0.75为良好）。对NSE不达标区域，需排查参数误差（如K因子是否未本地化）或模型适用性（如重力侵蚀占比＞30%时需切换至WEPP模型）。五、成果表达与应用（一）成果类型1.监测包括年度水土流失状况（面积、强度分级、变化趋势）、重点区域（如长江上游、黄河中游）专题分析、典型案例（如某小流域治理前后对比），附主要指标统计表（分省、分地貌类型）。2.专题图件：包括土壤流失模数图（分级：微度＜500，轻度500-2500，中度2500-5000，强烈5000-8000，极强烈8000-15000，剧烈＞15000，单位t/(km²·a)）、侵蚀类型分布图（面蚀、沟蚀、重力侵蚀）、驱动因子贡献图（各因子q值空间分布），比例尺1:10万（全国）、1:5万（重点流域）、1:1万（小流域）。3.数据库：包含地面观测数据（2000年至今）、遥感反演产品（2010年至今，30m分辨率）、模型模拟结果（2005年至今，年度更新），支持SQL查询、空间叠加分析及API接口调用。4.预警产品：针对强降雨过程（24小时降雨量≥50mm），结合前期土壤湿度（通过微波遥感反演）、地形条件，发布小流域尺度水土流失预警（黄色预警：流失模数可能达2500-5000；橙色预警：5000-8000；红色预警：＞8000），预警时效≤6小时。（二）应用方向1.生态修复：为国土绿化、小流域综合治理提供靶区识别（优先治理剧烈侵蚀区）、措施配置建议（如25°以上坡耕地退耕还林，沟蚀区布设谷坊）。2.国土空间规划：在"三区三线"划定中，将高强度侵蚀区纳入生态保护红线，限制开发建设活动；在城镇开发边界内，要求建设项目水土流失防治责任范围全覆盖。3.政策评估：量化"十三五"以来水土保持重点工程效益（如某工程区土壤流失模数下降40%），为"十四五"政策调整提供数据支撑。4.公众服务：通过国家水土保持监测网发布年度监测公报，开发手机APP提供县域水土流失查询、预警推送服务，提升公众生态保护意识。六、质量保障与持续改进建立"管理-技术-质控"三级质量保障体系。管理层面由省级水土保持监测中心统筹，负责制定年度监测计划、协调数据共享；技术层面依托科研院所组建专家团队，解决复杂技术问题（如模型参数优化、多源数据融合）；质控层面设立独立质控小组（成员占总人数15%），对数据采集、处理、成果编制全流程抽查（抽查比例≥20%），出具质控报告（含问题清单及整改要求）。定期开展技术培训（每半年1次），内容涵盖新仪器操作（如便携式激光测距仪）、新算法应用