Adjoint assimilation of altimetric, surface drifter, and hydrographic data in a quasi‐geostrophic model of the Azores Current

The flow characteristics in the region of the Azores Current are investigated by assimilating TOPEX/POSEIDON and ERS 1 altimeter data into the multilevel Harvard quasigeostrophic (QG) model with open boundaries (Miller et al., 1983) using an adjoint variational scheme (Moore, 1991). The study site lies in the path of the Azores Current, where a branch retroflects to the south in the vicinity of the Madeira Rise. The region was the site of an intensive field program in 1993, SEMAPHORE. We had two main aims in this adjoint assimilation project. The first was to see whether the adjoint method could be applied locally to optimize an initial guess field, derived from the continous assimilation of altimetry data using optimal interpolation (OI). The second aim was to assimilate a variety of different data sets and evaluate their importance in constraining our QG model. The adjoint assimilation of surface data was effective in optimizing the initial conditions from OI. After 20 iterations the cost function was generally reduced by 50–80%, depending on the chosen data constraints. The primary adjustment process was via the barotropic mode. Altimetry proved to be a good constraint on the variable flow field, in particular, for constraining the barotropic field. The excellent data quality of the TOPEX/POSEIDON (T/P) altimeter data provided smooth and reliable forcing; but for our mesoscale study in a region of long decorrelation times O(30 days), the spatial coverage from the combined T/P and ERS 1 data sets was more important for constraining the solution and providing stable flow at all levels. Surface drifters provided an excellent constraint on both the barotropic and baroclinic model fields. More importantly, the drifters provided a reliable measure of the mean field. Hydrographic data were also applied as a constraint; in general, hydrography provided a weak but effective constraint on the vertical Rossby modes in the model. Finally, forecasts run over a 2-month period indicate that the initial conditions optimized by the 20-day adjoint assimilation provide more stable, longer-term forecasts.

[1]  Optimal Interpolation in a model of the Azores Current in 1986–88 , 1994 .

[2]  J. Neumann,et al.  Numerical Integration of the Barotropic Vorticity Equation , 1950 .

[3]  Eric Greiner,et al.  Assimilation of Geosat Altimetric Data in a Nonlinear Reduced-Gravity Model of the Indian Ocean. Part 1: Adjoint Approach and Model-Data Consistency , 1994 .

[4]  Carl Wunsch,et al.  Assimilation of Sea Surface Topography into an Ocean Circulation Model Using a Steady-State Smoother , 1993 .

[5]  J. Schröter.  A sensitivity analysis of variational data assimilation , 1994 .

[6]  P. L. Traon,et al.  The eddy field associated with the Azores Front east of the Mid‐Atlantic Ridge as observed by the Geosat altimeter , 1994 .

[7]  Jens Schröter,et al.  Variational Assimilation of Geosat Data into an Eddy-resolving Model of the Gulf Stream Extension Area , 1993 .

[8]  M. Yaremchuk,et al.  Conductivity-temperature-depth data assimilation into a three-dimensional quasigeostrophic open ocean model , 1994 .

[9]  Dale B. Haidvogel,et al.  A baroclinic quasigeostrophic open ocean model , 1983 .

[10]  Allan R. Robinson,et al.  Assimilation of altimeter Eddy fields in a limited-area quasi-geostrophic model , 1987 .

[11]  Hans-Harald Hinrichsen,et al.  Can meddies be detected by satellite altimetry , 1991 .

[12]  W. J. Gould,et al.  Physical oceanography of the Azores front , 1985 .

[13]  P. Demey.  Synoptic estimates of an eddy field in the North Atlantic Current , 1992 .

[14]  B. Klein,et al.  On the origin of the Azores Current , 1989 .

[15]  Eric Dombrowsky,et al.  Continuous assimilation in an open domain of the northeast Atlantic: 1. Methodology and application to AthenA‐88 , 1992 .

[16]  F. Hernandez.  Comparaison et combinaison de données altimétriques et lagrangiennes. Application à la campagne SEMAPHORE-93 , 1995 .

[17]  W. Thacker,et al.  Fitting dynamics to data , 1988 .

[18]  C. Provost,et al.  General Circulation of the Mid-Latitude Ocean: Coupled Effects of Variable Wind Forcings and Bottom Topography Roughness on the Mean and Eddy Circulation. , 1989 .

[19]  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 .

[20]  Dale B. Haidvogel,et al.  The accuracy, efficiency, and stability of three numerical models with application to open ocean problems , 1980 .

[21]  H. Hurlburt,et al.  Dynamic transfer of simulated altimeter data into subsurface information by a numerical ocean model , 1986 .

[22]  W. Zenk,et al.  Large Lenses of Highly Saline Mediterranean Water , 1984 .

[23]  John Derber,et al.  A Global Oceanic Data Assimilation System , 1989 .

[24]  Andrew M. Moore,et al.  Data Assimilation in a Quasi-geostrophic Open-Ocean Model of the Gulf Stream Region Using the Adjoint Method , 1991 .

[25]  Wayne G. Leslie,et al.  Data assimilation and dynamical interpolation in GULFCAST experiments , 1989 .

[26]  U. Seiler.  Estimation of open boundary conditions with the adjoint method , 1993 .

[27]  E. Polak,et al.  Note sur la convergence de méthodes de directions conjuguées , 1969 .

[28]  C. Boissier,et al.  Spatial scales of mesoscale variability in the North Atlantic as deduced from Geosat data , 1990 .

[29]  F. L. Dimet,et al.  Variational algorithms for analysis and assimilation of meteorological observations: theoretical aspects , 1986 .

[30]  Margaret King Robinson,et al.  Atlas of North Atlantic-Indian Ocean Monthly Mean Temperatures and Mean Salinities of the Surface Layer. , 1979 .

[31]  K. Haines,et al.  A Direct Method for Assimilating Sea Surface Height Data into Ocean Models with Adjustments to the Deep Circulation , 1991 .