Abstract: The method of the present invention provides for processing seismic data over a subsurface area of interest. Downgoing seismic data, which may be VSP data, and surface seismic data are acquired over a subsurface area of interest. An inverse operator is determined from changes in signal characteristics, which may be first-arrival signals, between consecutive depth levels of the downgoing VSP data. The inverse operator is assigned to at least one surface seismic data level. Data levels may be in time or depth. Inverse operators may be interpolated for time or depth level operators for data samples between already determined operators. Inverse operators are applied to seismic data to restore attenuated signal components.