Abstract:Discontinuous Galerkin methods approximate the model solution as a linear combination of basis functions in each grid cell. For ensemble data assimilation this provides an intrinsic way of decomposing the background error covariance (a) by scale (b) allowing for scale‐dependent localisation. Furthermore, the basis functions resolve the solution at subgrid scales, permitting more effective assimilation of dense observation sets. 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.

A structurally informed data assimilation approach for nonlinear partial differential equations

A Structurally Informed Data Assimilation Approach for Nonlinear Partial Differential Equations

Adaptive Multifidelity Data Assimilation for Nonlinear Subsurface Flow Problems

A Closed-Form Nonlinear Data Assimilation Algorithm for Multi-Layer Flow Fields

Adaptive Multi-Fidelity Probabilistic Collocation-Based Kalman Filter for Subsurface Flow Data Assimilation: Numerical Modeling and Real-World Experiment

A Stochastic Collocation Based Kalman Filter for Data Assimilation

Data Assimilation for Unsaturated Flow Models with Restart Adaptive Probabilistic Collocation Based Kalman Filter

Data assimilation for nonlinear problems by ensemble Kalman filter with reparameterization

Statistical-Learning-Based Ensemble Data Assimilation Methods for Parameter Estimation in Hydrodynamic Models

Introduction to data assimilation for parameter estimation

A Modified Iterative Ensemble Kalman Filter Data Assimilation Method

Numerical Linear Algebra in Data Assimilation

Spatial estimation of unidirectional wave evolution based on ensemble data assimilation

Data Assimilation for Transient Flow in Geologic Formations Via Ensemble Kalman Filter

A Localized Ensemble Kalman Smoother

Preserving Nonlinear Constraints in Variational Flow Filtering Data Assimilation

Iterative ensemble smoothers for data assimilation in coupled nonlinear multiscale models

Importance Weighting in Hybrid Iterative Ensemble Smoothers for Data Assimilation

Learning Dynamical Characteristics with Neural Operators for Data Assimilation

Sequential data assimilation with multiple nonlinear models and applications to subsurface flow

Tailoring data assimilation to discontinuous Galerkin models