Structure of the transport uncertainty in mesoscale inversions of CO 2 sources and sinks using ensemble model simulations

T. Lauvaux, O. Pannekoucke, C. Sarrat, F. Chevallier, P. Ciais, J. Noilhan, P. J. Rayner

Research output: Contribution to journalArticle

53 Citations (Scopus)

Abstract

We study the characteristics of a statistical ensemble of mesoscale simulations in order to estimate the model error in the simulation of CO 2 concentrations. The ensemble consists of ten members and the reference simulation using the operationnal short range forecast PEARP, perturbed using the Singular Vector technique. We then used this ensemble of simulations as the initial and boundary conditions for the meso scale model (Méso-NH) simulations, which uses CO2 fluxes from the ISBA-A-gs land surface model. The final ensemble represents the model dependence to the boundary conditions, conserving the physical properties of the dynamical schemes, but excluding the intrinsic error of the model. First, the variance of our ensemble is estimated over the domain, with associated spatial and temporal correlations. Second, we extract the signal from noisy horizontal correlations, due to the limited size ensemble, using diffusion equation modelling. The computational cost of such ensemble limits the number of members (simulations) especially when running online the carbon flux and the atmospheric models. In the theory, 50 to 100 members would be required to explore the overall sensitivity of the ensemble. The present diffusion model allows us to extract a significant part of the noisy error, and makes this study feasable with a limited number of simulations. Finally, we compute the diagonal and non-diagonal terms of the observation error covariance matrix and introduced it into our CO2 flux matrix inversion for 18 days of the 2005 intensive campaign CERES over the South West of France. Variances are based on model-data mismatch to ensure we treat model bias as well as ensemble dispersion, whereas spatial and temporal covariances are estimated with our method. The horizontal structure of the ensemble variance mani-fests the discontinuities of the mesoscale structures during the day, but remains locally driven during the night. On the vertical, surface layer variance shows large correlations with the upper levels in the boundary layer (> 0.6), dropping to 0.4 with the lower levels of the free troposphere. Large temporal correlations were found during the afternoon (> 0.5 for several hours), reduced during the night. The diffusion equation model extracted relevant error covariance signals horizontally, with reduced correlations over mountain areas and during the night over the continent. The posterior error reduction on the inverted CO2 fluxes accounting for the model error correlations illustrates the predominance of the temporal over the spatial correlations when using tower-based CO2 concentration observations.

Original languageEnglish (US)
Pages (from-to)1089-1102
Number of pages14
JournalBiogeosciences
Volume6
Issue number6
DOIs
StatePublished - Jan 1 2009

Fingerprint

simulation models
uncertainty
simulation
boundary condition
inversion
Manis
matrix
carbon flux
extracts
land surface
surface layer
troposphere
discontinuity
physical property
boundary layer
physical properties
France
mountains
mountain
cost

All Science Journal Classification (ASJC) codes

  • Ecology, Evolution, Behavior and Systematics
  • Earth-Surface Processes

Cite this

