Estimating surface velocities from satellite data and numerical models: Implementation and testing of a new simple method

Abstract A simple method of fusing tracer observations and model outputs for computing surface velocities in the ocean is implemented and tested in the framework of the twin experiment approach. Synthetic data from realistic velocity outputs produced by the operational Mediterranean Forecasting System (MFS) are used. The method ( Piterbarg, 2009 ) allows to estimate a velocity field using two consecutive tracer snapshots. The focus is on testing realistic time intervals between snapshots and partial tracer observations. The considered configuration consists of a tracer patch released and advected by the current, and is motivated by the practical problem of estimating velocities and concentrations using satellite data in case of pollutant releases such as oil spills. An extensive set of experiments has been carried out, and the method performance has been quantified in terms of improvements in accuracy with respect to the model. The improvement ranges from values of approximately 80–90% for concentration and 50–60% for velocity in the case of almost perfect data, to values of 30–40% for realistic time intervals of the order of days and reduced tracer information, and values of 15–20% when only the boundary of the patch is observed. The results are found to be robust to flow variability and patch parameters.

[1]  Gregoire Mercier,et al.  DETECTION AND MAPPING OF OIL SLICKS IN THE SEA BY COMBINED USE OF HYPERSPECTRAL IMAGERY AND LASER-INDUCED FLUORESCENCE , 2006 .

[2]  A. H. Al-Rabeh,et al.  Modeling the fate and transport of Al-Ahmadi oil spill , 1992 .

[3]  M. McIntyre,et al.  On the Conservation and Impermeability Theorems for Potential Vorticity , 1990 .

[4]  Joseph H. LaCasce,et al.  Statistics from Lagrangian observations , 2008 .

[5]  G. Budillon,et al.  Tracers and Potential Vorticities in Ocean Dynamics , 2002 .

[6]  Gary B. Brassington,et al.  GODAE systems in operation , 2009 .

[7]  Patrick Pérez,et al.  Dense Estimation of Fluid Flows , 2002, IEEE Trans. Pattern Anal. Mach. Intell..

[8]  Nicolas Papadakis,et al.  A Variational Approach for Object Contour Tracking , 2005, VLSM.

[9]  Joachim Weickert,et al.  Scale Space and PDE Methods in Computer Vision, 5th International Conference, Scale-Space 2005, Hofgeismar, Germany, April 7-9, 2005, Proceedings , 2005, Scale-Space.

[10]  M. Astraldi,et al.  Temporal variability of currents in the eastern Ligurian Sea , 1990 .

[11]  Patrick Pérez,et al.  A multigrid approach for hierarchical motion estimation , 1998, Sixth International Conference on Computer Vision (IEEE Cat. No.98CH36271).

[12]  C. Wunsch,et al.  Determining the General Circulation of the Oceans: A Preliminary Discussion , 1977, Science.

[13]  F. L. Dimet,et al.  Coupling models and data: Which possibilities for remotely-sensed images? , 2004 .

[14]  S. Bauer.  Eddy-mean flow decomposition and eddy-diffusivity estimates in the tropical Pacific Ocean , 1998 .

[15]  Isabelle Herlin,et al.  Optical Flow and Phase Portrait Methods for Environmental Satellite Image Sequences , 1996, ECCV.

[16]  Isabelle Herlin,et al.  A generalized optical flow constraint and its physical interpretation , 2000, Proceedings IEEE Conference on Computer Vision and Pattern Recognition. CVPR 2000 (Cat. No.PR00662).

[17]  Kathryn A. Kelly,et al.  An Inverse Model for Near-Surface Velocity from Infrared Images , 1989 .

[18]  A. Mariano,et al.  Assimilation of drifter observations for the reconstruction of the Eulerian circulation field , 2003 .

[19]  V. Taillandier,et al.  Assimilation of Argo float positions in the north western Mediterranean Sea and impact on ocean circulation simulations , 2006 .

[20]  Henry Stommel,et al.  The beta spiral and the determination of the absolute velocity field from hydrographic station data , 1977 .

[21]  Leonid I. Piterbarg,et al.  Inversion of Upper Ocean Temperature Time Series for Entrainment, Advection, and Diffusivity , 2000 .

[22]  Philippe Courtier,et al.  Sea surface velocities from sea surface temperature image sequences: 1. Method and validation using primitive equation model output , 2000 .

[23]  R. Reynolds,et al.  Testing a Dynamical Model for Mid-Latitude Sea Surface Temperature Anomalies , 1983 .

[24]  A. H. Al-Rabeh,et al.  Gulfspill Version 2.0: a software package for oil spills in the Arabian Gulf , 2000, Environ. Model. Softw..

[25]  Marina Tonani,et al.  under a Creative Commons License. Ocean Science A high-resolution free-surface model of the Mediterranean Sea , 2007 .

[26]  Anne Cuzol,et al.  Vortex and Source Particles for Fluid Motion Estimation , 2005, Scale-Space.

[27]  Implementation of position assimilation for ARGO floats in a realistic Mediterranean Sea OPA model and twin experiment testing , 2006 .

