We present a scheme for multiscale phase inversion (MPI) of seismic data that is less sensitive than full-waveform inversion (FWI) to the unmodeled physics of wave propagation and to a poor starting model. To avoid cycle skipping, the multiscale strategy temporally integrates the traces several times, i.e., high-order integration, to produce low-boost seismograms that are used as input data for the initial iterations of MPI. As the iterations proceed, lower frequencies in the data are boosted by using integrated traces of lower order as the input data. The input data are also filtered into different narrow frequency bands for the MPI implementation. Numerical results with synthetic acoustic data indicate that, for the Marmousi model, MPI is more robust than conventional multiscale FWI when the initial model is moderately far from the true model. Results from synthetic viscoacoustic and elastic data indicate that MPI is less sensitive than FWI to some of the unmodeled physics. Inversion of marine data indicates that MPI is more robust and produces modestly more accurate results than FWI for this data set.