Surface consistent amplitude and deconvolution simultaneous joined inversion

A method and system for estimating a surface-consistent amplitude correction and deconvolution operator in a simultaneous joined inversion comprises a one-pass estimation using input trace data for generating the gain correction and deconvolution operator based on a least squares iteration method. A series of iterations (e.g. Gauss Seidel) are performed simultaneously and independently estimating amplitude scalars and autocorrelation spectra with a common convergence criterion. The gain and deconvolution operator can be used to correct the input trace data for pre-stack or stack imaging.