Joint CO 2 state and flux estimation with the 4D-Var system EURAD-IM

Abstract. This idealized regional atmospheric inversion study assesses the potential of the 4-dimensional variational (4D-Var) method to estimate CO 2 fluxes and the atmospheric CO 2 concentration state jointly. In order to distinguish and quantify the surface-atmosphere CO 2 fluxes, combining anthropogenic CO 2 emissions, photosynthesis, and respiration, we include uncertainties of initial values, which arise from highly uncertain surface fluxes and night-time transport. Therefor a new calculation of the background error standard deviation for the CO 2 fluxes was developed. To suppress spurious wiggles occurring from advection, an absolute monotone advection scheme with low numeric diffusion and its adjoint has been implemented. The inversion by the EURopean Air pollution Dispersion-Inverse Model (EURAD-IM) with 5 km resolution in Central Europe is validated by synthetic half hourly measurements from eleven concentration towers. A significant improvement of the analysis is shown if initial values and CO 2 fluxes are optimised jointly, compared to optimising CO 2 fluxes alone, without estimating uncertainty of atmospheric concentration. We find that joint estimation of carbon fluxes and initial states requires a careful balance of the background error covariance matrices but enables a more detailed analysis of atmospheric CO 2 and the surface-atmosphere fluxes.