Continuous assimilation of simulated Geosat altimetric sea level into an eddy‐resolving numerical ocean model: 2. Referenced sea level differences

Continuous assimilation of referenced altimetric sea level differences from the simulated Geosat Exact Repeat Mission (ERM) into a three-layer, wind-driven, quasi-geostrophic, eddy-resolving numerical ocean box model at mid-latitude was conducted utilizing an improved version of the optimal interpolation procedure found in Part I (White et al., this issue), forming the basis of an operational dynamical interpolation procedure to be used in subsequent gridding of real Geosat ERM altimetric sea level data. The major improvement over the dynamical interpolation method used in Part I was the inclusion of temporal decorrelation in the error covariance matrix and the introduction of noise into the system. As referenced altimetric sea level differences (i.e., differences referenced by the mean of the model) were observed from a control model integration in real model time by the simulated Geosat altimeter, attention focused upon the ability of the dynamical interpolation procedure to reconstruct a mesoscale eddy field resembling that of the control model integration. Results of this dynamical interpolation procedure were compared with those of a statistical (i.e., optimum) interpolation procedure in its ability to reconstruct the mesoscale eddy field in the western boundary extension current of the control model integration. A series of model integrations were made, with initial conditions independent from those of the control model integration (the latter taken to represent the observed mesoscale eddy field), where the dynamical interpolation of simulated referenced Geosat ERM altimetric sea level differences was found capable of reproducing the control model integration (i.e., explaining 80–90% of the control model variance) after only one repeat cycle (i.e., 17 days) of the ERM. This success was attributed to the use of time decorrelation in the error covariance (together with noise) to nudge the first-guess model upper level stream function field toward the corresponding observation field over a period of a few days, reducing the shock of inserting the observations into the model. The degree of similarity between the upper layer stream function of the updated model and the control model integration after one repeat cycle was, however, only marginally significantly greater than that achieved by statistical interpolation (i.e., explaining 65–90% of the control model variance), both interpolation procedures (i.e., dynamical and statistical) using covariance estimates that were determined a priori in both space and time. Two important advantages of the dynamical interpolation over the statistical interpolation became evident after three cycles (i.e., 51 days) of continuous assimilation of referenced simulated Geosat ERM altimetric sea level differences. First, mesoscale eddy activity in the middle and lower layers of the updated model integration began to resemble that of the model control integration (i.e., explaining 65% and 15% of the control model variance, respectively), without resorting to vertical statistical interpolation. Second, the spectral coherence between the upper layer stream function field of the updated model integration and that of the control model integration was symmetric with respect to direction of propagation, while the spectral coherence between the upper layer stream function field by statistical interpolation and that of the control model integration was not, with the latter robust enough only to detect the dominant westward propagating portion of the spectrum of mesoscale eddy activity and not the subdominant information. An important finding was the fact that both interpolation procedures were able to resolve wavelengths longer than the Nyquist (i.e., 280 km) dictated by the track separation distance, as expected, but both interpolation procedures were able to resolve periods (i.e., 12 days) much shorter than the Nyquist (i.e., 34 days) dictated by the exact repeat period.

[1]  Brian D. Beckley,et al.  Global mesoscale variability from collinear tracks of SEASAT altimeter data , 1983 .

[2]  W. White,et al.  Design of an Oceanographic Network in the Midlatitude North Pacific , 1979 .

[3]  Robert N. Miller Toward the Application of the Kalman Filter to Regional Open Ocean Modeling , 1986 .

[4]  W. R. Holland,et al.  Observed and modeled mesoscale variability near the Gulf Stream and Kuroshio Extension , 1986 .

[5]  H. Dantzler Geographic variations in intensity of the North Atlantic and North Pacific oceanic eddy fields , 1976 .

[6]  W. R. Holland,et al.  Recent long‐term moored instrument observations in the western North Pacific , 1982 .

[7]  J. Kindle Sampling strategies and model assimilation of altimetric data for ocean monitoring and prediction , 1986 .

[8]  William J. Schmitz,et al.  Mid‐latitude mesoscale variability , 1983 .

[9]  A. Robinson,et al.  A real-time dynamical forecast of ocean synoptic/mesoscale eddies , 1984, Nature.

[10]  Chester J. Koblinsky,et al.  Estimates of the geostrophic stream function in the western North Pacific from XBT surveys. [expendable bathythermograph , 1984 .

[11]  J. Thompson,et al.  Altimeter data and geoid error in mesoscale ocean prediction: Some results from a primitive equation model , 1986 .

[12]  A. Robinson,et al.  Estimation and prediction of oceanic Eddy fields , 1985 .

[13]  W. White,et al.  Time and Length Scales of Baroclinic Eddies in the Central North Pacific Ocean , 1974 .

[14]  W. R. Holland Simulation of Mesoscale Ocean Variability in Mid-Latitude Gyres , 1985 .

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

[16]  W. Emery On the Geographical Variability of the Upper Level Mean and Eddy Fields in the North Atlantic and North Pacific , 1983 .

[17]  W. R. Holland,et al.  A preliminary comparison of selected numerical eddy-resolving general circulation experiments with observations. , 1982 .

[18]  L. Gandin Objective Analysis of Meteorological Fields , 1963 .

[19]  W. R. Holland,et al.  On the Generation of Mesoscale Eddies and their Contribution to the OceanicGeneral Circulation. II. A Parameter Study , 1975 .

[20]  William R. Holland,et al.  On the Generation of Mesoscale Eddies and their Contribution to the OceanicGeneral Circulation. I. A Preliminary Numerical Experiment , 1975 .

[21]  A. Robinson,et al.  Dynamical Interpolation and Forecast of the Evolution of Mesoscale Features off Northern California , 1987 .

[22]  W. White,et al.  Annual Rossby Wave Characteristics in the California Current Region from the GEOSAT Exact Repeat Mission , 1990 .

[23]  W. R. Holland,et al.  The Stability of Ocean Currents in Eddy-Resolving General Circulation Models , 1978 .

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

[25]  William R. Holland,et al.  The Role of Mesoscale Eddies in the General Circulation of the Ocean—Numerical Experiments Using a Wind-Driven Quasi-Geostrophic Model , 1978 .

[26]  W. White,et al.  Stationary and Traveling Mesoscale Perturbations in the Kuroshio Extension Current , 1981 .

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

[28]  D. Webb,et al.  Assimilation of Altimeter Data into Ocean Models , 1986 .

[29]  W. Budgell,et al.  Ocean Data Assimilation and the Moan Filter: Spatial Regularity , 1987 .

[30]  A. Lorenc A Global Three-Dimensional Multivariate Statistical Interpolation Scheme , 1981 .

[31]  W. White,et al.  Horizontal Eddy Fluxes of Momentum and Kinetic Energy in the Near-Surface of the Kuroshio Extension , 1982 .