Full-waveform inversion (FWI) has become a popular method to retrieve high-resolution subsurface model parameters. It is a highly nonlinear optimization problem based on minimizing the misfit between the observed and predicted data. For intrinsically attenuating media, wave propagation experiences significant loss of energy. Thus, for better data fitting, it is sometimes crucial to consider attenuation in FWI. Viscoacoustic FWI aims at achieving a joint inversion of the velocity and attenuation models. However, multiparameter FWI imposes additional challenges including expanding the null space and facing parameter trade-off issues. Theoretically, an ideal way to mitigate the trade-off issue in multiparameter FWI is to apply the inverse Hessian operator to the parameter gradients. However, it is often not practical to calculate the full Hessian and its matrix inverse because this will be extremely expensive. To improve the computational efficiency and mitigate the trade-off issue, we have used an efficient wavefield inversion (EWI) method to invert for the velocity and the intrinsic attenuation. This approach is implemented in the frequency domain, and the velocity, in this case, is complex-valued in the viscoacoustic EWI. We evaluate a sequential update strategy for the velocity and the intrinsic attenuation, and we repeat the separate optimizations, which we refer to as outer iterations, until the convergence is achieved. Because viscoacoustic EWI is able to recover an accurate velocity model, the velocity update leakage to the model is largely reduced. We determine the effectiveness of this approach using synthetic data generated for the viscoacoustic Marmousi and Overthrust models. To further demonstrate the validity of our approach, we generate data in the time domain using a viscoelastic wave equation solver and obtain reasonable inversion results in the frequency domain using the viscoacoustic approximation.