07 Mar 2024
 | 07 Mar 2024
Status: this preprint is open for discussion.

Inverse modelling for surface methane flux estimation with 4DVar: impact of a computationally efficient representation of a non-diagonal B-matrix in INVICAT v4

Ross Noel Bannister and Chris Wilson

Abstract. Prior information is essential to most inverse problems and the surface flux estimation problem is no exception. The uncertainties of the prior fields, and their inter-correlations, should ideally be reflected in the a-priori error covariance matrix, often called B. The B-matrix, is however, difficult to quantify partly because it is typically a large matrix and partly because its numerical values are unknown.

We present a highly efficient method of representing the B-matrix to represent prior errors in the initial concentration and in the time sequence of surface fluxes for the 4DVar-based inverse modelling system (INVICAT) used to estimate the surface fluxes of methane. Our formulation is based on a spectral formulation of the square-root of B, which we believe has not been used in any such inverse modelling system before. It allows horizontal and vertical error correlations of the initial concentration, and horizontal and temporal error correlations of the flux to be represented. We provide full mathematical details. Our scheme allows the various correlation components to be switched on/off and for the respective length and timescales to be set in a way that is much more computationally efficient than representing such a B-matrix explicitly.

We test 14 configurations of the B-matrix (including the diagonal configuration) in a 100 day test assimilation of surface flask measurements of methane. We measure the performance of each by comparing the analysis to unassimilated observations held back for evaluation purposes. We find that the diagonal configuration is amongst the poorest performing choices of B. The best performing choice uses the spectral method. It does not include correlations for the initial concentration field, but does account for spatio-temporal correlations for the fluxes. These have the form of a SOAR (second order auto-regressive) function with a correlation length-scale of 600 km and a timescale of 3 months. Our results demonstrate the effectiveness of our method, which is applicable to very high resolution inverse modelling systems. We propose that potential biases in the prior initial condition field may be the reason for the poor performance when correlations in the prior initial concentration field are used.

Ross Noel Bannister and Chris Wilson

Status: open (until 08 May 2024)

Comment types: AC – author | RC – referee | CC – community | EC – editor | CEC – chief editor | : Report abuse
  • RC1: 'Comment on egusphere-2024-655', Anonymous Referee #1, 21 Mar 2024 reply
    • AC1: 'Reply on RC1', Ross Bannister, 01 Apr 2024 reply
Ross Noel Bannister and Chris Wilson

Model code and software

Inverse modelling for surface methane flux estimation with 4DVar: impact of a computationally efficient representation of a non-diagonal B-matrix Ross Bannister and Chris Wilson

Ross Noel Bannister and Chris Wilson


Total article views: 202 (including HTML, PDF, and XML)
HTML PDF XML Total BibTeX EndNote
144 46 12 202 4 4
  • HTML: 144
  • PDF: 46
  • XML: 12
  • Total: 202
  • BibTeX: 4
  • EndNote: 4
Views and downloads (calculated since 07 Mar 2024)
Cumulative views and downloads (calculated since 07 Mar 2024)

Viewed (geographical distribution)

Total article views: 207 (including HTML, PDF, and XML) Thereof 207 with geography defined and 0 with unknown origin.
Country # Views %
  • 1
Latest update: 12 Apr 2024
Short summary
Prior information is essential for the top-down estimation of CH4 surface fluxes. Errors in the prior are correlated in time/space, but accounting for correlations can be costly. We report on an efficient scheme to represent correlations in the inverse modelling system, INVICAT. The method is tested by assimilating CH4 observations using the scheme. Our findings show that accounting for spatio-temporal correlations improve CH4 flux estimates, demonstrating that the method should be further used.