基于水文年分层的泥沙评级曲线（SRC）区域化框架：以湄公河最大支流蒙河流域为例

《Journal of Hydrology: Regional Studies》：A hydrological-year stratified regionalization framework for sediment rating curves in the largest tributary of the Mekong River

　　

研究人员开发了一种面向数据稀缺流域（data-scarce catchments）的先进泥沙评级曲线（sediment rating curves, SRC）区域化框架。该框架将水文年分层（hydrological year stratification）与回归（regression）和相似性（similarity）方法相集成。在蒙河流域（Mun River Basin）的17个流域（catchments）中实现了模型性能提升。该框架为湄公河地区（Mekong region）的泥沙负荷（sediment load）估算提供了可移植工具（transferable tool）。

研究背景与问题：河流泥沙输移控制河床形态、水生栖息地稳定性和水资源可持续性。湄公河等大流域的泥沙通量关系三角洲形成、营养循环和农业生产力，但水电开发、土地利用变化和水资源调节改变泥沙连续性，使下游及三角洲泥沙输移下降。精确估算泥沙负荷对泥沙收支、水库淤积和流域及区域管理均十分重要。然而，许多流域泥沙资料稀缺，连续监测受后勤和经费限制，记录稀疏或不连续，限制了泥沙评级曲线（sediment rating curve, SRC；Q_s=aQ^b）的发展。SRC以经验幂律联系径流与悬移质泥沙负荷（suspended sediment load, Q_s），可在缺少直接测量处间接估算泥沙通量。SRC预测性能依赖数据覆盖和时间代表性；泥沙—流量关系具有非平稳性，受季风降雨、极端流量、泥沙供给和土地利用变化影响，多年单一SRC可能不能代表不同水文条件，在季风主导流域和数据稀缺流域尤为明显。数据稀缺时需通过物理或统计相似性把SRC参数从有资料流域转移至无资料流域。Heng and Suetsugi 在Lower Mekong Basin发现回归与流域相似性组合可提供可靠估计，但其框架依赖全时段数据，未显式考虑年际水文变异。研究空白在于：其SRC区域化方法未在蒙河流域系统评估；既有框架通常假设参数时间恒定，未联合考虑年际水文变异与空间流域异质性；何种年分类最能代表泥沙—流量关系、流域属性如何指导回归或相似性参数转移仍不清楚。为此，研究人员评估Heng and Suetsugi 方法在蒙河流域的性能，并开发与评估水文年分层区域化框架，以改善数据稀缺流域泥沙负荷估算；比较全时段、土地利用和流量分类，建立洪水、正常和干旱年区域回归方程，识别决定回归或流域相似性参数转移适用性的流域属性。其科学贡献在于联合考虑时间水文变异与空间流域异质性；蒙河流域作为代表性数据稀缺且水文变化大的湄公河支流，但更广泛应用前仍需验证。

关键技术方法：样本为蒙河流域17个流量—悬移质泥沙配对站点（2003—2022；流量/泥沙RID，降水TMD，土地利用SERVIR-Mekong，土壤FAO HWSD；缺测流量由分段校准水文模型重建）。研究人员以SRC Q_s=aQ^b幂函数/对数线性拟合为基础，比较流域相似性（CSI最小供体，单/多供体输出与参数平均）、回归（属性筛选、线性/指数/多项式）、空间邻近（质心距离与IDW）及其输出平均、参数平均和聚类阈值组合。年分层比较全时段、土地利用分类和流量分类；流量分类采用均值±1标准差统计指标及基于Gumbel、Log-Pearson III、Normal、Log-Normal拟合分布的2年/10年重现期频率分析。以留一交叉验证和PBIAS、NSE、RSR评价，满意阈值为RSR<0.70、NSE>0.50、|PBIAS|<55%。

