  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.
