1. A computer-implemented method for use in geophysical exploration, comprising: accessing a set of source time functions, a set of measured array data, and a subsurface attribute model, the subsurface attribute model defining a modeling grid; forward modeling the source time functions in the subsurface attribute model, including: generating a set of modeled receiver array data from modeled arrays gathered from the modeling grid; and differencing the measured receiver array data from the modeled array data to obtain a data residual for each gathered, modeled array; back propagating the data residual into the subsurface attribute model to produce a modeled residual wavefield; and updating the subsurface attribute model using the modeled residual wavefield; and iterating updates to the surface attribute model to convergence.