2 Citations (Scopus)
100 Downloads (Pure)

Abstract

Strong constraint 4D-Variational data assimilation (4D-Var) is a method used to create an initialisation for a numerical model, that best replicates subsequent observations of the system it aims to recreate. The method does not take into account the presence of errors in the model, using the model equations as a strong constraint. This paper gives a rigorous and quantitative analysis of the errors introduced into the initialisation through the use of finite difference schemes to numerically solve the model equations. The 1D linear advection equation together with circulant boundary conditions, are chosen as the model equations of interest as they are representative of the advective processes relevant to numerical weather prediction, where 4D-Var is widely used. We consider the deterministic error introduced by finite difference approximations in the form of numerical dissipation and numerical dispersion and identify the relationship between these properties and the error in the 4D-Var initialisation. In particular, we find that a solely numerically dispersive scheme has the potential to introduce destructive interference resulting in the loss of some wavenumber components in the initialisation. Bounds for the error in the initialisation due to finite difference approximations are determined with and without observation errors. The bounds are found to depend on the smoothness of the true initial condition we wish to recover and the numerically dissipative and dispersive properties of the scheme. Numerical results are presented to demonstrate the effectiveness of the bounds. These lead to the conclusion that there exists a critical number of discretisation points when considering full sets of observations, where the effects of both the considered numerical model error and observational errors on the initialisation are minimised. The numerically dissipative and dispersive properties of the finite difference schemes also have the potential to alter the properties of the noise found in observations. Correlated noise structures may be introduced into the 4D-Var initialisation as a result. We determine when this occurs for observational errors in the form of additive white noise and find that the effect is reduced through the use of numerically non-dissipative finite difference schemes.
Original languageEnglish
Pages (from-to)567-588
JournalJournal of Computational and Applied Mathematics
Volume290
DOIs
Publication statusPublished - 15 Dec 2015

Fingerprint

Data Assimilation
Model Error
Initialization
Numerical models
Finite Difference Scheme
Finite Difference Approximation
Numerical Dispersion
Numerical Weather Prediction
Advection Equation
Correlated Noise
Quantitative Analysis
White noise
Advection
Model
Dissipation
Smoothness
Linear equation
Initial conditions
Interference
Discretization

Keywords

  • data assimilation
  • numerical model error
  • observation error
  • deterministic error

Cite this

The effect of numerical model error on data assimilation. / Jenkins, Sian; Smith, Nathan; Budd, Christopher; Freitag, Melina.

In: Journal of Computational and Applied Mathematics, Vol. 290, 15.12.2015, p. 567-588.

Research output: Contribution to journalArticle

Jenkins, Sian ; Smith, Nathan ; Budd, Christopher ; Freitag, Melina. / The effect of numerical model error on data assimilation. In: Journal of Computational and Applied Mathematics. 2015 ; Vol. 290. pp. 567-588.
@article{e27bcade0ae344dbb8d2eeb5dbb7ef9f,
title = "The effect of numerical model error on data assimilation",
abstract = "Strong constraint 4D-Variational data assimilation (4D-Var) is a method used to create an initialisation for a numerical model, that best replicates subsequent observations of the system it aims to recreate. The method does not take into account the presence of errors in the model, using the model equations as a strong constraint. This paper gives a rigorous and quantitative analysis of the errors introduced into the initialisation through the use of finite difference schemes to numerically solve the model equations. The 1D linear advection equation together with circulant boundary conditions, are chosen as the model equations of interest as they are representative of the advective processes relevant to numerical weather prediction, where 4D-Var is widely used. We consider the deterministic error introduced by finite difference approximations in the form of numerical dissipation and numerical dispersion and identify the relationship between these properties and the error in the 4D-Var initialisation. In particular, we find that a solely numerically dispersive scheme has the potential to introduce destructive interference resulting in the loss of some wavenumber components in the initialisation. Bounds for the error in the initialisation due to finite difference approximations are determined with and without observation errors. The bounds are found to depend on the smoothness of the true initial condition we wish to recover and the numerically dissipative and dispersive properties of the scheme. Numerical results are presented to demonstrate the effectiveness of the bounds. These lead to the conclusion that there exists a critical number of discretisation points when considering full sets of observations, where the effects of both the considered numerical model error and observational errors on the initialisation are minimised. The numerically dissipative and dispersive properties of the finite difference schemes also have the potential to alter the properties of the noise found in observations. Correlated noise structures may be introduced into the 4D-Var initialisation as a result. We determine when this occurs for observational errors in the form of additive white noise and find that the effect is reduced through the use of numerically non-dissipative finite difference schemes.",
keywords = "data assimilation, numerical model error, observation error, deterministic error",
author = "Sian Jenkins and Nathan Smith and Christopher Budd and Melina Freitag",
year = "2015",
month = "12",
day = "15",
doi = "10.1016/j.cam.2015.05.020",
language = "English",
volume = "290",
pages = "567--588",
journal = "Journal of Computational and Applied Mathematics",
issn = "0377-0427",
publisher = "Elsevier",

}

