A method for use in seismic exploration of a region of the subsurface of the earth, wherein is provided a plurality of continuously recorded seismic traces collected in a seismic survey conducted proximate to the region of the subsurface of the earth, each of said traces being recorded by a receiver having a clock, the method being characterized by comprising the following steps: (a) identifying a plurality of proximate trace pairs, each proximate trace pair comprising a first trace recorded by a first receiver and a second trace recorded by a second receiver proximate to said first receiver, wherein the second receiver is proximate to the first receiver if it records signals from the same seismic event; (b) for each of said proximate trace pairs, determining a relative time shift between said first trace and said second trace; (c) inverting said relative time shifts determined in step (b) to obtain a clock drift estimate for each of the clocks of said first and second receivers by solving the equation: 1 − 1... 0 0 1 0 − 1... 0 ⋮ ⋮ 0 0... 1 − 1 Δt 1 Δt 2 ⋮ Δt N = Δ 1,2 Δ 1,3 ⋮ Δ N − 1, N for Δ t i, i = 1, N, where Δ t i is said clock drift estimate associated with an ith one of said clocks of said first and second receivers, where N is the number of said clocks of said first and second receivers, and, where Δ k,j is said relative time shift between a kth one of said first traces and a jth one of said second traces; and (d) using said clock drift estimates obtained in step (c) to time correct the first and second traces of each proximate trace pair, thereby forming a plurality of time-corrected seismic traces representative of the subsurface of the earth.. (a) identifying a plurality of proximate trace pairs, each proximate trace pair comprising a first trace recorded by a first receiver and a second trace recorded by a second receiver proximate to said first receiver, wherein the second receiver is proximate to the first receiver if it records signals from the same seismic event; (b) for each of said proximate trace pairs, determining a relative time shift between said first trace and said second trace; (c) inverting said relative time shifts determined in step (b) to obtain a clock drift estimate for each of the clocks of said first and second receivers by solving the equation: 1 − 1... 0 0 1 0 − 1... 0 ⋮ ⋮ 0 0... 1 − 1 Δt 1 Δt 2 ⋮ Δt N = Δ 1,2 Δ 1,3 ⋮ Δ N − 1, N for Δ t i, i = 1, N, where Δ t i is said clock drift estimate associated with an ith one of said clocks of said first and second receivers, where N is the number of said clocks of said first and second receivers, and, where Δ k,j is said relative time shift between a kth one of said first traces and a jth one of said second traces; and for Δ t i, i = 1, N, where Δ t i is said clock drift estimate associated with an ith one of said clocks of said first and second receivers, where N is the number of said clocks of said first and second receivers, and, where Δ k,j is said relative time shift between a kth one of said first traces and a jth one of said second traces; and (d) using said clock drift estimates obtained in step (c) to time correct the first and second traces of each proximate trace pair, thereby forming a plurality of time-corrected seismic traces representative of the subsurface of the earth..