The complexity of dynamical laws governing 3D atmospheric flows associated with incomplete and noisy observations make the recovery of atmospheric dynamics from satellite images sequences very difficult. In this paper, we face the challenging problem of estimating physical sound and time-consistent horizontal motion fields at various atmospheric depths for a whole image sequence. Based on a vertical decomposition of the atmosphere, we propose a dynamically consistent atmospheric motion

more »

... relying on a multi-layer dynamical model. This estimator is based on a weak constraint variational data assimilation scheme and is applied on noisy and incomplete pressure difference observations derived from satellite images. The dynamical model consists in a simplified vorticity-divergence form of a multi-layer shallow-water model. Average horizontal motion fields are estimated for each layer. The performance of the proposed technique is assessed on real world meteorological satellite image sequences.