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(x)=x+u(x,t)δt (Pannekoucke et al., 2014). Hence, it follows that the metric tensor field gx(t) evolves in time as g^x(t+δt)=Dx-TgD-1(x)(t)Dx-1, where Dx-1 is the gradient of the inverse deformation D-1 at x (Pannekoucke et al., 2014, see Eq. (35)). This can be formulated considering the diffusion tensor ν, defined by νx=12gx-1, as:

((1) )
νxb(t+δt)=(Dx-1)-1νD-1(x)b(t)(Dx-T)-1.

Eq. (1) can be simplified as follows. Since the derivative of the identiy D[D-1(x)]=x is DD-1(x)Dx-1=I, it results that (Dx-1)-1=DD-1(x). Hence, the dynamics of the diffusion tensor Eq. (1) writes

((2) )
νxb(t+δt)=DD-1(x)νD-1(x)b(t)(DD-1(x))T.

Considering Eq. (2), the Algorithm 2 of P16 now writes Algorithm 1.

Algorithm 1.

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

Require: Fields of νa and Va. δt=τ/N, t = 0

fork=1:Ndo

  1- Pure advection

D(x)=x+u(x,t)δt
νxb(t+δt)=DD-1(x)νD-1(x)b(t)(DD-1(x))T
νxb(t+δt)=Vb[D-1(x),t]

  2- Pure diffusion

νxb(t+δt)=νxb(t+δt)+2κ(x)δt
Vbx(t+δt)=νxb(t+δt)|νxb(t+δt)|1/2|νxb(t+δt)|1/2

  3- Update of the background statistics

Vbx(t)Vbx(t+δt)
νxb(t)νxb(t+δt)
tt+δt
end for

 Return fields νxb(τ) and Vbx(τ)

3.  

Eulerian dynamics of the diffusion tensor

This expression modifies the derivation of the Eulerian dynamics in Appendix D of P16, where this time

DD-1(x)=Dx-δtu(x,t)+o(δt)=I+δtu(x-δtu(x,t)+o(δt)),
leading to
((3) )
DD-1(x)=I+δtu(x,t)+o(δt).

With Eq. (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) )
tνb+uνb=νb(u)T+(u)νb.

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.