A method for estimating a time variant signal representing a seismic source obtains seismic data recorded by at least one receiver and generated by the seismic source, the recorded seismic data comprising direct arrivals and derives the time variant signal using an operator that relates the time variant signal to the acquired seismic data, the operator constrained such that the time variant signal is sparse in time.