Abstract: | An alternative method of solution for the linearized ‘theta‐based’ form of the Richards equation of unsaturated flow is developed in two spatial dimensions. The Laplace and Fourier transformations are employed to reduce the Richards equation to an ordinary differential equation in terms of a transformed moisture content and the transform variables, s and ξ. Separate analytic solutions to the transformed equation are developed for initial states which are either in equilibrium or dis‐equilibrium. The solutions are assembled into a finite layer formulation satisfying continuity of soil suction, thereby facilitating the analysis of horizontally stratified soil profiles. Solution techniques are outlined for various boundary conditions including prescribed constant moisture content, prescribed constant flux and flux as a function of moisture change. Example solutions are compared with linearized finite element solutions. The agreement is found to be good. An adaptation of the method for treating the quasilinearized Richards equation with variable diffusivity is also described. Comparisons of quasilinear solutions with some earlier semi‐analytical, finite element and finite difference results are also favourable. Copyright © 2001 John Wiley & Sons, Ltd. |