Lauvaux, T., Pannekoucke, O., Sarrat, C., Chevallier, F., Ciais, P., Noilhan, J., & Rayner, P. J. (2009). Structure of the transport uncertainty in mesoscale inversions of CO 2 sources and sinks using ensemble model simulations. Biogeosciences, 6(6), 1089-1102. https://doi.org/10.5194/bg-6-1089-2009
Lauvaux, T. ; Pannekoucke, O. ; Sarrat, C. ; Chevallier, F. ; Ciais, P. ; Noilhan, J. ; Rayner, P. J. / Structure of the transport uncertainty in mesoscale inversions of CO 2 sources and sinks using ensemble model simulations. In: Biogeosciences. 2009 ; Vol. 6, No. 6. pp. 1089-1102.
@article{d9345fff495747b89e2e1ea9f3b4c73b,
title = "Structure of the transport uncertainty in mesoscale inversions of CO 2 sources and sinks using ensemble model simulations",
abstract = "We study the characteristics of a statistical ensemble of mesoscale simulations in order to estimate the model error in the simulation of CO 2 concentrations. The ensemble consists of ten members and the reference simulation using the operationnal short range forecast PEARP, perturbed using the Singular Vector technique. We then used this ensemble of simulations as the initial and boundary conditions for the meso scale model (M{\'e}so-NH) simulations, which uses CO2 fluxes from the ISBA-A-gs land surface model. The final ensemble represents the model dependence to the boundary conditions, conserving the physical properties of the dynamical schemes, but excluding the intrinsic error of the model. First, the variance of our ensemble is estimated over the domain, with associated spatial and temporal correlations. Second, we extract the signal from noisy horizontal correlations, due to the limited size ensemble, using diffusion equation modelling. The computational cost of such ensemble limits the number of members (simulations) especially when running online the carbon flux and the atmospheric models. In the theory, 50 to 100 members would be required to explore the overall sensitivity of the ensemble. The present diffusion model allows us to extract a significant part of the noisy error, and makes this study feasable with a limited number of simulations. Finally, we compute the diagonal and non-diagonal terms of the observation error covariance matrix and introduced it into our CO2 flux matrix inversion for 18 days of the 2005 intensive campaign CERES over the South West of France. Variances are based on model-data mismatch to ensure we treat model bias as well as ensemble dispersion, whereas spatial and temporal covariances are estimated with our method. The horizontal structure of the ensemble variance mani-fests the discontinuities of the mesoscale structures during the day, but remains locally driven during the night. On the vertical, surface layer variance shows large correlations with the upper levels in the boundary layer (> 0.6), dropping to 0.4 with the lower levels of the free troposphere. Large temporal correlations were found during the afternoon (> 0.5 for several hours), reduced during the night. The diffusion equation model extracted relevant error covariance signals horizontally, with reduced correlations over mountain areas and during the night over the continent. The posterior error reduction on the inverted CO2 fluxes accounting for the model error correlations illustrates the predominance of the temporal over the spatial correlations when using tower-based CO2 concentration observations.",
author = "T. Lauvaux and O. Pannekoucke and C. Sarrat and F. Chevallier and P. Ciais and J. Noilhan and Rayner, {P. J.}",
year = "2009",
month = "1",
day = "1",
doi = "10.5194/bg-6-1089-2009",
language = "English (US)",
volume = "6",
pages = "1089--1102",
journal = "Biogeosciences",
issn = "1726-4170",
publisher = "European Geosciences Union",
number = "6",

}

Lauvaux, T, Pannekoucke, O, Sarrat, C, Chevallier, F, Ciais, P, Noilhan, J & Rayner, PJ 2009, 'Structure of the transport uncertainty in mesoscale inversions of CO 2 sources and sinks using ensemble model simulations', Biogeosciences, vol. 6, no. 6, pp. 1089-1102. https://doi.org/10.5194/bg-6-1089-2009

Structure of the transport uncertainty in mesoscale inversions of CO 2 sources and sinks using ensemble model simulations. / Lauvaux, T.; Pannekoucke, O.; Sarrat, C.; Chevallier, F.; Ciais, P.; Noilhan, J.; Rayner, P. J.

In: Biogeosciences, Vol. 6, No. 6, 01.01.2009, p. 1089-1102.

Research output: Contribution to journalArticle

TY - JOUR

T1 - Structure of the transport uncertainty in mesoscale inversions of CO 2 sources and sinks using ensemble model simulations

AU - Lauvaux, T.

AU - Pannekoucke, O.

AU - Sarrat, C.

AU - Chevallier, F.

AU - Ciais, P.

AU - Noilhan, J.

AU - Rayner, P. J.

PY - 2009/1/1

Y1 - 2009/1/1

