微地形调控下径流流态转换过程中的表面分形异质性

《Frontiers in Ecology and Evolution》：Surface fractal heterogeneity under microtopography during flow regime transition

　　

尽管已有研究探讨了微地形与水动力学的相互作用，但多数未能量化由微地形和持续降雨触发的流态转换过程中表面分形异质性的连续动态响应。研究人员调查了由空间异质性和水动力学诱导的土壤表面微结构，重点关注由平作（TFC）、横垄（THC）、人工挖穴（TAD）和锄耕（THE）四种地表粗糙度（SR）在可变降雨时长（RD）下的表现。结果表明，径流流速（V = 0.16–0.32 m/s）随RD先增后稳，与水力剪切应力（τ = 1.12–1.68 Pa）类似，而其随SR增大而减小，与τ不同。雷诺数（Re = 184–480）随RD准线性增大，弗劳德数（Fr = 1.3–2.7）亦如此，达西–魏斯巴赫摩擦系数（f = 0–2）则呈反向变化。表面微结构受流态转换主导。微团聚体含量与分形维数（D）与环境因子显著正相关（p < 0.05），大团聚体则相反。随机森林（RF）分析与结构方程模型（SEM）证实RD（0.36）和容重（BD，0.44）对D有正向效应，而Fr（?0.15）、f（?0.24*）和结皮厚度（TN，?0.14）为负向，模型拟合可接受（RMSEA、SRMR < 0.08，CFI、TLI > 0.9）。微地形控制下触发流态转换的阈值降雨时长可用于量化表面分形异质性。所建立的多因子因果框架为优化坡耕地抗蚀耕作实践提供了针对性理论支持。

研究背景与立项依据

坡耕地水蚀研究中，微地形（microtopography）作为基础地貌因子，通过重塑地表起伏调制溅蚀、片蚀与沉积过程。已有文献证实微地形与水动力耦合可调控产沙与颗粒分选，但多关注静态关联，忽视了在长历时降雨下由层流向量流过渡时表面分形异质性（surface fractal heterogeneity）的动态演化。分形维数（D，fractal dimension）可定量描述团聚体稳定性与沉积结皮（depositional crust）微结构，而流态转换（flow regime transition）由雷诺数（Re）、弗劳德数（Fr）及达西–魏斯巴赫摩擦系数（f）刻画。现有研究缺口在于：未将耕作微地形、水动力相变与结皮分形特征纳入统一因果路径。为此，研究人员以黄土高原壤土为对象，在模拟降雨下追踪四种耕作微地形中流态转换与表面分形异质性的连续响应，旨在为抗蚀耕作优化提供统计因果证据。

论文发表于《Frontiers in Ecology and Evolution》，研究意义在于首次将微地形控制的流态转换阈值、沉积结皮分形维数及多变量因果框架（RF+SEM）联结，弥补了动态微观结构演化机制的认识不足。

主要技术方法概括

研究人员采集陕西黄土高原耕层0–20 cm粉壤土，填筑5°坡度2.0 m×1.5 m径流槽，设平作、横垄、人工挖穴、锄耕四种微地形（n=4），以摆动喷头模拟60 mm/h降雨，分5、10、15、20、25、30 min六时段取样（每处理3重复）。微地形采用运动恢复结构（SfM）摄影测量生成0.5 mm分辨率DEM，计算表面粗糙度（SR）与凹陷蓄量（SDS）；水动力学用染色示踪测流速V，按τ=ρghs算剪切应力，由Re=Vh/ν_m、Fr=V/√(gh)、f=8ghJ/V²求无量纲指标；沉积结皮剥取后以蜡封法测容重（BD）与厚度（TN），湿筛得团聚体分级，由质量–粒径关系求D；统计上用单因素ANOVA、相关分析、随机森林回归（十折交叉验证、网格调参）筛选变量，再以最大似然SEM（216观测、19自由参数）分解RD→Fr→f→D等协变路径，以RMSEA/SRMR<0.08、CFI/TLI>0.9判拟合。

研究结果

3.1 径流动力学特征：随RD由5延至30 min，V（0.16–0.32 m/s）与τ（1.12–1.68 Pa）先升后稳，T_FC流速最高（0.28 m/s）、T_HC最低（0.19 m/s），与SR排序反比；Re（184–480）与Fr（1.3–2.7）准线性上升，f（0–2）振荡，表明层流向过渡湍流转换，平作面更早进入超临界流。相关分析示RD与Re（0.74*）、Fr（0.57）正关，与f（?0.72）负关。

3.2 微地形介导的沉积结皮异质性

3.2.1 结构属性：BD随RD增至T_HC最高1.68 g/cm³，T_FC最低1.45 g/cm³；TN在T_AD（2.7 mm）与T_HC（2.6 mm）波动上升，平作与锄耕略降。SDS以T_HC（9.36×10^-2 m³）最大，决定凹处滞沙。

3.2.2 团聚体组成：>0.25 mm大团聚体占比随RD下降，T_HC降幅60%最大，T_FC仅6%；<0.053 mm微团聚体随RD、SR、SDS正相关（p<0.05），大团聚体与SR、SDS负相关。

3.3 沉积过程诱发的表面分形动态：D在2.22–2.71间随RD先升后稳，增幅11%–16%，T_HC达2.71；D与微团聚体、BD极正关，与大团聚体、f极负关（p<0.001）。

3.4 水动力驱动的分形调控机制：RF中f与Fr联合重要性83%，RD居首；SEM示RD（0.36）、BD（0.44）正向直效于D，f（?0.24*）、Fr（?0.15）、TN（?0.14）负向；RD经Fr负向连f形成中介链，模型拟合达标。

讨论部分总结

微地形通过凹蓄延长汇流、糙度生涡降速，使T_HC30 min内无连续出流；Re/Fr在10–15 min跨临界，证实微地形尺度依赖地干预相变。结皮分形受“降雨时长—流态—摩擦”链控：长RD增Fr削弱黏性摩擦、降f，间接促D升，抵消黏滞抑制。局限为单一雨强（60 mm/h）、均质黄土、5°坡，未涵盖间歇雨与多质地。

结论翻译

本研究探明微地形扰动与降雨时长对四种耕作下沉积结皮分形异质性的交互效应，确定15 min降雨为层流向过渡流转换阈值；SR与SDS及RD决定BD、TN并提升<0.053 mm微团聚体比例；D与 environ 因子（p<0.05）及水动力参数（p<0.01）显著正关；RF与SEM示RD（β=0.36）、BD（β=0.44）促分形发展，Fr（β=?0.15）、f（β=?0.24*）、TN（β=?0.14）抑之。层流–湍流转换驱动的表面分形复杂性量化了微地形缓冲与异质性权衡，为多因子因果框架设计抗蚀耕作与精准保土提供依据。