Computer systems and methods are provided for seismic wavefield processing, including time domain reconstructed full wavefield inversion (TDRFWI) of seismic survey data. Suitable methods include forward propagating a source signal based on a model of a subsurface region, generating a residual based on the forward-propagated wavefield in comparison to field data, back propagating the residual to generate a reconstructed source perturbation, and adding the reconstructed source perturbation to the source signal to generate a reconstructed source wavefield. The reconstructed source perturbation can be forward propagated to generate a reconstructed wavefield perturbation, and added to the reconstructed wavefield perturbation to reconstruct the wavefield. The propagations can be performed in the time domain, and the model can be updated based on the reconstructed wavefield and reconstructed source perturbation to generate high fidelity subsurface images.