[論文レビュー] 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