In recent years discontinuous Galerkin (DG) methods have received increased interest from the geophysical community. In these methods the solution in each grid cell is approximated as a linear combination of basis functions. Ensemble data assimilation (DA) aims to approximate the true state by combining model outputs with observations using error statistics estimated from an ensemble of model runs. Ensemble data assimilation in geophysical models faces several well-documented issues. In this work we exploit the expansion of the solution in DG basis functions to address some of these issues. Specifically, it is investigated whether a DA-DG combination (a) mitigates the need for observation thinning, (b) reduces errors in the field's gradients, and (c) can be used to set up scale-dependent localisation. Numerical experiments are carried out using stochastically generated ensembles of model states, with different noise properties, and with Legendre polynomials as basis functions. It is found that strong reduction in the analysis error is achieved by using DA-DG and that the benefit increases with increasing DG order. This is especially the case when small scales dominate the background error. The DA improvement in the first derivative is, on the other hand, marginal. We think this to be a counter-effect of the power of DG to fit the observations closely, which can deteriorate the estimates of the derivatives. Applying optimal localisation to the different polynomial orders, thus exploiting their different spatial length, is beneficial: it results in a covariance matrix closer to the true covariance than the matrix obtained using traditional optimal localisation in state space.

Pasmans I., Chen Y., Carrassi A., Jones C.K.R.T. (2024). Tailoring data assimilation to discontinuous Galerkin models. QUARTERLY JOURNAL OF THE ROYAL METEOROLOGICAL SOCIETY, 150(762), 2820-2847 [10.1002/qj.4737].

Tailoring data assimilation to discontinuous Galerkin models

Carrassi A.;
2024

Abstract

In recent years discontinuous Galerkin (DG) methods have received increased interest from the geophysical community. In these methods the solution in each grid cell is approximated as a linear combination of basis functions. Ensemble data assimilation (DA) aims to approximate the true state by combining model outputs with observations using error statistics estimated from an ensemble of model runs. Ensemble data assimilation in geophysical models faces several well-documented issues. In this work we exploit the expansion of the solution in DG basis functions to address some of these issues. Specifically, it is investigated whether a DA-DG combination (a) mitigates the need for observation thinning, (b) reduces errors in the field's gradients, and (c) can be used to set up scale-dependent localisation. Numerical experiments are carried out using stochastically generated ensembles of model states, with different noise properties, and with Legendre polynomials as basis functions. It is found that strong reduction in the analysis error is achieved by using DA-DG and that the benefit increases with increasing DG order. This is especially the case when small scales dominate the background error. The DA improvement in the first derivative is, on the other hand, marginal. We think this to be a counter-effect of the power of DG to fit the observations closely, which can deteriorate the estimates of the derivatives. Applying optimal localisation to the different polynomial orders, thus exploiting their different spatial length, is beneficial: it results in a covariance matrix closer to the true covariance than the matrix obtained using traditional optimal localisation in state space.
2024
Pasmans I., Chen Y., Carrassi A., Jones C.K.R.T. (2024). Tailoring data assimilation to discontinuous Galerkin models. QUARTERLY JOURNAL OF THE ROYAL METEOROLOGICAL SOCIETY, 150(762), 2820-2847 [10.1002/qj.4737].
Pasmans I.; Chen Y.; Carrassi A.; Jones C.K.R.T.
File in questo prodotto:
File Dimensione Formato  
Quart J Royal Meteoro Soc - 2024 - Pasmans - Tailoring data assimilation to discontinuous Galerkin models.pdf

accesso aperto

Tipo: Versione (PDF) editoriale
Licenza: Licenza per Accesso Aperto. Creative Commons Attribuzione (CCBY)
Dimensione 15.85 MB
Formato Adobe PDF
15.85 MB Adobe PDF Visualizza/Apri

I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.

Utilizza questo identificativo per citare o creare un link a questo documento: https://hdl.handle.net/11585/970437
Citazioni
  • ???jsp.display-item.citation.pmc??? ND
  • Scopus 0
  • ???jsp.display-item.citation.isi??? 0
social impact