[论文解读] A unified kinematic wave theory for melt infiltration into firn

Mohammad Afzal Shadab, Anja Rutishauser+2 more|arXiv (Cornell University)|Mar 24, 2024

Landslides and related hazardsEnvironmental Science被引用 3

一句话总结

该论文通过建立描述质量与焓守恒的2×2双曲型偏微分方程组，发展了一套统一的融水入渗冰晶雪冰的运动波理论，能够为12种物理上相关的黎曼问题提供解析解。该理论可捕捉融化/冻结界面、冻结锋面、蓄水层及不透水冰透镜体的形成，为冰川水文与冰盖动力学的数值模型提供了精确基准。

ABSTRACT

Motivated by the refreezing of melt water in firn we revisit the one-dimensional percolation of liquid water and non-reactive gas in porous ice. We analyze the dynamics of infiltration in the absence of capillary forces and heat conduction to understand the coupling between advective heat and mass transport in firn. In this limit, we formulate a kinematic wave theory that results in a 2X2-system of hyperbolic partial differential equations (PDEs) corresponding to the conservation of composition and enthalpy. For simple initial conditions (Riemann problems) this system admits self-similar solutions that illuminate the structure of melting/refreezing fronts and analytical solutions are provided for 12 basic cases of physical relevance encountered in the literature. Further we develop an extended kinematic theory that encompasses the cases when the firn saturates completely to form a perched water table governed by elliptic PDE so that the model is no longer fully hyperbolic (local). These solutions provide benchmarks for numerical models of melt infiltration into firn. They also provide insight into important physical processes such as the formation of frozen fringes, the perching of meltwater on pre-existing low porosity layers and the conditions required for impermeable ice lens formation. Lastly, these analytic solutions can be utilized to improve and compare the performance of the firn hydrology, ice-sheet and Earth system models. Our analysis provides a theoretical framework to understand these important processes in firn which affect the partitioning between meltwater infiltration and surface runoff and therefore determine the surface mass loss from ice sheets and its contribution to sea level rise.

研究动机与目标

• 开发一种基于物理的、解析的融水入渗冰晶雪冰的理论框架，以考虑对流热与质量输运。
• 通过在完全饱和或非均质冰晶雪冰中引入椭圆行为，弥补经典运动波理论的局限性。
• 为关键物理过程（如冻结锋面形成与不透水冰透镜体发展）提供精确的解析解。
• 为验证和改进冰川水文、冰盖及地球系统模型提供基准解。
• 增进对融水在入渗与地表径流之间分配的理解，这直接影响表面质量平衡与海平面变化。

提出的方法

• 在忽略毛细力与导热作用的前提下，基于组分与焓守恒，构建2×2双曲型偏微分方程组。
• 应用特征线法求解黎曼问题，得到12种典型入渗情景的自相似解。
• 将运动波模型扩展至包含完全饱和区域的椭圆型偏微分方程行为，从而实现对蓄水层的建模。
• 引入三相区（冰、液态水、空气），并采用无量纲化的孔隙率、饱和度与焓变量，以统一描述。
• 推导并证明特征结构（特征值与特征向量）的性质，以确立波速的双曲性与单调性。
• 利用特征分解与通量梯度分析，验证波传播与前缘演化的物理解释一致性。

Figure 1 : The dependence of temperature and volume fractions on dimensional and dimensionless enthalpy and composition, ( $C,H$ ) and ( $\mathcal{C},\mathcal{H}$ ) respectively. Dimensional $C,H$ : (a) temperature and volume fractions of (b) water, (c) ice and (d) gas phases. Dimensionless $\mathca

实验结果

研究问题

• RQ1 在无毛细力与导热效应的冷冰晶雪冰中，融水入渗过程中的融化与再冻结界面如何演化？
• RQ2 在低孔隙率冰晶雪冰层中，蓄水层形成的解析条件是什么？
• RQ3 在何种条件下，由于完全饱和与孔隙率有限，会形成不透水冰透镜体？
• RQ4 如何将运动波模型扩展，以在完全饱和的椭圆型区域保持有效性？
• RQ5 在初始条件不同的分层冰晶雪冰中，湿润前缘与冻结前缘的相互作用存在何种精确解？

主要发现

• 该模型为12种基本物理情形的融水入渗（包括湿润与冻结前缘）提供了解析解，且明确依赖于初始饱和度与焓。
• 该理论证实，由于熔化潜热较大，再冻结对前缘速度的影响较小，与先前的野外观测一致。
• 冻结锋面与蓄水层的形成被解析地捕捉，其成因是低孔隙率层的孔隙率降低与通量不平衡。
• 当体积通量超过低孔隙率层的孔隙容量时，预测将形成不透水冰透镜体，导致完全饱和与水力屏障的产生。
• 扩展模型成功实现了在饱和区域从双曲型到椭圆型行为的过渡，解决了经典运动波理论的关键局限。
• 这些解可作为精确基准，用于验证数值模型，从而降低冰川水文与地球系统模拟中的不确定性。

Figure 2 : The dimensionless flux of composition or enthalpy in $\mathcal{C}\mathcal{H}$ phase space for $m=3$ and $n=2$ . In region 1 consisting of water and gas region ( $\mathcal{H}\leq 0$ ) as well as region 3 comprising of three phase region ( $0<\mathcal{H}<\mathcal{C}$ ), the flux of dimensio