Adaptive error estimation in linearized ocean general circulation models
Citable URI
http://hdl.handle.net/1912/4743Location
5°60°N132°252°E
North Pacific
DOI
10.1575/1912/4743Abstract
Data assimilation methods, such as the Kalman filter, are routinely used in oceanography.
The statistics of the model and measurement errors need to be specified a priori.
In this study we address the problem of estimating model and measurement error statistics
from observations. We start by testing the Myers and Tapley (1976, MT) method of
adaptive error estimation with lowdimensional models. We then apply the MT method
in the North Pacific (5°60° N, 132°252° E) to TOPEX/POSEIDON sea level anomaly
data, acoustic tomography data from the ATOC project, and the MIT General Circulation
Model (GCM). A reduced state linear model that describes large scale internal
(baroclinic) error dynamics is used. The MT method, closely related to the maximum likelihood
methods of Belanger (1974) and Dee (1995), is shown to be sensitive to the
initial guess for the error statistics and the type of observations. It does not provide
information about the uncertainty of the estimates nor does it provide information about
which structures of the error statistics can be estimated and which cannot.
A new offline approach is developed, the covariance matching approach (CMA),
where covariance matrices of modeldata residuals are "matched" to their theoretical
expectations using familiar least squares methods. This method uses observations directly
instead of the innovations sequence and is shown to be related to the MT method and the
method of Fu et al. (1993). The CMA is both a powerful diagnostic tool for addressing
theoretical questions and an efficient estimator for real data assimilation studies. It can
be extended to estimate other statistics of the errors, trends, annual cycles, etc.
Twin experiments using the same linearized MIT GCM suggest that altimetric data
are illsuited to the estimation of internal GCM errors, but that such estimates can in
theory be obtained using acoustic data. After removal of trends and annual cycles, the low
frequency /wavenumber (periods> 2 months, wavelengths> 16°) TOPEX/POSEIDON
sea level anomaly is of the order 6 cm2. The GCM explains about 40% of that variance.
By covariance matching, it is estimated that 60% of the GCMTOPEX/POSEIDON
residual variance is consistent with the reduced state linear model.
The CMA is then applied to TOPEX/POSEIDON sea level anomaly data and a
linearization of a global GFDL GCM. The linearization, done in Fukumori et al.(1999),
uses two vertical mode, the barotropic and the first baroclinic modes. We show that the CMA method can be used with a global model and a global data set, and that the
estimates of the error statistics are robust. We show that the fraction of the GCMTOPEX/
POSEIDON residual variance explained by the model error is larger than that
derived in Fukumori et al.(1999) with the method of Fu et al.(1993). Most of the model
error is explained by the barotropic mode. However, we find that impact of the change
in the error statistics on the data assimilation estimates is very small. This is explained
by the large representation error, i.e. the dominance of the mesoscale eddies in the TIP
signal, which are not part of the 20 by 10 GCM. Therefore, the impact of the observations
on the assimilation is very small even after the adjustment of the error statistics.
This work demonstrates that simultaneous estimation of the model and measurement
error statistics for data assimilation with global ocean data sets and linearized GCMs is
possible. However, the error covariance estimation problem is in general highly underdetermined,
much more so than the state estimation problem. In other words there exist
a very large number of statistical models that can be made consistent with the available
data. Therefore, methods for obtaining quantitative error estimates, powerful though
they may be, cannot replace physical insight. Used in the right context, as a tool for
guiding the choice of a small number of model error parameters, covariance matching can
be a useful addition to the repertory of tools available to oceanographers.
Description
Submitted in partial fulfillment of the requirements for the degree of Doctor of Philosophy at the Massachusetts Institute of Technology and the Woods Hole Oceanographic Institution June 1999
Collections
Related items
Showing items related by title, author, creator and subject.

Error and uncertainty in estimates of Reynolds stress using ADCP in an energetic ocean state
Rapo, Mark A. (Massachusetts Institute of Technology and Woods Hole Oceanographic Institution, 200602)The challenge of estimating the Reynolds stress in an energetic ocean environment derives from the turbulence process overlapping in frequency, or in wavenumber, with the wave process. It was surmised that they would not ... 
Production and analysis of a Southern Ocean state estimate
Mazloff, Matthew R. (Massachusetts Institute of Technology and Woods Hole Oceanographic Institution, 200609)A modern general circulation model of the Southern Ocean with onesixth of a degree resolution is optimized to the observed ocean in a weighted least squares sense. Convergence toward the state estimate solution is carried ... 
On minimising assignment errors and the trade‐off between false positives and negatives in parentage analysis
Harrison, Hugo B.; SaenzAgudelo, Pablo; Planes, Serge; Jones, Geoffrey P.; Berumen, Michael L. (201305)Genetic parentage analyses provide a practical means with which to identify parent–offspring relationships in the wild. In Harrison et al.'s study (2013a), we compare three methods of parentage analysis and showed that the ...