N2 - We study the characteristics of a statistical ensemble of mesoscale simulations in order to estimate the model error in the simulation of CO 2 concentrations. The ensemble consists of ten members and the reference simulation using the operationnal short range forecast PEARP, perturbed using the Singular Vector technique. We then used this ensemble of simulations as the initial and boundary conditions for the meso scale model (Méso-NH) simulations, which uses CO2 fluxes from the ISBA-A-gs land surface model. The final ensemble represents the model dependence to the boundary conditions, conserving the physical properties of the dynamical schemes, but excluding the intrinsic error of the model. First, the variance of our ensemble is estimated over the domain, with associated spatial and temporal correlations. Second, we extract the signal from noisy horizontal correlations, due to the limited size ensemble, using diffusion equation modelling. The computational cost of such ensemble limits the number of members (simulations) especially when running online the carbon flux and the atmospheric models. In the theory, 50 to 100 members would be required to explore the overall sensitivity of the ensemble. The present diffusion model allows us to extract a significant part of the noisy error, and makes this study feasable with a limited number of simulations. Finally, we compute the diagonal and non-diagonal terms of the observation error covariance matrix and introduced it into our CO2 flux matrix inversion for 18 days of the 2005 intensive campaign CERES over the South West of France. Variances are based on model-data mismatch to ensure we treat model bias as well as ensemble dispersion, whereas spatial and temporal covariances are estimated with our method. The horizontal structure of the ensemble variance mani-fests the discontinuities of the mesoscale structures during the day, but remains locally driven during the night. On the vertical, surface layer variance shows large correlations with the upper levels in the boundary layer (> 0.6), dropping to 0.4 with the lower levels of the free troposphere. Large temporal correlations were found during the afternoon (> 0.5 for several hours), reduced during the night. The diffusion equation model extracted relevant error covariance signals horizontally, with reduced correlations over mountain areas and during the night over the continent. The posterior error reduction on the inverted CO2 fluxes accounting for the model error correlations illustrates the predominance of the temporal over the spatial correlations when using tower-based CO2 concentration observations.

AB - We study the characteristics of a statistical ensemble of mesoscale simulations in order to estimate the model error in the simulation of CO 2 concentrations. The ensemble consists of ten members and the reference simulation using the operationnal short range forecast PEARP, perturbed using the Singular Vector technique. We then used this ensemble of simulations as the initial and boundary conditions for the meso scale model (Méso-NH) simulations, which uses CO2 fluxes from the ISBA-A-gs land surface model. The final ensemble represents the model dependence to the boundary conditions, conserving the physical properties of the dynamical schemes, but excluding the intrinsic error of the model. First, the variance of our ensemble is estimated over the domain, with associated spatial and temporal correlations. Second, we extract the signal from noisy horizontal correlations, due to the limited size ensemble, using diffusion equation modelling. The computational cost of such ensemble limits the number of members (simulations) especially when running online the carbon flux and the atmospheric models. In the theory, 50 to 100 members would be required to explore the overall sensitivity of the ensemble. The present diffusion model allows us to extract a significant part of the noisy error, and makes this study feasable with a limited number of simulations. Finally, we compute the diagonal and non-diagonal terms of the observation error covariance matrix and introduced it into our CO2 flux matrix inversion for 18 days of the 2005 intensive campaign CERES over the South West of France. Variances are based on model-data mismatch to ensure we treat model bias as well as ensemble dispersion, whereas spatial and temporal covariances are estimated with our method. The horizontal structure of the ensemble variance mani-fests the discontinuities of the mesoscale structures during the day, but remains locally driven during the night. On the vertical, surface layer variance shows large correlations with the upper levels in the boundary layer (> 0.6), dropping to 0.4 with the lower levels of the free troposphere. Large temporal correlations were found during the afternoon (> 0.5 for several hours), reduced during the night. The diffusion equation model extracted relevant error covariance signals horizontally, with reduced correlations over mountain areas and during the night over the continent. The posterior error reduction on the inverted CO2 fluxes accounting for the model error correlations illustrates the predominance of the temporal over the spatial correlations when using tower-based CO2 concentration observations.

UR - http://www.scopus.com/inward/record.url?scp=72049096959&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=72049096959&partnerID=8YFLogxK

U2 - 10.5194/bg-6-1089-2009

DO - 10.5194/bg-6-1089-2009

M3 - Article

AN - SCOPUS:72049096959

VL - 6

SP - 1089

EP - 1102

JO - Biogeosciences

JF - Biogeosciences

SN - 1726-4170

IS - 6

ER -