[28]  D. Durran.  Numerical methods for wave equations in geophysical fluid dynamics , 1999 .

[29]  Dirk Olbers,et al.  The Inference of North Atlantic Circulation Patterns From Climatological Hydrographic Data (Paper 5R0560) , 1985 .

[30]  Thomas M. Smith,et al.  An Improved In Situ and Satellite SST Analysis for Climate , 2002 .

[31]  William J. Emery,et al.  Computing Coastal Ocean Surface Currents From Infrared and Ocean Color Satellite Imagery , 2007, IEEE Transactions on Geoscience and Remote Sensing.

[32]  D. Olbers,et al.  Potential Vorticity Constraints on the Dynamics and Hydrography of the Southern Ocean , 1993 .

[33]  Joaquim Ballabrera-Poy,et al.  Tracking oceanic currents by singularity analysis of Microwave Sea Surface Temperature images , 2008 .

[34]  Obtaining Velocities from Tracer Distributions , 1984 .

[35]  V. Taillandier,et al.  A variational approach for the reconstruction of regional scale Eulerian velocity fields from Lagrangian data , 2006 .

[36]  Enrico Zambianchi,et al.  Surface circulation in the central Mediterranean Sea as deduced from Lagrangian drifters in the 1990s , 2007 .

[37]  Marina Tonani,et al.  Mediterranean Forecasting System: An improved assimilation scheme for sea‐level anomaly and its validation , 2005 .

[38]  Didier Dubois,et al.  Fuzzy information engineering: a guided tour of applications , 1997 .

[39]  Isabelle Herlin,et al.  Data Assimilation of Satellite Images Within an Oceanographic Circulation Model , 2006, 2006 IEEE International Conference on Acoustics Speech and Signal Processing Proceedings.

[40]  K. Speer,et al.  Lagrangian Eddy Scales in the Northern Atlantic Ocean , 2002 .

[41]  K. Ide,et al.  Lagrangian data assimilation for point vortex systems , 2002 .

[42]  L. Kantha,et al.  Numerical models of oceans and oceanic processes , 2000 .

[43]  G. Georgiou,et al.  Operational ocean forecasting in the Eastern Mediterranean: implementation and evaluation , 2006 .

[44]  L. Piterbarg,et al.  A New Method for Obtaining Velocity and Diffusivity from Time-Dependent Distributions of a Tracer via the Maximum Likelihood Estimator for the Advection-Diffusion Equation , 1997 .

[45]  G. Needler.  The absolute velocity as a function of conserved measurable quantities , 1985 .

[46]  A. Solberg,et al.  Oil spill detection by satellite remote sensing , 2005 .

[47]  Nikos Paragios,et al.  Variational, Geometric, and Level Set Methods in Computer Vision, Third International Workshop, VLSM 2005, Beijing, China, October 16, 2005, Proceedings , 2005, VLSM.

[48]  Roberto Cipolla,et al.  Computer Vision — ECCV '96 , 1996, Lecture Notes in Computer Science.

[49]  Fabrice Hernandez,et al.  Mapping mesoscale variability of the Azores Current using TOPEX/POSEIDON and ERS 1 altimetry, together with hydrographic and Lagrangian measurements , 1995 .

[50]  Isabelle Herlin,et al.  Estimation of a Motion Field on Satellite Images from a Simplified Ocean Circulation Model , 2006, 2006 International Conference on Image Processing.

[51]  Rosalia Santoleri,et al.  Lagrangian and Eulerian observations of the surface circulation in the Tyrrhenian Sea , 2010 .

[52]  Isabelle Herlin,et al.  Retrieving ocean surface current by 4-D variational assimilation of sea surface temperature images , 2008 .

[53]  Annalisa Griffa,et al.  Transport Properties in the Adriatic Sea as Deduced from Drifter Data , 2000 .

[54]  W. Emery,et al.  An objective method for computing advective surface velocities from sequential infrared satellite images , 1986 .

[55]  William J. Emery,et al.  Satellite-Image-derived Gulf Stream Currents Compared with Numerical Model Results , 1992 .

[56]  K. Ide,et al.  A Method for Assimilating Lagrangian Data into a Shallow-Water-Equation Ocean Model , 2006 .

[57]  I. Taupier-Letage,et al.  Surface circulation in the Eastern Mediterranean using drifters (2005–2007) , 2009 .

[58]  R. Courant,et al.  Über die partiellen Differenzengleichungen der mathematischen Physik , 1928 .

[59]  L. Piterbarg.  A simple method for computing velocities from tracer observations and a model output , 2009 .

[60]  Leonid I. Piterbarg,et al.  Inversion for heat anomaly transport from sea surface temperature time series in the northwest Pacific , 1995 .

[61]  R. Signell,et al.  Variational analysis of drifter positions and model outputs for the reconstruction of surface currents in the central Adriatic during fall 2002 , 2008 .

[62]  A. Bennett.  Inverse Methods in Physical Oceanography , 1992 .

[63]  C. Wunsch.  The North Atlantic general circulation west of 50°W determined by inverse methods , 1978 .