Abstract: High resolution three-dimensional (3D) spatial information of soil pollutants at contaminated site is extremely important for pollution risk assessment, remediation, and precise management and control of the site. This paper first introduced the main characteristics and affecting factors of pollutant variability and spatial distribution at contaminated site, and then systematically reviewed the commonly used 3D spatial interpolation method of site pollutant in existing studies and engineering applications. The shortcomings of the existing 3D interpolation method, and the research gaps between 3D spatial interpolation method and auxiliary information incorporation, as well as the uncertainty assessment were also identified. Finally, the research suggestions were proposed from the aspects of soil borehole sampling scheme design, auxiliary information acquisition, spatial interpolation model development, and technological specification of 3D interpolation application.