Verification and application results are presented for a new saturated-unsaturated 3D groundwater flow model (SU3D) developing in Dogan and Motz, which can be used to calculate the pressure distribution over the entire groundwater flow domain in response to rainfall and evapotranspiration. SU3D solves the nonlinear 3D, modified mixed form of the Richards equation continuously throughout the groundwater flow domain, including both the unsaturated and saturated zones. The block-centered finite-difference method, a modified Picard iteration scheme, and the preconditioned conjugate gradient method are used to solve the governing partial differential equations in the new model. SU3D can simulate evaporation from land surface and transpiration from the root zone. Potential evapotranspiration is partitioned into potential evaporation and potential transpiration as a function of the leaf area index; actual evaporation and actual transpiration are then calculated individually. As shown in this paper, SU3D has been verified by reproducing the results of five analytical and numerical studies and a 3D unconfined-aquifer pumping test from the published literature.