TY - JOUR

T1 - The effect of numerical model error on data assimilation

AU - Jenkins, Sian

AU - Smith, Nathan

AU - Budd, Christopher

AU - Freitag, Melina

PY - 2015/12/15

Y1 - 2015/12/15

N2 - Strong constraint 4D-Variational data assimilation (4D-Var) is a method used to create an initialisation for a numerical model, that best replicates subsequent observations of the system it aims to recreate. The method does not take into account the presence of errors in the model, using the model equations as a strong constraint. This paper gives a rigorous and quantitative analysis of the errors introduced into the initialisation through the use of finite difference schemes to numerically solve the model equations. The 1D linear advection equation together with circulant boundary conditions, are chosen as the model equations of interest as they are representative of the advective processes relevant to numerical weather prediction, where 4D-Var is widely used. We consider the deterministic error introduced by finite difference approximations in the form of numerical dissipation and numerical dispersion and identify the relationship between these properties and the error in the 4D-Var initialisation. In particular, we find that a solely numerically dispersive scheme has the potential to introduce destructive interference resulting in the loss of some wavenumber components in the initialisation. Bounds for the error in the initialisation due to finite difference approximations are determined with and without observation errors. The bounds are found to depend on the smoothness of the true initial condition we wish to recover and the numerically dissipative and dispersive properties of the scheme. Numerical results are presented to demonstrate the effectiveness of the bounds. These lead to the conclusion that there exists a critical number of discretisation points when considering full sets of observations, where the effects of both the considered numerical model error and observational errors on the initialisation are minimised. The numerically dissipative and dispersive properties of the finite difference schemes also have the potential to alter the properties of the noise found in observations. Correlated noise structures may be introduced into the 4D-Var initialisation as a result. We determine when this occurs for observational errors in the form of additive white noise and find that the effect is reduced through the use of numerically non-dissipative finite difference schemes.

AB - Strong constraint 4D-Variational data assimilation (4D-Var) is a method used to create an initialisation for a numerical model, that best replicates subsequent observations of the system it aims to recreate. The method does not take into account the presence of errors in the model, using the model equations as a strong constraint. This paper gives a rigorous and quantitative analysis of the errors introduced into the initialisation through the use of finite difference schemes to numerically solve the model equations. The 1D linear advection equation together with circulant boundary conditions, are chosen as the model equations of interest as they are representative of the advective processes relevant to numerical weather prediction, where 4D-Var is widely used. We consider the deterministic error introduced by finite difference approximations in the form of numerical dissipation and numerical dispersion and identify the relationship between these properties and the error in the 4D-Var initialisation. In particular, we find that a solely numerically dispersive scheme has the potential to introduce destructive interference resulting in the loss of some wavenumber components in the initialisation. Bounds for the error in the initialisation due to finite difference approximations are determined with and without observation errors. The bounds are found to depend on the smoothness of the true initial condition we wish to recover and the numerically dissipative and dispersive properties of the scheme. Numerical results are presented to demonstrate the effectiveness of the bounds. These lead to the conclusion that there exists a critical number of discretisation points when considering full sets of observations, where the effects of both the considered numerical model error and observational errors on the initialisation are minimised. The numerically dissipative and dispersive properties of the finite difference schemes also have the potential to alter the properties of the noise found in observations. Correlated noise structures may be introduced into the 4D-Var initialisation as a result. We determine when this occurs for observational errors in the form of additive white noise and find that the effect is reduced through the use of numerically non-dissipative finite difference schemes.

KW - data assimilation

KW - numerical model error

KW - observation error

KW - deterministic error

UR - http://dx.doi.org/10.1016/j.cam.2015.05.020

U2 - 10.1016/j.cam.2015.05.020

DO - 10.1016/j.cam.2015.05.020

M3 - Article

VL - 290

SP - 567

EP - 588

JO - Journal of Computational and Applied Mathematics

JF - Journal of Computational and Applied Mathematics

SN - 0377-0427

ER -