Analytical Solution for One‐Dimensional Steady‐ and Transient‐State Flow in Vertical Heterogeneous Unsaturated Soils.

Saved in:
Bibliographic Details
Title: Analytical Solution for One‐Dimensional Steady‐ and Transient‐State Flow in Vertical Heterogeneous Unsaturated Soils.
Authors: Liu, Shuang1 (AUTHOR), McCartney, John S.2 (AUTHOR), Xiao, Yang1 (AUTHOR) xiaoy@cqu.edu.cn
Source: Water Resources Research. Dec2025, Vol. 61 Issue 12, p1-22. 22p.
Subjects: One-dimensional flow, Hydraulic conductivity, Equilibrium, Waterlogging (Soils), Engineering models, Heterogeneity, Slope stability
Abstract: Exact solutions for one‐dimensional steady‐state and transient liquid flow toward a water table in heterogeneous unsaturated soils are critical in predicting saturation profiles in several real‐world applications including interpretation of climate change effects on the subsurface and impacts on slope stability. In this study, vertical heterogeneity in saturated hydraulic conductivity with depth is characterized by an exponential decay function. A steady‐state solution is derived based on Darcy's law and the water table depth, and two transient‐state solutions are obtained from Richards' equation using the Laplace transform and the modified Bessel equation under common upper boundary conditions, that is, flow rate and pressure head, following the initial steady‐state flow condition and a water table at a specified depth. The transient and steady‐state solutions are compared with numerical solutions obtained from a multi‐layered configuration and an analytical solution for homogeneous soils, demonstrating their reliability and efficacy. Various hypothetical heterogeneous soils with differing parameters are employed to assess the flow behaviors, illustrating that vertical heterogeneity impacts the pressure head profiles. More pronounced heterogeneous soils exhibit lower pressure head and effective saturation during flow compared to homogeneous soils, which depend on the air‐entry value and the water table depth, irrespective of the upper boundary conditions. The solution has been employed to predict volumetric water content profiles in the Loess Plateau of China and evaluate the effect of infiltration on the stability of a shallow landslide in Japan based on the infinite slope model and the suction stress concept. Key Points: Exact solutions are derived for 1D flow in heterogeneous unsaturated soils, considering depth‐dependent saturated hydraulic conductivityThe novel flow analytical solutions can be employed to assess shallow landslide mobilization, nitrate transport, and other surficial issuesThe analytical solution is a reliable benchmark for numerical models that handle seepage problems while considering spatial variability [ABSTRACT FROM AUTHOR]
Copyright of Water Resources Research is the property of Wiley-Blackwell and its content may not be copied or emailed to multiple sites without the copyright holder's express written permission. Additionally, content may not be used with any artificial intelligence tools or machine learning technologies. However, users may print, download, or email articles for individual use. This abstract may be abridged. No warranty is given about the accuracy of the copy. Users should refer to the original published version of the material for the full abstract. (Copyright applies to all Abstracts.)
Database: Engineering Source
Description
Abstract:Exact solutions for one‐dimensional steady‐state and transient liquid flow toward a water table in heterogeneous unsaturated soils are critical in predicting saturation profiles in several real‐world applications including interpretation of climate change effects on the subsurface and impacts on slope stability. In this study, vertical heterogeneity in saturated hydraulic conductivity with depth is characterized by an exponential decay function. A steady‐state solution is derived based on Darcy's law and the water table depth, and two transient‐state solutions are obtained from Richards' equation using the Laplace transform and the modified Bessel equation under common upper boundary conditions, that is, flow rate and pressure head, following the initial steady‐state flow condition and a water table at a specified depth. The transient and steady‐state solutions are compared with numerical solutions obtained from a multi‐layered configuration and an analytical solution for homogeneous soils, demonstrating their reliability and efficacy. Various hypothetical heterogeneous soils with differing parameters are employed to assess the flow behaviors, illustrating that vertical heterogeneity impacts the pressure head profiles. More pronounced heterogeneous soils exhibit lower pressure head and effective saturation during flow compared to homogeneous soils, which depend on the air‐entry value and the water table depth, irrespective of the upper boundary conditions. The solution has been employed to predict volumetric water content profiles in the Loess Plateau of China and evaluate the effect of infiltration on the stability of a shallow landslide in Japan based on the infinite slope model and the suction stress concept. Key Points: Exact solutions are derived for 1D flow in heterogeneous unsaturated soils, considering depth‐dependent saturated hydraulic conductivityThe novel flow analytical solutions can be employed to assess shallow landslide mobilization, nitrate transport, and other surficial issuesThe analytical solution is a reliable benchmark for numerical models that handle seepage problems while considering spatial variability [ABSTRACT FROM AUTHOR]
ISSN:00431397
DOI:10.1029/2024WR039104