A basically time-domain method for performing full wavefield inversion of seismic data to infer a subsurface physical property model (61), where however at least one quantity required for the inversion, such as the Hessian of the cost function, is computed in the frequency domain (64). The frequency-domain quantity or quantities may be obtained at only a few discrete frequencies (62), preferably low frequencies, and may be computed on a coarse spatial grid, thus saving computing time with minimal loss in accuracy. For example, the simulations of predicted data and the broadband gradient of the objective function may be computed in the time domain (67), and the Hessian matrix, approximated by its diagonal, may be computed in the frequency domain. It may be preferable to use time-domain and the frequency-domain solvers that employ different numerical schemes, such as finite-difference method, one-way wave equation, finite-element method (63).