Japan Geoscience Union Meeting 2025

Presentation information

[E] Oral

M (Multidisciplinary and Interdisciplinary) » M-IS Intersection

[M-IS01] Particulate Gravity Current

Fri. May 30, 2025 3:30 PM - 5:00 PM Exhibition Hall Special Setting (6) (Exhibition Hall 7&8, Makuhari Messe)

convener:Hajime Naruse(Department of Geology and Mineralogy, Graduate School of Science, Kyoto University), Yuichi Sakai(Faculty of Agriculture, Utsunomiya University), Hiroyuki A. Shimizu(Sabo and Landslide Technical Center), Takahiro Tanabe(National Research Institute for Earth Science and Disaster Resilience), Chairperson:Yuichi Sakai(Faculty of Agriculture, Utsunomiya University), Hiroyuki A. Shimizu(Sabo and Landslide Technical Center), Takahiro Tanabe(National Research Institute for Earth Science and Disaster Resilience)

4:30 PM - 4:45 PM

[MIS01-04] Numerical modeling of transient water table in shallow unconfined aquifers: A hyperbolic theory and well-balanced finite volume scheme

*Ying-Hsin Wu1 (1.Kyoto University)

Keywords:subsurface flow, shallow water table, unsteady Darcy's law, well-balanced finite volume method

We present a new methodology capable of modeling transient motion of shallow phreatic surface of groundwater in unconfined aquifers. This methodology is founded on a new and comprehensive theory for water table motion and a corresponding efficient numerical scheme. In the theoretical aspect, we derived a new set of governing equations constituted by a depth-averaged continuity equation and momentum equations based on unsteady Darcy's law. The derived governing equations are of the hyperbolic type and possess stiff terms in the momentum equations due to the inertia motion in a characteristic time scale that is relatively shorter than the time scale of seepage motion. To effectively solve the derived hyperbolic system with stiff terms, in the numerical aspect, we utilize f-wave propagation algorithm, an explicit finite volume method, that can ensure numerical convergence and well-balancing solutions when momentum is rapidly relaxing to an equilibrium of steady state. Verification is successfully performed by comparing the results with analytic solutions to the classic problem of multidimensional spreading of a groundwater mound. This study demonstrates that the proposed methodology can accurately and satisfactorily simulate the spatiotemporal distribution of shallow water table and its wetting front in unconfined aquifers.