Journal cover Journal topic
Nonlinear Processes in Geophysics An interactive open-access journal of the European Geosciences Union
Journal topic

Journal metrics

Journal metrics

  • IF value: 1.129 IF 1.129
  • IF 5-year value: 1.519 IF 5-year
  • CiteScore value: 1.54 CiteScore
  • SNIP value: 0.798 SNIP 0.798
  • SJR value: 0.610 SJR 0.610
  • IPP value: 1.41 IPP 1.41
  • h5-index value: 21 h5-index 21
  • Scimago H <br class='hide-on-tablet hide-on-mobile'>index value: 48 Scimago H
    index 48
Discussion papers
© Author(s) 2019. This work is distributed under
the Creative Commons Attribution 4.0 License.
© Author(s) 2019. This work is distributed under
the Creative Commons Attribution 4.0 License.

Research article 11 Mar 2019

Research article | 11 Mar 2019

Review status
This discussion paper is a preprint. It is a manuscript under review for the journal Nonlinear Processes in Geophysics (NPG).

Data assimilation using adaptive, non-conservative, moving mesh models

Ali Aydoğdu1, Alberto Carrassi1,2, Colin T. Guider3, Chris K. R. T. Jones3, and Pierre Rampal1 Ali Aydoğdu et al.
  • 1Nansen Environmental and Remote Sensing Center, Bergen, Norway
  • 2Geophysical Institute, University of Bergen, Norway
  • 3University of North Carolina, Chapel Hill, USA

Abstract. Numerical models solved on adaptive moving meshes have become increasingly prevalent in recent years. Motivating problems include the study of fluids in a Lagrangian frame and the presence of highly localized structures such as shock waves or interfaces. In the former case, Lagrangian solvers move the nodes of the mesh with the dynamical flow; in the latter, mesh resolution is increased in the proximity of the localized structure. Mesh adaptation can include remeshing, a procedure that adds or removes mesh nodes according to specific rules reflecting constraints in the numerical solver. In this case, the number of mesh nodes will change during the integration and, as a result, the dimension of the model’s state vector will not be conserved. This work presents a novel approach to the formulation of ensemble data assimilation for models with this underlying computational structure. The challenge lies in the fact that remeshing entails a different state space dimension across members of the ensemble, thus impeding the usual computation of consistent ensemble-based statistics. Our methodology adds one forward and one backward mapping step before and after the EnKF analysis respectively. This mapping takes all the ensemble members onto a fixed, uniform, reference mesh where the EnKF analysis can be performed. We consider a high- (HR) and a low-resolution (LR) fixed uniform reference mesh, whose resolutions are determined by the remeshing tolerances. This way the reference meshes embed the model numerical constraints and also are upper and lower uniform meshes bounding the resolutions of the individual ensemble meshes. Numerical experiments are carried out using 1D prototypical models: Burgers and Kuramoto-Sivashinsky equations, and both Eulerian and Lagrangian synthetic observations. While the HR strategy generally outperforms that of LR, their skill difference can be reduced substantially by an optimal tuning of the data assimilation parameters. The LR case is appealing in high-dimensions because of its lower computational burden. Lagrangian observations are shown to be very effective in that fewer of them are able to keep the analysis error at a level comparable to the more numerous observers for the Eulerian case. This study is motivated by the development of suitable EnKF strategies for 2D models of the sea-ice that are numerically solved on a Lagrangian mesh with remeshing.

Ali Aydoğdu et al.
Interactive discussion
Status: final response (author comments only)
Status: final response (author comments only)
AC: Author comment | RC: Referee comment | SC: Short comment | EC: Editor comment
[Login for Authors/Editors] [Subscribe to comment alert] Printer-friendly Version - Printer-friendly version Supplement - Supplement
Ali Aydoğdu et al.
Ali Aydoğdu et al.
Total article views: 297 (including HTML, PDF, and XML)
HTML PDF XML Total BibTeX EndNote
229 64 4 297 5 7
  • HTML: 229
  • PDF: 64
  • XML: 4
  • Total: 297
  • BibTeX: 5
  • EndNote: 7
Views and downloads (calculated since 11 Mar 2019)
Cumulative views and downloads (calculated since 11 Mar 2019)
Viewed (geographical distribution)  
Total article views: 196 (including HTML, PDF, and XML) Thereof 196 with geography defined and 0 with unknown origin.
Country # Views %
  • 1
No saved metrics found.
No discussed metrics found.
Latest update: 26 May 2019
Publications Copernicus
Short summary
Computational models involving adaptive meshes can both evolve dynamically and be remeshed. The remeshing means that the state vector dimension changes in time and across ensemble members, making EnKF unsuitable for the assimilation of observational data. We develop a modification in which the analysis is performed on a fixed uniform grid onto which the ensemble is mapped, with resolution related to the remeshing criteria. The approach is tested on two 1D models proving success.
Computational models involving adaptive meshes can both evolve dynamically and be remeshed. The...