In certain embodiments, a method includes generating a set of average velocity controls based on received seismic data, generating a depth to basement model based on received potential fields data, and generating an average velocity model using an interpolation model to interpolate the set of average velocity controls and the depth to basement model. The method may also include generating a variogram model based on the set of average velocity controls, and generating the average velocity model using the interpolation model to interpolate the average velocity controls, the variogram model, and the depth to basement model. The interpolation model may be Kriging with external drift. The external drift may be based on the depth to basement model. Additionally, the method includes generating a structural map of a subsurface formation using the average velocity model.