Appropriate vertical discretization of Richards' equation for two-dimensional watershed-scale modelling

A number of watershed-scale hydrological models include Richards' equation (RE) solutions, but the literature is sparse on information as to the appropriate application of RE at the watershed scale. In most published applications of RE in distributed watershed-scale hydrological modelling, coarse vertical resolutions are used to decrease the computational burden. Compared to point- or field-scale studies, application at the watershed scale is complicated by diverse runoff production mechanisms, groundwater effects on runoff production, rumon phenomena and heterogeneous watershed characteristics. An essential element of the numerical solution of RE is that the solution converges as the spatial resolution increases. Spatial convergence studies can be used to identify the proper resolution that accurately describes the solution with maximum computational efficiency, when using physically realistic parameter values. In this study, spatial convergence studies are conducted using the two-dimensional, distributed-parameter, gridded surface subsurface hydrological analysis (GSSHA) model, which solves RE to simulate vadose zone fluxes. Tests to determine if the required discretization is strongly a function of dominant runoff production mechanism are conducted using data from two very different watersheds, the Hortonian Goodwin Creek Experimental Watershed and the non-Hortonian Muddy Brook watershed. Total infiltration, stream flow and evapotranspiration for the entire simulation period are used to compute comparison statistics. The influences of upper and lower boundary conditions on the solution accuracy are also explored. Results indicate that to simulate hydrological fluxes accurately at both watersheds small vertical cell sizes, of the order of 1 cm, are required near the soil surface, but not throughout the soil column. The appropriate choice of approximations for calculating the near soil-surface unsaturated hydraulic conductivity can yield modest increases in the required cell size. Results for both watersheds are quite similar, even though the soils and runoff production mechanisms differ greatly between the two catchments. Copyright (C) 2003 John Wiley Sons, Ltd.