研究结果：
3.1 既有区域化方法在蒙河流域的应用：研究人员按Heng and Suetsugi 框架，以留一交叉验证将每个站点依次视为无资料站点，比较流域相似性、回归、空间邻近及平均/组合集成。单供体流域相似性最有效（10/17站满意），回归居中（9/17），空间邻近最差（7/17）；多供体加权不及单供体。回归策略2优于策略1，C_USLE、K_USLE×C_USLE、Sc、Sr和Q_a与SRC参数相关较强，AC和Pr较弱。流域相似性与回归组合最佳（12/17站满意）。以Sc=6.0%和K_USLE×C_USLE=0.05为阈值：Sc<6.0%且K_USLE×C_USLE>0.05多用回归，Sc>6.0%且K_USLE×C_USLE<0.05多用流域相似性。全时段区域化仅约70%站点满意，提示需年分层。

3.2 年分层方法评估：研究人员比较全时段、土地利用分类、统计指标流量分类和频率分析流量分类，并测试参数拟合、输出平均、参数平均。参数拟合策略变异性最低且总体优于平均策略。全时段参数拟合14/17站满意，但PBIAS中位-14.3%倾向低估；土地利用分类14/17，NSE中位0.77；统计指标流量分类12/17且部分站点负NSE；频率分析流量分类最佳，16/17站满意，中位RSR 0.25、NSE 0.73、PBIAS -2.1%。土地利用变化（2010、2017关键年）能反映下垫面变化，但流域尺度上年际水文变异对泥沙—流量关系影响更强；山区站点土地利用分类可能更好。E.22B因上游拦沙、水库调节和非平稳表现不佳。

3.3 通过水文年分层推进区域化：研究人员按2年/10年重现期把每年分为干旱、正常、洪水年，分别建立回归方程（参数a和b为Q_a、Sc、Sr、C_USLE×K_USLE的函数），并按前述阈值在缓坡流域用回归、陡坡流域用流域相似性。与Heng and Suetsugi原始方法比较，满意站点由12增至16；与Inseeyong et al.标准化SRC方法同为16站，但在E70、E87、E91、M2A、M4、M9、M104、M176、M184、M185等站更优；E.22B略降。时间序列显示分层框架更好再现峰值泥沙负荷和时序变化。

讨论部分总结：既有区域化的主要限制不仅是空间资料稀缺，还有泥沙—流量关系的非平稳性；全时段混合洪、平、枯观测，频率分类以站点特定重现期阈值减少混合，提高水文一致性和降低偏差。框架流程为流量分类→以Sc和K_USLE×C_USLE选择回归或相似性→按水文类估计SRC参数→应用于对应流量期；适用于资料有限流域，但需在蒙河以外验证。方法适用性与泥沙产生和输移过程有关，而非仅由地理邻近决定；低坡流域由流域尺度属性控制，回归可用；陡坡流域侵蚀和输移空间异质，相似性转移更稳健。区域回归限于蒙河条件（0.126USLE<0.214，5.8%
研究结论翻译：本研究通过首先评估既有SRC区域化方法、继而开发并评估蒙河流域水文年分层框架，实现了两个目标。回归与流域相似性的组合在既有方法中性能最佳；分层框架通过区分洪水、正常和干旱条件进一步改善预测性能。在17个有资料流域中，先进框架取得更高NSE、更低RSR和更小PBIAS绝对值，16个站满足满意性能标准。这些改进证实，将数据分为洪水、正常和干旱年等不同水文类别可改善SRC参数化，更好表征泥沙—流量关系的时间变异性。研究还识别了湄公河地区不同流域类型的最佳参数化策略：Sc<6.0%且K_USLE×C_USLE>0.05的流域最适合回归方法，Sc>6.0%且K_USLE×C_USLE<0.05的流域使用流域相似性表现更好。此外，为每种水文条件建立了区域回归模型，可直接用于无资料地点。尽管有效，这些区域模型受蒙河流域地形和水文条件限制；将方法转移到物理特征不同流域，即使在同一湄公河地区，也可能需要重校准属性阈值或开发适合当地条件的新回归模型。尽管如此，该框架为监测不足地区泥沙估算提供了实用且可调整的方法；其整合空间和时间动态的能力为流域管理中的泥沙预测提供了更可靠基础。