Abstract-In this paper, we propose a data assimilation method for consistently estimating the velocity fields from a whole image sequence depicting the evolution of sea surface temperature transported by oceanic surface flow. The estimator is conducted through an ensemble-based variational data assimilation, which is designed by combining the advantages of two approaches: the ensemble Kalman filter and the variational data assimilation. This idea allows us to obtain the optimal initial condition as well as the full system trajectory. In order to extract the velocity fields from fluid images, a surface quasi-geostrophic model representing the generic evolution of the temperature field of the flow, and the optical flow constraint equation derived from the image intensity constancy assumption, are involved in the assimilation context. Numerical experimental evaluation is presented on a synthetic fluid image sequence. The results indicate good performance and efficiency of the proposed estimator.