One of the main challenges in seismic imaging, especially of land data, is building the near surface velocity model, as it is normally characterized by very low velocity values with different types of local anomalies. Resolving the near surface velocity model using only refraction data is insufficient, as it does not provide the required lateral resolution. Using only reflection data is equally insufficient, as it does not provide the required vertical resolution.