Forward and inverse modeling of water flow in unsaturated soils with discontinuous hydraulic conductivities using physics-informed neural networks with domain decomposition

TL;DR

This paper explores the use of physics-informed neural networks (PINNs) with domain decomposition to model water flow in unsaturated soils, demonstrating comparable accuracy to traditional methods and addressing challenges like boundary condition estimation.

Contribution

The study introduces a PINNs approach with domain decomposition for modeling water flow in layered soils with discontinuous hydraulic conductivities, including boundary condition estimation from sparse data.

Findings

PINNs solutions are comparable to traditional numerical methods.

PINNs can estimate unknown boundary conditions from sparse measurements.

PINNs are sensitive to initialization and slower than traditional methods.

Abstract

Modeling water flow in unsaturated soils is vital for describing various hydrological and ecological phenomena. Soil water dynamics is described by well-established physical laws (Richardson-Richards equation (RRE)). Solving the RRE is difficult due to the inherent non-linearity of the processes, and various numerical methods have been proposed to solve the issue. However, applying the methods to practical situations is challenging because they require well-defined initial and boundary conditions. Here, we present a physics-informed neural networks (PINNs) method that approximates the solution to the RRE using neural networks while concurrently matching available soil moisture data. Although the ability of PINNs to solve partial differential equations, including the RRE, has been demonstrated previously, its potential applications and limitations are not fully known. This study conducted a comprehensive analysis of PINNs and carefully tested the accuracy of the solutions by comparing them with analytical solutions and accepted traditional numerical solutions. We demonstrated that the solutions by PINNs with adaptive activation functions are comparable with those by traditional methods. We showed that soil moisture dynamics in layered soils with discontinuous hydraulic conductivities are correctly simulated by PINNs with domain decomposition. We demonstrated that the unspecified upper boundary condition can be estimated from sparse soil moisture measurements. Nevertheless, there remain challenges that require further development. Chiefly, PINNs are sensitive to the initialization of NNs and are significantly slower than traditional numerical methods.