1. A method, stored in a memory and executing on a processor, for estimating a surface-consistent amplitude correction and a deconvolution operator, said method comprising: simultaneously utilizing a plurality of input seismic datasets, in a single pass, to calculate an amplitude value and an autocorrelation spectrum, within an associated time window, for each of said plurality of input seismic datasets; transforming said amplitudes and said autocorrelations to a log domain and summing for each surface-consistent component value to create an initial amplitude log scalar and an initial spectra log-frequency sum estimate; performing one or more iterations to simultaneously and independently estimate an improved amplitude log scalar and a spectra log-frequency value, with a shared solver, for each surface-consistent component until a shared predetermined stop condition for said estimating is met; and outputting one or more of an estimated surface-consistent amplitude correction and a deconvolution operator based on reconstructing an amplitude scalar and a autocorrelation spectrum, from said plurality of input seismic datasets, after said stop condition is met.