Corrigendum

This article refers to:

Parametric Kalman filter for chemical transport models

Pannekoucke, O., Ricci, S., Bathelemy, S. Ménard, R., and Thual, O. (2016). Parametric Kalman filter for chemical transport models. Tellus A, 68, 31547.

https://doi.org/10.3402/tellusa.v68.31547

1. Introduction

In our previous contribution Pannekoucke et al. (2016) (P16), an error has been made in the derivation of the error diffusion tensor dynamics Eq. (20). This involves an error in the Lagrangian dynamics of the uncertainty given in the Algorithm 2, while leaving the Eulerian dyanmics unchanged.

The corrigendum is organized as follows. The modification of the Lagrangian dynamics is presented in Section 2, where a new version of the Algorithm 2 of P16 is presented. The computation leading to the Eulerian dynamics is described in Section 3, which gives the dynamics presented in Eq. (26) of P16.

2. Lagrangian dynamics of the diffusion tensor

The advection over a small time step δt can be viewed as equivalent of the deformation of the error field ε^b by the transformation $D(\text{x})=\text{x}+\text{u}(\text{x},t)\delta t$ (Pannekoucke et al., 2014). Hence, it follows that the metric tensor field ${\text{g}}_{\text{x}}(t)$ evolves in time as ${\widehat{\text{g}}}_{\text{x}}(t+\delta t)=D_{\text{x}}^{-T}{\text{g}}_{D^{-1}(\text{x})}(t)D_{\text{x}}^{-1}$, where $D_{\text{x}}^{-1}$ is the gradient of the inverse deformation $D^{-1}$ at $\text{x}$ (Pannekoucke et al., 2014, see Eq. (35)). This can be formulated considering the diffusion tensor $\nu$, defined by ${\nu }_{\text{x}}=\frac{1}{2}{\text{g}}_{\text{x}}^{-1}$, as: (1) $${⁁\nu }_x^b(t+\delta t)={(D_x^{-1})}^{-1}{\nu }_{D^{-1}(x)}^b(t){(D_x^{-T})}^{-1}.$$(1)

Eq. (1)(1) $${⁁\nu }_x^b(t+\delta t)={(D_x^{-1})}^{-1}{\nu }_{D^{-1}(x)}^b(t){(D_x^{-T})}^{-1}.$$(1) can be simplified as follows. Since the derivative of the identiy $D[D^{-1}(x)]=x$ is $D_{D^{-1}(x)}{D}_x^{-1}=I$, it results that ${(D_x^{-1})}^{-1}=D_{D^{-1}(x)}$. Hence, the dynamics of the diffusion tensor Eq. (1)(1) $${⁁\nu }_x^b(t+\delta t)={(D_x^{-1})}^{-1}{\nu }_{D^{-1}(x)}^b(t){(D_x^{-T})}^{-1}.$$(1) writes (2) $${⁁\nu }_x^b(t+\delta t)=D_{D^{-1}(x)}{\nu }_{D^{-1}(x)}^b(t){(D_{D^{-1}(x)})}^T.$$(2)

Considering Eq. (2)(2) $${⁁\nu }_x^b(t+\delta t)=D_{D^{-1}(x)}{\nu }_{D^{-1}(x)}^b(t){(D_{D^{-1}(x)})}^T.$$(2) , the Algorithm 2 of P16 now writes Algorithm 1.

Algorithm 1.

Iteration process to forecast the background covariance matrix at time $t=\tau$ from the analysis covariance matrix given at time t = 0, under local homogenity assumption.

Require: Fields of ${\nu }^a$ and $V^a$. $\delta t=\tau /N$, t = 0

 for $k=1:N$ do

  1- Pure advection $$D(x)=x+\text{u}(x,t)\delta t$$ $${⁁\nu }_x^b(t+\delta t)=D_{D^{-1}(x)}{\nu }_{D^{-1}(x)}^b(t){(D_{D^{-1}(x)})}^T$$ $${⁁\nu }_x^b(t+\delta t)=V^b[D^{-1}(x),t]$$

  2- Pure diffusion $${\nu }_x^b(t+\delta t)={⁁\nu }_x^b(t+\delta t)+2\kappa (x)\delta t$$ $$V^b{}_x(t+\delta t)={⁁\nu }_x^b(t+\delta t)\frac{|{⁁\nu }_x^b(t+\delta t){|}^{1/2}}{|{\nu }_x^b(t+\delta t){|}^{1/2}}$$

  3- Update of the background statistics $$V^b{}_x(t)\leftarrow V^b{}_x(t+\delta t)$$ $${\nu }_x^b(t)\leftarrow {\nu }_x^b(t+\delta t)$$ $$t\leftarrow t+\delta t$$  end for

 Return fields ${\nu }_x^b(\tau )$ and $V^b{}_x(\tau )$

3. Eulerian dynamics of the diffusion tensor

This expression modifies the derivation of the Eulerian dynamics in Appendix D of P16, where this time $$\begin{array}{l}{D}_{D^{-1}(x)}=D_{x-\delta t\text{u}(x,t)+o(\delta t)}=I+\delta t\nabla \text{u}(x-\delta t\text{u}(x,t)+o(\delta t)),\end{array}$$ leading to (3) $$D_{D^{-1}(x)}=I+\delta t\nabla \text{u}(x,t)+o(\delta t).$$(3)

With Eq. (3)(3) $$D_{D^{-1}(x)}=I+\delta t\nabla \text{u}(x,t)+o(\delta t).$$(3) , the computation in P16 leading to Eq. (D1) applies, and Eq. (D1) is found again: for the advection process, the dynamics of the error diffusion tensor writes: (4) $${\partial }_t{\nu }^b+\text{u}\nabla {\nu }^b={\nu }^b{(\nabla \text{u})}^T+(\nabla \text{u}){\nu }^b.$$(4)

4. Conclusion

In this corrigendum, the dynamics of the metric tensor and of the diffusion tensor have been corrected.

Algorithm 2 in Pannekoucke et al. (2016) is now replaced by the present Algorithm 1.

This modification does not alter the numerical results presented in P16.