Development and Validation of a Three‐Dimensional Variably Saturated Flow Model for Future Water Resource Assessment at a Global Scale—Targeting Saturated Groundwater Flow in Plains

Abstract

A three-dimensional variably saturated flow model was developed for assessing future global water resources and parameterized for groundwater pumping. We applied this model to an actual watershed to verify its validity as an Earth System Model. For global applicability, the parameterization method for multi-layered groundwater pumping was developed and verified through comparison with observations and MODFLOW results. The parameterization proposed in this study is applicable even when multiple groundwater pumping wells are present within one horizontal computational grid and when the well spans multiple vertical grids. This method can be applied at the global scale without parameters such as the well radius, for which data may be difficult to obtain. The parameterization recreated seasonal and annual variations in the observed values. Furthermore, the results were comparable to those of MODFLOW. However, the calculation results were overestimated relative to the observed values. This overestimation was likely to be due to active groundwater pumping in the Central Valley before the start of the unsteady-state calculation. Therefore, the groundwater level at the beginning of the unsteady-state calculation was calculated using observed values, improving reproducibility. Furthermore, as observed groundwater levels are unlikely to be available at the global scale, steady-state calculations were conducted over 15 and 60 years considering groundwater pumping. However, the results were not as reproducible as those obtained using observed groundwater levels. These results suggest that the groundwater level set at the beginning of the calculation is important for global-scale groundwater flow calculation.