When present in the subsurface salt bodies impact the complexity of wave-equation-based seismic imaging techniques, such as least-squares reverse-time migration, and full-waveform inversion (FWI). Typically, the Born approximation used in every iteration of least-squaresbased inversions is incapable of handling the sharp, high-contrast boundaries of salt bodies.We develop a variance-based method for reconstruction of velocity models to resolve the imaging and inversion issues caused by salt bodies. Our main idea lies in retrieving useful information from independent updates corresponding to FWI at di↵erent frequencies. After several FWI iterations we compare the model updates by considering the variance distribution between them to identify locations most prone to cycle skipping. We interpolate velocities from the surrounding environment into these high-variance areas. This approach allows the model to gradually improve from identifying easily resolvable areas and extrapolating the model updates from those to the areas that are di cult to resolve at early FWI iterations. In numerical tests, our method demonstrates the ability to obtain convergent FWI results at higher frequencies.Downloaded 08/06/18 to 109.171.137.210. Redistribution subject to SEG license or copyright; see Terms of Use at http://library.seg.org/ de Hoop (2013), who used a total variation norm to invert a salt-a↵ected subsurface using FWI. Regularized inversion methods could be computationally expensive when multiple runs on a fine mesh are needed and these methods often require prior assumptions about the subsurface structure (Kadu et al., 2016).The second approach involves image or gradient manipulations with single iteration updates of FWI. Processed gradients can lead to shorter paths toward the global minimum without being trapped in the local minima (Alkhalifah, 2015b. Image-guided inversion (Ma et al., 2012), gradient optimization (Wu and Alkhalifah, 2016) or gradient conditioning through scattering angle-based filters (Alkhalifah, 2015a;Kazei et al., 2016) can serve the same purpose.The objective of FWI is to minimize both amplitude and phase di↵erences between observed and modeled seismic data. For successful inversion the Born approximation requires the initial velocity model to deliver mismatches in travel times less than half the period (Beydoun and Tarantola, 1988). Larger errors cause cycle-skipping artifacts on the wavepaths specific to each source-receiver pair. These artifacts in the image domain can be identified as repeated contrast velocity anomalies that, in turn, lead the inversion to convergence toward a local minimum on an objective function (Virieux and Operto, 2009).Inversion of mono-frequency data allows features to be resolved at a specific scale. The inversion results from di↵erent mono-frequencies therefore do not match exactly.When mono-frequency data are modeled, the phase of a harmonic plane wave at a given point on a wavepath depends on the frequency and travel time. Hence, at di↵erent frequencies, the phases of ...