Described is a finite-difference methodology for efficiently computing the response from a model subject to changes within sub-volumes with the sub-volume enclosed with an extrapolation surface within an injection boundary using Green's functions to update the injection boundary and transmitting an updated wavefield from the injection boundary back into the altered sub-volume to include higher order reflections from outside of the altered sub-space.