Traditional analytical solutions for the evaluation of pumping tests under non-steady state conditions are applicable when the aquifer can be simplified as constant horizontal strata of infinite extension and homogeneous soil. Nevertheless, when such conditions are not representative of the study area, analytical solutions present considerable limitations. This paper discusses the application of 3D numerical modeling as a powerful tool to evaluate complex pumping well systems when it is mandatory to represent topography and variable thickness stratigraphic conditions of the study area. Furthermore, application of the analytical solutions is highlighted as a previous analysis to characterize the hydraulic properties of the aquifer (permeability, transmissivity and storage coefficient) and be took into account in the numerical analyses. The practical application presented in this paper utilized data from the northeastern part of the Valley of Mexico, particularly from the area of the former Texcoco Lake, where the construction of a new airport was planned. At the end some conclusions are given about the design and analysis of pumping well systems.