When water resource or quality problems associated with phreatic water in geological settings are studied, information about water table elevations is often crucial. In those cases where the water table is a subdued replica of the ground surface, universal cokriging can be used to estimate water table elevations at un-sampled locations on the basis of water table and ground surface elevation measurements obtained at well locations. In this paper, universal cokriging, with ground surface elevation considered as a co-variate, was used to estimate water table elevations. Universal cokriging equations were derived, an iterative method for obtaining experimental variograms was established, and a case study of an initial groundwater flow simulation in the Xiuwu County, Henan Province, China, was presented. In the case study, the initial groundwater flow regime was represented both by universal cokriging with the ground surface elevation serving as a covariate and by universal kriging without the inclusion of ground surface elevation as a covariate. A comparison of the results from these two approaches shows that groundwater levels of phreatic water at locations without measurements in regions with unconsolidated porous media can be estimated more accurately by universal cokriging. When water resource or quality problems associated with phreatic water in geological settings are studied, information about water table elevations is often crucial. In those cases where the water table is a subdued replica of the ground surface, universal cokriging can be used to estimate water table elevations at un-sampled locations on the basis of water table and ground surface elevation measurements obtained at well locations. In this paper, universal cokriging, with ground surface elevation considered as a co-variate, was used to estimate water table elevations. Universal cokriging equations were derived, an iterative method for obtaining experimental variograms was established, and a case study of an initial groundwater flow simulation in the Xiuwu County, Henan Province, China, was presented. In the case study, the initial groundwater flow regime was exhibited both by universal cokriging with the ground surface elevation serving as a covariate and by universal kriging without the inclusion of ground surface elevation as a covariate. A comparison of the results from these two approaches shows that groundwater levels of phreatic water at locations without measurements in regions with unconsolidated porous media can be estimated more accurately by universal cokriging.