《European Journal of Soil Science》：Flood Modelling Evaluation in a Semiarid Catchment: The Role of a Soil Mineralogical Index From Satellite Imagery
编辑推荐：
摘要：地中海地区以偶尔发生的强降雨事件为特征。这些事件与岩性、地貌和土地利用因素共同影响流域的水沙动力学，可能引发山洪。在本研究中，研究人员对一次切断低压事件进行了水文和水力学分析。该事件发生在2019年9月11日至13日，影响了Rambla del Beal
摘要：地中海地区以偶尔发生的强降雨事件为特征。这些事件与岩性、地貌和土地利用因素共同影响流域的水沙动力学，可能引发山洪。在本研究中，研究人员对一次切断低压事件进行了水文和水力学分析。该事件发生在2019年9月11日至13日，影响了Rambla del Beal——一条位于西班牙东南部半干旱气候条件下的间歇性河流，该地区无任何流量数据可用。该流域上游以数十年积累的未经处理的采矿废物为特征。研究人员采用HEC-HMS模型进行水文分析，并应用HEC-RAS水力学模型模拟河道周围洪泛区。记录到总降雨量为226?mm，主要集中在6小时内，最大小时降雨量为76.71?mm。流量过程线显示峰值流量为171.73?m3/s。研究人员利用洪水中沉积的沉积物图评估了结果。该图通过Sentinel-2卫星影像的波段指数检测铁氧化物和氢氧化物生成。该指数被应用于洪水前和洪水后的影像。然后计算两个日期之间的指数差值。结果显示，在洪水期间发生沉积的区域呈现正值。因此，这些沉积物指示了最小淹没范围。此外，沉积物的空间分布反映了地质图上标注的冲积扇动力学的水文行为，尽管缺乏细节。通过历史航拍照片分析，确认了冲积扇沉积其沉积物的确切位置。由于农业和城市发展改变了地形，这些特征目前已不可见。总之，铁氧化物和氢氧化物矿物制图在验证受酸性矿山排水影响流域的水力学模型结果方面被证明是有价值的。
**论文解读文章**
**研究背景与目的**
地中海地区是全球变化最脆弱的区域之一，气候变化预测显示该地区极端事件（如强降雨）频率和强度将增加，进而导致洪水事件的增多。然而，大多数地中海河流为非永久性河流（intermittent rivers and ephemeral streams），其流量具有间歇性，且长期缺乏监测数据，使得对其水文响应的理解十分困难。西班牙东南部广泛分布着这类称为“ramblas”的间歇性河流，其洪水风险极高，且缺乏历史流量数据，传统水文模型难以应用。遥感技术为洪水建模提供了替代数据源，但光学卫星受云覆盖限制，而合成孔径雷达（SAR）影像在窄河段中分辨率不足且需专业处理。本研究旨在评估利用Sentinel-2卫星影像的土壤矿物学指数（铁氧化物/氢氧化物）来验证半干旱流域洪水建模结果的可行性，以解决缺乏监测数据时模型验证的难题。论文发表在《European Journal of Soil Science》。
**关键技术方法**
研究人员采用HEC-HMS 4.12水文模型进行降雨-径流模拟，基于邻近气象站（西班牙国家气象局AEMET）的日降雨数据，通过对比分析生成小时降雨过程线。利用Curve Number（CN）方法（基于土地利用、土壤水文组、坡度等，数据来源：西班牙国家地理研究所CNIG）估算净雨和洪水过程线。随后使用HEC-RAS v6.1二维水力学模型，基于2米分辨率数字高程模型（DEM，源于LiDAR点云）和曼宁糙率系数，模拟洪水淹没范围。模型验证采用两种方法：一是1956年历史航拍照片（美国军事飞行B系列）识别冲积扇地貌；二是基于Sentinel-2多光谱影像（波段B4红、B3绿）的比值指数（B4/B3）计算洪水前后铁氧化物/氢氧化物变化，以识别沉积区域。样本来源：2019年9月11-13日切断低压事件，研究区为Rambla del Beal流域（7.6 km2，平均坡度1.9%）。
**研究结果**
**3.1 水文结果**
**3.1.1 降雨估算**：通过相邻站点的降雨数据，得到2019年9月11日09:00至13日05:00累计降雨量226 mm，最大小时雨强出现在9月13日03:00（55.5 mm/h）和04:00（76.7 mm/h），降雨主要集中在12小时窗口内，前期降雨使土壤饱和，利于后期径流产生。
**3.1.2 洪水过程线模拟**：流域平均CN值为72.8，初始抽象量18.99 mm，滞时22.19分钟，响应迅速。上游峰值流量于9月13日05:00达到171.73 m3/s，河流在当日10:00后干涸；河口处峰值流量于07:00达到98.8 m3/s。
**3.2 洪水图**
**3.2.1 HEC-RAS模拟**：最大洪水范围出现在9月13日06:00，上游两城镇Beal和Llano del Beal受影响，水深多介于20-50 cm，但在两座公路桥处水深分别达9 m和4.7 m，导致水流侧向漫溢至冲积扇。下游左岸两处池塘（高约4 m）造成阻塞，水流向右岸农业用地扩散，水深超1 m。河口盐沼区域积水，洪水持续时间超过50小时。
**3.2.2 模型评估**：将1956年航拍照片与模拟洪水图对比，发现历史冲积扇和沉积痕迹与2019年洪水区域吻合，证实了人类活动（如修建池塘、改变河道）掩盖的地貌动力学。铁氧化物/氢氧化物指数差值图显示，正值区域（沉积）位于低洼、浅水、短时积水区，与模拟洪水边界高度一致，表明模型准确反映了2019年9月的淹没范围。
**讨论与结论**
卫星影像衍生的铁氧化物/氢氧化物指数与模拟洪水范围的高空间一致性，验证了该方法在酸性矿山排水（Acid Mine Drainage, AMD）活跃流域的有效性。研究区矿山废物集中于上游，强降雨期间被洪水搬运至下游低洼处（冲积扇和盐沼），导致B4/B3比值上升，指示最小淹没范围。研究还发现，HEC-RAS二维模型成功模拟了复杂冲积扇动力学，即使地形已被人类活动改变。但存在局限性：植被会干扰矿物光谱信号，需通过归一化植被指数（NDVI）掩膜；邻近流域洪水重叠导致河口沉积物来源难以区分；该方法仅反映最小洪水范围，非所有淹没区均有沉积。结论：基于Sentinel-2影像的沉积物分析是评估HEC-RAS模型洪水模拟结果的有效方法，尤其在缺乏地面监测数据的间歇性河流中。未来需进一步研究植被覆盖区域铁氧化物升高的其他来源。
