We introduce an effective iterative Least-Squares Wave-Equation Migration (LS-WEM) solution for broadband imaging. Least-Squares Migration (LSM) solutions are designed to produce images of the subsurface corrected for wavefield distortions caused by acquisition and propagation effects. They implicitly solve for the earth reflectivity by means of data residual reduction in an iterative fashion, which usually demands intensive computation. The LS-WEM is implemented using a visco-acoustic anisotropic one-way wave-equation wavefield propagator that is able to fully utilize both the broader seismic bandwidth and the high-resolution velocity information from Full Waveform Inversion (FWI). Our implementation combines the one-way extrapolator with fast linear inversion solvers into an efficient migration inversion system. Application to the 2D Sigsbee2b synthetic model improves the sub-salt illumination by balancing the image amplitudes and reducing the effects of the shadow zones, enhances temporal resolution by broadening the frequency spectrum, balances the wavenumber content and improves images of faults and dipping salt flanks. In addition, LS-WEM converges rapidly to the true solution, reducing the data residuals by 90% in only four iterations. Application to real 3D datasets from the Gulf of Mexico and the North Sea demonstrates high-resolution imaging with reduced acquisition footprint effects, improved spatial frequency content, and better structural imaging at all depths.