Summary This work introduces an iterative geostatistical seismic time-lapse inversion based on stochastic sequential simulation and co-simulation. The proposed approach inverts simultaneously all seismic vintages and partial angle stacks directly for differences of Ip and Is. The application to a challenging non-stationary environment shows promissing results.