In regions characterized by complex subsurface structure, wave-equation depth migration is a powerful tool for accurately imaging the earth’s interior. The quality of the final image greatly depends on the quality of the model which includes anisotropy parameters (Gray et al., 2001). In particular, it is important to construct subsurface velocity models using techniques that are consistent with the methods used for imaging. Generally speaking, there are two possible strategies for velocity estimation from surface seismic data in the context of wavefield-based imaging (Sava et al., 2010). One possibility is to formulate an objective function in the data space, prior to migration, by matching the recorded data with simulated data. Techniques in this category are known by the name of waveform inversion. Another possibility is to formulate an objective function in the image space, after migration, by measuring and correcting image features that indicate model inaccuracies. Techniques in this category are known as wave-equation migration velocity analysis (MVA).