基岩河道中磨蚀与冻融风化相互作用的数值模拟

《Water》：Numerical Simulation of the Interaction Between Abrasion and Freeze–Thaw Weathering in a Bedrock Channel

　　

基岩河道演化受机械下切与季节性变率下风化过程的共同影响。研究人员采用数值模拟考察床载（bedload）冲击磨蚀与冻融风化如何在流量与温度的季节强迫下相互作用，驱动基岩河道横断面调整。一组数值试验系统性改变初始冲积层厚度、基岩抗拉强度、高流量洪峰流量与冻融侵蚀深度。结果表明沉积物覆盖对磨蚀效率具强控制作用：薄冲积覆盖增强工具效应（tool effect）并促进以中心为主的下切，较厚覆盖则逐渐抑制全床面磨蚀并稳定横断面几何形态。较高基岩抗拉强度仅削弱磨蚀驱动下切，使抗蚀性较强基质中磨蚀-风化移除面积比降低。流量对下切呈阈值型影响，超过阈值后磨蚀急剧增强并在近河心产生更深下切。冻融风化移除集中于河岸与边缘，在中等规定侵蚀深度下促进侧向改造，但在最大深度下因碎屑产量增加而受强烈抑制。上述响应界定了三种特征区制：磨蚀主导型、相互作用主导型与冻融主导型，为解释基岩河道演化中耦合季节过程提供框架。

研究背景方面，移动沉积物的机械侵蚀是基岩河道下切的基本驱动力，盐跃-磨蚀（saltation–abrasion）框架描述了床载颗粒反复冲击对基岩表面的磨蚀，其效率受岩石强度、冲击频次与能量及冲积覆盖度控制，其中工具效应与覆盖效应相互竞争。野外观测表明基岩河道演化受风化弱化与机械侵蚀交互控制，岩性特征影响基岩对磨蚀与拔蚀（plucking）的响应，风化降低强度并促进裂隙发育从而提升对床载冲击的敏感性。高纬与高山河流具高泥沙输移与水文-泥沙强耦合特征，冻融风化通过冰扩张产生拉应力渐进破坏基岩，沿水位波动致间歇出露的河岸最为有效。既有基岩下切数值模型多未显式耦合冻融风化与床载盐跃-磨蚀，常低估寒区物理风化强度及其与机械下切的交互，因此研究人员开展本研究以填补该空白。

该研究发表于《Water》，研究人员以日本北海道Pekerebetsu河上游（花岗岩基岩、冲积覆盖较薄、河岸间歇出露弱固结基岩）为参照场点，建立二维数值模型显式耦合冻融风化与床载盐跃-磨蚀，目标为将冻融纳入二维盐跃-磨蚀框架、量化风化-磨蚀交互对横断面调整的控制、评估气候驱动冻融强度变化对寒区高沙输移河道下切与地貌的影响。

关键技术方法上，研究人员采用深度平均浅水方程（含连续与动量方程）求解二维流场，床面剪切应力由流速幅值结合Manning–Strickler关系评估，冲积层厚度由泥沙守恒方程与Meyer–Peter–Müller床载能力公式及Hasegawa横向床载公式更新；基岩下切采用宏观糙度盐跃-磨蚀-淤积（MRSAA）框架，冲积覆盖分数由局地沉积厚度与基岩宏观糙度高度归一化获得，磨蚀项与独立规定的冻融移除项相加；冻融深度由一维基岩热传导模型（Izumiyama等）依据气温与太阳辐射计算孔隙水冻融循环数，超临界循环数（标定4次）岩层破碎为细粒 wash load 移出系统；水文强迫采用理想双峰季节序列（夏季洪峰35 h、冬季低流1 h×50周期），边界为周期边界维持床载再循环，参数取自Pekerebetsu河野外资料与d4PDF（+4 K情景参考）及已发表实验，开展15组敏感试验与局部单因子敏感分析。

研究结果部分，3.1节初始冲积厚度影响表明，初始冲积厚度0.70–1.20 m变化调节基岩出露与泥沙活动性：薄覆盖（0.70–0.80 m）下河心高速流致低覆盖分数，磨蚀集中河心且冻融沿边缘，形成中心聚焦下切；增厚至1.20 m时全床持久覆盖屏蔽床载冲击，磨蚀大减、下切变浅、横断面 relief 减小，冻融仍活跃于边缘。3.2节基岩抗拉强度影响显示，强度由0.20增至1.00 MPa依式削弱磨蚀系数，河心下切深度与磨蚀面积递减，冻融移除面积相对稳定，故高强度下磨蚀-风化面积比下降、冻融相对重要性上升但非风化增效。3.3节流量影响指出，洪峰200–300 m3/s时中心覆盖未破、下切限于边缘；350 m3/s中心下切显现；400 m3/s中心深下切与relief大增，磨蚀面积随流量非线性上升而冻融面积几乎不变，证实流量经降低覆盖分数而非直接增磨蚀率来强化中心磨蚀。3.4节冻融深度影响中，5 mm时整体调整有限且风化移除面积早期快增后放缓（因边坡趋平减少冬夏出露差）；37 mm为现气候参考；50 mm时风化移除面积增大，且因河宽增加使中心覆盖分数降低间接增强磨蚀，但侧向范围未同比扩大。3.5节局部敏感分析（25周期）显示宏观糙度高度η对最大下切与磨蚀面积响应最强（±0.06 m变动引38%–136%下切变化），次之基岩抗拉强度σ_t，经验常数k_ab中等，水力糙度非单调弱影响，冻融临界循环数联动ε_ft仅微调风化面积；敏感排序η > σ_t > k_ab > 糙度 > 冻融参数。

讨论部分总结，研究人员指出周期边界下初始冲积厚度作长期上游供沙代理，薄沙促中心磨蚀、厚沙屏蔽转向边缘，与既往水槽及数值研究一致；浅水方程未解三维次流致近岸平滑；基岩与冲积糙度等值假设引入不确定；直道理想几何使冻融响应近线性flattening，自然弯曲河道或异；冻融深度按规定独立于σ_t故未反映强岩抗裂真实衰减，但分离两效率凸显形态反馈（宽槽降中心覆盖分数反促磨蚀）；水流-基岩热交换未耦合可能高估冬水位附近冻深；风化碎屑即产即冲无临时护甲与再淤积，或高估边缘下切；侵蚀区制图基于15组RBF插值分磨蚀主导（R_a/w>2）、交互主导（1–2）、冻融主导（<1），供河道管理区分桥基下切、季节边缘弱化等风险；气候变暖经改雪盖与冻频使寒区两过程平衡空间分异。

结论部分翻译：本研究通过数值模拟考察床载冲击磨蚀与冻融风化交互如何在流量与温度季节变率下控制基岩河道横断面演化。结果显示河道调整受沉积覆盖、基岩强度、流量与冻融侵蚀深度共同调节，薄覆盖与高流量增强输沙与冲击频次促河心磨蚀下切，厚覆盖与高强度抑制磨蚀并将调整推至边缘；冻融主要影响出露河岸，在磨蚀受抑时相对凸显。交互产生磨蚀主导、交互主导与冻融主导三区制。结论强调过程交互对基岩地貌的控制，并为寒区季节条件影响长期演化提供新见解。